Molecular phylogeny of Squaliformes and first occurrence of bioluminescence in sharks
© Straube et al. 2015
Received: 13 January 2015
Accepted: 4 August 2015
Published: 16 August 2015
Squaliform sharks represent approximately 27 % of extant shark diversity, comprising more than 130 species with a predominantly deep-dwelling lifestyle. Many Squaliform species are highly specialized, including some that are bioluminescent, a character that is reported exclusively from Squaliform sharks within Chondrichthyes. The interfamiliar relationships within the order are still not satisfactorily resolved. Herein we estimate the phylogenetic interrelationships of a generic level sampling of “squaloid” sharks and closely related taxa using aligned sequences derived from a targeted gene capture approach. The resulting phylogenetic estimate is further used to evaluate the age of first occurrence of bioluminescence in Squaliformes.
Our dataset comprised 172 putative ortholog exon sequences. Phylogenetic estimates result in a fully resolved tree supporting a monophyletic lineage of Squaliformes excluding Echinorhinus. Non-luminous Squalidae are inferred to be the sister to a clade comprising all remaining Squaliform families. Our results suggest that the origin of photophores is coincident with an elevated diversification rate and the splitting of families Dalatiidae, Etmopteridae, Oxynotidae and Somniosidae at the transition of the Lower to the Upper Cretaceous. The presence of luminous organs was confirmed for the Sleeper shark genus Zameus. These results indicate that bioluminescence in sharks is not restricted solely to the families Etmopteridae and Dalatiidae as previously believed.
The sister-clade to non-luminous Squalidae comprises five families. The presence of photophores is reported for extant members of three out of these five families based on results of this study, i.e. Lantern sharks (Etmopteridae), Kitefin sharks (Dalatiidae) and Sleeper sharks (Somniosidae). Our results suggest that the origin of luminous organs arose during the rapid diversification event that gave rise to the extant Squaliform families. These inferences are consistent with the idea of diversification of Squaliform sharks being associated with the emergence of new deep-sea habitats in the Lower Cretaceous, which may have been facilitated by the evolution of bioluminescence.
Squaliform sharks constitute a group of highly specialized species with a predominantly deep-dwelling lifestyle. They represent a substantial part of extant shark diversity (~27 % ) comprising 24 genera and more than 130 described species . Many Squaliform species are bioluminescent, a feature which appears to be exclusive within the Chondrichthyes. Currently, the families Echinorhinidae (Bramble - and Prickly sharks), Squalidae (Dogfish sharks), Centrophoridae (Gulper Sharks), Somniosidae (Sleeper sharks), Oxynotidae (Rough sharks), Dalatiidae (Kitefin sharks), and Etmopteridae (Lantern sharks) are discussed to form the Squaliformes. However, some previous morphological studies have suggested alternative intergeneric and interfamilial arrangements for the group [2–12].
The phylogenetic placement of Echinorhinidae has remained ambiguous in both morphological and molecular studies, either being included within Squaliformes, considered sister to Squaliformes, or placed in a separate group with Saw sharks (Pristiophoriformes) or Angel sharks (Squatiniformes). Further, recent molecular studies have recovered Squalidae, Centrophoridae, Dalatiidae, and Etmopteridae as monophyletic lineages within the Squaliformes, however, their interfamiliar relationships remain partially unresolved while the family Somniosidae appeared paraphyletic as Oxynotidae cluster within Somniosidae [2, 3, 8–24].
All of the molecular data sets examined to date have been based on the analysis of a single or few genes and none have recovered substantial support for branching events at the family level, likely due to limited phylogenetic signal supporting deeper nodes. Phylogenetic analyses based on morphological characters have not yielded consistent results either, e.g. [9, 10].
A dataset with strong phylogenetic signal is prerequisite for analyses of the evolution of taxa through time. So far, molecular clock analyses have delivered conflicting results concerning the origin and radiation ages of Squaliform sharks in general and the rise of families in particular [23, 24]. Molecular clocks are best calibrated using information from fossils or from vicariant biogeographic events. Squaliformes are well documented in the fossil record for sharks, which is largely comprised of teeth. Most Squaliform sharks display diagnostic clade specific dentitions pointing to high levels of trophic specialization and conservatism. A number of fossils can therefore be readily assigned to extant lineages such as the Gulper shark genus Centrophorus  or the Viper dogfish Trigonognathus , without the need to erect distinct genera for extinct forms whose phylogenetic affinities are unclear. According to , the fossil record of Squalidae extends back to the Upper Jurassic, while families Centrophoridae, Etmopteridae, Somniosidae, Oxynotidae, and Dalatiidae appeared rather instantaneously at the beginning of the Upper Cretaceous, which has been suggested to be a period of adaptive evolution in response to new ecological opportunities [23, 24]. The oldest Echinorhinid fossils are recorded from the Lower Cretaceous [25, 27] the evolution of bioluminescence in Kitefin (Dalatiidae) and Lantern sharks (Etmopteridae) appears to be correlated with the diversification of Squaliform sharks in the deep-sea [23, 24, 28, 29]. Surprisingly, it has not been clear at which point in their evolutionary trajectory, squaliform sharks first acquired photophores. Despite the fact that Shirai  had noted that all squaloid sharks except Echinorhinus, Centrophorus, Cirrhigaleus, Deania, Somniosus, and Squalus bear luminous organs, several recent studies suggested that photophores are only present in Etmopteridae and Dalatiidae [2, 23, 30, 31].
In this study, we estimate the phylogenetic interrelationships of Squaliform sharks by applying a gene capture approach that targets a large number of single-copy nuclear exons  to a generic level sampling of “squaloid” sharks and closely related taxa . We have used these data in conjunction with fossil calibration data, to estimate times of divergence and diversification rates among the extant lineages examined. We have also explored the potential role that bioluminescence may have had in promoting diversification in these animals, by reconstructing ancestral character states based on the inferred tree and the presence of photophores in extant forms.
Results and discussion
Molecular phylogeny of Squaliformes
On average, 200,000 of 352,605 possible basepairs, were sequenced per specimen (Additional file 1: Table S1). Characteristics of the raw dataset are given in Additional file 1: Table S2. Missing data were randomly distributed among specimens resulting in a large amount of incomplete sequences per captured locus and specimen.
MARE [33, 34] detected 174 phylogenetically informative loci in the raw dataset (Additional file 1: Figure S1). Re-blasting the full genome of C. milii against the 174 phylogenetically informative loci resulted in two potentially paraloguous loci (cds 1200 (unknown) and cds 1366 (LRP4)). Excluding these two loci and repeating the maximum likelihood analysis as described above did not affect the inferred tree topology.
This phylogenetic estimate reveals two major clades: the Squaliformes excluding Echinorhinidae and a clade containing Squatina, Pristiophoriformes, and Echinorhinus (Fig. 1). Within this clade, Echinorhinus is sister to Squatina and Pristiophoriformes. Results suggest that Echinorhinidae are not Squaliform sharks, but are the sister group to Angel- (Squatiniformes) and Saw sharks (Pristiophoriformes), as previously suggested by the analysis of mitochondrial data . Therefore, Squaliformes form a monophyletic group only, if Echinorhinus is excluded. This study does not support results from , suggesting Echinorhinus being the sistergroup to the remaining Squaliform lineages. The node time estimation for the Echinorhinus lineage suggests an Upper Jurassic splitting of the extant Echinorhinus lineage and the Squatina plus Pristiophoriformes clade. This dates the Echinorhinus lineage older than anticipated from the fossil record, which reports the oldest echinorhinid fossil from the early Cretaceous (Hauterivian) of southeastern France , while the oldest squatinids already appear in the Upper Jurassic .
Within the Squaliform clade, the first split separates Squalidae from the remaining families Centrophoridae, Etmopteridae, Dalatiidae, Somniosidae, and Oxynotidae. The genera Squalus and Cirrhigaleus appear as sister taxa. Centrophoridae split from Etmopteridae, Dalatiidae, Somniosidae, and Oxynotidae, where genera Deania and Centrophorus are sister. Dalatiidae are sister to a clade comprising Etmopteridae, Somniosidae, and Oxynotidae. There are two clades within the dalatiids, one comprising the Isistius and Dalatias lineages, the other Squaliolus and Euprotomicrus. As shown in Fig. 1, Somniosidae sensu stricto form two clearly distinct lineages that are sister to each other, one containing the genus Somniosus (Fig. 1), the other lineage contains all other remaining somniosid genera. Oxynotidae cluster within Somniosidae (Fig. 1). Within Etmopteridae, Trigonognathus is sister to a clade comprising Aculeola and Centroscyllium. Etmopterus is sister to this previously described clade forming four distinct lineages representing the subclades described in .
Oxynotus is inferred to be nested within Somniosidae, rendering the family Somniosidae paraphyletic (Fig. 1) in the current study. This result is repeatedly recovered in phylogenetic estimates based on DNA sequence data (both mitochondrial and nuclear) [19–24]. Given the consistency of the inferences from molecular data, it would be interesting to see if any anatomical features also support the link between Oxynotidae and Somniosidae. Oxynotus clusters with a group of otherwise morphologically similar species of somniosids, i.e. along with Zameus, Centroselachus, Scymnodon, and Centroscymnus. Our molecular results show that all five genera are closely related (Fig. 1). This is especially evident when comparing intergeneric diversity within Somniosidae with the large intrageneric sequence differences evident within the genus Etmopterus (Fig. 1). Moreover, there are limited morphological characters that can be used to differentiate some of these taxa [8, 39]. Together these results imply that assigning separate generic status to some species within Somniosidae may be an overrepresentation of the true diversity within the family.
Occurrence and significance of bioluminescence in Squaliform sharks
Node time estimates for major splitting events
95 % HPD
190 – 241.32
Middle Triassic to Lower Jurassic
Splitting of Squaliformes from the clade comprising Echinorhinus, Squatina, Pliotrema & Pristiophorus
153.85 – 203.99
Upper Triassic to Upper Jurassic
Clade comprising Echinorhinus, Squatina, Pliotrema & Pristiophorus
130 – 143.18
Split Centrophoridae from Dalatiidae, Etmopteridae, Oxynotidae & Somniosidae
113.94 – 137.88
Split Dalatiidae from Etmopteridae, Oxynotidae & Somniosidae
99.2 – 131.01
Transition Lower to Upper Cretaceous
Split Etmopteridae from Somniosidae & Oxynotidae
92.81 – 124.88
Transition Lower to Upper Cretaceous
Split Centrophorus from Deania
89 – 96.84
Split Somniosus from Oxynotidae & remaining Somniosidae
64.8 – 114.49
65 – 105.4
65 – 90.66
46.28 – 74.64
Upper Cretaceous to Palaeocene
Split Trigonognathus from Aculeola & Centroscyllium
44.5 – 76.86
Upper Cretaceous to Palaeocene
Radiation Somniosidae excluding Somniosus
24.46 – 63.94
Split Oxynotus from Scymnodon
15.32 – 47.11
We reconstructed ancestral character states in order to test the hypothesis that bioluminescence evolved in conjunction with the diversification of the Dalatiidae, Etmopteridae, Oxynotidae and Somniosidae. In the first analysis, we coded Dalatiidae and Etmopteridae as luminescent. Results from this analysis indicated that the common ancestor of families Dalatiidae, Etmopteridae, Oxynotidae, and Somniosidae was already likely carrying luminous organs. Interestingly, Somniosidae have been widely accepted as non-luminous [2, 23, 30, 31, 44]. However, Shirai  suggested that all Somniosidae are luminescent except for the genus Somniosus, which may have secondarily lost the ability to produce light.
Morphological data presented herein provide clear evidence that functional photophores are present within Somniosidae, at least within the genus Zameus (Fig. 3). All other inspected specimens showed no evidence of epidermal photophores. In light of this, the ancestral character state reconstruction was repeated incorporating results from the inspection of skin samples, i.e. coding the genus Zameus in addition to Etmopteridae and Dalatiidae as luminescent. Results from this analysis further increased the likelihood that the common ancestor of Dalatiidae, Etmopteridae and Somniosidae was luminescent (Fig. 2). The common ancestor of Centrophoridae, Etmopteridae, Dalatiidae, Oxynotidae, and Somniosidae is also implied to have been luminescent, but the likelihood is less compelling. A further analysis following  coding somniosid genera Centroselachus, Centroscymnus, Scymnodon, and Zameus as luminous further increases the likelihood so that the common ancestor of all Squaliformes except Squalidae may already have been luminescent (Additional file 1: Figure S8). This indicates that extant Centrophoridae may have secondarily lost their ability to emit light, i.e. that luminous organs may have already been present at the branching event giving rise to families Centrophoridae, Dalatiidae, Etmopteridae, Somniosidae, and Oxynotidae (Fig. 2). This suggests that luminescence evolved along and facilitated the Squaliform deep-sea radiation – a scenario that would be consistent with the elevated diversification rate detected for Etmopteridae, Somniosidae, and Oxynotidae. (Fig. 2, Additional file 1: Figures S8 and S9). We speculate that the common ancestor of families Dalatiidae, Etmopteridae, Oxynotidae, and Somniosidae was luminescent and used this to enhance camouflage by counterillumination as this is assumed to be the most basal function of shark bioluminescence [23, 28, 45, 47].
The occurrence of bioluminescence within the family Somniosidae is not surprising as especially the smaller sized genera (Centroselachus, Centroscymnus, Scymnodon, and Zameus) occur in sympatry with other luminous sharks such as etmopterids and dalatiids as well as a number of other luminescent deep-sea taxa including myctophid fishes which interestingly were estimated to have radiated in a similar time window . Results presented here lend further support to the hypothesis that bioluminescence in sharks evolved only once [29, 47]. Work in progress will allow identifying all luminous taxa within the family Somniosidae.
Our findings provide insights into the phylogeny of Squaliform sharks as well as the evolution of bioluminescence in the group. The radiation is estimated to have started in the Lower Cretaceous and continued through to the Upper Cretaceous. The initial elevated diversification rate is correlated with the likely first occurrence of luminous organs in sharks. The presence of photophores was confirmed for the genus Zameus in the family Somniosidae, implying that bioluminescence in sharks is not restricted to families Etmopteridae and Dalatiidae as is widely believed.
Targeted gene capturing
To ensure correct sample IDs of target samples, we either used genomic DNA of specimens which were previously analyzed in [21, 23, 52, 53] or generated NADH2 sequences as described in  and compared to the samples analysed in [21, 52, 53]. In the latter case, genomic DNA was extracted from collection material (tissues in GJPN tissue collection) already used in previous studies and stored in 95 % ethanol. Genomic DNA was obtained using the Promega Wizard ® DNA Purification System (Fisher Scientific). Total amounts of DNA were measured using a Qbit® Fluorometer (Life Technologies).
Subsequently, genomic DNA of the 28 target samples was sheared to approximately 500 bp using a Covaris® Sonicator. Sheared samples were used to prepare Illumina sequencing libraries following the protocol provided in . See Additional file 1: Table S1 for an overview of samples analysed.
We designed custom RNA bait libraries for targeting putatively single-copy orthologous genes based on sequences derived from seven shark species in , i.e. Chlamydoselachus anguineus, Etmopterus joungi, Isurus oxyrinchus, Orectolobus halei, Carcharhinus amblyrhynchos, Heterodontus portusjacksoni, and Squatina nebulosa. Each bait library comprised a pooled series of 120 bp baits designed for each target locus. As in , a 60 bp tiled overlap across baits was used to generate two-fold redundancy coverage for each target gene. When the length of the target gene was less than 120 bp, the sequence was extended in length to 120 bp by adding thymine nucleotides. The baits were manufactured by MyCroarray® (Ann Arbor, MI, USA).
Thereafter, gene capture was conducted by hybridization of target DNA to the baits. After hybridization, unbound and non-target DNA was washed away . The remaining library was enriched for target loci and was re-amplified to incorporate sample specific indices. Samples were pooled in equimolar ratios for sequencing. The pooled product was quantified using the CFX Connect Real-Time PCR system (Bio-Rad, Hercules, CA). Pooled sample was diluted to 8 pM and used for paired-end 150 bp or 250 bp sequencing on an Illumina MiSeq sequencing instrument (Illumina, Inc, San Diego, CA). Sequence reads associated with each sample were identified by their respective indices.
Alignment reconstruction of gene capture data
Adapters were trimmed from sequence reads using Trimgalore v0.3.7. [54, 55] and assembled de novo using ABySS ver1.3.5.  with a k-mer of 64. Assembled contigs were assigned to core ortholog groups using HaMStR . The core ortholog database consisted of profile hidden Markov models of orthologous sequence groups from model vertebrates . Any sequence that matched a core-ortholog pHMM was provisionally assigned to the corresponding orthologous group. In order to be retained in the final matrix, provisional sequence hits also had to satisfy a reciprocal best BLAST criterion when compared to Callorhinchus milii as the reference taxon. Orthologous exons were trimmed from non-target intron information and aligned with Mafft [58, 59]. Finally, all loci were concatenated.
Data analysis and phylogenetic reconstruction
Maximum likelihood (ML) trees were estimated using RAxML GUI [33, 60]. The initial ML analysis used the complete concatenated dataset (i.e. 1265 loci and 28 taxa) under GTR GAMMA using different partitioning schemes (Additional file 1: Table S4) using the automatic halt for bootstrapping . Squalomorph sharks are widely accepted as monophyletic [8, 19–22]. Within Squalomorphs, Hexanchiformes are considered to form the most basal lineage [8, 20–22], therefore, Hexanchus griseus was chosen as the outgroup taxon.
Subsequently, MARE ver.1.2 [33, 34] was used to examine the dataset for phylogenetically informative sites and taxa. MARE [33, 34] was designed to identify the most phylogenetically informative subset of sites contained in phylogenomic data sets. It is especially well-suited to analysis of data sets with a high proportion of missing data. MARE [33, 34] identified 174 maximally informative loci for our data set reducing the maximal total sequence length per specimen from 352,605 bp to 73,925 bp (Additional file 1: Figure S1).
As a further scan for paraloguous sequences, we re-blasted the full genome of Callorhinchus milii against the remaining 174 loci to check, if each sequence has only a single hit in the C. milii genome using customized Perl scripts (Additional file 1).
The reduced nuclear dataset comprising 172 concatenated nucleotide loci is deposited at . This data was analysed as described for the full dataset and additionally RaxML GUI [35, 60] was applied to the amino acid alignment comprising infomative loci using the best partitioning scheme suggested by PartitionFinder Protein v1.1.1 [62, 63]. Further, the reduced 172 loci DNA sequence data alignment was analysed with PartitionFinder v1.1.1 to look for best fitting partition schemes and models of molecular evolution [62, 63] to determine if different partition types influence the tree topology. We used the rcluster option with a rcluster percentage of 10  for the analysis. The Bayesian mixture model CAT  implemented in PhyloBayes 3.3f [36, 37] was used on the concatenated 172 amino acid alignment to partition sites into different rate categories using non-parametric modeling of site specific effects. This allowed us a topological comparison to the ML analysis under GTR GAMMA . Four independent chains were run in parallel. The tracefiles and treelists of all four chains were used to check for convergence . The analysis was stopped with a maximum difference of 0.16 and effective sample sizes exceeding 100, with the exception of the allocent statistic (see Additional file 1). A majority rule consensus tree was computed from 12997 input trees from each chain with a burn-in of 1000 trees and analyzing every second tree of the pooled trees. The consensus tree was rooted midpoint.
See Additional file 1: Table S4 and Figures S3 to S7 for a summary of partitioning schemes and phylogenetic analyses conducted. To ensure that the choice of a single outgroup does not have a negative effect such as long branch attraction on our phylogenetic analysis, we performed analysis not defining an outgroup, defining different outgroups as well as deleting Hexanchus griseus from the dataset and re-computing a phylogenetic estimate without defining an outgroup taxon (Additional file 1: Figure S2).
Node time estimation and diversification rate
Calibration points used for node time estimates of squaloid sharks
Minimum age (Ma)
Soft upper bound (Ma)
Root age (Squalomorphii)
Echinorhinidae, Squatinidae & Pristiophoridae
Trigonognathus, Aculeola & Centroscyllium
The clade comprising Echinorhinidae, Pristiophoriformes Pliotrema and Pristiophorus as well as Squatina was assumed to vary in age between 145 to 163 Ma (Upper Jurassic) based on articulated fossils of Squatinids at the lower and Echinorhinus sp. teeth at the upper end of the time frame [25, 27]. The minimum age of Squaliformes was calibrated to 130 Ma based on the fossil taxon Protosqualus with a soft upper bound at 163 Ma allowing for the possibility that Squaliform sharks were already present in the Upper Jurassic . Squaliform family-level diversity is assumed to have originated in the late Mesozoic (Upper Cretaceous) , while most extant genera likely originated in the Cenozoic. Fossil evidence was used to calibrate the minimum age of Centrophoridae, Etmopteridae and Dalatiidae to be 65 Ma (C/T boundary) with a soft upper bound of 100 Ma (beginning of the Upper Cretaceous). Further, the clade comprising Trigonognathus, Aculeola, and Centroscyllium was assumed to be of minimum age of 45 Ma and a lower bound of 100 Ma based on the fossil record of Trigonognathus virginiae  and the age estimate of extant Etmopteridae in .
All analyses assumed an exponential prior distribution for calibration points. Three independent runs were performed with a Markov Chain lasting 90 million generations each, sampling trees every 1000 generations. One run included the maximum likelihood inferred tree with the highest likelihood as a newick formatted starting tree. This BEAST input file is deposited in the Dryad data repository . Combined log files were analyzed in Tracer v.1.6  to check, if the effective sample sizes (ESS) of parameters represent the posterior distribution adequately; further trace and density plots were checked for convergence of the MCMC and posterior probability distributions in different runs. After defining a burn-in of 25 % of all sampled trees in each run, TreeAnnotator  was used to create a consensus tree which was visualized in FigTree v.1.4.0 .
We used the R  module MEDUSA (modeling evolutionary diversification using stepwise AIC)  implemented in the GEIGER package  to estimate changes in the diversification rate based on the consensus chronogram attained from the BEAST  analysis. Species richness values were obtained from .
Ancestral bioluminescence within Squaliformes
Ancestral character states of bioluminescence were reconstructed using Maximum Likelihood estimates implemented in the R  package GEIGER  and are based on the chronogram attained from the BEAST  analyses. In a first analysis, we coded only Dalatiidae and Etmopteridae as luminescent. Results from this analysis indicated that the common ancestor of families Dalatiidae, Etmopteridae, and Somniosidae was already likely luminescent. As an empirical test of this idea, we subsequently inspected the ventral surface area of Somniosidae and Oxynotidae specimens from the Bavarian State Collection of Zoology –Centroselachus crepidater (ZSM30842), Centroscymnus owstonii (ZSM36725), Oxynotus bruniensis (ZSM30862), and Zameus squamulosus (ZSM30966)– and the Zoological Museum Hamburg –Somniosus microcephalus (ZMH 123507), S. rostratus (ZMH 25751), Centroscymnus coelolepis (ZMH 119748), Centroscymnus owstonii (ZMH 104894), Centroselachus crepidater (ZMH 103185), Scymnodalatias sp. (ZMH 122774), Zameus squamulosus (ZMH 120262; ZMH 120485). When pigmentation was apparent, a 1 cm2 skin patch was excised from the ventral surface of the specimen and observed under a binocular microscope (Leica MZ6, Wetzlar, Germany). If photophores were observed, a picture was taken and analysed in Image J v. 1.46 using a random 1 × 1 mm counting frame to estimate photophore mean diameter, photophore density and proportion of the skin surface area occupied by photophores (PAP) following the method of . Thereafter, the ancestral character state reconstruction was repeated incorporating results from the inspected skin samples in a second analysis, and incorporating information on presence of luminous organs in somniosids following  in a third test. See Additional file 1 for documentation on R scripts used and the different photophore presence/ absence matrix (Additional file 1: Table S6).
Availability of supporting data
The data sets supporting the results of this article are available in the Dryad repository, http://datadryad.org/review?doi=doi:10.5061/dryad.n3581. See also Additional file 1.
This project was funded by the National Science Foundation (NSF), grant “Jaws and Backbone: Chondrichthyan Phylogeny and a Spine for the Vertebrate Tree of Life”; DEB-01132229 to GJPN. JMC is a postdoctoral researcher at Fonds National de la Recherche Scientifique (FNRS, Belgium). We would like to express our sincere thanks to Ting Kuang (Shanghai Ocean University, Shanghai) for her help with re-blasting, Elisabeth Rochel (CofC, Charleston) for help in the lab, Thomas Fussel (CofC, Charleston) and Adam Bazinet (University of Maryland, College Park) for analysis pipeline programming. Alan Pradel (MNHN, Paris), Jürgen Pollerspöck, Ulrich K Schliewen (ZSM, Munich), Alexander Cerwenka (ZSM, Munich), and Frederic Schedel (ZSM, Munich) are thanked for fruitful discussions. We would further like to acknowledge institutional support at the Bavarian State Collection of Zoology (ZSM, Munich, Dirk Neumann) as well as the Zoological Museum Hamburg (ZMH, Hamburg, Simon Weigmann).Two anonymous reviewers are acknowledged for their constructive criticism.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Pollerspöck J, Straube N: www.shark-references.com, World Wide Web electronic publication, 2015
- Ebert DA, Fowler S, Compagno LJV. Sharks of the World – A Fully Illustrated Guide. Plymouth: Wild Nature Press; 2013.Google Scholar
- Compagno LJV. Interrelationships of Living Elasmobranchs. In: Greenwood PH, editor. Interrelationships of Fishes. London: Academic; 1973. p. 15–61.Google Scholar
- Compagno LJV. Phyletic relationships of living sharks and rays. Am Zool. 1977;17:303–22.Google Scholar
- Maisey JG. An evaluation of jaw suspension in sharks. Am Mus Novi. 1980;2706:1–17.Google Scholar
- Maisey JG. Chondrichthyan phylogeny: a look at the evidence. J Vertebr Paleontol. 1984;4(3):359–71. doi:10.1080/02724634.1984.10012015.View ArticleGoogle Scholar
- Maisey JG. Higher elasmobranch phylogeny and biostratigraphy. Zool J Linnean Soc. 1984;82(1–2):33–54. doi:10.1111/j.1096-3642.1984.tb00534.x.View ArticleGoogle Scholar
- Shirai S. Squalean Phylogeny and Related Taxa. Hokkaido: University Press; 1992.Google Scholar
- Shirai S. Phylogenetic Interrelationships of Neoselachians (Chondrichthyes, Euselachii). In: Stiassny MLJ, Parenti LR, Johnson GD, editors. Interrelationships of Fishes. San Diego: Academic; 1996. p. 9–34.View ArticleGoogle Scholar
- Adnet S, Cappetta H. A palaeontological and phylogenetical analysis of squaliform sharks (Chondrichthyes: Squaliformes) based on dental characters. Lethaia. 2001;34:234–48.View ArticleGoogle Scholar
- de Carvalho MR. Higher-Level Elasmobranch Phylogeny, Basal Squaleans, and Paraphyly. In: Stiassny MLJ, Parenti LR, Johnson GD, editors. Interrelationships of Fishes. San Diego: Academic; 1996. p. 35–62.View ArticleGoogle Scholar
- de Carvalho MR, Maisey JG. The Phylogenetic Relationship of the Late Jurassic Shark Protospinax Woodward 1919 (Chondrichthyes: Elasmobranchii). In: Arratia G, Viohl G, editors. Mesozoic Fishes – Systematics and Paleoecology. München: Friedrich Pfeil Verlag; 1996. p. 9–46.Google Scholar
- Bass AJ, Compagno LJV, Heemstra PC. Squalidae. In: Smith MM, Heemstra PC, editors. Smith’s Sea Fishes. Berlin: Springer Verlag; 1986. p. 49–62.Google Scholar
- Bass AJ, Compagno LJV. Families Echinorhinidae, Proscylliidae, Scyliorhinidae. In: Smith MM, Heemstra PC, editors. Smith’s Sea Fishes. 63rd ed. Johannesburg: Macmillian; 1986. p. 87–95.Google Scholar
- Compagno LJV. FAO Species Catalogue. Vol. 4. Sharks of the World: An Annotated and Illustrated Catalogue of Shark Species Known to Date. Part 2. Carcharhiniformes. Rome: U.N. Food and Agriculture Organization; 1984. p. 251–655.Google Scholar
- Garrick JAF. Studies on New Zealand Elasmobranchii. Part VII. The Identity of Specimens of Centrophorus from New Zealand. Transactions of the Royal Society of New Zealand. 1959;86(1):127–41.Google Scholar
- Herman J. Die Selachier-Zähne aus der Maastricht Stufe von Hemmoor, Niederelbe (NW-Deutschland). Geologisches Jahrbuch Reihe A. 1982;61:29–159.Google Scholar
- Pfeil FH. Zahnmorphologische Untersuchungen an Rezenten und Fossilen Haien der Ordnungen Chlamydoselachiformes und Echinorhiniformes. Palaeo Ichthyologica. 1983;1:1–315.Google Scholar
- Maisey JG, Naylor GJP, Ward DJ. Mesozoic Elasmobranchs, Neoselachian Phylogeny and the Rise of Modern Elasmobranch Diversity. In: Arratia G, Tintori A, editors. Mesozoic Fishes 3-Systematics. Palaeoenvironments and Biodiversity. München: Verlag Dr. Friedrich Pfeil; 2004. p. 17–56.Google Scholar
- Velez-Zuazo X, Agnarsson I. Shark tales: a molecular species-level phylogeny of sharks (Selachimorpha, Chondrichthyes). Mol Phylogenet Evol. 2011;58:207–17.View ArticlePubMedGoogle Scholar
- Naylor GJP, Caira JN, Jensen K, Rosana KAM, Straube N, Lakner C. Elasmobranch Phylogeny: A Mitochondrial Estimate Based on 595 Species. In: Carrier JC, Musick JA, Heithaus MR, editors. Biology of Sharks and Their Relatives. Boca Raton: CRC Press; 2012. p. 31–56.View ArticleGoogle Scholar
- Naylor GJP, Ryburn JA, Ferigo O, Lopez A. Phylogenetic Relationships Among the Major Lineages of Modern Elasmobranchs. In: Hamlett WC, editor. Reproductive Biology and Phylogeny of Chondrichthyes: Sharks, Batoids and Chimaeras. Enfield: Science Publishers; 2005. p. 1–25.Google Scholar
- Straube N, Iglésias SP, Sellos DY, Kriwet J, Schliewen UK. Molecular Phylogeny and Node Time Estimation of Bioluminescent Lanternsharks (Elasmobranchii: Etmopteridae). Mol Phylogenet Evol. 2010;56:905–17.View ArticlePubMedGoogle Scholar
- Sorenson L, Santini F, Alfaro ME. The Effect of Habitat on Modern Shark Diversification. J Evol Biol. 2014;27(8):1536–48. doi:10.1111/jeb.12405.View ArticlePubMedGoogle Scholar
- Maisey JG. What is an ‘elasmobranch’? The impact of palaeontology in understanding elasmobranch phylogeny and evolution. J Fish Biol. 2012;80(5):918–51. doi:10.1111/j.1095-8649.2012.03245.x.View ArticlePubMedGoogle Scholar
- Cappetta H, Adnet S. Discovery of the recent genus Trigonognathus (Squaliformes: Etmopteridae) in the Lutetian of Landes (southwestern France). Remarks on the teeth of the recent species Trigonognathus kabeyai. Paläontol Z. 2001;74(4):575–81.View ArticleGoogle Scholar
- Adnet S, Guinot G, Cappetta H, Welcomme JL. Oldest evidence of bramble sharks (Elasmobranchii, Echonorhinidae) in the Lower Cretaceous of southeast France and the evolutionary history of orbitostylic sharks. Cretac Res. 2012;35:81–7.View ArticleGoogle Scholar
- Claes JM, Nilsson DE, Straube N, Collin SP, Mallefet J. Iso-luminance counterillumination drove bioluminescent shark radiation. Sci Rep. 2014;4:4328.PubMed CentralView ArticlePubMedGoogle Scholar
- Klug S, Kriwet J. Timing of deep-sea adaptation in dogfish sharks: insights from a supertree of extinct and extant taxa. Zool Scr. 2009;39:331–42. doi:10.1111/j.1463-6409.2010.00427.x.View ArticleGoogle Scholar
- Compagno LJV, Dando M, Fowler S. A field guide to the sharks of the world. London; Collins; 2005.Google Scholar
- Claes JM, Mallefet J. Bioluminescence of Sharks: First Synthesis. In: Meyer-Rochow VB, editor. Bioluminescence in Focus - A Collection of Illuminating Essays. Kerala: Research Signpost; 2009. p. 51–65.Google Scholar
- Li C, Hofreiter M, Straube N, Corrigan S, Naylor GJP. Capturing proteincoding genes across highly divergent species. Biotechniques. 2013;54:321–6. doi:10.2144/000114039.PubMedGoogle Scholar
- Meyer B, Meusemann K, Misof B: MARE. MAtrix REduction - A tool to select optimized data subsets from supermatrices for phylogenetic inference. [https://www.zfmk.de/de/forschung/forschungszentren-und-gruppen/mare]
- Misof B, Meyer B, von Reumont BM, Kück P, Misof K, Meusemann K. Selecting informative subsets of sparse supermatrices increases the chance to find correct trees. BMC Bioinformatics. 2013;14:348.PubMed CentralView ArticlePubMedGoogle Scholar
- Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3. open access.PubMed CentralView ArticlePubMedGoogle Scholar
- Lartillot N. PhyloBayes home page. [www.phylobayes.org].
- Lartillot N, Lepage T, Blanquart S. PhyloBayes 3: a Bayesian software package for phylogenetic reconstruction and molecular dating. Bioinformatics. 2009;25(17):2286–8. doi:10.1093/bioinformatics/btp368.View ArticlePubMedGoogle Scholar
- Straube N, Li C, Claes JM, Corrigan S, Naylor GJP (2015). Data from: Molecular Phylogeny of Squaliformes and first Occurrence of Bioluminescence in Sharks. http://datadryad.org/review?doi=doi:10.5061/dryad.n3581
- White WT, Vaz DFB, Ho H-C, Ebert DA, de Carvalho MR, Corrigan S, et al. Redescription of Scymnodon ichiharai Yano and Tanaka 1984 (Squaliformes: Somniosidae) from the western North Pacific, with comments on the definition of somniosid genera. Ichthyol Res. 2014. doi:10.1007/s10228-014-0430-y.Google Scholar
- Drummond AJ, Ho S, Philipps M, Rambaut A. Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006;4:88.View ArticleGoogle Scholar
- Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214.PubMed CentralView ArticlePubMedGoogle Scholar
- Harmon LJ, Weir JT, Brock CD, Glor RE, Challenger W. GEIGER: Investigating evolutionary radiations. Bioinformatics. 2009;24:129–31.View ArticleGoogle Scholar
- Alfaro ME, Santini F, Brock C, Alamillo H, Dornburg A, Rabosky DL, et al. Nine exceptional radiations plus high turnover explain species diversity in jawed vertebrates. PNAS. 2009;106(32):13410–4.PubMed CentralView ArticlePubMedGoogle Scholar
- Claes JM, Krönström J, Holmgren S, Mallefet J. GABA inhibition of luminescence from lantern shark (Etmopterus spinax) photophores. Comp Biochem Physiol C. 2011;153:231–6.Google Scholar
- Claes JM, Mallefet J. The lantern shark’s light switch: turning shallow water crypsis into midwater camouflage. Biol Lett. 2010;6:685–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Claes JM, Sato K, Mallefet J. Morphology and control of photogenic structures in a rare dwarf pelagic lantern shark (Etmopterus splendidus). J Exp Mar Biol Ecol. 2011;406:1–5.View ArticleGoogle Scholar
- Claes JM, Ho H-C, Mallefet J. Control of luminescence from pygmy shark (Squaliolus aliae) photophores. J Exp Biol. 2012;215(10):1691–9. doi:10.1242/jeb.066704.View ArticlePubMedGoogle Scholar
- Renwart M, Delroisse J, Claes JM, Mallefet J. Cytological changes during luminescence production in lanternshark (Etmopterus spinax Linnaeus, 1758) photophores. Zoomorphology. 2015;134:107–16.View ArticleGoogle Scholar
- Herring PJ, Morin JG. Bioluminescence of Fishes. In: Herring PJ, editor. Bioluminescence in Action. New York: Academic; 1978. p. 273–329.Google Scholar
- Claes JM, Dean MN, Nilsson DE, Hart NS, Mallefet J. A deepwater fish with ‘lightsabers’ – dorsal spine-associated luminescence in a counterilluminating lanternshark. Sci Rep. 2013;3:1308.PubMed CentralView ArticlePubMedGoogle Scholar
- Davis MP, Holcroft NI, Wiley EO, Sparks JS, Smith WL. Species specific bioluminescence facilitates speciation in the deep Sea. Mar Biol. doi: 10.1007/s00227-014-2406-x
- 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;367:1–262.View ArticleGoogle Scholar
- Straube N, White WT, Ho H-C, Rochel E, Corrigan S, Li C, Naylor GJP: A DNA sequence based identification checklist for Taiwanese chondrichthyans. In: Systematics and biodiversity of sharks, rays, and chimaeras (Chondrichthyes) of Taiwan. Edited by Carvalho MR de, Ebert DA, Ho H-C, White WT. Zootaxa 3752: 256–278Google Scholar
- Trimgalore [www.bioinformatics.babraham.ac.uk/projects/trim_galore]
- Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet journal. 2011;17:10–2.View ArticleGoogle Scholar
- Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJ, Birol I. ABySS: a parallel assembler for short read sequence data. Genome Res. 2009;19:1117–23.PubMed CentralView ArticlePubMedGoogle Scholar
- Ebersberger I, Strauss S, von Haeseler A. HaMStR: profile hidden markov model based search for orthologs in ESTs. BMC Evol Biol. 2009;9:157.PubMed CentralView ArticlePubMedGoogle Scholar
- Katoh K, Misawa K, Kuma K, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nuc Ac Res. 2002;30(14):3059–66. doi:10.1093/nar/gkf436. PMC 135756.View ArticleGoogle Scholar
- Katoh K, Kuma K, Toh H, Miyata T. MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nuc Ac Res. 2005;33(2):511–8. doi:10.1093/nar/gki198. PMC 548345.View ArticleGoogle Scholar
- Silvestro M. RAxMLGUI: a graphical front-end for RAxML. Org Divers Evol. 2012;12:335–7.View ArticleGoogle Scholar
- Pattengale ND, Alipour M, Bininda-Edmonds ORP, Moret BME, Stamatakis A. How many bootstrap replicates are necessary? J Comp Biol. 2010;17:337–54.View ArticleGoogle Scholar
- Lanfear R, Calcott B, Ho SYW, Guindon S. PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012;29(6):1695–701.View ArticlePubMedGoogle Scholar
- Lanfear R, Calcott B, Kainer D, Mayer C, Stamatakis A. Selecting optimal partitioning schemes for phylogenomic datasets: a comparison of clustering methods. BMC Evol Biol 2014, 14:82 doi:10.1186/1471-2148-14-82
- Lartillot N, Philippe H. A Bayesian mixture model for acrosssite heterogeneities in the amino-acid replacement process. Mol Biol Evol. 2004;21:1095–109.View ArticlePubMedGoogle Scholar
- Lartillot N, Brinkmann H, Hervé P. Suppression of long-branch attraction artefacts in the animal phylogeny using a site-heterogeneous model. BMC Evol Biol. 2007;7 Suppl 1:S4. doi:10.1186/1471-2148-7-S1-S4.PubMed CentralView ArticlePubMedGoogle Scholar
- Pradel A, Tafforeau P, Maisey JG, Janvier P. A New Paleozoic Symmoriiformes (Chondrichthyes) from the Late Carboniferous of Kansas (USA) and Cladistic Analysis of Early Chondrichthyans. PLoS ONE. 2011;6(9), e24938. doi:10.1371/journal.pone.0024938.PubMed CentralView ArticlePubMedGoogle Scholar
- Sansom IJ, Davies NS, Coates MI, Nicoll RS, Ritchie A. Chondrichthyan-like scales from the middle Ordovician of Australia. Palaeontology. 2012;55:243–7.View ArticleGoogle Scholar
- Tracer v.1.6 [http://tree.bio.ed.ac.uk/software/tracer/]
- Figtree v.1.6 [http://tree.bio.ed.ac.uk/software/figtree/]
- R Core Team: R: A Language and Environment for Statistical Computing [http://www.r-project.org]
- Ginter M, Hampe O, Duffin C. Handbook of Paleoichthyology, Vol. 3D, Chondrichthyes. Paleozoic Elasmobranchii: Teeth. Verlag Friedrich Pfeil: Munich; 2010.Google Scholar
- Cappetta H. Handbook of Paleoichthyology, Vol. 3B, Chondrichthyes II. Mesozoic and Cenozoic Elasmobranchii: teeth. Verlag Friedrich Pfeil: Munich; 2004.Google Scholar