Conserved but flexible modularity in the zebrafish skull: implications for craniofacial evolvability

Whole-body variational modularity in the zebrafish: an inside-out story of a model species

Actinopterygians are the most diversified clade of extant vertebrates. Their impressive morphological disparity bears witness to tremendous ecological diversity. Modularity, the organization of biological systems into quasi-independent anatomical/morphological units, is thought to increase evolvability of organisms and facilitate morphological diversification. Our study aims to quantify patterns of variational modularity in a model actinopterygian, the zebrafish (Danio rerio), using three-dimensional geometric morphometrics on osteological structures isolated from micro-CT scans. A total of 72 landmarks were digitized along cranial and postcranial ossified regions of 30 adult zebrafishes. Two methods were used to test modularity hypotheses, the covariance ratio and the distance matrix approach. We find strong support for two modules, one comprised paired fins and the other comprised median fins, that are best explained by functional properties of subcarangiform swimming. While the skull is tightly integrated with the rest of the body, its intrinsic integration is relatively weak supporting previous findings that the fish skull is a modular structure. Our results provide additional support for the recognition of similar hypotheses of modularity identified based on external morphology in various teleosts, and at least two variational modules are proposed. Thus, our results hint at the possibility that internal and external modularity patterns may be congruent.

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.

Flexible conservatism in the skull modularity of convergently evolved myrmecophagous placental mammals

Background

The skull of placental mammals constitutes one of the best studied systems for phenotypic modularity. Several studies have found strong evidence for the conserved presence of two- and six-module architectures, while the strength of trait correlations (integration) has been associated with major developmental processes such as somatic growth, muscle-bone interactions, and tooth eruption. Among placentals, ant- and termite-eating (myrmecophagy) represents an exemplar case of dietary convergence, accompanied by the selection of several cranial morphofunctional traits such as rostrum elongation, tooth loss, and mastication loss. Despite such drastic functional modifications, the covariance patterns of the skull of convergently evolved myrmecophagous placentals are yet to be studied in order to assess the potential consequences of this dietary shift on cranial modularity.

Results

Here, we performed a landmark-based morphometric analysis of cranial covariance patterns in 13 species of myrmecophagous placentals. Our analyses reveal that most myrmecophagous species present skulls divided into six to seven modules (depending on the confirmatory method used), with architectures similar to those of non-myrmecophagous placentals (therian six modules). Within-module integration is also similar to what was previously described for other placentals, suggesting that most covariance-generating processes are conserved across the clade. Nevertheless, we show that extreme rostrum elongation and tooth loss in myrmecophagid anteaters have resulted in a shift in intermodule correlations in the proximal region of the rostrum. Namely, the naso-frontal and maxillo-palatine regions are strongly correlated with the oro-nasal module, suggesting an integrated rostrum conserved from pre-natal developmental processes. In contrast, the similarly toothless pangolins show a weaker correlation between the anterior rostral modules, resembling the pattern of toothed placentals.

Conclusions

These results reveal that despite some integration shifts related to extreme functional and morphological features of myrmecophagous skulls, cranial modular architectures have conserved the typical mammalian scheme.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02030-9.

Dynamics of the Zebrafish Skeleton in Three Dimensions During Juvenile and Adult Development

Zebrafish are a valuable model for normal vertebrate skeletogenesis and the study of myriad bone disorders. Bones grow, ossify and change shape throughout the zebrafish lifetime, and 3D technologies allow us to examine skeletogenic processes in detail through late developmental stages. To facilitate analysis of shape, orientation and tissue density of skeletal elements throughout ontogeny and adulthood, we generated a high-resolution skeletal reference dataset of wild-type zebrafish development. Using microCT technology, we produced 3D models of the skeletons of individuals ranging from 12 to 25 mm standard length (SL). We analyzed the dynamics of skeletal density and volume as they increase during juvenile and adult growth. Our resource allows anatomical comparisons between meristic units within an individual-e.g., we show that the vertebral canal width increases posteriorly along the spine. Further, structures may be compared between individuals at different body sizes: we highlight the shape changes that the lower jaw undergoes as fish mature from juvenile to adult. We show that even reproductively mature adult zebrafish (17-25 mm SL) continue to undergo substantial changes in skeletal morphology and composition with continued adult growth. We provide a segmented model of the adult skull and a series of interactive 3D PDFs at a range of key stages. These resources allow changes in the skeleton to be assessed quantitatively and qualitatively through late stages of development, and can serve as anatomical references for both research and education.

Complex genetic architecture of three-dimensional craniofacial shape variation in domestic pigeons

Deciphering the genetic basis of vertebrate craniofacial variation is a longstanding biological problem with broad implications in evolution, development, and human pathology. One of the most stunning examples of craniofacial diversification is the adaptive radiation of birds, in which the beak serves essential roles in virtually every aspect of their life histories. The domestic pigeon (Columba livia) provides an exceptional opportunity to study the genetic underpinnings of craniofacial variation because of its unique balance of experimental accessibility and extraordinary phenotypic diversity within a single species. We used traditional and geometric morphometrics to quantify craniofacial variation in an F₂ laboratory cross derived from the straight-beaked Pomeranian Pouter and curved-beak Scandaroon pigeon breeds. Using a combination of genome-wide quantitative trait locus scans and multi-locus modeling, we identified a set of genetic loci associated with complex shape variation in the craniofacial skeleton, including beak shape, braincase shape, and mandible shape. Some of these loci control coordinated changes between different structures, while others explain variation in the size and shape of specific skull and jaw regions. We find that in domestic pigeons, a complex blend of both independent and coupled genetic effects underlie three-dimensional craniofacial morphology.

Variability and the primacy of the genotype

Cave animals and species flocks exhibit common evolutionary principles. In caves, all traits dependent on the information derived from light lose their biological function. Mutations destructive for such traits, but neutral for the organism as a whole, can persist and accumulate until a trait has vanished. Adaptive radiations start in ecosystems containing open niches. Here, selection on niche-specific traits, such as the viscerocranium in fish, is relaxed owing to the absence of competing species, and viscerocranial variability arises. It is transitorily high in recent and phylogenetically younger flocks, providing new phenotypes. It lessens and is completely lost after directional selection promotes the fixation of phenotypes that are best adapted. In cave animals and species flocks, single traits manifest phenotypic variability owing to relaxed selection. Like the eye in cave species, the viscerocranium can be classified a module, the development of which is encoded in gene regulatory networks. Mutations in these genes can result in new phenotypes. Regarding functionality, these mutations might be destructive and eliminated by selection, neutral and thus persisting, or beneficial and promoted to fixation by directional selection. Given the ancient heritage of teleostean fish, these gene regulatory networks might be prone to mutations at the same loci or to developmental reactions resulting in similar phenotypes in closely related or taxonomically and geographically distant species.

Identification of Individual Zebrafish (Danio rerio): A Refined Protocol for VIE Tagging Whilst Considering Animal Welfare and the Principles of the 3Rs

In aquatic ecology, studies have commonly employed a tagging technique known as visible implant elastomer (VIE). This method has not been widely adopted by the zebrafish research community and also lacks refinement with regard to animal welfare. The current paper introduces a new VIE tagging protocol, with the aim of improving existing tagging techniques by placing particular emphasis on the Three Rs. To improve animal welfare and fish survival, we added the use of an analgesic compound (lidocaine) through the marking procedure, followed by after-treatment with antiseptics (melaleuca, aloe vera, and PVP-I as active ingredients) to improve tissue regeneration and healing. The newly improved protocol has been quantitatively evaluated on different populations and age groups of zebrafish. This study will be useful to the scientific zebrafish community and to the wider field including biologist and aquarists, especially in consideration of animal welfare, where tagging techniques are considered as a potential noxious stimulus for fish.

The genetic basis of coordinated plasticity across functional units in a Lake Malawi cichlid mapping population

Adaptive radiations are often stereotypical, as populations repeatedly specialize along conserved environmental axes. Phenotypic plasticity may be similarly stereotypical, as individuals respond to environmental cues. These parallel patterns of variation, which are often consistent across traits, have led researchers to propose that plasticity can facilitate predictable patterns of evolution along environmental gradients. This "flexible stem" model of evolution raises questions about the genetic nature of plasticity, including how complex is the genetic basis for plasticity? Is plasticity across traits mediated by many distinct loci, or few "global" regulators? To address these questions, we reared a hybrid cichlid mapping population on alternate diet regimes mimicking an important environmental axis. We show that plasticity across an array of ecologically relevant traits is generally morphologically integrated, such that traits respond in a coordinated manner, especially those with overlapping function. Our genetic data are more ambiguous. While our mapping experiment provides little evidence for global genetic regulators of plasticity, these data do contain a genetic signal for the integration of plasticity across traits. Overall, our data suggest a compromise between genetic modularity, whereby plasticity may evolve independently across traits, and low level but widespread genetic integration, establishing the potential for plasticity to experience coordinated evolution.

Genetic Pleiotropy of Bone-Related Phenotypes: Insights from Osteoporosis

PURPOSE OF REVIEW

We summarize recent evidence on the shared genetics within and outside the musculoskeletal system (mostly related to bone density and osteoporosis).

RECENT FINDINGS

Osteoporosis is determined by an interplay between multiple genetic and environmental factors. Significant progress has been made regarding its genetic background revealing a number of robustly validated loci and respective pathways. However, pleiotropic factors affecting bone and other tissues are not well understood. The analytical methods proposed to test for potential associations between genetic variants and multiple phenotypes can be applied to bone-related data. A number of recent genetic studies have shown evidence of pleiotropy between bone density and other different phenotypes (traits, conditions, or diseases), within and outside the musculoskeletal system. Power benefits of combining correlated phenotypes, as well as unbiased discovery, make these studies promising. Studies in humans are supported by evidence from animal models. Drug development and repurposing should benefit from the pleiotropic approach. We believe that future studies should take into account shared genetics between the bone and related traits.

Ecomorphological divergence and habitat lability in the context of robust patterns of modularity in the cichlid feeding apparatus

BACKGROUND

Adaptive radiations are characterized by extreme and/or iterative phenotypic divergence; however, such variation does not accumulate evenly across an organism. Instead, it is often partitioned into sub-units, or modules, which can differentially respond to selection. While it is recognized that changing the pattern of modularity or the strength of covariation (integration) can influence the range or rate of morphological evolution, the relationship between shape variation and covariation remains unclear. For example, it is possible that rapid phenotypic change requires concomitant changes to the underlying covariance structure. Alternatively, repeated shifts between phenotypic states may be facilitated by a conserved covariance structure. Distinguishing between these scenarios will contribute to a better understanding of the factors that shape biodiversity. Here, we explore these questions using a diverse Lake Malawi cichlid species complex, Tropheops, that appears to partition habitat by depth.

RESULTS

We construct a phylogeny of Tropheops populations and use 3D geometric morphometrics to assess the shape of four bones involved in feeding (mandible, pharyngeal jaw, maxilla, pre-maxilla) in populations that inhabit deep versus shallow habitats. We next test numerous modularity hypotheses to understand whether fish at different depths are characterized by conserved or divergent patterns of modularity. We further examine rates of morphological evolution and disparity between habitats and among modules. Finally, we raise a single Tropheops species in environments mimicking deep or shallow habitats to discover whether plasticity can replicate the pattern of morphology, disparity, or modularity observed in natural populations.

CONCLUSIONS

Our data support the hypothesis that conserved patterns of modularity permit the evolution of divergent morphologies and may facilitate the repeated transitions between habitats. In addition, we find the lab-reared populations replicate many trends in the natural populations, which suggests that plasticity may be an important force in initiating depth transitions, priming the feeding apparatus for evolutionary change.

Skull morphology diverges between urban and rural populations of red foxes mirroring patterns of domestication and macroevolution

Human activity is drastically altering the habitat use of natural populations. This has been documented as a driver of phenotypic divergence in a number of wild animal populations. Here, we show that urban and rural populations of red foxes (Vulpes vulpes) from London and surrounding boroughs are divergent in skull traits. These changes are primarily found to be involved with snout length, with urban individuals tending to have shorter and wider muzzles relative to rural individuals, smaller braincases and reduced sexual dimorphism. Changes were widespread and related to muscle attachment sites and thus are likely driven by differing biomechanical demands of feeding or cognition between habitats. Through extensive sampling of the genus Vulpes, we found no support for phylogenetic effects on skull morphology, but patterns of divergence found between urban and rural habitats in V. vulpes quantitatively aligned with macroevolutionary divergence between species. The patterns of skull divergence between urban and rural habitats matched the description of morphological changes that can occur during domestication. Specifically, urban populations of foxes show variation consistent with 'domestication syndrome'. Therefore, we suggest that occurrences of phenotypic divergence in relation to human activity, while interesting themselves, also have the potential to inform us of the conditions and mechanisms that could initiate domestication. Finally, this also suggests that patterns of domestication may be developmentally biased towards larger patterns of interspecific divergence.

Kinematic integration during prey capture varies among individuals but not ecological contexts in bluegill sunfish, Lepomis macrochirus (Perciformes: Centrarchidae)

The general ability of components of an organism to work together to achieve a common goal has been termed integration and is often studied empirically by deconstructing organisms into component parts and quantifying covariation between them. Kinematic traits describing movement are useful for allowing organisms to respond to ecological contexts that vary over short time spans (milliseconds, minutes, etc.). Integration of these traits can contribute to the maintenance of the function of the whole organism, but it is unclear how modulation of component kinematic traits affects their integration. We examined the integration of swimming and feeding during capture of alternative prey types in bluegill sunfish (Lepomis macrochirus). Despite the expected modulation of kinematics, integration within individuals was inflexible across prey types, suggesting functional redundancy for solving a broad constraint. However, integration was variable among individuals, suggesting that individuals vary in their solutions for achieving whole-organism function and that this solution acts as a ‘top-down’ regulator of component traits, which provides insight into why kinematic variation is observed. Additionally, variation in kinematic integration among individuals could serve as an understudied target of environmental selection on prey capture, which is a necessary first step towards the observed divergence in integration among populations and species.

Does phenotypic plasticity initiate developmental bias?

The generation of variation is paramount for the action of natural selection. Although biologists are now moving beyond the idea that random mutation provides the sole source of variation for adaptive evolution, we still assume that variation occurs randomly. In this review, we discuss an alternative view for how phenotypic plasticity, which has become well accepted as a source of phenotypic variation within evolutionary biology, can generate nonrandom variation. Although phenotypic plasticity is often defined as a property of a genotype, we argue that it needs to be considered more explicitly as a property of developmental systems involving more than the genotype. We provide examples of where plasticity could be initiating developmental bias, either through direct active responses to similar stimuli across populations or as the result of programmed variation within developmental systems. Such biased variation can echo past adaptations that reflect the evolutionary history of a lineage but can also serve to initiate evolution when environments change. Such adaptive programs can remain latent for millions of years and allow development to harbor an array of complex adaptations that can initiate new bouts of evolution. Specifically, we address how ideas such as the flexible stem hypothesis and cryptic genetic variation overlap, how modularity among traits can direct the outcomes of plasticity, and how the structure of developmental signaling pathways is limited to a few outcomes. We highlight key questions throughout and conclude by providing suggestions for future research that can address how plasticity initiates and harbors developmental bias.

Bony Patchwork: Mosaic Patterns of Evolution in the Skull of Electric Fishes (Apteronotidae: Gymnotiformes)

Mosaic evolution refers to the pattern whereby different organismal traits exhibit differential rates of evolution typically due to reduced levels of trait covariation through deep time (i.e., modularity). These differences in rates can be attributed to variation in responses to selective pressures between individual traits. Differential responses to selective pressures also have the potential to facilitate functional specialization, allowing certain traits to track environmental stimuli more closely than others. The teleost skull is a multifunctional structure comprising a complex network of bones and thus an excellent system for which to study mosaic evolution. Here we construct an ultrametric phylogeny for a clade of Neotropical electric fishes (Apteronotidae: Gymnotiformes) and use three-dimensional geometric morphometrics to investigate patterns of mosaic evolution in the skull and jaws. We find strong support for a developmental, three-module hypothesis that consists of the face, braincase, and mandible, and we find that the mandible has evolved four times faster than its neighboring modules. We hypothesize that the functional specialization of the mandible in this group of fishes has allowed it to outpace the face and braincase and evolve in a more decoupled manner. We also hypothesize that this pattern of mosaicism may be widespread across other clades of teleost fishes.

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.

Conserved but flexible modularity in the zebrafish skull: implications for craniofacial evolvability

Morphological variation is the outward manifestation of development and provides fodder for adaptive evolution. Because of this contingency, evolution is often thought to be biased by developmental processes and functional interactions among structures, which are statistically detectable through forms of covariance among traits. This can take the form of substructures of integrated traits, termed modules, which together comprise patterns of variational modularity. While modularity is essential to an understanding of evolutionary potential, biologists currently have little understanding of its genetic basis and its temporal dynamics over generations. To address these open questions, we compared patterns of craniofacial modularity among laboratory strains, defined mutant lines and a wild population of zebrafish (Danio rerio). Our findings suggest that relatively simple genetic changes can have profound effects on covariance, without greatly affecting craniofacial shape. Moreover, we show that instead of completely deconstructing the covariance structure among sets of traits, mutations cause shifts among seemingly latent patterns of modularity suggesting that the skull may be predisposed towards a limited number of phenotypes. This new insight may serve to greatly increase the evolvability of a population by providing a range of 'preset' patterns of modularity that can appear readily and allow for rapid evolution.

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.

The modular organization of roe deer (Capreolus capreolus) body during ontogeny: the effects of sex and habitat

Background

As a small artiodactyl, the roe deer (Capreolus capreolus L.) is characterized by biological plasticity and great adaptability demonstrated by their survival under a wide variety of environmental conditions. In order to depict patterns of phenotypic variation of roe deer body this study aims to quantify variation during ontogenetic development and determine how sex-specific reproductive investment and non-uniform habitat differences relate to phenotypic variation and do these differential investments mold the patterns of phenotypic variation through modular organisation.

Results

Patterns of phenotypic correlation among body traits change during the ontogeny of roe deer, with differential influence of sex and habitat type. Modularity was found to be a feature of closed habitats with trunk+forelimbs+hindlimbs as the best supported integration/modularity hypothesis for both sexes. The indices of integration and evolvability vary with habitat type, age and sex where increased integration is followed by decreased evolvability.

Conclusion

This is the first study that quantifies patterns of correlation in the roe deer body and finds pronounced changes in correlation structure during ontogeny affected by sex and habitat type. The correlation structure of the roe deer body is developmentally written over the course of ontogeny but we do not exclude the influence of function on ontogenetic changes. Modularity arises with the onset of reproduction (subadults not being modular) and is differentially expressed in males and females from different habitats. Both adult males and females show modularity in primordial, closed habitats. Overall, all these findings are important as they provide support to the idea that modularity can evolve at the population level and change fast within a species.