Skip to main content

Advertisement

You are viewing the new article page. Let us know what you think. Return to old version

Review | Open | Published:

Creating diversity in mammalian facial morphology: a review of potential developmental mechanisms

Abstract

Mammals (class Mammalia) have evolved diverse craniofacial morphology to adapt to a wide range of ecological niches. However, the genetic and developmental mechanisms underlying the diversification of mammalian craniofacial morphology remain largely unknown. In this paper, we focus on the facial length and orofacial clefts of mammals and deduce potential mechanisms that produced diversity in mammalian facial morphology. Small-scale changes in facial morphology from the common ancestor, such as slight changes in facial length and the evolution of the midline cleft in some lineages of bats, could be attributed to heterochrony in facial bone ossification. In contrast, large-scale changes of facial morphology from the common ancestor, such as a truncated, widened face as well as the evolution of the bilateral cleft possessed by some bat species, could be brought about by changes in growth and patterning of the facial primordium (the facial processes) at the early stages of embryogenesis.

Morphological diversity in mammalian faces

Mammals (class Mammalia) are one of the major groups of vertebrates, containing over 5400 living species as well as abundant extinct species [1,2,3,4]. Living mammals consist of three major clades: monotremes (order Monotremata), marsupials (infraclass Marsupialia), and placentals (infraclass Placentalia; Fig. 1). Recent phylogenetics, including comparative phylogenomic studies, have lead to a general consensus concerning the deeper branches of the mammalian evolutionary tree, for example identifying four major clades within placentals: Xenarthra, Afrotheria, Laurasiatheria, and Euarchontoglires [5,6,7,8,9,10,11].

Fig. 1
figure1

Phylogeny of mammals and diversity of their craniofacial morphology. A, phylogeny of living mammals adapted from Tarver et al. [11], with phylogeny of Laurasiatheria based on Chen et al. [101]. B–M′, frontal and lateral views of mammalian heads: B and B′, platypus (Ornithorhynchus anatinus); C and C′, echidna (Tachyglossus aculeatus); D and D′, red kangaroo (Macropus rufus); E and E′, koala (Phascolarctos cinereus); F and F′, giant anteater (Myrmecophaga tridactyla); G and G′, Linnaeus’s two-toed sloth (Choloepus didactylus); H and H′, aardvark (Orycteropus afer); I and I′, rock hyrax (Procavia capensis); J and J′, domestic dog (Borzoi) (Canis lupus familiaris); K and K′, leopard (Panthera pardus); L and L′, common tree shrew (Tupaia glis); M and M′, common chimpanzee (Pan troglodytes)

Mammals have evolved diverse morphologies to adapt to a wide range of ecological niches [3, 4]. The morphological diversity of mammalian heads is especially remarkable, possibly due to the head’s fundamental role in sensing, communication, and feeding [12,13,14,15,16,17,18] (Fig. 1). For example, both long- and short-faced taxa are recognized in each mammalian group (Fig. 1). Craniofacial morphology in mammals has been quantitatively evaluated in each group by comparative morphological analyses, including modern geometric morphometrics (summarized in Table 1).

Table 1 Diversity of craniofacial morphology in mammals and recent studies evaluating this diversity using landmark-based geometric morphometrics

However, the genetic and developmental mechanisms underlying the diversification of mammalian craniofacial morphology remain largely unknown. In this review, we compiled the recent findings in the developmental genetics of mice, a model mammalian species, to attempt to deduce the potential diversification mechanisms of mammalian facial morphology. We also introduce the results of previous studies in which a strong correlation between the number of nucleotide tandem repeats within the Runx2 gene and the facial length in some placental mammals was reported. Finally, we focus on bats (order Chiroptera), which display a substantial degree of craniofacial diversity and discuss their potential as a model for understanding the evolution of mammalian craniofacial morphology.

Molecular and cellular mechanisms creating diversity in facial morphology uncovered by mouse transgenesis

Mouse transgenesis is a powerful tool to infer the function of genes related to vertebrate morphogenesis. We examine the phenotypes of transgenic mice to gain insights into the molecular and cellular mechanisms that produce morphological variation in mammalian faces. We focused on two developmental events: (1) growth and patterning of the facial primordium and (2) ossification of the facial bones that lead to a shortened face and the orofacial cleft (Table 2).

Table 2 The genes involved in shortening the face and making the orofacial cleft in mouse

Growth and patterning of the facial primordium

Formation of mammalian faces begins at the pharyngula stage of embryogenesis, through growth and fusion of the five facial processes: the frontonasal process (FNP), medial nasal processes (MNPs), lateral nasal processes (LNPs), maxillary processes (MAXs), and mandibular processes (MANs) [19]. In the facial development of mice, FNP first expands anteriorly in a nine-day-old embryo (E9.0). Subsequently, MNPs and LNPs start to bulge out from the FNP at E10.0. These two processes surround the nasal placodes, MNP surrounds its medial aspect, and LNP surrounds its lateral aspect. During the same embryonic stage, MAXs begin to bulge anteriorly covering the ventrolateral aspect of the FNP. MAXs and the FNP continue to grow and fuse to each other in later stages to form the upper jaw. Paired MANs begin to grow anteriorly at E9.0 and fuse to one another at the midline to form the mandible [19, 20].

The early patterning of the mammalian face is regulated by migration and proliferation of the neural crest-derived mesenchyme (ectomesenchyme hereafter) [19, 21]. Mice with genetic defects related to the migration or proliferation of the ectomesenchyme possess a shortened face [22,23,24,25] and/or cleft lip (CL) occasionally accompanying the cleft palate (CP) [19, 26,27,28].

Several major signaling pathways, including BMP, FGF, Shh, and Wnt signaling pathways, are associated with outgrowth and fusion of the facial processes [19]. Repression of the up-stream component genes of these signaling pathways (e.g., Bmp4, Fgf8, Shh, and Wnt3) leads to a truncated face [19, 22, 24, 29, 30]. Recent papers have reported that migration of ectomesenchyme in the heads of mouse embryos are directly regulated by Wnt5a, a ligand of non-canonical Wnt signaling pathway [22, 25, 31, 32]. Alteration of the level of neural crest-specific Wnt5a expression (by both knockout and over-expression) results in a widened, shortened face [25, 33]. In Wnt5a conditional knockout mice, the migration pattern of the ectomesenchyme that later occupies the internal space of the facial processes is altered from that in control wild type mice [25]. The change in the ectomesenchyme migration pattern was attributed to the disruption of the directionality of cell division [25]. The induction of the internal facial structures (e.g., cartilage, bones, sensory compartments, muscles, glands, and teeth) was not influenced, and the lower jaw’s volume in the Wnt5a conditional knockout mouse was almost equivalent to that of the control mouse [25]. These results suggest that Wnt5a could play a crucial role in generating a shortened, widened face (truncated face) as naturally seen in koalas, sloths, the great apes, and cats through regulating the ectomesenchyme’s migration pattern, which in turn governs growth and organization of the facial processes (Fig. 1).

Disruptions in the growth and fusion of the facial processes also cause CL with or without CP (collectively called ‘CL/P’) [26,27,28]. A fusion of the facial processes first occurs between LNP and MNP, followed by a fusion of LNP and MAX. Finally, the anterior ends of both MAX and MNP are fused to one another. Fusion of the facial processes is initiated by contact of the epithelium of each facial process through proper organization of the facial processes [19]. Subsequently, the epithelial seam between coadjacent facial processes disappears due to apoptosis. Fusion of the MNP and the MAX and fusion of the MNP and the LNP are defective in mutants of the genes (e.g., Bmp4, Bmpr1a, Tcfap2a, Sox11, and Wnt9b) that regulate apoptosis within the epithelium as well as outgrowth and organization of the facial processes. Failure of these facial processes fusing accompanies CL/P [26].

Ossification of the facial bones

The palate of mammals separates the oral cavity from the nasal cavity and is subdivided into the anterior bony hard palate (palatal bones) and posterior soft palate [34]. The formation of the palate (palatogenesis) proceeds in two steps, the primary and secondary palate formations. In mouse development, the primary palate is formed by the fusion of the MAXs and MNPs at E11.5. Subsequently, the secondary palate is formed through three consecutive events. First, a pair of palatal shelves is formed by an uplift of the tongue at E11.5. Second, at E14.5, each palatal shelf grows medially above the tongue through ‘palatal shelf elevation’ [34]. Third, the left and right palatal shelves meet and fuse at the midline at E15.0 with fusion completing at E17.0. Palatal bones (anterior premaxilla derived from the ectomesenchyme of the primary palate, and central maxilla and posterior palatine that are derived from the ectomesenchyme of the secondary palate) begin to form at E14.5.

In contrast to defects in facial process development that produce an extremely shortened face (see the previous section), defects in facial bone formation, which occur in later phases of facial development, lead to a shortened face with milder dysmorphology. For example, Fbln5 knockout mice exhibit decreased outgrowth of the premaxilla bones during postnatal stages, compared to control wild type mice [35]. Fibulin-5 is an extracellular matrix protein deposited as a fibrous matrix in neural crest-derived craniofacial suture mesenchyme and plays a role as a regulator of cellular function such as cell proliferation [35, 36]. While premaxilla-maxilla suture mesenchyme in Fbln5 knockout mice were capable of differentiating into osteoblasts, suture cells in the mutant were less proliferative, suggesting fibulin-5 is indispensable for the regulation of facial suture mesenchymal cell proliferation required for craniofacial skeletal morphogenesis [35]. External facial morphology of adult Fbln5 knockout mice is almost normal, although facial length is slightly shortened compared to the control [35].

Defected facial bone development also leads to a submucous cleft palate (SMCP). SMCP is a clinical subgroup of CP. While CP is characterized by the whole palate (including both bones and epithelium) separated at the midline, SMCP is characterized by incomplete fusion of left and right palatal bones at the midline without cleft formation in the oral epithelium covering the bones. In mouse transgenesis, SMCP is only observed in the region between left and right maxilla bones. Only two genes that cause SMCP have been reported to date, Bmpr1a and Tbx22. In Osr2-IresCre;Bmpr1af/f transgenic mice, Bmpr1a was specifically knocked out in the tissue constructing the secondary palate. Osr2, whose promoter sequence was used for tissue/time-specific Bmpr1a knockout, is uniquely expressed in secondary palate morphogenesis in mice (see [37] for detail). The tissue-specific inactivation of Bmpr1a causes reduction of mesenchymal condensation in the anterior part of the secondary palate which subsequently differentiates into the maxilla bones [38]. Expression of Runx2, Osterix, and Dlx5, genes encoding transcriptional factors for bone development, is severely down-regulated in the anteromedial part of the secondary palate of Osr2-IresCre;Bmpr1af/f transgenic mice. As a result, elongation of the maxilla bones toward the midline is blocked, resulting in a cleft between the left and right maxilla bones [38]. Tbx22 is a transcription factor required for palatal bone formation [39]. Tbx22 knockout embryos bear a CP or SMCP accompanied by delayed osteoblast differentiation and hypotrophic maxilla bones [39].

To our knowledge, elongation of the face in transgenic mice compared to wild type mice has not been reported to date. In fish and birds, longer and more pointed jaws or beaks are formed by up-regulation of calmodulin signaling [40,41,42,43]. In mammals, however, the function of calmodulin signaling in facial development is poorly understood. Runx2 may regulate facial length in mammals. We briefly review the correlation between facial length and the variation of glutamine/alanine tandem repeats within Runx2 in the next section.

The number of Runx2 tandem repeats and mammalian facial length

There are long- and short-faced taxa in each mammalian group, and both face types show a high degree of diversity and evolvability in facial length (Fig. 1). Runx2 (Runt-related transcription factor 2) is an important transcription factor protein that plays multiple roles in bone development (e.g., osteoblast differentiation) in vertebrates including mammals [44,45,46] (reviewed in [47]). Runx2 enhances early osteoblast differentiation but inhibits terminal osteoblast differentiation [48]. Therefore, up-regulation of Runx2 leads to accelerated (via early onset of osteoblast differentiation) and extended (via delayed termination of osteoblast differentiation) bone development, while down-regulation of Runx2 results in delayed, shortened bone development [48, 49].

The Runx2 protein contains a highly conserved RUNT DNA binding domain and a repetitive glutamine (Q) and alanine (A) domain [46, 50]. Changes to the tandem repeat glutamines to alanines ratio (QA ratio), calculated by dividing the number of consecutive glutamines by the number of consecutive alanines within Runx2, alter transcriptional activity of Runx2 and its target genes [49, 51].

The Runx2 QA tandem repeat ratio is correlated with facial length variation in carnivorans [49, 52, 53]. Species with higher QA ratios have longer faces [49] (Fig. 2). In contrast, a lower QA ratio leads to lower transcriptional activity of Runx2 and results in short-faced carnivorans [49] (Fig. 2). This suggests that the QA ratio is associated with allometric variation in carnivoran facial length and the timing of facial bone (e.g., premaxilla, maxilla, nasal, jugal, vomer, palatine, and dentary) ossification. A similar pattern has been reported in primates [54].

Fig. 2
figure2

Correlation between QA ratio of the transcription factor Runx2 and facial length in order Carnivora. The Siberian husky, a breed of the domestic dog (Canis lupus familiaris), has 20 glutamine- and 8 alanine-coding nucleotide sequences within the repetitive glutamine and alanine domain of Runx2. QA ratio, calculated by dividing total glutamine-coding sequences by total alanine-coding sequences, is 2.50 and results in a longer face. The domestic cat (Falis catus) has 21 glutamine- and 10 alanine-coding nucleotide sequences within the corresponding domain of Runx2. QA ratio is 2.10 and results in a shorter face

Conversely, there is no correlation between the Runx2 QA tandem repeat ratio and facial length in xenarthrans and afrotherians [55], and marsupials [51]. Although marsupials display variation in facial length roughly equivalent to that observed in placentals (Fig. 1), almost no variation is observed in the nucleotide sequence of glutamine/alanine repeats in Runx2 [51]. The extreme conservation of nucleotide sequence and the QA ratio in marsupials may heavily constrain the timing of facial bone ossification in marsupial species [51]. These results suggest that the variations of facial length in xenarthrans, afrotherians, and marsupials are brought about by distinct molecular mechanisms. For example, a missense mutation in the gene Bmp3 (that encodes a growth factor, Bone morphogenetic protein 3) causes brachycephaly (shortened head) in domestic dogs [56]. We recommend further research concerning the role of morphogenetic genes such as Bmp3 to improve our understanding of the mechanisms generating facial length variation in mammals other than carnivorans and primates.

Bats: a model for understanding the diversification of mammalian craniofacial morphology

As reviewed in section II, our understanding of mammalian facial development mechanisms has been informed by studies of laboratory mice. However, the developmental mechanisms that produce facial morphology in non-model, wild mammal species have been only partially understood, perhaps due to difficulties in obtaining embryonic materials for analyses. More is understood about the molecular and cellular mechanisms underlying diversification of facial (beak) morphology in non-model bird species thanks to a series of evo-devo studies of Darwin’s finches, one of the most famous examples of adaptive radiations in vertebrates [40, 57,58,59,60,61,62]. Although model mammals help us to understand the basic mechanisms of mammalian morphogenesis, studying non-model species is necessary to identify other molecular and cellular mechanisms that lead to the morphological evolution of this group of vertebrates (including humans). Here, we focus on bats as a potential model for understanding evolution of mammalian craniofacial morphology.

Bats (order Chiroptera) are the second largest group of mammals after rodents [2, 63]. More than 1300 extant bat species are known, classified into 20 families [63]. Recent molecular phylogenetic studies [64,65,66,67] identified two major clades within bats, the Yinpterochiroptera and Yangochiroptera (Fig. 3). Chiropterans are distributed worldwide in all but the coldest regions [63], probably facilitated by the evolution of flight [68,69,70,71,72,73,74,75,76,77,78,79,80].

Fig. 3
figure3

Phylogeny of bats and evolution of orofacial cleft. The basic framework of chiropteran phylogeny is based on Teeling et al. [65]. Phylogenetic relationships in the superfamily Vespertilionoidea (Natalidae, Molossidae, Miniopteridae, Vespertilionidae, and Cistugidae) adapted from Lack et al. [228]. The midline cleft is possessed by nine different families of bats. In Molossidae, at least two genera (Mormopterus, and Tadarida) bear the midline cleft. In Furipteridae, only the genus Furipterus bears the midline cleft. The bilateral cleft evolved only once in the common ancestor of Rhinolophidae and Hipposideridae. Character mapping was based on Orr et al. [88]

Although largely neglected by biologists, diversity in bat facial morphology is astonishing. This diversity reflects their adaption to various environments and highly impressed Ernst Haeckel, an influential comparative embryologist and an artist in nineteenth century [81] (Fig. 4). New World leaf-nosed bats (family Phyllostomidae) are especially known for their incredible facial diversity [82, 83]. Phyllostmid facial length is strongly correlated with diet [84,85,86]. For example, frugivorous species (e.g., the wrinkle-faced bat, Centurio senex) have a truncated, widened face that exerts a high bite force. In contrast, nectarivorous species (e.g., the mexican long-tongued bat, Choeronycteris mexicana) have a long, narrow face that helps them to insert their rostrum into flowers. However, the molecular and cellular mechanisms that regulate the facial length of bats and are responsible for generating existing diversity in craniofacial morphology are poorly understood.

Fig. 4
figure4

Diversity of craniofacial morphology in bats. Left, a picture drawn by Ernst Haeckel, an influential comparative embryologist and artist [81]. Right, the silhouettes of the bat species illustrated in the Haeckel’s picture: (1) lesser long-eared bat (Nyctophilus geoffroyi), frontal view of the head; (2) brown long-eared bat (Plecotus auratus), frontal view of the head; (3) brown long-eared bat, entire body; (4) lesser false vampire bat (Megaderma spasma), frontal view of the head; (5) big-eared woolly bat (Chrotopterus auritus), lateral view of the head; (6) Tomes’s sword-nosed bat (Lonchorhina aurita), caudo-lateral view of the head; (7) Tomes’s sword-nosed bat, frontal view of the head; (8) Mexican funnel-eared bat (Natalus stramineus), frontal view of the head; (9) Antillean ghost-faced bat (Mormoops blainvillei), frontal view of the head; (10) flower-faced bat (Anthops ornatus), high magnification of noseleaf; (11) greater spear-nosed bat (Phyllostomus hastatus), frontal view of the head; (12) thumbless bat (Furipterus horrens), frontal view of the head; (13) greater horseshoe bat (Rhinolophus ferrumequinum), frontal view of the head; (14) wrinkle-faced bat (Centurio senex), frontal view of the head; (I) spectral bat (Vampyrum spectrum), frontal view of the head

Bats have a unique morphological feature in the rostral part of the upper jaw, an orofacial cleft on the premaxilla and maxilla bones that is anatomically similar to that observed in humans with congenital anomalies [87, 88]. There are two types of chiropteran orofacial cleft, midline and bilateral clefts. The midline cleft is observed in nine families of bats: Megadermatidae, Emballonuridae, Molossidae, Miniopteridae, Vespertilliionidae, Cistugidae, Myzopodidae, Thyropteridae, Furipteridae [88] (Fig. 3). Midline clefts are U-shaped clefts present between two premaxilla bones that are highly reduced in size (Fig. 6). Each premaxilla bone bears two permanent incisors and is completely fused to the maxilla bone posteriorly. The inner space of the cleft is occupied with a robust, translucent, fibrous membrane. The bilateral cleft is only seen in Rhinolophidae and Hipposideridae [88] (Fig. 3). In this cleft type, the premaxilla bone, which bears a single diminutive incisor, is separated from the laterally located maxilla bone by a cleft. The cleft is filled with fibrous connective tissue. The posterior margin of the medially fused premaxilla bones is loosely connected to the maxilla bones with fibrous connective tissue.

Bat orofacial clefts may contribute to reduction of returning echolocation signal interference, modulation of nasal acoustic emissions, increasing oral gape to facilitate capture of large prey, reduction of overall weight, and increase of olfactory ability [88]. However, the molecular and cellular mechanisms underlying orofacial cleft development in bats and the degree to which development of the two cleft types is similar are currently unknown.

Few studies have investigated the molecular mechanisms related to craniofacial diversity in bats. One such study by Phillips et al. [89] focused on Pax9, a transcription factor that plays an important role in vertebrate craniofacial and dental development. The authors compared nucleotide sequences of the 3′ untranslated region (UTR) of Pax9 among phyllostomids, vespertilionids, and other mammalian orders and identified four Musashi-binding elements (MBE) within conserved regions of the 3′ UTR [89]. The number of MBEs in morphologically diverse phyllostomid bats varied but was invariant in morphologically similar vespertilionid bats with the exception of a Murina species [89]. Because the number of MBEs may affect the expression level of Pax9, the authors proposed that the evolution of Pax9 regulation may be a contributing mechanism to the radiation of craniofacial morphological diversity in bats [89]. Although this study provides valuable insight into a potential genetic mechanism underlying the evolution and diversification of craniofacial morphology in phyllostomid bats, our understanding of the fundamental facial development mechanisms is far from complete.

Because convergence or parallel evolution of morphological traits in vertebrates is often brought about by identical genetic mechanisms (e.g., [90,91,92,93]), common mechanisms might regulate facial length even in bats (superorder Laurasiatheria) and rodents (superorder Euarchontoglires; Table 2).

In mice, a shortened face without apparent facial bone defects is mainly brought about by a decrease in proliferation and differentiation of the ectomesenchyme which later differentiates into osteoblasts [35]. In addition, facial length variation observed in carnivorans and primates are correlated with the level of activity of Runx2, which influences facial bone development duration [49]. Therefore, facial length variation in bats could be attributed to differences in the duration of facial bone development among species. For example, nectarivorous bats (e.g., Choeronycteris mexicana) have a relatively longer face. In this case, the duration of facial bone development might be extended, giving facial bones time to enlarged, especially anteriorly (Fig. 5). Conversely, insectivorous or omnivorous bats (e.g., Macrophyllum macrophyllum) have a relatively shorter face. Here, the period of facial bone development may be shortened leading to earlier completion of facial bone growth and preventing further anterior elongation (Fig. 5). Indeed, heterochronic shift in formation and growth of the palatal bones may produce variations of craniofacial morphology in phyllostomid bats [94]. Sears supposed that the diversity of palate shapes along phyllostomids is the result of relatively subtle evolutionary changes in later rather than earlier developmental event. Although it is likely that Runx2 plays a crucial role in producing facial length diversity in carnivorans and primates [49, 52,53,54], its function in chiropteran craniofacial development has yet to be identified and warrants further investigations.

Fig. 5
figure5

Schematic diagram depicting potential developmental mechanisms regulating facial length in bats. Top row: A long face possessed by some nectarivorous bat species (e.g., Choeronycteris mexicana) is formed through extension of the period of facial bone (e.g., premaxilla, maxilla, and dentary) development. Middle row: A short face possessed by many insectivorous or omnivorous bat species (e.g., Macrophyllum macrophyllum) is formed through shortening the period of facial bone development. Bottom row: A truncated face possessed by some frugivorous bat species (e.g., Centurio senex) is formed through deficient outgrowth of the facial processes in pharyngula stages. This could be attributed to reduced proliferation as well as disrupted migration of cranial neural crest cells (ectomesenchyme) occupying the internal space of the facial processes. FNP, frontonasal process; MAN, mandibular process; MAX, maxillary process

The truncated face of Wnt5a conditional knockout mice is brought about by the disruption of ectomesenchyme migration within the facial processes [25]. Notably, some phyllostomid bats (e.g., Centurio senex) possess an extremely truncated face that shares multiple characteristics with Wnt5a knockout mice faces. Therefore, facial morphology in these bat species might be derived from changes in expression of the genes that control direction of migration of the ectomesenchyme through regulating the directionality of cell division within the facial processes (Fig. 5). It would be interesting to compare Wnt5a activity and expression pattern in facial ectomesenchyme among chiropteran species.

The orofacial clefts observed in bats are morphologically categorized as SMCP. They are probably brought about by changes in premaxilla and maxilla bone formation. As we introduced in section II, Osr2-IresCre;Bmpr1af/f mice have a cleft between paired maxilla bones [38]. If Bmpr1a expression is specifically inactivated in the primary palate region using a similar transgenic technique (e.g., using a promoter of the gene that is uniquely expressed in the primary palate in gene knockout), a cleft may appear between paired premaxilla bones that are derived from the ectomesenchyme distributed within the primary palate. Considering this, the midline cleft in bats, which is present between two premaxilla bones, could be explained by domain-specific repression or down-regulation of Bmpr1a in the ectomesenchyme within the primary palate (instead of the secondary palate) that later gives rise to the premaxilla bones (Fig. 6). Because Bmpr1a is a receptor of the growth factor, the down-regulation of Bmpr1a may decrease the degree of ossification of the premaxilla bone through heterochrony (shorter and/or delayed ossification of the bone compared to the ancestor) and may result in such a small-scale morphological change in the tip of the face.

Fig. 6
figure6

Schematic diagram depicting potential developmental mechanisms forming orofacial clefts in bats. Top row: Most bat species do not have an orofacial cleft in their skull. Palatal bones are composed of three sets of bones: the premaxilla (blue), the maxilla (pink), and the palatine (gray). Each is a paired structure originally, having left and right elements. In adults, these bones are fused to each other and form the palate. In bats, the sutures between the bones constructing their skull (including the palate region) usually become indistinct through complete fusion of the bones. Middle row: In bat species with the midline cleft, the growth of the premaxilla bones toward the midline is inhibited and this makes medially unfused premaxilla bones. Heterochrony in ossification of the premaxilla bone (shorter and/or delayed ossification of the bone compared to the ancestor) may result in such a small-scale morphological change in the tip of the face. Bottom row: In bat species with the bilateral cleft, the cleft is likely formed through three developmental steps: (1) the domain of the secondary palate expands antero-medially, possibly through changes in growth and patterning of the facial processes at the early stages of embryogenesis. This narrows the space for the primary palate (arrows in the left illustration). (2) The maxilla bones are elongated anteriorly (arrows in the central illustration) compared to in bats species without orofacial clefts as well as those with midline clefts, acquiring its anterior projection. Simultaneously, the position of the premaxilla bones is confined at the center of the tip of the face, due to reduction of the space for its lateral expansion. (3) The boundary between the (anterior) premaxilla and (posterior) maxilla is left as a joint connected through loose connective tissue (a white dashed line in the right illustration). The space between the (medial) premaxilla and the (lateral) maxilla bones is left as a cleft

The formation of the bilateral cleft could be much more complicated, perhaps associated with extensive alterations of the developmental program. The premaxilla bones are derived from the ectomesenchyme distributed within the primordium of the primary palate, while the maxilla bones are derived from that of the secondary palate. Therefore, in the facial development of bat species bearing the bilateral cleft, the relative position of the primary and secondary palates might be changed through alterations in formation and organization of the facial processes from those in bat species without orofacial cleft. We speculate that the bilateral cleft developed through the following three steps (Fig. 6). First, the ectomesenchyme occupying the secondary palate expanded its distribution antero-medially and restricted the space for primary palate development at the tip of the face. Second, the osteoblasts derived from the ectomesenchyme distributed within the anterior part of the secondary palate differentiated into bone and made anterior projection of the maxilla bones surrounding the premaxilla bone laterally. Thus, the position of the premaxilla bone became restricted at the center of the tip of the face. Third, inhibition of ossification at the suture between the medially positioned premaxilla and laterally positioned maxilla bones left the unossified area between the two bones as a cleft.

Orofacial clefts occur as a craniofacial anomaly in humans at a relatively high frequency (approximately 1 in 700 live births) [88]. Investigating the mechanisms behind orofacial cleft formation in bats may contribute not only to understanding the reason why this cranial feature, which usually occurs as a skeletal pathology in other mammals groups including humans, appears as a normal phenotype in bats, but also to developing novel therapies against human orofacial cleft.

In the last 15 years, several studies have described in detail the overall embryonic development [95,96,97,98,99,100] and specifically wing development of bat species where embryos could be obtained [68,69,70,71,72,73,74,75,76,77, 79, 80]. We believe that examination of bat facial development and its comparisons among the species provide profound insights into the molecular and cellular bases of craniofacial morphology diversification in mammals.

Conclusions

In this paper, we have reviewed recent advances in understanding how mammalian faces are formed and discussed how these data are being applied to make new hypotheses about the diversity creation in mammalian craniofacial morphology. Small-scale changes in facial morphology from the ancestor, such slight changes in facial length and the evolution of the midline cleft in some lineages of bats could be attributed to heterochrony in facial bone ossification. In contrast, large-scale changes in facial morphology from the ancestor, such as a truncated, widened faces, as well as the evolution of the bilateral cleft in some bat species, could be brought about by changes in growth and patterning of the facial primordium (the facial processes) at the early stages of embryogenesis. Significant work remains to be done to test these hypotheses.

Abbreviations

CL:

cleft lip

CP:

cleft palate

FNP:

frontonasal process

LNP:

lateral nasal process

MAN:

mandibular process

MAX:

maxillary process

MBE:

Musashi-binding elements

MNP:

medial nasal process

SMCP:

submucous cleft palate

References

  1. 1.

    Kemp TS. The origin and evolution of mammals. Oxford: Oxford University Press; 2005.

  2. 2.

    Wilson DE, Reeder DM. Mammal species of the world. 3rd ed. Baltimore: Johns Hopkins University Press; 2005.

  3. 3.

    Vaughan TA, Ryan JM, Czaplewski NJ. Mammalogy. 6th ed. Burlington: Jones & Bartlett Learning; 2013.

  4. 4.

    Feldhamer GA, Drickamer LC, Vessey SH, Merritt JF, Krajewski C. Mammalogy: adaptation, diversity, ecology. 4th ed. Baltimore: Johns Hopkins University Press; 2014.

  5. 5.

    Meredith RW, Janečka JE, Gatesy J, Ryder OA, Fisher CA, Teeling EC, Goodbla A, Eizirik E, Simão TL, Stadler T, Rabosky DL, Honeycutt RL, Flynn JJ, Ingram CM, Steiner C, Williams TL, Robinson TJ, Burk-Herrick A, Westerman M, Ayoub NA, Springer MS, Murphy WJ. Impacts of the cretaceous terrestrial revolution and KPg extinction on mammal diversification. Science. 2011;334:521–4.

  6. 6.

    Dos Reis M, Inoue J, Hasegawa M, Asher RJ, Donoghue PC, Yang Z. Phylogenomic datasets provide both precision and accuracy in estimating the timescale of placental mammal phylogeny. Proc Biol Sci. 2012;279:3491–500.

  7. 7.

    Song S, Liu L, Edwards SV, Wu S. Resolving conflict in eutherian mammal phylogeny using phylogenomics and the multispecies coalescent model. Proc Nat Acad Sci. 2012;109:14942–7.

  8. 8.

    O’Leary MA, Bloch JI, Flynn JJ, Gaudin TJ, Giallombardo A, Giannini NP, Goldberg SL, Kraatz BP, Luo ZX, Meng J, Ni X, Novacek MJ, Perini FA, Randall ZS, Rougier GW, Sargis EJ, Silcox MT, Simmons NB, Spaulding M, Velazco PM, Weksler M, Wible JR, Cirranello AL. The placental mammal ancestor and the post-K-Pg radiation of placentals. Science. 2013;339:662–7.

  9. 9.

    Mitchell KJ, Pratt RC, Watson LN, Gibb GC, Llamas B, Kasper M, Edson J, Hopwood B, Male D, Armstrong KN, Meyer M, Hofreiter M, Austin J, Donnellan SC, Lee MS, Phillips MJ, Cooper A. Molecular phylogeny, biogeography, and habitat preference evolution of Marsupials. Mol Biol Evol. 2014;31:2322–30.

  10. 10.

    Foley NM, Springer MS, Teeling EC. Mammal madness: is the mammal tree of life not yet resolved? Philos Trans R Sci. 2016;371:20150140.

  11. 11.

    Tarver JE, Dos Reis M, Mirarab S, Moran RJ, Parker S, O’Reilly JE, King BL, O’Connell MJ, Asher RJ, Warnow T, Peterson KJ, Donoghue PC, Pisani D. The interrelationships of placental mammals and the limits of phylogenetic inference. Genome Biol Evol. 2016;8:330–44.

  12. 12.

    Marcus E, Hingst-Zaher E, Zaher H. Application of landmark morphometrics to skulls representing the orders of living mammals. Hystrix. 2000;11:27–47.

  13. 13.

    Brugmann SA, Kim J, Helms JA. Looking different: understanding diversity in facial form. Am J Med Genet A. 2006;140:2521–9.

  14. 14.

    Goswami A. Cranial modularity shifts during mammalian evolution. Am Nat. 2006;168:270–80.

  15. 15.

    Porto A, de Oliveira FB, Shirai LT, De Conto V, Marroig G. The Evolution of modularity in the mammalian skull I: morphological integration patterns and magnitudes. Evol Biol. 2009;36:118–35.

  16. 16.

    Koyabu D, Maier W, Sánchez-Villagra MR. Paleontological and developmental evidence resolve the homology and dual embryonic origin of a mammalian skull bone, the interparietal. Proc Natl Acad Sci USA. 2012;109:14075–80.

  17. 17.

    Cardini A, Polly PD. Larger mammals have longer faces because of size-related constraints on skull form. Nat Commun. 2013;4:2458.

  18. 18.

    Koyabu D, Werneburg I, Morimoto N, Zollikofer CP, Forasiepi AM, Endo H, Kimura J, Ohdachi SD, Truong Son N, Sánchez-Villagra MR. Mammalian skull heterochrony reveals modular evolution and a link between cranial development and brain size. Nat Commun. 2014;5:3625.

  19. 19.

    Jiang R, Bush JO, Lidral AC. Development of the upper lip: morphogenetic and molecular mechanisms. Dev Dyn. 2006;235:1152–66.

  20. 20.

    Tehiler K. The house mouse: atlas of embryonic development. New York: Springer; 1989.

  21. 21.

    O’Rahilly R, Müller F. The development of the neural crest in the human. J Anat. 2007;211:335–51.

  22. 22.

    Reid BS, Yang H, Melvin VS, Taketo MM, Williams T. Ectodermal Wnt/β-catenin signaling shapes the mouse face. Dev Biol. 2011;349:261–9.

  23. 23.

    Medio M, Yeh E, Popelut A, Babajko S, Berdal A, Helms JA. Wnt/β-catenin signaling and Msx1 promote outgrowth of the maxillary prominences. Front Physiol. 2012;3:375.

  24. 24.

    Lipinski RJ, Holloway HT, O’Leary-Moore SK, Ament JJ, Pecevich SJ, Cofer GP, Budin F, Everson JL, Johnson GA, Sulik KK. Characterization of subtle brain abnormalities in a mouse model of hedgehog pathway antagonist-induced cleft lip and palate. PLoS ONE. 2014;9:3–9.

  25. 25.

    Kaucka M, Ivashkin E, Gyllborg D, Zikmund T, Tesarova M, Kaiser J, Xie M, Petersen J, Pachni V, Nicolis SK, Yu T, Sharpe P, Arenas E, Brismar H, Blom H, Clevers H, Suter U, Chagin AS, Fried K, Hellander A, Adameyko I. Analysis of neural crest-derived clones reveals novel aspects of facial development. Sci Adv. 2016;2:e1600060.

  26. 26.

    Juriloff DM, Harris MJ. Mouse genetic models of cleft lip with or without cleft palate. Birth Defects Res A Clin Mol Teratol. 2008;82:63–77.

  27. 27.

    Jugessur A, Farlie P, Kilpatrick N. The genetics of isolated orofacial clefts: From genotypes to subphenotypes. Oral Dis. 2009;15:437–53.

  28. 28.

    Dixon MJ, Marazita ML, Beaty TH, Murray JC. Cleft lip and palate: Synthesizing genetic and environmental influences. Nat Rev Genet. 2011;123:167–78.

  29. 29.

    Lan Y, Ryan RC, Zhang Z, Bullard AA, Bush JO, Maltby KM, Lidral AC, Jiang R. Expression of Wnt9b and activation of canonical Wnt signaling during midfacial morphogenesis in mice. Dev Dyn. 2006;235:1448–54.

  30. 30.

    Bonilla-Claudio M, Wang J, Bai Y, Klysik E, Selever J, Martin JF. Bmp signaling regulates a dose-dependent transcriptional program to control facial skeletal development. Development. 2012;139:709–19.

  31. 31.

    Bakker ERM, Raghoebir L, Franken PF, Helvensteijn W, van Gurp L, Meijlink F, van der Valk MA, Rottier RJ, Kuipers EJ, van Veelen W, Smits R. Induced Wnt5a expression perturbs embryonic outgrowth and intestinal elongation, but is well-tolerated in adult mice. Dev Biol. 2012;369:91–100.

  32. 32.

    Ho HY, Susman MW, Bikoff JB, Ryu YK, Jonas AM, Hu L, Kuruvilla R, Greenberg ME. Wnt5a-Ror-Dishevelled signaling constitutes a core developmental pathway that controls tissue morphogenesis. Proc Natl Acad Sci USA. 2012;109:4044–51.

  33. 33.

    Amerongen RV, Fuerer C, Mizutani M, Nusse R. Wnt5a can both activate and repress Wnt/β-catenin signaling during mouse embryonic development. Dev Biol. 2012;369:101–14.

  34. 34.

    Bush JO, Jiang R. Palatogenesis: morphogenetic and molecular mechanisms of secondary palate development. Development. 2012;139:828.

  35. 35.

    Noda K, Nakamura T, Komatsu Y. Fibulin-5 deficiency causes developmental defect of premaxillary bone in mice. Biochem Biophys Res Commun. 2015;36:1011–4.

  36. 36.

    Yanagisawa H, Schluterman MK, Breken RA. Fibulin-5, an integrin-binding matricellular protein: its function in development and disease. J Cell Commun Signal. 2009;3:337–47.

  37. 37.

    Lan Y, Ovitt CE, Cho E, Maltby KM, Wang Q, Jiang R. Odd-skipped related 2 (Osr2) encodes a key intrinsic regulator of secondary palate growth and morphogenesis. Development. 2004;121:3207–16.

  38. 38.

    Baek JA, Lan Y, Liu H, Maltby KM, Mishina Y, Jiang R. Bmpr1a signaling plays critical roles in palatal shelf growth and palatal bone formation. Dev Biol. 2011;350:520–31.

  39. 39.

    Pauws E, Hoshino A, Bentley L, Prajapati S, Keller C, Hammond P, Martinez-Barbera J, Moore GE, Stanier P. Tbx22 null mice have a submucous cleft palate due to reduced palatal bone formation and also display ankyloglossia and choanal atresia phenotypes. Hum Mol Genet. 2009;18:4171–9.

  40. 40.

    Abzhanov A, Kuo WP, Hartman C, Grant BR, Grant PR, Tabin CJ. The calmodulin pathway and evolution of elongated beak morphology in Darwin’s finches. Nature. 2006;442:563–7.

  41. 41.

    Parsons KJ, Albertson RC. Roles for Bmp4 and Cam1 in shaping the jaw: Evo-Devo and beyond. Annu Rev Genet. 2009;43:369–88.

  42. 42.

    Knief U, Schielzeth H, Kempenaers B, Ellegren H, Forstmeier W. QTL and quantitative genetic analysis of beak morphology reveals patterns of standing genetic variation in an Estrildid finch. Mol Ecol. 2012;21:3704–17.

  43. 43.

    Gunter H, Koppermann C, Meyer A. Revisiting de Beer’s textbook example of heterochrony and jaw elongation in fish: calmodulin expression reflects heterochronic growth, and underlies morphological innovation in the jaws of belonoid fishes. EvoDevo. 2014;5:8.

  44. 44.

    Suda T, Kobayashi K, Jimi E, Udagawa N, Takahashi N. The molecular basis of osteoblast differentiation and activation. Novartis Found Symp. 2001;232:235–50.

  45. 45.

    Komori T. Regulation of skeletal development by the Runx family of transcription factors. J Cell Biochem. 2005;95:445–53.

  46. 46.

    Schroeder TM, Jensen ED, Westendorf JJ. Runx2: A master organizer of gene transcription in developing and maturing osteoblasts. Birth Defects Res C Embryo Today. 2005;75:213–25.

  47. 47.

    Tokita M, Chaeychomsri W, Siruntawineti J. Skeletal gene expression in the temporal region of the reptilian embryos: implications for the evolution of reptilian skull morphology. Springerplus. 2013;2:336.

  48. 48.

    Komori T. Runx2, A multifunctional transcription factor in skeletal development. J Cell Biochem. 2002;87:1–8.

  49. 49.

    Sears KE, Goswami A, Flynn JJ, Niswander LA. The correlated evolution of Runx2 tandem repeats, transcriptional activity, and facial length in Carnivora. Evol Dev. 2007;9:555–65.

  50. 50.

    Ziros PG, Basdra EK, Papavassiliou AG. Runx2: of bone and stretch. Int J Biochem Cell Biol. 2008;40:1659–63.

  51. 51.

    Newton AH, Feigin CY, Pask AJ. RUNX2 repeat variation does not drive craniofacial diversity in marsupials. BMC Evol Biol. 2017;17:110.

  52. 52.

    Fondon JW 3rd, Garner HR. Molecular origins of rapid and continuous morphological evolution. Proc Natl Acad Sci USA. 2004;101:18058–63.

  53. 53.

    Fondon JW 3rd, Garner HR. Detection of length-dependent effects of tandem repeat alleles by 3-D geometric decomposition of craniofacial variation. Dev Genes Evol. 2007;217:79–85.

  54. 54.

    Ritzman TB, Banovich N, Buss KP, Guida J, Rubel MA, Pinney J, Khang B, Ravosa MJ, Stone AC. Facing the facts: the Runx2 gene is associated with variation in facial morphology in primates. J Hum Evol. 2017;111:139–51.

  55. 55.

    Pointer MA, Kamilar JM, Warmuth V, Chester SGB, Delsuc F, Mundy NI, AsherEmail RJ, Bradley BJ. RUNX2 tandem repeats and the evolution of facial length in placental mammals. BMC Evol Biol. 2012;12:103.

  56. 56.

    Schoenebeck JJ, Hutchinson SA, Byers A, Beale HC, Carrington B, Faden DL, Rimbault M, Decker B, Kidd JM, Sood R, Boyko AR, Fondon JW 3rd, Wayne RK, Bustamante CD, Ciruna B, Ostrander EA. Variation of BMP3 contributes to dog breed skull diversity. PLoS Genet. 2012;8:1–11.

  57. 57.

    Abzhanov A, Protas M, Grant BR, Grant PR, Tabin C. Bmp4 and morphological variation of beaks in Darwin’s Finches. Science. 2004;305:1462–5.

  58. 58.

    Campàs O, Mallarino R, Herrel A, Abzhanov A. Brenner MP. Scaling and shear transformations capture beak shape variation in Darwin’s finches. Proc Natl Acad Sci USA. 2010;107:3356–60.

  59. 59.

    Mallarino R, Grant PR, Grant BR, Herrel A, Kuo WP, Abzhanova A. Two developmental modules establish 3D beak-shape variation in Darwin’s finches. Proc Natl Acad Sci USA. 2011;108:4057–62.

  60. 60.

    Mallarino R, Campàs O, Fritz J, Burns K, Weeks OG, Brenner MP, Abzhanov A. Closely related bird species demonstrate flexibility between beak morphology and underlying developmental programs. Proc Natl Acad Sci USA. 2012;109:16222–7.

  61. 61.

    Fritz JA, Brancale J, Tokita M, Burns KJ, Hawkins MB, Abzhanov A, Brenner MP. Shared developmental programme strongly constrains beak shape diversity in songbirds. Nat Commun. 2014;5:3700.

  62. 62.

    Tokita M, Yano W, James HF, Abzhanov A. Cranial shape evolution in adaptive radiations of birds: comparative morphometrics of Darwin’s finches and Hawaiian honeycreepers. Philos Trans R Soc Lond B Biol Sci. 2017;372:20150481.

  63. 63.

    Fenton MB, Simmons NB. Bats a world of science and mystery. Chicago: University of Chicago Press; 2014.

  64. 64.

    Eick GN, Jacobs DS, Matthee CA. A nuclear DNA phylogenetic perspective on the evolution of echolocation and historical biogeography of extant bats (Chiroptera). Mol Biol Evol. 2005;22:1869–86.

  65. 65.

    Teeling EC, Springer MS, Madsen O, Bates P, O’brien SJ, Murphy WJ. A molecular phylogeny for bats illuminates biogeography and the fossil record. Science. 2005;307:580–4.

  66. 66.

    Tsagkogeorga G, Parker J, Stupka E, Cotton JA, Rossiter SJ. Phylogenomic analyses elucidate the evolutionary relationships of bats. Curr Biol. 2013;23:2262–7.

  67. 67.

    Shi JJ, Rabosky DL. Speciation dynamics during the global radiation of extant bats. Evolution. 2015;69:1528–45.

  68. 68.

    Chen CH, Cretekos CJ, Rasweiler JJ 4th, Behringer RR. Hoxd13 expression in the developing limbs of the short-tailed fruit bat, Carollia perspicillata. Evol Dev. 2005;7:130–41.

  69. 69.

    Sears KE, Behringer RR, Rasweiler JJ, Niswander LA. Development of bat flight: morphologic and molecular evolution of bat wing digits. Proc Natl Acad Sci USA. 2006;103:6581–6.

  70. 70.

    Weatherbee SD, Behringer RR, Rasweiler JJ, Niswander LA. Interdigital webbing retention in bat wings illustrates genetic changes underlying amniote limb diversification. Proc Natl Acad Sci. 2006;103:15103–7.

  71. 71.

    Cretekos CJ, Wang Y, Green ED, Martin JF, Rasweiler JJ, Behringer RR. Regulatory divergence modifies limb length between mammals. Genes Dev. 2008;22:141–51.

  72. 72.

    Hockman D, Cretekos CJ, Mason MK, Behringer RR, Jacobs DS, Illing N. A second wave of Sonic hedgehog expression during the development of the bat limb. Proc Natl Acad Sci USA. 2008;105:16982–7.

  73. 73.

    Wang Z, Dong D, Ru B, Young RL, Han N, Guo T, Zhang S. Digital gene expression tag profiling of bat digits provides robust candidates contributing to wing formation. BMC Genom. 2010;11:619.

  74. 74.

    Tokita M, Abe T, Suzuki K. The developmental basis of bat wing muscle. Nat Commun. 2012;3:1302.

  75. 75.

    Dai M, Wang Y, Fang L, Irwin DM, Zhu T, Zhang J, Zhang S, Wang Z. Differential expression of Meis2, Mab21l2 and Tbx3 during limb development associated with diversification of limb morphology in mammals. PLoS ONE. 2014;9:e106100.

  76. 76.

    Wang Z, Dai M, Wang Y, Cooper KL, Zhu T, Dong D, Zhang J, Zhang S. Unique expression patterns of multiple key genes associated with the evolution of mammalian flight. Proc Biol Sci. 2014;281:20133133.

  77. 77.

    Mason MK, Hockman D, Curry L, Cunningham TJ, Duester G, Logan M, Jacobs DS, Illing N. Retinoic acid-independent expression of Meis2 during autopod patterning in the developing bat and mouse limb. Evodevo. 2015;6:6.

  78. 78.

    Tokita M. How the pterosaur got its wing. Biol Rev. 2015;90:1163–78.

  79. 79.

    Booker BM, Friedrich T, Mason MK, VanderMeer JE, Zhao J, Eckalbar WL, Logan M, Illing N, Pollard KS, Ahituv N. Bat accelerated regions identify a bat forelimb specific enhancer in the HoxD locus. PLoS Genet. 2016;12:e1005738.

  80. 80.

    Eckalbar WL, Schlebusch SA, Mason MK, Gill Z, Parker AV, Booker BM, Nishizaki S, Muswamba-Nday C, Terhune E, Nevonen KA, Makki N, Friedrich T, VanderMeer JE, Pollard KS, Carbone L, Wall JD, Illing N, Ahituv N. Transcriptomic and epigenomic characterization of the developing bat wing. Nat Genet. 2016;48:528–36.

  81. 81.

    Haeckel E. Kunstformen der Natur. 1904.

  82. 82.

    Dumont ER, Dávalos LM, Goldberg A, Santana SE, Rex K, Voigt CC. Morphological innovation, diversification and invasion of a new adaptive zone. Proc Biol Sci. 2011;279:1797–805.

  83. 83.

    Baker RJ, Bininda-Emonds ORP, Mantilla-Meluk H, Porter CA, Van Den Bussche RA. Molecular timescale of diversification of feeding strategy and morphology in New World Leaf-Nosed Bats (Phyllostomidae): A phylogenetic perspective. In Evolutionary history of bats: fossils, molecules and morphology. Edited by Gunnell GF. Simmons NB. Cambridge: Cambridge University Press. 2012;385–409.

  84. 84.

    Dumont ER. The effect of food hardness on feeding behavior in frugivorous bats (Phyllostomidae): an experimental study. J Zool Lond. 1999;248:219–29.

  85. 85.

    Dumont ER. Patterns of diversity in cranial shape among plant-visiting bats. Acta Chiropterol. 2004;6:59–74.

  86. 86.

    Santana SE, Grosse IR, Dumont ER. Dietary hardness, loading behavior, and the evolution of skull form in bats. Evolution. 2012;66:2587–98.

  87. 87.

    Giannini NP, Simmons NB. The Chiropteran Premaxilla: a reanalysis of morphological variation and its phylogenetic interpretation. Am Museum Novit. 2007;3585:1–44.

  88. 88.

    Orr DJA, Teeling EC, Puechmaille SJ, Finarelli JA. Patterns of orofacial clefting in the facial morphology of bats: a possible naturally occurring model of cleft palate. J Anat. 2016;229:657–72.

  89. 89.

    Phillips CD, Butler B, Fondon JW, Mantilla-Meluk H, Baker RJ. Contrasting evolutionary dynamics of the developmental regulator PAX9, among bats, with evidence for a novel post-transcriptional regulatory mechanism. PLoS ONE. 2013;8:e57649.

  90. 90.

    Shapiro MD, Bell MA, Kingsley DM. Parallel genetic origins of pelvic reduction in vertebrates. Proc Natl Acad Sci USA. 2006;103:13753–8.

  91. 91.

    Kohlsdorf T, Cummings MP, Lynch VJ, Stopper GF, Takahashi K, Wagner GP. A molecular footprint of limb loss: sequence variation of the autopodial identity gene Hoxa-13. J Mol Evol. 2008;67:581–93.

  92. 92.

    Woltering JM, Vonk FJ, Müller H, Bardine N, Tuduce IL, de Bakker MA, Knöchel W, Sirbu IO, Durston AJ, Richardson MK. Axial patterning in snakes and caecilians: evidence for an alternative interpretation of the Hox code. Dev Biol. 2009;332:82–9.

  93. 93.

    Mallarino R, Henegar C, Mirasierra M, Manceau M, Schradin C, Vallejo M, Beronja S, Barsh GS, Hoekstra HE. Developmental mechanisms of stripe patterns in rodents. Nature. 2016;539:518–23.

  94. 94.

    Sears KE. Differences in growth generate the diverse palate shapes of new world leaf-nosed Bats (Order Chiroptera, Family Phyllostomidae). Evol Biol. 2014;41:12–21.

  95. 95.

    Cretekos CJ, Weatherbee SD, Chen CH, Badwaik NK, Niswander L, Behringer RR, Rasweiler JJ 4th. Embryonic staging system for the short-tailed fruit bat, Carollia perspicillata, a model organism for the mammalian order Chiroptera, based upon timed pregnancies in captive-bred animals. Dev Dyn. 2005;233:721–38.

  96. 96.

    Giannini N, Goswami A, Sánchez-Villagra MR. Development of integumentary structures in Rousettus amplexicaudatus (Mammalia: Chiroptera: Pteropodidae) during late-embryonic and fetal stages. J Mammal. 2006;87:993–1001.

  97. 97.

    Tokita M. Normal embryonic development of the Japanese pipistrelle, Pipistrellus abramus. Zoology. 2006;109:137–47.

  98. 98.

    Hockman D, Mason MK, Jacobs DS, Illing N. The role of early development in mammalian limb diversification: a descriptive comparison of early limb development between the natal long-fingered bat (Miniopterus natalensis) and the mouse (Mus musculus). Dev Dyn. 2009;238:965–79.

  99. 99.

    Wang Z, Han N, Racey PA, Ru B, He G. A comparative study of prenatal development in Miniopterus schreibersii fuliginosus, Hipposideros armiger and H. pratti. BMC Dev Biol. 2010;10:10.

  100. 100.

    Paksuz EP, Hayretdag S, Olgun K. Prenatal development in greater mouse- eared bat, Myotis myotis (Borkhausen, 1797) (Chiroptera, Vespertilionidae). Anat Histol Embryol. 2017;46:563–71.

  101. 101.

    Chen MY, Liang D, Zhang P. Phylogenomic Resolution of the Phylogeny of Laurasiatherian Mammals: Exploring Phylogenetic Signals within Coding and Noncoding Sequences. Genome Biol Evol. 2017;9:1998–2012.

  102. 102.

    De Moraes DA, Hingst-Zaher E, Maarcus L, Cerqueria R. A geometric morphometric analysis of cranial and mandibular shape variation of didelphid marsupials. Hyatrix It J Mamm. 2000;11:115–30.

  103. 103.

    Goswami A, Polly PD, Mock OB, Sánchez-Villagra MR. Shape, variance and integration during craniogenesis: contrasting marsupial and placental mammals. J Evol Biol. 2012;25:862–72.

  104. 104.

    Prevosti FJ, Turazzini GF, Ercoli MD, Hingst-Zaher E. Mandible shape in marsupial and placental carnivorous mammals: a morphological comparative study using geometric morphometrics. Zool J Linn Soc. 2012;164:836–55.

  105. 105.

    Bennett CV, Goswami A. Statistical support for the hypothesis of developmental constraint in marsupial skull evolution. BMC Biol. 2013;11:52.

  106. 106.

    De Iuliis G, Bargo MS. Vizcaíno SF. Variation in skull morphology and mastication in the fossil giant armadillos Pampatherium spp. And allied genera (Mammalia: Xenarthra: Pampatheriidae), with comments on their systematics and distribution. J Verterbr Paleontol. 2000;20:743-54.

  107. 107.

    Hautier L, Billet G, Eastwood B, Lane J. Patterns of morphological variation of extant sloth skulls and their implication for future conservation efforts. Anat Rec. 2014;297:979–1008.

  108. 108.

    Hossotani CMS, Ragusa-Netto J, e Luna HS. Skull morphometry and vault sutures of Myrmecophaga tridactyla and Tamandua tetradactyla. Iheringia, Série Zoologia. 2017;107:e2017038.

  109. 109.

    Panchetti F, Scalici M, Carpaneto GM, Gibertini G. Shape and size variations in the cranium of elephant-shrews: a morphometric contribution to a phylogenetic debate. Zoomorphology. 2008;127:69–82.

  110. 110.

    Scalici M, Panchetti F. Morphological cranial diversity contributes to phylogeny in soft-furred sengis (Afrotheria, Macroscelidea). Zoology. 2011;114:85–94.

  111. 111.

    Finlay S, Cooper N. Morphological diversity in tenrecs (Afrosoricida, Tenrecidae): comparing tenrec skull diversity to their closest relatives. PeerJ. 2015;30(3):e927.

  112. 112.

    Paunović G, Bogićević K, Urošević A. Mandible shape differentiation between Mammuthus trogontherii and M. primigenius and mandible shape ontogeny in M. primigenius specimens from Serbia: A preliminary explorative geometric morphometric study. Quant Int. 2017;443:212-20.

  113. 113.

    Barros H, Meirelles ACO, Luna FO, Marmontel M, Cordeiro-Esterla P, Santos N, Astúa D. Cranial and chromosomal geographic variation in manatees (Mammalia: Sirenia: Trichechidae) with the description of the Antillean manatee karyotype in Brazil. J Zool Syst Ecol Res. 2017;55:73–87.

  114. 114.

    Rohlf FJ, Loy A, Corti M. Morphometric analysis of old world Talpidae (Mammalia, Insectivora) using Partial-warp scores. Syst Biol. 1996;45:334–62.

  115. 115.

    Loy A, Capula M, Palombi A, Capanna E. Genetic and morphometric evidence of introgression between two species of moles (Insectivora: Talpa europaea and Talpa romana) in central Italy. J Zool Lond. 2001;254:229–38.

  116. 116.

    Barrow E, Macleod N. Shape variation in the mole dentary (Talpidae: Mammalia). Zool J Linn Soc. 2008;153:187–221.

  117. 117.

    Ventura J, López-Fuster MJ. Geometric morphometrics of the mandible in the Iberian desman, Galemys pyrenaicus (Mammalia: Soricomorpha): Is there a significant variation in form during post-weaning life? Mamm Biol. 2010;75:191–7.

  118. 118.

    Shchipanov NA, Voyta LL, Bobretsov AV, Kuprianova IF. Intra-species structuring in the common shrew Sorex araneus (Lipotyphla: Soricidae) in European Russia: morphometric variability could give evidence of limitation of interpopulation migration. Russian J Theriol. 2014;13:119–40.

  119. 119.

    Shchipanova NA, Sychevaa VB, Tumasyan FA. Morphometric distances and population structuring in the common shrew Sorex araneus L. (Lipotyphla: Soricidae). Biology Bulletin. 2016;43:437–49.

  120. 120.

    Bornholdt R, Oliveira LR, Fabián ME. Size and shape variability in the skull of Myotis nigricans (Schinz, 1821) (Chiroptera: Vespertilionidae) from two geographic areas in Brazil. Braz J Biol. 2008;68:623–31.

  121. 121.

    Sorensen DW, Butkus C, Cooper LN, Cretekos CJ, Rasweiler JJ 4th, Sears KE. Palate variation and evolution in New World leaf-nosed and Old World fruit bats (Order Chiroptera). Evol Biol. 2014;41:595–605.

  122. 122.

    Zuccarelli MD. Comparative morphometric analysis of captive vs. wild African lion (Panthera leo) skulls. Bios. 2004;75:131–8.

  123. 123.

    Slater GJ, Van Valkenburgh B. Allometry and performance: the evolution of skull form and function in felids. J Evol Biol. 2009;22:2278–87.

  124. 124.

    Schutz H, Polly PD, Krieger JD, Guralnick RP. Differential sexual dimorphism: size and shape in the cranium and pelvis of grey foxes (Urocyon). Biol J Linn Soc. 2009;96:339–53.

  125. 125.

    Goswami A, Polly PD. The influence of modularity on cranial morphological disparity in Carnivora and Primates (Mammalia). PLoS ONE. 2010;5:e9517.

  126. 126.

    Drake AG, Klingenberg CP. Large-scale diversification of skull shape in domestic dogs: disparity and modularity. Am Nat. 2010;175:289–301.

  127. 127.

    Figueirido B1, Serrano-Alarcón FJ, Slater GJ, Palmqvist P. Shape at the cross-roads: homoplasy and history in the evolution of the carnivoran skull towards herbivory. J Evol Biol. 2010;23:2579-94.

  128. 128.

    Milenković M, Šipetić VJ, Blagojević J, Tatović S, Vujošević M. Skull variation in Dinaric-Balkan and Carpathian gray wolf populations revealed by geometric morphometric approaches. J Mammal. 2010;91:376–86.

  129. 129.

    Tanner JB, Zelditch ML, Lundrigan BL, Holekamp KE. Ontogenetic change in skull morphology and mechanical advantage in the spotted hyena (Crocuta crocuta). J Morphol. 2010;271:353–65.

  130. 130.

    Goswami A, Milne N, Wroe S. Biting through constraints: cranial morphology, disparity and convergence across living and fossil carnivorous mammals. Proc Biol Sci. 2011;278:1831–9.

  131. 131.

    Drake AG. Dispelling dog dogma: an investigation of heterochrony in dogs using 3D geometric morphometric analysis of skull shape. Evol Dev. 2011;13:204–13.

  132. 132.

    Figueirido B, Palmqvist P, Pérez-Claros JA, Dong W. Cranial shape transformation in the evolution of the giant panda (Ailuropoda melanoleuca). Naturwissenschaften. 2011;98:107–16.

  133. 133.

    Asahara M. Shape variation in the skull and lower carnassial in a wild population of raccoon dog (Nyctereutes procyonoides). Zoolog Sci. 2013;30:205–10.

  134. 134.

    Figueirido B, Tseng ZJ, Martín-Serra A. Skull shape evolution in durophagous carnivorans. Evolution. 2013;67:1975–93.

  135. 135.

    Schoenebeck JJ, Ostrander EA. The genetics of canine skull shape variation. Genetics. 2013;193:317–25.

  136. 136.

    Piras P, Maiorino L, Teresi L, Meloro C, Lucci F, Kotsakis T, Raia P. Bite of the cats: relationships between functional integration and mechanical performance as revealed by mandible geometry. Syst Biol. 2013;62:878–900.

  137. 137.

    Asahara M. Shape variation in the skull within and between wild populations of the raccoon dog (Nyctereutes procyonoides) in Japan. Mammal Study. 2014;39:105–13.

  138. 138.

    Curtis AA, Van Valkenburgh B. Beyond the sniffer: frontal sinuses in Carnivora. Anat Rec (Hoboken). 2014;297:2047–64.

  139. 139.

    Meachen JA, Janowicz AC, Avery JE, Sadleir RW. Ecological changes in Coyotes (Canis latrans) in response to the ice age megafaunal extinctions. PLoS ONE. 2014;9:e116041.

  140. 140.

    Jones KE, Smaers JB, Goswami A. Impact of the terrestrial-aquatic transition on disparity and rates of evolution in the carnivoran skull. BMC Evol Biol. 2015;15:8.

  141. 141.

    Tseng ZJ, Flynn JJ. Are cranial biomechanical simulation data linked to known diets in extant taxa? A method for applying diet-biomechanics linkage models to infer feeding capability of extinct species. PLoS ONE. 2015;10:e0124020.

  142. 142.

    Lau AC, Asahara M, Han SY, Kimura J. Sexual dimorphism of the Eurasian otter (Lutra lutra) in South Korea: Craniodental geometric morphometry. J Vet Med Sci. 2016;78:1007–11.

  143. 143.

    de Moura Bubadué J, Cáceres N, Dos Santos Carvalho R, Meloro C. Ecogeographical Variation in Skull Shape of South-American Canids: Abiotic or Biotic Processes? Evol Biol. 2016;43:145–59.

  144. 144.

    Bertolini F, Gandolfi B, Kim ES, Haase B, Lyons LA, Rothschild MF. Evidence of selection signatures that shape the Persian cat breed. Mamm Genome. 2016;27:144–55.

  145. 145.

    Lau AC, Asahara M, Han SY, Kimura J. Geographic variation of craniodental morphology of the Eurasian otter (Lutra lutra) in East Asia. J Vet Med Sci. 2017;79:144–52.

  146. 146.

    Geiger M, Evin A, Sánchez-Villagra MR, Gascho D, Mainini C, Zollikofer CPE. Neomorphosis and heterochrony of skull shape in dog domestication. Sci Rep. 2017;7:13443.

  147. 147.

    Yoshida H, Shirakihara K, Shirakihara M, Takemura A. Geographic variation in the skull morphology of the finless porpoise Neophocaena phocaenoides in Japan waters. Fish Sci. 1995;61:555–8.

  148. 148.

    Monteiro-Filho ELA, Monteiro LR, Reis SF. Skull shape and size divergence in dolphins of the genus Sotalia: a tridimensional morphometric analysis. J Mammal. 2002;83:125–34.

  149. 149.

    Amaral AR, Coelho MM, Marugán-Lobón J, Rohlf FJ. Cranial shape differentiation in three closely related delphinid cetacean species: Insights into evolutionary history. Zoology. 2009;112:38–47.

  150. 150.

    Gutstein CS, Cozzuol MA, Vargas AO, Suárez MA, Schultz CL, Rubilar-Rogers DR. Patterns of skull variation of Brachydelphis (Cetacea, Odontoceti) from the Neogene of the southeastern pacific. J Mammal. 2009;90:504–19.

  151. 151.

    Galatius A, Berta A, Frandsen MS, Goodall RNP. Interspecific variation of ontogeny and skull shape among porpoises (Phocoenidae). J Morphol. 2011;272:136–48.

  152. 152.

    Loy A, Tamburelli A, Carlini R, Slice DE. Craniometric variation of some Mediterranean and Atlantic populations of Stenella coeruleoalba (Mammalia, Delphinidae): A three-dimensional geometric morphometric analysis. Marine Mammal Science. 2011;27:65–78.

  153. 153.

    Parés-Casanova PM. Fabre – Anatomia L. Size and shape variability in the skull of the bottlenose dolphin, Tursiops truncatus (Montagu, 1821). Anat Histol Embryol. 2013;42:379–83.

  154. 154.

    Del Castillo DL, Flores DA, Cappozzo HL. Ontogenetic development and sexual dimorphism of franciscana dolphin skull: A 3D geometric morphometric approach. J Morphol. 2014;00:1–10.

  155. 155.

    Owen J, Dobney K, Evin A, Cucchi T, Larson G, Vidarsdottir US. The zooarchaelogical application of quantifying cranial shape differences in wild boar and domestic pig (Sus scrofa) using 3D geometric morphometrics. J Archaeol Sci. 2014;43:159–67.

  156. 156.

    Gol’din P, Vishnyakova K. Habitat shapes skull profile of small cetaceans: evidence from geographical variation in Black Sea harbour porpoises (Phocoena phocoena relicta). Zoomorphology. 2016;135:387-93.

  157. 157.

    Evin A, Owen J, Larson G, Debiais-Thibaud M, Cucchi T, Vidarsdottir US, Dobney K. A test for paedomorphism in domestic pig cranial morphology. Biol Lett. 2017;13:20170321.

  158. 158.

    Endo H, Rashdi ABM, Yamagiwa D, Yamada J. Principal component analysis of the common tree shrew skulls from Laos, Thailand and Malaysia. National Science Museum monographs. 2000;18:171–5.

  159. 159.

    Yusoff AM, Kumaran JV, Mohd Tahir NFD, Mokhtar SI. Preliminary study of skull polymorphims of tupaia glis in Peninsular Malaysia by using morphoj. Jurnal Teknologi (Sciences & Engineering). 2015;72:67-9.

  160. 160.

    Sargis EJ, Woodman N, Morningstar NC, Bell TN, Olson LE. Skeletal variation and taxonomic boundaries among mainland and island populations of the common treeshrew (Mammalia: Scandentia: Tupaiidae). Biol J Linn Soc. 2017;120:286–312.

  161. 161.

    Corti M, Di Giuliomaria C, Verheyen W. Three-dimensional geometric morphometrics of the African genus Lophuromys (Rodenta Muridae). Hystrix It J Mamm. 2000;11.

  162. 162.

    Corti M, Aguilera M, Capanna E. Size and shape changes in the skull accompanying speciation of South American spiny rats (Rodentia: Proechimys spp.). J Zool Lond. 2001;253:537-47.

  163. 163.

    Klingenberg CP, Leamy LJ, Routman EJ, Cheverud JM. Genetic architecture of mandible shape in mice: effects of quantitative trait loci analyzed by geometric morphometrics. Genetics. 2001;157:785–802.

  164. 164.

    Dobigny G, Baylac M, Denys C. Geometric morphometrics, neural networks and diagnosis of sibling Taterillus species (Rodentia, Gerbillinae). Biol J Linn Soc. 2002;77:319–27.

  165. 165.

    Klingenberg CP, Mebus K, Auffray JC. Developmental integration in a complex morphological structure: how distinct are the modules in the mouse mandible? Evol Dev. 2003;5:522–31.

  166. 166.

    Klingenberg CP1, Leamy LJ, Cheverud JM. Integration and modularity of quantitative trait locus effects on geometric shape in the mouse mandible. Genetics. 2004;166:1909-21.

  167. 167.

    Zelditch ML1, Lundrigan BL, Garland T Jr. Developmental regulation of skull morphology. I. Ontogenetic dynamics of variance. Evol Dev. 2004;6:194-206.

  168. 168.

    Kawakami M, Yamamura K. Cranial bone morphometric study among mouse strains. BMC Evol Biol. 2008;8:73.

  169. 169.

    Fornel R, Cordeiro-Estrela P, De Freitas TR. Skull shape and size variation in Ctenomys minutus (Rodentia: Ctenomyidae) in geographical, chromosomal polymorphism, and environmental contexts. Biol J Linn Soc. 2010;101:705–20.

  170. 170.

    Boell L, Tautz D. Micro-evolutionary divergence patterns of mandible shapes in wild house mouse (Mus musculus) populations. BMC Evol Biol. 2011;11:306.

  171. 171.

    Cox PG, Fagan MJ, Rayfield EJ, Jeffery N. Finite element modelling of squirrel, guinea pig and rat skulls: using geometric morphometrics to assess sensitivity. J Anat. 2011;219:696–709.

  172. 172.

    Jamniczky HA, Hallgrímsson B. Modularity in the skull and cranial vasculature of laboratory mice: implications for the evolution of complex phenotypes. Evol Dev. 2011;13:28–37.

  173. 173.

    Burgio G, Baylac M, Heyer E, Montagutelli X. Exploration of the genetic organization of morphological modularity on the mouse mandible using a set of interspecific recombinant congenic strains between C57BL /6 and mice of the Mus spretus species. Genes Genomes Genetics. 2012;2-1257-68.

  174. 174.

    Hautier L, Lebrun R, Cox PG. Patterns of covariation in the masticatory apparatus of hystricognathous rodents: implications for evolution and diversification. J Morphol. 2012;273:1319–37.

  175. 175.

    Paradis MR, Raj MT, Boughner JC. Jaw growth in the absence of teeth: the developmental morphology of edentulous mandibles using the p63 mouse mutant. Evol Dev. 2013;15:268–79.

  176. 176.

    Gonzalez PN, Kristensen E, Morck DW, Boyd S, Hallgrímsson B. Effects of growth hormone on the ontogenetic allometry of craniofacial bones. Evol Dev. 2013;15:133–45.

  177. 177.

    Christians JK, de Zwaan DR, Fung SH. Pregnancy associated plasma protein A2 (PAPP-A2) affects bone size and shape and contributes to natural variation in postnatal growth in mice. PLoS ONE. 2013;8:e56260.

  178. 178.

    Casanovas-Vilar I, van Dam J. Conservatism and adaptability during squirrel radiation: what is mandible shape telling us? PLoS ONE. 2013;8:e61298.

  179. 179.

    Stolz JF, Gonçalves GL, Leipnitz L, Freitas TR. DNA-based and geometric morphometric analysis to validate species designation: a case study of the subterranean rodent Ctenomys bicolor. Genet Mol Res. 2013;12:5023–37.

  180. 180.

    Anderson PS, Renaud S, Rayfield EJ. Adaptive plasticity in the mouse mandible. BMC Evol Biol. 2014;14:85.

  181. 181.

    Parsons TE, Weinberg SM, Khaksarfard K, Howie RN, Elsalanty M, Yu JC, Cray JJ Jr. Craniofacial shape variation in Twist1 +/ mutant mice. Anat Rec (Hoboken). 2014;297:826–33.

  182. 182.

    Pallares LF, Harr B, Turner LM, Tautz D. Use of a natural hybrid zone for genomewide association mapping of craniofacial traits in the house mouse. Mol Ecol. 2014;23:5756–70.

  183. 183.

    Lu X, Ge D, Xia L, Huang C, Yang Q. Geometric morphometric study of the skull shape diversification in Sciuridae (Mammalia, Rodentia). Integr Zool. 2014;9:231–45.

  184. 184.

    Singh N, Dutka T, Devenney BM, Kawasaki K, Reeves RH, Richtsmeier JT. Acute upregulation of hedgehog signaling in mice causes differential effects on cranial morphology. Dis Model Mech. 2015;8:271–9.

  185. 185.

    Lalis A, Evin A, Janier M, Koivogui L, Denys C. Host evolution in Mastomys natalensis (Rodentia: Muridae): An integrative approach using geometric morphometrics and genetics. Integr Zool. 2015;10:505–14.

  186. 186.

    Walsh RE, Aprígio Assis AP, Patton JL, Marroig G, Dawson TE, Lacey EA. Morphological and dietary responses of chipmunks to a century of climate change. Glob Chang Biol. 2016;22:3233–52.

  187. 187.

    Pallares LF, Turner LM, Tautz D. Craniofacial shape transition across the house mouse hybrid zone: implications for the genetic architecture and evolution of between-species differences. Dev Genes Evol. 2016;226:173–86.

  188. 188.

    Quintela FM, Fornel R, Freitas TR. Geographic variation in skull shape of the water rat Scapteromys tumidus (Cricetidae, Sigmodontinae): isolation-by-distance plus environmental and geographic barrier effects? An Acad Bras Cienc. 2016;88:451–66.

  189. 189.

    Percival CJ, Liberton DK, Pardo-Manuel de Villena F, Spritz R, Marcucio R, Hallgrímsson B. Genetics of murine craniofacial morphology: diallel analysis of the eight founders of the Collaborative Cross. J Anat. 2016;228:96–112.

  190. 190.

    Vasil’ev AG, Bolshakov VN, Evdokimov NG, Sineva NV. Morphological diversity of mole vole mono- and polymorphic populations: Does Chernov’s “compensation principle” work within a population? Dokl Biol Sci. 2016;468:118-21.

  191. 191.

    Pavličev M, Mitteroecker P, Gonzalez PM, Rolian C, Jamniczky H, Villena FP, Marcucio R, Spritz R, Hallgrimsson B. Development Shapes a Consistent Inbreeding Effect in Mouse Crania of Different Line Crosses. J Exp Zool B Mol Dev Evol. 2016;326:474–88.

  192. 192.

    Singh N, Albert FW, Plyusnina I, Trut L, Pӓӓbo S, Harvati K. Facial shape differences between rats selected for tame and aggressive behaviors. PLoS ONE. 2017;12:e0175043.

  193. 193.

    Abramov SA, Lopatina NV, Litvinov YN. Cranial size and shape variation in isolated populations of the Olkhon mountain vole (Alticola olchonensis Litvinov, 1960). Zoology (Jena). 2017;123:91–100.

  194. 194.

    Caumul R, Polly PD. Phylogenetic and environmental components of morphological variation: skull, mandible, and molar shape in marmots (Marmota, Rodentia). Evolution. 2005;59:2460–72.

  195. 195.

    Ge D, Lv X, Xia L, Huang C, Yang Q. Geometric morphometric of postnatal size and shape changes in the cranium of cape hare (Lagomorpha, Leporidea, Lepus capensis). Acta Thariologica Sinica. 2012;32:12–24.

  196. 196.

    Kraatz BP, Sherratt E, Bumacod N, Wedel MJ. Ecological correlates to cranial morphology in Leporids (Mammalia, Lagomorpha). PeerJ. 2015;3:e844.

  197. 197.

    Ge D, Yao L, Xia L, Zhang Z, Yang Q. Geometric morphometric analysis of skull morphology reveals loss of phylogenetic signal at the generic level in extant lagomorphs (Mammalia: Lagomorpha). Contributions to Zoology. 2015;84:267–84.

  198. 198.

    Kraatz B, Sherratt E. Evolutionary morphology of the rabbit skull. PeerJ. 2016;4:e2453.

  199. 199.

    De León MSP, Zollikofer CPE. Neanderthal cranial ontogeny and its implications for late hominid diversity. Nature. 2001;412:534–8.

  200. 200.

    Singleton M. Patterns of cranial shape variation in the Papionini (Primates: Cercopithecinae). J Hum Evol. 2002;42:547–78.

  201. 201.

    Hallgrímsson B, Willmore K, Dorval C, Cooper DM. Craniofacial variability and modularity in macaques and mice. J Exp Zool B Mol Dev Evol. 2004;302:207–25.

  202. 202.

    Perez SI, Bernal V, Gonzalez PN. Differences between sliding semi-landmark methods in geometric morphometrics, with an application to human craniofacial and dental variation. J Anat. 2006;208:769–84.

  203. 203.

    Cardini A, Elton S. Variation in guenon skulls (I): species divergence, ecological and genetic differences. J Hum Evol. 2008;54:615–37.

  204. 204.

    Shirai LT, Marroig G. Skull modularity in neotropical marsupials and monkeys: size variation and evolutionary constraint and flexibility. J Exp Zool B Mol Dev Evol. 2010;314:663–83.

  205. 205.

    Bennett CV, Goswami A. Morphometric analysis of cranial shape in fossil and recent euprimates. Anat Res Int. 2012;2012:478903.

  206. 206.

    Hayashi K, Saitoh S, Mizoguchi I. Morphological analysis of the skeletal remains of Japanese females from the Ikenohata-Shichikencho site. Eur J Orthod. 2012;34:575–81.

  207. 207.

    Robinson C. Geometric morphometric analysis of mandibular shape diversity in Pan. J Hum Evol. 2012;63:191–204.

  208. 208.

    Neaux D, Guy F, Gilissen E, Coudyzer W, Vignaud P, Ducrocq S. Facial orientation and facial shape in extant great apes: a geometric morphometric analysis of covariation. PLoS ONE. 2013;8:e57026.

  209. 209.

    Ito T, Nishimura T, Takai M. Ecogeographical and phylogenetic effects on craniofacial variation in macaques. Am J Phys Anthropol. 2014;154:27–41.

  210. 210.

    Pérez-Claros JA, Jiménez-Arenas JM, Palmqvist P. Neurocranium versus Face: A Morphometric Approach with Classical Anthropometric Variables for Characterizing Patterns of Cranial Integration in Extant Hominoids and Extinct Hominins. PLoS ONE. 2015;10:e0131055.

  211. 211.

    Neaux D, Gilissen E, Coudyzer W, Guy F. Integration between the face and the mandible of Pongo and the evolution of the craniofacial morphology of orangutans. Am J Phys Anthropol. 2015;158:475–86.

  212. 212.

    Fleagle JG, Gilbert CC, Baden AL. Comparing primate crania: The importance of fossils. Am J Phys Anthropol. 2016;161:259–75.

  213. 213.

    Martinez-Maza C, Freidline SE, Strauss A, Nieto-Diaz M. Bone Growth Dynamics of the Facial Skeleton and Mandible in Gorilla gorilla and Pan troglodytes. Evol Biol. 2016;43:60–80.

  214. 214.

    Schroeder L, Scott JE, Garvin HM, Laird MF, Dembo M, Radovčić D, Berger LR, de Ruiter DJ, Ackermann RR. Skull diversity in the Homo lineage and the relative position of Homo naledi. J Hum Evol. 2017;104:124–35.

  215. 215.

    Liu W, Sun X, Braut A, Mishina Y, Behringer RR, Mina M, Martin JF. Distinct functions for Bmp signaling in lip and palate fusion in mice. Development. 2005;132:1453–61.

  216. 216.

    Metzis V, Courtney AD, Kerr MC, Ferguson C, Rondón Galeano MC, Parton RG, Wainwright BJ, Wicking C. Patched1 is required in neural crest cells for the prevention of orofacial clefts. Hum Mol Genet. 2013;22:5026–35.

  217. 217.

    Jeong J, Mao J, Tenzen T, Kottmann AH, McMahon AP. Hedgehog signaling in the neural crest cells regulates the patterning and growth of facial primordia. Genes Dev. 2004;18:937–51.

  218. 218.

    Green RM, Feng W, Phang T, Fish JL, Li H, Spritz RA, Marcucio RS, Hooper J, Jamniczky H, Hallgrímsson B, Williams T. Tfap2a-dependent changes in mouse facial morphology result in clefting that can be ameliorated by a reduction in Fgf8 gene dosage. Dis Model Mech. 2015;8:31–43.

  219. 219.

    Niemann S, Zhao C, Pascu F, Stahl U, Aulepp U, Niswander L, Weber JL, Müller U. Homozygous WNT3 mutation causes tetra-amelia in a large consanguineous family. Am J Hum Genet. 2004;74:558–63.

  220. 220.

    Juriloff DM, Harris MJ, McMahon AP, Carroll TJ, Lidral AC. Wnt9b is the mutated gene involved in multifactorial nonsyndromic cleft lip with or without cleft palate in A/WySn mice, as confirmed by a genetic complementation test. Birth Defects Res A Clin Mol Teratol. 2006;76:574–9.

  221. 221.

    Wu W, Gu S, Sun C, He W, Xie X, Li X, Ye W, Qin C, Chen Y, Xiao J, Liu C. Altered FGF Signaling Pathways Impair Cell Proliferation and Elevation of Palate Shelves. PLoS ONE. 2015;10:e0136951.

  222. 222.

    Brugmann SA, Allen NC, James AW, Mekonnen Z, Madan E, Helms JA. A primary cilia-dependent etiology for midline facial disorders. Hum Mol Genet. 2010;19:1577–92.

  223. 223.

    Satokata I, Maas R. Msx1 deficient mice exhibit cleft palate and abnormalities of craniofacial and tooth development. Nat Genet. 1994;6:348–56.

  224. 224.

    Zhang Z, Song Y, Zhao X, Zhang X, Fermin C, Chen Y. Rescue of cleft palate in Msx1-deficient mice by transgenic Bmp4 reveals a network of BMP and Shh signaling in the regulation of mammalian palatogenesis. Development. 2002;129:4135–46.

  225. 225.

    Meester-Smoor MA, Vermeij M, van Helmond MJ, Molijn AC, van Wely KH, Hekman AC, Vermey-Keers C, Riegman PH, Zwarthoff EC. Targeted disruption of the Mn1 oncogene results in severe defects in development of membranous bones of the cranial skeleton. Mol Cell Biol. 2005;25:4229–36.

  226. 226.

    Yu L, Gu S, Alappat S, Song Y, Yan M, Zhang X, Zhang G, Jiang Y, Zhang Z, Zhang Y, Chen Y. Shox2-deficient mice exhibit a rare type of incomplete clefting of the secondary palate. Development. 2005;132:4397–406.

  227. 227.

    Ito Y, Yeo JY, Chytil A, Han J, Bringas P Jr, Nakajima A, Shuler CF, Moses HL, Chai Y. Conditional inactivation of Tgfbr2 in cranial neural crest causes cleft palate and calvaria defects. Development. 2003;130:5269–80.

  228. 228.

    Lack JB, Van Den Bussche RA. Identifying the confounding factors in resolving phylogenetic relationships in Vespertilionidae. J Mammal. 2010;91:1435–48.

Download references

Authors’ contributions

All authors contributed to the writing and editing of the manuscript. Both authors read and approved the final manuscript.

Acknowledgements

We thank Dr. Shin-ichiro Kawada (National Museum of Nature and Science, Tokyo) and Dr. Hiroko Kudo-Hirotani (Kanagawa Prefectural Museum of Natural History) for arrangement of the stuffed specimens of mammals, Mr. Jason Preble for proofreading the early version of the manuscript, and all members of the laboratory of Vertebrate Diversity and Evolution, Toho University for helpful discussions.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

Not applicable.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Correspondence to Masayoshi Tokita.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Mammals
  • Craniofacial morphology
  • Diversity
  • Transgenic mice
  • Bats
  • Facial processes
  • Neural crest
  • Ectomesenchyme
  • Bone
  • Orofacial cleft