Introduction
Chondrichthyan fishes are composed of two subclasses, Holocephali and Elasmobranchii. Holocephalans are the more basal of the pair having first appeared more than 400 million years ago and include a single surviving order, Chimaeriformes, the chimaeras, with 39 extant species1. Elasmobranchs appeared approximately 350 million years ago and include more than 1000 species of sharks, skates, and rays2. Chondrichthyan fishes occupy a pivotal position at the base of the vertebrate phylogenetic tree. For research that includes an evolutionary component, representation of this diverse class affords a valuable perspective to evaluate all vertebrates.
Chondrichthyan fishes are circumglobal in distribution and occupy a wide range of ecological habitats. Their life history parameters are equally disparate but in general chondrichthyans are slow growing and late maturing fishes with an increased risk of extinction3–5. Fecundity is as few as 1 or 2 for viviparous species such as the sand tiger shark, Carcharias taurus6 and as high as 300 for the whale shark, Rhincodon typus7. They are of economic importance for fisheries as well as ecotourism. Management and assessment of stock is essential to ensure both ecotourism interests and food resources remain sustainable8. Management of fish populations has increasingly relied on molecular tools to investigate population structure, properly identify species, and compliance with fishing quotas9–13.
Elasmobranchs have been used as a model for biomedical research for more than 100 years. Elasmobranchs, like other cartilaginous fishes, exhibit many fundamental vertebrate characteristics, including a neural crest, jaws and teeth, an adaptive immune system, and a pressurized circulatory system. The skate is a powerful comparative model to study biological processes shared among jawed and limbed vertebrates such as development14–16, renal physiology17–20, immunology21–26, toxicology27, neurobiology28, and wound healing and regeneration29. They are the most ancient vertebrates to posses an adaptive immune system that generates antibodies using a V(D)J combinatorial mechanism30. Phylogenetically, cartilaginous fishes are the first vertebrates to possess a thymus, a central lymphoid organ that provides a microenvironment for the development of T cells31. The thymus shares a common organization with more derived vertebrates containing cortical and medullary regions32,33.
In addition to shared physiological characteristics, the diversity of specializations between species allows investigations of evolution within a single clade. For example, elasmobranchs use a plethora of reproductive strategies that span the full range of maternal investment from placental viviparity to strict lecitrophic oviparity. Besides sexual reproduction, captive elasmobranchs are capable of asexual parthenogenesis34–36. Of these reproductive mechanisms, the most tractable for research purposes is oviparity. Approximately 43% of chondrichthyans utilize oviparity including all Chimaeriformes, Heterodontiformes (bullhead sharks), Rajoidae (skates) and Scyliorhinidae (catsharks)37. Many species can be maintained in captivity and will breed and lay eggs throughout an annual season38. Artificial insemination has been reported for two oviparous species, the clearnose skate, Raja eglanteria39, and the cloudy catshark, Scyliorhinus torazame40. Additionally, sperm storage allows wild caught females to lay eggs for several years without requiring males or captive mating events41.
Leucoraja erinacea, the little skate, was chosen for a genome sequencing project to represent this clade of fishes because of their use as a biomedical model, experimental tractability, genome size, existing sequence data, and northeast regional distribution. The sequencing project is an ongoing effort of the North East Bioinformatics Collaborative (NEBC) of the North East Cyberinfrastructure Consortium (NECC), composed of the bioinformatics core facilities from Delaware, Maine, New Hampshire, Rhode Island, and Vermont funded by National Institutes of Health (NIH) Institutional Development Awards (IDeA) and/or National Science Foundation (NSF) Experimental Program to Stimulate Competitive Research (EPSCoR) programs.
Existing resources
There is a single order of holocephalans and 13 orders of elasmobranchs. The distribution of species in orders, families and genera is shown in Figure 1. The batoids are composed of 4 orders, Rajiiformes, Myliobatiformes, Torpidiformes, and Rhinopristiformes, and contain 54% of extant chondrichthyan species. Sharks are broadly divided into two super orders, Galeomorphii and Squalomorphii that together account for 43% of extant chondrichthyan species. The galean sharks include 4 orders: Heterdontiformes, Orectolobiformes, Lamniormes and Carcharhiniformes, and represent 30% of extant chondrichthyan species. Squalean sharks are composed of 4 orders: Squaliformes, Squatiniformes, Pristophoriformes, and Hexanchiformes, comprising 13% of extant chondrichthyan species. Among individual orders, Rajiiformes, the skates, have the most species (345) followed by Carcharhiniformes, the ground sharks (283) and Myliobatiformes (226)2. These ‘big three’ orders contain 854 species, 72% of extant chondrichthyans.
Figure 1. Species distribution within chondrichthyan orders.
There is a single order of Holocephalans, Chimaeriformes, and 13 orders of elasmobranchs. The distribution of chondrichthyan species in each of the 14 orders is shown relative to the total number of species, genera and families for the clade. The batoids are composed of 4 orders, Rajiiformes, Myliobatiformes, Torpidiformes, and Rhinopristiformes, and contain 54% of extant chondrichthyan species. Sharks are broadly divided into two super orders, Galeomorphii and Squalomorphii that together include the remaining 9 orders and 43% of extant chondrichthyan species.
Chondrichthyan conservation, management, and research all benefit from easily accessible and well-documented molecular resources. The organization of data and metadata in archival databases is critically important for efficient use of large and complex datasets. The International Nucleotide Sequence Database Collaboration (INSDC) is composed of three large public nucleotide repositories, DNA Data Bank of Japan (DDBJ), European Molecular Biology Laboratory-European Bioinformatics Institute (EMBL-EBI), and GenBank at the National Center for Biotechnology Information (NCBI). Recently, two new NCBI database projects were initiated to collect details of samples, BioSample, and project data, BioProject, and propagate the metadata to all associated database entries in an effort to expand the use of already existing and rapidly expanding molecular resources42. Figure 2 illustrates the relationship between BioProject, BioSample and the sequence data for SkateBase. Because the BioProject and BioSample databases were established in 2012, not all existing datasets have metadata or details of the biological source to populate a BioSample and BioProject entry. When available, BioProject and BioSample hyperlinks are included for Sequence Read Archive (SRA), Expressed Sequence Tag (EST) and Genome Survey Sequence (GSS) datasets in the tables below.
Figure 2. Representation of SkateBase data within the The National Center for Biotechnology Information (NCBI) databases.
A. The little skate genome project is represented as a BioProject entry that connects all samples and data thematically. A BioSample record describes the DNA sample that was used for genome sequencing that was generated from a single stage 32 skate embryo. The SRA catalogs the unassembled Illumina genome sequence data. The Whole Genome Shotgun (WGS) database contains the contiguous sequences from shotgun sequencing projects. The assembled and annotated mitochondrial genome was deposited in GenBank and subsequently included in the NCBI Reference Sequence Database (RefSeq). B. The project to characterize the embryonic transcriptomes of L. erinacea, C. milii and S. canicula is represented in a BioProject entry. Three BioSample entries, one for each species, lead to three SRA datasets. The transcriptome data is represented also in the Gene Expression Omnibus (GEO), a database of high-throughput functional genomic data derived from microarrays and next-generation sequencing technologies.
Table 1 is a summary of chondrichthyan sequence data in NCBI databases, UniProtKB, and the Protein Data Bank (PDB) with L. erinacea, Callorhinchus milii and Scyliorhinus canicula, the three species featured at SkateBase listed individually. The distribution of holocephalans and elasmobranchs in public databases is illustrated in Figure 3. Despite the majority of species belonging to Elasmobranchii, the GenBank, UniProtKB/TrEMBL, and Gene databases are dominated by chimaera data derived from the genome sequence of the elephant shark, C. milii43. Elasmobranch data predominates in UniProtKB/Swiss-Prot, PDB, BioProject and BioSample databases as well as the number of whole mitochondrial genomes (WMG) in GenBank. The EST and SRA databases are nearly equally split between the two subclasses.
Table 1. Chondrichthyan molecular sequence data in public databases.
| National Center for Biotechnology Information (NCBI) databases1 | | |
---|
| | | | | GenBank | | | | UniProtKB2 | |
---|
| Taxonomy | BioProject | BioSample | Gene | GenBank | WMG | EST | EST lib | GSS | GSS lib | WGS (Mbp) | GEO3 | SRA | Swiss- Prot | TrEMBL | PDB |
---|
Chondrichthyes | 7777 | 16 | 75 | 21069 | 55810 | 72 | 192948 | 33 | 28497 | 5 | 2492.3 | 3 | 22 | 276 | 26485* | 178 |
Holocephali | 7863 | 3 | 21 | 20201 | 39512 | 8 | 109965 | 6 | 27944 | 1 | 936.9 | 1 | 13 | 12 | 20170 | 0 |
C. milii
| 7868 | 3 | 21 | 20110 | 39232 | 1 | 109965 | 6 | 27944 | 1 | 936.9 | 1 | 13 | 3 | 19989 | 0 |
Elasmobranchii | 7778 | 13 | 54 | 868 | 16273 | 64 | 82983 | 27 | 553 | 4 | 1555.4 | 2 | 9 | 264 | 6299 | 178 |
L. erinacea
| 7782 | 3 | 7 | 13 | 284 | 1 | 31167 | 5 | 0 | 0 | 1555.4 | 1 | 2 | 6 | 123 | 0 |
S. canicula
| 7830 | 2 | 8 | 13 | 645 | 1 | 1600 | 7 | 0 | 0 | 0 | 1 | 1 | 38 | 283 | 1 |
Figure 3. Holocephalan and elasmobranch resources in public nucleotide and protein databases.
The distribution of data for Holocephalii (chimaeras) and elasmobranchii (sharks and rays) subclasses of chondrichthyan fishes does not always reflect their species distribution. The number of species represented in GenBank is representative of the actual species distribution but the amount of data in GenBank is not. Holocephalan data forms the majority of the NCBI Gene, GenBank, Genome Survey Sequence (GSS) and UniProt TrEMBL databases. The number of Short Reach Archive (SRA) experiments and EST sequences in nearly equal for each subclass and the remaining databases are primarily populated by elasmobranch data.
Chondrichthyan genomes
Currently there are multiple efforts to sequence an elasmobranch genome in various stages of completion (Table 2); however, only the skate genome project currently has data publically available. Efforts to sequence the whale shark are underway at the Georgia Aquarium and Emory University (personal communication, Alistair Dove, Georgia Aquarium). Genoscope leads a project to sequence the genome of another oviparous elasmobranch, the catshark, S. canicula. The current assembly is described in Table 2. A second version of the catshark genome with 200x coverage, including mate pair sequencing, is in progress (personal communication, Sylvie Mazan, French National Centre for Scientific Research). Among holocephalans, the genome of the elephant shark, C. milii, was first described in a 1.4x coverage assembly in 200644. With continued sequencing the assembly coverage is currently 19.25x and data has been made available through the project website (http://esharkgenome.imcb.a-star.edu.sg/) and Genbank43.
A powerful resource for characterizing genomes is large-insert clone libraries where each clone contains a large (~100kb) genomic region. Bacterial artificial chromosome (BAC) and P1-derived artificial chromosome (PAC) libraries are DNA constructs within a plasmid used to transform bacteria. As the bacteria grow the inserted DNA is amplified and subsequently isolated and sequenced. BACs are beneficial for genome sequencing projects because the insert size can be very large, nearly 350 kb, facilitating assembly post-sequencing. BAC/PAC libraries were built for several chondrichthyan species including the nurse shark, Ginglymostoma cirratum45; elephant shark, C. milli46; little skate, L. erinacea47; horn shark, Heterdontus francisci48; dogfish shark, Squalus acanthias49,50; and catshark, S. canicula51. These libraries were used to successfully characterize a handful of genomic regions such as little skate HoxA cluster47,52, horn shark Hox A and D clusters53, catshark HoxA, B and D clusters51,52, C. milii HoxA-D clusters54, immunoglobulin receptor IgW C regions30 and neurohypophysial gene loci46.
RNA databases
Transcriptome sequencing seeks to characterize all genes expressed in a tissue or set of tissues in a sample. Technologies to identify the complete RNA transcript sequence have developed from studies of a small number of transcripts to comprehensive characterizations. The application of large-scale cDNA cloning of Expressed Sequence Tags (ESTs) gave initial characterizations of 5-prime and/or 3-prime ends of transcripts in several elasmobranchs including L. erinacea and S. acanthias (Table 3). EST sequence data are available in the EST divisions of the GenBank, EMBL and DDBJ databases that make up the INSDC. cDNA clones and their sequences from these EST projects have enabled the complete characterization of the full-length cDNA sequence of several genes. In the last five years, high-throughput RNA sequencing (RNA-Seq) has been applied to comprehensively examine the complete sequence of transcripts in tissues of cartilaginous fishes. Among the most valuable RNA-Seq datasets are those from whole late-stage embryos following organogenesis. Our project has generated these datasets for L. erinacea, S. canicula and C. milii52. Public RNA-Seq data sets can be found in the NCBI Gene Expression Omnibus and Short Read Archive (SRA) databases or the EBI ArrayExpress and European Nucleic Acid (ENA) archives (Table 3 and Table 4).
Table 3. National Center for Biotechnology Information (NCBI) Expressed Sequence Tags (EST) and Genome Survey Sequences (GSS) databases (release 130101): Chondrichthyan sequence data.
BioSample | BioSample Description | Library ID | Organism | Sample age/sex | Sample type | ESTs | Facility1 | Date |
---|
Holocephali |
Chimaeriformes |
182978
| Whole-genome shotgun library of the elephant shark (aka elephant fish) | GSS: LIBGSS_009694 |
Callorhinchus milii
| - | testis | 27944 | IMCB | 2004 |
1000678
| Elephant shark full- length cDNA library from testis | EST: LIBEST_027873 |
Callorhinchus milii
| - | testis | 29234 | IMCB | 2012 |
1000677
| Elephant shark full- length cDNA library from spleen | EST: LIBEST_027872 |
Callorhinchus milii
| - | spleen | 16664 | IMCB | 2012 |
1000676
| Elephant shark full- length cDNA library from liver | EST: LIBEST_027871 |
Callorhinchus milii
| - | liver | 16573 | IMCB | 2012 |
1000675
| Elephant shark full- length cDNA library from kidney | EST: LIBEST_027870 |
Callorhinchus milii
| - | kidney | 19246 | IMCB | 2012 |
1000674
| Elephant shark full- length cDNA library from intestine | EST: LIBEST_027869 |
Callorhinchus milii
| - | intestine | 12146 | IMCB | 2012 |
1000673
| Elephant shark full- length cDNA library from gills | EST: LIBEST_027868 |
Callorhinchus milii
| - | gills | 16012 | IMCB | 2012 |
Elasmobranchii: Batoids (rays and skates) |
Torpediformes |
158311
| Torpedo marmorata electric organ | EST: LIBEST_003755 |
Torpedo marmorata
| - | electric organ | 8 | CNRS | 2000 |
158310
| Torpedo marmorata electric lobe | EST: LIBEST_003754 |
Torpedo marmorata
| - | electric lobe | 26 | CNRS | 2000 |
157461
| pFL61-TEL | EST: LIBEST_002905 |
Torpedo marmorata
| - | electric lobe | 2 | CNRS | 2000 |
157406
| pFL61-EL | EST: LIBEST_002849 |
Torpedo marmorata
| - | electric lobe | 5 | CNRS | 2000 |
154382
| Torpedo californica electric organ | EST: LIBEST_020696 |
Torpedo californica
| - | electric organ | 10185 | Children’s National Medical Center, USA | 2006 |
Rajiformes |
175126
| Little Skate Multiple Tissues, Normalized | EST: LIBEST_015890 |
Leucoraja erinacea
| adult | mixeda | 5698 | MDIBL | 2004 |
176484
| Little Skate Liver, Normalized | EST: LIBEST_017626 |
Leucoraja erinacea
| adult | liver | 6016 | MDIBL | 2005 |
165533
| Little Skate embryo cell line 1 (LEE-1): 5' sequences | EST: LIBEST_022984 |
Leucoraja erinacea
| embryonic cell line | stage 28 | 4825 | MDIBL | 2006 |
154366
| Little skate embryo tissues; 5' sequences | EST: LIBEST_020422 |
Leucoraja erinacea
| embryo | stage 19, 20, 25 | 5600 | MDIBL | 2006 |
166469
| Skate Multiple Tissues, Normalized | EST: LIBEST_023576 |
Leucoraja erinacea
| adult | mixeda | 9028 | MDIBL | 2008 |
Elasmonbranchii: Selachii (sharks)
|
Carcharhiniformes |
168576
| Dogfish testis - round spermatids zone (SSH) | EST: LIBEST_025578 |
Scyliorhinus canicula
| adult | testis | 20 | Caen University | 2009 |
168575
| Dogfish testis - spermatogonia zone (SSH) | EST: LIBEST_025577 |
Scyliorhinus canicula
| adult | testis | 12 | Caen University | 2010 |
222714
| Scyliorhinus canicula juvenile library | EST: LIBEST_026904 |
Scyliorhinus canicula
| juvenile | 5 days post-hatch | 56 | enoscope-CEA | 2011 |
222713
| Scyliorhinus canicula embryonic, stages 9–15 library | EST: LIBEST_026903 |
Scyliorhinus canicula
| embryo | stages 9–15 | 628 | Genoscope-CEA | 2011 |
222712
| Scyliorhinus canicula embryonic, stages 19–25 library | EST: LIBEST_026902 |
Scyliorhinus canicula
| embryo | stages 19–25 | 772 | Genoscope-CEA | 2011 |
222711
| Scyliorhinus canicula embryonic, stages 19–24 library | EST: LIBEST_026901 |
Scyliorhinus canicula
| embryo | stages 19–24 | 33 | Genoscope-CEA | 2011 |
222710
| Scyliorhinus canicula adult brain library | EST: LIBEST_026900 |
Scyliorhinus canicula
| adult | brain | 79 | Genoscope-CEA | 2011 |
699400
| cloudy catshark embryo cDNA library | EST: LIBEST_027410 |
Scyliorhinus torazame
| embryo | stage 31 | 2942 | RIKEN | 2011 |
Orectolobiformes |
183175
| GC__Ba | GSS: LIBGSS_009945 |
Ginglymostoma cirratum
| adult | red blood cells | 178 | University of Arizona | 2005 |
184343
| shark whole genome shotgun library 2 | GSS: LIBGSS_011249 |
Chiloscyllium plagiosum
| female | ventral fin | 177 | Tgen | 2008 |
184342
| shark whole genome shotgun library 1 | GSS: LIBGSS_011248 |
Chiloscyllium plagiosum
| female | ventral fin | 194 | Tgen | 2008 |
166749
| Shark liver regeneration | EST: LIBEST_023789 |
Chiloscyllium plagiosum
| adult | liver | 2103 | BGI | 2008 |
176026
| cDNA library of Shark hepatic regeneration tissues | EST: LIBEST_017019 |
Chiloscyllium plagiosum
| none | Hour 24 after 2/3 partial hepatectomy | 17 | CPU | 2005 |
254067
| Toll like receptor ligand induced Spleen | EST: LIBEST_027180 |
Chiloscyllium griseum
| male | spleen | 1051 | MVC | 2011 |
254066
| Spleen of Chiloscyllium griseum | EST: LIBEST_027179 |
Chiloscyllium griseum
| male | spleen | 1000 | MVC | 2011 |
1797282
| Suppressive subtractive hybridization library from peptidoglycan induced spleen of the shark | EST: LIBEST_028031 |
Chiloscyllium griseum
| male | spleen | 315 | MVC | 2012 |
Squaliformes |
175664
| Dogfish Shark Multiple Tissues, Normalized | EST: LIBEST_016552 |
Squalus acanthias
| adult | mixedb | 15078 | MDIBL | 2004 |
176998
| Dogfish Shark Embryo-derived Cell Line SAE, Normalized | EST: LIBEST_018195 |
Squalus acanthias
| embryonic cell line | embryo with external yolk sac | 5824 | MDIBL | 2005 |
154362
| Spiny dogfish shark rectal gland EST library | EST: LIBEST_020417 |
Squalus acanthias
| - | rectal gland | 5085 | MDIBL | 2006 |
150616
| Dogfish Shark Rectal Gland, Normalized | EST: LIBEST_020023 |
Squalus acanthias
| adult | rectal gland | 6575 | MDIBL | 2006 |
Hexanchiformes |
178140
| Hexanchus griseus DNA (Hunter C) | GSS: LIBGSS_003277 |
Hexanchus griseus
| - | - | 4 | HGMP-RC | 2001 |
Mitochondrial genomes
Individual mitochondrial genes such as cytochrome c oxidase subunit I (CO1 or COX1) and NADH-ubiquinone oxidoreductase chain 2 (NADH2 or MT-ND2) have been used extensively to construct molecular phylogenies55–57. The Fish barcode of life (FISH-BOL) a working group of the International Barcode of Life Project (iBOL), has CO1 barcodes for 54% of elasmobranchs and 62% of holocephalans (http://www.fishbol.org, accessed July 24, 2014). Recently, whole mitochondrial sequences are increasingly popular for their increased granularity when resolving branches of phylogenetic trees1. Whole mitochondrial genome sequences currently are available for 72 species of sharks, skates, rays and chimaeras. These sequences are accessible in the GenBank, EMBL and DDBJ databases summarized in Table 558.
Chondrichthyan Tree of Life
Currently, molecular data for cartilaginous fishes is being collected as part of the Chondrichthyan Tree of Life project (http://sharksrays.org). The project website currently includes 5 elements: 1) an interactive phylogenetic tree55; 2) scientific illustrations of specimens; 3) range information for all extant species; 4) interactive comparative anatomy through segmented CT scan data; and 5) DNA sequence for 1265 single copy orthologous genes59. Project data will be available in public databases as well as through the project website once collection and analysis is complete (personal communication, Gavin Naylor, Medical University of South Carolina).
Protein databases
Given the improved technologies to characterize full-length transcripts using RNA-Seq, there are increasingly more protein sequence data for chondrichthyans. The UniProt Consortium, consists of groups from the European Bioinformatics Institute (EBI), the Swiss Institute of Bioinformatics (SIB) and the Protein Information Resource (PIR). The consortium maintains the UniProt Knowledgebase (UniProtKB), a comprehensive and standardized catalogue of protein sequences and functional annotation knowledgebase60. Proteins with UniProtKB accessions are first automatically annotated, unreviewed UniProtKB/TrEMBL entries that progress to UniProtKB/Swiss-Prot entries following curator review. Among Chondrichthyes, there are 12 UniProtKB/Swiss-Prot and 20,170 UniProtKB/TrEMBL entries for holocephalans and 264 UniProtKB/Swiss-Prot and 6,299 UniProtKB/TrEMBL entries for elasmobranchs in Release 2014_07 of 09-Jul-2014 of the knowledgebase (Table 1). An unidentified fin sample accounts for 16 UniProtKB/TrEMBL entries that are not included in either Holocephali or Elasmobranchii. PDB, an archive of protein macromolecular structural data, has 178 entries for Chondrichthyes, all elasmobranchs61. Of these, 76% are derived from 2 species from a single family, Torpediniformes, the electric rays, and in total only 10 species are represented in PDB.
The distribution of data in NCBI databases, PDB, and UniProtKB for chondrichthyan orders is shown in Figure 4. When order Chimaeriformes is included (Figure 4A) the distributions are disproportionate due to the large volume of annotated sequence data from the elephant shark genome. The distributions are repeated exclusively for elasmobranchs. To understand if the data distribution is representative of the number of species in each order, a species distribution is included in each chart. A cladogram (Figure 4B) is linked to the chart legend and illustrates the phylogeny between chondrichthyan orders.
Figure 4. A survey of public data and phylogeny for chondrichthyan orders.
A. The 14 orders of chondrichthyan fish and their relative distribution in public nucleotide and protein databases for Chondrichthyes and Elasmobranchii are shown individually. The species distribution for each Order and GenBank are similar indicating sequence data has been collected for a broad range of chondrichthyans. For Chondrichthyes, the elephant shark genome project data contributes the majority of the data in NCBI Gene, GenBank, Genome Survey Sequence (GSS), and the Short Reach Archive (SRA) databases. The NCBI GSS, GSS libraries, and Protein Data Bank (PDB) are the least diverse with representation of 1–6 of the 14 Orders. The color of each Order as represented in the bar chart is included in the cladogram key with left to right in the bar chart corresponding with top to bottom in the cladogram. B. A cladogram of Chondrichthyes illustrates the phylogeny relationship between the 14 Orders. The color code associated with each Order appears consecutively in the bar chart.
SkateBase
SkateBase (http://skatebase.org) is the public portal for the little skate genome project and is a valuable collection of data and learning resources. The NEBC little skate genome project team hosted three week-long workshops and a mitochondrial genome annotation jamboree with the goal of using the project data to develop a bioinformatics aware workforce and foster collaborative and distributed big data research. The lecture materials and worked annotation examples are included at SkateBase for educational use (http://skatebase.org/workshops). The project vitae contains an overview and timeline of the genome project effort along with key personnel, project related publications and presentations, the curation team, and citation information for researchers utilizing the resource in their publication. A Gene Table currently represents manually curated genes derived from workshops and curriculum with extensive annotation evidence. The number of gene entries will continue to grow through usage and expansion of the SkateBase educational modules. Plans to update the annotation interface to enable community annotation by domain experts is planned for the future.
SkateBase provides links to web resources with chondrichthyan data including the Chondrichthyan Tree of Life, Elephant Shark Genome Project (http://esharkgenome.imcb.a-star.edu.sg), the first described genome for a chimaera43, and Vertebrate TimeCapsule, (http://transcriptome.cdb.riken.go.jp/vtcap), a project that aims to develop a gene database to represent evolution and development for vertebrates and currently includes transcriptome data for a hagfish (Eptatretus burger), shark (S. torazame) and birchir (Polypterus senegalus)62. SkateBase data is linked locally as well as from NCBI in the Gene Expression Omnibus (GSE26235), GenBank (AESE010000000) and Sequence Read Archive (SRA026856) to ensure convenient and easy access. A link to the American Elasmobranch Society (http://www.elasmo.org), a non-profit organization with the mission of advancing the scientific study of living and fossil sharks, skates, rays, and chimaeras and promoting education, conservation, and wise utilization of natural resources, connects domain scientists to the little skate genome project.
SkateBase data includes embryonic transcriptomes for three chondrichthyan species, a chimaera, C. milii, a shark, S. canicula and the little skate, L. erinacea as well as the first draft of the little skate genome. The assembled skate genome sequence gave a single high-coverage contiguous sequence that represented the entire length of the mitochondrial genome. The mitochondrial genome was subsequently annotated as part of a Jamboree in 201163. The annotated sequence is represented by the NCBI Reference Sequence (RefSeq) project, accession NC_016429, and provides extensive information for each gene.
Whole embryos were used to build the transcriptome libraries available at SkateBase35. Two C. milii embryos, stage 32, were combined and used to build a chimaera library. The transcriptome library for S. canicula was assembled from six pooled embryos, stages 24–30. The embryonic skate transcriptome library was assembled using six pooled embryos ranging in stage from 20–29. This combination of stages encompasses a large portion of the developmental period for these fishes and represents a catalog of genes important for organogenesis of all or part of every physiological system. Early developmental events are similar for nearly all elasmobranchs regardless of reproductive mode or adult body form enabling the data to be useful for more than just the specific species from where it was derived64. Since all three embryonic transcriptomes contain a similar stage embryo direct comparison for temporal expression patterns is possible. Skatebase includes tools for data investigation, SkateBLAST, a sequence retrieval tool, Skate Contig Lookup, and genome browsers for three skate whole mitochondrial sequences, L. erinacea, the thorny skate, Amblyraja radiata, and, the ocellate spot skate, Okamejei kenojei. Skatebase contains resources that can be used for teaching and research purposes. As an example, two use cases follow, one for sequence or homology based research and the other for education.
SkateBLAST
A common task for researchers is searching for genes of interest in a genome or transcriptome. Knowledge of the gene sequence at the DNA or RNA level is needed for many different studies, including phylogenetic analysis or designing primers for quantitative PCR gene expression studies. Here we describe the major steps necessary to identify relevant sequences for a gene of interest using the BLAST sequence similarity tool at SkateBase. SkateBase features a web interface to BLAST, named SkateBLAST, that builds upon the ViroBLAST package version 2.265, with custom modifications allowing parallel cluster-based execution of queries and enhanced display of results. The overall workflow consists of a) entering a query sequence and selecting the database to search; b) evaluating the alignments returned; c) retrieving the sequence from one of the SkateBLAST databases; and d) checking to make sure that the retrieved sequence aligns best to the query sequence. The following description provides a brief tutorial on the overall workflow while describing tools at SkateBase.
Figure 5 demonstrates the use of SkateBLAST to find expressed sequences for the gene, suppressor of cytokine signaling 6 (SOCS6). SOCS6 is a E3 ubiquitin ligase that interacts with c-KIT to suppress cellular proliferation through its SH2 domain66. The first step to identify SOCS6 in the skate transcriptome begins with entering the protein sequence for human SOCS6 that was obtained from UniProt and searching this sequence against the skate transcriptome using the tblastn program. The next step is to evaluate the alignments to determine which transcriptome sequences best represent SOCS6. When interpreting the pairwise alignments from SkateBlast as in any BLAST tool, it is important to examine: a) alignment statistics; b) alignment coverage; and c) presence of protein domains that you may expect to be conserved. The alignment statistics are reported to ascertain whether you would expect the given alignment by chance or not. There are three key alignment statistics, the expectation (E)-value, percent identity and alignment length. The E-value represents the probability that you would expect an alignment with that alignment score or better by random chance, thus the lower the E-value, the better the alignment. Conversely, the greater the percent identity (percent identical sequence) and alignment length, the more similar the two sequences are assumed to be. Alignment coverage with respect to the query or subject sequence (alignment length divided by the length of the query or subject sequence) can also be an important consideration, as low coverage suggests that important regions of one or both sequences may not be represented in the alignment. Finally, there may be particular sequence features, such as protein domains, that you would expect to find in the alignment. If those domains are missing, then it suggests that you have a partial or misleading alignment.
Figure 5. Example of using SkateBase and NCBI resources to find transcriptome data for SOCS6.
A. SkateBLAST query form showing the four steps to align the UniProt sequence for human SOCS6 (O14544) against the skate embryonic transcriptome using tblastn. Step 1 is to enter the sequence in FASTA format. The second step is to choose the tblastn program that will align the query protein sequence against translated sequences in all six possible reading frames. The third step is to select the embryonic transcriptome as the sequence database to search. The fourth step is to launch the search. B. The complete BLAST output can be accessed by clicking the “Inspect BLAST output” link at the top of the summary report page. This is necessary to examine the sequence alignments. C. Four important fields in the output should be examined carefully to interpret the alignments and determine which returned alignment best represents the skate ortholog to SOCS6. First, the alignment score, E-value, alignment length and percent identity can be used to interpret the overall alignment significance. Alignment coverage with respect to the query protein sequence and the subject transcriptome sequence can be interpreted by comparing the alignment coordinates to the length of the query protein sequence and length of the transcriptome sequence. In this example, the entire query protein sequence is covered by this transcriptome sequence. D. The SkateBase Contig Lookup tool can be used to retrieve the transcriptome sequence found in the SOCS6 tblastn search in FASTA format. Sequences from the skate genome assembly or the skate, S. canicula or C. milii transcriptome assemblies can be retrieved using this tool. E. Output from the NCBI ORF Finder tool showing a 536aa ORF in the skate transcriptome contig that best represents SOCS6 (left). Alignment from blastx search of the skate transcriptome sequence (contig 15542) against human UniProt using NCBI BLAST to validate that the contig aligned best to human SOSC6 rather than another human gene.
Once a transcriptome sequence of interest, such as contig15542, is identified in the SkateBLAST results, you must do a reciprocal search of that sequence against a database of protein sequences to confirm that the sequence aligns best to your gene of interest. You can retrieve the full sequence directly from the BLAST tool or using the Skate Contig Lookup tool (Figure 5D): a) specify the transcriptome that you had originally searched using SkateBLAST; b) enter the sequence identification or contig number is entered into the query box; and c) select the ‘GO’ button. The user can copy the returned sequence and use it for further exploration of sequence homology at NCBI or similar databases.
SkateBase classroom use case: teach concepts of gene and protein annotation
SkateBase includes valuable teaching resources derived from the project workshops on gene and protein annotation. Infrastructure for sequence annotation was developed and modules for use in teaching are available. Access to the teaching modules is through the Curator Access link from the homepage and permission is granted by request using the email link at the bottom of each page, info@SkateBase.org. Once successfully logged into the site, access to pre-computed blast results, guides and examples, annotation forms, and links to external tools helpful for sequence analysis are available. Gene annotation begins with a transcriptome contig identified through a SkateBlast search as illustrated above. The portion of the transcript that codes for protein is identified using an open reading frame or ORF finder tool. Annotation follows a workflow where complimentary sequences from the transcriptome and genome are aligned allowing annotation of both sequences using Sequence Ontology vocabulary67. The evidence is recorded in an annotation form that records information about the annotator and sequences and includes a comment box for questions and comments between students and teachers or curators and annotators. The annotation form records the pairwise alignment of the transcriptome and genome contigs, notes concerning mismatches or gaps, as well as output from the ORF tool. The untranslated regions (UTR) at the beginning and end of each sequence, 5’UTR and 3’ UTR regions, as well as the intron/exon structure for the genomic contig and CDS for the transcriptomic contig are recorded in the Gene Annotation Form. When completing the Gene Annotation Form, the appropriate activity must be selected and can be customized to specify the user’s course ID, institution or workshop title to track annotation history. Protein annotation uses the rapid annotation interface for proteins, RACE-P, developed by the PIR. A UniProt accession number is required to initiate a new annotation form. The form is composed of 6 blocks of information, protein information, gene information, a bibliography, Gene Ontology (GO), computational analysis using tools such as Pfam68, TMHMM69, SignalP70, COILS71, NetPhos72 and EMBOSS73, and protein family evidence.
Discussion
The volume of data in GenBank continues to grow exponentially, doubling nearly every 18 months. The first sequences for chondrichthyes appeared in 1983 and the overall data trend for chondrichthyans is similar to all of GenBank with three exceptions. First, the rate of increase is less than GenBank. Second, the number of sequences deposited during the first decade of the 21st century was nearly stagnant in comparison. Third, a large spike is observed in late 2012 attributed to the Elephant Shark Genome Project data (Figure 6). Molecular data is increasingly important for all aspects of research utilizing chondrichthyan fishes74. It can be a forensic tool to understand species when fins are landed without carcass and ensure protected species and quotas are respected75–79. For migrating species molecular data serves as a surrogate to classical tagging data to understand population structure and range80–83. In studies of evolution, molecular data provides estimates of divergence time and supplements morphological and ecological traits as the basis for a phylogeny. The benefits and uses of molecular data for these fishes are limited only by the amount of data available. SkateBase provides the only genomic data publically available for an elasmobranch in addition to embryonic transcriptomes, data tools, and educational resources.
Figure 6. GenBank and WGS data trends for Chondrichthyes and all taxa.
GenBank is the National Institutes of Health (NIH) genetic sequence database and together with the DNA Databank of Japan (DDBJ) and the European Molecular Biological Laboratory (EMBL) comprise the International Nucleotide Sequence Database Collaboration (INSDC). The cumulative base pair total for all taxa as well as chondrichthyan only data are given versus time for GenBank and Whole Genome Shotgun (WGS) data. The Elephant Shark Genome Project is responsible for the spike in chondrichthyan GenBank in 2011. The little skate and elephant shark genome projects are currently the only two WGS datasets (yellow line).
Sequencing projects require significant funding and personnel commitments but generate a large amount of information that can be translated to knowledge by domain experts. The efficiency of this process is affected most by allowing the scientific community to access the data. The value of data sharing can be measured by the number of publications that result from its distribution. To date, 19 publications in peer-reviewed journals have used data derived from SkateBase (http://skatebase.org/vitae). Molecular data are the means to investigate genes and develop reagents for gene expression studies by PCR or in situ hybridization. Small scale sequencing efforts that generate limited or fragmented data often get deposited to hard disks and remain ‘buried’ and out of reach. Efforts to deposit this data at public sequence repositories are encouraged to build the foundation of data required to describe this dynamic and ancient clade of fishes. We invite investigators to contact the authors in an effort to survey the volume of private data available for potential distribution through SkateBase.
The transcriptome data at SkateBase serves as a platform to teach molecular techniques, technologies, and bioinformatics in the context of studying elasmobranchs. As next generation sequencing (NGS) technologies evolve it is important for scientists and students to understand how the sequence was generated and caveats of workflow for each data type in order to recognize errors and customize analysis algorithms. The educational materials and infrastructure at SkateBase have been used by University of Delaware, Georgetown University, MDI Biological Laboratory, University of Maine at Machias, University of Rhode Island, and most recently the Virginia Institute of Marine Science to teach gene and protein annotation concepts. We invite and look forward to continued expansion of the SkateBase educational platform as we refine the infrastructure and expand the data available for investigation through continued sequencing efforts.
Author contributions
The little skate genome project is a collaborative effort with participation from all authors. JTW and BLK prepared the first draft of the manuscript. SWP oversaw the design and implementation of the SkateBase web portal and SkateBLAST tool. All authors were involved with the revision of the draft manuscript and have agreed to the final content.
Competing interests
The authors declare to have no competing interests.
Grant information
This work was supported by a re-entry career award to JTW, National Institute of General Medical Sciences (NIGMS) IDeA Networks of Biomedical Research Excellence (INBRE) 3P20GM103446-12S1. Skate genome sequencing was funded by National Institutes of Health (NIH) National Center for Research Resources (NCRR) ARRA Supplements to 5P20RR016463-12 (MDIBL), 5P20RR016472-12 (UD), and 5P20RR16462 (UVM). The North East Cyberinfrastructure Consortium (NECC) is funded by NIH NCRR grants 5P20RR016463-12 (MDIBL), 5P20RR016472-12 (UD), 5 P20 RR16462 (UVM), 5P20RR016457-11 (URI), and 5P20RR030360-03 (UNH) and NIH NIGMS grants 8 P20 GM103423-12 (MDIBL), 8P20GM103446-12 (UD), 8P20GM103449 (UVM), 8 P20 GM103430-11 (URI), and 8P20GM103506-03 (Dartmouth), and NSF Experimental Program to Stimulate Competitive Research (EPSCoR) grants EPS-0904155 (UM), EPS-081425 (UD), EPS-1101317 (UVM), EPS-1004057 (URI), and EPS-1101245 (UNH).
Acknowledgements
The authors thank the North East Bioinformatics Collaborative of the North East Cyberinfrastructure Consortium and Karl Steiner and Steven Stanhope of DE-INBRE, Patricia Hand of ME-INBRE, and Judith VanHouten of VT-INBRE. We thank the Delaware Biotechnology Institute (DBI) for hosting the website, database and NECC Shared Data Center. We thank Karol Miaskiewicz for systems administration support and Gang Li and Zhiwen Li for web development. We thank Bruce Kingham and the DBI Sequencing and Genotyping Center. Sylvie Mazan, Alistair Dove, Dave Ebert, Cecelia Arighi, Qinghua Wang, James Sulikowski and Eric Haenni contributed data, samples, artwork or editorial comments. We thank Gavin Naylor and the Chondrichthyan Tree of Life project team including Lindsay Marshall (illustrations), Jason Davies (database, computational work, and visualizations), Will White and Peter Last (maps and taxonomy), Shannon Corrigan and Lei Yang (gene capture data), and Callie Crawford and Thomas Fussell (CT scanning and anatomy) for providing a description of the project scope.
Faculty Opinions recommendedReferences
- 1.
Inoue JG, Miya M, Lam K, et al.:
Evolutionary origin and phylogeny of the modern holocephalans (Chondrichthyes: Chimaeriformes): a mitogenomic perspective.
Mol Biol Evol.
2010; 27(11): 2576–86. PubMed Abstract
| Publisher Full Text
- 2.
Ebert DA, Ho H, White WT, et al.:
Introduction to the systematics and biodiversity of sharks, rays, and chimaeras (Chondrichthyes) of Taiwan. 2013; 3752(1): 5–19. Publisher Full Text
- 3.
Stevens J, Bonfil R, Dulvy NK, et al.:
The effects of fishing on sharks, rays, and chimaeras (chondrichthyans), and the implications for marine ecosystems.
ICES J Mar Sci.
2000; 57(3): 476–94. Publisher Full Text
- 4.
García VB, Lucifora LO, Myers RA:
The importance of habitat and life history to extinction risk in sharks, skates, rays and chimaeras.
Proc Biol Sci.
2008; 275(1630): 83–9. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 5.
Field IC, Meekan MG, Buckworth RC, et al.:
Chapter 4. Susceptibility of sharks, rays and chimaeras to global extinction.
Adv Mar Biol.
2009; 56: 275–363. PubMed Abstract
| Publisher Full Text
- 6.
Springer S:
Oviphagous Embryos of the Sand Shark, Carcharias taurus.
Copeia.
1948; 1948(3): 153. Reference Source
- 7.
Joung SJ, Chen CT, Clark E, et al.:
The whale shark, Rhincodon typus, is a livebearer: 300 embryos found in one ‘megamamma’ supreme.
Environ Biol Fishes.
1996; 46(3): 219–23. Publisher Full Text
- 8.
Stevens J, Bonfil R, Dulvy NK, et al.:
The effects of fishing on sharks, rays, and chimaeras (chondrichthyans), and the implications for marine ecosystems.
ICES J Mar Sci.
2000; 57(3): 476–94. Publisher Full Text
- 9.
Liu SY, Chan CL, Lin O, et al.:
DNA barcoding of shark meats identify species composition and CITES-listed species from the markets in Taiwan.
PLoS One.
2013; 8(11): e79373. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 10.
Dudgeon CL, Blower DC, Broderick D, et al.:
A review of the application of molecular genetics for fisheries management and conservation of sharks and rays.
J Fish Biol.
2012; 80(5): 1789–843. PubMed Abstract
| Publisher Full Text
- 11.
Pinhal D, Shivji MS, Nachtigall PG, et al.:
A streamlined DNA tool for global identification of heavily exploited coastal shark species (genus Rhizoprionodon).
PLoS One.
2012; 7(4): e34797. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 12.
Barbuto M, Galimberti A, Ferri E, et al.:
DNA barcoding reveals fraudulent substitutions in shark seafood products: The Italian case of “palombo” (Mustelus spp.).
Food Res Int.
2010; 43(1): 376–81. Publisher Full Text
- 13.
Feldheim KA, Gruber SH, Dibattista JD, et al.:
Two decades of genetic profiling yields first evidence of natal philopatry and long-term fidelity to parturition sites in sharks.
Mol Ecol.
2014; 23(1): 110–7. PubMed Abstract
| Publisher Full Text
- 14.
Schneider I, Aneas I, Gehrke AR, et al.:
Appendage expression driven by the Hoxd Global Control Region is an ancient gnathostome feature.
Proc Natl Acad Sci U S A.
2011; 108(31): 12782–6. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 15.
Gillis JA, Dahn RD, Shubin NH:
Shared developmental mechanisms pattern the vertebrate gill arch and paired fin skeletons.
Proc Natl Acad Sci U S A.
2009; 106(14): 5720–4. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 16.
Gillis JA, Modrell MS, Baker CVH:
Developmental evidence for serial homology of the vertebrate jaw and gill arch skeleton.
Nat Commun.
2013; ; 4: 1436. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 17.
Hentschel DM, Bonventre JV:
Novel non-rodent models of kidney disease.
Curr Mol Med.
2005; 5(5): 537–46. PubMed Abstract
| Publisher Full Text
- 18.
Steele SL, Yancey PH, Wright PA:
Dogmas and controversies in the handling of nitrogenous wastes: osmoregulation during early embryonic development in the marine little skate Raja erinacea; response to changes in external salinity.
J Exp Biol.
2004; 207(Pt 12): 2021–31. PubMed Abstract
| Publisher Full Text
- 19.
Stolte H, Galaske RG, Eisenbach GM, et al.:
Renal tubule ion transport and collecting duct function in the elasmobranch little skate, Raja erinacea.
J Exp Zool.
1977; 199(3): 403–10. PubMed Abstract
| Publisher Full Text
- 20.
Evans DH:
A brief history of the study of fish osmoregulation: the central role of the Mt. Desert Island Biological Laboratory.
Front Physiol.
2010; 1: 13. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 21.
Wyffels JT, Luer CA, Walsh CJ, et al.:
In vivo exposure of clearnose skates, Raja eglanteria, to ionising X-radiation: acute effects on the peripheral blood, spleen, and epigonal and Leydig organs.
Fish Shellfish Immunol.
2007; 23(2): 401–18. PubMed Abstract
| Publisher Full Text
- 22.
Wyffels JT, Walsh CJ, Luer CA, et al.:
In vivo exposure of clearnose skates, Raja eglanteria, to ionizing X-radiation: acute effects on the thymus.
Dev Comp Immunol.
2005; 29(4): 315–31. PubMed Abstract
| Publisher Full Text
- 23.
Lutton BV, Callard IP:
Morphological relationships and leukocyte influence on steroid production in the epigonal organ-ovary complex of the skate, Leucoraja erinacea.
J Morphol.
2008; 269(5): 620–9. PubMed Abstract
| Publisher Full Text
- 24.
Walsh CJ, Luer CA, Bodine AB, et al.:
Elasmobranch immune cells as a source of novel tumor cell inhibitors: Implications for public health.
Integr Comp Biol.
2006; 46(6): 1072–81. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 25.
Eason DD, Litman RT, Luer CA, et al.:
Expression of individual immunoglobulin genes occurs in an unusual system consisting of multiple independent loci.
Eur J Immunol.
2004; 34(9): 2551–8. PubMed Abstract
| Publisher Full Text
- 26.
Anderson MK, Pant R, Miracle AL, et al.:
Evolutionary origins of lymphocytes: ensembles of T cell and B cell transcriptional regulators in a cartilaginous fish.
J Immunol.
2004; 172(10): 5851–60. PubMed Abstract
| Publisher Full Text
- 27.
Cai SY, Soroka CJ, Ballatori N, et al.:
Molecular characterization of a multidrug resistance-associated protein, Mrp2, from the little skate.
Am J Physiol Regul Integr Comp Physiol.
2003; 284(1): R125–30. PubMed Abstract
- 28.
Clusin W, Spray DC, Bennett MVL:
Activation of a voltage-insensitive conductance by inward calcium current.
Nature.
1975; 256(5516): 425–7. PubMed Abstract
| Publisher Full Text
- 29.
Elger M, Hentschel H, Litteral J, et al.:
Nephrogenesis is induced by partial nephrectomy in the elasmobranch Leucoraja erinacea.
J Am Soc Nephrol.
2003; 14(6): 1506–18. PubMed Abstract
| Publisher Full Text
- 30.
Zhang C, Du Pasquier L, Hsu E:
Shark IgW C region diversification through RNA processing and isotype switching.
J Immunol.
2013; 191(6): 3410–8. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 31.
Beard J:
The origin and histogenesis of the thymus in Raja batis.
Zool Jahrb Abt Anat Ontog Tiere.
1902; 17: 403–80. Reference Source
- 32.
Luer CA, Walsh CJ, Bodine AB, et al.:
The elasmobranch thymus: Anatomical, histological, and preliminary functional characterization.
J Exp Zool.
1995; 273(4): 342–54. Publisher Full Text
- 33.
Wyffels JT, Walsh CJ, Luer CA, et al.:
In vivo exposure of clearnose skates, Raja eglanteria, to ionizing X-radiation: acute effects on the thymus.
Dev Comp Immunol.
2005; 29(4): 315–31. PubMed Abstract
| Publisher Full Text
- 34.
Feldheim KA, Chapman DD, Sweet D, et al.:
Shark virgin birth produces multiple, viable offspring.
J Hered.
2010; 101(3): 374–7. PubMed Abstract
| Publisher Full Text
- 35.
Chapman DD, Shivji MS, Louis E, et al.:
Virgin birth in a hammerhead shark.
Biol Lett.
2007; 3(4): 425–7. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 36.
Robinson DP, Baverstock W, Al-Jaru A, et al.:
Annually recurring parthenogenesis in a zebra shark Stegostoma fasciatum.
J Fish Biol.
2011; 79(5): 1376–82. PubMed Abstract
| Publisher Full Text
- 37.
Compagno LJV:
Alternative life-history styles of cartilaginous fishes in time and space.
Environ Biol Fishes.
1990; 28(1–4): 33–75. Publisher Full Text
- 38.
Luer CA, Gilbert PW:
Mating behavior, egg deposition, incubation period, and hatching in the clearnose skate, Raja eglanteria.
Environ Biol Fishes.
1985; 13(3): 161–71. Publisher Full Text
- 39.
Luer CA, Walsh CJ, Bodine AB, et al.:
Normal embryonic development in the clearnose skate, Raja eglanteria, with experimental observations on artificial insemination.
Environ Biol Fishes.
2007; 80(2–3): 239–55. Publisher Full Text
- 40.
Masuda M, Izawa Y, Kametuta S, et al.:
Artificial insemination of the cloudy catshark.
J Japanese Assoc Zool Gard Aquariums.
2003; 44(2): 39–43. Reference Source
- 41.
Castro JI, Bubucis PM, Overstrom NA, et al.:
The Reproductive Biology of the Chain Dogfish, Scyliorhinus retifer.
Copeia.
1988; 1988(3): 740. Publisher Full Text
- 42.
Barrett T, Clark K, Gevorgyan R, et al.:
BioProject and BioSample databases at NCBI: facilitating capture and organization of metadata.
Nucleic Acids Res.
2012; 40(Database issue): D57–63. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 43.
Venkatesh B, Lee AP, Ravi V, et al.:
Elephant shark genome provides unique insights into gnathostome evolution.
Nature.
2014; 505(7482): 174–9. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 44.
Venkatesh B, Kirkness EF, Loh YH, et al.:
Survey sequencing and comparative analysis of the elephant shark (Callorhinchus milii) genome.
PLoS Biol.
2007; 5(4): e101. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 45.
Luo M, Kim H, Kudrna D, et al.:
Construction of a nurse shark (Ginglymostoma cirratum) bacterial artificial chromosome (BAC) library and a preliminary genome survey.
BMC Genomics.
2006; 7: 106. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 46.
Gwee PC, Tay BH, Brenner S, et al.:
Characterization of the neurohypophysial hormone gene loci in elephant shark and the Japanese lamprey: origin of the vertebrate neurohypophysial hormone genes.
BMC Evol Biol.
2009; 9: 47. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 47.
Mulley JF, Zhong YF, Holland PW:
Comparative genomics of chondrichthyan Hoxa clusters.
BMC Evol Biol.
2009; 9: 218. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 48.
Chiu C, Amemiya C, Dewar K, et al.:
Molecular evolution of the HoxA cluster in the three major gnathostome lineages.
Proc Natl Acad Sci U S A.
2002; 99(8): 5492–7. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 49.
Mattingly C, Parton A, Dowell L, et al.:
Cell and molecular biology of marine elasmobranchs: Squalus acanthias and Raja erinacea.
Zebrafish.
2004; 1(2): 111–20. PubMed Abstract
| Publisher Full Text
- 50.
Merson RR, Mattingly C, Planchart AJ:
Tandem duplication of aryl hydrocarbon receptor (AHR) genes in the genome of the spiny dogfish shark (Squalus acanthias).
Bull MDIBL.
2009; 48: 43–4. Reference Source
- 51.
Oulion S, Debiais-Thibaud M, d’Aubenton-Carafa Y, et al.:
Evolution of Hox gene clusters in gnathostomes: insights from a survey of a shark (Scyliorhinus canicula) transcriptome.
Mol Biol Evol.
2010; 27(12): 2829–38. PubMed Abstract
| Publisher Full Text
- 52.
King BL, Gillis JA, Carlisle HR, et al.:
A natural deletion of the HoxC cluster in elasmobranch fishes.
Science.
2011; 334(6062): 1517. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 53.
Kim CB, Amemiya C, Bailey W, et al.:
Hox cluster genomics in the horn shark, Heterodontus francisci.
Proc Natl Acad Sci U S A.
2000; 97(4): 1655–60. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 54.
Ravi V, Lam K, Tay BH, et al.:
Elephant shark (Callorhinchus milii) provides insights into the evolution of Hox gene clusters in gnathostomes.
Proc Natl Acad Sci U S A.
2009; 106(38): 16327–32. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 55.
Naylor GJP, Caira JN, Jensen K, et al.:
Elasmobranch Phylogeny: A mitochondrial estimate based on 595 species. In: Carrier JC, Musick JA, Heithaus MR, editors. The Biology of Sharks and Their Relatives. CRC Press, Taylor & Francis Group. 2012; 31–56. Reference Source
- 56.
Naylor GJP, Ryburn JA, Fedrigo O, et al.:
Phylogenetic relationships among the Major Lineages of Modern Elasmobranchs. 1995. Reference Source
- 57.
Naylor GJP, Caira JN, Jensen K, et al.:
A DNA Sequence–Based Approach To the Identification of Shark and Ray Species and Its Implications for Global Elasmobranch Diversity and Parasitology.
Bull Am Museum Nat Hist.
2012; 367: 1–262. Publisher Full Text
- 58.
Feijão PC, Neiva LS, de Azeredo-Espin AM, et al.:
AMiGA: the arthropodan mitochondrial genomes accessible database.
Bioinformatics.
2006; 22(7): 902–3. PubMed Abstract
| Publisher Full Text
- 59.
Li C, Hofreiter M, Straube N, et al.:
Capturing protein-coding genes across highly divergent species.
Biotechniques.
2013; 54(6): 321–6. PubMed Abstract
- 60.
UniProt Consortium. Activities at the Universal Protein Resource (UniProt).
Nucleic Acids Res.
2014; 42(Database issue): D191–8. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 61.
Berman HM, Westbrook J, Feng Z, et al.:
The Protein Data Bank.
Nucleic Acids Res.
2000; 28(1): 235–42. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 62.
Takechi M, Takeuchi M, Ota KG, et al.:
Overview of the transcriptome profiles identified in hagfish, shark, and bichir: current issues arising from some nonmodel vertebrate taxa.
J Exp Zool B Mol Dev Evol.
2011; 316(7): 526–46. PubMed Abstract
| Publisher Full Text
- 63.
Wang Q, Arighi CN, King BL, et al.:
Community annotation and bioinformatics workforce development in concert--Little Skate Genome Annotation Workshops and Jamborees.
Database (Oxford).
2012; 2012: bar064. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 64.
Jennifer TW:
Embryonic development of Chondrichthyan fishes – a review. In: W KY, Lue CA & Kapoor BG, editors. Development of Non-Teleost Fishes. Enfield: Science Publishers. 2009; 1–103. Publisher Full Text
- 65.
Deng W, Nickle DC, Learn GH, et al.:
ViroBLAST: a stand-alone BLAST web server for flexible queries of multiple databases and user’s datasets.
Bioinformatics.
2007; 23(17): 2334–6. PubMed Abstract
| Publisher Full Text
- 66.
Zadjali F, Pike AC, Vesterlund M, et al.:
Structural basis for c-KIT inhibition by the suppressor of cytokine signaling 6 (SOCS6) ubiquitin ligase.
J Biol Chem.
2011; 286(1): 480–90. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 67.
Mungall CJ, Batchelor C, Eilbeck K:
Evolution of the Sequence Ontology terms and relationships.
J Biomed Inform.
2011; 44(1): 87–93. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 68.
Finn RD, Bateman A, Clements J, et al.:
Pfam: the protein families database.
Nucleic Acids Res.
2014; 42(Database issue): D222–30. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 69.
Krogh A, Larsson B, von Heijne G, et al.:
Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes.
J Mol Biol.
2001; 305(3): 567–80. PubMed Abstract
| Publisher Full Text
- 70.
Petersen TN, Brunak S, von Heijne G, et al.:
SignalP 4.0: discriminating signal peptides from transmembrane regions.
Nat Methods.
2011; 8(10): 785–6. PubMed Abstract
| Publisher Full Text
- 71.
Lupas A, Van Dyke M, Stock J:
Predicting coiled coils from protein sequences.
Science.
1991; 252(5009): 1162–4. PubMed Abstract
| Publisher Full Text
- 72.
Blom N, Gammeltoft S, Brunak S:
Sequence and structure-based prediction of eukaryotic protein phosphorylation sites.
J Mol Biol.
1999; 294(5): 1351–62. PubMed Abstract
| Publisher Full Text
- 73.
Rice P, Longden I, Bleasby A:
EMBOSS: the European Molecular Biology Open Software Suite.
Trends Genet.
2000; 16(6): 276–7. PubMed Abstract
| Publisher Full Text
- 74.
Heist E:
Genetics of Sharks, Skates, and Rays. In: Musick J, Carrier J, Heithaus M, editors. Biology of Sharks and Their Relatives. CRC Press. 2004; 471–85. Publisher Full Text
- 75.
Chapman D, Pinhal D, Shivji M:
Tracking the fin trade: genetic stock identification in western Atlantic scalloped hammerhead sharks Sphyrna lewini.
Endanger Species Res.
2009; 9: 221–8. Publisher Full Text
- 76.
Clarke SC, Magnussen JE, Abercrombie DL, et al.:
Identification of shark species composition and proportion in the Hong Kong shark fin market based on molecular genetics and trade records.
Conserv Biol.
2006; 20(1): 201–11. PubMed Abstract
| Publisher Full Text
- 77.
Liu SY, Chan CL, Lin O, et al.:
DNA barcoding of shark meats identify species composition and CITES-listed species from the markets in Taiwan.
PLoS One.
2013; 8(11): e79373. PubMed Abstract
| Publisher Full Text
| Free Full Text
- 78.
Shivji MS:
DNA Forensic Applications in Shark Management and Conservation. In: Carrier JC., Musick JA., Heithaus MR., editors. Sharks and Their Relatives II Biodiversity, Adaptive Physiology, and Conservation. CRC Press. 2010; 593–610. Publisher Full Text
- 79.
Rodrigues-Filho LF, Pinhal D, Sodré D, et al.:
Applications and Impacts on Genetic Diversity. In: Caliskan M, editor. Analysis of Genetic Variation in Animals. InTech. 2012; 269–86. Publisher Full Text
- 80.
Chabot CL, Allen LG:
Global population structure of the tope (Galeorhinus galeus) inferred by mitochondrial control region sequence data.
Mol Ecol.
2009; 18(3): 545–52. PubMed Abstract
| Publisher Full Text
- 81.
Keeney DB, Heist EJ:
Worldwide phylogeography of the blacktip shark (Carcharhinus limbatus) inferred from mitochondrial DNA reveals isolation of western Atlantic populations coupled with recent Pacific dispersal.
Mol Ecol.
2006; 15(12): 3669–79. PubMed Abstract
| Publisher Full Text
- 82.
Schultz JK, Feldheim KA, Gruber SH, et al.:
Global phylogeography and seascape genetics of the lemon sharks (genus Negaprion).
Mol Ecol.
2008; 17(24): 5336–48. PubMed Abstract
| Publisher Full Text
- 83.
Ward RD, Holmes BH, White WT, et al.:
DNA barcoding Australasian chondrichthyans: results and potential uses in conservation.
Mar Freshw Res.
2008; 59(1): 57–71. Publisher Full Text
Comments on this article Comments (0)