Bracketing phenogenotypic limits of mammalian hybridization

Toward a universal measure of robustness across model organs and systems

The development of an individual must be capable of resisting the harmful effects of internal and external perturbations. This capacity, called robustness, can make the difference between normal variation and disease. Some systems and organs are more resilient in their capacity to correct the effects of internal disturbances such as mutations. Similarly, organs and organisms differ in their capacity to be resilient against external disturbances, such as changes in temperature. Furthermore, all developmental systems must be somewhat flexible to permit evolutionary change, and understanding robustness requires a comparative framework. Over the last decades, most research on developmental robustness has been focusing on specific model systems and organs. Hence, we lack tools that would allow cross-species and cross-organ comparisons. Here, we emphasize the need for a uniform framework to experimentally test and quantify robustness across study systems and suggest that the analysis of fluctuating asymmetry might be a powerful proxy to do so. Such a comparative framework will ultimately help to resolve why and how organs of the same and different species differ in their sensitivity to internal (e.g., mutations) and external (e.g., temperature) perturbations and at what level of biological organization buffering capacities exist and therefore create robustness of the developmental system.

Genomic evidence uncovers inbreeding and supports translocations in rescuing the genetic diversity of a landlocked seal population

Fragmentation of isolated populations increases the risk of inbreeding and loss of genetic diversity. The endemic Saimaa ringed seal (Pusa hispida saimensis) is one of the most endangered pinnipeds in the world with a population of only ~ 400 individuals. The current genetic diversity of this subspecies, isolated in Lake Saimaa in Finland for ca. 1000 generations, is alarmingly low. We performed whole-genome sequencing on Saimaa ringed seals (N = 30) and analyzed the level of homozygosity and genetic composition across the individual genomes. Our results show that the Saimaa ringed seal population has a high number of runs of homozygosity (RoH) compared with the neighboring Baltic ringed seal (Pusa hispida botnica) reference population (p < 0.001). There is also a tendency for stillborn seal pups to have more pronounced RoH. Since the population is divided into semi-isolated subpopulations within the Lake Saimaa exposing the population to deleterious genomic effects, our results support augmented gene flow as a genetic conservation action. Based on our results suggesting inbreeding depression in the population, we recommend Pihlajavesi as a potential source and Southern Saimaa as a potential recipient subpopulation for translocating individuals. The Saimaa ringed seal is a recognized subspecies and therefore translocations should be considered only within the lake to avoid an unpredictable risk of disease, the introduction of deleterious alleles, and severe ecological issues for the population.

Skull ecomorphological variation of narwhals (Monodon monoceros, Linnaeus 1758) and belugas (Delphinapterus leucas, Pallas 1776) reveals phenotype of their hybrids

Narwhals and belugas are toothed whales belonging to the Monodontidae. Belugas have a circumpolar Arctic and sub-Artic distribution while narwhals are restricted to the Atlantic Arctic. Their geographical ranges overlap during winter migrations in the Baffin Bay area (Canada/West Greenland) and successful interbreeding may occur. Here, we employed geometric morphometrics on museum specimens to explore the cranium and mandible morphology of a known hybrid (NHMD MCE 1356) and the cranium morphology of a putative hybrid (NHMD 1963.44.1.4) relative to skull morphological variation in the parental species. Specifically, we used 3D models of skulls from 69 belugas, 86 narwhals, and the two known/putative hybrids and 2D left hemi-mandibles from 20 belugas, 64 narwhals and the known hybrid. Skull shape analyses allowed clear discrimination between species. Narwhals are characterised by a relatively short rostrum and wide neurocranium while belugas show a more elongated and narrower cranium. Sexual size dimorphism was detected in narwhals, with males larger than females, but no sexual shape dimorphism was detected in either species (excluding presence/absence of tusks in narwhals). Morphological skull variation was also dependent on different allometric slopes between species and sexes in narwhals. Our analyses showed that the cranium of the known hybrid was phenotypically close to belugas but its 2D hemi-mandible had a narwhal shape and size morphology. Both cranium and mandible were strongly correlated, with the pattern of covariation being similar to belugas. The putative hybrid was a pure male narwhal with extruded teeth. Comparison of genomic DNA supported this result, and stable carbon and nitrogen isotope values suggested that the putative hybrid had a more benthic foraging strategy compared to narwhals. This work demonstrates that although the known hybrid could be discriminated from narwhals and belugas, detection of its affinities with these parental species was dependent on the part of the skull analysed.

Range-wide variation in grey seal (Halichoerus grypus) skull morphology

The large interspecific variation in marine mammal skull and dental morphology reflects ecological specialisations to foraging and communication. At the intraspecific level, the drivers of skull shape variation are less well understood, having implications for identifying putative local foraging adaptations and delineating populations and subspecies for taxonomy, systematics, management and conservation. Here, we assess the range-wide intraspecific variation in 71 grey seal skulls by 3D surface scanning, collection of cranial landmarks and geometric morphometric analysis. We find that skull shape differs slightly between populations in the Northwest Atlantic, Northeast Atlantic and Baltic Sea. However, there was a large shape overlap between populations and variation was substantially larger among animals within populations than between. We hypothesize that this pattern of intraspecific variation in grey seal skull shape results from balancing selection or phenotypic plasticity allowing for a remarkably generalist foraging behaviour. Moreover, the large overlap in skull shape between populations implies that the separate subspecies status of Atlantic and Baltic Sea grey seals is questionable from a morphological point of view.

An epithelial signalling centre in sharks supports homology of tooth morphogenesis in vertebrates

Development of tooth shape is regulated by the enamel knot signalling centre, at least in mammals. Fgf signalling regulates differential proliferation between the enamel knot and adjacent dental epithelia during tooth development, leading to formation of the dental cusp. The presence of an enamel knot in non-mammalian vertebrates is debated given differences in signalling. Here, we show the conservation and restriction of fgf3, fgf10, and shh to the sites of future dental cusps in the shark (Scyliorhinus canicula), whilst also highlighting striking differences between the shark and mouse. We reveal shifts in tooth size, shape, and cusp number following small molecule perturbations of canonical Wnt signalling. Resulting tooth phenotypes mirror observed effects in mammals, where canonical Wnt has been implicated as an upstream regulator of enamel knot signalling. In silico modelling of shark dental morphogenesis demonstrates how subtle changes in activatory and inhibitory signals can alter tooth shape, resembling developmental phenotypes and cusp shapes observed following experimental Wnt perturbation. Our results support the functional conservation of an enamel knot-like signalling centre throughout vertebrates and suggest that varied tooth types from sharks to mammals follow a similar developmental bauplan. Lineage-specific differences in signalling are not sufficient in refuting homology of this signalling centre, which is likely older than teeth themselves.

Genomic evaluation of hybridization in historic and modern North American Bison (Bison bison)

During the late nineteenth century North American bison underwent a significant population bottleneck resulting in a reduction in population size of over 99% and a species-level near-extinction event. Factors responsible for this destruction included indiscriminate killing, loss of access to suitable habitat, and diseases. At the nadir of this population crash, very few wild plains bison survived and were restricted to Yellowstone National Park, USA and a small number of wild wood bison remained in Wood Buffalo National Park, Canada. However, most surviving bison in the late 1800’s were maintained by cattle ranchers in private herds where hybridization between bison with various breeds of domestic cattle was often encouraged. Over the last 20 years, the legacy of this introgression has been identified using mitochondrial DNA and limited nuclear microsatellite analyses. However, no genome-wide assessment has been performed, and some herds were believed to be free of introgression based on current genetic testing strategies. Herein, we report detailed analyses using whole genome sequencing from nineteen modern and six historical bison, chosen to represent the major lineages of bison, to identify and quantitate signatures of nuclear introgression in their recent (within 200 years) history. Both low and high coverage genomes provided evidence for recent introgression, including animals from Yellowstone, Wind Cave, and Elk Island National Parks which were previously thought to be free from hybridization with domestic cattle. We employed multiple approaches, including one developed for this work, to identify putative cattle haplotypes in each bison genome. These regions vary greatly in size and frequency by sample and herd, though we detected domestic cattle introgression in all bison genomes tested. Since our sampling strategy spanned across the diversity of modern bison populations, these finding are best explained by multiple historical hybridization events between these two species with significant genetic recombination over the last 200 years. Our results demonstrate that whole genome sequencing approaches are required to accurately quantitate cattle introgression in bison.

Origin and expansion of the world's most widespread pinniped: range-wide population genomics of the harbour seal (Phoca vitulina)

The harbour seal (Phoca vitulina) is the most widely distributed pinniped, occupying a wide variety of habitats and climatic zones across the Northern Hemisphere. Intriguingly, the harbour seal is also one of the most philopatric seals, raising questions as to how it colonised virtually the whole of the Northern Hemisphere. To shed light on the origin, remarkable range expansion, population structure and genetic diversity of this species, we used genotyping-by-sequencing to analyse ~13,500 biallelic SNPs from 286 individuals sampled from 22 localities across the species' range. Our results point to a Northeast Pacific origin, colonisation of the North Atlantic via the Canadian Arctic, and subsequent stepping-stone range expansions across the North Atlantic from North America to Europe, accompanied by a successive loss of genetic diversity. Our analyses further revealed a deep divergence between modern North Pacific and North Atlantic harbour seals, with finer-scale genetic structure at regional and local scales consistent with strong philopatry. The study provides new insights into the harbour seal's remarkable ability to colonise and adapt to a wide range of habitats. Furthermore, it has implications for current harbour seal subspecies delineations and highlights the need for international and national red lists and management plans to ensure the protection of genetically and demographically isolated populations.

Wild versus lab house mice: Effects of age, diet, and genetics on molar geometry and topography

Molar morphology is shaped by phylogenetic history and adaptive processes related to food processing. Topographic parameters of the occlusal surface, such as sharpness and relief, can be especially informative regarding diet preferences of a species. The occlusal surface can however be deeply modified by wear throughout an animal's life, potentially obliterating other signals. Age being difficult to assess in wild populations, especially small rodents, experimental studies of wear through age in laboratory populations may constitute a powerful way to assess its impact on molar geometry and topography, and to validate descriptors of molar morphology that could mitigate this issue. Molar morphology was therefore quantified using 3D geometric morphometrics and topographic estimates in four groups of house mice: wild-trapped mice, lab-bred offspring of these wild mice, typical laboratory mice, and their hybrids. Three descriptors of the molar morphology were considered: the surface of the whole molar row, the surface of the first upper molar, and a truncated template of the first upper molar mimicking advanced wear. Increasing wear with age was demonstrated in the different groups, with a more pronounced effect in the wild-trapped population. The geometry of the molar row is not only modified by wear, but also by the relative position of the late developing molars on the jaw due to loading during mastication. As a consequence, the alignment of the molars is modified in wild mice, showing a qualitative difference between wild animals and their lab-bred offspring. Results obtained from the lab should thus be transferred with caution to the interpretation of differences in wild populations. Topographic estimates computed for the first upper molar seems to provide more stable parameters than those based on the whole molar row, because issues related to non-planar occlusal surface along the molar row are discarded. The truncated template was proven efficient in discarding the wear effect to focus on genetic differences, allowing an efficient characterization of the hybridization signature between wild and lab mice. Dominance of the wild phenotype for the first molar shape supports that the lab strain evolved in a context of relaxation of the selective pressures related to nutrition.

Effects of hybridization on pelvic morphology: A macaque model

Ancient DNA analyses have shown that interbreeding between hominin taxa occurred multiple times. Although admixture is often reflected in skeletal phenotype, the relationship between the two remains poorly understood, hampering interpretation of the hominin fossil record. Direct study of this relationship is often impossible due to the paucity of hominin fossils and difficulties retrieving ancient genetic material. Here, we use a sample of known ancestry hybrids between two closely related nonhuman primate taxa (Indian and Chinese Macaca mulatta) to investigate the effect of admixture on skeletal morphology. We focus on pelvic shape, which has potential fitness implications in hybrids, as mismatches between maternal pelvic and fetal cranial morphology are often fatal to mother and offspring. As the pelvis is also one of the skeletal regions that differs most between Homo sapiens and Neanderthals, investigating the pelvic consequences of interbreeding could be informative regarding the viability of their hybrids. We find that the effect of admixture in M. mulatta is small and proportional to the relatively small morphological difference between the parent taxa. Sexual dimorphism appears to be the main determinant of pelvic shape in M. mulatta. The lack of difference in pelvic shape between Chinese and Indian M. mulatta is in contrast to that between Neanderthals and H. sapiens, despite a similar split time (in generations) between the hybridizing pairs. Greater phenotypic divergence between hominins may relate to adaptations to disparate environments but may also highlight how the unique degree of cultural buffering in hominins allowed for greater neutral divergence. In contrast to some previous work identifying extreme morphologies in first- and second-generation hybrids, here the relationship between pelvic shape and admixture is linear. This linearity may be because most sampled animals have a multigenerational admixture history or because of relatively high constraints on the pelvis compared with other skeletal regions.

Influence of captivity and selection on limb long bone cross-sectional morphology of reindeer

The emergence of pastoralism and animal husbandry has been a critical point in the history of human evolution. Beyond profound behavioural changes in domesticated animals compared to wild ones, characterising the morphological changes associated with domestication process remains challenging. Because reindeer (Rangifer tarandus) can be considered to still be in the early phases of the domestication process, the study of modern populations provides a unique opportunity to examine the impact of captivity and selective breeding on skeletal changes. In this work, we investigated the morphological changes in long limb bone cross-sections using 137 wild and domestic reindeer individuals bred in free-range, in captivity or used for racing and pulling. The shape and shaft cortical thickness of the six long limb bones (i.e., humerus, radioulna, metacarpal, femur, tibia and metatarsal) were measured using a 2D-geometric morphometrics approach taking into account subspecies, sex, body mass and lifestyle differences. These bones are important to understanding functional morphological changes because they can provide information on feeding and locomotor behaviours, as well as on body propulsion and weight bearing. Apart from the effects of taxonomy, etho-ecology and sex, we have found that captivity and selection induced important variations in the size and body mass of modern reindeer. Our results also showed that patterns of variation in cortical bone thickness of long limb bone cross-sections were strongly impacted by body mass and human-imposed restrictions in roaming. This demonstrates that bone cross-sections can provide information on changes in locomotor, reproductive and feeding behaviours induced by the domestication process. These results are valuable not only for (paleo) biologists studying the impact of captivity and selection in ungulates but also for archaeologists exploring the origins of domestication and early herding strategies.

Predicting sample success for large-scale ancient DNA studies on marine mammals

In recent years, nonhuman ancient DNA studies have begun to focus on larger sample sizes and whole genomes, offering the potential to reveal exciting and hitherto unknown answers to ongoing biological and archaeological questions. However, one major limitation to such studies is the substantial financial and time investments still required during sample screening, due to uncertainty regarding successful sample selection. This study investigates the effect of a wide range of sample properties including latitude, sample age, skeletal element, collagen preservation, and context on endogenous content and DNA damage profiles for 317 ancient and historic pinniped samples collected from across the North Atlantic and surrounding regions. Using generalised linear and mixed-effect models, we found that a range of factors affected DNA preservation within each of the species under consideration. The most important findings were that endogenous content varied significantly within species according to context, the type of skeletal element, the collagen content and collection year. There also appears to be an effect of the sample's geographic origin, with samples from the Arctic generally showing higher endogenous content and lower damage rates. Both latitude and sample age were found to have significant relationships with damage levels, but only for walrus samples. Sex, ontogenetic age and extraction material preparation were not found to have any significant relationship with DNA preservation. Overall, skeletal element and sample context were found to be the most influential factors and should therefore be considered when selecting samples for large-scale ancient genome studies.

A universal power law for modelling the growth and form of teeth, claws, horns, thorns, beaks, and shells

Background

A major goal of evolutionary developmental biology is to discover general models and mechanisms that create the phenotypes of organisms. However, universal models of such fundamental growth and form are rare, presumably due to the limited number of physical laws and biological processes that influence growth. One such model is the logarithmic spiral, which has been purported to explain the growth of biological structures such as teeth, claws, horns, and beaks. However, the logarithmic spiral only describes the path of the structure through space, and cannot generate these shapes.

Results

Here we show a new universal model based on a power law between the radius of the structure and its length, which generates a shape called a ‘power cone’. We describe the underlying ‘power cascade’ model that explains the extreme diversity of tooth shapes in vertebrates, including humans, mammoths, sabre-toothed cats, tyrannosaurs and giant megalodon sharks. This model can be used to predict the age of mammals with ever-growing teeth, including elephants and rodents. We view this as the third general model of tooth development, along with the patterning cascade model for cusp number and spacing, and the inhibitory cascade model that predicts relative tooth size. Beyond the dentition, this new model also describes the growth of claws, horns, antlers and beaks of vertebrates, as well as the fangs and shells of invertebrates, and thorns and prickles of plants.

Conclusions

The power cone is generated when the radial power growth rate is unequal to the length power growth rate. The power cascade model operates independently of the logarithmic spiral and is present throughout diverse biological systems. The power cascade provides a mechanistic basis for the generation of these pointed structures across the tree of life.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-00990-w.

Identifying early stages of reindeer domestication in the archaeological record: a 3D morphological investigation on forelimb bones of modern populations from Fennoscandia

Reindeer herding probably developed during the Late Iron Age onwards and is still an important part of the subsistence and culture of many peoples in northern Eurasia. However, despite the importance of this husbandry in the history of these Arctic people, the period and place of the origin as well as the spread of domestic reindeer is still highly debated. Besides the existence of different breeding methods in these territories, identifying domesticated individuals in the archaeological record is complicated because reindeers are considered to still be in the early phases of the domestication process. Indeed, the traditional morphological markers used in zooarchaeology to decipher the domestication syndrome are hardly perceptible in these early stages. In this work, we propose solutions for identifying domestic reindeer bones using 3D geometric morphometrics on isolated elements from the long bones of the forelimb (i.e. humerus, radio-ulna and metacarpal). These bones are important to understand both the feeding behaviour and the mobility of reindeer, and the potential effect of load-carrying or draught in the case of domestic reindeer. We analysed 123 modern specimens from Fennoscandia, including the two interbreeding subspecies currently present in these territories: mountain reindeer (Rangifer tarandus tarandus) and forest reindeer (R.t. fennicus); and where the sex and the lifestyle were known (i.e. free-ranging, racing or draught and captive individuals). A good level of discrimination between the size and shape variables of the bones of the forelimb was found among both subspecies and sexes. Moreover, individuals bred in captivity had smaller bone elements and a thinner and more slender morphology than free-ranging individuals. This demonstrates that the long bones of the forelimb can provide information on changes in feeding and locomotor behaviour prompted by the domestication process, like control and/or reduction of mobility and food of individual reindeer by humans. This also demonstrates that analysis in 3D geometric morphometrics is useful in detecting reindeer incipient domestication markers. Our results can be used by archaeologists to trace the early stages of domestication from fossil reindeer remains, and aid in reconstructing the socio-economic changes of past Arctic populations over time.

Comparative analysis of squamate brains unveils multi-level variation in cerebellar architecture associated with locomotor specialization

Ecomorphological studies evaluating the impact of environmental and biological factors on the brain have so far focused on morphology or size measurements, and the ecological relevance of potential multi-level variations in brain architecture remains unclear in vertebrates. Here, we exploit the extraordinary ecomorphological diversity of squamates to assess brain phenotypic diversification with respect to locomotor specialization, by integrating single-cell distribution and transcriptomic data along with geometric morphometric, phylogenetic, and volumetric analysis of high-definition 3D models. We reveal significant changes in cerebellar shape and size as well as alternative spatial layouts of cortical neurons and dynamic gene expression that all correlate with locomotor behaviours. These findings show that locomotor mode is a strong predictor of cerebellar structure and pattern, suggesting that major behavioural transitions in squamates are evolutionarily correlated with mosaic brain changes. Furthermore, our study amplifies the concept of ‘cerebrotype’, initially proposed for vertebrate brain proportions, towards additional shape characters.

The cerebellum is critical in sensory-motor control and is structurally diverse across vertebrates. Here, the authors investigate the evolutionary relationship between locomotory mode and cerebellum architecture across squamates by integrating study of gene expression, cell distribution, and 3D morphology.

Unravelling the hybrid vigor in domestic equids: the effect of hybridization on bone shape variation and covariation

BACKGROUND

Hybridization has been widely practiced in plant and animal breeding as a means to enhance the quality and fitness of the organisms. In domestic equids, this hybrid vigor takes the form of improved physical and physiological characteristics, notably for strength or endurance. Because the offspring of horse and donkey is generally sterile, this widely recognized vigor is expressed in the first generation (F1). However, in the absence of recombination between the two parental genomes, F1 hybrids can be expected to be phenotypically intermediate between their parents which could potentially restrict the possibilities of an increase in overall fitness. In this study, we examine the morphology of the main limb bones of domestic horses, donkeys and their hybrids to investigate the phenotypic impact of hybridization on the locomotor system. We explore bone shape variation and covariation to gain insights into the morphological and functional expressions of the hybrid vigor commonly described in domestic equids.

RESULTS

Our data reveal the occurrence of transgressive effects on several bones in the F1 generation. The patterns of morphological integration further demonstrate that the developmental processes producing covariation are not disrupted by hybridization, contrary to functional ones.

CONCLUSIONS

These results suggest that an increase in overall fitness could be related to more flexibility in shape change in hybrids, except for the main forelimb long bones of which the morphology is strongly driven by muscle interactions. More broadly, this study illustrates the interest of investigating not only bone shape variation but also underlying processes, in order to contribute to better understanding how developmental and functional mechanisms are affected by hybridization.

Dental pulp exposure, periapical inflammation and suppurative osteomyelitis of the jaws in juvenile Baltic grey seals (Halichoerus grypus grypus) from the late 19th century

The systematic analysis of museum collections can provide important insights into the dental and skeletal pathology of wild mammals. Here we present a previously unreported type of dental defect and related skull pathology in five juvenile Baltic grey seals that had been collected in the course of a seal culling program along the Danish coast in 1889 and 1890. All five skulls exhibited openings into the pulp cavities at the crown tips of all (four animals) or two (one animal) canines as well as several incisors and (in one animal) also some anterior premolars. The affected teeth showed wide pulp cavities and thin dentin. Pulp exposure had caused infection, inflammation, and finally necrosis of the pulp. As was evidenced by the extensive radiolucency around the roots of the affected teeth, the inflammation had extended from the pulp into the periapical space, leading to apical periodontitis with extensive bone resorption. Further spreading of the inflammation into the surrounding bone regions had then caused suppurative osteomyelitis of the jaws. The postcanine teeth of the pathological individuals typically had dentin of normal thickness and, except for one specimen, did not exhibit pulp exposure. The condition may have been caused by a late onset of secondary and tertiary dentin formation that led to pulp exposure in anterior teeth exposed to intense wear. Future investigations could address a possible genetic causation of the condition in the studied grey seals.