Skip to main content

Global shifts in species richness have shaped carpet shark evolution

Abstract

Background

The evolutionary processes that shape patterns of species richness in marine ecosystems are complex and may differ between organismal groups. There has been considerable interest in understanding the evolutionary processes that led to marine species richness being concentrated in specific geographical locations. In this study we focus on the evolutionary history of a group of small-to-medium sized sharks known as carpet sharks. While a few carpet shark species are widespread, the majority of carpet shark species richness is contained within a biodiversity hotspot at the boundary of the Indian and Pacific oceans. We address the significance of this biodiversity hotspot in carpet shark evolution and speciation by leveraging a rich fossil record and molecular phylogenetics to examine the prehistoric distribution of carpet sharks.

Results

We find that carpet sharks species richness was greatest in shallow seas connected to the Atlantic Ocean during the Late Cretaceous, but that there was a subsequent loss of biodiversity in Atlantic waters. Fossil evidence from sites in close geographic proximity to the current center of carpet shark diversity are generally restricted to younger geologic strata.

Conclusions

From this data we conclude that (1) center of carpet shark biodiversity has shifted during the last 100 million years, (2) carpet sharks have repeatedly dispersed to nascent habitat (including to their current center of diversity), and (3) the current center of carpet shark biodiversity conserves lineages that have been extirpated from this prehistoric range and is a source of new carpet shark species. Our findings provide insights into the roles of marine biodiversity hotspots for higher-tropic level predators and the methods applied here can be used for additional studies of shark evolution.

Peer Review reports

Background

Geographical regions that contain a rich taxonomic diversity, across multiple taxonomically distant groups, with many endemic species are known as biodiversity hotspots [1]. The assembly of diverse species communities in marine biodiversity hotspots has provided insights into the evolution of marine species and ecosystems [2,3,4]. The region at the boundary of the Indian and Pacific oceans represents a large marine biodiversity hotspot, housing a rich diversity of shallow-water marine and reef species [2, 4,5,6]. There has been considerable interest in understanding the evolutionary and biogeographical processes that have shaped species diversity and composition in this region [2,3,4, 6]. This region appears to have been a cradle of biodiversity, by facilitating speciation during the Neogene [2,3,4, 6]. However, some of the lineages predate the formation of the biodiversity hotspot at the boundary of the Indian and Pacific oceans, suggesting these lineages immigrated to the region and then speciated [2]. Phylogenetic analysis and fossil evidence suggests that similar marine biodiversity hotspots were present near what is now southwest Europe during the Paleogene and that many of these species-groups present there immigrated into the current center of biodiversity during the Neogene [2]. Much of this work focused on coastal plants or small marine animals with rich fossil records [2]. However, recent studies suggest that marine predators may have followed a similar evolutionary path, diversifying elsewhere and subsequently arriving at the boundary of the Indian and Pacific oceans, where they then radiated [7, 8]. Here we examine the evolutionary history of carpet sharks (Chondrichthyes: Orectolobiformes) to assess their prehistoric distribution and estimate the age of species divergence times to provide insights into the colonization of a biodiversity hotspot by small-to-medium sized marine predators.

Carpet sharks represent a monophyletic group that encompasses a rich diversity of extant [7, 9,10,11,12] and fossil species [13, 14]. Extant carpet shark species (45 extant species within 7 families) constitute mostly small-to-medium-sized sharks (< 3.2 m long) and inhabit warm to temperate waters of continental and insular shelves [9, 15]. The only exception, the pelagic whale shark (Rhincodon typus), is the largest fish in the world [10, 11, 16]. Most carpet shark species have teeth that can rotate to accommodate either soft or hard-bodied prey [17] and this tooth morphology is unique among sharks [13, 18] (E. Manning pers. comm.). This has allowed many carpet shark species to prey on both small fishes and cephalopods (soft-bodied) as well as invertebrates with hard shells (hard-bodied), taking advantage of a variety of dietary resources [19, 20]. However, some groups of carpet sharks have become dietary specialists. The whale shark is largely planktivorous [19] and has reduced teeth [13, 14, 18]. The wobbegongs (Orectolobidae) and blind sharks (Brachilaridae) are higher-level predators that consume bony fishes, cephalopods, and smaller sharks [21, 22]. The wobbegong’s teeth are modified for a diet of soft-bodied prey, while blind sharks retain the characteristic rotating teeth [17]. The greatest species richness of carpet sharks is centered within the eastern Indian and western Pacific oceans, with a total of 41 species of carpet sharks are present in this region, 39 species that are endemic to one, or both of these oceans [10, 23, 24]. The richness of carpet shark species within a marine biodiversity hotspot and a rich fossil record made carpet sharks ideal to study the role of biodiversity hotspots in the speciation of marine predators.

In order to describe the evolutionary history of carpet sharks, we first needed estimate species divergence times. The process of estimating species divergence times from molecular sequence data requires external information to calibrate the rate at which new substitutions arise and are fixed. Fossils provide the most reliable source of information to calibrate DNA substitution rates and estimate divergence times [25]. The fossilized remains of shark teeth are commonly encountered in paleontological studies [13, 14], because a typical shark loses tens of thousands of teeth during a lifetime [26]. A few studies have used fossil shark evidence described in the primary literature to obtain divergence calibration times [8, 27,28,29,30,31,32]. The results of these studies suggest the literature contains useful and well-founded information on fossilized shark teeth that can be used to calibrate molecular phylogenies; however, these studies highlight issues around selecting the optimal fossil to calibrate the rate of DNA substitutions among many potentially informative fossils. These studies either selected the oldest fossil within a specious and evolutionarily old clade [27,28,29,30,31,32] or hand selected a few fossils representing different chronologically distant constraints [32]. The first approach excludes many potentially informative fossils and the later could lead to inconsistencies between studies. Here we employ a simple and repeatable method for quickly screening hundreds of fossil descriptions from the primary literature to find the optimal fossils for calibrating when to species of sharks diverged in geological time, while avoiding the time-consuming re-examination of fossils, and avoiding the assumption that all fossils are correctly identified within the literature. We then used these fossils to build a time-calibrated phylogenetic tree, comparing, in three cases, estimates of species divergence times to expected divergence times to evaluate our approach. Once we had a phylogenetic tree representing both species relationships and species divergence times, we needed to identify prehistoric centers of biodiversity. We did this by classifying to family the hundreds of fossils evaluated for the molecular calibration. From this data we identified regions that contained the greatest family level diversity at different points in the past.

Results

We complied 345 fossil records, each representing a unique combination of locality, age, and species (338 published records, 7 unpublished specimens; a full account of each record and its sources are included in the supplementary data). This was an attempt at creating an exhaustive inventory of fossil carpet shark records and only duplicate records were excluded. Species attributed to multiple carpet shark families were present in the Late Cretaceous deposits of North America, Europe, and Western Asia (Fig. 1A). Fossil descriptions from Paleogene of North America were dominated reports attributed to the family Ginglymostomatidae, with fossils attributed to the Brachaeluridae and Orectolobidae from Virginia and Alabama respectively, being the two exceptions, suggesting a decline in diversity near the end of the Cretaceous period. Despite a decline in North American reports, fossils from the Paleogene of Europe were attributed to species in multiple families. We also noted reports of fossils attributed to species from multiple families from the Paleogene deposits of northwestern and west-central Africa. Fossil descriptions from central or east Asia were few and included two Late Cretaceous reports from central Asia and two Paleogene reports, one from southern Asia, and the other from east Asia. Reports of fossils attributed to family Ginglymostomatidae are common throughout the Late Cretaceous and Paleogene fossil record. With the exception of two fossils described from southeast Asia attributed to the family Orectolobidae, all fossils obtained from Neogene deposits were attributed as species belonging to the Ginglymostomatidae. Fossils attributed to the pelagic family Rhincodontidae were limited to the Paleogene or Neogene deposits, across multiple continents (Fig. 1B). Collectively, centers of carpet shark diversity appear to have shifted while members of the Ginglymostomatidae and Rhincodontidae became and remained established throughout the world.

Fig. 1
figure 1

Geographical and chronological distribution of orectolobiform fossils attributed to extant families from the Late Cretaceous to recent. Fossils older than 100 MYA, fossils not attributed to a family, and suspect records (fossil identifiers 220, 273, 294) were excluded

We constructed a phylogenetic tree representing all seven carpet shark families and 13 genera. We did not capture the full species diversity within two families, Hemiscylliidae and Parascylliidae, but we did not exclude any major groups or species with an atypical biogeographical range. This tree was largely consistent with other phylogenetic results. Bootstrap support was greater than or equal to 75 % for most nodes in this tree. Nodes with lower bootstrap support included those node connecting species within families (Hemiscylliidae and Orectolobidae) and the node connecting Stegostomatidae and Rhincodontidae. Of the 345 fossils descriptions obtained for this study, 298 could readily be assigned to a species divergence with our phylogenetic tree. From each of the 298 records we attempted to build 298 time-calibrated phylogenetic trees using penalized likelihood [33, 34]. From these trees we extracted node ages for the three nodes representing divergences at different time frames (Table 1). We then determined the node age range and excluded any tree, along with its associated fossil record, that was a statistical outlier. This left us with a population of trees and their associated fossil records. Of these remaining records we selected the oldest fossil and used its minimum geological age as our calibration point for that node, leaving us with four fossil records assigned to four nodes. A final calibrated phylogenetic tree was produced using these four fossils. Mean node ages were reached in two identical runs, suggesting a stable result given our fossil data. We then examined three nodes in the tree where we had an expectation of node ages, finding that the mean node ages fit within our expectations and compared them to single point calibrations (Table 1). Node ages based on single point calibrations were consistent with our expected divergence times in two of three instances and all three predicted node ages were younger than ages obtained from multiple point calibrations. By “weeding out” uninformative or problematic fossils through a systematic process, we produced an age-calibrated phylogenetic tree of carpet sharks that was consistent with our expectation of divergence times at three different nodes. Therefore, at least with regards to carpet sharks, our method of culling and selecting fossils appears valid and provides improvement over single point calibrations.

Table 1 Expected and observed species divergence times in MYA

We arrived at a Triassic or Jurassic origin of carpet sharks (Fig. 2) with a mean age of 203.7 MYA. Extant families of orectolobiform sharks diverged from one another prior to the Cretaceous–Paleogene (K–Pg) boundary, except for Rhincodontidae and Stegostomatidae, which diverged near the time of the K–Pg boundary. There was subsequent speciation following the K–Pg boundary, including a radiation event within the Orectolobidae during the Neogene.

Fig. 2
figure 2

Time-calibrated phylogeny of carpet sharks based on mitochondrial DNA. Red arrows denote fossil-calibrated nodes used in the analysis. Blue bars and adjacent labels denote the confidence interval obtained using the program MCMCTree to sample the node ages. Values in italics represent bootstrap support for nodes. Colored boxes at tree tips represent modern ranges for species (blue = eastern Pacific Ocean, green = Atlantic Ocean, yellow = south western Indian Ocean, blue = north western Indian Ocean, black = central Indian Ocean, red = eastern Indian and western Pacific oceans north of Weber’s line, and Purple = eastern Indian and western Pacific oceans south of Weber’s line. Geologic periods are noted using gray and white boxes and associated text. Major clades of sharks are described by their common names and representative body form to the right of the tree. Globe created using rnaturalearth (https://github.com/ropensci/rnaturalearth)

Discussion

Today, the greatest species richness of carpet sharks is centered within the eastern Indian and western Pacific oceans, a region that has been shaped by the collision of crustal plates during the Neogene (represented by red and purple regions in Fig. 2) [7, 35, 36]. A total of 41 species of orectolobiform sharks are present in this region including 39 species that are endemic to one or both of these oceans [10, 23, 24]. We were only able to identify a few fossils described from sites in close geographical proximity to this current center of carpet shark diversity. However, we identified locations distant to the current center of biodiversity, that contained diverse assemblages of fossil carpet shark species, representing multiple families, during the Late Cretaceous and Paleogene. Based on the data assembled here, we propose that the center of carpet shark diversity shifted to the eastern Indian and western Pacific oceans during the Cenozoic.

Fossil descriptions indicate two regions were host to diverse assemblages of carpet shark species during the Late Cretaceous. Fossil from western Eurasia suggest that elements of the Tethys seaway hosted a diverse assemblage of carpet shark species, representing multiple families [37]. Fossil evidence also suggests a diverse assemblage of carpet shark species were present in the western interior seaway and inundated Atlantic coastal plain of North America during the same time period [38,39,40]. Late Cretaceous shallow seas in these regions likely provided suitable habitat for carpet sharks [39, 40] and facilitated speciation. The end of the Cretaceous was marked by a decline in carpet shark diversity in North America, while diversity remained high in what is now Europe. Fossil evidence also shows a diverse assemblage of carpet shark species in west central and north Africa during the Paleogene, also likely facilitated by suitable habitat that was available at this time [41]. While fossil data is often patchy and does not fully capture the timing and extent of invasion and speciation within a region, the data clearly demonstrates that a diverse carpet shark assemblage previously occupied the Tethyan seaway and shallow seas connected to the Atlantic basin. Currently, the Atlantic basin is home to two carpet shark species, the nurse shark (G. cirratum) and the pelagic whale shark (R. typus). We conclude that the Atlantic basin was likely a cradle of carpet shark biodiversity, but that there was subsequent carpet shark extinction during the Paleogene. Collectively, the data suggests carpet sharks are capable of long-distance dispersal and colonizing nascent habitat.

The eastern Indian and western Pacific oceans have conserved carpet shark lineages, while promoting speciation within some of these lineages. Wobbegongs and blind and collared sharks are restricted to the eastern Indian and western Pacific oceans, but fossil and phylogenetic evidence suggest these groups likely evolved elsewhere [8]. Therefore, the eastern Indian and western Pacific oceans have acted as refuge following extinction of these groups in other regions. Within the eastern Indian and western Pacific oceans, the waters off Australia and New Guinea, comprising the Sahul shelf (purple region in Fig. 2), appear to have played a particularly important role in the evolution of carpet sharks [7, 8]. The waters over the Sahul shelf are host to most species of wobbegongs. Only two closely related species (O. japonicus and O. leptolineatus), along with Eucrossorhinus dasypogon, occur outside of Sahul shelf, to the north on the Sunda shelf [9, 42] and our phylogenetic inference suggests these species dispersed out of the Sahul shelf. Therefore, we conclude that these wobbegong species, until recently, may have been restricted to the Sahul shelf. Other carpet shark groups, including blind, epaulette (a highly specious lineage within the family Hemiscylliidae), and collared (in part) sharks are also restricted to the Sahul shelf and adjacent islands [8, 9, 15]. Therefore, we conclude that a relatively small geographical region has played a large role in conserving carpet shark diversity. Previous phylogenetic studies have also highlighted the Sahul shelf as a cradle of biodiversity [7, 8]. Within this region, both wobbegongs and epaulette sharks have speciated extensively, likely a result of allopatric speciation within the region [7, 8].

While it appears that dispersal, radiation, and extinction have shaped wobbegong, blind, and collared shark lineages, three groups, nurse, zebra, bamboo, and whale sharks, appear to defy this pattern. We noted that fossil nurse shark remains are present throughout the world during the Late Cretaceous, Paleogene, and Neogene. Furthermore, nurse shark remains have been identified after other carpet shark species disappear from the fossil record. While extant nurse shark diversity is small, the three recognized species collectively occupy a wide geographical region. Collectively, this suggests that nurse sharks acquired one or more traits that allowed them to persist in environments that led to the localized extinction of other carpet shark groups. The whale shark, and allied extinct species, shifted to planktonic feeding [13, 14]. Whale shark teeth have been found in many fossil sites including Africa, Eurasia, and the New World. The single extant whale shark species is found throughout the world [20]. The zebra shark ranges widely in the Indian Ocean, unfortunately, fossils records were few and we cannot comment at this time about its prehistoric distribution. Finally, the fossil remains of bamboo sharks have been recovered in Europe and North America. While extant species are absent from North America, the whitespotted (C. plagiosum) and slender (C. indicum) bamboo sharks are found through much of the Indian Ocean. Thus, these sharks appear to have escaped the evolutionary pressures that have shaped the remainder of carpet shark groups.

Conclusions

The teeth of sharks, and closely related rays, are abundant in the fossil record. When used in combination with molecular phylogenetics, these fossil teeth present an opportunity to examine the prehistoric distribution of sharks and rays. We utilized this data here to explore the evolutionary and biogeographic processes that shaped carpet shark species richness. We suggest our study on carpet sharks could serve as a template for future studies utilizing the rich fossil evidence left by shark and ray species.

Methods

Molecular phylogenetics

Molecular [29, 43] and morphological studies [24] suggest carpet sharks form a clade within the superorder Galeomorphii. Therefore, the first step was to build a phylogenetic tree representing extant species diversity within carpet sharks. We used DNA sequence data from the mitochondrial gene NADH2, which provided us with data from all seven families and 13 genera of carpet sharks. We did not capture the full species diversity within the families Hemiscylliidae and Parascylliidae; however, we did not exclude any species having an atypical range compared to closely related species. DNA sequences were aligned as amino acid sequences using the multiple sequence aligner Muscle (v3.7; BLOSSUM62 substitution matrix) [44] and visualized in Geneious (v2019.2.1) for a manual check of alignment quality before back-translating to nucleotide. PartitionFinder2 (v2.1.1) [45] was used to identify the optimal model of sequence evolution, using Akaike information criterion modified for small sample size [46], for each codon position. Each codon position was separated into a partition and the GTR + Γ + I model was applied to each partition. Free-parameters were estimated for each codon position independently. The optimal phylogenetic arrangement was reached under a Maximum Likelihood (ML) optimality criterion using RAxML (v8.2.12) [47]. We validated our tree in two ways. First, the species arrangements in our ML tree were compared with previous findings. Second, we assessed support among 1000 bootstrap replicates. An ML tree was calculated for each bootstrap replicate under the same model of sequence evolution, again allowing free-parameters of the model to be calculated for each codon position.

Fossil data collection and selection

We obtained fossil data by reviewing the relevant primary literature. Because some geographic regions were under-represented in our data, we made a special effort to seek out published reports from these regions. We also visited a large collection of fossil shark remains with well documented origins (G. Hubbell coll.) to supplement the published data, particularly in regards to filling-in data gaps in under-represented regions. Most fossils were described using geologic age and were converted to coarse age ranges in MYA following Walker et al. [48]. It was appropriate to use coarse age ranges as fossil teeth often aggregate in lag deposits that can prevent assignment to a narrow age range. We obtained 345 fossil records in total.

Carpet sharks have teeth that rotate to accommodate soft or hard-bodied prey [17] and the associated tooth morphology is unique among sharks [13, 18] (E. Manning pers. comm.). Within carpet sharks, the upper margin of the labial tooth face differs between families [18]. These differences facilitate the family-level classification of teeth and suggests most fossil descriptions likely represent accurate fossil classification to an extant family and not an extinct sister lineage. The loss of diagnostic tooth morphology in a clade could be problematic for incorporating fossil remains. The wobbegongs (Orectolobidae) and the whale shark (Rhincodontidae) have atypical novel tooth morphology associated with novel diets [9, 21, 49]. Wobbegong teeth superficially resemble teeth from angel sharks (Squatiniformes: Squatinidae), but characters along the upper margin can be used to assign teeth to Orectolobidae. Whale shark teeth are greatly reduced, [18] resembling teeth from other filter-feeding sharks, but root structure differs among filter-feeding shark species [13]. Given the unique morphology, we expect that teeth described in the literature were in most cases correctly identified and each fossil was assigned to a node in our tree based on its described taxonomic classification (Fig. 2). Fossils with no family level assignment reported and fossils from the genus Pararhincodon, which appears to represent an extinct clade, were left unassigned (47 records). Our phylogenetic reconstruction suggested the family Ginglymostomatidae is polyphyletic; therefore, fossils assigned to the genus Pseudoginglymostoma were grouped with fossils from the family Stegostomatidae and were considered for calibration of node 5.

We developed a bioinformatics process to rapidly evaluate each fossil record individually and identify problematic fossils. First, each fossil record was used to convert our ML tree into a chronogram based on its node assignment and age range. Because we had to repeat the process of node calibration for each of the 298 fossil records, we used a ML method described by Sanderson [33] to quickly estimate a chronogram, rather than a more computationally expensive approach. This method used an optimized ML approach to apply a relaxed molecular clock, implemented in TreePL [34]. A script was written to prepare an input file for each fossil record (298 total), with each fossil assigned to a node (nodes labeled 2–7; Table 2). TreePL was run with each input file and then node ages were collected from each output file using the R package ape (R v3.6.1; ape v5.3) [50].

Table 2 Nodes used to calibrate the phylogenetic tree using treePL

Once we generated a rate-calibrated molecular phylogeny using data from each fossil record, we needed to exclude any potentially problematic or uninformative fossil records based on the resulting time-calibrated phylogeny. To do this, we excluded any fossil that provided a node age that was outside the interquartile region at one or more of three nodes in the phylogenetic tree. These three nodes were distributed across the phylogenetic tree and where we had an expectation of node age (Table 1). We examined median node ages to determine if they were appropriate (Table 1). The remaining 110 fossils were assigned to nodes 2–6, with no remaining fossils assigned to node 7. Finally, we removed any fossils that were assigned to a node having three or fewer fossils. This resulted in removing two fossil records assigned to node 3. 108 fossil records remained.

Fossil calibration

After we evaluated each fossil record individually and removed outliers, we identified the minimum age for the oldest fossil at each node (node 2 > 170 MYA, node 4 > 94 MYA, node 5 > 66 MYA, and node 6 > 113 MYA). Based on these four minimum calibrations of divergence times, we estimated median node ages and 95 % confidence intervals for all nodes within our phylogenetic tree using MCMCTree (v4.9j; 500,000 generations, sampling every 50 generations with a burn-in of 50,000 generations) [51]. Baseml (v4.7a) [51] was used to estimate substitution rates (0.290696) following Stein et al. [12]. Root age of our analysis was constrained to < 334 MYA, and an independent rates model was used. Mate selection, reproductive strategies, and age at maturity vary between species in this group [5, 20, 24, 52, 53]; therefore, we believed it was best to use the independent rates option available in MCMCTree. We examined results using Tracer (v.1.7.1) [54] and then replicated the analysis to check for median node age convergence.

Availability of data and materials

No public or private databases were used in this study to obtain fossil data; sequence data was obtained from https://www.ncbi.nlm.gov/genbank/. Fossil data was obtained from the primary literature or from accession data from the Gordon Hubbell collection with permission. Raw data is available from the figshare data repository https://doi.org/10.6084/m9.figshare.12595169. This repository includes all fossil records gathered from the literature and the Hubbell collection, aligned sequences, sequence partitions, ML tree with bootstrap values, individual calibrated phylogenetic trees from implementation of TreePL, summary data on node ages from each individual calibrated trees (both in total and with outliers removed), baseml control file, MCMCTree control file, MCMCTree calibration tree, calibrated phylogenetic tree based on multiple node calibration points, convergence plot for calibrated phylogenies, biogeographic model comparisons, ancestral range estimations, fossil data and node associations, node definition table, and range data for extant species.

References

  1. Reid WV. Biodiversity hotspots. Trend Ecol Evol. 1998;13:275–80.

    Article  CAS  Google Scholar 

  2. Renema W, Bellwood DR, Braga JC, Bromfield K, Hall R, Johnson KG, Lunt P, Meyer CP, McMonagle LB, Morley RJ, et al. Hopping hotspots: global shifts in marine biodiversity. Science. 2008;321:654–7.

    Article  CAS  PubMed  Google Scholar 

  3. Palumbi SR. Molecular biogeography of the Pacific. Coral Reefs. 1997;16:S47–52.

    Article  Google Scholar 

  4. Read CI, Bellwood DR, Herwerden L. Ancient origins of Indo-Pacific coral reef fish biodiversity: a case study of the leopard wrasses (Labridae: Macropharyngodon). Mol Phylo Evol. 2006;38:808–19.

    Article  CAS  Google Scholar 

  5. Stehli FG, Wells JW. Diversity and age patterns in hermatypic corals. Syst Biol. 1971;20:115–26.

    Google Scholar 

  6. Bellwood DR, Hughes TP. Regional-scale assembly rules and biodiversity of coral reefs. Science. 2001;292:1532–5.

    Article  CAS  PubMed  Google Scholar 

  7. Corrigan S, Beheregaray LB. A recent shark radiation: molecular phylogeny, biogeography and speciation of wobbegong sharks (Family: Orectolobidae). Mol Phylo Evol. 2009;52:205–16.

    Article  CAS  Google Scholar 

  8. Dudgeon CL, Corrigan S, Yang L, Allen GR, Erdmann MV, Fahmi, Sugeha HY, White WT, Naylor GJP. Walking, swimming or hitching a ride? Phylogenetics and biogeography of the walking shark genus Hemiscyllium. Marine Freshwater Res. 2020. https://doi.org/10.1071/MF19163.

  9. Last PR, Stevens JD. Sharks and Rays of Australia, Second Edition. Harvard University Press, MA, 2009.

  10. Weigmann S. Annotated checklist of the living sharks, batoids and chimaeras (Chondrichthyes) of the world, with a focus on biogeographical diversity. J Fish Biol. 2016;88:837–1037.

    Article  CAS  PubMed  Google Scholar 

  11. Fricke R, Eschmeyer WN, Van der Laan R. Eschmeyer’s Catalog of Fishes. 2020. Accessed 22 July 2020 online at http://researcharchive.calacademy.org/research/‌ichthyology/‌catalog/fishcatmain.asp

  12. Castro JI. The Sharks of North America. New York: Oxford University Press; 2011.

    Google Scholar 

  13. Cappetta H. "Chondrichthyes II, Mesozoic and Cenozoic Elasmobranchii, Volume 3B”. In: Schultze HP, editor. Handbook of Paleoichthyology. Stuttgart: Gustav Fischer; 1987.

    Google Scholar 

  14. Cappetta H. “Chondrichthyes, Mesozoic and Cenozoic Elasmobranchii: teeth, Volume 3E.” In: Schultze HP, editor. Handbook of Paleoichthyology. Stuttgart: Gustav Fischer; 2012.

    Google Scholar 

  15. Allen GR, Erdmann MV, White WT, Fahmi, Dudgeon CL. Review of the bamboo shark genus Hemiscyllium (Orectolobiformes: Hemiscyllidae). J Ocean Sci Found. 2016;23:51–97.

    Google Scholar 

  16. Rowat D, Brooks KS. A review of the biology, fisheries and conservation of the whale shark Rhincodon typus. J Fish Biol. 2012;80:1019–56.

    Article  CAS  PubMed  Google Scholar 

  17. Ramsay JB, Wilga CD. Morphology and mechanics of the teeth and jaws of white-spotted bamboo sharks (Chiloscyllium plagiosum). J Morph. 2007;268:664–82.

    Article  PubMed  Google Scholar 

  18. Herman JM, Hovestadt-Euler M, Hovestadt DC. Contributions to the study of the comparative morphology of teeth and other relevant ichthyodorulites in living supraspecific taxa of Chondrichthyan fishes. Part A: Selachii. No. 4: Order: Orectolobiformes, Families: Brachaeluridae, Ginglymostomatidae, Hemiscylliidae, Orectolobidae, Parascylliidae, Rhiniodontidae, Stegostomatidae. Order: Pristiophoriformes–Family: Pristiophoridae, Order: Squatiniformes–Family: Squatinidae. Bull L’Institut Royal Sci Nature Belgique. 1992;67:193–254.

    Google Scholar 

  19. Cortes E. Standardized diet composition and trophic levels of sharks. J Marine Sci. 1999;56:707–17.

    Google Scholar 

  20. Castro JI. The biology of the nurse shark, Ginglymostoma cirratum, off the Florida east coast and the Bahama islands. Environ Biol Fishes. 2000;58:1–22.

    Article  Google Scholar 

  21. Huveneers C, Otway NM, Gibbs SE, Harcourt RG. Quantitative diet assessment of wobbegong sharks (genus Orectolobus) in New South Wales, Australia. ICES J Mar Sci. 2007;64:1272–81.

    Article  Google Scholar 

  22. Kyne PM, Compagno LJV, Stead J, Jackson MV, Bennett MB. Distribution, habitat and biology of a rare and threatened eastern Australian endemic shark: Colclough’s shark, Brachaelurus colcloughi Ogilby, 1908. Marine Freshwater Res. 2011;62:540–7.

    Article  CAS  Google Scholar 

  23. Springer VG. Pacific Plate Biogeography, with Special Reference to Shorefishes. Smithsonian Contributions to Zoology No. 367. 1982.

  24. Goto T. Comparative anatomy, phylogeny and cladistic classification of the Order Orectolobiformes (Chondrichthyes, Elasmobranchii). Mem Grad School of Fisheries Sci, Hokkaido Uni. 2001;48:1–100.

    Google Scholar 

  25. Forest F. Calibrating the tree of life: fossils, molecules and evolutionary timescales. Ann Bot. 2009;104:789–94.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Luer CA, Blum PC, Gilbert PW. Rate of tooth replacement in the nurse shark, Ginglymostoma cirratum. Copeia. 1990;1:182–91.

    Article  Google Scholar 

  27. Heinicke MP, Naylor GJP, Hedges SB. Cartilaginous fishes (Chondricthyes). In: Hedges SB, Kumar S, editors. The Timetree of Life. Oxford: Oxford University Press; 2013. p. 320–7.

    Google Scholar 

  28. Puckridge M, Last PR, White WT, Andreakis N. Phylogeography of the Indo-West Pacific maskrays (Dasyatidae, Neotrygon): a complex example of chondrichthyan radiation in the Cenozoic. Ecol Evol. 2013;3:217–32.

    Article  PubMed  Google Scholar 

  29. Sorenson L, Santini F, Alfaro ME. The effect of habitat on modern shark diversification. J Evol Biol. 2014;27:1536–48.

    Article  CAS  PubMed  Google Scholar 

  30. Straube N, Li C, Claes JM, Corrigan S, Naylor GJP. Molecular phylogeny of Squaliformes and first occurrence of bioluminescence in sharks. BMC Evol Biol. 2015;15:162.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Stein RW, Mull CG, Kuhn TS, Aschliman NC, Davidson LNK, Joy JB, Smith GJ, Dulvy NK, Mooers AO. Global priorities for conserving the evolutionary history of sharks, rays and chimeras. Nat Ecol Evol. 2018;2:288–98.

    Article  PubMed  Google Scholar 

  32. Swift DG, Dunning LT, Igea J, Brooks EJ, Jones CS, Noble LR, Ciezarek A, Humble E, Savolainen V. Evidence of positive selection associated with placental loss in tiger sharks. BMC Evol Biol. 2016;16:1–10.

    Article  CAS  Google Scholar 

  33. Sanderson MJ. Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach. Mol Biol Evol. 2002;19:101–9.

    Article  CAS  PubMed  Google Scholar 

  34. Smith AA, O’Meara BC. treePL: divergence time estimation using penalized likelihood for larger phylogenies. Bioinformatics. 2012;28:2689–90.

    Article  CAS  PubMed  Google Scholar 

  35. Hall R. Cenozoic geological and plate tectonic evolution of SE Asia and the SW Pacific: computer-based reconstructions, model and animations. J Asian Earth Sci. 2002;20:353–431.

    Article  Google Scholar 

  36. Hall R. Australia–SE Asia collision: plate tectonics and crustal flow. In: Hall R, Cottam MA, Wilson MEJ, Eds. The SE Asian Gateway: History and Tectonics of the Australia–Asia Collision. Geol. Soc. Lond. Spec. Pub. 2011;355:75–109.

  37. Seton M, Muller RD, Zahirovic S, Gaina C, Torsvik T, Shepard G, Talsma A, Gurnis M, Turner M, Maus S, Chandler M. Global continental and ocean basin reconstructions since 200 Ma. Earth Sci Rev. 2012;113:212–70.

    Article  Google Scholar 

  38. Haq BU. Cretaceous eustasy revisited. Glob Planet Change. 2014;113:44–58.

    Article  Google Scholar 

  39. Murray AM, Cook TD. Overview of the Late Cretaceous fishes of the northern western interior seaway. In: Khosla A, Lucas SG, ed. Cretaceous period: biotic diversity and biogeography. New Mexico Museum of Natural History and Science Bulletin 2016;71: pp255-261.

  40. Nicholls EL, Russell AP. Paleobiogeography of the Cretaceous western interior seaway of North America: the vertebrate evidence. Paleogeogr Palaeoclimatol Palaeoecol. 1990;79:149–69.

    Article  Google Scholar 

  41. Sole F, Noiret C, Desmares D, Adnet S, Taverne L, Putter TD, Mees F, Yans J, Louwye S, Folie A, et al. Reassessment of historical sections from the Paleogene marine margin of the Congo Basin reveals an almost complete absence of Danian deposits. Geosci Front. 2019;10:1039–63.

    Article  Google Scholar 

  42. Ebert DA, Fowler S, Compagno L. Sharks of the World, a Fully Illustrated Guide. NH: Wild Nature Press; 2013.

    Google Scholar 

  43. Naylor GJP, Caira JN, Jensen K, Rosana KAM, White WT, Last PR. A DNA sequence-based approach to the identification of shark and ray species and its implications for global elasmobranch diversity and parasitology. Bull Am Mus Nat Hist. 2012, 1–262.

  44. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucl Acids Res. 2004;32:1792–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Lanfear R, Fradsen PB, Wright AM, Senfeld T, Calcott B. PartitionFinder2: new methods for selecting partitioned models of evolution for molecular and morphological analyses. Mol Biol Evol. 2017;34:772–3.

    CAS  PubMed  Google Scholar 

  46. Sullivan J, Joyce P. Model selection in phylogenetics. Ann Rev Ecol Evol Syst. 2014;36:445–66.

    Article  Google Scholar 

  47. Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of larger phylogenies. Bioinf. 2014;30:1312–1213.

    CAS  Google Scholar 

  48. Walker JD, Geisman JW, Bowring SA, Babcock LE. Geologic Time Scale v. 5.0, Geological Society of America. 2018.

  49. Borell A, Aguilar A, Gazo A, Kumarran RP, Cardona L. Stable isotope profiles in whale shark (Rhincodon typus) suggest segregation and dissimilarities in the diet depending on sex and size. Environ Biol Fishes. 2011;92:559–67.

    Article  Google Scholar 

  50. Paradis E, Claude J, Strimmer K. APE: analysis of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–90.

    Article  CAS  PubMed  Google Scholar 

  51. Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24:1586–91.

    Article  CAS  PubMed  Google Scholar 

  52. Pratt HL, Carrier JC. A review of elasmobranch reproductive behavior with a case study on the nurse shark, Ginglymostoma cirratum. Environ Biol Fishes. 2001;60:157–88.

    Article  Google Scholar 

  53. Nakaya K, White WT, Ho HC. Discovery of a new mode of oviparous reproduction in sharks and its evolutionary implications. Sci Reports. 2020;10:12280.

    CAS  Google Scholar 

  54. Rambaut A, Drummond AJ, Xie D, Baele G, Suchard MA. Posterior summarization in Bayesian phylogenetics using Tracer 1.7. Syst Biol. 2018;67:901–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. del Moral-Flores LF, Ramírez-Antonia E, Angulo A, Pérez-Ponce de León G. Ginglymostoma unami sp. nov. (Chondrichthyes: Orectolobiformes: Ginglymostomatidae): una especie nueva de tiburón gata del Pacifíco oriental tropical. Rev Mex Biodivers. 2015;86:48–58.

  56. O’dea A, Lessios HA, Coates AG, Eytan RI, Restrepo-Moreno SA, Cione AL, Collins LS, de Queiroz A, Farris DW, Norris RD, et al. Formation of the isthmus of panama. Sci Adv. 2016;2:e1600883.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Maisey JG. What is an ‘elasmobranch’? The impact of palaeontology in understanding elasmobranch phylogeny and evolution. J Fish Biol. 2012;80:918–51.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Gordon Hubbell (Gainesville, Florida) for providing fossil data and review of an earlier draft of this manuscript, Leila Alickovic (Virginia Commonwealth University) for help with coding and processing fossil and biogeographic data, Jürgen Pollerspöck (Shark-References.com) for his help in locating obscure bibliographic records, Earl Manning (Davenport, Iowa) for sharing his knowledge of fossil carpet shark teeth, and Connie Steen (ANAMAR) and Gregg Poulakis (Florida Fish and Wildlife Conservation Commission) for providing reviews of earlier drafts of this manuscript. High Performance Computing resources provided by the High Performance Research Computing (HPRC) Core Facility at Virginia Commonwealth University (https://chipc.vcu.edu) were used for conducting the research reported in this work.

Funding

This work was supported by the Virginia Commonwealth University, Department of Life Sciences and ANAMAR Environmental Consulting. Virginia Commonwealth University provided support through startup funds to BMB. ANAMAR Environmental Consulting provided JCS with time to assemble the data set. This study was initiated by the authors and neither organization influenced the study design or interpretation of the results.

Author information

Authors and Affiliations

Authors

Contributions

JCS aggregated and evaluated fossil and range data used in this study. BMB conducted phylogenetic and statistical analyses, evaluated phylogenetic results, and prepared the figures. BMB and JCS interpreted the results and prepared the manuscript. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Bret M. Boyd.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors have no competing interests to declare.

Additional information

Publisher’s Note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Boyd, B.M., Seitz, J.C. Global shifts in species richness have shaped carpet shark evolution. BMC Ecol Evo 21, 192 (2021). https://doi.org/10.1186/s12862-021-01922-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-021-01922-6

Keywords