The genetic basis of a complex functional system

Many-to-many mapping: A simulation study of how the number of traits and tasks affect the evolution of form and function

Many-to-many mapping of form-to-function posits that multiple morphological and physiological traits affect the performance of multiple tasks in an organism, and that redundancy and multitasking occur simultaneously to shape the evolution of an organism's phenotype. Many-to-many mapping is expected to be ubiquitous in nature, yet little is known about how it influences the evolution of organismal phenotype. The F-matrix is a powerful tool to study these issues because it describes how multiple traits affect multiple tasks. We undertook a simulation study using the F-matrix to test how the number of traits and the number of tasks affect trait integration and evolvability, as well as the relationships among tasks. We found that as the number of traits and/or tasks increases, the relationships between the tasks and the integration between the traits become weaker, and that the evolvability of the traits increases, all resulting in a system that is freer to evolve. We also found that as the number of traits increases, performance tradeoffs tend to become weaker, but only to a point. Our work shows that it is important to consider not only multiple traits, but also the multitude of tasks that those traits carry out when studying form-function relationships. We suggest that evolution of these relationships follows functional lines of least resistance, which are less defined in more complex systems, resulting in a mechanism for diversification.

Genetic architecture of adaptive radiation across two trophic levels

Evolution of trophic diversity is a hallmark of adaptive radiation. Yet, transitions between carnivory and herbivory are rare in young adaptive radiations. Haplochromine cichlid fish of the African Great Lakes are exceptional in this regard. Lake Victoria was colonized by an insectivorous generalist and in less than 20 000 years, several clades of specialized herbivores evolved. Carnivorous versus herbivorous lifestyles in cichlids require many different adaptations in functional morphology, physiology and behaviour. Ecological transitions in either direction thus require many traits to change in a concerted fashion, which could be facilitated if genomic regions underlying these traits were physically linked or pleiotropic. However, linkage/pleiotropy could also constrain evolvability. To investigate components of the genetic architecture of a suite of traits that distinguish invertivores from algae scrapers, we performed quantitative trait locus (QTL) mapping using a second-generation hybrid cross. While we found indications of linkage/pleiotropy within trait complexes, QTLs for distinct traits were distributed across several unlinked genomic regions. Thus, a mixture of independently segregating variation and some pleiotropy may underpin the rapid trophic transitions. We argue that the emergence and maintenance of associations between the different genomic regions underpinning co-adapted traits that evolved and persist against some gene flow required reproductive isolation.

Assessing the Levels of Functional Adaptation: Finite Element Analysis Reveals Species, Hybrid, and Sexual Variation in the Biomechanics of African Cichlid Mandibles

To understand how adaptive divergence emerges it is essential to examine the function of phenotypic traits along a continuum. For vertebrates, the mandible provides a key link with foraging and other important activities which has made it highly relevant for investigations of biomechanical change. Variation in mandible shape is known to correspond with ecology but its function is often only investigated between distinct species. However, for such divergence to occur and be maintained selection likely draws from many sources of biomechanical variation. African cichlids represent an exemplar model for understanding how such processes unfold with mandible variation existing between species, sexes, and is likely generated in nature by the potential for hybridization. We explored such mandible variation through a finite element modelling approach and predicted that hybrids and females would have reduced functional capabilities, the former in line with disruptive selection and the latter due to potential trade-offs incurred by maternal mouthbrooding in Malawian haplochromines. We revealed evidence of structural adaptations between Tropheops ‘Red Cheek’ and Labeotrophues fuelleborni that impacted the dispersion of mechanical stress in ways that matched the foraging of these species. Also, hybrids showed higher stresses relative to both species across the mandible. Sexual dimorphism in stress handling was evident despite minor differences in shape with males showing enhanced load resistance. However, in hybrids it appeared that males were disadvantaged relative to females, and displayed asymmetry in load handling. Together, these results show evidence of species and sex based biomechanical variation, that could be targeted by divergent selection.

OUP accepted manuscript

Whole genome sequences are beginning to revolutionize our understanding of phylogenetic relationships. Yet, even whole genome sequences can fail to resolve the evolutionary history of the most rapidly radiating lineages, where incomplete lineage sorting, standing genetic variation, introgression, and other factors obscure the phylogenetic history of the group. To overcome such challenges, one emerging strategy is to integrate results across different methods. Most such approaches have been implemented on reduced representation genomic data sets, but whole genomes should provide the maximum possible evidence approach. Here, we test the ability of single nucleotide polymorphisms extracted from whole genome resequencing data, implemented in an integrative genomic approach, to resolve key nodes in the phylogeny of the mbuna, rock-dwelling cichlid fishes of Lake Malaŵi, which epitomize the phylogenetic intractability that often accompanies explosive lineage diversification. This monophyletic radiation has diversified at an unparalleled rate into several hundred species in less than 2 million years. Using an array of phylogenomic methods, we consistently recovered four major clades of mbuna, but a large basal polytomy among them. Although introgression between clades apparently contributed to the challenge of phylogenetic reconstruction, reduction of the data set to nonintrogressed sites still did not help to resolve the basal polytomy. On the other hand, relationships among six congeneric species pairs were resolved without ambiguity, even in one case where existing data led us to predict that resolution would be difficult. We conclude that the bursts of diversification at the earliest stages of the mbuna radiation may be phylogenetically unresolvable, but other regions of the tree are phylogenetically clearly supported. Integration of multiple phylogenomic approaches will continue to increase confidence in relationships inferred from these and other whole-genome data sets. [Incomplete lineage sorting; introgression; linkage disequilibrium; multispecies coalescence; rapid radiation; soft polytomy.]

Ciliary Rootlet Coiled-Coil 2 (crocc2) Is Associated with Evolutionary Divergence and Plasticity of Cichlid Jaw Shape

Cichlid fishes exhibit rapid, extensive, and replicative adaptive radiation in feeding morphology. Plasticity of the cichlid jaw has also been well documented, and this combination of iterative evolution and developmental plasticity has led to the proposition that the cichlid feeding apparatus represents a morphological "flexible stem." Under this scenario, the fixation of environmentally sensitive genetic variation drives evolutionary divergence along a phenotypic axis established by the initial plastic response. Thus, if plasticity is predictable then so too should be the evolutionary response. We set out to explore these ideas at the molecular level by identifying genes that underlie both the evolution and plasticity of the cichlid jaw. As a first step, we fine-mapped an environment-specific quantitative trait loci for lower jaw shape in cichlids, and identified a nonsynonymous mutation in the ciliary rootlet coiled-coil 2 (crocc2), which encodes a major structural component of the primary cilium. Given that primary cilia play key roles in skeletal mechanosensing, we reasoned that this gene may confer its effects by regulating the sensitivity of bone to respond to mechanical input. Using both cichlids and zebrafish, we confirmed this prediction through a series of experiments targeting multiple levels of biological organization. Taken together, our results implicate crocc2 as a novel mediator of bone formation, plasticity, and evolution.

Surprising spatiotemporal stability of a multi-peak fitness landscape revealed by independent field experiments measuring hybrid fitness

The effect of the environment on fitness in natural populations is a fundamental question in evolutionary biology. However, experimental manipulations of both environment and phenotype at the same time are rare. Thus, the relative importance of the competitive environment versus intrinsic organismal performance in shaping the location, height, and fluidity of fitness peaks and valleys remains largely unknown. Here, we experimentally tested the effect of competitor frequency on the complex fitness landscape driving adaptive radiation of a generalist and two trophic specialist pupfishes, a scale-eater and molluscivore, endemic to hypersaline lakes on San Salvador Island (SSI), Bahamas. We manipulated phenotypes, by generating 3407 F4/F5 lab-reared hybrids, and competitive environment, by altering the frequency of rare transgressive hybrids between field enclosures in two independent lake populations. We then tracked hybrid survival and growth rates across these four field enclosures for 3-11 months. In contrast to competitive speciation theory, we found no evidence that the frequency of hybrid phenotypes affected their survival. Instead, we observed a strikingly similar fitness landscape to a previous independent field experiment, each supporting multiple fitness peaks for generalist and molluscivore phenotypes and a large fitness valley isolating the divergent scale-eater phenotype. These features of the fitness landscape were stable across manipulated competitive environments, multivariate trait axes, and spatiotemporal heterogeneity. We suggest that absolute performance constraints and divergent gene regulatory networks shape macroevolutionary (interspecific) fitness landscapes in addition to microevolutionary (intraspecific) competitive dynamics. This interplay between organism and environment underlies static and dynamic features of the adaptive landscape.

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.

Rapid generation of ecologically relevant behavioral novelty in experimental cichlid hybrids

The East African cichlid radiations are characterized by repeated and rapid diversification into many distinct species with different ecological specializations and by a history of hybridization events between nonsister species. Such hybridization might provide important fuel for adaptive radiation. Interspecific hybrids can have extreme trait values or novel trait combinations and such transgressive phenotypes may allow some hybrids to explore ecological niches neither of the parental species could tap into. Here, we investigate the potential of second-generation (F2) hybrids between two generalist cichlid species from Lake Malawi to exploit a resource neither parental species is specialized on: feeding by sifting sand. Some of the F2 hybrids phenotypically resembled fish of species that are specialized on sand sifting. We combined experimental behavioral and morphometric approaches to test whether the F2 hybrids are transgressive in both morphology and behavior related to sand sifting. We then performed a quantitative trait loci (QTL) analysis using RADseq markers to investigate the genetic architecture of morphological and behavioral traits. We show that transgression is present in several morphological traits, that novel trait combinations occur, and we observe transgressive trait values in sand sifting behavior in some of the F2 hybrids. Moreover, we find QTLs for morphology and for sand sifting behavior, suggesting the existence of some loci with moderate to large effects. We demonstrate that hybridization has the potential to rapidly generate novel and ecologically relevant phenotypes that may be suited to a niche neither of the parental species occupies. Interspecific hybridization may thereby contribute to the rapid generation of ecological diversity in cichlid radiations.

Evolution of a soft-tissue foraging adaptation in African cichlids: Roles for novelty, convergence, and constraint

Understanding the origins of biodiversity demands consideration of both extrinsic (e.g., ecological opportunity) and intrinsic (e.g., developmental constraint) factors. Here, we use a combination of phylogenetic and genetic tools to address the origin of novelty in African cichlids. In particular, we focus on an extreme hypertrophied snout that is structurally integrated with the upper jaw. We show that this bizarre trait has evolved independently in at least two distinct and ecologically successful cichlid clades. We find that snout dimensions are decoupled both phenotypically and genetically, which has enabled it to evolve independently in multiple directions. Further, patterns of variation among species and within a genetic mapping pedigree suggest that relative to snout length, depth is under greater genetic and/or developmental constraint. Models of evolution suggest that snout shape is under selection for feeding behavior, with snout depth being important for algae scraping and snout length for sand sifting. Indeed, the deep snout of some algivores is achieved via an expansion of the intermaxillary ligament, which is important for jaw stability and may increase feeding performance. Overall, our data imply that the evolution of exaggerated snout depth required overcoming a genetic/developmental constraint, which led to expanded ecological opportunity via foraging adaptation.

Functional coupling in the evolution of suction feeding and gill ventilation of sculpins (Perciformes: Cottoidei)

Suction feeding and gill ventilation in teleosts are functionally coupled, meaning that there is an overlap in the structures involved with both functions. Functional coupling is one type of morphological integration, a term that broadly refers to any covariation, correlation, or coordination among structures. Suction feeding and gill ventilation exhibit other types of morphological integration, including functional coordination (a tendency of structures to work together to perform a function) and evolutionary integration (a tendency of structures to covary in size or shape across evolutionary history). Functional coupling, functional coordination, and evolutionary integration have each been proposed to limit morphological diversification to some extent. Yet teleosts show extraordinary cranial diversity, suggesting that there are mechanisms within some teleost clades that promote morphological diversification, even within the highly integrated suction feeding and gill ventilatory systems. To investigate this, we quantified evolutionary integration among four mechanical units associated with suction feeding and gill ventilation in a diverse clade of benthic, primarily suction-feeding fishes (Cottoidei; sculpins and relatives). We reconstructed cottoid phylogeny using molecular data from 108 species, and obtained 24 linear measurements of four mechanical units (jaws, hyoid, opercular bones, and branchiostegal rays) from micro-CT reconstructions of 44 cottoids and 1 outgroup taxon. We tested for evolutionary correlation and covariation among the four mechanical units using phylogenetically corrected principal component analysis to reduce the dimensionality of measurements for each unit, followed by correlating phylogenetically independent contrasts and computing phylogenetic generalized least squares models from the first principle component axis of each of the four mechanical units. The jaws, opercular bones, and branchiostegal rays show evolutionary integration, but the hyoid is not positively integrated with these units. To examine these results in an ecomorphological context, we used published ecological data in phylogenetic ANOVA models to demonstrate that the jaw is larger in fishes that eat elusive or grasping prey (e.g., prey that can easily escape or cling to the substrate) and that the hyoid is smaller in intertidal and hypoxia-tolerant sculpins. Within Cottoidei, the relatively independent evolution of the hyoid likely has reduced limitations on morphological evolution within the highly morphologically integrated suction feeding and gill ventilatory systems.

Pleiotropic jaw morphology links the evolution of mechanical modularity and functional feeding convergence in Lake Malawi cichlids

Complexity in how mechanistic variation translates into ecological novelty could be critical to organismal diversification. For instance, when multiple distinct morphologies can generate the same mechanical or functional phenotype, this could mitigate trade-offs and/or provide alternative ways to meet the same ecological challenge. To investigate how this type of complexity shapes diversity in a classic adaptive radiation, we tested several evolutionary consequences of the anterior jaw four-bar linkage for Lake Malawi cichlid trophic diversification. Using a novel phylogenetic framework, we demonstrated that different mechanical outputs of the same four jaw elements are evolutionarily associated with both jaw protrusion distance and jaw protrusion angle. However, these two functional aspects of jaw protrusion have evolved independently. Additionally, although four-bar morphology showed little evidence for attraction to optima, there was substantial evidence of adaptive peaks for emergent four-bar linkage mechanics and jaw protrusion abilities among Malawi feeding guilds. Finally, we highlighted a clear case of two cichlid species that have -independently evolved to graze algae in less than 2 Myr and have converged on similar jaw protrusion abilities as well as four-bar linkage mechanics, but have evolved these similarities via non-convergent four-bar morphologies.

Repeated Evolution Versus Common Ancestry: Sex Chromosome Evolution in the Haplochromine Cichlid Pseudocrenilabrus philander

Why sex chromosomes turn over and remain undifferentiated in some taxa, whereas they degenerate in others, is still an area of ongoing research. The recurrent occurrence of homologous and homomorphic sex chromosomes in distantly related taxa suggests their independent evolution or continued recombination since their first emergence. Fishes display a great diversity of sex-determining systems. Here, we focus on sex chromosome evolution in haplochromines, the most species-rich lineage of cichlid fishes. We investigate sex-specific signatures in the Pseudocrenilabrus philander species complex, which belongs to a haplochromine genus found in many river systems and ichthyogeographic regions in northern, eastern, central, and southern Africa. Using whole-genome sequencing and population genetic, phylogenetic, and read-coverage analyses, we show that one population of P. philander has an XX-XY sex-determining system on LG7 with a large region of suppressed recombination. However, in a second bottlenecked population, we did not find any sign of a sex chromosome. Interestingly, LG7 also carries an XX-XY system in the phylogenetically more derived Lake Malawi haplochromine cichlids. Although the genomic regions determining sex are the same in Lake Malawi cichlids and P. philander, we did not find evidence for shared ancestry, suggesting that LG7 evolved as sex chromosome at least twice in haplochromine cichlids. Hence, our work provides further evidence for the labile nature of sex determination in fishes and supports the hypothesis that the same genomic regions can repeatedly and rapidly be recruited as sex chromosomes in more distantly related lineages.

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.

The Integrated Genomic Architecture and Evolution of Dental Divergence in East African Cichlid Fishes (Haplochromis chilotes x H. nyererei)

The independent evolution of the two toothed jaws of cichlid fishes is thought to have promoted their unparalleled ecological divergence and species richness. However, dental divergence in cichlids could exhibit substantial genetic covariance and this could dictate how traits like tooth numbers evolve in different African Lakes and on their two jaws. To test this hypothesis, we used a hybrid mapping cross of two trophically divergent Lake Victoria species (Haplochromis chilotes × Haplochromis nyererei) to examine genomic regions associated with cichlid tooth diversity. Surprisingly, a similar genomic region was found to be associated with oral jaw tooth numbers in cichlids from both Lake Malawi and Lake Victoria. Likewise, this same genomic location was associated with variation in pharyngeal jaw tooth numbers. Similar relationships between tooth numbers on the two jaws in both our Victoria hybrid population and across the phylogenetic diversity of Malawi cichlids additionally suggests that tooth numbers on the two jaws of haplochromine cichlids might generally coevolve owing to shared genetic underpinnings. Integrated, rather than independent, genomic architectures could be key to the incomparable evolutionary divergence and convergence in cichlid tooth numbers.

Genetic dissection of adaptive form and function in rapidly speciating cichlid fishes

Genes of major phenotypic effects and strong genetic correlations can facilitate adaptation, direct selective responses, and potentially lead to phenotypic convergence. However, the preponderance of this type of genetic architecture in repeatedly evolved adaptations remains unknown. Using hybrids between Haplochromis chilotes (thick-lipped) and Pundamilia nyererei (thin-lipped) we investigated the genetics underlying hypertrophied lips and elongated heads, traits that evolved repeatedly in cichlids. At least 25 loci of small-to-moderate and mainly additive effects were detected. Phenotypic variation in lip and head morphology was largely independent. Although several QTL overlapped for lip and head morphology traits, they were often of opposite effects. The distribution of effect signs suggests strong selection on lips. The fitness implications of several detected loci were demonstrated using a laboratory assay testing for the association between genotype and variation in foraging performance. The persistence of low fitness alleles in head morphology appears to be maintained through antagonistic pleiotropy/close linkage with positive-effect lip morphology alleles. Rather than being based on few major loci with strong positive genetic correlations, our results indicate that the evolution of the Lake Victoria thick-lipped ecomorph is the result of selection on numerous loci distributed throughout the genome.

Mechanical Transgressive Segregation and the Rapid Origin of Trophic Novelty

Hybrid phenotypes are often intermediate between those of parental species. However, hybridization can generate novel phenotypes when traits are complex. For instance, even when the morphologies of individual musculo-skeletal components do not segregate outside the parental range in hybrid offspring, complex functional systems can exhibit emergent phenotypes whose mechanics exceed the parental values. To determine if transgression in mechanics could facilitate divergence during an adaptive radiation, we examined three functional systems in the trophic apparatus of Lake Malawi cichlid fishes. We conducted a simulation study of hybridization between species pairs whose morphology for three functional systems was empirically measured, to determine how the evolutionary divergence of parental species influences the frequency that hybridization could produce mechanics that transgress the parental range. Our simulations suggest that the complex mechanical systems of the cichlid trophic apparatus commonly exhibit greater transgression between more recently diverged cichlid species. Because (1) all three mechanical systems produce hybrids with transgressive mechanics in Lake Malawi cichlids, (2) hybridization is common, and (3) single hybrid crosses often recapitulate a substantial diversity of mechanics, we conclude that mechanical transgressive segregation could play an important role in the rapid accumulation of phenotypic variation in adaptive radiations.

Variations on a theme: Genomics of sex determination in the cichlid fish Astatotilapia burtoni

Background

Sex chromosomes change more frequently in fish than in mammals or birds. However, certain chromosomes or genes are repeatedly used as sex determinants in different members of the teleostean lineage. East African cichlids are an enigmatic model system in evolutionary biology representing some of the most diverse extant vertebrate adaptive radiations. How sex is determined and if different sex-determining mechanisms contribute to speciation is unknown for almost all of the over 1,500 cichlid species of the Great Lakes. Here, we investigated the genetic basis of sex determination in a cichlid from Lake Tanganyika, Astatotilapia burtoni, a member of the most species-rich cichlid lineage, the haplochromines.

Results

We used RAD-sequencing of crosses for two populations of A. burtoni, a lab strain and fish caught at the south of Lake Tanganyika. Using association mapping and comparative genomics, we confirmed male heterogamety in A. burtoni and identified different sex chromosomes (LG5 and LG18) in the two populations of the same species. LG5, the sex chromosome of the lab strain, is a fusion chromosome in A. burtoni. Wnt4 is located on this chromosome, representing the best candidate identified so far for the master sex-determining gene in our lab strain of A. burtoni.

Conclusions

Cichlids exemplify the high turnover rate of sex chromosomes in fish with two different chromosomes, LG5 and LG18, containing major sex-determining loci in the two populations of A. burtoni examined here. However, they also illustrate that particular chromosomes are more likely to be used as sex chromosomes. Chromosome 5 is such a chromosome, which has evolved several times as a sex chromosome, both in haplochromine cichlids from all Great Lakes and also in other teleost fishes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3178-0) contains supplementary material, which is available to authorized users.

Foraging environment determines the genetic architecture and evolutionary potential of trophic morphology in cichlid fishes

Phenotypic plasticity allows organisms to change their phenotype in response to shifts in the environment. While a central topic in current discussions of evolutionary potential, a comprehensive understanding of the genetic underpinnings of plasticity is lacking in systems undergoing adaptive diversification. Here, we investigate the genetic basis of phenotypic plasticity in a textbook adaptive radiation, Lake Malawi cichlid fishes. Specifically, we crossed two divergent species to generate an F₃ hybrid mapping population. At early juvenile stages, hybrid families were split and reared in alternate foraging environments that mimicked benthic/scraping or limnetic/sucking modes of feeding. These alternate treatments produced a variation in morphology that was broadly similar to the major axis of divergence among Malawi cichlids, providing support for the flexible stem theory of adaptive radiation. Next, we found that the genetic architecture of several morphological traits was highly sensitive to the environment. In particular, of 22 significant quantitative trait loci (QTL), only one was shared between the environments. In addition, we identified QTL acting across environments with alternate alleles being differentially sensitive to the environment. Thus, our data suggest that while plasticity is largely determined by loci specific to a given environment, it may also be influenced by loci operating across environments. Finally, our mapping data provide evidence for the evolution of plasticity via genetic assimilation at an important regulatory locus, ptch1. In all, our data address long-standing discussions about the genetic basis and evolution of plasticity. They also underscore the importance of the environment in affecting developmental outcomes, genetic architectures, morphological diversity and evolutionary potential.

Do relaxed selection and habitat temperature facilitate biased mitogenomic introgression in a narrowly endemic fish?

Introgression might be exceptionally common during the evolution of narrowly endemic species. For instance, in the springs of the small and isolated Cuatro Ciénegas Valley, the mitogenome of the cichlid fish Herichthys cyanoguttatus could be rapidly introgressing into populations of the trophically polymorphic H. minckleyi. We used a combination of genetic and environmental data to examine the factors associated with this mitochondrial introgression. A reduced representation library of over 6220 single nucleotide polymorphisms (SNPs) from the nuclear genome showed that mitochondrial introgression into H. minckleyi is biased relative to the amount of nuclear introgression. SNP assignment probabilities also indicated that cichlids with more hybrid ancestry are not more commonly female providing no support for asymmetric backcrossing or hybrid‐induced sex‐ratio distortion in generating the bias in mitochondrial introgression. Smaller effective population size in H. minckleyi inferred from the SNPs coupled with sequences of all 13 mitochondrial proteins suggests that relaxed selection on the mitogenome could be facilitating the introgression of “H. cyanoguttatus” haplotypes. Additionally, we showed that springs with colder temperatures had greater amounts of mitochondrial introgression from H. cyanoguttatus. Relaxed selection in H. minckleyi coupled with temperature‐related molecular adaptation could be facilitating mitogenomic introgression into H. minckleyi.

Coevolutionary patterning of teeth and taste buds

Teeth and taste buds are iteratively patterned structures that line the oro-pharynx of vertebrates. Biologists do not fully understand how teeth and taste buds develop from undifferentiated epithelium or how variation in organ density is regulated. These organs are typically studied independently because of their separate anatomical location in mammals: teeth on the jaw margin and taste buds on the tongue. However, in many aquatic animals like bony fishes, teeth and taste buds are colocalized one next to the other. Using genetic mapping in cichlid fishes, we identified shared loci controlling a positive correlation between tooth and taste bud densities. Genome intervals contained candidate genes expressed in tooth and taste bud fields. sfrp5 and bmper, notable for roles in Wingless (Wnt) and bone morphogenetic protein (BMP) signaling, were differentially expressed across cichlid species with divergent tooth and taste bud density, and were expressed in the development of both organs in mice. Synexpression analysis and chemical manipulation of Wnt, BMP, and Hedgehog (Hh) pathways suggest that a common cichlid oral lamina is competent to form teeth or taste buds. Wnt signaling couples tooth and taste bud density and BMP and Hh mediate distinct organ identity. Synthesizing data from fish and mouse, we suggest that the Wnt-BMP-Hh regulatory hierarchy that configures teeth and taste buds on mammalian jaws and tongues may be an evolutionary remnant inherited from ancestors wherein these organs were copatterned from common epithelium.

Soft Bodies, Hard Jaws: An Introduction to the Symposium, with Rotifers as Models of Jaw Diversity

Jaws have evolved numerous times in the animal kingdom and they display a wide variety of structural, compositional, and functional characteristics that reflect their polyphyletic origins. Among soft-bodied invertebrates, jaws are known from annelids, chaetognaths, flatworms, gnathostomulids, micrognathozoans, mollusks, rotifers, and several ecdysozoans. Depending on the taxon, jaws may function in the capture of prey (e.g., chaetognaths and flatworms), processing of prey (e.g., gnathostomulids and onychophorans), or both (e.g., rotifers). Although structural diversity among invertebrates’ jaws is becoming better characterized with the use of electron microscopy, many details remain poorly described, including neuromuscular control, elemental composition, and physical characteristics, such as hardness and resistance to wear. Unfortunately, absence of relevant data has impeded understanding of their functional diversity and evolutionary origins. With this symposium, we bring together researchers of disparately jawed taxa to draw structural and mechanistic comparisons among species to determine their commonalities. Additionally, we show that rotifers’ jaws, which are perhaps the best-characterized jaws among invertebrates, are still enigmatic with regard to their origins and mechanics. Nevertheless, technologies such as energy dispersive X-ray spectroscopy (EDX) and 3D modeling are being used to characterize their chemical composition and to develop physical models that allow exploration of their mechanical properties, respectively. We predict that these methods can also be used to develop biomimetic and bioinspired constructs based on the full range of the complexity of jaws, and that such constructs also can be developed from other invertebrate taxa. These approaches may also shed light on common developmental and physiological processes that facilitate the evolution of invertebrates’ jaws.

Nested Levels of Adaptive Divergence: The Genetic Basis of Craniofacial Divergence and Ecological Sexual Dimorphism

Exemplary systems for adaptive divergence are often characterized by their large degrees of phenotypic variation. This variation represents the outcome of generations of diversifying selection. However, adaptive radiations can also contain a hierarchy of differentiation nested within them where species display only subtle phenotypic differences that still have substantial effects on ecology, function, and ultimately fitness. Sexual dimorphisms are also common in species displaying adaptive divergence and can be the result of differential selection between sexes that produce ecological differences between sexes. Understanding the genetic basis of subtle variation (between certain species or sexes) is therefore important for understanding the process of adaptive divergence. Using cichlids from the dramatic adaptive radiation of Lake Malawi, we focus on understanding the genetic basis of two aspects of relatively subtle phenotypic variation. This included a morphometric comparison of the patterns of craniofacial divergence between two ecologically similar species in relation to the larger adaptive radiation of Malawi, and male-female morphological divergence between their F2 hybrids. We then genetically map craniofacial traits within the context of sex and locate several regions of the genome that contribute to variation in craniofacial shape that is relevant to sexual dimorphism within species and subtle divergence between closely related species, and possibly to craniofacial divergence in the Malawi radiation as a whole. To enhance our search for candidate genes we take advantage of population genomic data and a genetic map that is anchored to the cichlid genome to determine which genes within our QTL regions are associated with SNPs that are alternatively fixed between species. This study provides a holistic understanding of the genetic underpinnings of adaptive divergence in craniofacial shape.

The life history of whole-organism performance

For almost 40 years, studies of whole-organism performance have formed a cornerstone of evolutionary physiology. Although its utility as a heuristic guide is beyond question, and we have learned much about morphological evolution from its application, the ecomorphological paradigm has frequently been applied to performance evolution in ways that range from unsatisfactory to inappropriate. More importantly, the standard ecomorphological paradigm does not account for tradeoffs among performance and other traits, nor between performance traits that are mediated by resource allocation. A revised paradigm that includes such tradeoffs, and the possible ways that performance and fitness-enhancing traits might affect each other, could potentially revivify the study of phenotypic evolution and make important inroads into understanding the relationships between morphology and performance and between performance and Darwinian fitness. We describe such a paradigm, and discuss the various ways that performance and key life-history traits might interact with and affect each other. We emphasize both the proximate mechanisms potentially linking such traits, and the likely ultimate factors driving those linkages, as well as the evolutionary implications for the overall, multivariate phenotype. Finally, we highlight several research directions that will shed light on the evolution and ecology of whole-organism performance and related life-history traits.

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

Background

Adaptation to a new environment can be facilitated by co-inheritance of a suite of phenotypes that are all advantageous in the new habitat. Although experimental evidence demonstrates that multiple phenotypes often map to the same genomic regions, it is challenging to determine whether phenotypes are associated due to pleiotropic effects of a single gene or to multiple tightly linked genes. In the threespine stickleback fish (Gasterosteus aculeatus), multiple phenotypes are associated with a genomic region around the gene Ectodysplasin (Eda), but only the presence of bony lateral plates has been conclusively shown to be caused by Eda.

Results

Here, we ask whether pleiotropy or linkage is responsible for the association between lateral plates and the distribution of the neuromasts of the lateral line. We first identify a strong correlation between plate appearance and changes in the spatial distribution of neuromasts through development. We then use an Eda transgene to induce the formation of ectopic plates in low plated fish, which also results in alterations to neuromast distribution. Our results also show that other loci may modify the effects of Eda on plate formation and neuromast distribution.

Conclusions

Together, these results demonstrate that Eda has pleiotropic effects on at least two phenotypes, highlighting the role of pleiotropy in the genetic basis of adaptation.

Complex craniofacial changes in blind cave-dwelling fish are mediated by genetically symmetric and asymmetric loci

The genetic regulators of regressive craniofacial morphologies are poorly understood. To shed light on this problem, we examined the freshwater fish Astyanax mexicanus, a species with surface-dwelling and multiple independent eyeless cave-dwelling forms. Changes affecting the skull in cavefish include morphological alterations to the intramembranous circumorbital bones encircling the eye. Many of these modifications, however, have evolved separately from eye loss, such as fragmentation of the third suborbital bone. To understand the genetic architecture of these eye-independent craniofacial alterations, we developed and scored 33 phenotypes in the context of an F2 hybrid mapping pedigree bred from Pachón cavefish and surface fish. We discovered several individuals exhibiting dramatic left-right differences in bone formation, such as extensive fragmentation on the right side only. This observation, along with well-known eye size asymmetry in natural cave-dwelling animals, led us to further evaluate left-right genetic differences for the craniofacial complex. We discovered three phenotypes, inclusive of bone fragmentation and fusion, which demonstrated a directional heritable basis only on one side. Interestingly, the overall areas of affected bones were genetically symmetric. Phenotypic effect plots of these novel craniofacial QTL revealed that cave alleles are associated with abnormal conditions such as bony fusion and fragmentation. Moreover, many linked loci overlapped with other cave-associated traits, suggesting regressive craniofacial changes may evolve through linkage or as antagonistic pleiotropic consequences of cave-associated adaptations. These novel findings illuminate significant craniofacial changes accompanying evolution in complete darkness and reveal complex changes to the skull differentially influenced by genetic changes affecting the left and right sides.

Construction of chromosome markers from the Lake Victoria cichlid Paralabidochromis chilotes and their application to comparative mapping

Cichlid fishes in the African Great Lakes are known as a spectacular example of adaptive radiation in vertebrates. Four linkage maps have been constructed to identify the genes responsible for adaptation and speciation, and the genetic linkages of those genes are assumed to play an important role during adaptive evolution. However, it is difficult to analyze such linkages because the linkage groups of one species do not match well with those of the other species. Chromosome markers are a powerful tool for the direct identification of linkage homology between different species. We used information about the linkage map of the Lake Malawi cichlid (Labeotropheus fuelleborni/Metriaclima zebra) to isolate bacterial artificial chromosome (BAC) clones from the BAC library of Paralabidochromis chilotes, Lake Victoria. We identified 18 of 22 P. chilotes chromosomes by single- and multi-color BAC fluorescence in situ hybridization using 19 BAC clones. Comparative mapping with the chromosome markers of P. chilotes in Astatotilapia burtoni (2n = 40) from Lake Tanganyika revealed the chromosome rearrangements that have occurred in this lineage. These chromosome markers will be useful for delineating the process of genome and chromosome evolution in African species.

Introgressive hybridization in a trophically polymorphic cichlid

Trophically polymorphic species could represent lineages that are rapidly diverging along an ecological axis or could phenotypically mark the collapse of species through introgressive hybridization. We investigated patterns of introgression between the trophically polymorphic cichlid fish Herichthys minckleyi and its relative H. cyanoguttatus using a combination of population genetics and species tree analyses. We first examined the distribution of mitochondrial haplotypes within the alternative H. minckleyi pharyngeal jaw morphotypes that are endemic to the small desert valley of Cuatro Ciénegas. We recovered two clusters of mitochondrial haplotypes. The first contained a number of slightly differentiated cytochrome b (cytb) haplotypes that showed some phylogeographic signal and were present in both jaw morphotypes. The other haplotype was monomorphic, highly differentiated from the other cluster, present in equal frequencies in the morphotypes, and identical to H. cyanoguttatus haplotypes found outside Cuatro Ciénegas. Then, we investigated whether H. minckleyi individuals with the H. cyanoguttatus cytb were more evolutionarily similar to H. cyanoguttatus or other H. minckleyi using a species tree analysis of 84 nuclear loci. Both H. minckleyi pharyngeal morphotypes, regardless of their cytb haplotype, were quite distinct from H. cyanoguttatus. However, hybridization could be blurring subdivision within H. minckleyi as the alternative jaw morphotypes were not genetically distinct from one another. Accounting for introgression from H. cyanoguttatus will be essential to understand the evolution of the trophically polymorphic cichlid H. minckleyi.

Lake Malawi cichlid evolution along a benthic/limnetic axis

Divergence along a benthic to limnetic habitat axis is ubiquitous in aquatic systems. However, this type of habitat divergence has largely been examined in low diversity, high latitude lake systems. In this study, we examined the importance of benthic and limnetic divergence within the incredibly species-rich radiation of Lake Malawi cichlid fishes. Using novel phylogenetic reconstructions, we provided a series of hypotheses regarding the evolutionary relationships among 24 benthic and limnetic species that suggests divergence along this axis has occurred multiple times within Lake Malawi cichlids. Because pectoral fin morphology is often associated with divergence along this habitat axis in other fish groups, we investigated divergence in pectoral fin muscles in these benthic and limnetic cichlid species. We showed that the eight pectoral fin muscles and fin area generally tended to evolve in a tightly correlated manner in the Lake Malawi cichlids. Additionally, we found that larger pectoral fin muscles are strongly associated with the independent evolution of the benthic feeding habit across this group of fish. Evolutionary specialization along a benthic/limnetic axis has occurred multiple times within this tropical lake radiation and has produced repeated convergent matching between exploitation of water column habitats and locomotory morphology.

A hybrid genetic linkage map of two ecologically and morphologically divergent Midas cichlid fishes (Amphilophus spp.) obtained by massively parallel DNA sequencing (ddRADSeq)

Cichlid fishes are an excellent model system for studying speciation and the formation of adaptive radiations because of their tremendous species richness and astonishing phenotypic diversity. Most research has focused on African rift lake fishes, although Neotropical cichlid species display much variability as well. Almost one dozen species of the Midas cichlid species complex (Amphilophus spp.) have been described so far and have formed repeated adaptive radiations in several Nicaraguan crater lakes. Here we apply double-digest restriction-site associated DNA sequencing to obtain a high-density linkage map of an interspecific cross between the benthic Amphilophus astorquii and the limnetic Amphilophus zaliosus, which are sympatric species endemic to Crater Lake Apoyo, Nicaragua. A total of 755 RAD markers were genotyped in 343 F₂ hybrids. The map resolved 25 linkage groups and spans a total distance of 1427 cM with an average marker spacing distance of 1.95 cM, almost matching the total number of chromosomes (n = 24) in these species. Regions of segregation distortion were identified in five linkage groups. Based on the pedigree of parents to F₂ offspring, we calculated a genome-wide mutation rate of 6.6 × 10⁻⁸ mutations per nucleotide per generation. This genetic map will facilitate the mapping of ecomorphologically relevant adaptive traits in the repeated phenotypes that evolved within the Midas cichlid lineage and, as the first linkage map of a Neotropical cichlid, facilitate comparative genomic analyses between African cichlids, Neotropical cichlids and other teleost fishes.

Evolution of cichlid vision via trans-regulatory divergence

BACKGROUND

Phenotypic evolution may occur through mutations that affect either the structure or expression of protein-coding genes. Although the evolution of color vision has historically been attributed to structural mutations within the opsin genes, recent research has shown that opsin regulatory mutations can also tune photoreceptor sensitivity and color vision. Visual sensitivity in African cichlid fishes varies as a result of the differential expression of seven opsin genes. We crossed cichlid species that express different opsin gene sets and scanned their genome for expression Quantitative Trait Loci (eQTL) responsible for these differences. Our results shed light on the role that different structural, cis-, and trans-regulatory mutations play in the evolution of color vision.

RESULTS

We identified 11 eQTL that contribute to the divergent expression of five opsin genes. On three linkage groups, several eQTL formed regulatory "hotspots" associated with the expression of multiple opsins. Importantly, however, the majority of the eQTL we identified (8/11 or 73%) occur on linkage groups located trans to the opsin genes, suggesting that cichlid color vision has evolved primarily via trans-regulatory divergence. By modeling the impact of just two of these trans-regulatory eQTL, we show that opsin regulatory mutations can alter cichlid photoreceptor sensitivity and color vision at least as much as opsin structural mutations can.

CONCLUSIONS

Combined with previous work, we demonstrate that the evolution of cichlid color vision results from the interplay of structural, cis-, and especially trans-regulatory loci. Although there are numerous examples of structural and cis-regulatory mutations that contribute to phenotypic evolution, our results suggest that trans-regulatory mutations could contribute to phenotypic divergence more commonly than previously expected, especially in systems like color vision, where compensatory changes in the expression of multiple genes are required in order to produce functional phenotypes.

Genetic interactions controlling sex and color establish the potential for sexual conflict in Lake Malawi cichlid fishes

Sex-determining systems may evolve rapidly and contribute to lineage diversification. In fact, recent work has suggested an integral role of sex chromosome evolution in models of speciation. We use quantitative trait loci analysis of restriction site-associated DNA -tag single nucleotide polymorphisms to identify multiple loci responsible for sex determination and reproductively adaptive color phenotypes in Lake Malawi cichlids. We detect a complex epistatic sex system consisting of a major female heterogametic ZW locus on chromosome 5, two separate male heterogametic XY loci on chromosome 7, and two additional interacting loci on chromosomes 3 and 20. Our data support the known chromosomal linkage between orange blotch color and ZW, as well as novel genetic associations between male blue nuptial color and two sex determining regions (an XY and ZW locus). These results provide further empirical evidence for a complex antagonistic sex-color system in this species flock and suggest a possible role for, and effect of, polygenic sex-determining systems in rapid evolutionary diversification.