Hostname: page-component-cd9895bd7-mkpzs Total loading time: 0 Render date: 2024-12-25T18:32:29.748Z Has data issue: false hasContentIssue false

Diversification trajectories and paleobiogeography of Neogene chondrichthyans from Europe

Published online by Cambridge University Press:  08 February 2023

Jaime A. Villafaña*
Affiliation:
Vienna Doctoral School of Ecology and Evolution, Vienna, Austria; Centro de Investigación en Recursos Naturales y Sustentabilidad, Universidad Bernardo O'Higgins, Avenida Viel 1497, 8370993 Santiago, Chile. E-mail: [email protected]
Marcelo M. Rivadeneira
Affiliation:
Laboratorio de Paleobiología, Centro de Estudios Avanzados en Zonas Áridas (CEAZA), Avenida Bernardo Ossandón 877, 1781681, Coquimbo, Chile; Departamento de Biología Marina, Facultad de Ciencias del Mar, Universidad Católica del Norte, Larrondo 1281, Coquimbo, Chile; Departamento de Biología, Universidad de La Serena, Avenida Raul Bitrán 1305, La Serena, Chile. E-mail: [email protected]
Catalina Pimiento
Affiliation:
Paleontological Institute and Museum, University of Zurich, CH-8006 Zurich, Switzerland; Department of Biosciences, Swansea University, Swansea SA28PP, United Kingdom; Smithsonian Tropical Research Institute, Balboa, Panama. E-mail: [email protected]
Jürgen Kriwet*
Affiliation:
Department of Palaeontology, University of Vienna, Althanstraße 14, Geocenter, 1090 Vienna, Austria. E-mail: [email protected]
*
*Corresponding author.
*Corresponding author.

Abstract

Despite the rich fossil record of Neogene chondrichthyans (chimaeras, sharks, rays, and skates) from Europe, little is known about the macroevolutionary processes that generated their current diversity and geographical distribution. We compiled 4368 Neogene occurrences comprising 102 genera, 41 families, and 12 orders from four European regions (Atlantic, Mediterranean, North Sea, and Paratethys) and evaluated their diversification trajectories and paleobiogeographic patterns. In all regions analyzed, we found that generic richness increased during the early Miocene, then decreased sharply during the middle Miocene in the Paratethys, and moderately during the late Miocene and Pliocene in the Mediterranean and North Seas. Origination rates display the most significant pulses in the early Miocene in all regions. Extinction rate pulses varied across regions, with the Paratethys displaying the most significant pulses during the late Miocene and the Mediterranean and North Seas during the late Miocene and early Pliocene. Overall, up to 27% and 56% of the European Neogene genera are now globally and regionally extinct, respectively. The observed pulses of origination and extinction in the different regions coincide with warming and cooling events that occurred during the Neogene globally and regionally. Our study reveals complex diversity dynamics of Neogene chondrichthyans from Europe and their distinct biogeographic composition despite the multiple marine passages that connected the different marine regions during this time.

Type
Articles
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © The Author(s), 2023. Published by Cambridge University Press on behalf of The Paleontological Society

Introduction

The formation of the Mediterranean Sea was the result of intense geographic and oceanographic changes that took place over the last 150 Myr (Rögl Reference Rögl1999). During the Mesozoic, the Tethys Ocean separated the continents of Laurasia and Gondwana until the end of the Eocene, when continental drift and Alpine orogeny caused it to vanish (Rögl Reference Rögl1999; Berra and Angiolini Reference Berra, Angiolini, Marlow, Kendall and Yose2014). By the Oligocene, Europe's marine areas reorganized into three distinct biogeographic provinces: Paratethys, Mediterranean, and Atlantic boreal (Steininger and Wessely Reference Steininger and Wessely2000). The Paratethys formed a separate branch of the former Tethys Ocean in the north, and a proto–Mediterranean Sea formed in the south (Rögl Reference Rögl1998). As a consequence of these tectonic changes, Europe was largely covered by parts of the eastern Atlantic Ocean, the Mediterranean Sea, the Paratethys, and the North Sea, which together formed shallow epicontinental seas during the Neogene (23 to 2.6 Ma; Rögl Reference Rögl1998).

The connectivity of Europe's main marine regions resulted in major environmental and biotic changes during the Neogene. The Mediterranean and Paratethys Seas experienced a series of connections and disconnections through sea passages from the Miocene to Pliocene (Rögl Reference Rögl1999). During the Burdigalian (20.43–15.97 Ma), the Paratethys was connected to the Atlantic Ocean, the North Sea, and the Mediterranean Sea, enabling faunal exchanges (Kroh Reference Kroh2007; Reinecke et al. Reference Reinecke, Louwye, Havekost and Moths2011). At the end of the Burdigalian, the western part of the Paratethys fell dry, but the western and central Paratethys remained under marine conditions and connected to the Mediterranean Sea (Rögl Reference Rögl1999). From the Serravalian (13.82–11.63 Ma) onward, the Paratethys became gradually isolated from the Mediterranean Sea, and at the end of the Serravalian, the uplift of the Carpathian Mountains separated them (Harzhauser and Kowalke Reference Harzhauser and Kowalke2002; Harzhauser and Piller Reference Harzhauser and Piller2007). A gradual decrease of temperature and salinity caused a continuous development of endemic faunas in the Central Paratethys (Kroh Reference Kroh2007), and the closure of sea passages triggered the final isolation of the Paratethys and its transition from a marine to a freshwater environment. The Mediterranean Sea was affected by a gradual decrease of temperature and desiccation events during the Tortonian and Messinian (11.63–5.33 Ma; Rögl Reference Rögl1999), although marine organisms persisted throughout the three stages of the Messinian salinity crisis (e.g., Carnevale et al. Reference Carnevale, Gennari, Lozar, Natalicchio, Pellegrino and Dela Pierre2019). The temperature drop continued through the Pliocene, reaching its lowest levels during the beginning of the Pleistocene glaciations (Westerhold et al. Reference Westerhold, Marwan, Drury, Liebrand, Agnini, Anagnostou, Barnet, Bohaty, De Vleeschouwer, Florindo, Frederichs, Hodell, Holbourn, Kroon, Lauretano, Littler, Lourens, Lyle, Pälike, Röhl, Tian, Wilkens, Wilson and Zachos2020). The waxing and waning of large high-latitude ice sheets caused intense fluctuations of sea levels during the Quaternary (Lambeck and Chappell Reference Lambeck and Chappell2001), including in the Mediterranean Sea (Brunović et al. Reference Brunović, Miko, Hasan, Gapatheodorou, Ilijanić, Miserocchi, Correggiaria, Geraga, Channel and Adriatic Sea2020). As a result, marine faunas in Europe faced large and persistent environmental changes throughout the Neogene, which likely influenced their diversity and distribution.

Indeed, previous studies have regarded the intense climatic and oceanographic events of the European Neogene as the cause for the faunistic changes that occurred in that region during this time (Harzhauser et al. Reference Harzhauser, Mandic and Zuschin2003; Moissette et al. Reference Moissette, Dulai and Müller2006; Harzhauser and Piller Reference Harzhauser and Piller2007; Kroh Reference Kroh2007; Piller et al. Reference Piller, Harzhauser and Mandic2007; Borgh et al. Reference Borgh, Vasiliev, Stoica, Knežević, Matenco, Krijgsman, Rundic and Cloetingh2013). For instance, it has been proposed that the high diversity of marine taxa reached at the beginning of the middle Miocene was a result of temperature increase and favorable oceanographic conditions (Kroh Reference Kroh2007; Schwarzhans Reference Schwarzhans2010). However, these studies predominantly focused on invertebrates, and therefore, understanding of the response of marine communities to the large environmental changes that took place in the last 23 Myr in Europe remains elusive, as it lacks data from vertebrate communities. Despite the large amount of paleontological information that has been accumulated on Neogene chondrichthyans (chimaeras, sharks, rays, skates) from Europe in the last decades (e.g., Marsili Reference Marsili2008; Antunes and Balbino Reference Antunes and Balbino2010; Bor et al. Reference Bor, Reinecke and Verschueren2012; Cappetta Reference Cappetta and Schultze2012; Schultz Reference Schultz2013; Marramà et al. Reference Marramà, Schultz and Kriwet2019), studies of their diversification trajectories during this time of great environmental change are scarce and based on few localities (Kriwet and Klug Reference Kriwet and Klug2008; Reinecke et al. Reference Reinecke, Louwye, Havekost and Moths2011; Fuchs Reference Fuchs2015; Schwab Reference Schwab2015; Villafaña et al. Reference Villafaña, Marramà, Klug, Pollerspöck, Balsberger, Rivadeneira and Kriwet2020). Here, we synthesize the rich fossil record of Neogene chondrichthyans from Europe in order to reconstruct their diversification trajectories and paleobiogeographic dynamics across different geographic areas from the past to the present. Our results reveal the distinct biogeographic composition of chondrichthyan faunas during the Neogene of Europe and the potential links between diversification trajectories and global and regional climatic changes. As such, this study advances our understanding of the long-term, regional responses of marine communities to major environmental perturbations.

Materials and Methods

Data

We gathered chondrichthyan occurrences at the genus level from the Neogene (23–2.6 Ma) of Europe based on a comprehensive literature quest that consisted of searching for the terms “chondrichthyans”, “fossil”, “Neogene”, and “Europe” in Google Scholar (https://scholar.google.com). This resulted in a list of 122 journal articles, unpublished theses, conference abstracts, and books. This information was complemented with data downloaded from the Paleobiology Database (https://paleobiodb.org) and from museum online collection databases (Supplementary Table S1). Additionally, collections housed in the Natural History Museum of Vienna and the State Museum of Natural History of Stuttgart, Germany, were examined. In total, we collected 4368 occurrences (Fig. 1, Supplementary Table S1), which we assigned to four regions based on the paleogeographic reconstructions proposed by Rögl (Reference Rögl1999): Atlantic (n = 433), Mediterranean Sea (n = 750), North Sea (n = 563), and Paratethys (n = 2622). The regional stratigraphic stages (i.e., Paratethys Sea) were updated based on more recent studies (Grunert et al. Reference Grunert, Soliman, Ćorić, Scholger, Harzhauser and Piller2010; Heckeberg et al. Reference Heckeberg, Pippèrr, Läuchli, Heimann and Reichenbacher2010; Hohenegger et al. Reference Hohenegger, Ćorić and Wagreich2014; Kováč et al. Reference Kováč, Halásová, Hudáčková, Holcová, Hyžný, Jamrich and Ruman2018). Ambiguous records with unclear taxonomic names (i.e., nonvalid synonyms) or localities (i.e., assigned only to the country level) were excluded from the database. All taxonomic names were updated according to the most recent taxonomic reviews (Cappetta Reference Cappetta and Schultze2012; Pollerspöck and Straube Reference Pollerspöck and Straube2021).

Figure 1. Fossil occurrences of Neogene chondrichthyans from the Atlantic, Mediterranean, North Sea, and Paratethys regions.

It has been demonstrated that deep-water chondrichthyans are differently affected by abiotic factors than shallow-water taxa (e.g., Kriwet and Benton Reference Kriwet and Benton2004). However, we did not attempt to analyze deep- (<500 m) and shallow-water (>500 m) associations separately, because although several deep-water chondrichthyan faunas have been reported from the Miocene of the Paratethys (e.g., Underwood and Schlögl Reference Underwood and Schlögl2013) and the northern margin of the Tethys (e.g., Cigala-Fulgosi Reference Cigala-Fulgosi1996; Adnet Reference Adnet2006), up to now, no deep-water chondrichthyan assemblages from other regions analyzed here have been unambiguously identified. Additionally, dental morphologies of many modern deep-water chondrichthyans are still poorly known, so it is unclear whether taxa (except for most squaliforms) identified in deep-water settings can be associated with such habitat, or if they are representatives of shallow-water taxa but were mixed with deeper-water sediments taphonomically. Consequently, analyzing deep- and shallow-water taxa separately would likely introduce unintended artifacts.

Analyses

To reconstruct diversification trajectories (i.e., genus richness, origination, and extinction rates), we used the first and last appearance of each genus based on its occurrence, with occurrences distributed in 1 Ma time bins in each region (i.e., Mediterranean, North Sea, and Paratethys). We excluded the Atlantic region from the analysis due to its low number of occurrences (n = 433; 10% of the total diversity) and the lack of information on the stratigraphic ages of many localities. Genus richness was calculated per time bin using two approaches: (1) “boundary-crossers” (i.e., number of taxa that cross the boundary of the interval; Foote Reference Foote2000); and (2) shareholder quorum subsampling (SQS; i.e., fixed coverage of the frequency curve of genus occurrences; Alroy Reference Alroy, Alroy and Hunt2010). We used 1000 iterations and quorums of 0.4, 0.6, and 0.8. This last approach was implemented because it accounts for differences in sampling effort, unlike the other approaches, which are prone to sampling biases (Alroy Reference Alroy, Alroy and Hunt2010). However, a strong positive correlation between genus richness (estimated using the boundary-crosser method) and the subsampled genus richness (estimated using the SQS method) might suggest that sampling bias is relatively systematic in time (Supplementary Table S2). Origination and extinction rates were estimated as described in the per capita rates of Foote (Reference Foote1999). Origination rates (Eq. 1) are estimated as:

(1)$$p = -\ln [ {N_{bt}/( {N_{bt} + N_b} ) } ] $$

and extinction rates (Eq. 2) as:

(2)$$q = -\ln [ {N_{bt}/( {N_{bt} + N_t} ) } ] $$

where N bt is the number of taxa crossing both bottom and top interval boundaries, N b is the number of taxa crossing the bottom boundary, and N t is number of taxa crossing the top boundary. Singletons were excluded from the origination and extinction rate estimates as proposed by Foote (Reference Foote1999).

To reconstruct the overall extinction magnitude of European chondrichthyans from the Neogene to the Recent, we estimated the Lyellian percentages, that is, the proportion of Neogene genera still living in each region today. Information on present-day distributions of each genus was obtained from the Ocean Biogeographic Information System (2021) and FishBase (Froese and Pauly Reference Froese and Pauly2021). Paratethys records were compared against the present-day Mediterranean occurrences. Differences in the proportions of extinct genera across regions were compared using χ2 tests. We further made comparisons at higher taxonomic levels (i.e., order and family). We did so by only using the Mediterranean fauna as a comparative region because of its high diversity of chondrichthyans (86 species; Cariani et al. Reference Cariani, Messinetti, Ferrari, Arculeo, Bonello, Bonnici, Cannas, Carbonara, Cau, Charilaou, El Ouamari, Fiorentino, Follesa, Garofalo, Golani, Guarniero, Hanner, Hemida, Kada, Lo Brutto, Mancusi, Morey, Schembri, Serena, Sion, Stagoni, Tursi, Vrgoc, Steinke and Tinti2017; Ebert and Dando Reference Ebert and Dando2020) and because its fossil chondrichthyan record was intensively studied in the past (e.g., Cappetta Reference Cappetta and Schultze2012 and references therein; this study). All analyses were made in R (R Core Team 2021) using the Divdyn library (Kocsis et al. Reference Kocsis, Reddin, Alroy and Kiessling2019).

Results

Taxonomic Composition

The Neogene chondrichthyan fauna comprises 102 genera (Supplementary Table S1) representing four superorders, 12 orders, and 41 families (Table 1). At the superorder level, galeomorphs were the most dominant (i.e., >50%) in the four regions followed by squalomorphs (11–26%). Holocephalans were the least abundant (4%, 5 out of 102), with the Atlantic region displaying the highest number of genera (7%, 4 out of 54). Fossils of this group have not been recorded from the Mediterranean region so far. At the order level, carcharhiniforms and lamniforms were the most abundant groups in the four regions (24–31% and 19–29% of genera, respectively; Table 1, Supplementary Tables S3–S6). Within batoids, the order Myliobatiformes (11–17%) was the best-represented group. At the family level, carchahinids were the most abundant shark group in the Atlantic (11%, 6 out of 54), Mediterranean (10%, 7 out of 70), and Paratethys (10%, 7 out of 69) regions, whereas lamnids (11%, 7 out of 62) were the most common shark family in the North Sea (Table 1). Among batoids, the family Myliobatidae was the most abundant in the Atlantic (7%, 4 out of 54), Mediterranean (6%, 4 out of 70), and Paratethys (6%, 4 of 69) regions, whereas rajids (6%, 4 out of 62) were the most abundant batoids in the North Sea region (Table 1).

Table 1. Faunal composition of Neogene European chondrichthyans at superorder, order, and family levels.

Diversification Trajectories

The diversity analyses showed marked differences between regions (Fig. 2, Table 2). In the Paratethys, the number of occurrences and genus richness increased and reached maximum values around the early Miocene, but subsequently decreased toward the late Miocene (Fig. 2A,B). In the Mediterranean (Fig. 2E,F) and North Seas (Fig. 2I,J), the number of occurrences and genera increased during the early Miocene, reaching maximum values around the middle Miocene, and then decreasing toward the early Pliocene. The only significant correlation between regions in terms of genus richness was found to be between the Mediterranean and North Sea regions (r = 0.80, p < 0.05; Table 2).

Figure 2. Diversification trajectories of European chondrichthyans. A, E, I, Number of occurrences; B, F, J, standing generic richness; C, G, K, origination rates; and D, H, L, extinction rates. The confidence intervals around the lines were excluded for clarity. Abbreviations: early Miocene (eM), middle Miocene (mM), late Miocene (lM) and Pliocene (Pl). Quorums used to estimate genus richness based on SQS are represented with black lines

Table 2. Spearman moment-product correlation for diversity trajectories between paired regions.a

a Significant values (p < 0.05) are in bold.

In terms of origination rates (Fig. 2C,G,K), the Paratethys and the North Sea showed a main pulse during the early Miocene (Fig. 2C,K), whereas in the Mediterranean, two main pulses were observed in the early and middle Miocene (Fig. 2G). There was no significant correlation in origination rates between regions (Table 2). Extinction rates presented substantial pulses during the entire time studied in all regions (Fig. 2D,H,L), with the highest peaks taking place in the Paratethys during the late Miocene (Fig. 2D). In the Mediterranean, there was one main extinction pulse during the Pliocene (Fig. 2H). In the North Sea, there were two main extinction pulses in the late Miocene and Pliocene (Fig. 2L). Extinction rates were significantly correlated between the Mediterranean and North Sea region (r = 0.49, p < 0.05; Table 2) and the Mediterranean and the Paratethys (r = 0.59, p < 0.05; Table 2).

The trajectories observed seem to be independent of the method used (Supplementary Table S2). For all regions, the boundary-crosser method was significantly correlated with at least one of the quorums used to estimate the taxonomic richness based on SQS. For instance, in the Paratethys region, the taxonomic richness based on SQS (quota = 0.8) significantly correlates with the boundary-crosser method results (r = 0.97, p < 0.05) (Supplementary Table S2).

Biogeography

The comparison between the fossil and current distributions (Fig. 3, Supplementary Tables S3–S6) at the genus level showed that 27% (28 out of 102; see “combined” bar in Fig. 3) of genera from the Neogene of Europe are now globally extinct (e.g., Megascyliorhinus, Otodus, and Striatolamia). The proportion of globally extinct genera was only different between the Paratethys (25%, 18 out of 71; Fig. 3) and Atlantic faunas (15%, 8 out of 54; χ2 = 3.84, df = 1, p = 0.049; Table 3). The proportion of regionally extinct genera was much higher (40–56%; Fig. 3), but there were no differences between regions (p > 0.05 in all cases; Table 3). The biogeographic distribution of European chondrichthyans among the different regions from the Neogene to the Recent is shown in Supplementary Tables S3–S6.

Figure 3. Percentages of extinct/living chondrichthyan genera at each region, at global scale (i.e., the taxon is no longer present in the global ocean) and regional scale (i.e., the taxon was extirpated from the region but is still living elsewhere).

Table 3. Differences in the proportions of extinct genera at global and regional scales between regions using χ2 test.a

a Significant values (p < 0.05) are in bold.

At the order level, 12 out of 14 extant orders (86%) of chondrichthyans are present in the Neogene of Europe. Using the current diversity of the Mediterranean Sea for comparison, 12 extant orders are shared with the Neogene. Two of the shared orders are absent from the Mediterranean Sea today. As such, orectilobiforms and pristiophoriforms are found in the Neogene, but are absent today in the Mediterranean Sea (Supplementary Table S7). Carcharhiniforms, lamniforms, myliobatiforms, and squaliforms display the highest number of genera both in the Neogene and today (between 11% and 17% of the Neogene generic diversity and between 16% and 18% of today's diversity in the Mediterranean Sea; Supplementary Table S7). Although chimaeriforms and torpediniforms have not been recovered from the Neogene, they are present today; however, these are the least genus-rich chondrichthyan orders in the Mediterranean Sea today (1% and 2% of the total diversity, respectively; Supplementary Table S7).

At the family level, 40 out of 66 extant families of chondrichthyans were present in the Neogene of Europe (61%). Comparisons based on the Mediterranean Sea show that there are 36 chondrichthyan families shared between the Neogene and today. From these, 27 families (75%) are present today, with 24 (67%) present in both the Neogene and today (Supplementary Table S7). Carcharhinidae was the most genus-rich family in the Neogene, comprising 10% of the diversity (7 out of 70 fossil genera), whereas Rajidae and Dasyatidae are the most diverse families in the Mediterranean Sea today, comprising 11% and 9% of the current diversity (5 and 4 out of 45, respectively; Supplementary Table S7). The family Lamnidae represents one of the second most genus-rich groups in the Neogene (6%, 4 out of 70) and today (7%, 3 out of 45). Fossil specimens of aetobatids, chlamydoselachids, ginglymostomatids, hemigaleids, mitsukurinids, otondontids, plesiobatids, pristiophorids, and rhinids have been recovered from the Neogene, but are absent today in the Mediterranean Sea. Chimaerids, oxynotids, and torpedinids represent the only families not recovered from Neogene localities, but are present today in the Mediterranean region (Supplementary Table S7).

Discussion

Diversification Trajectories

Our analyses revealed marked differences in Neogene chondrichthyan diversification trajectories between the studied regions (Fig. 2; Supplementary Fig. S1), likely reflecting the impact of paleoenvironmental conditions that occurred at different spatial scales across Europe and the globe (Fig. 4). For instance, the peak in diversity observed in the Paratethys Sea in the early Miocene (~17.5 Ma), which is evidenced by the high genus richness and origination rates of the time, was coeval with warm temperate conditions during the Eggenburgian (~20.4–18.3 Ma; Nebelsick Reference Nebelsick1992), which may have also promoted elevated invertebrate origination rates (Kroh Reference Kroh2007). Similarly, in the North Sea, the highest number of occurrences and genera took place in the middle Miocene, while origination rates peaked in the early Miocene (19.5 Ma; Fig. 2). Around this time, specifically between ~19 and 14.5 Ma, the North Sea experienced a warming event due to the incursion of warm-temperate Atlantic waters through a southwest-oriented sea passage (Gürs Reference Gürs2001) that has been associated with increased species richness in vertebrates and invertebrates (e.g., Gürs and Janssen Reference Gürs and Janssen2002; Kowalewski et al. Reference Kowalewski, Gürs, Nebelsick, Oschmann, Piller and Hoffmeister2002; Moths et al. Reference Moths, Albrecht and Stein2010; Schwarzhans Reference Schwarzhans2010; Reinecke et al. Reference Reinecke, Louwye, Havekost and Moths2011). In the case of the Mediterranean Sea, there is a steep increase in generic diversity in the early and middle Miocene and two origination peaks at 18.5 and 14.5 Ma. These peaks coincide with the two intervals when the Mediterranean Sea and the Indo-Pacific were connected (23–18 Ma and ~16–15 Ma; Rögl Reference Rögl1999), which resulted in a temperature increase in the region (Harzhauser et al. Reference Harzhauser, Kroh, Mandic, Piller, Göhlich, Reuter and Berning2007; Vertino et al. Reference Vertino, Stolarski, Bosellini, Taviani, Goffredo and Dubinsky2014). Taken together, these results suggest that increases in Neogene chondrichthyan diversity coincided with regional warming events during the early and middle Miocene. In current marine systems, species richness tends to be higher in warm areas (e.g., the tropics) than in colder environments (Hillebrand Reference Hillebrand2004; Kinlock et al. Reference Kinlock, Prowant, Herstoff, Foley, Akin-Fajiye, Bender, Umarani, Ryu, Şen and Gurevitch2018). One of the possible explanations of this temperature dependence for diversity is the kinetic energy or temperature hypothesis, which postulates that high temperatures increase metabolic rates, promoting higher rates of speciation, ultimately leading to greater diversity (Tittensor et al. Reference Tittensor, Mora, Jetz, Lotze, Ricard, Berghe and Worm2010).

Figure 4. Climatic, oceanographic and biotic events through the Neogene and their relation with the origination (blue arrows) and extinction (red arrows) of Mediterranean (M), North Sea (NS), and Paratethys (P) chondrichthyans. 1: Westerhold et al. (Reference Westerhold, Marwan, Drury, Liebrand, Agnini, Anagnostou, Barnet, Bohaty, De Vleeschouwer, Florindo, Frederichs, Hodell, Holbourn, Kroon, Lauretano, Littler, Lourens, Lyle, Pälike, Röhl, Tian, Wilkens, Wilson and Zachos2020); 2: Abreu and Haddad (Reference Abreu, Haddad, Graciansky, Hardenbol, Jacquin and Vail1998); 3: Nebelsick (Reference Nebelsick1992); 4: Kroh (Reference Kroh2007); 5: Gürs and Janssen (Reference Gürs and Janssen2002); 6: Rögl (Reference Rögl1999); 7: Harzhauser et al. (Reference Harzhauser, Kroh, Mandic, Piller, Göhlich, Reuter and Berning2007); 8: Reinecke et al. (Reference Reinecke, Louwye, Havekost and Moths2011); 9: Harzhauser et al. (Reference Harzhauser, Mandic and Zuschin2003); 10: Harzhauser and Piller (Reference Harzhauser and Piller2007), 11: Piller et al. (Reference Piller, Harzhauser and Mandic2007); 12: Bannikov et al. (Reference Bannikov, Schwarzhans and Carnevale2018); 13: Hohenegger et al. (Reference Hohenegger, Ćorić and Wagreich2014); 14: Krijgsman et al. (Reference Krijgsman, Garcés, Agustı, Raffi, Taberner and Zachariasse2000); 15: Krijgsman et al. (Reference Krijgsman, Stoica, Vasiliev and Popov2010); 16: Todd et al. (Reference Todd, Jackson, Johnson, Fortunato, Heitz, Alvarez and Jung2002); 17: Rivadeneira and Marquet (Reference Rivadeneira and Marquet2007); 18: Villafaña and Rivadeneira (Reference Villafaña and Rivadeneira2014). CP, Central Paratethys; NS, North Sea.

The extinction peaks in the different regions also coincide with regional and global climatic and oceanographic events (Fig. 4). The highest extinction peak in the Paratethys at 11.5 Ma coincides with the isolation of the Central Paratethys Sea from all surrounding marine environments during the Sarmatian/Pannonian (11.6 Ma; Rögl Reference Rögl1999). This isolation caused a gradual change from marine to freshwater conditions, triggering the regional Sarmatian–Pannonian extinction event that resulted in the disappearance of more than 90% of gastropods (Harzhauser and Piller Reference Harzhauser and Piller2007; Borgh et al. Reference Borgh, Vasiliev, Stoica, Knežević, Matenco, Krijgsman, Rundic and Cloetingh2013). The second-highest extinction peak in the Paratethys (14.5 Ma) coincides with an intense cooling event (14.25 Ma; Abreu and Haddad Reference Abreu, Haddad, Graciansky, Hardenbol, Jacquin and Vail1998) that could have been responsible for the extinction of mollusks, bony fishes, and foraminifera in the middle Badenian (Harzhauser and Piller Reference Harzhauser and Piller2007; Piller et al. Reference Piller, Harzhauser and Mandic2007; Borgh et al. Reference Borgh, Vasiliev, Stoica, Knežević, Matenco, Krijgsman, Rundic and Cloetingh2013; Hohenegger et al. Reference Hohenegger, Ćorić and Wagreich2014; Bannikov et al. Reference Bannikov, Schwarzhans and Carnevale2018). In the North Sea region, the first main extinction peak took place at 8.5 Ma, after the connection to the Atlantic Sea was closed (14.5 Ma; Gürs Reference Gürs2001) and a major ice sheet on Antarctica was reestablished (14 Ma; Westerhold et al. Reference Westerhold, Marwan, Drury, Liebrand, Agnini, Anagnostou, Barnet, Bohaty, De Vleeschouwer, Florindo, Frederichs, Hodell, Holbourn, Kroon, Lauretano, Littler, Lourens, Lyle, Pälike, Röhl, Tian, Wilkens, Wilson and Zachos2020). Around 4.5 Ma, there was another extinction peak both in the North Sea and Mediterranean regions, coinciding with the onset of the Northern Hemisphere glaciation (Westerhold et al. Reference Westerhold, Marwan, Drury, Liebrand, Agnini, Anagnostou, Barnet, Bohaty, De Vleeschouwer, Florindo, Frederichs, Hodell, Holbourn, Kroon, Lauretano, Littler, Lourens, Lyle, Pälike, Röhl, Tian, Wilkens, Wilson and Zachos2020). Our results indicate high extinction rates in European chondrichthyans mostly from the middle Miocene onward that coincided with a global decline in sea temperature linked to the polar glaciation (Westerhold et al. Reference Westerhold, Marwan, Drury, Liebrand, Agnini, Anagnostou, Barnet, Bohaty, De Vleeschouwer, Florindo, Frederichs, Hodell, Holbourn, Kroon, Lauretano, Littler, Lourens, Lyle, Pälike, Röhl, Tian, Wilkens, Wilson and Zachos2020). Indeed, cooling events in Cenozoic have also been linked to the global decline in diversity of lamniform sharks (Condamine et al. Reference Condamine, Romieu and Guinot2019). As such, the chondrichthyan extinctions in the Neogene of Europe could be related to cooling events.

Overall, our results indicate that chondrichthyan diversification trajectories during the Neogene of Europe coincided with regional and global warming and cooling events (Fig. 4). In line with this observation, sea temperature has been shown to be a major driver of chondrichthyan species richness today (Tittensor et al. Reference Tittensor, Mora, Jetz, Lotze, Ricard, Berghe and Worm2010; Guisande et al. Reference Guisande, Patti, Vaamonde, Manjarrés-Hernández, Pelayo-Villamil, García-Roselló, González-Dacosta, Heine and Granado-Lorencio2013). Nevertheless, to be able to postulate a causal link between ocean temperature and peaks in chondrichthyan origination and extinction rates, future research would need to analyze the effect of sea temperature in origination and extinction rates and their interaction with the ecological and life-history traits of each taxon.

Chondrichthyan Faunas from Europe: Neogene versus Recent

Twenty-seven percent of the genera that inhabited the marine environments of Europe during the Neogene are now globally extinct (Fig. 3), with as much as 56% being regionally extinct (i.e., in the North Sea; Supplementary Table S5). The level of extinction found at the generic level in the Neogene of Europe is much higher than in other areas. For instance, only 3% of the genera from the Neogene of the Pacific of South America (Villafaña and Rivadeneira Reference Villafaña and Rivadeneira2014, Reference Villafaña and Rivadeneira2018) are now globally extinct, and 34% are regionally extinct. Similarly, only 13% of the genera from tropical America (Carrillo-Briceño et al. Reference Carrillo-Briceño, Carrillo, Aguilera and Sanchez-Villagra2018) are now globally extinct. We attribute these differences in extinction intensities to the large destruction of marine habitats in Europe associated with the drying of the Paratethys and the fluctuation of sea levels in the Mediterranean Sea (Rögl Reference Rögl1999; Brunović et al. Reference Brunović, Miko, Hasan, Gapatheodorou, Ilijanić, Miserocchi, Correggiaria, Geraga, Channel and Adriatic Sea2020). In contrast, despite the significant oceanic changes associated with the rise of the Isthmus of Panama and the closure of the Central American Seaway (Klaus et al. Reference Klaus, Lutz, Mcneill, Budd, Johnson and Ishman2011; Montes et al. Reference Montes, Cardona, Jaramillo, Pardo, Silva, Valencia, Ayala, Pérez-Angel, Rodriguez-Parra, Ramirez and Niño2015; O'Dea et al. Reference O'Dea, Lessios, Coates, Eytan, Restrepo-Moreno, Cione, Collins, Queiroz, Farris, Norris, Stallard, Woodburne, Aguilera, Aubry, Berggren, Budd, Cozzuol, Coppard, Duque-Caro, Finnegan, Gasparini, Grossman, Johnson, Keiwin, Knowlton, Leigh, Leonard-Pingel, Marko, Pyenson, Rachello-Dolmen, Soibelzon, Soibelzon, Todd, Vermeij and Jackson2016), the loss of epicontinental seas is not recorded in the Pacific of South America or in tropical America (Miller et al. Reference Miller, Kominz, Browning, Wright, Mountain, Katz, Sugarman, Cramer, Christie-Blick and Pekar2005; Le Roux et al. Reference Roux, Achurra, Henríquez, Carreño, Rivera, Suárez, Ishman, Pyenson and Gutstein2016). Therefore, the intense oceanographic events ultimately resulting in the vanishing of the Paratethys may have had a more significant effect on chondrichthyan faunas in Europe relative to those from the Americas. Indeed, habitat loss is linked to chondrichthyan extinctions in both the fossil record and modern seas (Cione et al. Reference Cione, Mennucci, Santalucita and Hospitaleche2007; Pimiento et al. Reference Pimiento, Griffin, Clements, Silvestro, Varela, Uhenm and Jaramillo2017; Dulvy et al. Reference Dulvy, Pacoureau, Rigby, Pollom, Jabado, Ebert, Finucci, Pollock, Cheok, Derrick, Herman, Sherman, VanderWright, Lawson, Walls, Carlson, Charvet, Bineesh, Fernando, Ralph, Matsushiba, Hilton-Taylor, Fordham and Simpfendorfer2021). Nevertheless, to further unveil the mechanisms of chondrichthyan extinctions in Europe, and specifically the role of habitat loss, future studies should analyze how sea-level changes affect species diversification in the region.

As expected, the observed changes in the generic composition vary with taxonomic level. As such, the proportion of genera across orders is similar between the past to the present, with 12 extant orders being found both in the Neogene of Europe and today in the Mediterranean Sea. However, Torpediniformes and Chimaeriformes are absent from the Neogene and Orectolobiformes and Pristiophoriformes are found in the Neogene but are absent today in the Mediterranean Sea (Supplementary Table S7). Indeed, it has been shown that the chondrichthyan fossil record is largely conserved at the order level relative to today (Pimiento and Benton Reference Pimiento and Benton2020; Paillard et al. Reference Paillard, Shimada and Pimiento2021). At the family level there are larger differences, with nine families (Supplementary Table S7) present in the Neogene of the Mediterranean Sea but absent today, and three families (Chimaeridae, Oxynotidae and Torpedinidae) found in the Mediterranean Sea today, but absent from the Neogene. Our results therefore suggest that Neogene chondrichthyan faunas from Europe suffered a significant taxonomic loss relative to other regions, as well as a possible rearrangement of taxonomic composition at the order and family levels.

Conclusions and Perspectives

Our analysis provides first steps toward synthesizing and understanding the macroevolutionary diversification trajectories and paleobiogeographic changes of chondrichthyans during the Neogene in Europe. Increases in number of genera and peaks of origination seem to be related to warming events, whereas the highest extinction peaks appear to be related to cooling events. The biogeographic comparison between past and present shows that global and regional extinctions of chondrichthyans in Europe were much higher than in other regions. The present study nevertheless should be considered an initial step toward understanding Neogene diversity patterns of Europe; additional paleontological studies are needed to provide additional data for further refinements of the present analyses, including (1) more precise dating of fossil inventories; (2) more specific information related to the collection method used; (3) inclusion of ecological and life-history traits of each taxon in order to understand the biogeographic dynamics through the time; (4) establishment of the correct ecology for extinct taxa based on detailed comparisons with living taxa to better distinguish between deep- (<500) and shallow-water (>500 m) chondrichthyan taxa to determine the influence of abiotic factors on their diversity fluctuations and extinction risk; and (5) more robust tests of the role of environmental drivers based on new paleoceanographic reconstructions.

Acknowledgments

This research was supported by Becas de Doctorado en el Extranjero, Becas Chile (grant 72160271) and a SYNTHESYS grant (DE-TAF). We deeply appreciate the comments and help provided by J. Pollerspöck (Stephansposching, Germany), S. Stumpf and I. Feichtinger (Vienna, Austria). For their kindness in providing access to paleontological collections, we would like to thank E. Maxwell (Staatliches Museum fur Naturkunde Stuttgart) and U. Göhlich (Naturhistorisches Museum Wien). This project was partially funded by ANID/FONDECYT no. 1200843, ANID/CENTROS REGIONALES R20F0008 (CLAP), and Millennium Science Initiative Nucleus UPWELL (NCN19–153) (to M.M.R.). C.P. is funded by a PRIMA grant from the Swiss National Science Foundation (no. 185798). This research was funded in part by the Austrian Science Fund (FWF) [P 33820] to J.K. CC BY or equivalent license is applied to the Author Accepted Manuscript (AAM) arising from this submission, in accordance with the grant's open access conditions.

Declaration of Competing Interests

The authors have declared that no competing interests exist.

Data Availability Statement

Data available from the Dryad Digital Repository: https://doi.org/10.5061/dryad.34tmpg4pd.

References

Literature Cited

Abreu, V. S., and Haddad, G. A.. 1998. Glacioeustatic fluctuations: the mechanism linking stable isotope events and sequence stratigraphy from the Early Oligocene to Middle Miocene. In Graciansky, C. P., Hardenbol, J., Jacquin, T., and Vail, P. R., eds. Mesozoic and Cenozoic sequence stratigraphy of European basins. SEPM Special Publication 60:245260.Google Scholar
Adnet, S. 2006. Nouvelles faunes de sélaciens (Elasmobranchii, Neoselachii) de l’Éocène des Landes (Sud-Ouest, France). Implication dans les connaissances des communautés d'eaux profondes. Palaeo Ichthyologica 10:1128.Google Scholar
Alroy, J. 2010. Fair sampling of taxonomic richness and unbiased estimation of origination and extinction rates. In Alroy, J. and Hunt, G., eds. Quantitative methods in paleobiology. Paleontological Society Papers 16:5580.Google Scholar
Antunes, M. T., and Balbino, A. C.. 2010. Uppermost Miocene lamniform selachians (Pisces) from the Alvalade basin (Portugal). Ciencias da Terra 15:141154.Google Scholar
Bannikov, A. F., Schwarzhans, W., and Carnevale, G.. 2018. Neogene Paratethyan croakers (Teleostei, Sciaenidae). Rivista Italiana di Paleontologia e Stratigrafia 124:535571.Google Scholar
Berra, F., and Angiolini, L.. 2014. The evolution of the Tethys region throughout the Phanerozoic: a brief tectonic reconstruction. In Marlow, L., Kendall, C., and Yose, L., eds. Petroleum systems of the Tethyan region. American Association of Petroleum Geologists Memoir 106:127.Google Scholar
Bor, T., Reinecke, T., and Verschueren, S.. 2012. Miocene Chondrichthyes from Winterswijk-Miste, the Netherlands. Palaeontos 21:1136.Google Scholar
Borgh, M., Vasiliev, I., Stoica, M., Knežević, S., Matenco, L., Krijgsman, W., Rundic, L., and Cloetingh, S.. 2013. The isolation of the Pannonian basin (Central Paratethys): new constraints from magnetostratigraphy and biostratigraphy. Global and Planetary Change 103:99118.CrossRefGoogle Scholar
Brunović, D., Miko, S., Hasan, O., Gapatheodorou, G., Ilijanić, N., Miserocchi, S., Correggiaria, A., and Geraga, M.. 2020. Late Pleistocene and Holocene paleoenvironmental reconstruction of a drowned karst isolation basin (Channel, Lošinj, Adriatic Sea, NE). Palaeogeography, Palaeoclimatology, Palaeoecology 544:109587.Google Scholar
Cappetta, H. 2012. Chondrichthyes: Mesozoic and Cenozoic Elasmobranchii: teeth. Pp. 1512 in Schultze, H. P., ed. Handbook of paleoichthyology, Vol. 3E. Verlag F. Pfeil, Munich.Google Scholar
Cariani, A., Messinetti, S., Ferrari, A., Arculeo, M., Bonello, J. J., Bonnici, L., Cannas, R., Carbonara, P., Cau, A., Charilaou, C., El Ouamari, N., Fiorentino, F., Follesa, M. C., Garofalo, G., Golani, D., Guarniero, I., Hanner, R., Hemida, F., Kada, O., Lo Brutto, S., Mancusi, C., Morey, G., Schembri, P. J., Serena, F., Sion, L., Stagoni, M., Tursi, A., Vrgoc, N., Steinke, D., and Tinti, F.. 2017. Improving the conservation of Mediterranean chondrichthyans: the ELASMOMED DNA barcode reference library. PLoS ONE 12: e0170244.CrossRefGoogle ScholarPubMed
Carnevale, G., Gennari, R., Lozar, F., Natalicchio, M., Pellegrino, L., and Dela Pierre, F.. 2019. Living in a deep desiccated Mediterranean Sea: an overview of the Italian fossil record of the Messinian salinity crisis. Bollettino della Società Paleontologica Italiana 58:109140Google Scholar
Carrillo-Briceño, J. D., Carrillo, J. D., Aguilera, O. A., and Sanchez-Villagra, M. R.. 2018. Shark and ray diversity in the Tropical America (Neotropics)—an examination of environmental and historical factors affecting diversity. PeerJ 6:e5313.CrossRefGoogle ScholarPubMed
Cigala-Fulgosi, F. 1996. Rare oceanic deep water squaloid sharks from the Lower Pliocene of the Northern Apennines (Parma province, Italy). Bollettino della Società Paleontologica Italiana 34:301322.Google Scholar
Cione, A. L., Mennucci, J. A., Santalucita, F., and Hospitaleche, C. A.. 2007. Local extinction of sharks of genus Carcharias Rafinesque, 1810 (Elasmobranchii, Odontaspididae) in the eastern Pacific Ocean. Andean Geology 34:139145.CrossRefGoogle Scholar
Condamine, F.L., Romieu, J., and Guinot, G.. 2019. Climate cooling and clade competition likely drove the decline of lamniform sharks. Proceedings of the National Academy of Sciences USA 116:2058420590.Google ScholarPubMed
Dulvy, N. K., Pacoureau, N., Rigby, C. L., Pollom, R. A., Jabado, R. W., Ebert, D. A., Finucci, B., Pollock, C. M., Cheok, J., Derrick, D. H., Herman, K. B., Sherman, C. S., VanderWright, W. J., Lawson, J. M., Walls, R. H. L., Carlson, J. K., Charvet, P., Bineesh, K. K., Fernando, D., Ralph, G. M., Matsushiba, J. H., Hilton-Taylor, C., Fordham, S. V., and Simpfendorfer, C. A. 2021. Overfishing drives over one-third of all sharks and rays toward a global extinction crisis. Current Biology 31:47734787.CrossRefGoogle Scholar
Ebert, D. A., and Dando, M.. 2020. Field guide to sharks, rays and chimaeras of Europe and the Mediterranean. Princeton University Press, Princeton, N.J.Google Scholar
Foote, M. 1999. Morphological diversity in the evolutionary radiation of Paleozoic and post-Paleozoic crinoids. Paleobiology 25:1115.CrossRefGoogle Scholar
Foote, M. 2000. Origination and extinction components of taxonomic diversity: general problems. Paleobiology 26:74102.CrossRefGoogle Scholar
Froese, R., and Pauly, D.. 2021. eds. FishBase. www.fishbase.org, November 2021 version.Google Scholar
Fuchs, I. 2015. Palaeoecology and diversity of early Miocene Elasmobranchii of the Paratethys of Austria. Bachelor's thesis, University of Vienna, Vienna.Google Scholar
Grunert, P., Soliman, A., Ćorić, S., Scholger, R., Harzhauser, M., and Piller, W. E.. 2010. Stratigraphic re-evaluation of the stratotype for the regional Ottnangian stage (Central Paratethys, middle Burdigalian). Newsletters on Stratigraphy 44:116.CrossRefGoogle Scholar
Guisande, C., Patti, B., Vaamonde, A., Manjarrés-Hernández, A., Pelayo-Villamil, P., García-Roselló, E., González-Dacosta, J., Heine, J., and Granado-Lorencio, C.. 2013. Factors affecting species richness of marine elasmobranchs. Biodiversity and Conservation 22:17031714.CrossRefGoogle Scholar
Gürs, K. 2001. Neues zur Fauna und Entwicklung des Nordseebeckens im Mittel- und Obermioän. Meyniana 53:5174.Google Scholar
Gürs, K., and Janssen, A. W.. 2002. Revised pteropod biostratigraphy for the Miocene of the North Sea Basin. Pp. 117131 in Proceedings of the Eighth Biannual Meeting of the RCNNS/RCNPS. Landesamt für Natur und Umwelt des Landes Schleswig-Holstein, Flintbek, Germany.Google Scholar
Harzhauser, M., and Kowalke, T.. 2002. Sarmatian (Late Middle Miocene) gastropod assemblages of the central Paratethys. Facies 46:5782.CrossRefGoogle Scholar
Harzhauser, M., and Piller, W. E.. 2007. Benchmark data of a changing sea—palaeogeography, palaeobiogeography and events in the Central Paratethys during the Miocene. Palaeogeography, Palaeoclimatology, Palaeoecology 253:831.CrossRefGoogle Scholar
Harzhauser, M., Mandic, O., and Zuschin, M.. 2003. Changes in Paratethyan marine molluscs at the Early/Middle Miocene transition-diversity, palaeogeography and palaeoclimate. Acta Geologica Polonica 53:323339.Google Scholar
Harzhauser, M., Kroh, A., Mandic, O., Piller, W. E., Göhlich, U., Reuter, M., and Berning, B.. 2007. Biogeographic responses to geodynamics: a key study all around the Oligo–Miocene Tethyan Seaway. Zoologischer Anzeiger—A Journal of Comparative Zoology 246:241256.CrossRefGoogle Scholar
Heckeberg, N., Pippèrr, M., Läuchli, B., Heimann, F. U., and Reichenbacher, B.. 2010. The Upper Marine Molasse (Burdigalian, Ottnangian) in Southwest Germany—facies interpretation and a new lithostratigraphic terminology. Zeitschrift der Deutschen Gesellschaft für Geowissenschaften 161:285302.CrossRefGoogle Scholar
Hillebrand, H. 2004. On the generality of the latitudinal diversity gradient. American Naturalist 163:192211.CrossRefGoogle ScholarPubMed
Hohenegger, J., Ćorić, S., and Wagreich, M.. 2014. Timing of the middle Miocene Badenian stage of the central Paratethys. Geologica Carpathica 65:5566.CrossRefGoogle Scholar
Kinlock, N. L., Prowant, L., Herstoff, E. M., Foley, C. M., Akin-Fajiye, M., Bender, N., Umarani, M., Ryu, H. Y., Şen, B., and Gurevitch, J.. 2018. Explaining global variation in the latitudinal diversity gradient: meta-analysis confirms known patterns and uncovers new ones. Global Ecology and Biogeography 27:125141.CrossRefGoogle Scholar
Klaus, J. S., Lutz, B. P., Mcneill, D. F., Budd, A. F., Johnson, K. G., and Ishman, S. E.. 2011. Rise and fall of Pliocene free-living corals in the Caribbean. Geology 39:375378.CrossRefGoogle Scholar
Kocsis, A. T., Reddin, C. J., Alroy, J., Kiessling, W.. 2019. The R package divDyn for quantifying diversity dynamics using fossil sampling data. Methods in Ecology and Evolution 10:735743.CrossRefGoogle Scholar
Kowalewski, M., Gürs, K., Nebelsick, J. H., Oschmann, W., Piller, W. E., and Hoffmeister, A. P.. 2002. Multivariate hierarchical analyses of Miocene mollusk assemblages of Europe: paleogeographic, paleoecological, and biostratigraphic implications. Geological Society of America Bulletin 114:239256.2.0.CO;2>CrossRefGoogle Scholar
Kriwet, J., and Benton, M. J.. 2004. Neoselachian (Chondrichthyes, Elasmobranchii) diversity across the Cretaceous–Tertiary boundary. Palaeogeography, Palaeoclimatology, Palaeoecology 214:181194.CrossRefGoogle Scholar
Kriwet, J., and Klug, S.. 2008. Diversity and biogeography patterns of Late Jurassic neoselachians (Chondrichthyes: Elasmobranchii). Geological Society of London Special Publication 295:5570.CrossRefGoogle Scholar
Kováč, M., Halásová, E. V. A., Hudáčková, N., Holcová, K., Hyžný, M., Jamrich, M., and Ruman, A.. 2018. Towards better correlation of the Central Paratethys regional time scale with the standard geological time scale of the Miocene Epoch. Geologica Carpathica 69:283300.CrossRefGoogle Scholar
Krijgsman, W., Garcés, M., Agustı, J., Raffi, I., Taberner, C., and Zachariasse, W. J.. 2000. The “Tortonian salinity crisis” of the eastern Betics (Spain). Earth and Planetary Science Letters 181:497511.CrossRefGoogle Scholar
Krijgsman, W., Stoica, M., Vasiliev, I., and Popov, V. V.. 2010. Rise and fall of the Paratethys Sea during the Messinian Salinity Crisis. Earth and Planetary Science Letters 290:183191.CrossRefGoogle Scholar
Kroh, A. 2007. Climate changes in the Early to Middle Miocene of the Central Paratethys and the origin of its echinoderm fauna. Palaeogeography, Palaeoclimatology, Palaeoecology 253:169207.CrossRefGoogle Scholar
Lambeck, K., and Chappell, J.. 2001. Sea level change through the last glacial cycle. Science 292:679686.CrossRefGoogle ScholarPubMed
Roux, Le, J. P., Achurra, L., Henríquez, A., Carreño, C., Rivera, H., Suárez, M. E., Ishman, S. E., Pyenson, N. D., and Gutstein, C. S.. 2016. Oroclinal bending of the Juan Fernández Ridge suggested by geohistory analysis of the Bahía Inglesa Formation, north-central Chile. Sedimentary Geology 333:3249.CrossRefGoogle Scholar
Marramà, G., Schultz, O., and Kriwet, J.. 2019. A new Miocene skate from the central Paratethys (Upper Austria): the first unambiguous skeletal record for the Rajiformes (Chondrichthyes: Batomorphii). Journal of Systematic Palaeontology 17:937960.CrossRefGoogle ScholarPubMed
Marsili, S. 2008. Systematic, paleoecologic and paleobiogeographic analysis of the Plio–Pleistocene Mediterranean elasmobranch fauna. Atti Della Società Toscana Di Scienze Naturali Memorie, serie A 113:8188.Google Scholar
Miller, K. G., Kominz, M. A., Browning, J. V., Wright, J. D., Mountain, G. S., Katz, M. E., Sugarman, P. J., Cramer, B. S., Christie-Blick, N., and Pekar, S. F.. 2005. The Phanerozoic record of global sea-level change. Science 310:12931298.CrossRefGoogle ScholarPubMed
Moissette, P., Dulai, A., and Müller, P.. 2006. Bryozoan faunas in the Middle Miocene of Hungary: biodiversity and biogeography. Palaeogeography, Palaeoclimatology, Palaeoecology 233:300314.CrossRefGoogle Scholar
Montes, C., Cardona, A., Jaramillo, C., Pardo, A., Silva, J. C., Valencia, V., Ayala, C., Pérez-Angel, L. C., Rodriguez-Parra, L. A., Ramirez, V., and Niño, H.. 2015. Middle Miocene closure of the Central American seaway. Science 348:226229.CrossRefGoogle ScholarPubMed
Moths, H., Albrecht, F., and Stein, G.. 2010. Die Molluskenfauna (Hemmoorium, Untermiozän) aus der Kiesgrube Krinke bei Werder (Nordwest-Niedersachsen). Palaeofocus 3:1155.Google Scholar
Nebelsick, J. H. 1992. Component analysis of sediment composition in Early Miocene temperate carbonates from the Austrian Paratethys. Palaeogeography, Palaeoclimatology, Palaeoecology 91:59–69CrossRefGoogle Scholar
[OBIS] Ocean Biogeographic Information System. 2021. Home page. www.obis.org, accessed 1 November 2021.Google Scholar
O'Dea, A., Lessios, H. A., Coates, A. G., Eytan, R. I., Restrepo-Moreno, S. A., Cione, A. L., Collins, L. S., Queiroz, A. D, Farris, D. W., Norris, R. D., Stallard, R. F., Woodburne, M. O., Aguilera, O., Aubry, M., Berggren, W. A., Budd, A. F., Cozzuol, M. A., Coppard, S. E., Duque-Caro, H., Finnegan, S., Gasparini, G., Grossman, E. L., Johnson, K. G., Keiwin, L. D., Knowlton, N., Leigh, E. G., Leonard-Pingel, J. S., Marko, P. B., Pyenson, N. D., Rachello-Dolmen, P. G., Soibelzon, E., Soibelzon, L, Todd, J. A., Vermeij, G. J., and Jackson, J. B.. 2016. Formation of the Isthmus of Panama. Science Advances 2(8):e1600883.CrossRefGoogle ScholarPubMed
Paillard, A., Shimada, K., and Pimiento, C.. 2021. The fossil record of extant elasmobranchs. Journal of Fish Biology 98:445455.CrossRefGoogle ScholarPubMed
Piller, W. E., Harzhauser, M., and Mandic, O.. 2007. Miocene Central Paratethys stratigraphy—current status and future directions. Stratigraphy 4:13.Google Scholar
Pimiento, C., and Benton, M. J.. 2020. The impact of the Pull of the Recent on extant elasmobranchs. Palaeontology 63:369374.CrossRefGoogle Scholar
Pimiento, C., Griffin, J. N., Clements, C. F., Silvestro, D., Varela, S., Uhenm, M. D., and Jaramillo, C.. 2017. The Pliocene marine megafauna extinction and its impact on functional diversity. Nature Ecology and Evolution 1:11001106.CrossRefGoogle ScholarPubMed
Pollerspöck, J., and Straube, N.. 2021. SharkReferences home page. www.shark-references.com, accessed 1 November 2021.Google Scholar
R Core Team. 2021. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna. http://www.R-project.org.Google Scholar
Reinecke, T., Louwye, S., Havekost, U., and Moths, H.. 2011. The Elasmobranch fauna of the Late Burdigalian, Miocene, at Werder-Uesen, Lower Saxony, Germany and its relationships with Early Miocene faunas in the North Atlantic, central Paratethys and Mediterranean. Palaeontos 20:1170.Google Scholar
Rivadeneira, M. M., and Marquet, P. A.. 2007. Selective extinction of late Neogene bivalves on the temperate Pacific coast of South America. Paleobiology 33:455468.CrossRefGoogle Scholar
Rögl, F. 1998. Palaeogeographic considerations for Mediterranean and Paratethys seaways (Oligocene to Miocene). Annalen Naturhistorischen Museums in Wien 99A:279310.Google Scholar
Rögl, F. 1999. Mediterranean and Paratethys. Facts and hypotheses of an Oligocene to Miocene paleogeography (short overview). Geologica Carpathica 50:339349.Google Scholar
Schultz, O. 2013. Catalogus Fossilium Austriae. Vol. 3, Pisces. Verlag der Österreichischen Akademie der Wissenschaften, Vienna.CrossRefGoogle Scholar
Schwab, J. 2015. Paleoecology and diversity of Middle Miocene Elasmobranchii of the Paratethys in the Vienna Basin. Bachelor's thesis, University of Vienna, Vienna.Google Scholar
Schwarzhans, W. 2010. The otoliths from the Miocene of the North Sea Basin. Backhuys, Weikersheim, Germany.Google Scholar
Steininger, F. F., and Wessely, G.. 2000. From the Tethyan Ocean to the Paratethys Sea: Oligocene to Neogene stratigraphy, paleogeography and paleobiogeography of the circum-Mediterranean region and the Oligocene to Neogene Basin evolution in Austria. Mitteilungen der Österreichischen Geologischen Gesellschaft 92:95116.Google Scholar
Tittensor, D.P., Mora, C., Jetz, W., Lotze, H. K., Ricard, D., Berghe, E. V., and Worm, B.. 2010. Global patterns and predictors of marine biodiversity across taxa. Nature 466:10981101.CrossRefGoogle ScholarPubMed
Todd, J. A., Jackson, J. B., Johnson, K. G., Fortunato, H. M., Heitz, A., Alvarez, M., and Jung, P.. 2002. The ecology of extinction: molluscan feeding and faunal turnover in the Caribbean Neogene. Proceedings of the Royal Society of London B 269:571577.CrossRefGoogle ScholarPubMed
Underwood, C. J., and Schlögl, J.. 2013. Deep-water chondrichthyans from the Early Miocene of the Vienna Basin (Central Paratethys, Slovakia). Acta Palaeontologica Polonica 58:487509.Google Scholar
Vertino, A., Stolarski, J., Bosellini, F. R., and Taviani, M.. 2014. Mediterranean corals through time: from Miocene to Present. Pp. 257274 in Goffredo, S. and Dubinsky, Z., eds. The Mediterranean Sea. Springer, Dordrecht, Netherlands.CrossRefGoogle Scholar
Villafaña, J. A., and Rivadeneira, M. M.. 2014. Rise and fall in diversity of Neogene marine vertebrates on the temperate Pacific coast of South America. Paleobiology 40:659674.CrossRefGoogle Scholar
Villafaña, J. A., and Rivadeneira, M. M.. 2018. The modulating role of traits on the biogeographic dynamics of chondrichthyans from the Neogene to the present. Paleobiology 44:251262.CrossRefGoogle Scholar
Villafaña, J. A., Marramà, G., Klug, S., Pollerspöck, J., Balsberger, M., Rivadeneira, M. M., and Kriwet, J.. 2020. Sharks, rays and skates (Chondrichthyes, Elasmobranchii) from the Upper Marine Molasse (middle Burdigalian, early Miocene) of the Simssee area (Bavaria, Germany), with comments on palaeogeographic and ecological patterns. Paläontologische Zeitschrift 94:725757.CrossRefGoogle ScholarPubMed
Westerhold, T., Marwan, N., Drury, A. J., Liebrand, D., Agnini, C., Anagnostou, E., Barnet, J. S. K., Bohaty, S. M., De Vleeschouwer, D., Florindo, F., Frederichs, T., Hodell, D. A., Holbourn, A. E., Kroon, D., Lauretano, V., Littler, K., Lourens, L. J., Lyle, M., Pälike, H., Röhl, U., Tian, J., Wilkens, R. H., Wilson, P. A., and Zachos, J.. 2020. An astronomically dated record of Earth's climate and its predictability over the last 66 million years. Science 369:13831387.CrossRefGoogle ScholarPubMed
Figure 0

Figure 1. Fossil occurrences of Neogene chondrichthyans from the Atlantic, Mediterranean, North Sea, and Paratethys regions.

Figure 1

Table 1. Faunal composition of Neogene European chondrichthyans at superorder, order, and family levels.

Figure 2

Figure 2. Diversification trajectories of European chondrichthyans. A, E, I, Number of occurrences; B, F, J, standing generic richness; C, G, K, origination rates; and D, H, L, extinction rates. The confidence intervals around the lines were excluded for clarity. Abbreviations: early Miocene (eM), middle Miocene (mM), late Miocene (lM) and Pliocene (Pl). Quorums used to estimate genus richness based on SQS are represented with black lines

Figure 3

Table 2. Spearman moment-product correlation for diversity trajectories between paired regions.a

Figure 4

Figure 3. Percentages of extinct/living chondrichthyan genera at each region, at global scale (i.e., the taxon is no longer present in the global ocean) and regional scale (i.e., the taxon was extirpated from the region but is still living elsewhere).

Figure 5

Table 3. Differences in the proportions of extinct genera at global and regional scales between regions using χ2 test.a

Figure 6

Figure 4. Climatic, oceanographic and biotic events through the Neogene and their relation with the origination (blue arrows) and extinction (red arrows) of Mediterranean (M), North Sea (NS), and Paratethys (P) chondrichthyans. 1: Westerhold et al. (2020); 2: Abreu and Haddad (1998); 3: Nebelsick (1992); 4: Kroh (2007); 5: Gürs and Janssen (2002); 6: Rögl (1999); 7: Harzhauser et al. (2007); 8: Reinecke et al. (2011); 9: Harzhauser et al. (2003); 10: Harzhauser and Piller (2007), 11: Piller et al. (2007); 12: Bannikov et al. (2018); 13: Hohenegger et al. (2014); 14: Krijgsman et al. (2000); 15: Krijgsman et al. (2010); 16: Todd et al. (2002); 17: Rivadeneira and Marquet (2007); 18: Villafaña and Rivadeneira (2014). CP, Central Paratethys; NS, North Sea.