Molecular systematics and global phylogeography of angel sharks (genus Squatina)

https://doi.org/10.1016/j.ympev.2009.07.029Get rights and content

Abstract

Angel sharks of the genus Squatina represent a group comprising 22 extant benthic species inhabiting continental shelves and upper slopes. In the present study, a comprehensive phylogenetic reconstruction of 17 Squatina species based on two mitochondrial markers (COI and 16S rRNA) is provided. The phylogenetic reconstructions are used to test biogeographic patterns. In addition, a molecular clock analysis is conducted to estimate divergence times of the emerged clades. All analyses show Squatina to be monophyletic. Four geographic clades are recognized, of which the Europe–North Africa–Asia clade is probably a result of the Tethys Sea closure. A second sister group relationship emerged in the analyses, including S. californica (eastern North Pacific) and S. dumeril (western North Atlantic), probably related to the rise of the Panamanian isthmus. The molecular clock analysis show that both lineage divergences coincide with the estimated time of these two geological events.

Introduction

The angel shark genus Squatina (Chondrichthyes: Elasmobranchii: Neoselachii: Squalea: Squatinidae) comprises 22 extant, morphologically homogenous, benthic species (cf. Vooren and Da Silva, 1992, Compagno et al., 2005, Castro-Aguirre et al., 2006, Last and White, 2008), which inhabit continental shelves and upper slopes down to 500 m (Compagno et al., 2005). They are moderately-sized (total length about 1–2 m) and globally distributed in temperate to tropical seas (Compagno et al., 2005, Last and White, 2008). While some species occur over a wide geographic range, the majority are restricted to a smaller area (Compagno et al., 2005). Restriction in geographic range might be as a result of the behavior of Squatina species, which are ambush predators with a corresponding stationary bottom-dwelling habit (Compagno et al., 2005). Thus, trans-ocean migration is extremely unlikely, even though large-scale coastal migratory patterns have been reported in species such as S. squatina (Wheeler et al., 1975) and S. californica (Kato et al., 1967, Eschmeyer et al., 1983, Natanson and Cailliet, 1986, Compagno et al., 2005).

Based on their distribution, eight distinct putative zoogeographic species groups can be distinguished (distribution areas based on Compagno et al., 2005, Last and White, 2008): (1) Eastern North Atlantic–Mediterranean–North Africa (S. aculeata, sympatric with S. oculata and S. squatina in the Mediterranean), (2) southwest Indian Ocean, South Africa (S. africana), (3) western North Pacific, i.e. Asian species (S. formosa, S. japonica, S. nebulosa and S. tergocellatoides), (4) western South Pacific and eastern Indian Ocean, i.e. Australian species (S. australis along the south coast of Australia, partly sympatric with S. albipunctata along south-eastern Australia and with S. tergocellata along south-western Australia; S. pseudocellata occurs along north-western Australia) and S. legnota, which is recorded from Indonesia, (5) eastern North Pacific (S. californica), (6) eastern South Pacific (S. armata), (7) western North Atlantic (S. dumeril, sympatric with S. heteroptera and S. mexicana in the Gulf of Mexico) and (8) western South Atlantic (S. argentina, sympatric with S. guggenheim, S. occulta and S. punctata). These species groupings can be used to formulate and test biogeographic hypotheses, such as the influence of the rise of the Panamanian isthmus and the Tethys Sea closure on speciation.

Interspecific relationships of angel sharks have not been investigated except for a study by Furtado-Neto and Carr (2002) on three Brazilian species. The position of the squatinids within the elasmobranchs (sharks, skates and rays) based on morphological characters is still widely debated (cf. Compagno, 1973, Maisey, 1984, Thies and Reif, 1985, Seret, 1986, Shirai, 1992, de Carvalho, 1996). Recent mtDNA and nDNA studies by Douady et al. (2003) and Winchell et al. (2004) show closer relationships between Squatiniformes (angel sharks), Pristiophoriformes (saw sharks) and Squaliformes (dogfish sharks), which is congruent with morphological data.

In the present paper we present the first comprehensive phylogeography of the squatinid sharks based on two mitochondrial markers, cytochrome oxidase subunit I (COI) and 16S rRNA, using different tree reconstruction methods. In addition, fossil and geological calibration points are used to estimate divergence times among the Squatina species and to test biogeographic patterns of inferred species groups.

Section snippets

Material examined

In this study, 17 of the 22 described Squatina species from different localities of all geographic species groups (see Section 1) were analyzed for molecular purposes (Table 1, Fig. 1). Tissue samples could not be obtained from Squatina argentina, S. heteroptera, S. mexicana, S. nebulosa or S. punctata.

In addition, 22 cytochrome oxidase subunit I (COI) sequences of four Australian and two Asian species of Squatina from Ward et al., 2005, Ward et al., 2008 were used as reference data for

Molecular phylogeny and geographic patterns

All three phylogenetic methods showed Squatina to be monophyletic and yielded essentially the same topology (Fig. 2). Four major clades were found in all analyses (COI, 16S rRNA and concatenated sequences). These clades correspond to geographic regions and comprise (1) the European and North African species S. aculeata, S. oculata and S. squatina, together with the Asian species S. formosa, S. japonica, S. legnota and S. tergocellatoides (clade 1, blue), (2) the South African species S. africana

Reliability of divergence time estimates

Estimation of divergence times using fossil taxa is, with good reason, hotly debated and frequently criticized (e.g. Heads, 2005; Pulquerio and Nichols, 2007). Squatina cranei was used, as it is the oldest known, verified fossil squatinid, to determine the tmrca of Squalus and Squatina. However, it cannot be excluded that the squatinid lineage is even older depending on the exact position of the fossil within Squatina (see e.g. Magallon, 2004), resulting in generally older mean node ages and

Acknowledgments

This study would not have been possible without the following colleagues, who provided tissue material and additional sequences: Julie Carter (NOAA, Charleston, USA), Sebastian Hernandez (Chile) and Nicolas Straube (Berlin/Munich), Mark McGrouther (Australian Museum, Sydney, Australia), Andrew Bentley (Kansas University, Lawrence, USA), Samuel Iglesias (Muséum national d’Histoire naturelle, Paris, France), Gavin Naylor (Florida State University, Tallahassee, USA), Juan M. Diaz de Astarloa

References (63)

  • C.J. Winchell et al.

    Phylogeny of elasmobranchs based on LSU and SSU ribosomal RNA genes

    Mol. Phylogenet. Evol.

    (2004)
  • J.L. Castro-Aguirre et al.

    Dos nuevas especies del genero Squatina (Chondrichthyes: Squatinidae) del Golfo de Mexico

    Rev. Biol. Trop.

    (2006)
  • L.J.V. Compagno

    Interrelationships of living elasmobranchs

  • L.J.V. Compagno et al.

    A Field Guide to the Sharks of the World

    (2005)
  • A. de Queiroz

    The resurrection of oceanic dispersal in historical biogeography

    Trends Ecol. Evol.

    (2005)
  • S. Doiron et al.

    A comparative mitogenomic analysis of the potential adaptive value of arctic charr mtDNA introgression in brook charr populations (Salvelinus fontinalis Mitchill)

    Mol. Biol. Evol.

    (2002)
  • D.J. Douady et al.

    Molecular phylogenetic evidence refuting the hypothesis of Batoidea (rays and skates) as derived sharks

    Mol. Phylogenet. Evol.

    (2003)
  • A.J. Drummond et al.

    BEAST: Bayesian evolutionary analysis by sampling trees

    BMC Evol. Biol.

    (2007)
  • J.A. Endler

    Geographic Variation, Speciation and Clines

    (1977)
  • W.N. Eschmeyer et al.

    A Field Guide to Pacific Coast Fishes of North America

    (1983)
  • O. Folmer et al.

    DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates

    Mol. Mar. Biol. Biotechnol.

    (1994)
  • Furtado-Neto, M., Carr, S., 2002. Molecular genetics of some Brazilian sharks. The IUCN/SSC Shark Specialist Group,...
  • W.S. Grant et al.

    Shallow population histories in deep evolutionary lineages of marine fishes: insights from sardines and anchovies and lessons for conservation

    J. Hered.

    (1998)
  • M. Heads

    Dating nodes on molecular phylogenies: a critique of molecular biogeography

    Cladistics

    (2005)
  • T. Hrbek et al.

    Closing of the Tethys Sea and the phylogeny of the Eurasian killifishes (Cyprinodontiformes: Cyprinodontidae)

    J. Evol. Biol.

    (2003)
  • J.P. Huelsenbeck et al.

    Phylogeny estimation and hypothesis testing using maximum likelihood

    Annu. Rev. Ecol. Syst.

    (1997)
  • J.P. Huelsenbeck et al.

    A likelihood-ratio test of monophyly

    Syst. Biol.

    (1996)
  • G. Jobb et al.

    TREEFINDER: a powerful graphical analysis environment for molecular phylogenetics

    BMC Evol. Biol.

    (2004)
  • D.S. Jordan

    The law of geminate species

    Am. Nat.

    (1908)
  • S. Kato et al.

    Field guide to eastern Pacific and Hawaiian sharks

    US Fish and Wildl. Serv., Circ.

    (1967)
  • M.X. Kirby et al.

    Lower Miocene stratigraphy along the Panama Canal and its bearing on the Central American Peninsula

    PLoS ONE

    (2008)
  • Cited by (0)

    View full text