Oldest skeleton of a fossil flying squirrel casts new light on the phylogeny of the group
Abstract
Flying squirrels are the only group of gliding mammals with a remarkable diversity and wide geographical range. However, their evolutionary story is not well known. Thus far, identification of extinct flying squirrels has been exclusively based on dental features, which, contrary to certain postcranial characters, are not unique to them. Therefore, fossils attributed to this clade may indeed belong to other squirrel groups. Here we report the oldest fossil skeleton of a flying squirrel (11.6 Ma) that displays the gliding-related diagnostic features shared by extant forms and allows for a recalibration of the divergence time between tree and flying squirrels. Our phylogenetic analyses combining morphological and molecular data generally support older dates than previous molecular estimates (~23 Ma), being congruent with the inclusion of some of the earliest fossils (~36 Ma) into this clade. They also show that flying squirrels experienced little morphological change for almost 12 million years.
https://doi.org/10.7554/eLife.39270.001eLife digest
Mammals can walk, hop, swim and fly; a few, like marsupial sugar gliders or colugos, can even glide. With 52 species scattered across the Northern hemisphere, flying squirrels are by far the most successful group that adopted this way of going airborne. To drift from tree to tree, these small animals pack their own ‘parachute’: a membrane draping between their lower limbs and the long cartilage rods that extend from their wrists. Tiny specialized wrist bones, which are unique to flying squirrels, help to support the cartilaginous extensions.
The origin of flying squirrels is a point of contention: while most genetic studies point towards the group splitting from tree squirrels about 23 million years ago, the oldest remains – mostly cheek teeth – suggest the animals were already soaring through forests 36 million years ago. However, recent studies show that the dental features used to distinguish between gliding and non-gliding squirrels may actually be shared by the two groups.
In 2002, the digging of a dump site in Barcelona unearthed a peculiar skeleton: first a tail and two thigh bones, big enough that the researchers thought it could be the fossil of a small primate. In fact, and much to the disappointment of paleoprimatologists, further excavating revealed that it was a rodent. As the specimen – nearly an entire skeleton – was being prepared, paleontologists insisted that all the ‘dirt’ attached to the bones had to be carefully screen-washed. From the mud emerged the minuscule specialized wrist bones: the primate-turned-rodent was in fact Miopetaurista neogrivensis, an extinct flying squirrel.
Here, Casanovas-Vilar et al. describe the 11.6 million years old fossil, the oldest ever found. The wrist bones reveal that the animal belongs to the group of flying squirrels that have large sizes. Evolutionary analyses that combined molecular and paleontological data demonstrated that flying squirrels evolved from tree squirrels as far back as 31 to 25 million years ago, and possibly even earlier. In addition, the results show that Miopetaurista is closely related to Petaurista, a modern group of giant flying squirrels. In fact, their skeletons are so similar that the large species that currently inhabit the tropical and subtropical forests of Asia could be considered living fossils.
Molecular and paleontological data are often at odds, but this fossil shows that they can be reconciled and combined to retrace history. Discovering older fossils, or even transitional forms, could help to retrace how flying squirrels took a leap from the rest of their evolutionary tree.
https://doi.org/10.7554/eLife.39270.002Introduction
Flying squirrels (Sciurinae, Pteromyini) are the only group of gliding mammals to have achieved a significant diversity (52 species in 15 genera) and wide geographical distribution across Eurasia and North America (Koprowski et al., 2016). They have been classically regarded as a distinct subfamily among the Sciuridae (McKenna and Bell, 1997; McLaughlin, 1984; Simpson, 1945), and even sometimes considered a separate family derived from a different group than the remaining sciurids (De Bruijn and Ünay, 1989; Forsyth Major, 1893; Mein, 1970). The fact that presumed fossil flying squirrels are at least as old as (or maybe even older than) the oldest tree squirrels (36.6 – 35.8 Ma) may support the latter hypothesis. However, flying squirrels are currently recognized as a monophyletic clade, as supported by a set of synapomorphies in the wrist (Thorington, 1984). The carpal anatomy of flying squirrels is unique, being related to the structures that support the patagium and their particular gliding position, which is different from that of all other gliding mammals (Thorington, 1984; Thorington and Darrow, 2000). Molecular phylogenies indicate that flying squirrels (tribe Pteromyini) are nested within tree squirrels (subfamily Sciurinae) and likely diverged as recently as the latest Oligocene–early Miocene (23 ± 2.1 Ma) (Fabre et al., 2012; Mercer and Roth, 2003; Steppan et al., 2004). Notwithstanding, the pteromyin fossil record suggests a much older split. Indeed, one of the earliest sciurids, Hesperopetes thoringtoni from the late Eocene (36.6 – 35.8 Ma) of North America, has been related to the lineage leading to flying squirrels according to dental morphology (Emry and Korth, 2007). In the light of molecular results, it was conceded that Hesperopetes unlikely represented a pteromyin and was not assigned to any squirrel subfamily (Emry and Korth, 2007). On the other hand, this genus appears to have been closely related to Oligopetes (Emry and Korth, 2007), an earliest Oligocene (ca. 34 – 31 Ma) purported flying squirrel from Europe and Pakistan (Cuenca Bescós and Canudo, 1992; De Bruijn and Ünay, 1989; Heissig, 1979; Marivaux and Welcomme, 2003). Hesperopetes is last recorded during the earliest Oligocene (Orellan; Korth, 2017), coinciding with the oldest record of Sciurion (Bell, 2004), yet another alleged flying squirrel. Isolated cheek teeth are the only material available for all these taxa, which have been related to flying squirrels exclusively based on dental morphology. In fact, the whole fossil record of flying squirrels almost exclusively consists of isolated cheek teeth and a few mandibular and maxillary fragments. Unfortunately, dental features commonly used to recognize flying squirrels are not unique but also present in other sciurids (Thorington et al., 2005), so it is uncertain if any of the extinct ‘flying’ squirrels belonged to this group. Furthermore, if any of the oldest (late Eocene–early Oligocene) forms truly represented a pteromyin this would imply a discrepancy of more than 10 Myr between molecular and paleontological data. Contrary to dental material, postcranial remains do show diagnostic characters of the pteromyins (Thorington, 1984; Thorington et al., 2005; Thorington and Darrow, 2000). Therefore, they are of utmost importance to clarify the assignment of extinct ‘flying’ squirrels and calibrate their divergence date from other sciurids. Yet, these have not been described and are rarely preserved in the fossil record. Here we report a remarkably complete skeleton of a Miocene squirrel that displays the gliding-related diagnostic features shared by extant pteromyins and allows for a recalibration of the time of origin and diversification of the group. The fossil record of ‘flying’ squirrels is further discussed in the light of this new finding and the results of our phylogenetic analyses.
Results
Recovered material and specific attribution
The described partial skeleton (IPS56468; Figure 1, Videos 1, 3D model in Supplementary file 1) was recovered at Abocador de Can Mata site ACM/C5-D1 (els Hostalets de Pierola, Catalonia, Spain; see Materials and Methods), with an estimated age of 11.63 Ma (Alba et al., 2017). The recovered remains were found partly articulated (Figure 2—figure supplement 1) and comprise more than 80 complete and fragmentary bones including the skull (Figure 6—figure supplement 1) and elements of the fore- and hindlimbs (Figures 2–5, Figure 2, Table 1). Additional material, including a second cranium (Figure 6—figure supplement 2), has been recovered from the same horizon and other roughly coeval ACM localities (Table 2). The specimens are assigned to Miopetaurista neogrivensis based on diagnostic cheek tooth morphology (Figure 3; for detailed description and comparisons of cheek teeth morphology see Appendix 3.1). In the ACM localities a second genus of ‘flying’ squirrel, Albanensia, is recorded, but Miopetaurista is clearly distinguished by its larger size, and several morphological features. The diagnostic characters of M. neogrivensis comprise: its large size; the presence of a complete entolophid and the frequent occurrence of a short mesolophid in the lower molars; and the large mesostyle in the P4 (Casanovas-Vilar et al., 2015; Mein, 1970). Miopetaurista neogrivensis has only been reported from La Grive L5 (type locality) and L3 in France, from Bellestar (Seu d’Urgell Basin, also in Catalonia), and from several sites from the Vallès-Penedès Basin (Casanovas-Vilar et al., 2015). This species is extremely rare, being represented by just a few isolated cheek teeth in most of the Vallès-Penedès sites.
Extant and fossil flying squirrels have been classified into different groups according to the complexity of dental morphology (Mein, 1970). The cheek teeth of Miopetaurista show a simple occlusal pattern, with enamel wrinkling only in the lower molars and no additional lophules (Figure 3, Appendix 3.1). This pattern clearly differs from the more complex one of other large-sized flying squirrels, such as Aeretes and Petaurista (Mein, 1970; Thorington et al., 2002). Therefore, Miopetaurista has been included within the group that comprises Aeromys and the small-sized flying squirrels, which do show simple dental patterns (Mein, 1970). However, our phylogenetic analyses (see below) show that M. neogrivensis is the sister taxon of extant Petaurista, a genus that would belong to a completely different group according to dental classification (Mein, 1970). Considering dental morphology Petaurista is assigned to a group characterized by its complex dental pattern with additional transverse lophules which would also comprise the genera Aeretes, Belomys, Eupetaurus and Trogopterus, among others (Mein, 1970). This clearly illustrates that dental characters, although useful to diagnose the different species and genera, should not warrant high consideration for disentangling the phylogenetic relationships between flying squirrels.
Morphological description and comparisons
Among the recovered postcranial material, a complete scapholunate and the dorsal end of the pisiform (Figure 4 and Video 2) are the most diagnostic elements of pteromyins, because they form the functional complex associated with the extension of the gliding membrane (Thorington, 1984; Thorington et al., 2005; Thorington et al., 2002; Thorington and Darrow, 2000). The styliform cartilage, which supports the patagium in all members of the group, attaches to the pisiform and is extended when the wrist is radially abducted and dorsiflexed (Thorington, 1984; Thorington and Darrow, 2000). The pisiform of M. neogrivensis displays an elevated process for the articulation with the scapholunate (Figure 4 and Video 2). This is characteristic of pteromyins, serving as a stabilizer of the styliform cartilage, whereas in other squirrels this bone articulates only with the triquetrum and the distal end of the ulna. Moreover, in the scapholunate of M. neogrivensis, the articular surface for the radius is much more convex than in tree squirrels, thus resembling the flying squirrel condition, which enables a greater radial abduction. Therefore, the proximal wrist joint morphology of M. neogrivensis indicates that this species belongs to the pteromyin clade and provides the oldest evidence of gliding locomotion in sciurids (see also Appendix 3.3). The latter is further confirmed by other postcranial adaptations shared with extant pteromyins (Thorington et al., 2005), including the elongated and slender limb bones with reduced muscular attachments (Figure 5 and Figure 5—figure supplement 1), which enhance joint extension during gliding (Thorington et al., 2005), as well as the elongated lumbar vertebrae (Figure 2; see Appendix 3.3 for a detailed description and comparisons of the postcranial elements). The elongation of lumbar vertebrae and limbs determines the size and shape of the patagium and dictates important aerodynamic features, such as the decreased wing loading of flying squirrels (Thorington and Heaney, 1981; Thorington and Santana, 2007).
Based on morphological (Thorington et al., 2002) and molecular data (Mercer and Roth, 2003), flying squirrels are divided into two distinct subtribes: Pteromyina, comprising large-sized forms, and Glaucomyina, for the small-sized ones (Thorington and Hoffmann, 2005). The skeleton of M. neogrivensis morphologically resembles that of the Pteromyina, further being comparable in size to their largest representatives. Body mass was estimated by means of an allometric regression of body mass vs. skull length in extant sciurids (see Materials and Methods), resulting in 1339 g (50% confidence intervals 1116 – 1606 g thus being in the range of most species of the extant giant flying squirrel Petaurista (about 1200 – 2000 g; Thorington et al., 2012). The long bones are almost indistinguishable of Petaurista. The skull, which was virtually reconstructed from two well-preserved specimens (Figure 6, Figure 6—figure supplement 1-2, Video 3, Table 3), is strikingly similar in size and morphology to that of the other large-sized flying squirrels, particularly Aeromys and Petaurista (for a detailed morphological description of the skull and comparisons see Appendix 3.2). These genera are characterized by their short and wide rostrum, moderately inflated bullae and relatively wide posterior region of the skull. Other morphological similarities include the robust and long postorbital process that partially encloses the orbit, the well-developed jugal process in the zygomatic arch and the presence of two septa in the tympanic cavity (Video 4; Appendix 3.2). Most of the smaller flying squirrels show more elongate muzzles, slender or shorter postorbital processes and, in some cases, a higher number of transbullar septa. The proximal carpal bones of M. neogrivensis not only unambiguously indicate that it is a flying squirrel, but also allow assigning it to the Pteromyina (Figure 4 and Video 2). The pisiform displays a distinct spur (triquetral process) that fits between the palmar surfaces of the scapholunate and the triquetrum. This process is completely lacking in the Glaucomyina (Thorington et al., 2002; Thorington and Darrow, 2000) (Figure 4 and Video 2; see also Appendix 3.3). Both subtribes are also distinguished by the origin of the tibiocarpalis muscle, which runs from the ankle to the tip of the styliform cartilage, defining the edge of the patagium. In the Glaucomyina the tibiocarpalis originates from a tuberosity on the distal tibia which is lacking in M. neogrivensis and the Pteromyina (Thorington et al., 2002). In the latter, the tibiocarpalis originates from the metatarsals instead.
Phylogenetic results
The assignment of M. neogrivensis to the subtribe Pteromyina is further confirmed by a total evidence phylogenetic analysis combining morphological and molecular data (see Materials and Methods). This analysis relied on 35 extant species of sciurids plus Aplodontia used as outgroup (both with molecular and morphological data) and two fossils (M. neogrivensis and the oldest-known tree squirrel, Douglassciurus jeffersoni, both represented only by morphological data; see also Appendices 1.1 and 2). Our analysis strongly supports M. neogrivensis as the sister taxon of the Petaurista spp. clade and indicates a divergence date between flying and tree squirrels ranging from the late Eocene to the late Oligocene (95% highest posterior density [HPD] interval 36.5 – 24.9 Ma; Figure 7). Such range is congruent with previous molecular estimates (Fabre et al., 2012; Mercer and Roth, 2003) but also supports older dates, as old as the oldest records of purported pteromyins (Emry and Korth, 2007). The Pteromyina and Glaucomyina would have diverged between the late Oligocene–early Miocene (HPD interval 27.1 – 18.1 Ma; Figure 7). We independently tested the estimates of Pteromyini/Sciurini divergence and the onset of Pteromyini crown diversification by means of a node dating analysis (see Materials and Methods) of the extant Sciurinae using two different calibration points, one for each tribe (Figure 7—figure supplement 1, Table 4; see also Appendix 1.2). Estimates for many nodes are somewhat younger using this alternative approach, but mostly overlap with the younger half of HPD intervals of the total evidence analysis (Figure 7). The age range for the Sciurini/Pteromyini split spans from the mid Oligocene to early Miocene (HPD interval 30.6 – 17.4 Ma; Figure 7—figure supplement 1). Therefore, divergence estimates derived using independent methods overlap for the late Oligocene.
Discussion
Divergence date between flying and tree squirrels
The two independent phylogenetic analyses estimate divergence dates ranging from 36.5 to 17.4 Ma (late Eocene–early Miocene; Figure 7 and Figure 7—figure supplement 1). These broad ranges are congruent with a previous molecular phylogenetic analysis of the sciurids (Mercer and Roth, 2003), which placed the split between tree and flying squirrels near the Oligocene/Miocene boundary (23 ± 2.1 Ma), although allowing for substantially older dates. Total evidence analysis in particular provides older estimates (36.5 – 24.9 Ma) than node dating (30.6 – 17.4 Ma), but these still marginally overlap with those of previous molecular results. The different results may arise from the different selection of calibration points (see Table 4). In that study (Mercer and Roth, 2003) the age of Douglassciurus (ca. 36 Ma) was assigned to the base of the sciurid crown radiation (i.e. the origin of extant major clades) to calibrate the phylogenetic tree whereas the root of our phylogenetic trees is calibrated using multiple points (see Materials and Methods).
A few molecular phylogenetic studies (Montgelard et al., 2008; Tapaltsyan et al., 2015), generally dealing with the whole rodent order and including only a few sciurid genera, have provided older ages for the divergence of flying squirrels dating back to the late Eocene and earliest Oligocene (34.5 – 30.9 Ma). Again, this may be attributed to the different selection of calibration points, which in most cases consider paleontological data from other rodent groups. Finally, rodent diversification has been analyzed using a molecular supermatrix that included 98% of extant squirrel genera (Fabre et al., 2012). This study found a late Oligocene divergence date between pteromyins and sciurins, which is perfectly congruent with our results.
Total evidence 95% highest posterior density (HPD) interval for the pteromyin divergence is very broad (36.5–24.9 Ma; Figure 7), being consistent with the range of Oligopetes (Cuenca Bescós and Canudo, 1992; De Bruijn and Ünay, 1989; Heissig, 1979; Marivaux and Welcomme, 2003), the oldest records of Sciurion (Bell, 2004) and even the earliest (late Eocene) occurrences of Hesperopetes (Korth, 2017) (Figure 8). In contrast, node dating analysis estimates a younger HPD interval that would exclude the older records of these genera (30.6 – 17.4 Ma; Figure 7—figure supplement 1). According to our results, late Eocene to early Oligocene alleged flying squirrels might indeed belong to this group, but this is only supported by total evidence analysis (Figure 7). Therefore, our analyses are not conclusive to this regard and further fossil data are required to elucidate the phylogenetic position of these older ‘flying’ squirrels. Although allowing for older and younger ages, the independently derived estimates are generally older than previous molecular results (Mercer and Roth, 2003) and overlap for the late Oligocene, which should be considered the most likely time for pteromyin divergence.
Flying squirrel radiation and biogeography in relation to Miocene climatic changes
Except for the partial skeleton of Miopetaurista neogrivensis described here, no diagnostic postcranial material has been recovered for any other extinct ‘flying’ squirrel, so their assignment to the pteromyins is doubtful. Here we discuss the fossil record of purported pteromyins in the light of our phylogenetic results, but in these instances our conclusions depend of the correctness of the tribe assignment.
Our analyses indicate that the initial diversification of flying squirrels into the two extant subtribes (Glaucomyina and Pteromyina) occurred between the latest Oligocene and the middle Miocene, with estimates derived from both methods overlapping for the early Miocene (Figure 7 and Figure 7—figure supplement 1). Such time interval coincides with the almost simultaneous earliest records of new genera of ‘flying’ squirrels in both North America and Europe. In Europe, these include the small-sized genus Blackia, which is first recorded near the Oligocene/Miocene transition (biozone MP30) (Engesser and Storch, 2008) at around 23.3 – 23.0 Ma (Figure 8). According to dental morphology this genus appears to be closely related to the older North American Sciurion (Skwara, 1986), so in case their pteromyin affinities were confirmed this would argue for a North American origin of the group. Besides Blackia three additional genera are recorded in the European early Miocene (De Bruijn, 1999): Aliveria, Miopetaurista and Neopetes (Figure 8). These occurrences date back to biozone MN3, corresponding to the early Miocene, yielding a relative age of 20.0 – 18.0 Ma. A species of the extant genus Hylopetes has been erected based on material from site of Oberdorf (Austria), correlated to the MN4 (18.0 – 17.0 Ma) (De Bruijn, 1998). Such occurrence would pull back the range of this genus into the early Miocene and has been taken for granted in some recent molecular studies (Lu et al., 2013) even though it blatantly disagrees with previous molecular results (Mercer and Roth, 2003). Indeed, some paleontologists have argued that the characters justifying the ascription of the Austrian and other material to Hylopetes are symplesiomorphies shared with many other flying squirrel genera (e.g. Petinomys, Glaucomys) and assign this material to the extinct ‘flying’ squirrel genera Neopetes and Pliopetes (Casanovas-Vilar et al., 2015; Daxner-Höck, 2004). In North America, aside from Sciurion the larger-sized genus Petauristodon is also present from the late Oligocene to the late Miocene (Goodwin, 2008) (Figure 8). While Sciurion is already known since the early Oligocene (33.5-33.0 Ma) (Bell, 2004), the oldest record of Petauristodon is a single molar dated between 25.9 and 23.8 Ma from the John Day Formation of Oregon (Korth and Samuels, 2015). Subsequent records date back to the early Miocene (ca. 19 Ma; Goodwin, 2008). Interestingly, the specimens ascribed to Sciurion and Petauristodon were previously referred to the European genera Blackia and Miopetaurista, respectively (Goodwin, 2008). Finally, in Asia the genera Meinia (Qiu, 1981), Parapetaurista (Qiu, 1981) and Shuanggouia (Qiu and Liu, 1986) have been reported from the early Miocene (Shanwangian) Shanwang and Xiacaowan formations of North and Eastern China, with an estimated age of ca. 18 – 16 Ma or slightly younger (Qiu et al., 2013b) (Figure 8). Even though fossil evidence can only justify the ascription of Miopetaurista to the pteromyins, it is worth noting that Europe records highest ‘flying’ squirrel diversity at the time (Figure 8), which has led to the suggestion that the group may well have originated there and immediately dispersed into Asia and North America near the Oligocene/Miocene transition (Lu et al., 2013). Notwithstanding, this hypothesis is challenged by the fossil record of ‘flying’ squirrels since the oldest occurrences are in North America. If certain Oligocene forms can be ultimately assigned to the pteromyins this would confirm an opposed model with an early origin and initial diversification in North America and a later dispersal into Eurasia by the latest Oligocene.
Initial flying squirrel diversification and dispersal across the Northern Hemisphere coincided with a period of high mean global temperatures during the early Miocene that peaked between 17 to 15 Ma, during the so-called Mid-Miocene Climatic Optimum (Zachos et al., 2001) (Figure 8). Humid warm-temperate broadleaf and mixed forests, resembling those existing in the southeastern coast of Asia, characterized the mid latitudes in Eurasia and North America (Pound et al., 2012). These forests provided a suitable habitat for flying squirrels and would have contributed to their initial radiation and dispersal into different continents. Our phylogenetic analyses show that some extant genera, such as Petaurista, Pteromys and Glaucomys, diverged approximately at that time (Figure 7 and Figure 7—figure supplement 1), thus agreeing with previous molecular results (Arbogast et al., 2017; Mercer and Roth, 2003).
The warm early Miocene phase was followed by a gradual cooling and the reestablishment of permanent major ice sheets on Antarctica by about 10 Ma (Zachos et al., 2001). Warm-temperate forests were still dominant in the mid latitudes during the middle and the beginning of the late Miocene. However, their distribution became more restricted, particularly in Western North America and Central Asia, where they began to be replaced by cooler and drier biomes (Pound et al., 2012). ‘Flying’ squirrels are particularly diverse in Europe during this interval, with as many as five different species co-occurring in a single site (Casanovas-Vilar et al., 2015). In addition, certain genera such as Miopetaurista, Pliopetaurista and Albanensia, are widely distributed (Casanovas-Vilar et al., 2015) (Figure 8). Miopetaurista and Pliopetaurista are also know from the beginning of the late Miocene of Amuwusu (Inner Mongolia, China) (Qiu et al., 2013b), whereas Albanensia or a closely related form would occur in the latest Miocene of Shihuiba (Yunnan, China) (Qiu, 2002). Our phylogenetic analyses show that most extant Asian flying squirrel genera diverged during the interval from 15 to 10 Ma (Figure 7 and Figure 7—figure supplement 1), thus coinciding with this time of high diversity and geographic range extension. On the contrary, ‘flying’ squirrel diversity in North America stayed at very low levels, with only the genera Petauristodon and Sciurion known from the middle and the beginning of the late Miocene (Goodwin, 2008) (Figure 8). Open habitats (mixed scrubland-grassland) housing a remarkable diversity of grazing mammals became widespread in central North America at that time, even though C4-dominated grasslands did not spread until the late Miocene (ca. 7 Ma) (Strömberg, 2011). These unsuitable habitats would have hampered the radiation of flying squirrels there.
Geographical range contraction and origins of extant flying squirrel genera
Throughout the later Miocene, warm-temperate forests continued reducing their extension and intermingled with drier and cooler biomes in Eurasia (Pound et al., 2012). In the Mediterranean regions, open habitats corresponding to woodlands and scrublands, already occurred in the southern Iberian Peninsula and Turkey since the middle Miocene (Fauquette et al., 2007; Strömberg, 2011), but would increase their extension at that time. Furthermore, the characteristic Mediterranean rainfall seasonality (summer drought) appeared during the Pliocene (Suc, 1984). In Central Europe, deciduous forests increasingly replaced warm-temperate mixed ones, which became restricted to coastal areas of the Mediterranean and the Paratethys (Mosbrugger et al., 2005). In Asia, major physiographical changes such as the uplift of the Himalaya-Tibetan Plateau affected atmospheric circulation and resulted in increased aridification of large areas of the continent (Zhisheng et al., 2001). Carbon stable isotope records of soil carbonates and ungulate tooth enamel indicate that C4 grasses rise to dominance in ecosystems between 8 – 7 Ma in Pakistan (Cerling et al., 1993; Strömberg, 2011) and later in China (Strömberg, 2011). All these environmental changes had profound effects in flying squirrel diversity and biogeography. The group became increasingly rarer over much of its former range and several genera (Forsythia, Albanensia, Miopetaurista, Blackia) disappeared during the late Miocene and the Pliocene (De Bruijn, 1999) (Figure 8). In Asia, fossil occurrences during the late Miocene are mostly confined to the southeast (Yunnan province in China, Thailand) (Lu et al., 2013). In North America, only the species Miopetaurista webbi is known by very scarce remains from the latest Miocene and the Pliocene of Florida (Goodwin, 2008); Robertson, 1976), being the only record of this genus outside Eurasia (Figure 8). It is not surprising that this last record of a large-sized flying squirrel in North America comes from Florida, an area that is still characterized by a humid subtropical climate with abundant densely forested areas. The extant American flying squirrel genus Glaucomys, which today inhabits temperate deciduous forests and boreal coniferous forests, is recorded for the first time already in the Pleistocene (Ruez, 2001) although it would have diverged significantly earlier according to molecular results (Figure 7 and Figure 7—figure supplement 1; Arbogast et al., 2017; Mercer and Roth, 2003).
The Pleistocene records of flying squirrels mostly correspond to extant genera and species (Lu et al., 2013) (Figure 8). This is again congruent with our results, which show that most extant species had already diverged during the Pliocene and some even at the latest Miocene (Figure 7 and Figure 7—figure supplement 1). During the Pleistocene, Glaucomys is the only flying squirrel known from North America. In Asia, the genera Aretes, Belomys, Petaurista, Pteromys and Trogopterus have been recorded from several sites of south and eastern China (Lu et al., 2013). In addition, Belomys has also been reported from Thailand and Pteromys from Japan (Lu et al., 2013). The Pleistocene records of all these genera are located within their current geographical range and generally correspond to extant species. In the European Pleistocene, flying squirrels are represented by the genera Neopetes, Petauria and Pliopetaurista (Figure 8). Quite surprisingly, there are no fossil records of Pteromys, the only flying squirrel genus still extant in Europe. Neopetes comprises species formerly included in the extant genus Hylopetes (Jackson and Thorington, 2012). Petauria is a large-sized flying squirrel known from middle Pleistocene fissure fillings and cave deposits of Germany and Poland. The available cheek teeth show several striking morphological similarities with the extant Petaurista, including: smooth enamel with numerous lophules, particularly in the basin of the lower cheek teeth; prominent mesostylid in the lower molars; and presence of a well-developed postero-lingual re-entrant fold in the upper molars. We agree with the opinion of some authors that this genus is a junior subjective synonym of Petaurista (Jackson and Thorington, 2012; Thorington et al., 2005). Therefore, the geographical range of Petaurista during the middle Pleistocene included Europe. Whether the genus originated in Europe or Asia cannot be resolved, since its closest relative, Miopetaurista, shows a similarly broad geographical range. It is worth noting that Miopetaurista is remarkably similar to Petaurista, to the point that their postcranial skeleton is virtually identical, even in specificities such as the more reduced lateral epicondylar ridge of the humerus or the wider patellar groove in the femur as compared to other (generally smaller) flying squirrels (for more detailed comparisons see Appendix 3.3). Cranial morphology evidences a close affinity between both genera, with details such as the development of the postorbital processes or the short and wide rostrum, being surprisingly similar. The only differences lie in cheek tooth morphology. Miopetaurista shows a relatively simple morphology, with faint enamel in the upper cheek teeth and no additional longitudinal lophules. In contrast, Petaurista presents relatively higher-crowned teeth and a much more complex morphology with additional longitudinal lophules (Mein, 1970; Thorington et al., 2002). In addition, the upper cheek teeth show a well-defined distolingual flexus, a feature that only occurs in the M3 in Miopetaurista. Dental morphological differences apart, it is worth noting that large-sized flying squirrels must be regarded as a very conservative group, having experienced little morphological changes since the late middle Miocene.
Concerning the genus Petaurista, our analyses recognize the three main species groups (Figure 7 and Figure 7—figure supplement 1) that had already been recognized by previous molecular analyses (Li et al., 2013). The first group solely includes Petaurista leucogenys, a species endemic to Japan. It is found to have diverged during the late Miocene (ca. 11 – 6 Ma), significantly earlier than the remaining species. The second clade includes Petaurista petaurista and Petaurista philippensis, the most widely distributed species in the genus, together with Petaurista xanthotis, endemic to south-central China and Tibet. Finally, the third clade comprises species mostly occurring in southern and central China (Petaurista albiventer, Petaurista alborufus, Petaurista yunanensis) and Indonesia (Petaurista elegans), with Petaurista hainana, which is endemic to Hainan Island (southeastern China). These two Petaurista subclades diverged during the late Miocene and estimates for the divergence of most species range from the late Pliocene to the early Pleistocene, that is coinciding with global cooling after 3 Ma and the start of northern hemisphere glaciation (Zachos et al., 2001). Other flying squirrel genera, such as Hylopetes, Petinomys and Glaucomys show a similar pattern, with estimates for the divergence of many extant species ranging from the late Pliocene to the early Pleistocene.
In the present, flying squirrels are widely distributed across the northern hemisphere, but only one genus lives in North America (Glaucomys) and a single species (Pteromys volans) occurs in Europe. Very few species are endemic to central and northern Asia. In contrast, flying squirrels are diverse in the tropical and subtropical forests of the Indo-Malayan region, which apparently acted both as refugia and diversification center (also for tree squirrels) since the late Miocene (Lu et al., 2013; Mercer and Roth, 2003).
Conclusions
Miopetaurista neogrivensis is the oldest unquestionable flying squirrel and dates back to the middle/late Miocene boundary (11.6 Ma). Its diagnostic wrist anatomy indicates that the two subtribes of flying squirrels had already diverged at that time. Moreover, this new fossil allows for a recalibration of flying squirrel time of origin and diversification, generally providing somewhat older estimates than previous molecular analyses. These differ according to the phylogenetic method used, total evidence analysis estimates an interval of 36.6 – 24.9 Ma while node dating results in a younger estimate of 30.6 – 17.4 Ma. Therefore, we cannot rule out that at least some of the oldest (ca. 36 Ma) fossils tentatively identified as flying squirrels may indeed belong to this group. However, the estimates of both independent phylogenetic approaches overlap for the late Oligocene (31 – 25 Ma), which should be considered the most likely interval for flying squirrel divergence. The two flying squirrel subtribes are found to have diverged during the early Miocene (22 – 18 Ma) while most extant genera would do so during the Miocene, although they are not recorded until the Pleistocene. Miopetaurista neogrivensis is estimated to have diverged from Petaurista spp., its sister taxon, between 18.8 – 12.4 Ma, the oldest boundary overlapping with the earliest record of the genus Miopetaurista (18 – 17 Ma). Perhaps not surprisingly, the skeletons of both genera show little differences. Sciurids are often regarded as a morphologically conservative group and flying squirrels are no exception having experienced few morphological changes for almost 12 million years.
Important information on data availability
3D surface models of the described material of M. neogrivensis (reconstructed skeleton, reconstructed skull, carpal bones) are available at MorphoBank https://morphobank.org/index.php/Projects/ProjectOverview/project_id/3108. All other data files are provided as supplementary data.
Materials and methods
Provenance and chronology
Request a detailed protocolThe partial skeleton of Miopetaurista neogrivensis (IPS56468) and all other described material (Figure 2—figure supplement 1-2) is housed in the Institut Català de Paleontologia Miquel Crusafont (ICP) in Sabadell (Barcelona, Catalonia, Spain). IPS56468 was found in 2008 during paleontological surveillance of excavation works at the Can Mata landfill (Abocador de Can Mata [ACM], els Hostalets de Pierola, Catalonia, Spain). It was unearthed from a rich fossiliferous horizon ACM/C5-D1 (ACM Cell 5, sector D, locality 1) in two different blocks with some elements in anatomical connection (Figure 2; Figure 7—figure supplement 1). The excavation of ACM/C5-D1 provided additional material of M. neogrivensis (Table 2). A partial, dorsoventrally crushed skull (IPS88677) was also recovered from stratigraphically close locality ACM/C8-Af (ACM Cell 8, sector A, locality f) (Alba et al., 2017).
The ACM series is located in the Vallès-Penedès Basin, an elongated half-graben filled mostly by continental deposits during the Miocene, which are rich in vertebrate fossils (Figure 9). The ACM composite series ranges from ca. 12.6 to 11.5 Ma, the age of the paleontological localities being well constrained thanks to high-resolution litho-, bio- and magnetostratigraphical data (Alba et al., 2017). ACM/C5-D1 is correlated to chron C5r.2n (11.657 – 11.592 Ma) (Alba et al., 2017), and its interpolated age is 11.64 Ma (Alba et al., 2017).
Anatomical terminology
Request a detailed protocolDental terminology, abreviations and measurement methods for sciurid cheek teeth follow Casanovas-Vilar et al., 2015 and references therein.
Three-dimensional data acquisition and reconstruction
Request a detailed protocolA virtual three-dimensional (3D) cranial reconstruction of M. neogrivensis was performed based on both IPS56468h and IPS88677 (Figure 6—figure supplement 1-2, Video 3). The specimens were analyzed separately by microfocus X-ray computed tomography (μCT) at the Multidisciplinary Laboratory of the ‘Abdus Salam’ International Centre of Theoretical Physics (Trieste, Italy), using a system specifically designed for the study of archaeological and paleontological materials (Table 5). Raw data from each scanning were imported (as stack of TIFF 8-bit files) to Avizo 7.0 and Rhinoceros 5.0 for segmentation, repositioning, mirroring and visualization. Each cranial bone or bone fragment was segmented virtually removing the surrounding matrix using semiautomatic thresholding tools and obtaining individual 3D digital models. Up to 64 3D models were generated for both specimens prior to repositioning and mirroring them to assemble an almost complete skull. Reconstruction primarily relied upon IPS56468h, but used 3D bone models from IPS88677 for elements that were particularly damaged or missing from the former (Table 3). Due to poor preservation of some fragments, it was necessary to import some of the 3D models to Rhinoceros 5.0 to repair the meshes, split the model to keep only the well-preserved regions, and use, if available, the mirrored region from the other side of the same skull, or alternatively take it from IPS88677 (mirrored if necessary). The 3D models were repositioned using Avizo 7.0 based on bilateral symmetry and fracture congruence. An almost complete 3D virtual skull model was finally assembled using a total of 41 3D models of cranial bones (Figure 6, Table 3). Besides cranial restoration, the two mandibles belonging to IPS56468 were scanned using a 3D desktop laser scanner (NextEngine) at high definition and with a dimensional accuracy of 0.13 mm (Macro Mode). Mandibles were exported to Rhinoceros 5.0 to repair meshes and were repositioned and joined to the resulting skull model using Avizo 7.0.
In addition, the posterior region of the skull of the extant flying squirrel Petaurista petaurista was scanned to compare its inner ear morphology to that of M. neogrivensis (Video 4). The specimen is kept at the Naturalis Biodiversity Center (NBC, Leiden, the Netherlands) and was scanned at this institution with a µCT scanner, using a Skyscan system model 1172 (Bruker company, Belgium) (Table 5).
The 3D surface of the postcranial skeleton of M. neogrivensis was digitized by means of photogrammetry using a variable number of sets of two-dimensional images. Within each set, individual images were captured with a minimum overlap of 66% using a digital single-lens reflex camera with a 50 mm f/2 lens. Aperture and focus were constant to ensure consistency between sets of photographs. Each set was then aligned and scaled to create a 3D model using the photogrammetric software Agisoft PhotoScan Professional. Individual surface models were exported to the 3D software Pixologic ZBrush version 4R7, which was used to merge the different meshes and improve the quality of the final merged 3D surface. Missing areas were modeled according to mirrored bones when available. Missing bone elements and broken areas in the available material were reconstructed using P. petaurista as reference. The morphology of the axial skeleton follows that of P. petaurista but has been modified according to the proportions of the recovered vertebrae.
Finally, carpal bones of extant Sciurus vulgaris, P. petaurista and Hylopetes sagitta, as well as those from M. neogrivensis were digitized for visualization and comparison. All extant specimens are stored in the collections of the NBC and were scanned there using a µCT scanner Skyscan system model 1172 (Table 5). Raw CT data from each scanning were imported directly (as stack of TIFF 16-bit files) to Avizo 7.0.
Body mass estimation
Request a detailed protocolBody mass (BM, in g) was estimated using an allometric regression of BM vs. skull length (SL, in mm), which is the estimator most tightly correlated with BM in extant rodents (Bertrand et al., 2016). SL was measured in M. neogrivensis from the 3D reconstruction based on crania IPS56468h and IPS88677 simultaneously (69.80 mm). The allometric regression was performed using log-transformed (ln) mean species data, as it is customary when predicting BM (Ruff, 2003). Sexes were treated separately (when possible) to avoid potentially confounding effects of body size dimorphism. The regression was computed for sciurids only to avoid biases due to allometric grade shifts among rodent families. SL sex-species means were calculated from published individual values (Bertrand et al., 2016) or measured by one of the authors (I.C.V.). BM average data were computed from individual published values (Bertrand et al., 2016) or taken from the literature (Hayssen, 2008; Thorington and Heaney, 1981; Zahler, 2001). The ordinary least-squares regression method, selected as the most suitable for prediction (Ruff, 2003; Smith, 1994), yielded the following equation with SPSS v. 17: ln BM = 4.369 ln SL – 11.384 (N = 33, p<0.001, r = 0.972, SEE = 0.256). The logarithmic detransformation bias was corrected using the quasimaximum likelihood estimator (Ruff, 2003; Smith, 1993): QMLE = 1.033. 50% confidence intervals (CI) for the prediction were computed.
Phylogenetic analyses and divergence dating
Request a detailed protocolTo infer the phylogenetic placement of M. neogrivensis and assess its impact on the dating estimates of the phylogeny of flying squirrels, we performed a new phylogenetic analysis of tribe Pteromyini by means of two alternative approaches: total evidence (Ronquist et al., 2012) and node dating (Ronquist et al., 2016). These approaches relied on the combination of morphological and molecular datasets, or on molecular datasets alone, respectively.
For the molecular dataset, we retrieved from GenBank the combination of the four genes more frequently used in previous phylogenetic studies of the Sciuridae (12S, 16S, cytochrome b, and irbp) (Mercer and Roth, 2003), which were downloaded for all Sciurinae species available in GenBank (Supplementary file 2). This dataset comprises 58 of the 89 extant Sciurinae species (Koprowski et al., 2016) and all genera but Biswamoyopterus. The Pteromyini are represented by 29 species, comprising 65% of extant diversity. We also obtained the sequence for Aplodontia rufa, sole extant member of the Aplodontiidae (the sister group of the Sciuridae; Fabre et al., 2012; Huchon et al., 2002), which was used as outgroup in the total evidence analysis. Each gene was aligned using two procedures: ribosomal coding genes were aligned by means of MAFFT v6 (Katoh et al., 2002) and protein coding genes were aligned using the translation alignment algorithm implemented in TranslatorX (using MAFFT to align proteins) (Abascal et al., 2010). Poorly-aligned regions in the ribosomal coding genes were eliminated with Gblocks (Castresana, 2000) under low stringency options (Talavera and Castresana, 2007). The morphological dataset consisted of 105 characters comprising the dentition, skull and postcranial skeleton of 36 extant taxa plus two fossil squirrels: M. neogrivensis and Douglassciurus jeffersoni (Appendix 2; Supplementary file 3). Douglassciurus is chosen because it is the earliest squirrel represented by abundant postcranial material and is currently recognized as an outgroup to all other sciurids (Appendix 1.1). The character list is mostly based in Thorington et al., 2002 with numerous additions of diagnostic pteromyin characters from the limb bones after Thorington et al., 2005. Finally, a set of new characters was also added to resolve the relationships of A. rufa and D. jeffersoni with other taxa (Appendix 2, Supplementary file 3).
For the total evidence dating we used MrBayes v. 3.2.6 (Ronquist et al., 2012). The analysis included 36 species for which we had both molecular and morphological data plus two fossil species with only morphological data. The molecular data consisted in a total of 3345 base pairs (bps) distributed in each gene as follows: 12S (505 bps), 16S (521 bps), cytb (1140 pbs), irbp (1179 pbs). The protein coding genes were split in codon positions and the most appropriate partitioning scheme and the model of molecular evolution for each partition were estimated by means of the software Partitionfinder (Lanfear et al., 2012) (Table 6). The five partitions derived from Partitionfinder were then concatenated with the previously described morphological dataset. We modeled node ages and tree topology using a fossilized birth and death process (Gavryushkina et al., 2017; Heath et al., 2014; Stadler, 2010; Zhang et al., 2016) with broad priors on speciation (exp[10]), extinction (beta[1,1]), and fossilization (beta[1,1]) (Pyron, 2017; Zhang et al., 2016).The analysis also relied on a relaxed-clock with the independent gamma-rates (IGR) model (Ronquist et al., 2012), with a broad prior for the variance increase parameter (exp[10]; Pyron, 2017). Following Pyron (Pyron, 2017) we used 1/mean of the root age and exp(1/mean of the root age) to generate a wide clock rate prior compatible with both morphological and molecular clocks. For the root age we used the interval 52.7–45.7 Ma as a uniform prior. The lower bound of this prior was informed by the most recent age estimate (MRE) of the oldest Aplodontiidae (Table 4; see also Appendix 1.1) and we calculated the upper bound using the algorithm proposed by Hedman, 2010, which calculates the probability distribution of the age of a clade, given the ages of the oldest fossil representatives of the outgroups of that clade. We used the upper limit of the 95% confidence interval for the age of the root as a plausible upper bound for the root. This was informed by the MRE for the minimum age of the oldest Gliridae (47.4 Ma, outgroup to Aplodontiidae+ Sciuridae (Blanga-Kanfi et al., 2009; Fabre et al., 2012; Huchon et al., 2002) and the MRE for the minimum age of Erlianomys combinatus, the oldest Myodonta (53.9 Ma, outgroup to Gliridae +Aplodontiidae + Sciuridae; Blanga-Kanfi et al., 2009; Fabre et al., 2012; Huchon et al., 2002); Table 4, Appendix 1.1). The maximum age oldest rodent was used as the maximum possible age of the root (56 Ma; Table 4, Appendix 1.1).
Partitions and nucleotide substitution models for the molecular data were estimated by means of Partitionfinder v. 2.1.1 (Lanfear et al., 2012) with linked branch lengths, a Bayesian Information Criterion (BIC) model of selection and a greedy search algorithm (Table 6). For the morphological data (only variable characters included), we used a k-state Markov (Mkv) model (Lewis, 2001) with a rate variation modeled by means of a discrete gamma model. We distinguished between ordered and unordered characters.
We also used an alternative Bayesian program, BEAST v. 1.8.4 (Drummond et al., 2012), to provide an independent estimate of the time of the Pteromyini/Sciurini divergence and the onset of Pteromyini crown diversification by means of node dating. This analysis included all species of Sciurinae available in GenBank (58 taxa) and consisted in 3225 pbs distributed in each gene as follows: 12S (383 bps), 16S (523 bps), cytb (1140 pbs), irbp (1179 pbs). Protein coding genes were subsequently split in codon positions and we used Partitionfinder to determine the best set of molecular partitions and models of molecular evolution (Table 6). We estimated the tree in time units using an uncorrelated lognormal clock applied to each of the five partitions derived from Partitionfinder and two different calibration points (one located in each tribe; Table 4). Miopetaurista neogrivensis from ACM/C5-D1 (11.6 Ma) provides a minimum age for the divergence between Petaurista and the remaining Pteromyina, while the oldest record of the genus Sciurus (dating back to 13.6–10.3 Ma) gives minimum age for the Sciurus/Tamiasciurus split (Table 4, Appendix 1.2). Soft maxima for both calibrations were calculated as described in Hedman, 2010. The soft maximum for the calibration point in the Pteromyina was informed by the MRE of the oldest Heteroxerus, (25.0 Ma), the earliest Xerinae (sister group to Sciurinae), and by the oldest estimate of D. jeffersoni (35.8 Ma), basal to crown Sciuridae (Table 4, Appendix 1.2). For the calibration point in the Sciurini, the soft maximum was informed by the age of M. neogrivensis and the ages of oldest Heteroxerus and D. jeffersoni stated above. Following the ‘consistent approach’ (Hedman, 2010) the minimum age of the oldest Aplodontiidae (45.7 Ma) is taken as the maximum possible age of both calibration points.
Both MrBayes and BEAST analyses relied on four independent runs of 25,000,000 generations, each sampled every 4000 generations. Each of the runs in MrBayes relied on four Markov chain Monte Carlo (MCMC) chains.
Appendix 1
Comments on the fossils used to calibrate the phylogenetic trees
1.1 Total evidence phylogenetic analysis
The prior assigned to root age relied on the following fossils: the oldest Rodentia, the oldest Aplodontiidae, the oldest Gliridae and the oldest Myodonta (see also Table 4). The oldest rodent known to date is Acritoparamys atavus (Reithroparamyinae, Ischyromidae) from the latest Paleocene (early Clarkforkian, 56.0 – 55.8 Ma) of Montana (Janis et al., 2008; Korth, 1994). Spurimus selbyi is the oldest aplodontid (Flynn and Jacobs, 2008; Korth, 1994), and Eogliravus wildi, the earliest glirid (Daams and Bruijn, 1995; Dawson, 2003; Hartenberger , 1971). Spurimus is a basal member of the aplodontiids, the sister group of sciurids (Fabre et al., 2012; Huchon et al., 2002), dating back to the middle Eocene (early Uintan, 46.3 – 45.7 Ma) of Wyoming (Flynn and Jacobs, 2008). Eogliravus wildi is a stem glirid known from isolated cheek teeth from France (Hartenberger J-L and Hartenberger, 1971) and an exceptionally preserved specimen (articulated skeleton with pelage impression) from the Messel Lagerstätte in Germany (Storch and Seiffert, 2007). Eogliravus and Gliravus are distinguished from all other glirids by their plesiomorphic protrogomorphous zygomasseteric anatomy which they share with basal rodents (Daams and Bruijn, 1995; Storch and Seiffert, 2007; Vianey-Liaud, 1985). The French fossils date back to the earliest Eocene (biozone MP10) (Daams and Bruijn, 1995; Dawson, 2003), yielding an age range of 50.7 – 47.4 Ma (age boundaries after Vanderberghe et al., 2012). Finally, Erlianomys combinatus, recovered at the lower part of the Arshanto Fromation of Inner Mongolia (China) (Li and Meng, 2010; Wu et al., 2012), represents the earliest record of myodonts, the mouse-like rodent clade (Blanga-Kanfi et al., 2009; Fabre et al., 2012; Huchon et al., 2002) (sister group to Aplodontiidae + Sciuridae + Gliridae; Blanga-Kanfi et al., 2009; Fabre et al., 2012; Huchon et al., 2002). This occurrence has been relatively dated by means of bio- and magnetostratigraphy, yielding an age range of 54.8 – 53.9 Ma (early Eocene) (Wang et al., 2010). All these data were considered in the calculation of root age and a uniform prior of 52.7 – 45.7 was finally used.
The oldest record of Douglassciurus jeffersoni, one of the oldest sciurids, is included as tip date. Douglassciurus jeffersoni dates back to the late Eocene (middle Chadronian, 36.6 – 35.8 Ma) of Wyoming and is the oldest-known tree squirrel (Emry and Korth, 1996; Goodwin, 2008), being somewhat younger than the oldest purported flying squirrel, Hesperopetes thoringtoni (Emry and Korth, 2007). Although it is not a Sciurinae, the postcranial skeleton of Douglassciurus is surprisingly similar to that of extant tree squirrels such as Sciurus (Emry and Thorington, 1982; Emry and Thorington, 1984), but the archaic protrogomorphous zygomasseteric structure of the skull (Emry and Korth, 1996; Emry and Thorington, 1982) is shared with aplodontiids and other primitive rodent groups such as ischyromids. Consequently, Douglassciurus is generally considered the most basal sciurid, since all other extinct genera, including other stem taxa are sciuromorphous or show a modified protrogomorphous zygomasseteric structure (Black, 1963; Goodwin, 2008; Korth and Emry, 1991). This basal position is confirmed by our results (see Figure 7). Finally, the record of Miopetaurista neogrivensis from ACM/C5-D1 (11.63 Ma) is also considered as tip date in this analysis.
1.2 Node dating phylogenetic analysis
Node dating for the Sciurinae (Figure 7—figure supplement 1) considered two different calibration points: one for the divergence between Sciurus and Tamiasciurus, and another between Petaurista and the remaining Pteromyina, that is, the large-sized flying squirrel clade (Table 4). The minimum age for the split between Sciurus and Tamiasciurus is given by the oldest record of the former genus, Sciurus olsoni (Emry et al., 2005), dating back to the middle–late Miocene (Clarendonian, 13.6 – 10.3 Ma) from the Truckee Formation of Nevada (USA) (Emry et al., 2005; Goodwin, 2008; Woodburne, 2004). Protosciurus and Miosciurus, which are recorded from the Oligocene of North America, have been considered sciurins (Black, 1963; Goodwin, 2008), but their archaic zygomasseteric structure indicates that they should be better regarded as stem sciurids. Miopetaurista neogrivensis from ACM/C5-D1 (11.6 Ma) provides a minimum age for the divergence of Petaurista from the remaining Pteromyina, because Miopetaurista is recognized as its sister taxon by total evidence analysis (Figure 7).
The soft maximum for the calibration point in the Pteromyina was informed by the most recent age estimate of the earliest Xerinae, the sister group to the Sciurinae (Mercer and Roth, 2003; Steppan et al., 2004). Xerines are the most diverse subfamily of extant squirrels (Koprowski et al., 2016; Thorington et al., 2012), including the speciose Holarctic ground squirrels (Marmotini) as well as the African ground (Xerini) and tree squirrels (Protoxerini). Ground squirrels (sensu lato) have a rich fossil record and include some of the oldest-known squirrels, such as Palaeosciurus goti from the early Oligocene of France (Vianey-Liaud, 1974; Vianey-Liaud, 1975), the earliest squirrel to be recorded in Eurasia. A high diversity of ground squirrels is known from the late Oligocene to the middle Miocene of North America, including Nototamias, Palaearctomys and Protospermophilus, among other genera (Goodwin, 2008; Korth, 1994). While some authors have included these genera in the marmotin tribe (Black, 1963; Korth, 1994), recent revisions regard them as stem ground squirrels of unknown affinities (Goodwin, 2008), because they lack the diagnostic features characterizing extant subfamilies. The same can be argued about the genus Palaeosciurus, which is recognized as a ‘ground squirrel’ on the basis of limb proportions (Vianey-Liaud, 1974; Vianey-Liaud, 1975), but does not show the diagnostic features of extant groups. The genus Miospermophilus is tentatively included into the marmotini by some authors (Goodwin, 2008; Korth, 1994), which (if correct) would push the record of the subfamily as far back as 30.0–28.0 Ma (early Arikarean, late Oligocene) (Goodwin, 2008). However, the skull of this genus, which would unambiguously show if it is a marmotin or not, is unknown. Xerins are the most basal tribe of xerines and include just six extant species. Their oldest representative is Heteroxerus, which is first recorded in the late Oligocene (Chattian, MP28) of Spain (Álvarez-Sierra et al., 1999), with an estimated age of 26.4–25.0 Ma (Barberà et al., 2001). This genus is abundantly represented in the European Miocene by mandibular and cranial fragments as well as isolated cheek teeth (De Bruijn, 1999). The soft maximum for this calibration point is also informed by the minimum age of the oldest estimate of D. jeffersoni, basal to crown Sciuridae (Table 4). For the calibration point in the sciurins, the soft maximum was informed by the age of M. neogrivensis and the ages of oldest Heteroxerus and D. jeffersoni stated above. The minimum age of the oldest aplodontiid (45.7 Ma) is taken as the maximum possible age of both calibration points.
Appendix 2
Morphological characters for phylogenetic analyses
Phylogenetic analyses consider 36 different species including representatives of all extant flying squirrel genera but Biswamayopterus, which is known by very few specimens (Koprowski et al., 2016; Thorington et al., 2012). In addition, species of the tree squirrels Ratufa, Sciurus, Callosciurus and Protoxerus are included. Finally the outgroups include D. jeffersoni, the oldest-known squirrel for which well-preserved cranial and postcranial material is available (Emry and Korth, 1996; Emry and Thorington, 1982), and Aplodontia rufa. The latter is the only living aplodontiid, which has been long recognized as the sister group of sciurids (Simpson, 1945), as confirmed by recent molecular phylogenetic analyses (Adkins et al., 2001; Fabre et al., 2012; Huchon et al., 2002; Montgelard et al., 2008; Montgelard et al., 2002). The morphological character matrix comprises 105 characters: 45 dental (1–45), 26 cranial (46–71), 31 postcranial (75–105) and three characters related to the muscles that support the patagium (72–74). Dental and cranial characters are based on Thorington et al., 2002 with ten new characters added in order to resolve the relationships of A. rufa and D. jeffersoni to other taxa. Most of these new characters are taken from Korth and Emry, 1991. The inclusion of A. rufa as outgroup implied some modifications and additions in character states of the matrix of Thorington et al., 2002. The postcranial characters are mainly based in Thorington et al., 2005. For D. jeffersoni the character list from Thorington et al., 2002 was completed with data from Emry and Korth, 1996. Dental and cranial character states for most extant taxa were directly recorded after the study of the material housed within the collections of the Naturalis Biodiversity Center of Leiden (Aeromys tephromelas; Aplodontia rufa; Belomys pearsonii; Callosciurus notatus; Eupetaurus cinereus; Glaucomys sabrinus; Glaucomys volans; Hylopetes alboniger; Hylopetes phayrei; Hylopetes sagitta; Iomys horsfieldii; Petaurista alborufus; Petaurista elegans; Petaurista petaurista; Petaurista philippensis; Petinomys setosus; Protoxerus stangeri; Pteromys volans; Ratufa bicolor; Ratufa affinis; Sciurus carolinensis; Sciurus vulgaris), the Muséum National d’Histoire Naturelle of Paris (Iomys horsfieldii, Petaurillus kinlochii, Pteromyscus pulverulentus, Trogopterus xantiphes), and the Institut Català de Paleontologia Miquel Crusafont (Callosciurus notatus; Glaucomys sabrinus). Similarly, postcranial character states were directly recorded for Aplodontia rufa, Callosciurus notatus, Glaucomys sabrinus, Glaucomys volans, Hylopetes sagitta, Petaurista petaurista, Ratufa affinis, Sciurus carolinensis and Sciurus vulgaris. Characters on the muscles that support the patagium and characters for all other extant species in the analyses were taken from the literature. Postcranial material is generally rare in museum collections of extant rodents, and generally only pelts and skulls are available for most species. For certain genera such as Petaurista or Hylopetes only the postcranial material of a single species could be studied. Therefore, for most flying squirrel species character states are based in Thorington et al., 2005. We have assumed that postcranial character states are the same for all the species within a given genus, which is apparently reasonable since the diagnoses for most squirrel genera and species take into account cranial and dental traits (Ellerman, 1940; Moore, 1959). Below we list the morphological characters used in the analyses. Character treatment is indicated with the abbreviations ‘O’ for ordered characters and ‘U’ for unordered characters. See Supplementary file 3 for the complete morphological character matrix.
Hypselodont cheek teeth (O): (0) absent; (1) present.
P3 (O): (0) absent; (1) very small to small; (2) intermediate to large.
P4 size relative to M1 size (O): (0) smaller; (1) approximately equal; (2) larger.
Crenulation in the upper cheek teeth (O): (0) none; (1) low; (2) heavy.
Enamel surface in the upper cheek teeth (U): (0) smooth; (1) folds: lumpy appearance; (2) pitting; (3) well-defined lophules.
Color of the upper incisors (U): (0) deep orange; (1) yellowish orange; (2) yellow; (3) variable in color.
Longitudinal lophs in the upper cheek teeth (O): (0) none; (1) one; (2) two or more.
Hypocone (U): (0) absent; (1) protocone and hypocone separated from each other by a small groove; (2) large protocone connected by a thin loph tow a smaller hypocone, but both cusps moderately separated; (3) large hypocone well separated from the protocone.
Mesostyle on the P4 (U): (0) absent; (1) variable; (2) present; (3) present and very broad.
Mesostyle on the M1 and M2 (U): (0) absent; (1) variable; (2) present; (3) present and very broad.
Mesostyle on the M3 (U): (0) absent; (1) variable; (2) present; (3) present and very broad.
Mesostyle position in P4 (U): (0) N/A; (1) connecting with the paracone; (2) connecting with the metacone; (3) isolated.
Mesostyle position on the M1/M2 (U): (0) N/A; (1) connecting with the paracone; (2) connecting with the metacone; (3) isolated.
Mesostyle position on the M3 (U): (0) N/A; (1) connecting with the paracone; (2) connecting with the metacone; (3) isolated.
Mesoloph on the P4 (O): (0) obsolete; (1) short; (2) long.
Mesoloph on the M1/M2 (O): (0) obsolete; (1) short; (2) long.
Mesoloph on the M3 (O): (0) obsolete; (1) short; (2) long.
Divided parastyle on the P4 (O): (0) absent; (1) present.
Divided parastyle on the M1/M2 (O): (0) absent; (1) present.
Protoconule on the P4 (O): (0) not evident; (1) incorporated into the protoloph but evident; (3) pinched off from the protoloph.
Protoconule on the M1/M2 (O): (0) not evident; (1) incorporated into the protoloph but evident; (3) pinched off from the protoloph.
Protoconule on the M3 (O): (0) not evident; (1) incorporated into the protoloph but evident; (3) pinched off from the protoloph.
Metaconule on the P4 (O): (0) not evident; (1) incorporated into the metaloph but evident; (3) pinched off from the metaloph.
Metaconule on the M1/M2 (O): (0) not evident; (1) incorporated into the metaloph but evident; (3) pinched off from the metaloph.
Metaconule on the M3 (U): (0) not evident; (1) incorporated into the metaloph but evident; (3) pinched off from the metaloph.
Metaloph connection on the P4 (U): (0) absent; (1) connected to the protocone; (2) connected to the hypocone.
Metaloph connection on the M1/M2 (U): (0) absent; (1) connected to the protocone; (2) connected to the hypocone.
Metaloph connection on the M3 (U): (0) absent; (1) connected to the protocone; (2) connected to the hypocone.
Protoloph-protocone connection on the P4 (U): (0) absent; (1) weak; (2) variable; (3) strong.
Protoloph-protocone connection on the M1/M2 (U): (0) absent; (1) weak; (2) variable; (3) strong.
Protoloph-protocone connection on the M3 (U): (0) absent; (1) weak; (2) variable; (3) strong.
Metaloph connection to protocone/hypocone on P4 (U): (0) N/A; (1) weak; (2) variable; (3) strong.
Metaloph connection to protocone/hypocone on M1/M2 (U): (0) N/A; (1) weak; (2) variable; (3) strong.
Metaloph connection to protocone/hypocone on M3 (U): (0) N/A; (1) weak; (2) variable; (3) strong.
Distolingual flexus (O): (0) absent; (1) present.
Lingual cingulum (O): (0) absent; (1) present.
Metacone on M3 (O): (0) absent; (1) present.
Hypocone on M3 (O): (0) absent; (1) present.
Anterolophid and metalophid development in the lower molars (U): (0) anterolophid and metalophid absent; (1) anterolophid and metalophid complete and joining the metaconid; (2) metalophid joining anterolophid before reaching the metaconid; (3) anterolophid reaching the metaconid and protolophid terminating in the talonid basin.
Mesostylid development in all lower cheek teeth (O): (0) absent; (1) present; (2) present and very broad.
Anterosinusid (O): (0) obsolete; (1) present.
Hypoconulid on the p4 (O): (0) obsolete; (1) present.
Hypoconulid on the m1/m2 (O): (0) obsolete; (1) present.
Hypoconulid on the m3 (O): (0) obsolete; (1) present.
Complexity of the lower cheek teeth pattern (O): (0) simple; (1) somewhat complex; (2) complex.
Zygomasseteric structure (O): (0) protrogomorphous; (1) sciuromorphous.
Orientation of the zygomatic plate (O): (0) completely horizontal; (1) rising only slightly above the level of the cheek teeth; (2) anterodrosally tilted.
Shape of the infraorbital foramen (O): (0) round; (1) laterally compressed.
Position of the infraorbital foramen (O): (0) high; (2) low.
Orientation of angular process in the mandible (O): (0) parallel to the rest of mandibular processes; (1) bent medially.
End of masseteric fossa in the mandible (O): (0) below m1; (1) below p4.
Pterygoid ridges (O): (0) short; (1) elongated.
Basioccipital orientation (O): (0) vertical; (1) tilted ventrally.
Length of the rostrum (O): (0) short; (1) intermediate; (2) long.
Width of the rostrum (O): (0) wide; (1) intermediate; (2) narrow.
Length of incisive foramina (O): (0) short; (1) intermediate; (2) long.
Masseteric tubercles (O): (0) none; (1) not prominent; (2) somewhat prominent; (3) prominent.
Shape of masseteric tubercles (O): (0) N/A; (1) flat; (2) bulbous; (3) bulbous and tine-like.
Openness of the distal half of the orbit (O): (0) mostly open; (1) somewhat open; (2) slightly open.
Postorbital process of the frontal (O): (0) absent; (1) short; (2) long.
Dorsal process of the jugal (O): (0) absent; (1) short; (2) long.
Orientation of the pterygoids (O): (0) parallel; (1) pinched in.
Zygomaticomandibularis fossa (O): (0) deep, (1) shallow; (2) obsolete
Stapedial artery (O): (0) enclosed within a bony tube; (1) no stapedial artery (i.e., lack of stapedial and sphenofrontal foramina).
Shape of the auditory bulla (O): (0) flask-shapped; (1) oval.
Number of septa per bulla (U): (0) zero to one; (1) two to three; (2) more than three, not arranged in a honeycomb pattern; (3) more than three and arranged in a honeycomb pattern; (4) more than three, arranged as a sliver.
Width of the bulla (O): (0) narrow; (1) intermediate; (2) wide.
Inflation of the bulla (O): (0) inflated; (1) flat.
Distance between bullae (O): (0) narrow; (1) intermediate; (2) wide.
Size of the bony process from the opening of ear canal (O): (0) small; (1) large.
Mastoid process inflated (O): (0) no; (1) yes.
Origin of semitendinosus III muscle (O): (0) inseparable from origin of semitendinosus II; (1) caudal vertebrae 5–8; (2) caudal vertebrae 9–15.
Semitendinosus III discrete or joined (O): (0) inseparable from semitendinosus I and II; (1) muscular attachment to semitendinosus I, II; (2) tendinous attachment to semitendinosus I, II; (3) separate, but lies just caudal to the other heads of semitendinosus; (4) well separated from the rest of musculature, deeply embedded in the skin of the uropatagium.
Insertion of the tibiocarpalis muscle (U): (0) absent; (1) distal tibial tuberosity; (2) metatarsal I; (3) metatarsal II.
Elevated scapholunate tuberosity on the pisiform (O): (0) absent; (1) present.
Triquetral process on the pisiform (O): (0) absent; (1) present.
Twisting of the pisiform (O): (0) yes; (1) no.
Enlarged surface for pisiform-triquetrum articulation on the ulnar flange (O): obsolete; (1) present, not prominent; (2) present, prominent.
Enlarged surface for pisiform-triquetrum articulation on the radial flange (O): obsolete; (1) present, not prominent; (2) present, prominent.
Articulation between metacarpal II and central (O): (0) yes; (1) no.
Scapholunate-capitate articulation (O): (0) no; (1) yes.
Centrale-hamate articulation (O): (0) no; (1) yes.
Lesser multangular-scapholunate articulation (O): (0) no; (1) yes.
Size of the triquetral facet on the scapholunate (O): (0) short; (1) long.
Shape of the pisiform facet on scapholunate (O): (0) circular; (1) grooved.
Proximal falciform process on scapholunate (U): (0) obsolete; (1) present, prominent; (2) present, not prominent.
Distal falciform process on scapholunate (U): (0) obsolete; (1) present, prominent; (2) present, not prominent.
Angle between proximal falciform process and body of the scapholunate (O): (0) perpendicular; (1) acutely angular.
Distal tibial tuberosity (O): (0) absent; (1) present.
Patagium (O): (0) absent; (1) present.
Fused scapholunate (O): (0) no; (1) yes.
Deltoid ridge on humerus (U): (0) very broad and prominent; (1) broad to moderately broad; (2) very narrow and knife-like.
Lateral supracondylar ridge on humerus (U): (0) large, extending to the lateral edge of the epicondyle; (1) reduced to greatly reduced, not extending to the lateral edge of the epicondyle.
Medial epicondyle of the humerus (U): (0) elongated; (1) reduced.
Radio-ulnar joint (O): (0) moveable (connected by ligaments or diarthrosis); (1) fused by a prominent syndesmosis.
Process that supports the distal compartment for muscle abductor pollicis longus in the distal radius (U): (0) long; (1) short.
Olecranon process of ulna (U): (0) long; (1) short.
Distal end of ulna (U): (0) robust and broad; (1) delicate and narrow.
Styliform process of ulna (U): (0) similar in size to styliform process of radius; (1) more prominent than styliform process in radius.
Third trochanter of femur (U): (0) prominent (may be ridge-like); (1) moderately to highly reduced.
Position of the third trochanter of femur (U): (0) distal to the lesser trochanter; (1) opposite and proximal to the lesser trochanter.
Popliteal fossa of tibia (U): (0) deep; (1) intermediate; (2) shallow.
Anterior (cnemial) crest of the tibia (U): (0) prominent; (1) intermediate; (2) faint.
Pit on the plantar surface of the astragalus between the calcaneal articular surfaces (U): (0) present; (1) absent.
Position of the sustentacular process of the calcaneus (U): (0) opposite to the peroneal process; (1) slightly distal to the peroneal process; (2) clearly more distal to the peroneal process.
Appendix 3
Detailed morphological descriptions and comparisons
3.1 Cheek tooth morphology and species attribution
IPS56468. The cheek teeth of the associated skull (Figure 5—figure supplement 1h-i) are moderately worn, so certain details cannot be appreciated, especially in the lower series. The enamel in the upper cheek teeth is smooth and the main transverse ridges are thick and simple. The P3 is small and rounded, presenting a single anterolingual cusp. The P4 shows a well-developed parastyle forming a broad platform that makes this premolar longer that wide (see Supplementary file 4). The parastyle is as developed as the main labial cusps (paracone and metacone). The protoloph and the metaloph in the P4–M2 are higher than the remaining transverse ridges and parallel, reaching the protocone separately. The protoloph shows a thickening near its lingual end, presumably corresponding to the protoconule, in the P4 and M3, but this conule is missing in the M1 and M2. The metaconule is absent in all upper cheek teeth. The metaloph shows a short distal spur in the P4, whereas a similarly-developed mesial spur is present in the protoloph of the M3. The protocone is longitudinally elongated and continuous with the ectoloph, whereas the hypocone is highly reduced and cannot be distinguished from this ridge. The mesostyle is present in the P4 and M3 but not in the M1 and M2. The M3 is relatively short and shows only two main cusps, protocone and paracone, and a wide valley delimited by a continuous ridge. This molar presents a short, thin and arched metaloph which forms a reduced distolingual valley. This valley is closed by a small style connected to the metaloph by means of a short spur.
In contrast to the upper cheek teeth, the enamel in the lower ones shows extensive meandriform crenulations that mostly affect the valleys but also extend to the main ridges. There are three well-defined cusps (protoconid, metaconid and hypoconid) and three smaller ones (entoconid, mesoconid and mesostylid) surrounding a relatively shallow central valley. The mesostylid is placed quite mesially, next to the metaconid. The hypoconulid is absent. The ectolophid is M-shaped and the sinusid reduced and open. The anterosinusid is highly reduced and closed by an anterior cingulid descending from the anterolophid. The metalophid is short in the p4–m2, not reaching the metaconid. The p4 is markedly longer and narrower than the molars and presents a small and rounded anteroconulid. It shows a short, distally-directed mesolophid placed just distal to the protoconid. The entolophid is thick but it is interrupted before reaching the entoconid. In contrast, the entolophid is thick and complete in the m1–m3. In these molars the entolophid is transverse, merging the entoconid with the ectolophid between the mesoconid and the hypoconid. In addition, in the m1 the mesostylid shows a labially-directed ridge pointing towards an opposite ridge that departs from the ectolophid. These two ridges are interpreted as an incomplete mesolophid and cannot be discerned in the m2. The m3 shows a complex morphology, with a low metalophid that bifurcates into two lower meandriform ridges that join the mesostylid and the mesoconid, respectively. A vestigial mesolophid is also present in the m3 between the mesoconid and the protoconid.
Additional material. Additional remains from ACM/C5-D1 include a maxillary fragment with the broken P4–M2 (IPS43724), two partial mandibles with p4–m3 (IPS43677, IPS85340), 37 isolated cheek teeth and additional cheek teeth fragments (see Table 2 and Supplementary file 4). In addition, a partial cranium (IPS88677) has been recovered in the approximately stratigraphically equivalent site ACM/C8-Af (Alba et al., 2017) (Figure 6—figure supplement 2). A highly fragmentary third cranium (IPS85410) comes from locality ACM/C6-A5, with an estimated age of also 11.6 Ma (Alba et al., 2017). IPS85410 includes only part of the frontal, parietal and maxillary with damaged P4–M3. Finally, another partial mandible with p4–m3 (IPS87560) has been recovered in sector ACM/C8-B (Figure 3f–g). Age range for the localities in this sector of the ACM series is from 11.6 to 11.7 Ma (Alba et al., 2017), so the age of all additional specimens is very close to that of the partial skeleton IPS56468.
Overall the size and morphology of these additional specimens is very similar to that of IPS56468. However, this additional material shows that some occlusal features are variable, such as the development of the metalophid, the mesolophid and the entolophid in the lower cheek teeth. The metalophid is generally short and does not reach the metaconid, however, this ridge is complete in one isolated p4 (IPS77873) and two m1/m2 (IPS43675, IPS77879). The entolophid can be partial in some premolars, but this ridge is always thick and complete in the molars. A well-defined short to medium-length mesolophid can be recognized in several individuals, such as in IPS43677. However, this ridge (if present) is usually too low to be clearly distinguished from the abundant crenulations that occur in the central valley. Three isolated dp4 have been recovered. These are much smaller but morphologically similar to the p4, although the anteroconulid is more reduced and situated closer to the protoconid. As far as the upper cheek teeth are concerned, the P4 is more variable. A metaconule is generally present, but the protoconule is distinct only in a few specimens. The number of distally-directed spurs in the metalophule varies between two and three, with one specimen (IPS43675) showing a mesially-directed spur as well. The parastyle of the P4 is invariably defined by a single cusp, and a mesostyle is present in all instances. The M1 and M2 do not show conules in the protoloph and the metaloph, but the later ridge may develop a short distally-directed spur as in the P4. The recovered isolated M3 show a well-defined protoconule, a small mesostyle and a small distolingual flexus. Even though the enamel is smooth in the upper cheek teeth, a few specimens (IPS77859, IPS77862, IPS77863) present weak crenulations in the protoloph and metaloph.
Discussion. The diagnoses of fossil ‘flying’ squirrel genera and species rely entirely on cheek tooth morphology and dimensions. The partial skeleton IPS56468 and additional material described here (Table 2) are attributed to the genus Miopetaurista. In the ACM localities a second genus of ‘flying’ squirrel, Albanensia, is recorded, but Miopetaurista is clearly distinguished by its larger size, absence of enamel crenulations in the upper cheek teeth, and protoloph and metaloph of the P4–M2 parallel instead of lingually convergent (Daxner-Höck and Mein, 1975; De Bruijn, 1999; Mein, 1970). In addition, in Miopetaurista the protoloph and metaloph rarely present conules, the hypoconulid is not apparent in the lower cheek teeth, and the anterosinusid is always present, although reduced (Daxner-Höck and Mein, 1975; De Bruijn, 1999; Mein, 1970).
Miopetaurista ranges from the early Miocene to the Pliocene of Eurasia and it includes up to ten different species, although some of them are poorly known and the validity of others is questionable (Casanovas-Vilar et al., 2015; De Bruijn, 1999; De Bruijn, 1995). The specimens described here are outstanding because of their large size (Supplementary file 4), clearly above the size range of most species of the genus, particularly the older ones such as Miopetaurista diescalidus, Miopetaurista dehmi and Miopetaurista lappi. They further differ from these species in several morphological features including the absence of mesoloph in the upper cheek teeth, the different shape of the P4, and the presence of a complete entolophid in the lower cheek teeth. The roughly coeval Miopetaurista gaillardi is also smaller, but shows a similar morphology, including the presence of a complete entolophid in some lower molars. However, this species is distinguished by its shorter P4 as result of the less developed parastyle. Miopetaurista gibberosa, the type species of the genus, is only known from a single mandible from the middle Miocene (MN5) of Göriach (Austria) (Daxner-Höck and Höck, 2015). The lower cheek teeth of this species are smaller than the described material and in addition they show a highly reduced entolophid. The size of the described specimens is only comparable to those of Miopetaurista neogrivensis, Miopetaurista crusafonti and Miopetaurista thaleri. The latter species is slightly larger and further displays a number of morphological differences (Mein, 1970), such as the presence of a complete mesolophid that joins the entoconid in the lower molars. The described material is well above the size range of M. crusafonti, which is recorded in other Vallès-Penedès sites of similar age (Casanovas-Vilar et al., 2015; Casanovas-Vilar et al., 2010). However, the described molars always show a complete entolophid, whereas this ridge is always partial in M. crusafonti. Furthermore, in the studied material the mesolophids are better developed than in M. crusafonti and the P4 is less triangular (Casanovas-Vilar et al., 2015). The described specimens fit perfectly in size and morphology with M. neogrivensis, which is present elsewhere in the Vallès-Penedès Basin (Aldana Carrasco, 1991; Aldana Carrasco, 1992; Casanovas-Vilar et al., 2015). This species is only surpassed in size by M. thaleri and shows a more complex occlusal morphology than M. crusafonti.
Mein, 1970 described a number of dental features in extant and fossil ‘flying’ squirrels that were used to distinguish three different groups mainly on the basis of enamel wrinkling and the presence of additional lophules. According to this author Miopetaurista (= Cryptopterus in Mein (1970) and Daxner-Höck and Mein, 1975) would belong to the second group, characterized by its smooth enamel and absence of lophules. This group would also include the extant genera Glaucomys, Eoglaucomys and Iomys. Other common characters of the group are the absence of metaloph on M3, the absence of mesostyle, the presence of anterosinusid in the lower molars and the absence of anteroconid. Although the mesostyle is generally absent in the upper molars of Miopetaurista, it frequently occurs in the P4. On the other hand, the anteroconid is absent from the lower molars but present in the p4 and dp4. These dental characters have proven to be inadequate to distinguish between clades of flying squirrels, or even to separate flying squirrels from other squirrels (Thorington et al., 2002; Thorington et al., 2005). Petaurista is characterized by its smooth enamel but shows a much more complex morphology, with additional longitudinal lophules and would belong to a completely different group than Miopetaurista according to dental classification (Mein, 1970). This contradicts our phylogenetic results (see main text). Nevertheless, there is a remarkable dental similarity between Miopetaurista and Petaurista, namely the presence of a distolingual flexus. This occurs in all upper cheek teeth in Petaurista, whereas in Miopetaurista it only occurs in the M3. Such feature is only observed in two other genera from the subtribe Pteromyina, Aeretes and Eupetaurus, being absent in the Glaucomyina.
3.2 Skull morphology
The reconstruction of skull morphology is based on the almost complete cranium and mandible of IPS56468 from ACM/C5-D1 (Figure 3a–e, Figure 6—figure supplement 1), as well as on the additional information provided by the partial cranium IPS88677 from ACM/C8-Af (Figure 6—figure supplement 2) and the partial mandible IPS87560 from sector ACM/C8-B (Figure 3f–g). The cranium IPS56468 is laterally crushed in an oblique angle, with the rostrum deviated towards the left side. The zygomatic arch is complete but crushed, with the jugal and other fragments pressed against the orbit. On the other hand, IPS88677 is dorsoventraly crushed and most of the premaxillary and frontal bones are not preserved. The virtual reconstruction of the skull (see Materials and Methods; Figure 6 and Video 3) relied primarily on IPS56468h, but used some elements from IPS88677 when the bone was particularly damaged or missing in the former specimen (zygomatic plate, jugal, pterygoid, alisphenoid; Figure 6—figure supplement 1-2, Table 3). The description of cranial morphology relies on the original specimens but also considers aspects of the virtual reconstruction concerning overall shape and proportions of certain structures (i.e., rostrum breadth and height, separation between bullae). The skull was directly compared to that of extant flying squirrels housed within the collections of the Naturalis Biodiversity Center of Leiden (Aeromys tephromelas; Belomys pearsonii; Eupetaurus cinereus; Glaucomys sabrinus; Glaucomys volans; Hylopetes alboniger; Hylopetes phayrei; Hylopetes platyurus; Hylopetes sagitta; Iomys horsfieldii; Petaurista alborufus; Petaurista elegans; Petaurista petaurista; Petaurista philippensis; Petinomys hageni; Petinomys setosus; Petinomys vordermanni; Pteromys momonga; Pteromys volans) and the Muséum National d’Histoire Naturelle of Paris (Iomys horsfieldii, Petaurillus kinlochii, Pteromyscus pulverulentus, Trogopterus xantiphes). Only Eoglaucomys and the poorly-known genera Aeretes and Biswamayopterus could not be directly compared. Comparisons with Aeretes melanopterus are based on published data (Thorington et al., 2002; Thorington et al., 2012). In addition, the posterior region of the skull of P. petaurista was scanned in order to compare its inner ear morphology to that of M. neogrivensis (see Materials and Methods, Video 4). Several cranial and mandibular measurements were taken in the specimens of M. neogrivensis and some comparative material (see Supplementary file 5–6). Measurements were also taken from the virtually reconstructed cranium of M. neogrivensis using Avizo 7.0.
Cranium. The cranium is relatively short and wide, as in most large-sized flying squirrels but Eupetaurus and a few species of Petaurista, such as Petaurista alborufus (Figure 6). The upper incisors are orthodont, as in most flying squirrels except Pteromys and Aeretes, which show opisthodont incisors. The incisors are not as wide and robust as in other large-sized flying squirrels. The rostrum is similar in length to that of Petaurista, Aeretes and Aeromys. It is also relatively wide and high, although not as much as in Aeretes. The nasals are not preserved in any of the specimens, but their suture with the premaxilla is almost complete in IPS88677 (Figure 6—figure supplement 2), showing that they were most likely wide, resembling those in Aeretes and Petaurista. The incisive foramina are slightly longer than in most flying squirrels, approaching the condition seen in Aeromys and Eupetaurus. The masseteric tubercles are somewhat prominent and bulbous, being very similar in size and shape to those of Petaurista. In addition, there is a prominent semicircular scar for the origin of the superficial masseter (Ball and Roth, 1995; Thorington and Darrow, 1996) in front of the cheek tooth row. The palate is not particularly wide, being markedly narrower than in Glaucomys, Pteromys and most Hylopetes and Petinomys species. A short posterior nasal spine is present at the suture between the palatines. This is highly variable in flying squirrels, but the spine is generally reduced in Aeromys, Hylopetes and Glaucomys. The pterygoid ridges diverge distally, as in Aeretes, instead of being parallel as in most flying squirrels or slightly convergent as in Pteromys and Trogopterus. The auditory bullae are inflated and similar in size to those of Petaurista and Aeromys and separated by a wide distance as in Petaurista and most Hylopetes species. This is in marked contrast with Aeromys and, especially, Aeretes, which show a narrower space between the bullae. A μCT scan of the middle ear cavity of IPS88677 shows the presence of two ventral transbullar septa (Video 4), while this feature is not preserved in IPS56468h. Septated bullae occur in the tympanic cavity of certain rodents (Aplodontiidae, Sciuridae, Geomyidae and some Arvicolinae), although their function is unclear (Mason, 2015). Several authors have stressed the taxonomic importance of the presence or absence of these structures and the number of septa (Moore, 1959; Pfaff et al., 2015; Thorington et al., 2002). Most squirrels show transbullar septa, and their number typically ranges from one to three, but some flying squirrels show many more, sometimes even arranged in a honeycomb (as in Petinomys) or in a sliver pattern (as in Trogopterus and Pteromyscus) (Moore, 1959; Pfaff et al., 2015; Thorington et al., 2002). Miopetaurista neogrivensis shows just two septa, as in most flying squirrels including Petaurista (Video 4). The external auditory meatus is large and rounded. The bony process from the opening ear canal is large. As in all flying squirrels but Iomys and Petaurillus, the mastoid process is not inflated.
In dorsal view the postorbital process of the frontal is robust and long, partly enclosing the distal half of the orbit, thus approaching the morphology of large-sized flying squirrels such as Petaurista and Aeromys. Small-sized flying squirrels usually show more slender and shorter postorbital processes. The frontal is depressed on its anterior part. There is a conspicuous semicircular notch just in front of the postorbital processes. This notch is also present in many flying squirrels, including Aeromys and Petaurista, but it is more pronounced in Miopetaurista, thus approaching the morphology of Pteromys, Glaucomys and most Hylopetes species. The zygomatic plate of the premaxillary bone is similar in size and morphology to that of most flying squirrels, except for Pteromys and Petaurillus, which show wider plates bordered by more marked ridges. The infraorbital foramen is relatively wide. The zygomatic arches are wide and approximately parallel to the longitudinal axis of the skull. The jugal is deeper than in most flying squirrels and shows a dorsal process pointing towards the postorbital process of the frontal. This process is also observed in many flying squirrel genera, including Petaurista, Aeromys and Aeretes, but it is more pointed than in Miopetaurista. The origin of the deep masseter on the jugal bone (Ball and Roth, 1995; Thorington and Darrow, 1996) defines a marked scar that is larger and extends somewhat more dorsally than in any other studied squirrel. The cranial vault is relatively flat as in most large-sized flying squirrels, except for Aeromys tephromelas and some Petaurista species. Small-sized flying squirrels generally show more convex cranial vaults, sometimes markedly inflated as in Iomys. The parietal is wide and shows well defined temporal ridges (marking the insertion of the temporalis muscle (Ball and Roth, 1995; Thorington and Darrow, 1996) that approach each other towards the posterior part of the skull. Even though this occurs to some degree in certain flying squirrels (Pteromys, Hylopetes), the ridges come so close only in Eupetaurus.
Mandible. Both hemimandibles of IPS56468 are damaged in their posterior part, with the mandibular process partially broken and displaced (Figure 3a–e). The articular process is completely missing in the left hemimandible (IPS56468j), which also lacks the tip of the coronoid process. On the contrary, both processes are preserved in the right hemimandible (IPS56468i), although detached and moved from their anatomical location. The partial hemimandible IPS87560 preserves the articular process in its original position as well as the angular process (Figure 3f–g). A complete mandible was assembled using three-dimensional models of these three specimens digitized by means of photogrammetry (see Materials and Methods; Figure 1, Videos 1 and 3D model in Supplementary file 1). However, the description of mandibular morphology relies on the original specimens, the three-dimensional model being only considered for assessing the shape of the mandibular ramus between the coronoid and articular processes.
The mandible is relatively deep and short, with the lower incisor projecting slightly above the occlusal plane of the cheek teeth. The corpus is strong and deep. The symphysis is strong and the diastema is relatively short. The masseteric ridges converge at the level of the distal root of the m1. The point where the ridges meet does not form a bump or a tubercle as in other sciurids. The masseteric fossa is quite deep. The coronoid process is robust and hook-shaped, resembling that of Aeretes, Aeromys, Petaurista, Belomys, Petinomys and some Hylopetes species such as Hylopetes sagitta. The coronoid is also noticeably higher than the articular process as in all these genera, but rises in an angle closer to that of Hylopetes and Petinomys and not as steep as in Petaurista, Aeretes, Belomys or Trogopterus. IPS56468i partially preserves the outline of the ramus between the coronoid and articular processes. These two processes are not as separated as in Petaurista, Aeretes, Belomys or Trogopterus, being closer as in Aeromys, Glaucomys and some species of Hylopetes and Petinomys. IPS87560 preserves the articular process in its original position and shows that it is long and projects distally. Aeretes, Glaucomys and Hylopetes show similarly long articular processes, while most large-sized flying squirrels such as Aeromys, Eupetaurus and Petaurista show shorter processes. The articular condyle is longitudinally elongated. The angular process is of intermediate length and distally directed, not extending far ventrally. As in Petaurista, Aeromys and Hylopetes, for example, the angular process does not reach farther posteriorly than the articular one. The inner surface of the angular process shows a marked fossa for the insertion of the medial pterygoid partially divided by a marked longitudinal ridge. The ventral edge of the angular process presents a well-defined area for the insertion of the superficial masseter bordered by faint ridges. Overall, mandibular morphology and robusticity recall Aeromys and Hylopetes, being slightly more gracile and elongated than in Aeretes and Petaurista, for example.
The skull of Miopetaurista resembles that of other large-sized flying squirrels, particularly Aeromys and Petaurista, which are characterized by their short and wide rostrum, moderately inflated bullae and relatively wide posterior region of the skull. Other details are strikingly similar to these extant genera, including the robust and long postorbital process that partially encloses the orbit, the well-developed jugal process in the zygomatic arch, or the presence of two septa in the tympanic cavity. Most of the smaller flying squirrels show more elongate muzzles, slender or shorter postorbital processes and, in some cases, a higher number of transbullar septa. However, the mandible of Miopetaurista is somewhat slenderer than in most large-sized flying squirrels, showing a more elongated articular process resembling that of Hylopetes and other small-sized flying squirrels instead of the lower and shorter process of Petaurista. The characters of the skull agree with our phylogenetic results (see main text and Figure 7) in showing a close relationship between Miopetaursita and extant large-sized flying squirrels, particularly Petaurista. However, we must note that none of these characters distinguishes flying squirrels exclusively. Indeed, overall skull morphology of some large-sized tree squirrels such as Ratufa and Protoxerus (Ellerman, 1940; Moore, 1959) is similar to that of large-sized flying squirrels.
3.3 Postcranial morphology
The partial skeleton IPS56468 comprises up to 80 different elements (Table 1) including the skull, most of the long bones of the fore- and hindlimbs, elements of the wrist and ankle, metapodials and phalanges, fragmentary pelves and clavicles, a few vertebrae from different regions and scarce rib fragments. Some of these elements, particularly those from the hindlimb, are remarkably well preserved and were even found in anatomical connection (Figure 2 and Figure 2—figure supplement 1). Here we briefly describe the postcranial skeleton, outlining the distinguishing morphological features of flying squirrels (Thorington et al., 2005), which are mostly related to gliding locomotion. The skeletal elements of M. neogrivensis were compared with those of extant squirrels (Supplementary file 7) in the collections of the Institut Català de Paleontologia Miquel Crusafont (Glaucomys sabrinus, Sciurus niger) and the Naturalis Biodiversity Center (Urocitellus parryii, Xerus erythropus, Callosciurus prevostii, Sciurus vulgaris, Hylopetes sagitta and Petaurista petaurista).
The recovered long bones include the two humeri and femora, as well as the right tibia. Only the distal half of the left tibia and that of the right radius have been preserved, the distal epiphysis of the radius further being damaged. Complete long bones of flying squirrels are much longer and more slender than those of tree and ground squirrels (Samuels and Van Valkenburgh, 2008; Thorington et al., 2005; Thorington and Heaney, 1981; Thorington and Santana, 2007), and this is also the case for M. neogrivensis (Figure 5—figure supplement 1). Moreover, several additional features can be recognized in the proximal and distal ends of these bones as described below.
Humerus. On the proximal half of the humerus, the deltoid ridge is not as prominent as in tree and specially ground squirrels (Figure 5—figure supplement 1a-e). In flying squirrels this ridge is also sharper and much thinner than in ground squirrels and most tree squirrels. The deltoid ridge in M. neogrivensis is less prominent than in Hylopetes, being roughly comparable to that of Petaurista. The reduced deltoid ridge of flying squirrels probably allows increased flexibility at the shoulder joint, as opposed to the condition of tree and ground squirrels, which serves to increase the power of forelimb flexion at the expense of flexibility at this joint (Thorington et al., 2005). The distal end of the humerus shows additional differences between the three kinds of squirrels. The epicondylar ridge is large in tree and particularly in ground squirrels, extending to the lateral edge of the epicondyle (Figure 5—figure supplement 1f-g). In most flying squirrels, this ridge is reduced and does not reach the epicondyle (Figure 5—figure supplement 1h-j). In large-sized flying squirrels and M. neogrivensis this ridge is even more reduced than in small-sized flying squirrels such as Hylopetes and Eoglaucomys. Flying squirrels and M. neogrivensis are further characterized by a small medial epicondyle that contrasts with the much larger one of ground squirrels and the elongated epicondyle of tree squirrels (Figure 5—figure supplement 1f-j). The reduction of the lateral epicondylar ridge in flying squirrels represents a reduction in the leverage of the extensor muscles for flexion of the elbow joint, thus facilitating elbow joint extension for gliding (Thorington et al., 2005). Finally, the trochlea is not sharply angled as in tree and ground squirrels, being flatter as in flying squirrels (Figure 5—figure supplement 1f-j).
Radius. The radius is not well preserved, so that the very few features that distinguish the radius of flying squirrel from that of other squirrels cannot be observed.
Carpal bones. Carpals are the most diagnostic elements of pteromyins and can be used to unambiguously diagnose the group because of their role in the extension of the patagium (Thorington, 1984; Thorington et al., 2005; Thorington et al., 2002; Thorington and Darrow, 2000). The matrix removed during the preparation of the two blocks containing the skeleton of M. neogrivensis was screen-washed, allowing for the recovery of small elements, including the complete right scapholunate and the dorsal end of the pisiform (Figure 4 and Video 2). The pisiform shows three processes for articulation with the ulna, triquetrum and scapholunate. The presence of a radial thumb-like process that articulates with the scapholunate is a pteromyin synapomorphy not observed in other squirrels. In M. neogrivensis, as well as in other flying squirrels, the scapholunate process of the pisiform is well developed and fits into a distinct facet of the former bone (Figure 4 and Video 2). The patagium in flying squirrels is supported by the styliform cartilage, which is strongly bound to the palmar end of the pisiform, and extended and raised by means of the abductor of the thumb on the opposite side of the wrist (Thorington, 1984; Thorington et al., 1998; Thorington and Darrow, 2000). The articulation of the pisiform with both the triquetrum and the scapholunate causes the proximal row of carpal bones to function as a single unit, providing a more stable base for the styliform cartilage and the wing tip during gliding (Thorington, 1984; Thorington et al., 1998; Thorington and Darrow, 2000). On the other hand, the articular surface between the radius and scapholunate is more spherical in flying squirrels and M. neogrivensis, whereas it is planoconvex in tree squirrels (Thorington, 1984; Thorington and Darrow, 2000). Maximum extension of the wing tip of the patagium is achieved when the hand is dorsiflexed and radially deviated (Thorington, 1984; Thorington et al., 1998; Thorington and Darrow, 2000). Dorsiflexion is required to elevate the tip of the patagium at the wrist relative to the rest of the membrane and it would likely function to reduce induced drag when gliding (Thorington et al., 1998). The greater convexity and extension of the dorsal articular surface of the scapholunate in the pteromyins allow for a greater angular movement at the proximal wrist joint, which is required to bring the hands to gliding position. Besides the scapholunate and ulnar process, the pisiform of M. neogrivensis shows a distinct triquetral process (Figure 4 and Video 2). This process is diagnostic of the Pteromyina (i.e., the large-sized flying squirrel clade, see main text) and fits between the palmar surfaces of the scapholunate and triquetrum (Thorington et al., 2002; Thorington and Darrow, 2000). Therefore, the proximal carpal bones of M. neogrivensis not only unambiguously indicate that it is a flying squirrel, but also allow assigning it to one of the two extant flying squirrel clades (Figure 7).
Femur. The proximal end of the femur shows reduced lesser trochanter and third trochanter, as seen in flying squirrels (Figure 5—figure supplement 1k-o). The lesser trochanter is larger in most tree squirrels, extending medially as far or farther than the femoral head. In marmotin ground squirrels, the lesser trochanter is even shorter than in flying squirrels (Figure 5—figure supplement 1k), but it is more robust and posteromedially directed, being presumably linked to their digging activities and burrowing habits. In flying squirrels, the third trochanter is markedly more reduced than in other squirrels and it is placed more proximally. It has been proposed that the reduction of the third trochanter and the trochanteric ridge in flying squirrels may increase flexibility at the hip joint (Thorington et al., 2005), as described for the case of the shoulder joint. Although the distal epiphysis shows fewer diagnostic characters, M. neogrivensis resembles large-sized flying squirrels (Petaurista, Aeromys) in showing a broader patellar groove (Figure 5f–j).
Tibia. The popliteal fossa is shallow and not limited by low ridges as in tree and small flying squirrels (Figure 5—figure supplement 1t-s). This fossa is generally deeper in tree squirrels, and markedly pronounced in ground squirrels. The anterior crest is less pronounced than in tree and ground squirrels (Figure 5—figure supplement 1t-s). The anterior surface of the distal epiphysis does not present the tuberosity for the origin of the tibiocarpalis muscle that characterizes the Glaucomyina (i.e., the small-sized flying squirrel clade, see main text; Figure 5h–j; Thorington et al., 2002). The tibiocarpalis runs from the ankle to the tip of the styliform cartilage and defines the margin of the flying membrane. In the Pteromyina, and presumably in M. neogrivensis as well, this muscle originates from the dorsal surface of the first or second metatarsals (Thorington et al., 2002). This further supports the inclusion of M. neogrivensis among the latter subtribe.
Astragalus. The trochlear surface of the astragalus is said to be more V-shaped in flying squirrels than in other squirrels (Thorington et al., 2005). In the extant sciurids we examined, the shape of the trochlea is similar in tree and flying squirrels, although in small-sized flying squirrels (Hylopetes and Glaucomys) the angle between the lateral and medial malleolar surfaces in the trochlea is sharper. The plantar surface of the astragalus shows a well-defined sulcus tali between the posterior and the middle calcaneal articular surfaces (Figure 2k–l). This sulcus accommodates the astragalocalcaneal interosseus ligament and synovial tissue of the joints. In flying squirrels the two calcaneal articular surfaces are enlarged and occlude the sulcus to the point that it is barely appreciated (Ginot et al., 2016; Thorington et al., 2005). This derived morphology has been related to a greater mobility of the ankle joint associated with extreme plantar inversion of the foot during gliding (Thorington et al., 2005). A greatly reduced sulcus is also present in the tree squirrels of the tribes Sciurini (Sciurus, Tamiasciurus) and Protoxerini (Funisciurus, Paraxerus), presumably permitting the reversal of the foot during headfirst descent of trees (Ginot et al., 2016; Thorington et al., 2005). On the contrary, this sulcus is really prominent in Ratufa (tribe Ratufini) and at least in some Callosciurini (Callosciurus, Sundasciurus) (Thorington et al., 2005). However, all these tree squirrels are perfectly capable of rotating their feet. A similarly developed sulcus is present in Douglassciurus jeffersoni and has been interpreted as a plesiomorphic feature unrelated to tree climbing (Emry and Thorington, 1982; Emry and Thorington, 1984). On the other hand, in marmotine ground squirrels, which do not climb trees, the sulcus is very narrow or completely missing (Thorington et al., 2005). Miopetaurista neogrivensis shows a deep groove comparable to that of Ratufa. Interestingly, the astragalus of the early Miocene North American putative flying squirrel Petauristodon pattersoni presents a similarly pronounced sulcus (Pratt and Morgan, 1989) which has been previously taken as a proof for the exclusion of this genus from pteromyins (Thorington et al., 2005). Clearly, the presence of an occluded sulcus in the plantar surface of the astragalus is not diagnostic of flying squirrels and is not required for foot inversion.
Calcaneus. The sustentacular process of the calcaneus is clearly larger than the fibular process and it is placed slightly more distally (Figure 2m–n), as observed in tree and flying squirrels (Ginot et al., 2016; Thorington et al., 2005). In ground squirrels both processes are aligned opposite to one another and the sustentacular process is only slightly larger than the fibular one, resulting in a cruciform calcaneus. As in flying squirrels and some tree squirrels of the tribes Protoxerini, Ratufini and Callosciurini, M. neogrivensis shows a pronounced calcaneal sulcus between the two facets for articulation with the astragalus. The morphology of the calcaneus is diagnostic for ground squirrels but does not appear to be of use in discriminating between tree and flying squirrels.
Vertebrae. Several vertebrae of M. neogrivensis have been recovered (Figure 2g–h, Table 1), including a partial axis, a partial cervical (C3), four partial thoracics (T1?–T4?), four lumbars in anatomical connection (L3I–L5; Figure 2g–h) and four caudals that correspond to the mid part of the tail. Small-sized flying squirrels of the genus Glaucomys have more elongated lumbar and mid-caudal vertebrae as compared to similarly sized tree squirrels (Microsciurus) and show a reduction of the transverse processes in lumbar vertebrae (Thorington and Santana, 2007). The lumbar vertebrae of Petaurista are also slenderer than in tree squirrels and present reduced transverse processes oriented parallel to the vertebral body. The four recovered lumbar vertebrae of M. neogrivensis show the characteristic morphology of Petaurista. The elongation of lumbar vertebrae and limbs determines the size and shape of the patagium, which in turn dictates important aerodynamic features such as wing loading and aspect ratio (Thorington and Heaney, 1981; Thorington and Santana, 2007). The elongation of these elements results in a decrease in wing loading of flying squirrels.
Clavicle. Two clavicles have been recovered, the right one missing only the articulation with the scapula. The clavicle is straight and slender as in arboreal and flying squirrels, while in ground squirrels this bone is conspicuously robust.
The remaining recovered bones (Table 1) include a partial pelvis (Figure 2i–j), some rib fragments and two sternebrae. In addition, several elements of the ankle and foot have been recovered, including the navicular, two cuneiforms, and several metatarsals, with left metatarsals 2 – 4 in anatomical connection. Most phalanges were recovered after screen-washing and are difficult to assign to ray. Four proximal hand phalanges as well as four proximal foot phalanges have been recovered. Seven intermediate phalanges and several phalangeal fragments are also attributed to the foot. Metapodials as well as proximal and intermediate phalanges are slender as in tree and flying squirrels. Six distal phalanges (tentatively attributed to the foot) have also been recovered. These are markedly sharp and curved as in most tree and flying squirrels, instead of blunt as in many ground squirrels.
Data availability
3D surface models of the described material of M. neogrivensis (reconstructed skeleton, reconstructed skull, carpal bones) are available at MorphoBank with identifier http://morphobank.org/permalink/?P3108. GenBank accession numbers for all the extant taxa considered in the phylogenetic analyses are provided in Supplementary File 2. All other data files are included in the supplementary information.
References
-
TranslatorX: multiple alignment of nucleotide sequences guided by amino acid translationsNucleic Acids Research 38:W7–W13.https://doi.org/10.1093/nar/gkq291
-
Molecular phylogeny and divergence time estimates for Major rodent groups: evidence from multiple genesMolecular Biology and Evolution 18:777–791.https://doi.org/10.1093/oxfordjournals.molbev.a003860
-
BookRoedores Esciuromorfos Del Neógeno De CataluñaBellaterra: Universitat Autònoma de Barcelona.
-
Los esciurópteros del mioceno de la cuenca del Vallès-Penedès (Cataluña, España)Geogaceta 11:114–116.
-
The late oligocene rodent faunas of canales (MP28) and parrales (MP29) from the Loranca Basin, province of Cuenca, SpainRevista Española De Paleontología 14:93–116.
-
Genetic data reveal a cryptic species of new world flying squirrel: glaucomys oregonensisJournal of Mammalogy 98:1027–1041.https://doi.org/10.1093/jmammal/gyx055
-
Jaw muscles of new world squirrelsJournal of Morphology 224:265–291.https://doi.org/10.1002/jmor.1052240303
-
A complete terrestrial oligocene magnetobiostratigraphy from the Ebro Basin, SpainEarth and Planetary Science Letters 187:1–16.https://doi.org/10.1016/S0012-821X(01)00270-9
-
BookAplodontid, Castorid, Zapodid and Geomyoid Rodents of the Rodent Hill Locality, Cypress Hills Formation, South-West Saskatchewan (Masters Thesis)Saskatoon: University of Saskatchewan.
-
Cranial dimensions as estimators of body mass and locomotor habits in extant and fossil rodentsJournal of Vertebrate Paleontology 36:e1014905.https://doi.org/10.1080/02724634.2015.1014905
-
A review of the North American tertiary SciuridaeBulletin of the Museum of Comparative Zoology at Harvard College 130:109–248.
-
Rodents and lagomorphs from the middle miocene hominoid-bearing site of barranc de can Vila 1 (els hostalets de pierola, Catalonia, Spain)Neues Jahrbuch Für Geologie Und Paläontologie - Abhandlungen 257:297–315.https://doi.org/10.1127/0077-7749/2010/0077
-
Late miocene flying squirrels from can llobateres 1 (Vallès-Penedès basin, Catalonia): systematics and palaeobiogeographyPalaeobiodiversity and Palaeoenvironments 95:353–372.https://doi.org/10.1007/s12549-015-0192-1
-
The miocene mammal record of the Vallès-Penedès basin (Catalonia)Comptes Rendus Palevol 15:791–812.https://doi.org/10.1016/j.crpv.2015.07.004
-
Selection of conserved blocks from multiple alignments for their use in phylogenetic analysisMolecular Biology and Evolution 17:540–552.https://doi.org/10.1093/oxfordjournals.molbev.a026334
-
Los Sciuridae (Rodentia, mammalia) del oligoceno inferior de montalbán y olalla (Teruel, españa). Discusión sobre el origen de los esciúridosBoletín De La Sociedad Española De Historia Natural Sección Geológica 87:155–169.
-
A classification of the Gliridae (Rodentia) on the basis of dental morphologyHystrix : Rivista Di Teriologia 6:3–50.
-
Distribution and Migration of Tertiary Mammals in Eurasia. Volume in Honour of Hansde Bruijn-Deinsea97–126, Paleogene rodents of Eurasia, Distribution and Migration of Tertiary Mammals in Eurasia. Volume in Honour of Hansde Bruijn-Deinsea, 10.
-
Taxonomische probleme um das GenusMiopetaurista kretzoi, 1962 (Fam. Sciuridae)Paläontologische Zeitschrift 49:75–77.https://doi.org/10.1007/BF02988067
-
Biber und ein zwerghamster aus mataschen (Unter-Pannonium, Steirisches Becken)Joannea Geologie und Paläontologie 5:19–33.
-
BookCatalogus Fossilium Austriae Band 4 - Rodentia Neogenica (First edition)Springer Verlag der Österreichischen Akademie der Wissenschaften.
-
Petauristinae (Mammalia, Rodentia) from the Oligocene of Spain, Belgium and Turkish ThraceNatural History Museum of Los Angeles Scholarly Series 33:139–146.
-
BookSciuridae, petauristidae and eomyidae (Rodentia, mammalia)In: Schmidt-Kittler N, editors. The Vertebrate Locality Maramena (Macedonia, Greece) at the Turolian-Ruscinian Boundary (Neogene). Munich: Verlag Dr. Friedrich Pfeil. pp. 87–102.
-
Vertebrates from the early miocene lignite deposits of opencast mine oberdorf (Western styrian basin, Austria).6 Rodentia 1 (Mammalia). Ann Naturhist Mus Wien 99 A:99–137.
-
BookSuperfamily SciuroideaIn: Rössner G. E, Heissig K, editors. The Miocene Land Mammals of Europe. Munich: Verlag Dr. Friedrich Pfie. pp. 271–280.
-
Bayesian phylogenetics with BEAUti and the BEAST 1.7Molecular Biology and Evolution 29:1969–1973.https://doi.org/10.1093/molbev/mss075
-
BookThe Families and Genera of Living Rodents with a List of Named Forms (1758-1936) (Reprint edition)British Museum.
-
Descriptive and comparative osteology of the oldest fossil squirrel protosciurus (Rodentia: Sciuridae). Smithson. ContribPaleobiology 47:1–35.
-
BookThe tree squirrel sciurus (Sciuridae, Rodentia) as a living fossilIn: Elredge N, Stanley S. M, editors. Living Fossils. New York: Springer-Verlag. pp. 23–31.https://doi.org/10.1007/978-1-4613-8271-3_3
-
A tree squirrel (Rodentia, Sciuridae, sciurini) from the late miocene (Clarendonian) of NevadaJournal of Vertebrate Paleontology 25:228–235.https://doi.org/10.1671/0272-4634(2005)025[0228:ATSRSS]2.0.CO;2
-
A new genus of squirrel (Rodentia, Sciuridae) from the mid-Cenozoic of north americaJournal of Vertebrate Paleontology 27:693–698.https://doi.org/10.1671/0272-4634(2007)27[693:ANGOSR]2.0.CO;2
-
Latest oligocene didelphimorphia, rodentia and lagomorpha (Mammalia) from Oberleichtersbach, rhön mountains, GermanyCour Forsch-Inst Senckenberg 260:185–251.
-
A glimpse on the pattern of rodent diversification: a phylogenetic approachBMC Evolutionary Biology 12:88.https://doi.org/10.1186/1471-2148-12-88
-
BookLatitudinal climatic gradients in the Western European and Mediterranean regions from the Mid-Mioce ne (c. 15 Ma) to the Mid-Pliocene (c. 3.5 Ma) as quantified from pollen dataIn: Williams M, Haywood A. M, Gregory F. M, Schmidt D. N, editors. Deep-Time Perspectives on Climate Change: Marrying the Signal From Computer Models and Biological Proxies, the Micropalaeontological Society, Special Publications. London: The Geological Society. pp. 481–502.
-
BookEvolution of Tertiary Mammals of North America: Volume 2, Small Mammals, Xenarthrans, and Marine MammalsIn: Janis C. M, Gunnell G. F, Uhen M. D, editors. Aplodontoidea. New York: Cambridge University Press. pp. 377–390.
-
On some Miocene squirrels, with remarks on the dentition and classification of the SciurinaeProc. Zool. Soc. Lond 1893:179–215.
-
BookEvolution of Tertiary Mammals of North America: Volume 2, Small Mammals, Xenarthrans, and Marine MammalsIn: Janis C. M, Gunnell G. F, Uhen M. D, editors. Sciuridae. New York: Cambridge University Press. pp. 355–376.
-
Patterns of body and tail length and body mass in SciuridaeJournal of Mammalogy 89:852–873.https://doi.org/10.1644/07-MAMM-A-217.1
-
The fossilized birth-death process for coherent calibration of divergence-time estimatesProceedings of the National Academy of Sciences 111:E2957–E2966.https://doi.org/10.1073/pnas.1319091111
-
Die frühesten flughörnchen und primitive ailuravinae (Rodentia, mamm) aus dem süddeutschen oligozän.Mitteilungen der Bayerischen Staatssammlung für Paläontologie und Histor. Geologie. 19:139–169.
-
BookGliding Mammals: Taxonomy of Living and Extinct Species, Smithsonian Contributions to ZoologyWashington: Smithsonian Institution Scholarly Press.
-
BookGlires SummaryIn: Christine M, Gunnell G. F, Uhen M. D, editors. Evolution of Tertiary Mammals of North America - Volume 2: Small Mammals, Xenarthrans, and Marine Mammals. New York: Cambridge University Press. pp. 263–292.
-
MAFFT: a novel method for rapid multiple sequence alignment based on fast fourier transformNucleic Acids Research 30:3059–3066.https://doi.org/10.1093/nar/gkf436
-
BookFamily Sciuridae (tree, flying and ground squirrels, chipmunks, marmots and praire dogs)In: Wilson D. E, Lacher T. E, Mittermeier R. A, editors. Handbook of the Mammals of the World. 6. Lagomorphs and Rodents I, Handbook of Mammals of the World. Barcelona: Lynx Edicions. pp. 648–837.
-
The skull of cedromus and a review of the cedromurinae (Rodentia, Sciuridae)Journal of Paleontology 65:984–994.https://doi.org/10.1017/S0022336000033291
-
First record of Hesperopetes Emry and Korth, 2007 (Rodentia, Sciuridae), from the Orellan (early Oligocene)Journal of Vertebrate Paleontology 37:e1318890.https://doi.org/10.1080/02724634.2017.1318890
-
Partitionfinder: combined selection of partitioning schemes and substitution models for phylogenetic analysesMolecular Biology and Evolution 29:1695–1701.https://doi.org/10.1093/molbev/mss020
-
Erlianomys combinatus, a primitive myodont rodent from the eocene arshanto formation, Nuhetingboerhe, nei mongol, ChinaVertebrata PalAsiatica 48:133–144.
-
BookEvolution of the Rodents: Advances in Phylogeny, Functional Morphology, and Development, Cambridge Studies in Morphology and MoleculesIn: Cox P. G, Hautier L, editors. Functional Morphology of Rodent Middle Ears (Eds. United Kingdom: Cambridge University Press. pp. 373–404.
-
BookProtrogomorph, Sciuromorph, Castorimorph, Myomorph (Geomyoid, Anomaluroid, Pedetoid and Ctenodactyloid) RodentsIn: Anderson S, Jones J. K, editors. Orders and Families of Recent Mammals of the World. New York: John Wiley & Sons Inc. pp. 267–288.
-
Molecular systematics of sciurognathi (rodentia): the mitochondrial cytochrome b and 12S rRNA genes support the anomaluroidea (Pedetidae and Anomaluridae)Molecular Phylogenetics and Evolution 22:220–233.https://doi.org/10.1006/mpev.2001.1056
-
Relationships among the living squirrels of the sciurinaeBulletin of the American Museum of Natural History 118:157–205.
-
New Sciuridae (Mammalia: rodentia) from the early miocene thomas farm local fauna, FloridaJournal of Vertebrate Paleontology 9:89–100.https://doi.org/10.1080/02724634.1989.10011741
-
A new sciuroptere from middle Miocene of Linqu, ShandongVertebrata PalAsiatica 19:228–238.
-
The aragonian vertebrate fauna of Xiacaowan, Jiangsu -5. Sciuridae (Rodentia, Mammalia)Vertebrata PalAsiatica 24:195–209.
-
Sciurids from the Late Miocene Lufeng hominoid loality, YunnanVertebrata PalAsiatica 40:177–193.
-
BookNeogene faunal succession and biochronology of central Nei Mongol (Inner Mongolia) (Eds).In: Wang X, Flynn L. J, Fortelius M, editors. Fossil Mammals of Asia: Neogene Biostratigraphy and Chronology. New York: Columbia University Press. pp. 155–186.
-
BookNeogene land mammal stages/ages of China. Toward the goal to establish an Asian land mammal stage/age schemeIn: Flynn L. J, Fortelius M, editors. Fossil Mammals of Asia: Neogene Biostratigraphy and Chronology. . New York: Columbia University Press. pp. 29–90.
-
Latest Pliocene mammals from Haile XV A, Alachua County, FloridaBulletin of the Florida Museum of Natural History. Biological Sciences 20:111–186.
-
Closing the gap between rocks and clocks using total-evidence datingPhilosophical Transactions of the Royal Society B: Biological Sciences 371:20150136.https://doi.org/10.1098/rstb.2015.0136
-
Early irvingtonian (Latest pliocene) rodents from inglis 1C, Citrus county, FloridaJournal of Vertebrate Paleontology 21:153–171.https://doi.org/10.1671/0272-4634(2001)021[0153:EILPRF]2.0.CO;2
-
Long bone articular and diaphyseal structure in old world monkeys and apes. II: estimation of body massAmerican Journal of Physical Anthropology 120:16–37.https://doi.org/10.1002/ajpa.10118
-
Skeletal indicators of locomotor adaptations in living and extinct rodentsJournal of Morphology 269:1387–1411.https://doi.org/10.1002/jmor.10662
-
The principles of classification and a classification of mammalsBulletin of the American Museum of Natural History 85:1–350.
-
A new “flying squirrel” (Rodentia: Sciuridae) from the Early Miocene of southwestern SaskatchewanJournal of Vertebrate Paleontology 6:290–294.https://doi.org/10.1080/02724634.1986.10011625
-
Bias in equations used to estimate fossil primate body massJournal of Human Evolution 25:31–41.https://doi.org/10.1006/jhev.1993.1036
-
Regression models for prediction equationsJournal of Human Evolution 26:239–244.https://doi.org/10.1006/jhev.1994.1013
-
A proposal for a standard terminology of anatomical notation and orientation in fossil vertebrate dentitionsJournal of Vertebrate Paleontology 23:1–12.https://doi.org/10.1671/0272-4634(2003)23[1:APFAST]2.0.CO;2
-
Sampling-through-time in birth-death treesJournal of Theoretical Biology 267:396–404.https://doi.org/10.1016/j.jtbi.2010.09.010
-
Nuclear DNA phylogeny of the squirrels (Mammalia: rodentia) and the evolution of arboreality from c-myc and RAG1Molecular Phylogenetics and Evolution 30:703–719.https://doi.org/10.1016/S1055-7903(03)00204-5
-
Extraordinarily preserved specimen of the oldest known glirid from the middle eocene of messel (Rodentia)Journal of Vertebrate Paleontology 27:189–194.https://doi.org/10.1671/0272-4634(2007)27[189:EPSOTO]2.0.CO;2
-
Evolution of grasses and grassland ecosystemsAnnual Review of Earth and Planetary Sciences 39:517–544.https://doi.org/10.1146/annurev-earth-040809-152402
-
Body Proportions and Gliding Adaptations of Flying Squirrels (Petauristinae)Journal of Mammalogy 62:101–114.https://doi.org/10.2307/1380481
-
Flying squirrels are monophyleticScience 225:1048–1050.https://doi.org/10.1126/science.225.4666.1048
-
Jaw muscles of old world squirrelsJournal of Morphology 230:145–165.https://doi.org/10.1002/(SICI)1097-4687(199611)230:2<145::AID-JMOR3>3.0.CO;2-G
-
Wing tip anatomy and aerodynamics in flying squirrelsJournal of Mammalogy 79:245–250.https://doi.org/10.2307/1382860
-
Phylogenies of flying squirrels (Pteromyinae)Journal of Mammalian Evolution 9:99–135.https://doi.org/10.1023/A:1021335912016
-
Mammal Species of the World. a Taxonomic and Geographic Reference754–818, Family Sciuridae, Mammal Species of the World. a Taxonomic and Geographic Reference, Baltimore, The Johns Hopkins University Press.
-
The difficulties of identifying flying squirrels (Sciuridae: pteromyini) in the fossil recordJournal of Vertebrate Paleontology 25:950–961.https://doi.org/10.1671/0272-4634(2005)025[0950:TDOIFS]2.0.CO;2
-
How to make a flying squirrel: Glaucomys anatomy in phylogenetic perspectiveJournal of Mammalogy 88:882–896.https://doi.org/10.1644/06-MAMM-S-325R2.1
-
BookThe Paleogene perioIn: Gradstein F. M, Ogg J. G, Schmitz M, Ogg G, editors. The Geologic Time Scale . Amsterdam: Elsevier. pp. 855–921.https://doi.org/10.1016/B978-0-444-59425-9.00028-7
-
Palaeosciurus goti nov. sp., écureuil terrestre de l’Oligocène moyen du Quercy: donnés nouvelles sur l’apparition des Sciuridés en EuropeAnnales De Paléontologie pp. 103–122.
-
Problèmes Actuels Paléontol. Évolution Vertébrés765–776, Caractéristiques évolutives des rongeurs de l’Oligocène d’Europe Occidentale, Problèmes Actuels Paléontol. Évolution Vertébrés, 218, Colloq. Int. CNRS.
-
BookPossible evolutionary relationships among Eocene and lower Oligocene rodents of Asia, Europe and North AmericaIn: Luckett W. P, Hartenberger J. -L, editors. Multidisciplinary Analysis, NATO ASI Series A. New York: Plenum Press. pp. 277–309.
-
BookLate Cretaceous and Cenozoic Mammals of North America: Biostratigraphy and GeochronologyNew York: Columbia University Press.https://doi.org/10.7312/wood13040
-
The woolly flying squirrel and gliding: does size matterActa Theriol Warsz 46:429–435.https://doi.org/10.4098/AT.arch.01-46
-
Total-Evidence dating under the fossilized Birth-Death processSystematic Biology 65:228–249.https://doi.org/10.1093/sysbio/syv080
Article and author information
Author details
Funding
Agencia Estatatal de Investigación and European Regional Development Fund (CGL2017-82654-P)
- Isaac Casanovas-Vilar
- Josep Fortuny
Generalitat de Catalunya (CERCA Programme)
- Isaac Casanovas-Vilar
- Joan Garcia-Porta
- Josep Fortuny
- Marina Querejeta
- Sergio Llácer
- Josep M Robles
- David M Alba
Handel T. Martin Endowment
- Óscar Sanisidro
Agencia Estatal de Investigación and European Regional Development Fund (CGL2016-76431-P)
- Isaac Casanovas-Vilar
- Josep M Robles
- David M Alba
Agencia Estatal de Investigación and European Regional Development Fund (RYC-2013-12470)
- Isaac Casanovas-Vilar
Agencia Estatal de Investigación and European Regional Development Fund (FJCI-2014-20380)
- Joan Garcia-Porta
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This research has been funded by the Spanish Agencia Estatal de Investigación and the European Regional Development Fund of the European Union (CGL2016-76431-P, CGL2017-82654-P, RYC-2013 – 12470 to IC-V. FJCI-2014 – 20380 fellowship to JG-P.), by the Generalitat de Catalunya (CERCA Programme), and by the Handel T Martin Endowment to O S. We further acknowledge the support of the Servei d’Arqueologia i Paleontologia of the Generalitat de Catalunya. Fieldwork was defrayed by CESPA Gestión de Residuos, S A U. We thank M Valls (ICP) for the excellent preparation of the specimens, which was partly funded by the Department of Cultura of the Generalitat de Catalunya (2014/100609). The authors thank L W van den Hoek Ostende and W Renema for allowing access to the Naturalis Biodiversity Center (Leiden, the Netherlands) µCT facility, S van der Mije and P Kaminga for their assistance during the study of the squirrel collections kept at this institution. S Moyà-Solà (ICP), L J Flynn and D Pilbeam (Harvard University) provided valuable discussions on the fossil material and comments on the manuscript. We thank J Madurell-Malapeira (ICP) for his assessment in the preparation of some of the figures and J X Samuels (East Tennessee State University) for providing valuable data regarding the North American ‘flying squirrel’ record. Finally, we thank the editors and reviewers of this work for their valuable and constructive comments.
Copyright
© 2018, Casanovas-Vilar et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 9,449
- views
-
- 558
- downloads
-
- 17
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
Fossils found in Spain have shed light on the evolution of flying squirrels.
-
- Evolutionary Biology
Lineages of rod-shaped bacteria such as Escherichia coli exhibit a temporal decline in elongation rate in a manner comparable to cellular or biological aging. The effect results from the production of asymmetrical daughters, one with a lower elongation rate, by the division of a mother cell. The slower daughter compared to the faster daughter, denoted respectively as the old and new daughters, has more aggregates of damaged proteins and fewer expressed gene products. We have examined further the degree of asymmetry by measuring the density of ribosomes between old and new daughters and between their poles. We found that ribosomes were denser in the new daughter and also in the new pole of the daughters. These ribosome patterns match the ones we previously found for expressed gene products. This outcome suggests that the asymmetry is not likely to result from properties unique to the gene expressed in our previous study, but rather from a more fundamental upstream process affecting the distribution of ribosomal abundance. Because damage aggregates and ribosomes are both more abundant at the poles of E. coli cells, we suggest that competition for space between the two could explain the reduced ribosomal density in old daughters. Using published values for aggregate sizes and the relationship between ribosomal number and elongation rates, we show that the aggregate volumes could in principle displace quantitatively the amount of ribosomes needed to reduce the elongation rate of the old daughters.