Rapid mitochondrial genome sequencing based on Oxford Nanopore Sequencing and a proxy for vertebrate species identification

Comparative analysis of bottom trawl and nanopore sequencing in fish biodiversity assessment: The sylt outer reef example

The assessment of fish diversity is crucial for effective conservation and management strategies, especially in ecologically sensitive regions such as marine protected areas. This study contrasts the effectiveness of environmental DNA (eDNA) metabarcoding analysis employing Nanopore technology with compare beam trawl surveys at the Sylt Outer Reef, a Natura 2000 site in the North Sea, Germany. Out of the 17 fish species caught in a bottom trawl (using a 3m beam trawl), 14 were also identified through eDNA extracted from water samples. The three species not detected in the eDNA results were absent because they lacked representation in public DNA databases. The eDNA method detected twice as many fish species as the beam trawl, totalling 36 species, of which 14 were also detected by the trawl. Additionally, the selection of primers (Mifish) facilitated the identification of one marine mammal species, the harbour porpoise. In conclusion, the findings underscore the potential of eDNA coupled with MinION sequencing (Long read technology) as a robust tool for biodiversity assessment, surpassing traditional methods in detecting species richness.

Exploring Barbronia species diversity and phylogenetic relationship within Suborder Erpobdelliformes (Clitellata: Annelida)

Background

Barbronia, a genus of freshwater macrophagous leeches, belongs to Erpobdelliformes (Salifidae: Clitellata: Annelida), and B. weberi, a well-known leech within this genus, has a worldwide distribution. However, the systematics of Barbronia have not yet been adequately investigated, primarily due to a few molecular markers, and only 20 Barbronia sequences available in the GenBank database. This gap significantly limits our understanding of the Barbronia species identification, as well as the phylogenetic placement of the genus Barbronia within Salifidae.

Methods

Next-generation sequencing (NGS) was used to simultaneously capture the entire mitochondrial genome and the full-length 18S/28S rDNA sequences. The species boundary of Barbronia species was estimated using bGMYC and bPTP methods, based on all available Barbronia COI sequences. Uncorrected COI p-distance was calculated in MEGA. A molecular data matrix consisting of four loci (COI, 12S, 18S, and 28S rDNA) for outgroups (three Haemopis leeches) and 49 erpobdellid leeches, representing eight genera within the Suborder Erpobdelliformes was aligned using MAFFT and LocARNA. This matrix was used to reconstruct the phylogenetic relationship of Barbronia via Bayesian inference (BI) and the maximum likelihood (ML) method.

Results

The full lengths of the mitochondrial genome, 18S and 28S rDNAs of B. cf. gwalagwalensis, are 14847 bp, 1876 bp 1876 bp, and 2863 bp, respectively. Both bGMYC and bPTP results based on COI data are generally congruent, suggesting that the previously proposed taxa (B. arcana, B. weberi formosana, and B. wuttkei or Erpobdella wuttkei) are synonyms of B. weberi. The specimens listed in the B. gwalagwalensis group, however, are split into at least two Primary Species Hypotheses (PSHs). The p-distance of the first PSH is less than 1.3% but increased to 4.5% when including the secondary PSH (i.e., B. cf. gwalagwalensis). In comparison, the interspecific p-distance between the B. weberi group and the B. gwalagwalensis group ranged from 6.4% to 8.7%, and the intraspecific p-distance within the B. weberi group is less than 0.8%. Considering the species delimitation results and the sufficient large p-distance, the specimen sampled in China is treated as B. cf. gwalagwalensis. The monophyly of the four Erpobdelliformes families Salifidae, Orobdellidae, Gastrostomobdellidae sensu stricto and Erpobdellidae is well supported in ML and BI analysis based on a data of four markers. Within the Salifidae, a well-supported Barbronia is closely related to a clade containing Odontobdella and Mimobdella, and these three genera are sister to a clade consisted of Salifa and Linta. According to the results of this study, the strategy of simultaneous obtaining both whole mitochondria and nuclear markers from extensively sampled Salifids species using NGS is expected to fathom both the species diversity of B. gwalagwalensis and the evolutionary relationship of Salifidae.

Rapid and collective determination of the complete "hot-spring frog" mitochondrial genome containing long repeat regions using Nanopore sequencing

The mitochondrial genome (mt-genome) is one of the promising molecular markers for phylogenetics and population genetics. Recently, various mt-genomes have been determined rapidly by using massively parallel sequencers. However, the control region (CR, also called D-loop) in mt-genomes remain difficult to precisely determine due to the presence of repeat regions. Here, using Nanopore sequencing, we succeeded in rapid and collective determination of complete mt-genome of the hot-spring frog, Buergeria japonica, and found that its mt-genome size was 22,274 bp including CR (6,929 bp) with two types of tandem repeat motifs forming repeat regions. Comparison of assembly strategies revealed that the long- and short-read data combined together enabled efficient determination of the CR, but the short-read data alone did not. The B. japonica CR was longer than that of a congenic species inhabiting cooler climate areas, Buergeria buergeri, because of the long repeat regions in the former. During the thermal adaptation of B. japonica, the longer repeat regions in its CR may have accumulated within a period after divergence from B. buergeri.

Genome skimming with nanopore sequencing precisely determines global and transposon DNA methylation in vertebrates

Genome skimming is defined as low-pass sequencing below 0.05× coverage and is typically used for mitochondrial genome recovery and species identification. Long-read nanopore sequencers enable simultaneous reading of both DNA sequence and methylation and can multiplex samples for low-cost genome skimming. Here I present nanopore sequencing as a highly precise platform for global DNA methylation and transposon assessment. At coverage of just 0.001×, or 30 Mb of reads, accuracy is sub-1%. Biological and technical replicates validate high precision. Skimming 40 vertebrate species reveals conserved patterns of global methylation consistent with whole-genome bisulfite sequencing and an average mapping rate >97%. Genome size directly correlates to global DNA methylation, explaining 39% of its variance. Accurate SINE and LINE transposon methylation in both the mouse and primates can be obtained with just 0.0001× coverage, or 3 Mb of reads. Sample multiplexing, field portability, and the low price of this instrument combine to make genome skimming for DNA methylation an accessible method for epigenetic assessment from ecology to epidemiology and for low-resource groups.

How genomics can help biodiversity conservation

The availability of public genomic resources can greatly assist biodiversity assessment, conservation, and restoration efforts by providing evidence for scientifically informed management decisions. Here we survey the main approaches and applications in biodiversity and conservation genomics, considering practical factors, such as cost, time, prerequisite skills, and current shortcomings of applications. Most approaches perform best in combination with reference genomes from the target species or closely related species. We review case studies to illustrate how reference genomes can facilitate biodiversity research and conservation across the tree of life. We conclude that the time is ripe to view reference genomes as fundamental resources and to integrate their use as a best practice in conservation genomics.

Genome Skimming with Nanopore Sequencing Precisely Determines Global and Transposon DNA Methylation in Vertebrates

Genome skimming is defined as low-pass sequencing below 0.05X coverage and is typically used for mitochondrial genome recovery and species identification. Long read nanopore sequencers enable simultaneous reading of both DNA sequence and methylation and can multiplex samples for low-cost genome skimming. Here I present nanopore sequencing as a highly precise platform for global DNA methylation and transposon assessment. At coverage of just 0.001X, or 30 Mb of reads, accuracy is sub-1%. Biological and technical replicates validate high precision. Skimming 40 vertebrate species reveals conserved patterns of global methylation consistent with whole genome bisulfite sequencing and an average mapping rate above 97%. Genome size directly correlates to global DNA methylation, explaining 44% of its variance. Accurate SINE and LINE transposon methylation in both mouse and primates can be obtained with just 0.0001X coverage, or 3 Mb of reads. Sample multiplexing, field portability, and the low price of this instrument combine to make genome skimming for DNA methylation an accessible method for epigenetic assessment from ecology to epidemiology, and by low resource groups.

Forensic Genetic Analyses of Melanistic Iguanas Highlight the Need to Monitor the Iguanid Trade

Simple Summary

The illegal pet trade remains an ongoing, substantial threat to wild populations, especially small insular populations, and can even lead to extinction. Fraudulent activity within the global reptile trade is known to occur, but its identification through forensic applications depends on knowledge of diversity within wild populations. In this study, we assessed the geographic origin of melanistic iguanas (Iguana iguana), which are only found in nations that have never authorized legal export of live animals. Analysis of genetic data from two pet iguanas in the USA flag these as originating from Saba or Montserrat, from which no export permits have ever been issued, confirming their illegal origin. Despite the international trade in I. iguana, in which tens- if not hundreds of thousands of specimens are traded each year, only a handful of individuals have been genetically assessed. Our work highlights the utility of applying forensic genetic techniques to this trade in order to track and discourage illegal activity.

Abstract

Lizards within the Iguana iguana species complex are among the most common reptilian pets, with the widest natural geographic range among iguanids. Deep phylogenetic divergence distinguishes multiple mitochondrial clades, and several taxonomic changes have recently been proposed. These small populations, typically island endemics, are threatened by numerous factors, including the international pet trade. Recent investigations reveal the absence of required CITES permits for lawful export of animals, providing evidence of ongoing illegal trade. Additional monitoring of trade in iguanas can be achieved through the application of forensic molecular techniques. In this study, two captive melanistic iguanas were genotyped for molecular markers for which geographic distributions of alleles have been established. Mitochondrial sequencing indicates that both animals carry a haplotype known to originate from the islands of Saba and Montserrat, populations taxonomically proposed to be Iguana melanoderma. Genotypes at 15 microsatellite loci are equally consistent with this origin, given the results of a principal component analysis. This first forensic genetic assessment within the extensive I. iguana pet trade highlights the presence of illegal activity. The need for additional forensic assessments of pet-trade iguanas is evident, especially given that their value is driven by variety and rarity, which is further intensified by recent taxonomic changes.

Evolution of the Noncoding Features of Sea Snake Mitochondrial Genomes within Elapidae

Mitochondrial genomes of four elapid snakes (three marine species [Emydocephalus ijimae, Hydrophis ornatus, and Hydrophis melanocephalus], and one terrestrial species [Sinomicrurus japonicus]) were completely sequenced by a combination of Sanger sequencing, next-generation sequencing and Nanopore sequencing. Nanopore sequencing was especially effective in accurately reading through long tandem repeats in these genomes. This led us to show that major noncoding regions in the mitochondrial genomes of those three sea snakes contain considerably long tandem duplications, unlike the mitochondrial genomes previously reported for same and other sea snake species. We also found a transposition of the light-strand replication origin within a tRNA gene cluster for the three sea snakes. This change can be explained by the Tandem Duplication—Random Loss model, which was further supported by remnant intervening sequences between tRNA genes. Mitochondrial genomes of true snakes (Alethinophidia) have been shown to contain duplicate major noncoding regions, each of which includes the control region necessary for regulating the heavy-strand replication and transcription from both strands. However, the control region completely disappeared from one of the two major noncoding regions for two Hydrophis sea snakes, posing evolutionary questions on the roles of duplicate control regions in snake mitochondrial genomes. The timing and molecular mechanisms for these changes are discussed based on the elapid phylogeny.

Utilisation of Oxford Nanopore sequencing to generate six complete gastropod mitochondrial genomes as part of a biodiversity curriculum

High-throughput sequencing has enabled genome skimming approaches to produce complete mitochondrial genomes (mitogenomes) for species identification and phylogenomics purposes. In particular, the portable sequencing device from Oxford Nanopore Technologies (ONT) has the potential to facilitate hands-on training from sampling to sequencing and interpretation of mitogenomes. In this study, we present the results from sampling and sequencing of six gastropod mitogenomes (Aplysia argus, Cellana orientalis, Cellana toreuma, Conus ebraeus, Conus miles and Tylothais aculeata) from a graduate level biodiversity course. The students were able to produce mitogenomes from sampling to annotation using existing protocols and programs. Approximately 4 Gb of sequence was produced from 16 Flongle and one MinION flow cells, averaging 235 Mb and N50 = 4.4 kb per flow cell. Five of the six 14.1-18 kb mitogenomes were circlised containing all 13 core protein coding genes. Additional Illumina sequencing revealed that the ONT assemblies spanned over highly AT rich sequences in the control region that were otherwise missing in Illumina-assembled mitogenomes, but still contained a base error of one every 70.8-346.7 bp under the fast mode basecalling with the majority occurring at homopolymer regions. Our findings suggest that the portable MinION device can be used to rapidly produce low-cost mitogenomes onsite and tailored to genomics-based training in biodiversity research.

Are we there yet? Benchmarking low-coverage nanopore long-read sequencing for the assembling of mitochondrial genomes using the vulnerable silky shark Carcharhinus falciformis

Background

Whole mitochondrial genomes are quickly becoming markers of choice for the exploration of within-species genealogical and among-species phylogenetic relationships. Most often, ‘primer walking’ or ‘long PCR’ strategies plus Sanger sequencing or low-pass whole genome sequencing using Illumina short reads are used for the assembling of mitochondrial chromosomes. In this study, we first confirmed that mitochondrial genomes can be sequenced from long reads using nanopore sequencing data exclusively. Next, we examined the accuracy of the long-reads assembled mitochondrial chromosomes when comparing them to a ‘gold’ standard reference mitochondrial chromosome assembled using Illumina short-reads sequencing.

Results

Using a specialized bioinformatics tool, we first produced a short-reads mitochondrial genome assembly for the silky shark C. falciformis with an average base coverage of 9.8x. The complete mitochondrial genome of C. falciformis was 16,705 bp in length and 934 bp shorter than a previously assembled genome (17,639 bp in length) that used bioinformatics tools not specialized for the assembly of mitochondrial chromosomes. Next, low-pass whole genome sequencing using a MinION ONT pocket-sized platform plus customized de-novo and reference-based workflows assembled and circularized a highly accurate mitochondrial genome in the silky shark Carcharhinus falciformis. Indels at the flanks of homopolymer regions explained most of the dissimilarities observed between the ‘gold’ standard reference mitochondrial genome (assembled using Illumina short reads) and each of the long-reads mitochondrial genome assemblies. Although not completely accurate, mitophylogenomics and barcoding analyses (using entire mitogenomes and the D-Loop/Control Region, respectively) suggest that long-reads assembled mitochondrial genomes are reliable for identifying a sequenced individual, such as C. falciformis, and separating the same individual from others belonging to closely related congeneric species.

Conclusions

This study confirms that mitochondrial genomes can be sequenced from long-reads nanopore sequencing data exclusively. With further development, nanopore technology can be used to quickly test in situ mislabeling in the shark fin fishing industry and thus, improve surveillance protocols, law enforcement, and the regulation of this fishery. This study will also assist with the transferring of high-throughput sequencing technology to middle- and low-income countries so that international scientists can explore population genomics in sharks using inclusive research strategies. Lastly, we recommend assembling mitochondrial genomes using specialized assemblers instead of other assemblers developed for bacterial and/or nuclear genomes.

A benchmarking of human mitochondrial DNA haplogroup classifiers from whole-genome and whole-exome sequence data

The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. We also assessed the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short-read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

Tracing Viral Transmission and Evolution of Bovine Leukemia Virus through Long Read Oxford Nanopore Sequencing of the Proviral Genome

Bovine leukemia virus (BLV) causes Enzootic Bovine Leukosis (EBL), a persistent life-long disease resulting in immune dysfunction and shortened lifespan in infected cattle, severely impacting the profitability of the US dairy industry. Our group has found that 94% of dairy farms in the United States are infected with BLV with an average in-herd prevalence of 46%. This is partly due to the lack of clinical presentation during the early stages of primary infection and the elusive nature of BLV transmission. This study sought to validate a near-complete genomic sequencing approach for reliability and accuracy before determining its efficacy in characterizing the sequence identity of BLV proviral genomes collected from a pilot study made up of 14 animals from one commercial dairy herd. These BLV-infected animals were comprised of seven adult dam/daughter pairs that tested positive by ELISA and qPCR. The results demonstrate sequence identity or divergence of the BLV genome from the same samples tested in two independent laboratories, suggesting both vertical and horizontal transmission in this dairy herd. This study supports the use of Oxford Nanopore sequencing for the identification of viral SNPs that can be used for retrospective genetic contact tracing of BLV transmission.

Genomic skimming and nanopore sequencing uncover cryptic hybridization in one of world's most threatened primates

The Brazilian buffy-tufted-ear marmoset (Callithrix aurita), one of the world's most endangered primates, is threatened by anthropogenic hybridization with exotic, invasive marmoset species. As there are few genetic data available for C. aurita, we developed a PCR-free protocol with minimal technical requirements to rapidly generate genomic data with genomic skimming and portable nanopore sequencing. With this direct DNA sequencing approach, we successfully determined the complete mitogenome of a marmoset that we initially identified as C. aurita. The obtained nanopore-assembled sequence was highly concordant with a Sanger sequenced version of the same mitogenome. Phylogenetic analyses unexpectedly revealed that our specimen was a cryptic hybrid, with a C. aurita phenotype and C. penicillata mitogenome lineage. We also used publicly available mitogenome data to determine diversity estimates for C. aurita and three other marmoset species. Mitogenomics holds great potential to address deficiencies in genomic data for endangered, non-model species such as C. aurita. However, we discuss why mitogenomic approaches should be used in conjunction with other data for marmoset species identification. Finally, we discuss the utility and implications of our results and genomic skimming/nanopore approach for conservation and evolutionary studies of C. aurita and other marmosets.

A Rapid and Cost-Effective Identification of Invertebrate Pests at the Borders Using MinION Sequencing of DNA Barcodes

The rapid and accurate identification of invertebrate pests detected at the border is a challenging task. Current diagnostic methods used at the borders are mainly based on time consuming visual and microscopic examinations. Here, we demonstrate a rapid in-house workflow for DNA extraction, PCR amplification of the barcode region of the mitochondrial cytochrome oxidase subunit I (COI) gene and Oxford Nanopore Technologies (ONT) MinION sequencing of amplified products multiplexed after barcoding on ONT Flongle flow cells. A side-by-side comparison was conducted of DNA barcode sequencing-based identification and morphological identification of both large (>0.5 mm in length) and small (<0.5 mm in length) invertebrate specimens intercepted at the Australian border. DNA barcode sequencing results supported the morphological identification in most cases and enabled immature stages of invertebrates and their eggs to be identified more confidently. Results also showed that sequencing the COI barcode region using the ONT rapid sequencing principle is a cost-effective and field-adaptable approach for the rapid and accurate identification of invertebrate pests. Overall, the results suggest that MinION sequencing of DNA barcodes offers a complementary tool to the existing morphological diagnostic approaches and provides rapid, accurate, reliable and defendable evidence for identifying invertebrate pests at the border.

Yes, we can use it: a formal test on the accuracy of low-pass nanopore long-read sequencing for mitophylogenomics and barcoding research using the Caribbean spiny lobster Panulirus argus

Background

Whole mitogenomes or short fragments (i.e., 300–700 bp of the cox1 gene) are the markers of choice for revealing within- and among-species genealogies. Protocols for sequencing and assembling mitogenomes include ‘primer walking’ or ‘long PCR’ followed by Sanger sequencing or Illumina short-read low-coverage whole genome (LC-WGS) sequencing with or without prior enrichment of mitochondrial DNA. The aforementioned strategies assemble complete and accurate mitochondrial genomes but are time consuming and/or expensive. In this study, I first tested whether mitogenomes can be sequenced from long-read nanopore sequencing data exclusively. Second, I explored the accuracy of the long-read assembled genomes by comparing them to a ‘gold’ standard reference mitogenome retrieved from the same individual using Illumina sequencing. Third and lastly, I tested if the long-read assemblies are useful for mitophylogenomics and barcoding research. To accomplish these goals, I used the Caribbean spiny lobster Panulirus argus, an ecologically relevant species in shallow water coral reefs and target of the most lucrative fishery in the greater Caribbean region.

Results

LC-WGS using a MinION ONT device and various de-novo and reference-based assembly pipelines retrieved a complete and highly accurate mitogenome for the Caribbean spiny lobster Panulirus argus. Discordance between each of the long-read assemblies and the reference mitogenome was mostly due to indels at the flanks of homopolymer regions. Although not ‘perfect’, phylogenetic analyses using entire mitogenomes or a fragment of the cox1 gene demonstrated that mitogenomes assembled using long reads reliably identify the sequenced specimen as belonging to P. argus and distinguish it from other related species in the same genus, family, and superorder.

Conclusions

This study serves as a proof-of-concept for the future implementation of in-situ surveillance protocols using the MinION to detect mislabeling in P. argus across its supply chain. Mislabeling detection will improve fishery management in this overexploited lobster. This study will additionally aid in decreasing costs for exploring meta-population connectivity in the Caribbean spiny lobster and will aid with the transfer of genomics technology to low-income countries.