Evidence of positive selection associated with placental loss in tiger sharks
© The Author(s). 2016
Received: 16 March 2016
Accepted: 2 June 2016
Published: 14 June 2016
All vertebrates initially feed their offspring using yolk reserves. In some live-bearing species these yolk reserves may be supplemented with extra nutrition via a placenta. Sharks belonging to the Carcharhinidae family are all live-bearing, and with the exception of the tiger shark (Galeocerdo cuvier), develop placental connections after exhausting yolk reserves. Phylogenetic relationships suggest the lack of placenta in tiger sharks is due to secondary loss. This represents a dramatic shift in reproductive strategy, and is likely to have left a molecular footprint of positive selection within the genome.
We sequenced the transcriptome of the tiger shark and eight other live-bearing shark species. From this data we constructed a time-calibrated phylogenetic tree estimating the tiger shark lineage diverged from the placental carcharhinids approximately 94 million years ago. Along the tiger shark lineage, we identified five genes exhibiting a signature of positive selection. Four of these genes have functions likely associated with brain development (YWHAE and ARL6IP5) and sexual reproduction (VAMP4 and TCTEX1D2).
Our results indicate the loss of placenta in tiger sharks may be associated with subsequent adaptive changes in brain development and sperm production.
KeywordsReproduction Viviparous Placenta Transcriptome RNA-Seq Positive selection Elasmobranchs Carcharhinids Galeocerdo
Aristotle was the first to record some animals give birth to live young (i.e. viviparity) whilst others lay eggs (i.e. oviparity) . Viviparity may offer selective advantages to parents and offspring, such as enhanced survival of offspring, compensation for low fecundity, amplification of reproductive niches to reduce competition, exploitation of pelagic niches, colonisation of new habitats, and possibly increased energetic efficiency; disadvantages may include reduced fecundity, cost to the female, and risk of brood loss through maternal death . Viviparity is thought to have first evolved over 350 million years ago (mya), and is an unprecedented example of convergent evolution having independently evolved at least 150 times in mammals, reptiles, amphibians, and fishes [3, 4]. Among these viviparous organisms there are differences in the method of foetal nutrition, with the supply of nutrients via yolk in eggs or yolk-sacs considered to be ancestral . In fact, yolk-sac placentation is the most common type of placentation in vertebrates .
Oviparity also has its own benefits given it has been retained by the majority of vertebrates. Although viviparity has evolved many times, it creates post-fertilization opportunities for genomic conflicts absent in oviparous species [6, 7]. Furthermore, transitions in reproductive modes require numerous morphological, physiological, and behavioural adaptations. These can be associated with variation in geographic distribution and environmental conditions, as documented in reptiles [8, 9] and amphibians .
Sharks are considered the most enduring of fishes having survived the mass extinction events of the last 420 million years. Approximately 70 % of these species give birth to live young and approximately 30 % lay eggs [11, 12]. Sharks are among the first vertebrates to display viviparity . They are also members of a lineage which is a sister group to all teleosts and tetrapods, thereby ideally placed for comparative studies of these taxa. To date, however, there have been few studies of sharks at genome level [13–15]; one of the reasons for this is the large size of shark genomes (up to 34 picograms per haploid genome; >30 Gb)  in comparison to bony fish models (e.g. zebrafish 1.4 Gb ).
During the early stages of development all viviparous sharks are nourished by a foetal yolk-sac . Viviparous sharks also display a diversity of embryonic nourishment derived from maternally obtained nutrients, known as matrotrophic nutrition. Types of matrotrophic nutrition in sharks include mucus produced by the uterus (mucoid histotrophy), supply of unfertilised eggs (oophagy), and direct exchange between maternal and foetal tissues via a placenta (placental viviparity). Placental connections develop in sharks when the empty yolk-sac morphs and attaches to the uterus wall . In viviparous shark species that do not develop placental connections (i.e. non-placental), the empty yolk-sac is reabsorbed into the developing embryo .
Placenta are thought to have one evolutionary origin in sharks and are restricted to five families from the Carcharhiniformes order (i.e., Carcharhinidae, Sphyrnidae, Hemigaleidae, Leptochariidae, and Triakidae) . Two of these families, Carcharhinidae and Triakidae, contain both placental and non-placental species, with the lack of placenta thought to be due to secondary loss . Additional nutrition via a placenta is suggested to increase the embryonic development rate of energetically expensive tissues such as the brain, and there is evidence associating increased reproductive nutrition in sharks with larger brain sizes, relative to body mass .
Here we focused on one shark family, the carcharhinids (Order: Carcharhiniformes, Family: Carcharhinidae), which all display placental viviparity, apart from the tiger shark (Galeocerdo cuvier) . As the sister groups to Galeocerdo are placental, it is most likely this lineage has undergone placental loss. Recent research has identified distinct embryonic nutrition displayed by tiger sharks where the egg case housing the embryo becomes filled with approximately one litre of energy-rich, yellowish fluid .
In this study, we reconstructed a phylogenetic relationship for the tiger shark, six other carcharhinid sharks, and two outgroup species. We used fossil data to produce a time-calibrated phylogeny to estimate when the tiger shark lineage diverged from the placental carcharhinids. We also aimed to identify orthologous genes among the transcriptomes of the nine species and determine if there is evidence of genes evolving under positive selection, and possibly associated with the loss of placenta in the tiger shark lineage. Based on the proposed link between embryonic nutrition and brain size, we hypothesised to detect in the tiger shark lineage evidence of positive selection in genes encoding proteins associated with brain development, as well as sexual reproduction.
Transcriptome statistics for the nine viviparous shark species sampled here
Number of filtered reads
Number of transcripts
Atlantic sharpnose shark (Rhizoprionodon terraenovae)
Blacknose shark (Carcharhinus acronotus)
Blue shark (Prionace glauca)
Bull shark (Carcharhinus leucas)
Caribbean reef shark (Carcharhinus perezii)
Dusky smoothhound (Mustelus canis insularis).
Lemon shark (Negaprion brevirostris)
Sand tiger shark (Carcharias taurus)
Tiger shark (Galeocerdo cuvier)
Transcriptome assembly and orthologue identification
We used Trinity (version 2013-05-08) to assemble species specific transcriptomes . The mean number of transcripts was 109,709 (Standard Deviation: 31,801) and the mean N50 value was 1835 (Standard Deviation: 348). We then identified from the open reading frames (ORFs) of assembled transcripts 3,215 putative orthologous sequences using a reciprocal best-hit Blast search. High confidence alignments were generated for these orthologous sequences using multiple aligners. Alignments were subsequently filtered again to remove low confidence codon alignments, finally resulting in 1,197 orthologues for further analysis.
Phylogenetic tree construction and positive selection analyses
The 1,197 orthologue alignments were concatenated (1,101,288 bp) and used to construct a phylogenetic tree (Additional file 1: Figure S1) using RAxML (version 8.0.0) . The 1,197 orthologues were then analysed with CodeML, a program of PAML (version 4.7) . The CodeML one ratio model (M0) found no orthologue exhibiting signatures of positive selection across the entire sequence. Therefore, we used an additional test, a comparison between the neutral model 7 (M7) and the non-neutral model 8 (M8) to identify specific regions of genes that may be evolving under positive selection. Sequences with significant M8:M7 likelihood ratios provide evidence of positive selection; therefore these sequences were further analysed using the Bayes Empirical Bayes (BEB) method to identify specific codons which may be under positive selection . Of the 1,197 orthologues across nine species analysed by CodeML and BEB, 95 orthologues (Additional file 2: Table S1) were found to have specific codon sites showing signatures of positive selection (Benjamini-Hochberg corrected p value < 0.05).
Fossil calibration used for calibrating the shark phylogeny. Minimum and upper bound fossils with estimated ages and references for nodes 1, 2, 3, 5 and 8
Fossil (minimum age)
Lower Bound Fossil Age (mya)
Fossil (maximum age)
Upper Bound Fossil Age (mya)
Lamniformes – Carcharhiniformes
Underwood & Ward, 2004 
Delsate & Lepage, 1990 
Triakidae – Carcharhinidae
Underwood et al. 1999 
Underwood & Ward, 2004 
Galeocerdo – Carcharhinidae
Noubhani & Cappetta, 1997 
Underwood et al. 1999 
Negaprion – Carcharhinus/Prionace (Node 5)
Probst, 1879 
Noubhani & Cappetta, 1997 
Prionace – Carcharhinus perezii (Node 8)
Landini, 1977 
Probst, 1879 
Using 1,007,817 bp of transcriptome data per species we reconstructed phylogenetic relationships for the species sampled. We included the dusky smoothhound and sand tiger shark as outgroups; the placental dusky smoothhound belongs to the only other shark family with both placental and non-placental species (Triakidae), whereas the sand tiger shark is one of many non-placental species making up a sister order to the Carcharhiniformes. The Carcharhiniformes and Lamniformes appear to have diverged before the evolution of placenta in sharks.
Inferred divergence times for the Carcharhinidae – Galeocerdo node in shark phylogenetic trees. Mean age and 95 % highest posterior density credibility intervals (HPD CI) are provided for this study and Sorenson et al. 
Sorenson et al. 2014 
Mean Age (mya)
95 % HPD CI (mya)
Mean Age (mya)
95 % HPD CI (mya)
Carcharhinidae – Galeocerdo
Genes evolving under positive selection in the tiger shark
Positive selection on genes linked to sexual reproduction
Tiger sharks have a gestation period of approximately 15-16 months, several months longer than many placental carcharhinids, and follow a triennial reproductive cycle, with an estimated one year of sexual inactivity . As male tiger sharks generally reproduce annually, changes to female reproductive behaviour may have led to an increase in sperm competition and post-copulatory sexual selection.
Here, we detected the signature of positive selection in the TCTEX1D2 and VAMP4 genes. TCTEX1D2 encodes a dynein-2 light chain protein required for cilia function and found in the flagellum of sperm in a variety of taxa, including humans, mice, teleosts and sea urchins [38–40]. Sperm motility is reduced and spermatogenesis disrupted in mice lacking TCTEX1D2 as a result of increased apoptosis in male germ cells . Deletion of VAMP4 causes a significant increase in sperm head abnormalities in mice, resulting from aberrant acrosome formation . This has implications for sperm morphology and hydrodynamics, therefore positive selection in TCTEX1D2 and VAMP4 may alter sperm count and motility, which could be a response to changes in sperm competition in tiger sharks induced by the loss of placenta.
It is assumed when there is strong sperm competition an individual can gain a competitive advantage by increasing sperm production ; however, there is a trade-off between sperm size and number . Sharks display extensive variation in sperm morphology which is thought to demonstrate variation in the intensity of post-copulatory sexual selection [44–46]. Tanaka et al.  showed tiger sharks have the shortest sperm flagella length and total length of 27 sharks across seven orders. Therefore, evidence of positive selection in TCTEX1D2 and VAMP4 in tiger sharks could reflect the evolution of an increased sperm count and shorter sperm in response to increased sperm competition and post-copulatory sexual selection.
Positive selection on genes linked to brain development
Signatures of positive selection were detected in ARL6IP5, a transmembrane protein inhibiting EAAC1, the latter being a transporter of the excitatory neurotransmitter glutamate . EAAC1 is associated with neuron development in both vertebrates and invertebrates [48, 49]. When deprived of oxygen, levels of glutamate in neurons increase, causing neuronal death and potentially brain damage . Therefore, the ability to maintain glutamate levels below damaging thresholds may have enabled the exploitation of marginal habitats. Penetrating marginal habitats could be advantageous for these predators by allowing expansion into sub-optimal habitats where prey may take refuge . Contemporary evidence for this is demonstrated by the near-global distribution of tiger sharks in coastal and pelagic habitats around the world, as well as their diverse diet . Tiger sharks also spend considerable time in shallow seagrass and neritic habitats [53, 54], where dissolved oxygen concentrations fluctuate diurnally due to high productivity and demand .
The brain is one of the most energetically expensive organs to develop and maintain, and non-placental sharks are proposed to generally have smaller brains relative to body mass, compared to placental species [19, 56]. Animals need to maintain a balance between maintenance of the brain and other organs when exposed to hypoxic conditions . There is evidence of intra and interspecific variation in brain size of fishes exposed to different concentrations of dissolved oxygen, with larger brains evident in species inhabiting well-oxygenated waters [57–59].
Signatures of positive selection were also detected in YWHAE, a member of the 14-3-3 protein family . These proteins are expressed ubiquitously, particularly in the brain, and are highly conserved across animals [60, 61]. The 14-3-3 proteins are vital for differentiation of neurons in Drosophila, while mice lacking YWHAE have restricted brain development and neuronal migration [60, 62]. In humans, YWHAE is absent in sufferers of Miller-Dieker syndrome, which is characterised by severe mental disability . This suggests YWHAE is vital for brain development in humans and other animals. Positive selection in YWHAE could thus reflect a reduction in brain size in response to hypoxia experienced in sub-optimal habitats occupied by tiger sharks.
Positive selection in the TCTEX1D2, VAMP4, ARL6IP5, and YWHAE genes and their functional associations suggest they may have played a role in adaptation following the loss of placenta in tiger sharks. Noteworthy, changes in gene expression could have also been involved; however, as the time since species divergence increases, so do gene expression differences . Given the tiger shark and placental carcharhinids diverged c. 94 mya, the magnitude of gene expression changes is considerable, likely erasing any potential signal of gene expression changes associated with the loss of placenta. Furthermore, controlling the variety of factors (e.g. environment, sex, life history stage) potentially influencing gene expression would require a completely different sampling design to the one we employ here. Investigations of gene expression would also require uterus or yolk-sac tissue which are difficult to sample in a non-lethal way; hence, we sampled white muscle tissue instead. Furthermore, positive selection is intermittent in nature and signatures of it can be lost over time due to recombination and accumulating neutral substitutions [63, 64], thus some signatures of positive selection in other genes may no longer be detectable.
Previous studies have shown positive selection can act strongly on genes that greatly affect an individual’s fitness (e.g. sexual reproduction and sensory perception) [65, 66]. Thus, there could be additional factors explaining the signatures of positive selection exhibited in these four genes in the tiger shark lineage. We have here only one point of comparison, these genes might affect other traits of the tiger shark and may not be linked simply to the loss of placenta. Their functional associations, however, suggest these genes are good candidates for further study using additional genomic techniques.
Other studies have used genomics to test for signatures of positive selection associated with placental evolution in mammals [67, 68]. Evidence of positive selection in this group was found across 1,254 genes by Crosley et al. (2013) and in approximately 300 genes by Elliot & Crespi (2015). Signatures of positive selection were detected in a greater number of genes in these studies compared with our study; however, both Crosley et al. and Elliot & Crespi used a larger set of closely related species, and therefore looked at more genes than us, i.e. 16,578 genes in Elliot & Crespi (2013) and approximately 18,000 genes in Crosley et al. (2015) versus 1,197 genes here. In addition, Crosley et al. and Elliot & Crespi tested for positive selection along multiple lineages, whereas we tested for positive selection along the tiger shark lineage only. Interestingly, signatures of positive selection were also detected in genes associated with brain development and sexual reproduction in Elliot & Crespi and Crosley et al., respectively [67, 68], suggesting similar adaptations associated with placentation may have occurred in both mammals and sharks.
We were limited in the outgroups we could sample. We included as an outgroup one placental species (i.e. dusky smoothhound) from the other shark family (i.e. Triakidae) in which species may have also lost placenta. Unfortunately, we could not obtain samples from non-placental houndshark species, and so could not compare signatures of positive selection in genes along two independent lineages where placenta may have been lost. Also, we were unable to sample additional carcharhiniform species, but had lamniforms instead (which do not develop placenta). We identified more orthologous genes in the sand tiger transcriptome compared with other potential outgroup species, and so we used it as the second outgroup taxon.
Finally, another caveat of our study is the sampling of a single individual for each species during a short time frame. Consequently, the sequence data for each species lack individual, gender, ontogenetic, and temporal variation.
We suggest the tiger shark lineage split from the placental carcharhinids approximately 94 mya. We also propose at least four genes associated with brain development and sperm production have been evolving under positive selection in the tiger shark lineage, potentially reflecting adaptation following placental loss. Future work should utilise supplementary genomic techniques to investigate similar changes in reproductive traits of additional shark species.
In total, one individual from six carcharhinid species were sampled, the non-placental tiger shark and five placental carcharhinids: Atlantic sharpnose (Rhizoprionodon terraenovae), blacknose (Carcharhinus acronotus), bull (Carcharhinus leucas), Caribbean reef (Carcharhinus perezii), and lemon (Negaprion brevirostris) sharks. A placental dusky smoothhound individual was also sampled. Collection of white muscle samples was carried out off the coast of Eleuthera, the Bahamas (N 24° 50’ 05”: W 076° 19’ 32”) in a two-week period covering January and February 2014. The seven shark species were caught using 400 m stationary longlines with 30 to 33 non-offset, 16/0 circular hooks (Mustad, Gjövik, Norway). Hooks were spaced five metres apart and were baited with Atlantic bonito (Sarda sarda). The longlines were left in the water for approximately 90 min. White muscle samples were collected using a biopsy punch from an area adjacent to the dorsal fin. Samples were immediately placed in RNAlater (Sigma, St. Louis, MO, USA) and stored at 4 °C for 24 h, before storing at −20 °C. Samples were stored in an icebox for 16 h during transport back to the UK.
Sample collection was carried out under Cape Eleuthera Institute animal care protocols developed within the guidelines of the Association for the Study of Animal Behaviour and the Animal Behaviour Society . All sample collection activities were approved by The Bahamas Department of Marine Resources under research permits MAF/FIS/17 and MAF/FIS/34.
RNA extraction and purification
Muscle tissue samples were removed from RNAlater and individually left to dry in a Petri dish for five minutes before being cut up finely with a scalpel and individually homogenised in 300 μl of Buffer RLT (Qiagen, Hilden, Germany) using a PowerGen 120 Homogeniser (Fisher Scientific, Loughborough, UK). Total RNA was extracted using the RNeasy fibrous tissue mini kit (Qiagen, Hilden, Germany), following the manufacturers protocol.
RNA quantity was assessed using a NanoDrop 2000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA) and RNA quality assessed by electrophoresis on 1 % TAE-agarose gels. RNA samples were subjected to a cleaning and concentration phase using the RNeasy cleanup kit (Qiagen, Hilden, Germany), before being further assessed for quality and quantity using a Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). At least 500 ng of RNA with an RNA integrity number (RIN) of 6.0 or above, were stored at -80 °C prior to being sent for sequencing. In addition to the seven species sampled, filtered paired-end reads from the white muscle transcriptomes of the placental blue shark (Prionace glauca) (Order: Carcharhiniformes, Family: Carcharhinidae) and sand tiger shark were procured. For each species, a single individual was sampled and a transcriptome sequenced. The reads were generated following the same methods used for the other species in this study.
Sequencing, quality control, and de novo assembly
A normalised cDNA library was synthesised and sequenced for each species by BGI Tech Solutions (Hong Kong) using the Illumina TruSeq kit (Illumina, San Diego, CA, USA). These libraries were later sequenced using RNA-Seq and an Illumina HiSeq 2000 system following standard protocol. cDNA libraries were normalised in order to maximise transcript coverage and to sequence as much of the complete transcriptomes as possible. BGI Tech carried out an initial round of quality control, trimming raw reads, primers and adaptor sequences, as well as removing low quality reads (Phred quality < 20). FastQC (version 0.10.1) further assessed the filtered reads before assembly . The transcriptomes were de novo assembled using Trinity (version 2013-05-08) with default parameters .
Identification and alignment of orthologous sequences
The transcriptomes of all nine species were clustered separately using the Uclust function of Usearch (version 7.0.1090) with a similarity threshold of 98 % to breakdown and remove putative splice variants . The Trinity ORF predictor was employed using TransDecoder (version 2013-05-08) to calculate the longest and most probable translated region for each sequence and to remove multiple ORF sequences. A custom Perl script was used to remove all stop codons from the end of ORF sequences.
Pairwise reciprocal Blast searches were employed using Blastn (version 2.2.29) to identify putative orthologous sequences shared between each pair of species . Pairwise putative orthologues were then collated using R (version 3.1.1) in order to find sequences shared between all nine species. The multiple sequence alignment program M-Coffee (version 11.00.9103146) was used to align orthologous sequences between all species using Mafft, Muscle, T-Coffee, and Kalign methods . Gblocks was also used to remove poorly aligned and divergent regions to further reduce alignment errors and gaps . The Gblocks parameters used were: minimum of seven sequences for a conserved position, minimum of seven sequences for a flank position, maximum of six contiguous non-conserved positions, minimum block length of nine, and 50 % or more of sequences with a gap were treated as a gap position. M-Coffee translated nucleotide sequences into amino acid sequences, which were aligned for all species and then back-translated to nucleotide sequences. Alignments were then trimmed and graded. Only alignments with quality grades of nine, the highest score, were retained for further analysis.
Phylogenetic tree construction and positive selection analyses
The aligned orthologues were concatenated and used to construct a phylogenetic tree using RAxML (version 8.0.0) . We utilised the nucleotide substitution model GTRGAMMA, as determined by jModelTest (version 2.1.4) using the lowest value of the Akaike Information Criterion [75, 76]. Bootstrap values were calculated using 1,000 replicates. CodeML calculated the number of substitutions which alter the amino acid sequence (nonsynonymous (dN )) and the number of substitutions which do not alter the amino acid sequence (synonymous (dS)) . The ratio of these substitutions can be used to detect genes exhibiting signatures of positive selection. Positive selection, the favouring of distinct phenotypes, can be indicated by a dN / dS ratio (ω) > 1. ω < 1 is indicative of negative selection, the removal of deleterious alleles, and ω = 1 can indicate neutral selection, drift of alleles not affecting an individual’s ability to pass on their genes . It is suggested positive selection can only be detected if the average ω across all codon sites, calculated by the one ratio model (M0), is greater than 1. This is conservative, however, considering most codon sites will be highly conserved to maintain protein structure and function . Also, estimates of ω can be artificially decreased due to partial sequences produced by high-throughput sequencing technologies .
A maximum-likelihood site test based on a comparison between the neutral model 7 (M7) and the non-neutral model 8 (M8) was therefore employed. M7 assumes a β distribution of ω between 0 and 1, not allowing ω > 1 at any sites . M8 also assumes a β distribution of ω but allows an additional category of sites were ω > 1 . A CodeML likelihood ratio test (LRT) was used to test if M8 fits the data significantly better than M7. The natural log likelihood (lnL) values of M7 were contrasted with those of M8. The lnL ratios were then compared to a chi-squared distribution with two degrees of freedom. False positive results are possible when implementing positive selection analyses on a genomic scale, therefore the Benjamini-Hochberg false discovery rate (FDR) correction was applied to all p values [81, 82]. The BEB method was employed for sequences with significant M8:M7 likelihood ratios and used to identify codon sites exhibiting signatures of positive selection . Codon sites were considered to be showing evidence of positive selection if the probability of ω > 1 was more than 95 %.
We concatenated the 1,102 orthologues not showing signatures of positive selection for each species to construct a phylogeny in RAxML . We used the nucleotide substitution model GTRGAMMA, as determined by jModelTest . Bootstrap values for this phylogeny were calculated using 1,000 replicates. Fossil data were used to produce a time-calibrated phylogenetic tree (Table 2). This was accomplished using four independent runs of the MCMCTree function of PAML with 50,000 iterations, a burn in of 10,000 iterations, a sample frequency of three, an independent molecular clock, and the nucleotide substitution model HKY85, determined using jModelTest .
The branch-site test was utilised to detect signatures of positive selection at specific codon sites in the tiger shark lineage. The branch-site test is considered more powerful than the site test as signals of positive selection are not averaged over all branches of the phylogeny . For the branch-site test, an alternative hypothesis was contrasted to the null hypothesis using an LRT where the lnL values were compared to a chi-squared distribution with one degree of freedom. Estimates of ω were not determined using M0 as this model has been shown to be unreliable when detecting positive selection in specific branches [78, 83]. The branch-site test has also been known to experience convergence problems when calculating likelihoods, leading to artificial lnL ratios ; thus, three independent runs of this model were performed for both the alternative and null hypotheses, with the highest lnL values kept to calculate the lnL ratios. The Benjamini-Hochberg FDR correction was applied to all p values . The BEB method was employed to identify codon sites exhibiting signatures of positive selection . Codon sites were considered to be showing evidence of positive selection if the probability of ω > 1 was more than 95 %.
BEB, Bayes Empirical Bayes; bp, base pairs; dN, number of nonsynonymous substitutions; dS, number of synonymous substitutions; FDR, false discovery rate; HPD CI, 95 % highest posterior density credibility interval; lnL, natural log likelihood; LRT, likelihood ratio test; M0, one ratio model; M7, neutral model 7; M8, non-neutral model 8; mya, million years ago; ORFs, open reading frames; RIN, RNA integrity number; TAE, Tris base, acetic acid and EDTA; ω, dN/dS ratio
We thank staff at the Cape Eleuthera Institute for assistance in the field, Dominique Barthelemy and Jean Goasdoue for providing samples, and Helen Hipperson for assistance in the lab.
Availability of data and materials
The dataset of Illumina filtered paired-end reads supporting the results of this article are available on the NCBI Sequence Read Archive. Accession numbers: SAMN03333345 (Rhizoprionodon terraenovae), SAMN03333346 (Carcharhinus acronotus), SAMN03333347 (Prionace glauca), SAMN03333348 (Carcharhinus leucas), SAMN03333349 (Carcharhinus perezi), SAMN03333350 (Mustelus canis insularis), SAMN03333351 (Negaprion brevirostris), SAMN03333352 (Carcharias taurus) and SAMN03333353 (Galeocerdo cuvier).
DGS carried out the research and wrote the manuscript. LTD helped with transcriptome analyses. JI helped with phylogenetic analyses. EJB managed sample collection in the Bahamas. CSJ and LRN provided samples and advice regarding ecological context. AC helped with transcriptome assembly and orthologue detection. EH helped with RNA extraction. VS supervised the research. All authors commented on the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Open AccessThis 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.
- Blackburn DG. Convergent evolution of viviparity, matrotrophy, and specialisations for fetal nutrition in reptiles and other vertebrates. Am Zool. 1992;32:313–21.View ArticleGoogle Scholar
- Wourms JP, Lombardi J. Reflections on the Evolution of Piscine Viviparity. Am Zool. 1992;293:276–93.View ArticleGoogle Scholar
- Grogan ED, Lund R. Superfoetative viviparity in a Carboniferous chondrichthyan and reproduction in early gnathostomes. Zool J Linn Soc. 2011;161:587–94.View ArticleGoogle Scholar
- Blackburn DG: Evolution of Vertebrate Viviparity and Specializations for Fetal Nutrition: A Quantitative and Qualitative Analysis. J Morphol 2014: In Press.Google Scholar
- Mossman HW. Vertebrate Fetal Membranes: Comparative Ontogeny and Morphology; Evolution; Phylogenetic Significance; Basic Functions; Research Opportunities. New Brunswick, NJ: Rutgers University Press; 1987.View ArticleGoogle Scholar
- Zeh DW, Zeh JA. Reproductive mode and speciation: the viviparity-driven conflict hypothesis. BioEssays. 2000;22:938–46.View ArticlePubMedGoogle Scholar
- Crespi B, Semeniuk C. Parent-Offspring Conflict in the Evolution of Vertebrate Reproductive Mode. Am Nat. 2004;163:635–53.View ArticlePubMedGoogle Scholar
- Shine R. Does Viviparity Evolve in Cold Climate Reptiles Because Pregnant Females Maintain Stable (Not High) Body Temperatures? Evolution. 2004;58:1809–18.View ArticlePubMedGoogle Scholar
- Lambert SM, Wiens JJ. Evolution of viviparity: A phylogenetic test of the cold-climate hypothesis in phrynosomatid lizards. Evolution. 2013;67:2614–30.View ArticlePubMedGoogle Scholar
- Velo-Antón G, Santos X, Sanmartín-Villar I, Cordero-Rivera A, Buckley D. Intraspecific variation in clutch size and maternal investment in pueriparous and larviparous Salamandra salamandra females. Evol Ecol. 2015;29:185–204.View ArticleGoogle Scholar
- Dulvy NK, Reynolds JD. Evolutionary transitions among egg-laying, live-bearing and maternal inputs in sharks and rays. Proc R Soc B Biol Sci. 1997;264:1309–15.View ArticleGoogle Scholar
- Conrath CL, Musick JA. Reproductive Biology of Elasmobranchs. In: Carrier JC, Musick JA, Heithaus MR, editors. Biology of Sharks & their Relatives. London, UK: CRC Press; 2012.Google Scholar
- Venkatesh B, Kirkness EF, Loh Y, Halpern AL, Lee AP, Johnson J, Dandona N, Viswanathan LD, Tay A, Venter JC, Strausberg RL, Brenner S. Survey Sequencing and Comparative Analysis of the Elephant Shark (Callorhinchus milii) Genome. PLOS Biol. 2007;5:932–46.View ArticleGoogle Scholar
- Richards VP, Suzuki H, Stanhope MJ, Shivji MS. Characterization of the heart transcriptome of the white shark (Carcharodon carcharias). BMC Genomics. 2013;14:1–27.View ArticleGoogle Scholar
- Venkatesh B, Lee AP, Ravi V, Maurya AK, Lian MM, Swann JB, Ohta Y, Flajnik MF, Sutoh Y, Kasahara M, Hoon S, Gangu V, Roy SW, Irimia M, Korzh V, Kondrychyn I, Lim ZW, Tay B, Tohari S, Kong KW, Ho S, Lorente-Galdos B, Quilez J, Marques-Bonet T, Raney BJ, Ingham PW, Tay A, Hillier LW, Minx P, Boehm T, et al. Elephant shark genome provides unique insights into gnathostome evolution. Nature. 2014;505:174–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Stingo V, Rocco L. Selachian cytogenetics: a review. Genetica. 2001;111:329–47.View ArticlePubMedGoogle Scholar
- Howe K, Clark MD, Torroja CF, Torrance J, Berthelot C, Muffato M, Collins JE, Humphray S, Mclaren K, Matthews L, Mclaren S, Sealy I, Caccamo M, White S, Chow W, Kilian B, Churcher C, Scott C, Barrett JC, Koch R, Zhou Y, Gu Y, Yen J, Vogel J, Eyre T, Redmond S, Banerjee R, Chi J, Fu B, Langley E, et al. The zebrafish reference genome sequence and its relationship to the human genome. Nature. 2013;496:2–8.View ArticleGoogle Scholar
- Hamlett W, Hysell M. Uterine specializations in elasmobranchs. J Exp Zool. 1998;282:438–59.View ArticleGoogle Scholar
- Mull CG, Yopak KE, Dulvy NK. Does more maternal investment mean a larger brain? Evolutionary relationships between reproductive mode and brain size in chondrichthyans. Mar Freshw Res. 2011;62:567–75.View ArticleGoogle Scholar
- Compagno L, Dando M, Fowler S. Sharks of the World. London, UK: Princeton Press; 2005.Google Scholar
- Castro JI, Sato K, Bodine AB. A novel mode of embryonic nutrition in the tiger shark, Galeocerdo cuvier. Mar Biol Res. 2016;1000:1–6.Google Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang Z. PAML 4: Phylogenetic Analysis by Maximum Likelihood. Mol Biol Evol. 2007;24:1586–91.View ArticlePubMedGoogle Scholar
- Yang Z, Wong WSW, Nielsen R. Bayes Empirical Bayes Inference of Amino Acid Sites Under Positive Selection. Mol Biol Evol. 2005;22:11007–118.Google Scholar
- Probst J. Beiträge zur Kenntniss der fossilen Fische aus der Molasse von Baltringen. Hayfische. Jahreshefte des Vereins für vaterländische Naturkd Württemb. 1879;35:127–91.Google Scholar
- Landini W. Revisione degli “ittiodontoliti pliocenici” della collezione Lawley. Palaeontogr Ital. 1977;70:92–134.Google Scholar
- Delsate D, Lepage J. Decouverte d’une faune originale d'elasmobranches dans les phosphates du Toarcien Lorrain (couche á Coeloceras crassum). Bull l’Académie Société Lorraine des Sci. 1990;29:153–61.Google Scholar
- Noubhani A, Cappetta H. Les Orectolobiformes, Carchariniformes, et Myliobatiformes (Elasmobranchii, Neoselachii) des Bassins a phosphate du Maroc (Maastrichtian-Lutétien basal). Palaeo Ichthyol. 1997;8:1–327.Google Scholar
- Underwood CJ, Mitchell SF, Veltkamp KJ. Shark and ray teeth from the Hauterivian (Lower Cretaceous) of north-east England. Palaeontology. 1999;42:287–302.View ArticleGoogle Scholar
- Underwood CJ, Ward DJ. Neoselachian sharks and rays from the British Bathonian (Middle Jurassic). Palaeontology. 2004;47:447–501.View ArticleGoogle Scholar
- Gradstein FM, Ogg JG, Schmitz MD, Ogg GM. A Geologic Time Scale 2012. Boston: Elsevier; 2012.Google Scholar
- Sorenson L, Santini F, Alfaro ME. The effect of habitat on modern shark diversification. J Evol Biol. 2014;27:1536–48.View ArticlePubMedGoogle Scholar
- Zangerl R. Handbook of Paleoichthyology 3. Stuttgart, Germany: Gustav Fischer; 1981.Google Scholar
- Grogan ED, Lund R, Greenfest-Allen E. The Origin and Relationships of Early Chondricthyans. In: Carrier JC, Musick JA, Heithaus MR, editors. Biology of Sharks & their Relatives. London, UK: CRC Press; 2012.Google Scholar
- Motta PJ, Huber DR. Prey Capture Behavior and Feeding Mechanics of Elasmobranchs. In: Carrier J, Musick JA, Heithaus M, editors. Biology of Sharks & their Relatives. London, UK: CRC Press; 2012.Google Scholar
- Whitney NM, Crow GL. Reproductive biology of the tiger shark (Galeocerdo cuvier) in Hawaii. Mar Biol. 2007;151:63–70.View ArticleGoogle Scholar
- Inaba K, Kagami O, Ogawa K. Tctex2-related outer arm dynein light chain is phosphorylated at activation of sperm motility. Biochem Biophys Res Commun. 1999;256:177–83.View ArticlePubMedGoogle Scholar
- Pazour GJ, Koutoulis A, Benashski SE, Bethany L, Sheng H, Patel-king RS, King SM, Witman GB. LC2, the Chlamydomonas Homologue of the t Complex-encoded Protein Tctex2, Is Essential for Outer Dynein Arm Assembly. Mol Biol Cell. 1999;10:3507–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Asante D, Stevenson NL, Stephens DJ. Subunit composition of the human cytoplasmic dynein-2 complex. J Cell Sci. 2014;127:4774–87.View ArticlePubMedPubMed CentralGoogle Scholar
- Rashid S, Grzmil P, Drenckhahn J-D, Meinhardt A, Adham I, Engel W, Neesen J. Disruption of the murine dynein light chain gene Tcte3-3 results in asthenozoospermia. Reproduction. 2010;139:99–111.View ArticlePubMedGoogle Scholar
- Guo X, Shen J, Xia Z, Zhang R, Zhang P, Zhao C, Xing J, Chen L, Chen W, Lin M, Huo R, Su B, Zhou Z, Sha J. Proteomic Analysis of Proteins Involved in Spermiogenesis in Mouse. J Proteome Res. 2010;9:1246–56.View ArticlePubMedGoogle Scholar
- Parker GA. Sperm competition and its evolutionary consequences in the insects. Biol Rev. 1970;45:525–67.View ArticleGoogle Scholar
- Tanaka S, Kurokawa H, Hara M. Comparative morphology of the sperm in chondrichthyan fishes. In: Jamieson BGM, Ausio J, Justine JL, editors. Advances in Spermatozoal Phylogeny and Taxonomy. Paris, France: Mémoires du Muséum National d’Histoire Naturelle; 1995. p. 313–20.Google Scholar
- Hunter FM, Birkhead TR. Sperm viability and sperm competition in insects. Curr Biol. 2002;12:121–3.View ArticlePubMedGoogle Scholar
- Fitzpatrick JL, Baer B. Polyandry reduces sperm length variation in social insects. Evolution. 2011;65:3006–12.View ArticlePubMedGoogle Scholar
- Choi BY, Kim JH, Kim HJ, Lee BE, Kim IY, Sohn M, Suh SW. EAAC1 gene deletion increases neuronal death and blood brain barrier disruption after transient cerebral ischemia in female mice. Int J Mol Sci. 2014;15:19444–57.View ArticlePubMedPubMed CentralGoogle Scholar
- Maleszka R, Helliwell P, Kucharski R. Pharmacological interference with glutamate re-uptake impairs long-term memory in the honeybee, Apis mellifera. Behav Brain Res. 2000;115:49–53.View ArticlePubMedGoogle Scholar
- Katagiri H, Tanaka K, Manabe T. Requirement of appropriate glutamate concentrations in the synaptic cleft for hippocampal LTP induction. Eur J Neurosci. 2001;14:547–53.View ArticlePubMedGoogle Scholar
- Takahashi M, Billups B, Rossi D, Sarantis M, Hamann M, Attwell D. The role of glutamate transporters in glutamate homeostasis in the brain. J Exp Biol. 1997;200:401–9.PubMedGoogle Scholar
- Jorgensen SJ, Klimley AP, Muhlia-Melo AF. Scalloped hammerhead shark Sphyrna lewini, utilizes deep-water, hypoxic zone in the gulf of California. J Fish Biol. 2009;74:1682–7.View ArticlePubMedGoogle Scholar
- Papastamatiou YP, Wetherbee BM, Lowe CG, Crow GL. Distribution and diet of four species of carcharhinid shark in the Hawaiian Islands: Evidence for resource partitioning and competitive exclusion. Mar Ecol Prog Ser. 2006;320:239–51.View ArticleGoogle Scholar
- Afonso AS, Hazin FH V. Vertical Movement Patterns and Ontogenetic Niche Expansion in the Tiger Shark, Galeocerdo cuvier. PloS One. 2015;10:e0116720.View ArticlePubMedPubMed CentralGoogle Scholar
- Heithaus MR, Dill LM, Marshall GJ, Buhleier B. Habitat use and foraging behavior of tiger sharks (Galeocerdo cuvier) in a seagrass ecosystem. Mar Biol. 2002;140:237–48.View ArticleGoogle Scholar
- Ziegler S, Benner R. Ecosystem metabolism in a subtropical, seagrass-dominated lagoon. Mar Ecol Prog Ser. 1998;173:1–12.View ArticleGoogle Scholar
- Elliot MG, Crespi BJ. Placental invasiveness and brain-body allometry in eutherian mammals. J Evol Biol. 2008;21:1763–78.View ArticlePubMedGoogle Scholar
- Chapman LJ, Hulen KG. Implications of hypoxia for the brain size and gill morphometry of mormyrid fishes. J Zool. 2001;254:461–72.View ArticleGoogle Scholar
- Chapman LJ, Galis F, Shinn J. Phenotypic plasticity and the possible role of genetic assimilation: Hypoxia-induced trade-offs in the morphological traits of an African cichlid. Ecol Lett. 2000;3:387–93.View ArticleGoogle Scholar
- Chapman L, Albert J, Galis F. Developmental Plasticity, Genetic Differentiation, and Hypoxia-induced Trade-offs in an African Cichlid Fish. Open Evol J. 2008;2:75–88.View ArticleGoogle Scholar
- Toyo-oka K, Shionoya A, Gambello MJ, Cardoso C, Leventer R, Ward HL, Ayala R, Tsai L, Dobyns W, Ledbetter D, Hirotsune S, Wynshaw-boris A. 14-3-3 ε is important for neuronal migration by binding to NUDEL: a molecular explanation for Miller-Dieker syndrome. Nat Genet. 2003;34:274–85.View ArticlePubMedGoogle Scholar
- Tzivion G, Avruch J. 14-3-3 Proteins: Active Cofactors in Cellular Regulation By Serine/Threonine Phosphorylation. J Biol Chem. 2002;277:3061–4.View ArticlePubMedGoogle Scholar
- Chang H, Rubin G. 14-3-3 epsilon positively regulates Ras-mediated signaling in Drosophila. Genes Dev. 1997;11:1132–9.View ArticlePubMedGoogle Scholar
- Reich DE, Goldstein DB. Microsatellites: Evolution and Applications. Oxford, UK: Oxford University Press; 1999.Google Scholar
- O’Neill MJ, Lawton BR, Mateos M, Carone DM, Ferreri GC, Hrbek T, Meredith RW, Reznick DN, O’Neill RJ. Ancient and continuing Darwinian selection on insulin-like growth factor II in placental fishes. Proc Natl Acad Sci U S A. 2007;104:12404–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Duret L. Neutral theory: the null hypothesis of molecular evolution. Nat Educ. 2008;1:218.Google Scholar
- Kosiol C, Vinar T, da Fonseca RR, Hubisz MJ, Bustamante CD, Nielsen R, Siepel A. Patterns of Positive Selection in Six Mammalian Genomes. PLoS Genet. 2008;4:e1000144.View ArticlePubMedPubMed CentralGoogle Scholar
- Crosley EJ, Elliot MG, Christians JK, Crespi BJ. Placental invasion, preeclampsia risk and adaptive molecular evolution at the origin of the great apes: Evidence from genome-wide analyses. Placenta. 2013;34:127–32.View ArticlePubMedGoogle Scholar
- Elliot MG, Crespi BJ, Elliot MG. Genetic recapitulation of human pre-eclampsia risk during convergent evolution of reduced placental invasiveness in eutherian mammals. Philos Trans R Soc B Biol Sci. 2015;370:20140069.View ArticleGoogle Scholar
- Rollin B, Kessel M. Guidelines for the treatment of animals in behavioural research and teaching. Anim Behav 1998;55:251-57.Google Scholar
- Andrews S: FastQC, a quality control tool for high throughput sequence data. Available at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc. 2010. Accessed May 2014.
- Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.View ArticlePubMedGoogle Scholar
- Altshul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic Local Alignment Search tool. J Mol Biol. 1990;215:403–10.View ArticleGoogle Scholar
- Wallace IM, O’Sullivan O, Higgins DG, Notredame C. M-Coffee: combining multiple sequence alignment methods with T-Coffee. Nucleic Acids Res. 2006;34:1692–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Castresana J. Selection of Conserved Blocks from Multiple Alignments for Their Use in Phylogenetic Analysis. Mol Biol Evol. 2000;17:540–52.View ArticlePubMedGoogle Scholar
- Akaike H. A new look at the statistical model identification. IEEE Trans Automat Contr. 1974;19:716–23.View ArticleGoogle Scholar
- Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9:722.View ArticleGoogle Scholar
- Nielsen R. Molecular signatures of natural selection. Annu Rev Genet. 2005;39:197–218.View ArticlePubMedGoogle Scholar
- Yang Z, dos Reis M. Statistical properties of the branch-site test of positive selection. Mol Biol Evol. 2011;28:1217–28.View ArticlePubMedGoogle Scholar
- Brieuc MS, Naish K. Detecting signatures of positive selection in partial sequences generated on a large scale: pitfalls, procedures and resources. Mol Ecol Resour. 2011;11:172–83.View ArticlePubMedGoogle Scholar
- Metzger KJ, Thomas MA. Evidence of positive selection at codon sites localized in extracellular domains of mammalian CC motif chemokine receptor proteins. BMC Evol Biol. 2010;10:139–48.View ArticlePubMedPubMed CentralGoogle Scholar
- Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57:289–300.Google Scholar
- Roux J, Privman E, Moretti S, Daub JT, Robinson-Rechavi M, Keller L. Patterns of positive selection in seven ant genomes. Mol Biol Evol Adv Access. 2014;31:1661–85.View ArticleGoogle Scholar
- Zhang J, Nielsen R, Yang Z. More genes underwent positive selection in chimpanzee evolution than in human evolution. PNAS. 2007;104:7489–94.View ArticlePubMedPubMed CentralGoogle Scholar