The complete mitochondrial genome sequence of the world's largest fish, the whale shark (Rhincodon typus), and its comparison with those of related shark species

Characterization of the complete mitochondrial genome of the brazilian cownose ray Rhinoptera brasiliensis (Myliobatiformes, Rhinopteridae) in the western Atlantic and its phylogenetic implications

BACKGROUND

The Brazilian cownose ray, Rhinoptera brasiliensis has undergone a global population reduction and is currently classified by IUCN as Vulnerable. This species is sometimes confused with Rhinoptera bonasus, the only external diagnostic characteristic to distinguish between both species is the number of rows of tooth plates. Both cownose rays overlap geographically from Rio de Janeiro to the western North Atlantic. This calls for a more comprehensive phylogenetic assessment using mitochondria DNA genomes to better understand the relationships and delimitation of these two species.

METHODS AND RESULTS

The mitochondrial genome sequences of R. brasiliensis was obtained by next-generation sequencing. The length of the mitochondrial genome was 17,759 bp containing 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a non-coding control region (D-loop). Each PCG was initiated by an authoritative ATG codon, except for COX1 initiated by a GTG codon. Most of the PCGs were terminated by a complete codon (TAA/TAG), while an incomplete termination codon (TA/T) was found in five out of the 13 PCGs. The phylogenetic analysis showed that R. brasiliensis was closely related to R. steindachneri whereas the reported mitogenome as R. steindachneri (GenBank accession number KM364982), differs from multiple mitocondrial DNA sequences of R. steindachneri and is nearly identical to that of R. javanica.

CONCLUSION

The new mitogenome determined in this study provides new insight into the phylogenetic relationships in Rhinoptera, while providing new molecular data that can be applied to population genetic studies.

Characterization of mitochondrial genome of Indian Ocean blue-spotted maskray, Neotrygon indica and its phylogenetic relationship within Dasyatidae Family

The present study characterized complete mitochondrial genome of Blue-spotted maskray, Neotrygon indica and studied the evolutionary relationship of the species within the Dasyatidae family. The total length of the mitogenome was 17,974 bp including 37 genes and a non-coding control region. The average frequency of nucleotides in protein-coding genes was A: 29.1 %, T: 30.2 %, G: 13.0 % and C: 27.7 % with AT content of 59.3 %. The values of AT and GC skewness were -0.018 and -0.338, respectively. Comparative analyses showed a large number of average synonymous substitutions per synonymous site (Ks) in gene NADH4 (5.07) followed by NADH5 (4.72). High values of average number of non-synonymous substitutions per non-synonymous site (Ka) were observed in genes ATPase8 (0.54) and NADH2 (0.44). Genes NADH4L and NADH2 showed high interspecific genetic distance values of 0.224 ± 0.001 and 0.213 ± 0.002, respectively. Heat map analysis showed variation in codon usage among different species of the Dasyatidae family. The phylogenetic tree showed a sister relationship between the Dasyatinae and the Neotrygoninae subfamilies. Neotrygon indica formed as a sister species to the clade consisting of N. varidens and N. orientalis. Based on the present results, Neotrygon indica could have diverged from the common ancestor of the two latter in the Plio-Pleistocene. The present study showed distinct characteristics of N. indica from its congeners through comparative mitogenomics.

Pieces in a global puzzle: Population genetics at two whale shark aggregations in the western Indian Ocean

The whale shark Rhincodon typus is found throughout the world's tropical and warm-temperate ocean basins. Despite their broad physical distribution, research on the species has been concentrated at a few aggregation sites. Comparing DNA sequences from sharks at different sites can provide a demographically neutral understanding of the whale shark's global ecology. Here, we created genetic profiles for 84 whale sharks from the Saudi Arabian Red Sea and 72 individuals from the coast of Tanzania using a combination of microsatellite and mitochondrial sequences. These two sites, separated by approximately 4500 km (shortest over-water distance), exhibit markedly different population demographics and behavioral ecologies. Eleven microsatellite DNA markers revealed that the two aggregation sites have similar levels of allelic richness and appear to be derived from the same source population. We sequenced the mitochondrial control region to produce multiple global haplotype networks (based on different alignment methodologies) that were broadly similar to each other in terms of population structure but suggested different demographic histories. Data from both microsatellite and mitochondrial markers demonstrated the stability of genetic diversity within the Saudi Arabian aggregation site throughout the sampling period. These results contrast previously measured declines in diversity at Ningaloo Reef, Western Australia. Mapping the geographic distribution of whale shark lineages provides insight into the species' connectivity and can be used to direct management efforts at both local and global scales. Similarly, understanding historical fluctuations in whale shark abundance provides a baseline by which to assess current trends. Continued development of new sequencing methods and the incorporation of genomic data could lead to considerable advances in the scientific understanding of whale shark population ecology and corresponding improvements to conservation policy.

The complete mitochondrial genome of the endemic Iberian pygmy skate Neoraja iberica Stehmann, Séret, Costa, & Baro 2008 (Elasmobranchii, Rajidae)

Skates, Chondrichthyes fishes from order Rajiformes, are the most species-rich group of all Batoidea. However, their phylogenetic relationships and systematics is still a highly discussed and controversial subject. The use of complete mitogenome has shown to be a promising tool to fill this gap of knowledge. Here, the complete mitogenome of the Iberian pygmy skate Neoraja iberica (Stehmann, Séret, Costa & Baro 2008) was sequenced and assembled. The mitogenome is 16,723 bp long and its gene content (i.e. 13 protein-coding genes, 22 transfer RNA, and 2 ribosomal RNA genes) and arrangement are the expected for Batoidea. Phylogenetic reconstructions, including 89 Rajiformes and two outgroup Rhinopristiformes, recovered family Rajidae as monophyletic, and further divided in the monophyletic tribe Rajini, sister to tribes Amblyrajini and Rostrorajini. The newly sequenced N. iberica mitogenome is the first representative of the tribe Rostrorajini.

Mitochondrial Genome Structures and Phylogenetic Analyses of Two Tropical Characidae Fishes

The Characidae family contains the largest number of tropical fish species. Morphological similarities make species identification difficult within this family. Here, the complete mitogenomes of two Characidae fish were determined and comparatively analyzed with those of nine other Characidae fish species. The two newly sequenced complete mitogenomes are circular DNA molecules with sizes of 16,701 bp (Hyphessobrycon amandae; MT484069) and 16,710 bp (Hemigrammus erythrozonus; MT484070); both have a highly conserved structure typical of Characidae, with the start codon ATN (ATG/ATT) and stop codon TAR (TAA/TAG) or an incomplete T--/TA-. Most protein-coding genes of the 11 Characidae mitogenomes showed significant codon usage bias, and the protein-coding gene cox1 was found to be a comparatively slow-evolving gene. Phylogenetic analyses via the maximum likelihood and Bayesian inference methods confirmed that H. amandae and H. erythrozonus belong to the family Characidae. In all Characidae species studied, one genus was well supported; whereas other two genera showed marked differentiation. These findings provide a phylogenetic basis for improved classification of the family Characidae. Determining the mitogenomes of H. erythrozonus and H. amandae improves our understanding of the phylogeny and evolution of fish species.

Chromosome-level genome assembly of Lethenteron reissneri provides insights into lamprey evolution

The reissner lamprey Lethenteron reissneri, belonging to the class Cyclostomata, serves as a bridge between invertebrates and jawed vertebrates, and is considered the sister group of jawed vertebrates. However, despite this evolutionary significance, the genetic mechanisms underlying the adaptive evolution of the lamprey lineage remain unclear. Here, we assembled a 1.06 Gb chromosome-level draft genome of L. reissneri, with 72 chromosomes (ranging in length from 4.5 Mb to 25.9 Mb) and a scaffold N50 length of 13.23 Mb. Genome quality comparisons revealed that the reissner lamprey genome has higher completeness and contiguity than the previously published sea lamprey and Japanese lamprey genomes. Moreover, reissner lamprey, sea lamprey, and Japanese lamprey species share similar transposable element profiles and Hox gene cluster compositions, suggesting that a burst of transposable element activity and whole genome duplication occurred before their divergence. Additionally, the Lip gene copy numbers, which have been studied for their functions in the host defence system, were found to be expanded uniquely in lamprey lineages, suggesting key roles for these genes in lamprey evolution and adaptation. We also identified two neural-related genes, Nrn1 and Unc13a, with copy number expansions in jawed vertebrates, which may be functionally relevant to the origin of lamprey brains. Hence, this study not only provides the first chromosome-level reference genome for Cyclostomata, but also highlights features of the unique biology and adaptive evolution of the lamprey lineage.

Complete mitochondrial genome and assembled DNA barcoding analysis of Lutjanus fulgens (Valenciennes, 1830) and its comparison with other Lutjanus species

Lutjanus fulgens (Valenciennes, 1830) is a teleost species classified under the family Lutjanidae which is a native of the Eastern Atlantic Ocean. Though highly commercialized due to its abundance and good taste, the production output has declined in recent years. This is an indication of the need for effective management and conservation measures. However, accurate species identification will ensure strategic management and conservation measure. DNA-based species identification has proven its reliability in this regard via precise species identification. Several researchers have confirmed the accuracy of DNAbarcode as a species identification tool as well as species phylogeny analysis based on both the complete mitogenome and COI gene. Currently, nine specimens of L. fulgens were sampled from Ghana and subjected to DNA-based analysis, namely, complete mitochondrial DNAand COI gene (DNA barcoding) analyses. The mitogenomic result revealed that L. fulgens is made up of a 16,500 base pairs (bp) mtDNA which consists of 22 transfer RNAs, 13 protein-coding genes, and two ribosomal RNAs (GenBank Accession Number: MN398650). Furthermore, a sequence polymorphism analysis of the COIgene (MN986442-MN986450) detected two haplotypes. These haplotypes were both collected from the same fish landing site which suggests a possible cryptic linage diversity in the L. fulgens population at Vodza. According to the phylogeny examination, a close taxonomic relationship exists between L. fulgens and Lutjanus buccanella caused by a recent evolution termed as sympatric speciation. This study serves as a novel study for this species, building the foundation for future molecular-based study for this species and as a DNA barcode reference data.

Rapid detection of CITES-listed shark fin species by loop-mediated isothermal amplification assay with potential for field use

Shark fin is a delicacy in many Asian countries. Overexploitation of sharks for shark fin trading has led to a drastic reduction in shark population. To monitor international trade of shark fin products and protect the endangered species from further population decline, we present rapid, user-friendly and sensitive diagnostic loop-mediated isothermal amplification (LAMP) and effective polymerase chain reaction (PCR) assays for all twelve CITES-listed shark species. Species-specific LAMP and PCR primers were designed based on cytochrome oxidase I (COI) and NADH2 regions. Our LAMP and PCR assays have been tested on 291 samples from 93 shark and related species. Target shark species could be differentiated from non-target species within three hours from DNA extraction to LAMP assay. The LAMP assay reported here is a simple and robust solution for on-site detection of CITES-listed shark species with shark fin products.

Morphological characterization and phylogenetic relationships of Indochinese box turtles-The Cuora galbinifrons complex

The members of the Indochinese box turtle complex, namely Cuora galbinifrons, Cuora bourreti, and Cuora picturata, rank the most critically endangered turtle species on earth after more than three decades of over-harvesting for food, traditional Chinese medicine, and pet markets. Despite advances in molecular biology, species boundaries and phylogenetic relationships, the status of the C. galbinifrons complex remains unresolved due to the small number of specimens observed and collected in the field. In this study, we present analyses of morphologic characters as well as mitochondrial and nuclear DNA data to reconstruct the species boundaries and systematic relationships within the C. galbinifrons complex. Based on principal component analysis (PCA) and statistical analysis, we found that phenotypic traits partially overlapped among galbinifrons, bourreti, and picturata, and that galbinifrons and bourreti might be only subspecifically distinct. Moreover, we used the mitochondrial genome, COI, and nuclear gene Rag1 under the maximum likelihood criteria and Bayesian inference criteria to elucidate whether C. galbinifrons could be divided into three separate species or subspecies. We found strong support for a sister relationship between picturata and the other two species, and consequently, we recommend maintaining picturata as a full species, and classifying bourreti and galbinifrons as subspecies of C. galbinifrons. These findings provide evidence for a better understanding of the evolutionary histories of these critically endangered turtles.

Sequencing and characterization of the complete mitochondrial genome of Japanese Swellshark (Cephalloscyllium umbratile)

To further comprehend the genome features of Cephalloscyllium umbratile (Carcharhiniformes), an endangered species, the complete mitochondrial DNA (mtDNA) was firstly sequenced and annotated. The full-length mtDNA of C. umbratile was 16,697 bp and contained ribosomal RNA (rRNA) genes, 13 protein-coding genes (PCGs), 23 transfer RNA (tRNA) genes, and a major non-coding control region. Each PCG was initiated by an authoritative ATN codon, except for COX1 initiated by a GTG codon. Seven of 13 PCGs had a typical TAA termination codon, while others terminated with a single T or TA. Moreover, the relative synonymous codon usage of the 13 PCGs was consistent with that of other published Carcharhiniformes. All tRNA genes had typical clover-leaf secondary structures, except for tRNA-Ser (GCT), which lacked the dihydrouridine 'DHU' arm. Furthermore, the analysis of the average Ka/Ks in the 13 PCGs of three Carcharhiniformes species indicated a strong purifying selection within this group. In addition, phylogenetic analysis revealed that C. umbratile was closely related to Glyphis glyphis and Glyphis garricki. Our data supply a useful resource for further studies on genetic diversity and population structure of C. umbratile.

Draft sequencing and assembly of the genome of the world's largest fish, the whale shark: Rhincodon typus Smith 1828

BACKGROUND

The whale shark (Rhincodon typus) has by far the largest body size of any elasmobranch (shark or ray) species. Therefore, it is also the largest extant species of the paraphyletic assemblage commonly referred to as fishes. As both a phenotypic extreme and a member of the group Chondrichthyes - the sister group to the remaining gnathostomes, which includes all tetrapods and therefore also humans - its genome is of substantial comparative interest. Whale sharks are also listed as an endangered species on the International Union for Conservation of Nature's Red List of threatened species and are of growing popularity as both a target of ecotourism and as a charismatic conservation ambassador for the pelagic ecosystem. A genome map for this species would aid in defining effective conservation units and understanding global population structure.

RESULTS

We characterised the nuclear genome of the whale shark using next generation sequencing (454, Illumina) and de novo assembly and annotation methods, based on material collected from the Georgia Aquarium. The data set consisted of 878,654,233 reads, which yielded a draft assembly of 1,213,200 contigs and 997,976 scaffolds. The estimated genome size was 3.44Gb. As expected, the proteome of the whale shark was most closely related to the only other complete genome of a cartilaginous fish, the holocephalan elephant shark. The whale shark contained a novel Toll-like-receptor (TLR) protein with sequence similarity to both the TLR4 and TLR13 proteins of mammals and TLR21 of teleosts. The data are publicly available on GenBank, FigShare, and from the NCBI Short Read Archive under accession number SRP044374.

CONCLUSIONS

This represents the first shotgun elasmobranch genome and will aid studies of molecular systematics, biogeography, genetic differentiation, and conservation genetics in this and other shark species, as well as providing comparative data for studies of evolutionary biology and immunology across the jawed vertebrate lineages.

Complete mitochondrial DNA genome of bonnethead shark, Sphyrna tiburo, and phylogenetic relationships among main superorders of modern elasmobranchs

Elasmobranchs are one of the most diverse groups in the marine realm represented by 18 orders, 55 families and about 1200 species reported, but also one of the most vulnerable to exploitation and to climate change. Phylogenetic relationships among main orders have been controversial since the emergence of the Hypnosqualean hypothesis by Shirai (1992) that considered batoids as a sister group of sharks. The use of the complete mitochondrial DNA (mtDNA) may shed light to further validate this hypothesis by increasing the number of informative characters. We report the mtDNA genome of the bonnethead shark Sphyrna tiburo, and compare it with mitogenomes of other 48 species to assess phylogenetic relationships. The mtDNA genome of S. tiburo, is quite similar in size to that of congeneric species but also similar to the reported mtDNA genome of other Carcharhinidae species. Like most vertebrate mitochondrial genomes, it contained 13 protein coding genes, two rRNA genes and 22 tRNA genes and the control region of 1086 bp (D-loop). The Bayesian analysis of the 49 mitogenomes supported the view that sharks and batoids are separate groups.

PCR-free shotgun sequencing of the stone loach mitochondrial genome (Barbatula barbatula)

The complete mitochondrial genome of the stone loach Barbatula barbatula (Linnaeus, 1758) (Actinopterygii: Cypriniformes: Nemacheilidae) has been sequenced using a genome-skimming approach on an Illumina Hiseq 2500 platform. The mitochondrial genome of B. barbatula was determined to be 16,630 bp long and presents an organization typical of vertebrate mitogenomes. The mean coverage was 82× with a minimum coverage of 33× for the control region and 52× for the remaining part of the genome. A phylogenetic analysis of the Nemacheilidae family shows the monophyly of the Barbatula genus with strong support.