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

New mitochondrial gene order arrangements and evolutionary implications in the class Octocorallia

The complete mitochondrial genomes of octocorals typically range from 18.5 kb to 20.5 kb in length and include 14 protein-coding genes (PCGs), two ribosomal RNA genes and one tRNA. To date, seven different gene orders (A-G) have been described, yet comprehensive investigations of the actual number of arrangements, as well as comparative analyses and evolutionary reconstructions of mitochondrial genome evolution within the whole class Octocorallia, have been often overlooked. Here, we considered the complete mitochondrial genomes available for octocorals and explored their structure and gene order variability. Our results updated the actual number of mitochondrial gene order arrangements so far known for octocorals from 7 to 14 and allowed us to explore and preliminarily discuss the role of some of the structural and functional factors in the mitogenomes. We performed comparative mitogenomic analyses on the existing and novel octocoral gene orders, considering different mitogenomic structural features such as genome size, GC percentage, AT and GC skewness. The mitochondrial gene order history mapped on a recently published nuclear loci phylogeny showed that the most common rearrangement events in octocorals are inversions, inverted transpositions and transpositions. Furthermore, gene order rearrangement events were restricted only to some regions of the tree. Overall, different rearrangement events arose independently and from the ancestral and most common gene order, instead of being derived from other rearranged orders. Finally, our data demonstrate how the study of mitochondrial gene orders can be used to explore the evolution of octocorals and in some cases can be used to assess the phylogenetic placement of certain taxa.

Endemic Radiation of African Moonfish, Selene dorsalis (Gill 1863), in the Eastern Atlantic: Mitogenomic Characterization and Phylogenetic Implications of Carangids (Teleostei: Carangiformes)

This study offers an in-depth analysis of the mitochondrial genome of Selene dorsalis (Gill 1863), a species native to the Eastern Atlantic Ocean. The circular mitochondrial DNA molecule measures 16,541 base pairs and comprises 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA genes, and a control region (CR). The nucleotide composition exhibits a notable adenine-thymine (AT) bias, accounting for 53.13%, which aligns with other species in the Carangidae family. Most PCGs initiate with the ATG codon, with the exception of Cytochrome C oxidase subunit I, which starts with GTG. Analysis of relative synonymous codon usage reveals that leucine and serine are the most prevalent amino acids in the mitochondrial genome of S. dorsalis and its congeners (S. vomer and S. setapinnis). All tRNAs display the typical cloverleaf structure, though tRNA Serine (S1) lacks a dihydrouracil arm. Pairwise comparisons of synonymous and nonsynonymous substitutions for all PCGs yielded values below '1', indicating strong purifying selection. The CR spans 847 bp, representing 5.12% of the mitochondrial genome, and is characterized by high AT content (62.81%). It is situated between tRNA-Pro (TGG) and tRNA-Phe (GAA). The CR contains conserved sequence blocks, with CSB-1 being the longest at 22 bp and CSB-D the shortest at 18 bp. Phylogenetic analysis, using Bayesian and Maximum-likelihood trees constructed from concatenated PCGs across 72 species, successfully differentiates S. dorsalis from other carangids. This study also explores how ocean currents and gyres might influence lineage diversification and parapatric speciation of Selene species between the Atlantic and Pacific Oceans. These results highlight the importance of the mitochondrial genome in elucidating the structural organization and evolutionary dynamics of S. dorsalis and its relatives within marine ecosystems.

The first complete mitochondrional genome of Anopheles gibbinsi using a skimming sequencing approach

Mosquitoes belonging to the genus Anopheles are the only vectors of human malaria. Anopheles gibbinsi has been linked to malaria transmission in Kenya, with recent collections in Zambia reporting the mosquito species exhibiting zoophilic and exophilic behavioral patterns with occasional contact with humans. Given the paucity of genetic data, and challenges to identification and molecular taxonomy of the mosquitoes belonging to the Anopheles genus; we report the first complete mitochondrial genome of An. gibbinsi using a genome skimming approach. An Illumina Novaseq 6000 platform was used for sequencing, the length of the mitochondrial genome was 15401 bp, with 78.5% AT content comprised of 37 genes. Phylogenetic analysis by maximum likelihood using concatenation of the 13 protein coding genes demonstrated that An. marshallii was the closest relative based on existing sequence data. This study demonstrates that the skimming approach is an inexpensive and efficient approach for mosquito species identification and concurrent taxonomic rectification, which may be a useful alternative for generating reference sequence data for evolutionary studies among the Culicidae.

A CRISPR-based strategy for targeted sequencing in biodiversity science

Many applications in molecular ecology require the ability to match specific DNA sequences from single- or mixed-species samples with a diagnostic reference library. Widely used methods for DNA barcoding and metabarcoding employ PCR and amplicon sequencing to identify taxa based on target sequences, but the target-specific enrichment capabilities of CRISPR-Cas systems may offer advantages in some applications. We identified 54,837 CRISPR-Cas guide RNAs that may be useful for enriching chloroplast DNA across phylogenetically diverse plant species. We tested a subset of 17 guide RNAs in vitro to enrich plant DNA strands ranging in size from diagnostic DNA barcodes of 1,428 bp to entire chloroplast genomes of 121,284 bp. We used an Oxford Nanopore sequencer to evaluate sequencing success based on both single- and mixed-species samples, which yielded mean chloroplast sequence lengths of 2,530-11,367 bp, depending on the experiment. In comparison to mixed-species experiments, single-species experiments yielded more on-target sequence reads and greater mean pairwise identity between contigs and the plant species' reference genomes. But nevertheless, these mixed-species experiments yielded sufficient data to provide ≥48-fold increase in sequence length and better estimates of relative abundance for a commercially prepared mixture of plant species compared to DNA metabarcoding based on the chloroplast trnL-P6 marker. Prior work developed CRISPR-based enrichment protocols for long-read sequencing and our experiments pioneered its use for plant DNA barcoding and chloroplast assemblies that may have advantages over workflows that require PCR and short-read sequencing. Future work would benefit from continuing to develop in vitro and in silico methods for CRISPR-based analyses of mixed-species samples, especially when the appropriate reference genomes for contig assembly cannot be known a priori.

New species of redbait from the Philippines (Teleostei, Emmelichthyidae, Emmelichthys)

We describe a new species of redbait in the genus Emmelichthys collected from fish markets on Panay and Cebu islands in the Visayas region of the Philippines. The species is externally similar to E.struhsakeri but is diagnosable by two prominent fleshy papillae associated with the cleithrum and fewer pectoral-fin rays (18-19 vs. 19-21) and gill rakers (30-33 vs. 34-41). Additionally, mitochondrial DNA differentiates this taxon from other species of Emmelichthys. We generate mitochondrial genomes for two of the three type specimens and several other emmelichthyids to place the new taxon in a phylogenetic context. Analysis of the protein-coding mitochondrial loci calls into question the monophyly of two emmelichthyid genera (Emmelichthys and Erythrocles) and highlights the need for subsequent analyses targeting the intrarelationships of the Emmelichthyidae.

Insights into the nuclear and mitochondrial genome of the Lemon shark Negaprion brevirostris using low-coverage sequencing: Genome size, repetitive elements, mitochondrial genome, and phylogenetic placement

The Lemon shark Negaprion brevirostris is an important species experiencing conservation issues that is in need of genomic resources. Herein, we conducted a genome survey sequencing in N. brevirostris and determined genome size, explored repetitive elements, assembled and annotated the 45S rRNA DNA operon, and assembled and described in detail the mitochondrial genome. Lastly, the phylogenetic position of N. brevirostris in the family Carcharhinidae was examined using translated protein coding genes. The estimated haploid genome size ranged between 2.29 and 2.58 Gbp using a k-mer analysis, which is slightly below the genome size estimated for other sharks belonging to the family Carcharhinidae. Using a k-mer analysis, approx. 64-71 % of the genome of N. brevirostris was composed of repetitive elements. A relatively large proportion of the 'repeatome' could not be annotated. Taking into account only annotated repetitive elements, Class I - Long Interspersed Nuclear Element (LINE) were the most abundant repetitive elements followed by Class I - Penelope and Satellite DNA. The nuclear ribosomal operon was fully assembled. The AT-rich complete mitochondrial genome was 16,703 bp long and encoded 13 protein coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Negaprion brevirostris is closely related to the genera Carcharhinus, Glyphis and Lamiopsis in the family Carcharinidae. This new genomic resources will aid with the development of conservation plans for this large coastal shark.

Mitogenomic Characterization and Phylogenetic Placement of African Hind, Cephalopholis taeniops: Shedding Light on the Evolution of Groupers (Serranidae: Epinephelinae)

The global exploration of evolutionary trends in groupers, based on mitogenomes, is currently underway. This research extensively investigates the structure of and variations in Cephalopholis species mitogenomes, along with their phylogenetic relationships, focusing specifically on Cephalopholis taeniops from the Eastern Atlantic Ocean. The generated mitogenome spans 16,572 base pairs and exhibits a gene order analogous to that of the ancestral teleost's, featuring 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an AT-rich control region. The mitogenome of C. taeniops displays an AT bias (54.99%), aligning with related species. The majority of PCGs in the mitogenome initiate with the start codon ATG, with the exceptions being COI (GTG) and atp6 (TTG). The relative synonymous codon usage analysis revealed the maximum abundance of leucine, proline, serine, and threonine. The nonsynonymous/synonymous ratios were <1, which indicates a strong negative selection among all PCGs of the Cephalopholis species. In C. taeniops, the prevalent transfer RNAs display conventional cloverleaf secondary structures, except for tRNA-serine (GCT), which lacks a dihydrouracil (DHU) stem. A comparative examination of conserved domains and sequence blocks across various Cephalopholis species indicates noteworthy variations in length and nucleotide diversity. Maximum likelihood, neighbor-joining, and Bayesian phylogenetic analyses, employing the concatenated PCGs and a combination of PCGs + rRNAs, distinctly separate all Cephalopholis species, including C. taeniops. Overall, these findings deepen our understanding of evolutionary relationships among serranid groupers, emphasizing the significance of structural considerations in mitogenomic analyses.

Plumbing the depths with environmental DNA (eDNA): Metabarcoding reveals biodiversity zonation at 45-60 m on mesophotic coral reefs

Mesophotic coral ecosystems (MCEs) are tropical reefs found at depths of ~30-150 m, below the region most heavily impacted by heat stress and other disturbances. Hence, MCEs may serve as potential refugia for threatened shallow reefs, but they also harbour depth-endemic fauna distinct from shallow reefs. Previous studies have characterized biodiversity patterns along depth gradients, but focussed primarily on conspicuous taxa (fishes, corals, etc.). Environmental DNA (eDNA) metabarcoding offers a more holistic approach to assess biodiversity patterns across the tree of life. Here, we use three metabarcoding assays targeting fishes (16S rRNA), eukaryotes (18S rDNA) and metazoans (COI) to assess biodiversity change from the surface to ~90 m depth across 15-m intervals at three sites within the Hawaiian Archipelago. We observed significant community differences between most depth zones, with distinct zonation centred at 45-60 m for eukaryotes and metazoans, but not for fishes. This finding may be attributable to the higher mobility of reef fishes, although methodological limitations are likely a contributing factor. The possibility for MCEs to serve as refugia is not excluded for fishes, but invertebrate communities >45 m are distinct, indicating limited connectivity for the majority of reef fauna. This study provides a new approach for surveying biodiversity on MCEs, revealing patterns in a much broader context than the limited-taxon studies that comprise the bulk of our present knowledge.

Biodiversity of Philippine marine fishes: A DNA barcode reference library based on voucher specimens

Accurate identification of fishes is essential for understanding their biology and to ensure food safety for consumers. DNA barcoding is an important tool because it can verify identifications of both whole and processed fishes that have had key morphological characters removed (e.g., filets, fish meal); however, DNA reference libraries are incomplete, and public repositories for sequence data contain incorrectly identified sequences. During a nine-year sampling program in the Philippines, a global biodiversity hotspot for marine fishes, we developed a verified reference library of cytochrome c oxidase subunit I (COI) sequences for 2,525 specimens representing 984 species. Specimens were primarily purchased from markets, with additional diversity collected using rotenone or fishing gear. Species identifications were verified based on taxonomic, phenotypic, and genotypic data, and sequences are associated with voucher specimens, live-color photographs, and genetic samples catalogued at Smithsonian Institution, National Museum of Natural History. The Biodiversity of Philippine Marine Fishes dataset is released herein to increase knowledge of species diversity and distributions and to facilitate accurate identification of market fishes.

Decoding dissolved information: environmental DNA sequencing at global scale to monitor a changing ocean

The use of environmental DNA (eDNA) technology for environmental monitoring is rapidly expanding, with applications for fisheries, coral reefs, harmful algal blooms, invasive and endangered species, and biodiversity monitoring. By enabling detection of species over space and time, eDNA fulfills a fundamental need of environmental surveys. Traditional surveys are expensive, require significant capital expenditure, and can be destructive; eDNA offers promise for cheaper, less invasive, and higher-resolution (i.e. genetic) assessments of environments and stocks. However, challenges in quantification, detection limits, biobanking capacity, reference databases, and data management and integration remain significant hurdles to efficient eDNA monitoring at global and decadal scale. Here, we consider the current state of eDNA technology and its suitability for the problems for which it is being used. We explore the current best practices, the logistical and social challenges that prevent eDNA from widespread adoption and benefit, and the emerging technologies that may address those challenges.