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

A comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important 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.

Morphology and phylogenetic significance of the pelvic articular region in elasmobranchs (Chondrichthyes)

The morphology of paired fins is commonly overlooked in morphological studies, particularly the pelvic girdle and fins. Consequently, previous phylogenetic studies incorporating morphological data used few skeletal characters from this complex. In this paper, the phylogenetic significance of pelvic articular characters for elasmobranchs is discussed in light of the morphological variation observed in 130 species, the most comprehensive study exploring the morphology of the pelvic girdle done so far. The 10 morphological characters proposed herein for the pelvic articulation were incorporated into a molecular matrix of NADH2 sequences and submitted to an analysis of maximum parsimony employing extended implied weighting. The most stable tree was selected based on the distortion coefficients, SPR distances (subtree pruning and regrafting) and fit values. Some of the striking synapomorphies recovered within elasmobranchs include the presence of an articular surface for the first enlarged pelvic radial supporting Elasmobranchii and the pelvic articular region for the basipterygium extending from the posterolatral margin of the pelvic girdle over its lateral surface in Echinorhinus + Hexanchiformes. Additionally, the proposed characters and their distributions are discussed considering the relationships recovered and also compared with previous morphological and molecular phylogenetic hypotheses.

Characterization of the complete mitochondrial genomes of two Critically Endangered wedgefishes: Rhynchobatus djiddensis and Rhynchobatus australiae

We present the complete mitochondrial genomes of the Critically Endangered whitespotted wedgefish, Rhynchobatus djiddensis (Forsskål, 1775), and bottlenose wedgefish, Rhynchobatus australiae (Whitley, 1939), with the R. djiddensis mitogenome documented for the first time. The genomes for R. djiddensis and R. australiae are 16,799 and 16,805 bp in length, respectively. Both comprise 13 protein-coding regions, 22 tRNA genes, two rRNA genes, and a non-coding control region. All protein-coding regions consistently start with the ATG start codon; however, the alternative start codon GTG is observed at the start of the COX1 gene. NADH2, COX2, and NADH4 have incomplete stop codons: T or TA, and tRNA^Leu and tRNA^Ser , have atypical codons: UAA, UGA, GCU, and UAG. The phylogenetic analysis places R. djiddensis and R. australiae within the Rhynchobatus genus, separate from other families in the order Rhinopristiformes. We also highlight the most variable gene regions to expedite future primer design, of which NADH2 was the most variable (4.5%) when taking gene length into account. These molecular resources could promote the taxonomic resolution of the whitespotted wedgefish species complex and aid in the genetic characterization of populations of these and related species.

Molecular identification and first demographic insights of sharks based on artisanal fisheries bycatch in the Pacific Coast of Colombia: implications for conservation

The Pacific coast of Colombia is characterized by mangrove ecosystems which play a crucial role as possible nurseries for juvenile sharks. However, trophic food webs from coastal ecosystems are heavily disturbed by increased fishing pressure, which affects numerous shark species. In this region of the Eastern Tropical Pacific (ETP), fisheries' data from coastal areas are scarce and unspecific, as most sharks from artisanal fisheries are landed decapitated and finless, making their morphological identification difficult. For the establishment and implementation of effective regional conservation and management policies, information on the diversity and population dynamics of shark species is crucial. We therefore sequenced the mitochondrial NADH2 gene of 696 samples taken from fishermen's landings of shark's bycatch along the Colombian north Pacific coast. We were able to identify 14 species of sharks, two of the most abundant species were Sphyrna lewini and Carcharhinus falciformis, both evaluated on IUCN the Red List of Threatened species (Critically Endangered and Vulnerable) and CITES regulated. We found low genetic diversity in the sampled area increasing the concern for both species in the region, even more considering that the majority of individuals were juveniles. Our results showed the importance of genetic markers for first population genetic insights as a complementary tool during the decision-making process in management plans. For this specific region, strategies such as the delimitation of conservation priority areas or the regulation of fishing gears could help improve the sustainability of shark populations in the Colombian Pacific.

Novel duplication remnant in the first complete mitogenome of Hemitriakis japanica and the unique phylogenetic position of family Triakidae

In this study, we firstly determined the complete mitogenome of the Japanese topeshark (Hemitriakis japonica), which belong to the family Triakidae and was assessed as Endangered A2d on the IUCN Red List in 2021. The mitogenome is 17,301 bp long, has a high AT content (60.0%), and contains 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, a control region and specially a 594 bp-long non-coding region between Cytb gene and tRNA-Thr gene. The novel non-coding region share high sequence similarity with segments of the former and latter genes, so it was recognized as a duplication remnant. In addition, the Cytb gene and tRNA-Thr gene tandemly duplicated twice while accompanied by being deleted once at least. This is the first report of mitogenomic gene-arrangement in Triakidae. The phylogenetic trees were constructed using Bayesian inference (BI) and maximum likelihood (ML) methods based on the mitogenomic data of 51 shark species and two outgroups. In summary, basing on a novel type of gene rearrangements in houndshark mitogenome, the possibly rearranged process was analyzed and contributed further insight of shark mitogenomes evolution and phylogeny.

Complete Mitochondrial DNA Genome of Nine Species of Sharks and Rays and Their Phylogenetic Placement among Modern Elasmobranchs

Chondrichthyes occupy a key position in the phylogeny of vertebrates. The complete sequence of the mitochondrial genome (mitogenome) of four species of sharks and five species of rays was obtained by whole genome sequencing (DNA-seq) in the Illumina HiSeq2500 platform. The arrangement and features of the genes in the assembled mitogenomes were identical to those found in vertebrates. Both Maximum Likelihood (ML) and Bayesian Inference (BI) analyses were used to reconstruct the phylogenetic relationships among 172 species (including 163 mitogenomes retrieved from GenBank) based on the concatenated dataset of 13 individual protein coding genes. Both ML and BI analyses did not support the “Hypnosqualea” hypothesis and confirmed the monophyly of sharks and rays. The broad notion in shark phylogeny, namely the division of sharks into Galeomorphii and Squalomorphii and the monophyly of the eight shark orders, was also supported. The phylogenetic placement of all nine species sequenced in this study produced high statistical support values. The present study expands our knowledge on the systematics, genetic differentiation, and conservation genetics of the species studied, and contributes to our understanding of the evolutionary history of Chondrichthyes.

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.

Multiple Evolutionary Origins of Ubiquitous Cu2+ and Zn2+ Binding in the S100 Protein Family

The S100 proteins are a large family of signaling proteins that play critical roles in biology and disease. Many S100 proteins bind Zn2+, Cu2+, and/or Mn2+ as part of their biological functions; however, the evolutionary origins of binding remain obscure. One key question is whether divalent transition metal binding is ancestral, or instead arose independently on multiple lineages. To tackle this question, we combined phylogenetics with biophysical characterization of modern S100 proteins. We demonstrate an earlier origin for established S100 subfamilies than previously believed, and reveal that transition metal binding is widely distributed across the tree. Using isothermal titration calorimetry, we found that Cu2+ and Zn2+ binding are common features of the family: the full breadth of human S100 paralogs-as well as two early-branching S100 proteins found in the tunicate Oikopleura dioica-bind these metals with μM affinity and stoichiometries ranging from 1:1 to 3:1 (metal:protein). While binding is consistent across the tree, structural responses to binding are quite variable. Further, mutational analysis and structural modeling revealed that transition metal binding occurs at different sites in different S100 proteins. This is consistent with multiple origins of transition metal binding over the evolution of this protein family. Our work reveals an evolutionary pattern in which the overall phenotype of binding is a constant feature of S100 proteins, even while the site and mechanism of binding is evolutionarily labile.

Complete mitochondrial genome of the porbeagle shark, Lamna nasus (Chondrichthyes, Lamnidae)

We report for the first time, the complete mitochondrial genome sequence of the porbeagle shark, Lamna nasus, from a specimen collected from offshore waters of New England, USA in the western North Atlantic Ocean. The genome structure of this species is similar to the other reported shark mitogenomes. The genome sequence has a total length of 16,697 bases; similar in size to the mtDNA genomes reported for other lamnid species. A Bayesian phylogenetic tree was reconstructed for the Lamnidae family using mitogenome sequences available in the Genbank.