Constraint and diversification of developmental trajectories in cichlid facial morphologies

From hatching to juvenile: Larval development of Vieja fenestrata (Teleostei: Cichlidae)

This study delves into the early development of Vieja fenestrata (Cichlidae), with a specific focus on the description of external morphological and morphometric changes, and growth patterns from hatching to the loss of larval characters under controlled laboratory conditions at a temperature of 28°C. Asynchronous hatching was observed between 58 and 60 h postfertilization, with the posterior body emerging first. Over 14 days, significant morphological, physiological, and behavioral changes were observed, revealing a complex developmental trajectory. The initial developmental phases were characterized by rapid vascularization, fin differentiation, and heightened activity, and the subsequent days witnessed the flexion of the notochord, emergence of swim bladder functionality, and transition to exogenous feeding. Maturation progressed with the absorption of the yolk sac, regression of cement glands, and fin ray development, culminating in metamorphosis by 14 days post-hatching. Throughout this period, evolving pigmentation patterns and structural adaptations highlight the species' adaptive strategies. During the larval period of V. fenestrata, substantial changes in morphological proportions were observed. Before the inflection, tail length, trunk length, and body depth had negative allometric growth, and head length, eye diameter, and snout length had positive allometric growth. After the inflection, body depth and snout length showed positive allometric growth; head length and trunk length exhibited isometric growth, whereas tail length and eye diameter demonstrated negative allometric growth. These findings contribute insights into the intricate developmental dynamics of V. fenestrata. Moreover, further research may explore these developmental dynamics' ecological and evolutionary implications.

Heterochronic shift in gene expression leads to ontogenetic morphological divergence between two closely related polyploid species

Heterochrony-alteration to the rate or timing of development-is an important mechanism of trait differentiation associated with speciation. Heterochrony may explain the morphological divergence between two polyploid species, June sucker (Chasmistes liorus) and Utah sucker (Catostomus ardens). The larvae of both species have terminal mouths; however, as adults, June sucker and Utah sucker develop subterminal and ventral mouths, respectively. We document a difference in the timing of shape development and a corresponding change in the timing of gene expression, suggesting the distinctive mouth morphology in June suckers may result from paedomorphosis. Specifically, adult June suckers exhibit an intermediate mouth morphology between the larval (terminal) and ancestral (ventral) states. Endemic and sympatric Chasmistes/Catostomus pairs in two other lakes also are morphologically divergent, but genetically similar. These species pairs could have resulted from the differential expression of genes and corresponding divergence in trait development. Paedomorphosis may lead to adaptive diversification in Catostomids.

Epigenetics and the evolution of form: Experimental manipulation of a chromatin modification causes species-specific changes to the craniofacial skeleton

A central question in biology is the molecular origins of phenotypic diversity. While genetic changes are key to the genotype-phenotype relationship, alterations to chromatin structure and the physical packaging of histone proteins may also be important drivers of vertebrate divergence. We investigate the impact of such an epigenetic mechanism, histone acetylation, within a textbook example of an adaptive radiation. Cichlids of Lake Malawi have adapted diverse craniofacial structures, and here we investigate how histone acetylation influences morphological variation in these fishes. Specifically, we assessed the effect of inhibiting histone deacetylation using the drug trichostatin A (TSA) on developing facial structures. We examined this during three critical developmental windows in two cichlid species with alternate adult morphologies. Exposure to TSA during neural crest cell (NCC) migration and as postmigratory NCCs proliferate in the pharyngeal arches resulted in significant changes in lateral and ventral shape in Maylandia, but not in Tropheops. This included an overall shortening of the head, widening of the lower jaw, and steeper craniofacial profile, all of which are paedomorphic morphologies. In contrast, treatment with TSA during early chondrogenesis did not result in significant morphological changes in either species. Together, these data suggest a sensitivity to epigenetic alterations that are both time- and species-dependent. We find that morphologies are due to nonautonomous or potentially indirect effects on NCC development, including in part a global developmental delay. Our research bolsters the understanding that proper histone acetylation is essential for early craniofacial development and identifies a species-specific robustness to developmental change. Overall, this study demonstrates how epigenetic regulation may play an important role in both generating and buffering morphological variation.

Evo-Devo Mechanobiology: The Missing Link

While the modern framework of evolutionary development (evo-devo) has been decidedly genetic, historic analyses have also considered the importance of mechanics in the evolution of form. With the aid of recent technological advancements in both quantifying and perturbing changes in the molecular and mechanical effectors of organismal shape, how molecular and genetic cues regulate the biophysical aspects of morphogenesis is becoming increasingly well studied. As a result, this is an opportune time to consider how the tissue-scale mechanics that underlie morphogenesis are acted upon through evolution to establish morphological diversity. Such a focus will enable a field of evo-devo mechanobiology that will serve to better elucidate the opaque relations between genes and forms by articulating intermediary physical mechanisms. Here, we review how the evolution of shape is measured and related to genetics, how recent strides have been made in the dissection of developmental tissue mechanics, and how we expect these areas to coalesce in evo-devo studies in the future.

Genes, Morphology, Performance, and Fitness: Quantifying Organismal Performance to Understand Adaptive Evolution

To understand the complexities of morphological evolution, we must understand the relationships between genes, morphology, performance, and fitness in complex traits. Genomicists have made tremendous progress in finding the genetic basis of many phenotypes, including a myriad of morphological characters. Similarly, field biologists have greatly advanced our understanding of the relationship between performance and fitness in natural populations. However, the connection from morphology to performance has primarily been studied at the interspecific level, meaning that in most cases we lack a mechanistic understanding of how evolutionarily relevant variation among individuals affects organismal performance. Therefore, functional morphologists need methods that will allow for the analysis of fine-grained intraspecific variation in order to close the path from genes to fitness. We suggest three methodological areas that we believe are well suited for this research program and provide examples of how each can be applied within fish model systems to build our understanding of microevolutionary processes. Specifically, we believe that structural equation modeling, biological robotics, and simultaneous multi-modal functional data acquisition will open up fruitful collaborations among biomechanists, evolutionary biologists, and field biologists. It is only through the combined efforts of all three fields that we will understand the connection between evolution (acting at the level of genes) and natural selection (acting on fitness).

Zebrafish crocc2 mutants exhibit divergent craniofacial shape, misregulated variability, and aberrant cartilage morphogenesis

BACKGROUND

Phenotypic variation is of paramount importance in development, evolution, and human health; however, the molecular mechanisms that influence organ shape and shape variability are not well understood. During craniofacial development, the behavior of skeletal precursors is regulated by both biochemical and environmental inputs, and the primary cilia play critical roles in transducing both types of signals. Here, we examine a gene that encodes a key constituent of the ciliary rootlets, crocc2, and its role in cartilage morphogenesis in larval zebrafish.

RESULTS

Geometric morphometric analysis of crocc2 mutants revealed altered craniofacial shapes and expanded variation. At the cellular level, we observed altered chondrocyte shapes and planar cell polarity across multiple stages in crocc2 mutants. Notably, cellular defects were specific to areas that experience direct mechanical input. Cartilage cell number, apoptosis, and bone patterning were not affected in crocc2 mutants.

CONCLUSIONS

Whereas "regulatory" genes are widely implicated in patterning the craniofacial skeleton, genes that encode "structural" aspects of the cell are increasingly implicated in shaping the face. Our results add crocc2 to this list, and demonstrate that it affects craniofacial geometry and canalizes phenotypic variation. We propose that it does so via mechanosensing, possibly through the ciliary rootlet. If true, this would implicate a new organelle in skeletal development and evolution.

Morphological and temporal variation in early embryogenesis contributes to species divergence in Malawi cichlid fishes

The cichlid fishes comprise the largest extant vertebrate family and are the quintessential example of rapid "explosive" adaptive radiations and phenotypic diversification. Despite low genetic divergence, East African cichlids harbor a spectacular intra- and interspecific morphological diversity, including the hyper-variable, neural crest (NC)-derived traits such as coloration and craniofacial skeleton. Although the genetic and developmental basis of these phenotypes has been investigated, understanding of when, and specifically how early, in ontogeny species-specific differences emerge, remains limited. Since adult traits often originate during embryonic development, the processes of embryogenesis could serve as a potential source of species-specific variation. Consequently, we designed a staging system by which we compare the features of embryogenesis between three Malawi cichlid species-Astatotilapia calliptera, Tropheops sp. 'mauve' and Rhamphochromis sp. "chilingali"-representing a wide spectrum of variation in pigmentation and craniofacial morphologies. Our results showed fundamental differences in multiple aspects of embryogenesis that could underlie interspecific divergence in adult adaptive traits. First, we identified variation in the somite number and signatures of temporal variation, or heterochrony, in the rates of somite formation. The heterochrony was also evident within and between species throughout ontogeny, up to the juvenile stages. Finally, the identified interspecific differences in the development of pigmentation and craniofacial cartilages, present at the earliest stages of their overt formation, provide compelling evidence that the species-specific trajectories begin divergence during early embryogenesis, potentially during somitogenesis and NC development. Altogether, our results expand our understanding of fundamental cichlid biology and provide new insights into the developmental origins of vertebrate morphological diversity.

East African cichlid fishes

Cichlid fishes are a very diverse and species-rich family of teleost fishes that inhabit lakes and rivers of India, Africa, and South and Central America. Research has largely focused on East African cichlids of the Rift Lakes Tanganyika, Malawi, and Victoria that constitute the biodiversity hotspots of cichlid fishes. Here, we give an overview of the study system, research questions, and methodologies. Research on cichlid fishes spans many disciplines including ecology, evolution, physiology, genetics, development, and behavioral biology. In this review, we focus on a range of organismal traits, including coloration phenotypes, trophic adaptations, appendages like fins and scales, sensory systems, sex, brains, and behaviors. Moreover, we discuss studies on cichlid phylogenies, plasticity, and general evolutionary patterns, ranging from convergence to speciation rates and the proximate and ultimate mechanisms underlying these processes. From a methodological viewpoint, the last decade has brought great advances in cichlid fish research, particularly through the advent of affordable deep sequencing and advances in genetic manipulations. The ability to integrate across traits and research disciplines, ranging from developmental biology to ecology and evolution, makes cichlid fishes a fascinating research system.

The Evolution of Feeding Mechanics in the Danioninae, or Why Giant Danios Don't Suck Like Zebrafish

By linking anatomical structure to mechanical performance we can improve our understanding of how selection shapes morphology. Here we examined the functional morphology of feeding in fishes of the subfamily Danioninae (order Cypriniformes) to determine aspects of cranial evolution connected with their trophic diversification. The Danioninae comprise three major lineages and each employs a different feeding strategy. We gathered data on skull form and function from species in each clade, then assessed their evolutionary dynamics using phylogenetic-comparative methods. Differences between clades are strongly associated with differences in jaw protrusion. The paedomorphic Danionella clade does not use jaw protrusion at all, members of the Danio clade use jaw protrusion for suction production and prey capture, and members of the sister clade to Danio (e.g., Devario and Microdevario) use jaw protrusion to retain prey after capture. The shape of the premaxillary bone is a major determinant of protrusion ability, and premaxilla morphology in each of these lineages is consistent with their protrusion strategies. Premaxilla shapes have evolved rapidly, which indicates that they have been subjected to strong selection. We compared premaxilla development in giant danio (Devario aequipinnatus) and zebrafish (Danio rerio) and discuss a developmental mechanism that could shift danionine fishes between the feeding strategies employed by these species and their respective clades. We also identified a highly integrated evolutionary module that has been an important factor in the evolution of trophic mechanics within the Danioninae.

Morphometric and Genetic Description of Trophic Adaptations in Cichlid Fishes

Since Darwin, biologists have sought to understand the evolution and origins of phenotypic adaptations. The skull is particularly diverse due to intense natural selection on feeding biomechanics. We investigated the genetic and molecular origins of trophic adaptation using Lake Malawi cichlids, which have undergone an exemplary evolutionary radiation. We analyzed morphological differences in the lateral and ventral head shape among an insectivore that eats by suction feeding, an obligate biting herbivore, and their F2 hybrids. We identified variation in a series of morphological traits—including mandible width, mandible length, and buccal length—that directly affect feeding kinematics and function. Using quantitative trait loci (QTL) mapping, we found that many genes of small effects influence these craniofacial adaptations. Intervals for some traits were enriched in genes related to potassium transport and sensory systems, the latter suggesting co-evolution of feeding structures and sensory adaptations for foraging. Despite these indications of co-evolution of structures, morphological traits did not show covariation. Furthermore, phenotypes largely mapped to distinct genetic intervals, suggesting that a common genetic basis does not generate coordinated changes in shape. Together, these suggest that craniofacial traits are mostly inherited as separate modules, which confers a high potential for the evolution of morphological diversity. Though these traits are not restricted by genetic pleiotropy, functional demands of feeding and sensory structures likely introduce constraints on variation. In all, we provide insights into the quantitative genetic basis of trophic adaptation, identify mechanisms that influence the direction of morphological evolution, and provide molecular inroads to craniofacial variation.

Early shape divergence of developmental trajectories in the jaw of galeomorph sharks

Background

The onset of morphological differences between related groups can be tracked at early stages during embryological development. This is expressed in functional traits that start with minor variations, but eventually diverge to defined specific morphologies. Several processes during this period, like proliferation, remodelling, and apoptosis for instance, can account for the variability observed between related groups. Morphological divergence through development is often associated with the hourglass model, in which early stages display higher variability and reach a conserved point with reduced variability from which divergence occurs again to the final phenotype.

Results

Here we explored the patterns of developmental shape changes in the lower jaw of two shark species, the bamboo shark (Chiloscyllium punctatum) and the catshark (Scyliorhinus canicula). These two species present marked differences in their foraging behaviour, which is reflected in their adult jaw morphology. By tracing the developmental sequence of the cartilage condensation, we identified the onset of cartilage for both species at around stage 31. Other structures that developed later without a noticeable anlage were the labial cartilages, which appear at around stage 33. We observed that the lower jaw displays striking differences in shape from the earliest moments, without any overlap in shape through the compared stages.

Conclusions

The differences observed are also reflected in the functional variation in feeding mechanism between both species. Likewise, the trajectory analysis shows that the main differences are in the magnitude of the shape change through time. Both species follow a unique trajectory, which is explained by the timing between stages.

Time to synchronize our clocks: Connecting developmental mechanisms and evolutionary consequences of heterochrony

Heterochrony, defined as a change in the timing of developmental events altering the course of evolution, was first recognized by Ernst Haeckel in 1866. Haeckel's original definition was meant to explain the observed parallels between ontogeny and phylogeny, but the interpretation of his work became a source of controversy over time. Heterochrony took its modern meaning following the now classical work in the 1970-80s by Steven J. Gould, Pere Alberch, and co-workers. Predicted and described heterochronic scenarios emphasize the many ways in which developmental changes can influence evolution. However, while important examples of heterochrony detected with comparative morphological methods have multiplied, the more mechanistic understanding of this phenomenon lagged conspicuously behind. Considering the rapid progress in imaging and molecular tools available now for developmental biologists, this review aims to stress the need to take heterochrony research to the next level. It is time to synchronize the different levels of heterochrony research into a single analysis flow: from studies on organismal-level morphology to cells to molecules and genes, using complementary techniques. To illustrate how to achieve a more comprehensive understanding of phyletic morphological diversification associated with heterochrony, we discuss several recent case studies at various phylogenetic scales that combine morphological, cellular, and molecular analyses. Such a synergistic approach offers to more fully integrate phylogenetic and ontogenetic dimensions of the fascinating evolutionary phenomenon of heterochrony.

Developmental origins of the crocodylian skull table and platyrostral face

The dorsoventrally flattened skull typifies extant Crocodylia perhaps more than any other anatomical feature and is generally considered an adaptation for semi-aquatic feeding. Although the evolutionary origins of caniofacial flattening have been extensively studied, the developmental origins have yet to be explored. To understand how the skull table and platyrostral snout develop, we quantified embryonic development and post-hatching growth (ontogeny) of the crocodylian skull in lateral view using geometric morphometrics. Our dataset (n = 103) includes all but one extant genus and all of the major ecomorphs, including the extremely slender-snouted Gavialis and Tomistoma. Our analysis reveals that the embryonic development of the flattened skull is remarkably similar across ecomorphs, including the presence of a conserved initial embryonic skull shape, similar to prior analysis of dorsal snout shape. Although differences during posthatching ontogeny are recovered among ecomorphs, embryonic patterns are not distinct, revealing an important shift in developmental rate near hatching. In particular, the flattened skull table is achieved by the end of embryonic development with no changes after hatching. Further, the rotation of skull roof and facial bones during development is critical for the stereotypical flatness of the crocodylian skull. Our results suggest selection on hatchling performance and constraints on embryonic skull shape may have been important in this pattern of developmental conservation. The appearance of aspects of cranial flatness among Jurassic stem crocodylians suggests key aspects of these cranial developmental patterns may have been conserved for over 200 million years.

Comparative scale morphology in the adaptive radiation of cichlid fishes (Perciformes: Cichlidae) from Lake Tanganyika

The morphology of fish scales has been investigated for > 200 years, but research on evolutionary patterns of scale morphology is scarce. Here, we study scale morphology and its evolution in the adaptive radiation of cichlid fishes from Lake Tanganyika, which are known for their exceptional diversity in habitat use, feeding ecology and morphology. Based on a geometric morphometric approach on eight scales per specimen (covering different body regions), we quantify scale types and morphology across nearly all ~240 species of the cichlid adaptive radiation in Lake Tanganyika. We first show that scale type, shape and ctenii coverage vary along the body, which is probably attributable to adaptations to different functional demands on the respective scales. Our comparative analyses reveal that flank scale size is tightly linked to phylogeny, whereas scale shape and ctenii coverage can be explained only in part by phylogenetic history and/or our proxy for ecology (stable isotopes and body shape), suggesting an additional adaptive component. We also show that our measured scale characteristics can help to assign an individual scale to a taxonomic group or ecotype. Thus, our data may serve as a valuable resource for taxonomic studies and to interpret fossil finds.

Natural toxic impact and thyroid signalling interplay orchestrates riverine adaptive divergence of salmonid fish

Adaptive radiation in fishes has been actively investigated over the last decades. Along with numerous well-studied cases of lacustrine radiation, some examples of riverine sympatric divergence have been recently discovered. In contrast to the lakes, the riverine conditions do not provide evident stability in the ecological gradients. Consequently, external factors triggering the radiation, as well as developmental mechanisms underpinning it, remain unclear. Herein, we present the comprehensive study of external and internal drivers of the riverine adaptive divergence of the salmonid fish Salvelinus malma. In the Kamchatka River, north-east Asia, this species splits in the reproductively isolated morphs that drastically differ in ecology and morphology: the benthivorous Dolly Varden (DV) and the piscivorous stone charr (SC). To understand why and how these morphs originated, we performed a series of field and experimental work, including common-garden rearing, comparative ontogenetic, physiological and endocrinological analyses, hormonal 'engineering' of phenotypes and acute toxicological tests. We revealed that the type of spawning ground acts as the decisive factor driving the radiation of S. malma. In contrast to DV spawning in the leaf krummholz zone, SC reproduces in the zone of coniferous forest, which litter has a toxic impact on developing fishes. SC enhances resistance to the toxicants via metabolism acceleration provided by the elevated thyroid hormone expenditure. These physiological changes lead to the multiple heterochronies resulting in a specific morphology and ecology of SC. Salvelinus malma represents a notable example of how the thyroid axis contributes to the generation of diverse phenotypic outcomes underlying the riverine sympatric divergence. Our findings, along with the paleoecology data concerning spruce forest distribution during the Pleistocene, provide an opportunity to reconstruct a scenario of S. malma divergence. Taken together, obtained results with the data of the role of thyroid hormones in the ontogeny and diversification of fishes contribute a resource to consider the thyroid axis as a prime director orchestrating the phenotypic plasticity promoting evolutionary diversification under the changing environmental conditions.

Endochondral growth zone pattern and activity in the zebrafish pharyngeal skeleton

BACKGROUND

Endochondral ossification is a major bone forming mechanism in vertebrates, defects in which can result in skeletal dysplasia or craniofacial anomalies in humans. The zebrafish holds great potential to advance our understanding of endochondral growth zone development and genetics, yet several important aspects of its biology remain unexplored. Here we provide a comprehensive description of endochondral growth zones in the pharyngeal skeleton, including their developmental progression, cellular activity, and adult fates.

RESULTS

Postembryonic growth of the pharyngeal skeleton is supported by endochondral growth zones located either at skeletal epiphyses or synchondroses. Col2a1a and col10a1a in situ hybridization and anti-PCNA immunostaining identify resting-, hypertrophic- and proliferative zones, respectively, in pharyngeal synchondroses. Cellular hypertrophy and matrix deposition contribute little, if at all, to axial growth in most skeletal elements. Zebrafish endochondral growth zones develop during metamorphosis and arrest in adults.

CONCLUSIONS

Two endochondral growth zone configurations in the zebrafish pharyngeal skeleton produce either unidirectional (epiphyses) or bidirectional (synchondroses) growth. Cell proliferation drives endochondral growth and its modulation, in contrast to mammalian long bones in which bone length depends more on cell enlargement during hypertrophy and intramembranous ossification is the default mechanism of bone growth in zebrafish adults.

The Genomic Substrate for Adaptive Radiation: Copy Number Variation across 12 Tribes of African Cichlid Species

The initial sequencing of five cichlid genomes revealed an accumulation of genetic variation, including extensive copy number variation in cichlid lineages particularly those that have undergone dramatic evolutionary radiation. Gene duplication has the potential to generate substantial molecular substrate for the origin of evolutionary novelty. We use array-based comparative heterologous genomic hybridization to identify copy number variation events (CNVEs) for 168 samples representing 53 cichlid species including the 5 species for which full genome sequence is available. We identify an average of 50-100 CNVEs per individual. For those species represented by multiple samples, we identify 150-200 total CNVEs suggesting a substantial amount of intraspecific variation. For these species, only ∼10% of the detected CNVEs are fixed. Hierarchical clustering of species according to CNVE data recapitulates phylogenetic relationships fairly well at both the tribe and radiation level. Although CNVEs are detected on all linkage groups, they tend to cluster in "hotspots" and are likely to contain and be flanked by transposable elements. Furthermore, we show that CNVEs impact functional categories of genes with potential roles in adaptive phenotypes that could reasonably promote divergence and speciation in the cichlid clade. These data contribute to a more complete understanding of the molecular basis for adaptive natural selection, speciation, and evolutionary radiation.

Divergence in larval jaw gene expression reflects differential trophic adaptation in haplochromine cichlids prior to foraging

Background

Understanding how variation in gene expression contributes to morphological diversity is a major goal in evolutionary biology. Cichlid fishes from the East African Great lakes exhibit striking diversity in trophic adaptations predicated on the functional modularity of their two sets of jaws (oral and pharyngeal). However, the transcriptional basis of this modularity is not so well understood, as no studies thus far have directly compared the expression of genes in the oral and pharyngeal jaws. Nor is it well understood how gene expression may have contributed to the parallel evolution of trophic morphologies across the replicate cichlid adaptive radiations in Lake Tanganyika, Malawi and Victoria.

Results

We set out to investigate the role of gene expression divergence in cichlid fishes from these three lakes adapted to herbivorous and carnivorous trophic niches. We focused on the development stage prior to the onset of exogenous feeding that is critical for understanding patterns of gene expression after oral and pharyngeal jaw skeletogenesis, anticipating environmental cues. This framework permitted us for the first time to test for signatures of gene expression underlying jaw modularity in convergent eco-morphologies across three independent adaptive radiations. We validated a set of reference genes, with stable expression between the two jaw types and across species, which can be important for future studies of gene expression in cichlid jaws. Next we found evidence of modular and non-modular gene expression between the two jaws, across different trophic niches and lakes. For instance, prdm1a, a skeletogenic gene with modular anterior-posterior expression, displayed higher pharyngeal jaw expression and modular expression pattern only in carnivorous species. Furthermore, we found the expression of genes in cichlids jaws from the youngest Lake Victoria to exhibit low modularity compared to the older lakes.

Conclusion

Overall, our results provide cross-species transcriptional comparisons of modularly-regulated skeletogenic genes in the two jaw types, implicating expression differences which might contribute to the formation of divergent trophic morphologies at the stage of larval independence prior to foraging.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1483-3) contains supplementary material, which is available to authorized users.

Extending the Geometric Approach for Studying Biomechanical Motions

Whether it is swimming, walking, eating, or jumping, motions are a fundamental way in which organisms interact with their environment. Understanding how morphology contributes to motion is a primary focus of kinematic research and is necessary for gaining insights into the evolution of functional systems. However, an element that is largely missing from traditional analyses of motion is the spatial context in which they occur. We explore an application of geometric morphometrics (GM) for analyzing and comparing motions to evaluate the outputs of biomechanical linkage models. We focus on a common model for oral jaw mechanics of perciform fishes, the fourbar linkage, using GM to summarize motion as a trajectory of shape change. Two traits derived from trajectories capture the total kinesis generated by a linkage (trajectory length) and the kinematic asynchrony (KA) of its mobile components (trajectory nonlinearity). Oral jaw fourbar data from two subfamilies of Malagasy cichlids were used to generate form–function landscapes, describing broad features of kinematic diversity. Our results suggest that kinesis and KA have complex relationships with fourbar morphology, each displaying a pattern in which different shapes possess equivalent kinematic trait values, known as many-to-one mapping of form-to-function. Additionally, we highlight the observation that KA captures temporal differences in the activation of motion components, a feature of kinesis that has long been appreciated but was difficult to measure. The methods used here to study fourbar linkages can also be applied to more complex biomechanical models and broadly to motions of live organisms. We suggest that they provide a suitable alternative to traditional approaches for evaluating linkage function and kinematics.

Heterochronic shifts and conserved embryonic shape underlie crocodylian craniofacial disparity and convergence

The distinctive anatomy of the crocodylian skull is intimately linked with dietary ecology, resulting in repeated convergence on blunt- and slender-snouted ecomorphs. These evolutionary shifts depend upon modifications of the developmental processes which direct growth and morphogenesis. Here we examine the evolution of cranial ontogenetic trajectories to shed light on the mechanisms underlying convergent snout evolution. We use geometric morphometrics to quantify skeletogenesis in an evolutionary context and reconstruct ancestral patterns of ontogenetic allometry to understand the developmental drivers of craniofacial diversity within Crocodylia. Our analyses uncovered a conserved embryonic region of morphospace (CER) shared by all non-gavialid crocodylians regardless of their eventual adult ecomorph. This observation suggests the presence of conserved developmental processes during early development (before Ferguson stage 20) across most of Crocodylia. Ancestral state reconstruction of ontogenetic trajectories revealed heterochrony, developmental constraint, and developmental systems drift have all played essential roles in the evolution of ecomorphs. Based on these observations, we conclude that two separate, but interconnected, developmental programmes controlling craniofacial morphogenesis and growth enabled the evolutionary plasticity of skull shape in crocodylians.

How a growing organismal perspective is adding new depth to integrative studies of morphological evolution

Over the past half century, the field of Evolutionary Developmental Biology, or Evo-devo, has integrated diverse fields of biology into a more synthetic understanding of morphological diversity. This has resulted in numerous insights into how development can evolve and reciprocally influence morphological evolution, as well as generated several novel theoretical areas. Although comparative by default, there remains a great gap in our understanding of adaptive morphological diversification and how developmental mechanisms influence the shape and pattern of phenotypic variation. Herein we highlight areas of research that are in the process of filling this void, and areas, if investigated more fully, that will add new insights into the diversification of morphology. At the centre of our discussion is an explicit awareness of organismal biology. Here we discuss an organismal framework that is supported by three distinct pillars. First, there is a need for Evo-devo to adopt a high-resolution phylogenetic approach in the study of morphological variation and its developmental underpinnings. Secondly, we propose that to understand the dynamic nature of morphological evolution, investigators need to give more explicit attention to the processes that generate evolutionarily relevant variation at the population level. Finally, we emphasize the need to address more thoroughly the processes that structure variation at micro- and macroevolutionary scales including modularity, morphological integration, constraint, and plasticity. We illustrate the power of these three pillars using numerous examples from both invertebrates and vertebrates to emphasize that many of these approaches are already present within the field, but have yet to be formally integrated into many research programs. We feel that the most exciting new insights will come where the traditional experimental approaches to Evo-devo are integrated more thoroughly with the principles of this organismal framework.

Changes in ontogenetic patterns facilitate diversification in skull shape of Australian agamid lizards

BACKGROUND

Morphological diversity among closely related animals can be the result of differing growth patterns. The Australian radiation of agamid lizards (Amphibolurinae) exhibits great ecological and morphological diversity, which they have achieved on a continent-wide scale, in a relatively short period of time (30 million years). Amphibolurines therefore make an ideal study group for examining ontogenetic allometry. We used two-dimensional landmark-based geometric morphometric methods to characterise the postnatal growth patterns in cranial shape of 18 species of amphibolurine lizards and investigate the associations between cranial morphology, and life habit and phylogeny.

RESULTS

For most amphibolurine species, juveniles share a similar cranial phenotype, but by adulthood crania are more disparate in shape and occupy different sub-spaces of the total shape space. To achieve this disparity, crania do not follow a common post-natal growth pattern; there are differences among species in both the direction and magnitude of change in morphospace. We found that these growth patterns among the amphibolurines are significantly associated with ecological life habits. The clade Ctenophorus includes species that undergo small magnitudes of shape change during growth. They have dorsoventrally deep, blunt-snouted skulls (associated with terrestrial lifestyles), and also dorsoventrally shallow skulls (associated with saxicolous lifestyles). The sister clade to Ctenophorus, which includes the bearded dragon (Pogona), frill-neck lizard (Chlamydosaurus), and long-nosed dragon (Gowidon), exhibit broad and robust post-orbital regions and differing snout lengths (mainly associated with scansorial lifestyles).

CONCLUSIONS

Australian agamids show great variability in the timing of development and divergence of growth trajectories which results in a diversity of adult cranial shapes. Phylogenetic signal in cranial morphology appears to be largely overwritten by signals that reflect life habit. This knowledge about growth patterns and skull shape diversity in agamid lizards will be valuable for placing phylogenetic, functional and ecological studies in a morphological context.

Muscle-induced loading as an important source of variation in craniofacial skeletal shape

The shape of the craniofacial skeleton is constantly changing through ontogeny and reflects a balance between developmental patterning and mechanical-load-induced remodeling. Muscles are a major contributor to producing the mechanical environment that is crucial for "normal" skull development. Here, we use an F₅ hybrid population of Lake Malawi cichlids to characterize the strength and types of associations between craniofacial bones and muscles. We focus on four bones/bone complexes, with different developmental origins, alongside four muscles with distinct functions. We used micro-computed tomography to extract 3D information on bones and muscles. 3D geometric morphometrics and volumetric measurements were used to characterize bone and muscle shape, respectively. Linear regressions were performed to test for associations between bone shape and muscle volume. We identified three types of associations between muscles and bones: weak, strong direct (i.e., muscles insert directly onto bone), and strong indirect (i.e., bone is influenced by muscles without a direct connection). In addition, we show that although the shape of some bones is relatively robust to muscle-induced mechanical stimulus, others appear to be highly sensitive to muscular input. Our results imply that the roles for muscular input on skeletal shape extend beyond specific points of origin or insertion and hold significant potential to influence broader patterns of craniofacial geometry. Thus, changes in the loading environment, either as a normal course of ontogeny or if an organism is exposed to a novel environment, may have pronounced effects on skeletal shape via near and far-ranging effects of muscular loading.

Modeling craniofacial development reveals spatiotemporal constraints on robust patterning of the mandibular arch

How does pattern formation occur accurately when confronted with tissue growth and stochastic fluctuations (noise) in gene expression? Dorso-ventral (D-V) patterning of the mandibular arch specifies upper versus lower jaw skeletal elements through a combination of Bone morphogenetic protein (Bmp), Endothelin-1 (Edn1), and Notch signaling, and this system is highly robust. We combine NanoString experiments of early D-V gene expression with live imaging of arch development in zebrafish to construct a computational model of the D-V mandibular patterning network. The model recapitulates published genetic perturbations in arch development. Patterning is most sensitive to changes in Bmp signaling, and the temporal order of gene expression modulates the response of the patterning network to noise. Thus, our integrated systems biology approach reveals non-intuitive features of the complex signaling system crucial for craniofacial development, including novel insights into roles of gene expression timing and stochasticity in signaling and gene regulation.

Experimental evidence of the role of heterochrony in evolution of the Mesoamerican cichlids pigment patterns

The Mesoamerican cichlids display a spectacular diversity of pigment patterns, which serve a variety of functions and serve as a strong selective trait for this lineage. The development and variation of coloration in the Mesoamerican cichlids have been detailed by several groups. In particular, Říčan, Musilová, Muška, and Novák () and Říčan, Piálek, Dragová, and Novák () determined homology of pattern and revealed four alternative types of coloration and their ontogeny. In this work, this group posed an "ontogenetic timing hypothesis" proposing heterochronic shifts underlying major transitions in the evolution of the Mesoamerican cichlids. Here, we experimentally test this hypothesis by experimentally altering timing of pigment pattern formation in the convict cichlid Amatitlania nigrofasciata, a member of the Mesoamerican cichlids, via manipulations of thyroid hormone (TH) function. The response of different pigment cell lineages to TH-perturbations revealed that the transition from larval to juvenile coloration in the convict cichlid is under the control of TH-signaling. Importantly, hormonally induced changes in the timing of pigment cell lineages' development resulted in shifts of coloration ontogeny type observed between lineages and led to the appearance of phenotypes mimicking those in phylogenetically close and distant species. Thus, our findings support the hypothesis that simple changes in ontogenetic timing underlies species specific patterns in pigmentation and provide new perspectives for studying the role of endocrine signaling in the evolution of cichlids.

Feeding ecology underlies the evolution of cichlid jaw mobility

The fish feeding apparatus is among the most diverse functional systems in vertebrates. While morphological and mechanical variations of feeding systems are well studied, we know far less about the diversity of the motions that they produce. We explored patterns of feeding movements in African cichlids from Lakes Malawi and Tanganyika, asking whether the degree of kinesis is associated with dietary habits of species. We used geometric morphometrics to measure feeding kinesis as trajectories of shape change, based on 326 high-speed videos in 56 species. Cranial morphology was significantly related to feeding movements, both of which were distributed along a dietary axis associated with prey evasiveness. Small-mouthed cichlids that feed by scraping algae and detritus from rocks had low kinesis strikes, while large-mouthed species that eat large, evasive prey (fishes and shrimps) generated the greatest kinesis. Despite having higher overall kinesis, comparisons of trajectory shape (linearity) revealed that cichlids that eat mobile prey also displayed more kinematically conserved, or efficient, feeding motions. Our work indicates that prey evasiveness is strongly related to the evolution of cichlid jaw mobility, suggesting that this same relationship may explain the origins and diversity of highly kinetic jaws that characterize the super-radiation of spiny-rayed fishes.

Genetic analyses in Lake Malawi cichlids identify new roles for Fgf signaling in scale shape variation

Elasmoid scales are the most common epithelial appendage among vertebrates, however an understanding of the genetic mechanisms that underlie variation in scale shape is lacking. Using an F2 mapping cross between morphologically distinct cichlid species, we identified >40 QTL for scale shape at different body positions. We show that while certain regions of the genome regulate variation in multiple scales, most are specific to scales at distinct positions. This suggests a degree of regional modularity in scale development. We also identified a single QTL for variation in scale shape disparity across the body. Finally, we screened a QTL hotspot for candidate loci, and identified the Fgf receptor fgfr1b as a prime target. Quantitative rtPCR and small molecule manipulation support a role for Fgf signaling in shaping cichlid scales. While Fgfs have previously been implicated in scale loss, these data reveal new roles for the pathway in scale shape variation.

Juvenile divergence in adaptive traits among seven sympatric fish ecomorphs arises before moving to different lacustrine habitats

Identifying the mechanisms initiating sympatric diversification in vertebrates has remained a conceptual challenge. Here, we analyse an assemblage of sympatric charr (Salvelinus malma) morphs from landlocked Lake Kronotskoe basin as a model to uncover the divergence pathways in freshwater fishes during the early life history stages. All morphs have distinct developmental biology, but a similar developmental rate retardation compared to the ancestor. Our study reveals that adult morphological differences, which acquire functionality at maturation, originate in the early juvenile stages due to heterochrony in skeletogenesis and allometric changes triggered by variation in metabolic activity. The craniofacial differences among the morphs result from asynchronous development of several skeletal modules. The accelerated ossification of teeth-armed bones occurs in predatory feeding morphs, whereas cranial cover ossification is promoted in benthivorous morphs. These contrasting growth patterns have led to seven phenotypes that span a range far beyond the ancestral variability. The most distinct morphs are a riverine spawning, epilimnetic predator and a lacustrine spawning, profundal benthic feeder. Taken together, we argue that the adaptive morphological differentiation in these sympatric freshwater fishes is driven by diverging patterns in ossification rate and metabolic activity against a background of uneven somatic growth. This divergence is primarily associated with basic environmental differences on the nursery grounds that might be unrelated to resource use. This nonheritable phenotype divergence is then exposed to natural selection that could result in further adaptive genetic changes.

The Role of Alternative Splicing and Differential Gene Expression in Cichlid Adaptive Radiation

Species diverge eco-morphologically through the continuous action of natural selection on functionally important structures, producing alternative adaptive morphologies. In cichlid fishes, the oral and pharyngeal jaws are such key structures. Adaptive variation in jaw morphology contributes to trophic specialization, which is hypothesized to fuel their rapid speciation in the East African Great Lakes. Much is known about the genes involved in cichlid jaw and craniofacial development. However, it is still unclear what salient sources of variation gave rise to trophic-niche specialization, facilitating adaptive radiation. Here, we explore two sources of transcriptional variation that may underlie species-specific disparities in jaw morphology. Using whole transcriptome RNA-sequencing, we analyze differences in gene expression and alternative splicing, at the end of postlarval development, in fully functional jaws of six species of cichlids from the Lake Tanganyika tribe Tropheini. Our data reveal a surprisingly high degree of alternative splicing events compared with gene expression differences among species and trophic types. This suggests that differential trophic adaptation of the jaw apparatus may have been shaped by transcriptional rewiring of splicing as well as gene expression variation during the rapid radiation of the Tropheini. Specifically, genes undergoing splicing across most species were found to be enriched for pharyngeal jaw gene ontology terms. Overall, jaw transcriptional patterns at postlarval developmental stage were highly dynamic and species-specific. In conclusion, this work indicates that shifts in alternative splicing could have played a more important role in cichlid adaptive radiation, and possibly adaptive radiation in general, than currently recognized.

Baby fish working out: an epigenetic source of adaptive variation in the cichlid jaw

Understanding the developmental processes that underlie the production of adaptive variation (i.e. the 'arrival of the fittest') is a major goal of evolutionary biology. While most evo-devo studies focus on the genetic underpinnings of adaptive phenotypic variation, factors beyond changes in nucleotide sequence can also play a major role in shaping developmental outcomes. Here, we document a vigorous but enigmatic gaping behaviour during the early development of Lake Malawi cichlid larvae. The onset of the behaviour precedes the formation of bone, and we predicted that it might influence craniofacial shape by affecting the mechanical environment in which bone develops. Consistent with this, we found that both natural variation and experimental manipulation of this behaviour induced differential skeletal development that foreshadows adaptive variation in adult trophic morphology. In fact, the magnitude of difference in skeletal morphology induced by these simple shifts in behaviour was similar to those predicted to be caused by genetic factors. Finally, we demonstrate that this mechanical-load-induced shift in skeletal development is associated with differences in ptch1 expression, a gene previously implicated in mediating between-species differences in skeletal shape. Our results underscore the complexity of development, and the importance of epigenetic (sensu Waddington) mechanisms in determining adaptive phenotypic variation.

The Cyprinodon variegatus genome reveals gene expression changes underlying differences in skull morphology among closely related species

BACKGROUND

Understanding the genetic and developmental origins of phenotypic novelty is central to the study of biological diversity. In this study we identify modifications to the expression of genes at four developmental stages that may underlie jaw morphological differences among three closely related species of pupfish (genus Cyprinodon) from San Salvador Island, Bahamas. Pupfishes on San Salvador Island are trophically differentiated and include two endemic species that have evolved jaw morphologies unlike that of any other species in the genus Cyprinodon.

RESULTS

We find that gene expression differs significantly across recently diverged species of pupfish. Genes such as Bmp4 and calmodulin, previously implicated in jaw diversification in African cichlid fishes and Galapagos finches, were not found to be differentially expressed among species of pupfish. Instead we find multiple growth factors and cytokine/chemokine genes to be differentially expressed among these pupfish taxa. These include both genes and pathways known to affect craniofacial development, such as Wnt signaling, as well as novel genes and pathways not previously implicated in craniofacial development. These data highlight both shared and potentially unique sources of jaw diversity in pupfish and those identified in other evolutionary model systems such as Galapagos finches and African cichlids.

CONCLUSIONS

We identify modifications to the expression of genes involved in Wnt signaling, Igf signaling, and the inflammation response as promising avenues for future research. Our project provides insight into the magnitude of gene expression changes contributing to the evolution of morphological novelties, such as jaw structure, in recently diverged pupfish species.

Cranial asymmetry arises later in the life history of the blind Mexican cavefish, Astyanax mexicanus

As a consequence of adaptation to the cave environment, the blind Mexican cavefish, Astyanax mexicanus, has evolved several cranial aberrations including changes to bone sizes, shapes and presence of numerous lateral asymmetries. Prior studies of cranial asymmetry in cavefish focused strictly on adult specimens. Thus, the extent to which these asymmetries emerge in adulthood, or earlier in the life history of cavefish, was unknown. We performed a geometric morphometric analysis of shape variation in the chondrocranium and osteocranium across life history in two distinct cavefish populations and surface-dwelling fish. The cartilaginous skull in juveniles was bilaterally symmetric and chondrocranial shape was conserved in all three populations. In contrast, bony skull shapes segregated into significantly distinct groups in adults. Cavefish demonstrated significant asymmetry for the bones surrounding the collapsed eye orbit, and the opercle bone posterior to the eye orbit. Interestingly, we discovered that cavefish also exhibit directional “bends” in skull shape, almost always biased to the left. In sum, this work reveals that asymmetric craniofacial aberrations emerge later in the cavefish life history. These abnormalities may mirror asymmetries in the lateral line sensory system, reflect a ‘handedness’ in cavefish swimming behavior, or evolve through neutral processes.

Signalling pathways in trophic skeletal development and morphogenesis: Insights from studies on teleost fish

During the development of the vertebrate feeding apparatus, a variety of complicated cellular and molecular processes participate in the formation and integration of individual skeletal elements. The molecular mechanisms regulating the formation of skeletal primordia and their development into specific morphological structures are tightly controlled by a set of interconnected signalling pathways. Some of these pathways, such as Bmp, Hedgehog, Notch and Wnt, are long known for their pivotal roles in craniofacial skeletogenesis. Studies addressing the functional details of their components and downstream targets, the mechanisms of their interactions with other signals as well as their potential roles in adaptive morphological divergence, are currently attracting considerable attention. An increasing number of signalling pathways that had previously been described in different biological contexts have been shown to be important in the regulation of jaw skeletal development and morphogenesis. In this review, I provide an overview of signalling pathways involved in trophic skeletogenesis emphasizing studies of the most species-rich group of vertebrates, the teleost fish, which through their evolutionary history have undergone repeated episodes of spectacular trophic diversification.

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.

Cichlid fishes as a model to understand normal and clinical craniofacial variation

We have made great strides towards understanding the etiology of craniofacial disorders, especially for 'simple' Mendelian traits. However, the facial skeleton is a complex trait, and the full spectrum of genetic, developmental, and environmental factors that contribute to its final geometry remain unresolved. Forward genetic screens are constrained with respect to complex traits due to the types of genes and alleles commonly identified, developmental pleiotropy, and limited information about the impact of environmental interactions. Here, we discuss how studies in an evolutionary model - African cichlid fishes - can complement traditional approaches to understand the genetic and developmental origins of complex shape. Cichlids exhibit an unparalleled range of natural craniofacial morphologies that model normal human variation, and in certain instances mimic human facial dysmorphologies. Moreover, the evolutionary history and genomic architecture of cichlids make them an ideal system to identify the genetic basis of these phenotypes via quantitative trait loci (QTL) mapping and population genomics. Given the molecular conservation of developmental genes and pathways, insights from cichlids are applicable to human facial variation and disease. We review recent work in this system, which has identified lbh as a novel regulator of neural crest cell migration, determined the Wnt and Hedgehog pathways mediate species-specific bone morphologies, and examined how plastic responses to diet modulate adult facial shapes. These studies have not only revealed new roles for existing pathways in craniofacial development, but have identified new genes and mechanisms involved in shaping the craniofacial skeleton. In all, we suggest that combining work in traditional laboratory and evolutionary models offers significant potential to provide a more complete and comprehensive picture of the myriad factors that are involved in the development of complex traits.

What shapes the oral jaws? Accommodation of complex dentition correlates with premaxillary but not mandibular shape

BACKGROUND

Teeth are integrated into the vertebrate oral jaws to provide a functional unit for feeding, however little is known about how this integration occurs during growth and development. The purpose of this study is to identify the ontogenetic changes in oral jaw shape that are associated with the transition of the oral dentition from unicuspid teeth to multicuspid teeth. Here, we compare the shape of the occluding upper (premaxilla) and lower (mandible) jaws of the toothed Mexican tetra (Astyanax mexicanus) and the toothless (oral teeth present, pharyngeal teeth absent) zebrafish (Danio rerio) over development. Gross morphology combined with morphometric analyses were used to analyse shape changes of the occluding oral jaws in each species. Histological analyses were also used to examine the development of the mandibular symphysis.

RESULTS

The occluding edge of the premaxilla is the first region to ossify in the Mexican tetra, but the last to ossify in zebrafish. Morphometric analyses revealed that the early shape of the premaxillae (in fish younger than 8mm SL) is the same in each species but that the premaxilla shape changes significantly at larger sizes. These changes are apparent in the tooth bearing region of the Mexican tetra. The rostral region of the mandible also houses teeth, however ossification and shape in this region were surprisingly similar between species despite differences in the presence of oral dentition. Furthermore, we found that the mandibular symphysis of the Mexican tetra is composed of interdigitating bone, while the symphyseal region of the zebrafish is composed of fibrous connective tissue.

CONCLUSIONS

These differences in the jaw skeleton have likely evolved due to different feeding strategies utilised by each species. Our results show that premaxillae shape correlates strongly with the development of complex dentitions unlike in the mandible. This study provides important insights into the relationship between jaw and tooth development in bony fishes and suggests that these mechanisms may be similar amongst vertebrates.

Investigation of the effects of estrogen on skeletal gene expression during zebrafish larval head development

The development of craniofacial skeletal structures requires well-orchestrated tissue interactions controlled by distinct molecular signals. Disruptions in normal function of these molecular signals have been associated with a wide range of craniofacial malformations. A pathway mediated by estrogens is one of those molecular signals that plays role in formation of bone and cartilage including craniofacial skeletogenesis. Studies in zebrafish have shown that while higher concentrations of 17-β estradiol (E 2) cause severe craniofacial defects, treatment with lower concentrations result in subtle changes in head morphology characterized with shorter snouts and flatter faces. The molecular basis for these morphological changes, particularly the subtle skeletal effects mediated by lower E 2 concentrations, remains unexplored. In the present study we address these effects at a molecular level by quantitative expression analysis of sets of candidate genes in developing heads of zebrafish larvae treated with two different E 2 concentrations. To this end, we first validated three suitable reference genes, ppia2, rpl8 and tbp, to permit sensitive quantitative real-time PCR analysis. Next, we profiled the expression of 28 skeletogenesis-associated genes that potentially respond to estrogen signals and play role in craniofacial development. We found E 2 mediated differential expression of genes involved in extracellular matrix (ECM) remodelling, mmp2/9/13, sparc and timp2a, as well as components of skeletogenic pathways, bmp2a, erf, ptch1/2, rankl, rarab and sfrp1a. Furthermore, we identified a co-expressed network of genes, including cpn1, dnajc3, esr1, lman1, rrbp1a, ssr1 and tram1 with a stronger inductive response to a lower dose of E 2 during larval head development.

Developmental basis of phenotypic integration in two Lake Malawi cichlids

Background

Cichlid fishes from the Rift Lakes of East Africa have undergone the most spectacular adaptive radiations in vertebrate history. Eco-morphological adaptations in lakes Victoria, Malawi and Tanganyika have resulted in a vast array of skull shapes and sizes, yet primary axes of morphological variation are conserved in all three radiations, prominently including the size of the preorbital region of the skull. This conserved pattern suggests that development may constrain the trajectories of cichlid head morphological evolution.

Results

Here, we (1) present a comparative analysis of adult head morphology in two sand-dweller cichlids from Lake Malawi with preorbital size differences representative of the main axis of variation among the three lakes and (2) analyze the ontogeny of shape and size differences by focusing on known developmental modules throughout the head. We find that (1) developmental differences between the two species correlate with known developmental modules; (2) differences in embryonic cartilage development result in phenotypically integrated changes among all bones derived from a single cartilage, while differences in dermal bone development tend to influence isolated regions within a bone; and lastly (3) species-specific morphologies appear in the embryo as subtle differences, which become progressively amplified throughout ontogeny. We propose that this amplification takes place at skeletal growth zones, the locations and shapes of which are patterned during embryogenesis.

Conclusions

This study is the most anatomically comprehensive analysis of the developmental differences underlying cichlid skull evolution in the Rift Lakes of East Africa. The scale of our analysis reveals previously unnoticed correlations between developmental modules and patterns of phenotypic integration. We propose that the primary axes of morphological variation among East African cichlid adaptive radiations are constrained by the hierarchical modularity of the teleost head skeleton.

Differential expression of the aryl hydrocarbon receptor pathway associates with craniofacial polymorphism in sympatric Arctic charr

Background

The developmental basis of craniofacial morphology hinges on interactions of numerous signalling systems. Extensive craniofacial variation in the polymorphic Arctic charr, a member of the salmonid family, from Lake Thingvallavatn (Iceland), offers opportunities to find and study such signalling pathways and their key regulators, thereby shedding light on the developmental pathways, and the genetics of trophic divergence.

Results

To identify genes involved in the craniofacial differences between benthic and limnetic Arctic charr, we used transcriptome data from different morphs, spanning early development, together with data on craniofacial expression patterns and skeletogenesis in model vertebrate species. Out of 20 genes identified, 7 showed lower gene expression in benthic than in limnetic charr morphs. We had previously identified a conserved gene network involved in extracellular matrix (ECM) organization and skeletogenesis, showing higher expression in developing craniofacial elements of benthic than in limnetic Arctic charr morphs. The present study adds a second set of genes constituting an expanded gene network with strong, benthic–limnetic differential expression. To identify putative upstream regulators, we performed knowledge-based motif enrichment analyses on the regulatory sequences of the identified genes which yielded potential binding sites for a set of known transcription factors (TFs). Of the 8 TFs that we examined using qPCR, two (Ahr2b and Ap2) were found to be differentially expressed between benthic and limnetic charr. Expression analysis of several known AhR targets indicated higher activity of the AhR pathway during craniofacial development in benthic charr morphotypes.

Conclusion

These results suggest a key role of the aryl hydrocarbon receptor (AhR) pathway in the observed craniofacial differences between distinct charr morphotypes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0022-6) contains supplementary material, which is available to authorized users.