A dated molecular phylogeny of manta and devil rays (Mobulidae) based on mitogenome and nuclear sequences

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.

The complete mitochondrial genome of the Reef Manta Ray, Mobula alfredi, from Hawaii

We provide the complete mitochondrial genome of the reef manta ray, Mobula alfredi, using an ezRAD approach. The total length of the mitogenome was 18,166 bp and contained 13 protein-coding genes, 22 transfer RNAs genes, two ribosomal RNA genes, and one non-coding control region. The gene organization and length are similar to other Mobula species. This reference mitogenome that includes the control region is expected to be a valuable resource for molecular-based species identification, population genomics, and phylogeography.

Tracing Patterns and Biodiversity Aspects of the Overlooked Skates and Rays (Subclass Elasmobranchii, Superorder Batoidea) in Greece

Species belonging to the superorder Batoidea have been poorly assessed due to the lack of information on their life history aspects and their limited economic value. This work presents an overview of skates and rays inhabiting the marine Greek waters and reports biodiversity issues requiring resolution for conservation purposes. Overall, 30 species from nine families and 16 genera have been documented within the past 22 years, based on the available literature and technical reports from research surveys of the Hellenic Centre for Marine Research (HCMR). However, 28 species are currently confirmed, since the presence of two rajids has not been sufficiently demonstrated and hence is considered as doubtful. Recent changes in nomenclature allowed us to replace old with new names in four species on the Greek list; patterns in the frequency of occurrence were observed and species were assigned into five categories; and diversity and misidentification issues were reported by family. Although Greek waters are oligotrophic and not considered a biodiversity hot spot for elasmobranchs, a high number of batoid species is documented in the area; therefore, the need to reinforce knowledge on biological aspects of skates and rays, define their status and identify their main threats is essential.

Skimming for barcodes: rapid production of mitochondrial genome and nuclear ribosomal repeat reference markers through shallow shotgun sequencing

DNA barcoding is critical to conservation and biodiversity research, yet public reference databases are incomplete. Existing barcode databases are biased toward cytochrome oxidase subunit I (COI) and frequently lack associated voucher specimens or geospatial metadata, which can hinder reliable species assignments. The emergence of metabarcoding approaches such as environmental DNA (eDNA) has necessitated multiple marker techniques combined with barcode reference databases backed by voucher specimens. Reference barcodes have traditionally been generated by Sanger sequencing, however sequencing multiple markers is costly for large numbers of specimens, requires multiple separate PCR reactions, and limits resulting sequences to targeted regions. High-throughput sequencing techniques such as genome skimming enable assembly of complete mitogenomes, which contain the most commonly used barcoding loci (e.g., COI, 12S, 16S), as well as nuclear ribosomal repeat regions (e.g., ITS1&2, 18S). We evaluated the feasibility of genome skimming to generate barcode references databases for marine fishes by assembling complete mitogenomes and nuclear ribosomal repeats. We tested genome skimming across a taxonomically diverse selection of 12 marine fish species from the collections of the National Museum of Natural History, Smithsonian Institution. We generated two sequencing libraries per species to test the impact of shearing method (enzymatic or mechanical), extraction method (kit-based or automated), and input DNA concentration. We produced complete mitogenomes for all non-chondrichthyans (11/12 species) and assembled nuclear ribosomal repeats (18S-ITS1-5.8S-ITS2-28S) for all taxa. The quality and completeness of mitogenome assemblies was not impacted by shearing method, extraction method or input DNA concentration. Our results reaffirm that genome skimming is an efficient and (at scale) cost-effective method to generate all mitochondrial and common nuclear DNA barcoding loci for multiple species simultaneously, which has great potential to scale for future projects and facilitate completing barcode reference databases for marine fishes.

The distribution of manta rays in the western North Atlantic Ocean off the eastern United States

In 2018, the giant manta ray was listed as threatened under the U.S. Endangered Species Act. We integrated decades of sightings and survey effort data from multiple sources in a comprehensive species distribution modeling (SDM) framework to evaluate the distribution of giant manta rays off the eastern United States, including the Gulf of Mexico. Manta rays were most commonly detected at productive nearshore and shelf-edge upwelling zones at surface thermal frontal boundaries within a temperature range of approximately 20–30 °C. SDMs predicted highest nearshore occurrence off northeastern Florida during April, with the distribution extending northward along the shelf-edge as temperatures warm, leading to higher occurrences north of Cape Hatteras, North Carolina from June to October, and then south of Savannah, Georgia from November to March as temperatures cool. In the Gulf of Mexico, the highest nearshore occurrence was predicted around the Mississippi River delta from April to June and again from October to November. SDM predictions will allow resource managers to more effectively protect manta rays from fisheries bycatch, boat strikes, oil and gas activities, contaminants and pollutants, and other threats.

Checklist of marine and estuarine fishes from the AlaskaYukon Border, Beaufort Sea, to Cabo San Lucas, Mexico

This paper is a checklist of the fishes that have been documented, through both published and unpublished sources, in marine and estuarine waters, and out 200 miles, from the United States-Canadian border on the Beaufort Sea to Cabo San Lucas, Mexico. A minimum of 241 families and 1,644 species are known within this range, including both native and nonnative species. For each of these species, we include maximum size, geographic and depth ranges, whether it is native or nonnative, as well as a brief mention of any taxonomic issues.

Phylogenomics and species delimitation for effective conservation of manta and devil rays

Practical biodiversity conservation relies on delineation of biologically meaningful units. Manta and devil rays (Mobulidae) are threatened worldwide, yet morphological similarities and a succession of recent taxonomic changes impede the development of an effective conservation strategy. Here, we generate genome-wide single nucleotide polymorphism (SNP) data from a geographically and taxonomically representative set of manta and devil ray samples to reconstruct phylogenetic relationships and evaluate species boundaries under the general lineage concept. We show that nominal species units supported by alternative data sources constitute independently evolving lineages, and find robust evidence for a putative new species of manta ray in the Gulf of Mexico. Additionally, we uncover substantial incomplete lineage sorting indicating that rapid speciation together with standing variation in ancestral populations has driven phylogenetic uncertainty within Mobulidae. Finally, we detect cryptic diversity in geographically distinct populations, demonstrating that management below the species level may be warranted in certain species. Overall, our study provides a framework for molecular genetic species delimitation that is relevant to wide-ranging taxa of conservation concern, and highlights the potential for genomic data to support effective management, conservation and law enforcement strategies.

Kinematics of swimming of the manta ray: three-dimensional analysis of open-water maneuverability

For aquatic animals, turning maneuvers represent a locomotor activity that may not be confined to a single coordinate plane, making analysis difficult, particularly in the field. To measure turning performance in a three-dimensional space for the manta ray (Mobula birostris), a large open-water swimmer, scaled stereo video recordings were collected. Movements of the cephalic lobes, eye and tail base were tracked to obtain three-dimensional coordinates. A mathematical analysis was performed on the coordinate data to calculate the turning rate and curvature (1/turning radius) as a function of time by numerically estimating the derivative of manta trajectories through three-dimensional space. Principal component analysis was used to project the three-dimensional trajectory onto the two-dimensional turn. Smoothing splines were applied to these turns. These are flexible models that minimize a cost function with a parameter controlling the balance between data fidelity and regularity of the derivative. Data for 30 sequences of rays performing slow, steady turns showed the highest 20% of values for the turning rate and smallest 20% of turn radii were 42.65±16.66 deg s^-1 and 2.05±1.26 m, respectively. Such turning maneuvers fall within the range of performance exhibited by swimmers with rigid bodies.

Complete mitochondrial genome of the Devil Ray, Mobula thurstoni (Lloyd, 1908) (Myliobatiformes: Myliobatidae)

The Devil Ray (Mobula thurstoni) is a species with global distribution and is an important species in conservation terms, here we present its complete mitochondrial genome assembled with Illumina sequencing data. The circular genome was 17,610 bp in length, and consists of 13 protein-coding, two ribosomal RNAs (rRNAs), and 22 transfer RNA (tRNA) genes. Base composition is 30.7% A, 29.1% T, 26.5% C, and 13.7% G, and 40.2% GC content. Protein-coding genes present two start codon (ATG and GUG) and seven stop codon (UAA, AUA, UUU, UUA, AAU, CCU, and UAG). The control region possesses the highest A + T (66.6%) content among all mitochondrial regions. These data would contribute to the evolutionary studies of this genus, where there has been recent reclassification.

DNA analysis of traded shark fins and mobulid gill plates reveals a high proportion of species of conservation concern

Continuously increasing demand for plant and animal products causes unsustainable depletion of biological resources. It is estimated that one-quarter of sharks and rays are threatened worldwide and although the global fin trade is widely recognized as a major driver, demand for meat, liver oil, and gill plates also represents a significant threat. This study used DNA barcoding and 16 S rRNA sequencing as a method to identify shark and ray species from dried fins and gill plates, obtained in Canada, China, and Sri Lanka. 129 fins and gill plates were analysed and searches on BOLD produced matches to 20 species of sharks and five species of rays or – in two cases – to a species pair. Twelve of the species found are listed or have been approved for listing in 2017 in the appendices of the Convention on International Trade in Endangered Species of Fauna and Flora (CITES), including the whale shark (Rhincodon typus), which was surprisingly found among both shark fin and gill plate samples. More than half of identified species fall under the IUCN Red List categories ‘Endangered’ and ‘Vulnerable’, raising further concerns about the impacts of this trade on the sustainability of these low productivity species.

Sympathy for the devil: a conservation strategy for devil and manta rays

BACKGROUND

International trade for luxury products, medicines, and tonics poses a threat to both terrestrial and marine wildlife. The demand for and consumption of gill plates (known as Peng Yu Sai, "Fish Gill of Mobulid Ray") from devil and manta rays (subfamily Mobulinae, collectively referred to as mobulids) poses a significant threat to these marine fishes because of their extremely low productivity. The demand for these gill plates has driven an international trade supplied by largely unmonitored and unregulated catches from target and incidental fisheries around the world. Scientific research, conservation campaigns, and legal protections for devil rays have lagged behind those for manta rays despite similar threats across all mobulids.

METHODS

To investigate the difference in attention given to devil rays and manta rays, we examined trends in the scientific literature and updated species distribution maps for all mobulids. Using available information on target and incidental fisheries, and gathering information on fishing and trade regulations (at international, national, and territorial levels), we examined how threats and protective measures overlap with species distribution. We then used a species conservation planning approach to develop the Global Devil and Manta Ray Conservation Strategy, specifying a vision, goals, objectives, and actions to advance the knowledge and protection of both devil and manta rays.

RESULTS AND DISCUSSION

Our literature review revealed that there had been nearly 2.5-times more "manta"-titled publications, than "mobula" or "devil ray"-titled publications over the past 4.5 years (January 2012-June 2016). The majority of these recent publications were reports on occurrence of mobulid species. These publications contributed to updated Area of Occupancy and Extent of Occurrence maps which showed expanded distributions for most mobulid species and overlap between the two genera. While several international protections have recently expanded to include all mobulids, there remains a greater number of national, state, and territory-level protections for manta rays compared to devil rays. We hypothesize that there are fewer scientific publications and regulatory protections for devil rays due primarily to perceptions of charisma that favour manta rays. We suggest that the well-established species conservation framework used here offers an objective solution to close this gap. To advance the goals of the conservation strategy we highlight opportunities for parity in protection and suggest solutions to help reduce target and bycatch fisheries.

Comparative analysis of the mitochondrial genomes of Colletotrichum gloeosporioides sensu lato: insights into the evolution of a fungal species complex interacting with diverse plants

BACKGROUND

The fungal species complex Colletotrichum gloeosporioides sensu lato contains over 20 plant-interacting species. These species exhibit different life styles (e.g., endophytes, foliar and fruit pathogens) and show considerable variation in host and tissue adaptation strategies. Accurate species delimitation in C. gloeosporioides s.l. is very challenging due to nascent lineage boundaries and phenotypic plasticity, which strongly impedes studies of the complex's host-interaction biology. In this study, we first sequenced and compared nine mitogenomes belonging to four C. gloeosporioides s.l. species lineages (C. gloeosporioides, C. fructicola, C. aenigma, and C. siamense s.l.), and evaluated the usefulness of mitogenome sequence in complementing prevailing nuclear markers for species delimitation.

RESULTS

The C. gloeosporioides s.l. mitogenomes ranged between 52,671 and 58,666 bp in size, and each contained an identical set of genes transcribed in the same direction. Compared with previously reported Colletotrichum mitogenomes, these mitogenomes were uniquely featured by: (1) significantly larger genome size due to richer intron content and longer intergenic region; (2) striking GC content elevation at the intergenic region; and (3) considerable intron content variation among different species lineages. Compared with nuclear DNA markers commonly used in phylogeny, the mitogenome nucleotide diversity was extremely low, yet the mitogenome alignment contained the highest number of parsimony informative sites, which allowed the generation of a high-resolution phylogeny recognizing all taxonomic lineages, including ones belonging to the very nascent C. siamense s.l. complex. The tree topology was highly congruent with the phylogeny based on nuclear marker concatenation except for lineages within C. siamense s.l. Further comparative phylogenetic analysis indicated that lineage-specific rapid divergence of GS and SOD2 markers confounded concatenation-based species relationship inference.

CONCLUSIONS

This study sheds light on the evolution of C. gloeosporioides s.l. mitogenomes and demonstrates that mitogenome sequence can complement prevailing nuclear markers in improving species delimitation accuracy. The mitogenome sequences reported will be valuable resources for further genetic studies with C. gloeosporioides s.l. and other Colletotrichum species.

Improving the Conservation of Mediterranean Chondrichthyans: The ELASMOMED DNA Barcode Reference Library

Cartilaginous fish are particularly vulnerable to anthropogenic stressors and environmental change because of their K-selected reproductive strategy. Accurate data from scientific surveys and landings are essential to assess conservation status and to develop robust protection and management plans. Currently available data are often incomplete or incorrect as a result of inaccurate species identifications, due to a high level of morphological stasis, especially among closely related taxa. Moreover, several diagnostic characters clearly visible in adult specimens are less evident in juveniles. Here we present results generated by the ELASMOMED Consortium, a regional network aiming to sample and DNA-barcode the Mediterranean Chondrichthyans with the ultimate goal to provide a comprehensive DNA barcode reference library. This library will support and improve the molecular taxonomy of this group and the effectiveness of management and conservation measures. We successfully barcoded 882 individuals belonging to 42 species (17 sharks, 24 batoids and one chimaera), including four endemic and several threatened ones. Morphological misidentifications were found across most orders, further confirming the need for a comprehensive DNA barcoding library as a valuable tool for the reliable identification of specimens in support of taxonomist who are reviewing current identification keys. Despite low intraspecific variation among their barcode sequences and reduced samples size, five species showed preliminary evidence of phylogeographic structure. Overall, the ELASMOMED initiative further emphasizes the key role accurate DNA barcoding libraries play in establishing reliable diagnostic species specific features in otherwise taxonomically problematic groups for biodiversity management and conservation actions.

DNA barcoding of Mobulid Ray Gill Rakers for Implementing CITES on Elasmobranch in China

The Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES) has been counted on for conserving threatened marine fish since it regulates the commercial international trade of these species. Implementation of the international treaty for Mantas included on CITES Appendix II is challenging due to insufficient information on species identification and markets management. To fill the gap in such aspects, we identified five species of Mobulid rays (Mobula spps. and Manta spp) by using COI and NADH2 mtDNA markers in dried ray gill rakers from Chinese markets, namely, Mobula japonica (representing 54.8% of the sample set), M. tarapacana (14.4%), M. kuhlii (13.3%), M. thurstoni (6.4%), along with Manta birostris (11.2%; CITES Appendix II). The utilization and conservation statuses of these species were discussed. Based on combination of DNA barcodes and key morphological characters, we developed a three-step process for identifying the gill rakers of Mobulid rays which has been adopted by frontline enforcement in China. We hope that our work can serve as a foundation and basis to reinforce objectives of international treaties, regulation of consumer-driven markets, regional cooperation, and national fishery management on endangered elasmobranchs in China as well as related countries.

A potential third Manta Ray species near the Yucatán Peninsula? Evidence for a recently diverged and novel genetic Manta group from the Gulf of Mexico

We present genetic and morphometric support for a third, distinct, and recently diverged group of Manta ray that appears resident to the Yucatán coastal waters of the Gulf of Mexico. Individuals of the genus Manta from Isla Holbox are markedly different from the other described manta rays in their morphology, habitat preference, and genetic makeup. Herein referred to as the Yucatán Manta Ray, these individuals form two genetically distinct groups: (1) a group of mtDNA haplotypes divergent (0.78%) from the currently recognized Manta birostris and M. alfredi species, and (2) a group possessing mtDNA haplotypes of M. birostris and highly similar haplotypes. The latter suggests the potential for either introgressive hybridization between Yucatán Manta Rays and M. birostris, or the retention of ancestral M. birostris signatures among Yucatán Manta Rays. Divergence of the genetically distinct Yucatán Manta Ray from M. birostris appears quite recent (<100,000 YBP) following fit to an Isolation-with-Migration model, with additional support for asymmetrical gene flow from M. birostris into the Yucatán Manta Ray. Formal naming of the Yucatán Manta Ray cannot yet be assigned until an in-depth taxonomic study and further confirmation of the genetic identity of existing type specimens has been performed.

Growth, productivity, and relative extinction risk of a data-sparse devil ray

Devil rays (Mobula spp.) face intensifying fishing pressure to meet the ongoing international demand for gill plates. The paucity of information on growth, mortality, and fishing effort for devil rays make quantifying population growth rates and extinction risk challenging. Furthermore, unlike manta rays (Manta spp.), devil rays have not been listed on CITES. Here, we use a published size-at-age dataset for the Spinetail Devil Ray (Mobula japanica), to estimate somatic growth rates, age at maturity, maximum age, and natural and fishing mortality. We then estimate a plausible distribution of the maximum intrinsic population growth rate (r_max) and compare it to 95 other chondrichthyans. We find evidence that larger devil ray species have low somatic growth rate, low annual reproductive output, and low maximum population growth rates, suggesting they have low productivity. Fishing rates of a small-scale artisanal Mexican fishery were comparable to our estimate of r_max, and therefore probably unsustainable. Devil ray r_max is very similar to that of manta rays, indicating devil rays can potentially be driven to local extinction at low levels of fishing mortality and that a similar degree of protection for both groups is warranted.

Deep-water feeding and behavioral plasticity in Manta birostris revealed by archival tags and submersible observations

Foraging drives many fundamental aspects of ecology, and an understanding of foraging behavior aids in the conservation of threatened species by identifying critical habitats and spatial patterns relevant to management. The world's largest ray, the oceanic manta (Manta birostris) is poorly studied and threatened globally by targeted fisheries and incidental capture. Very little information is available on the natural history, ecology and behavior of the species, complicating management efforts. This study provides the first data on the diving behavior of the species based on data returned from six tagged individuals, and an opportunistic observation from a submersible of a manta foraging at depth. Pop-off archival satellite tags deployed on mantas at the Revillagigedo Archipelago, Mexico recorded seasonal shifts in diving behavior, likely related to changes in the location and availability of zooplankton prey. Across seasons, mantas spent a large proportion of their time centered around the upper limit of the thermocline, where zooplankton often aggregate. Tag data reveal a gradual activity shift from surface waters to 100-150m across the tagging period, possibly indicating a change in foraging behavior from targeting surface-associated zooplankton to vertical migrators. The depth ranges accessed by mantas in this study carry variable bycatch risks from different fishing gear types. Consequently, region-specific data on diving behavior can help inform local management strategies that reduce or mitigate bycatch of this vulnerable species.

Comparative anatomy of the extraocular muscles in four Myliobatoidei rays (Batoidea, Myliobatiformes)

Extraocular muscles are classically grouped as four rectus and two oblique muscles. However, their description and potential associations with species behavior are limited. The objective was to characterize extraocular muscles in four Myliobatoidei rays from diverse habitats with divergent behaviors. Heads (10 per species) of Dasyatis hypostigma, Gymnura altavela, Mobula thurstoni and Pteroplatytrygon violacea were decalcified and dissected to characterize and describe extraocular muscles. Principal component analysis (PCA) was used to evaluate relationships between muscle length and species; for P. violacea, D. hypostigma and G. altavela, these were qualitatively and quantitatively consistent with the general pattern of extraocular muscles in vertebrates. In contrast, for M. thurstoni, the two oblique muscles were completely fused and there was a seventh extraocular muscle, named m. lateral rectus β (both were apparently novel findings in this species). There were also significant differences in eye disposition in the chondrocranium. The PCA axis 1 (rectus muscles) and PCA axis 2 (oblique muscles) accounted for 98.47% of data variability. Extraocular muscles had significant differences in length and important anatomical differences among sampled species that facilitated grouping species according to their life history. In conclusion, extraocular muscles are not uniform in all vertebrate species, thereby providing another basis for comparative studies.

The phylogenetic position of the giant devil ray Mobula mobular (Bonnaterre, 1788) (Myliobatiformes, Myliobatidae) inferred from the mitochondrial genome

The giant devil ray, Mobula mobular, is a member of one of the most distinct groups of cartilaginous fishes, the Mobulidae (manta and devil rays), and is the only mobulid assessed as Endangered due its restricted distribution, high bycatch mortality and suspected population decline. The complete mitochondrial genome is 18 913 base pairs in length and comprises 2 rRNAs, 13 protein-coding genes, 22 tRNAs and 2 non-coding regions. Comparison with the partial mitogenome of M. japanica suggests a sister-cryptic species complex and two different taxonomic units. However, the limited divergence within the species (>99.9% genetic identity) may be the result of a geographically and numerically restricted population of M. mobular within the Mediterranean Sea.