Comparing the performance of multiple mitochondrial genes in the analysis of Australian freshwater fishes

Identification of non-model mammal species using the MinION DNA sequencer from Oxford Nanopore

Background

The Neotropics harbors the largest species richness of the planet; however, even in well-studied groups, there are potentially hundreds of species that lack a formal description, and likewise, many already described taxa are difficult to identify using morphology. Specifically in small mammals, complex morphological diagnoses have been facilitated by the use of molecular data, particularly from mitochondrial sequences, to obtain accurate species identifications. Obtaining mitochondrial markers implies the use of PCR and specific primers, which are largely absent for non-model organisms. Oxford Nanopore Technologies (ONT) is a new alternative for sequencing the entire mitochondrial genome without the need for specific primers. Only a limited number of studies have employed exclusively ONT long-reads to assemble mitochondrial genomes, and few studies have yet evaluated the usefulness of such reads in multiple non-model organisms.

Methods

We implemented fieldwork to collect small mammals, including rodents, bats, and marsupials, in five localities in the northern extreme of the Cordillera Central of Colombia. DNA samples were sequenced using the MinION device and Flongle flow cells. Shotgun-sequenced data was used to reconstruct the mitochondrial genome of all the samples. In parallel, using a customized computational pipeline, species-level identifications were obtained based on sequencing raw reads (Whole Genome Sequencing). ONT-based identifications were corroborated using traditional morphological characters and phylogenetic analyses.

Results

A total of 24 individuals from 18 species were collected, morphologically identified, and deposited in the biological collection of Universidad EAFIT. Our different computational pipelines were able to reconstruct mitochondrial genomes from exclusively ONT reads. We obtained three new mitochondrial genomes and eight new molecular mitochondrial sequences for six species. Our species identification pipeline was able to obtain accurate species identifications for up to 75% of the individuals in as little as 5 s. Finally, our phylogenetic analyses corroborated the identifications from our automated species identification pipeline and revealed important contributions to the knowledge of the diversity of Neotropical small mammals.

Discussion

This study was able to evaluate different pipelines to reconstruct mitochondrial genomes from non-model organisms, using exclusively ONT reads, benchmarking these protocols on a multi-species dataset. The proposed methodology can be applied by non-expert taxonomists and has the potential to be implemented in real-time, without the need to euthanize the organisms and under field conditions. Therefore, it stands as a relevant tool to help increase the available data for non-model organisms, and the rate at which researchers can characterize life specially in highly biodiverse places as the Neotropics.

Quality of DNA extracted from freshwater fish scales and mucus and its application in genetic diversity studies of Perca fluviatilis and Rutilus rutilus

Studies on genetic diversity require biological material containing a reliable source of DNA that can be extracted and analyzed. Recently, non-invasive sampling has become a preferred sampling method of biological material. The suitability of a less invasive approach that involves obtaining samples by swabbing the fish skin (including live, non-anesthetized fish) should be considered. In this study, we compared the efficiency of DNA extraction, amplification, and sequencing of mtDNA fragments of two fish species Perca fluviatilis and Rutilus rutilus based on DNA collected from the scales and mucus using the modified Aljanabi and Martinez method. The results revealed a higher quality of DNA extracted from the mucus; however, the mean DNA concentration obtained from the scales of both fish species was higher. We verified the method suitable for amplification and sequencing of mtDNA fragments of both fish species using newly designed markers (D-loop, ATP6) and examined the potential risk of intraspecific cross-contamination. The DNA sequence alignment analysis revealed identical sequences attributed to the same individual when DNA, extracted from two different sources (scales and mucus), was used. We demonstrated that the quantity and quality of DNA extracted from the scales and mucus using the proposed method were high enough to carry out genetic diversity studies based on sampling of live fish with the possibility to release it after collecting samples.

Cenozoic colonisation of the Indian Ocean region by the Australian freshwater originating glassperch family Ambassidae (Teleostei)

We examined the phylogeny and biogeography of the glassperch family Ambassidae (Teleostei), which is widely distributed in the freshwater, brackish and marine coastal habitats across the Indo-West Pacific region. We first built a comprehensive time-calibrated phylogeny of Ambassidae using five genes. We then used this tree to reconstruct the evolution of the salinity preference and ancestral areas. Our results indicate that the two largest genera of Ambassidae, Ambassis and Parambassis, are each not monophyletic. The most recent common ancestor of Ambassidae was freshwater adapted and lived in Australia about 56 million years ago. Three independent freshwater-to-marine transitions are inferred, but no marine-to-freshwater ones. To explain the distribution of ambassids, we hypothesise two long-distance marine dispersal events from Australia. A first event was towards Southeast Asia during the early Cenozoic, followed by a second one towards Africa during mid-Cenozoic. The phylogenetic signal associated with the salinity adaptation of these events was not detected, possibly because of the selective extinction of intermediate marine lineages. The Ambassidae shares two characteristics with other freshwater fish groups distributed in continental regions surrounding the Indian Ocean: They are too young to support the hypothesis that their distribution is the result of the fragmentation of Gondwana, but they did not retain the phylogenetic signal of their marine dispersal.

Ancient DNA from the extinct New Zealand grayling (Prototroctes oxyrhynchus) reveals evidence for Miocene marine dispersal

The evolutionary history of Southern Hemisphere graylings (Retropinnidae) in New Zealand (NZ), including their relationship to the Australian grayling, is poorly understood. The NZ grayling (Prototroctes oxyrhynchus) is the only known fish in NZ to have gone extinct since human arrival there. Despite its historical abundance, only 23 wet and dried, formalin-fixed specimens exist in museums. We used high-throughput DNA sequencing to generate mitogenomes from formalin-fixed P. oxyrhynchus specimens, and analysed these in a temporal phylogenetic framework of retropinnids and osmerids. We recovered a strong sister-relationship between NZ and Australian grayling (P. mareana), with a common ancestor ~13.8 Mya [95% highest posterior density (HPD): 6.1–23.2 Mya], after the height of Oligocene marine inundation in NZ. Our temporal phylogenetic analysis suggests a single marine dispersal between NZ and Australia, although the direction of dispersal is equivocal, followed by divergence into genetically and morphologically distinguishable species through isolation by distance. This study provides further insights into the possible extinction drivers of the NZ grayling, informs discussion regarding reintroduction of Prototroctes to NZ and highlights how advances in palaeogenetics can be used to test evolutionary hypotheses in fish, which, until relatively recently, have been comparatively neglected in ancient-DNA research.

Use of complete mitochondrial genome sequences to identify barcoding markers for groups with low genetic distance

Complete mitochondrial sequences can be rapidly obtained and are widely available, providing a great source of species information and allowing for the discovery of new specific molecular markers. However, for some taxonomic groups, traditional approaches for species delimitation are impaired by the low genetic distance values. In these cases, other species-level markers are used. For Prochilodus, which includes important neotropical fish species, species-level delimitation usually results in poor phylogenetic resolution when using mitochondrial COI/cytB genes as barcoding markers because of low genetic variability and low species-level resolution. Thus, in this study, we developed an approach to design and validate new barcoding markers with high species-level resolution obtained from the D-loop region, using Prochilodus spp. as a model. For the new barcoding marker validation, the amplicon region was used to infer the phylogenetic relationships of Prochilodus spp. through three distinct methods: Bayesian inference (BI), Neighbor-Joining method (NJ), and Maximum Likelihood method (ML). The phylogenetic relationships of Prochilodus spp. revealed high resolution at species-level, nonoverlapping clades, and high branch support. The genetic distance results allied to two different clustering methods (Bayesian Poisson tree processes and automatic barcode gap discovery) revealed the existence of a barcoding gap, thus, validating the use of the barcoding markers designed in this study. The approach proposed here may, therefore, be expanded to other taxa to access and validate new barcoding markers with higher resolution at the species level.

Multiple molecular markers reinforce the systematic framework of unique Australian cave fishes (Milyeringa : Gobioidei)

Australia was once thought to be a biodiversity desert when considering the subterranean world; however, recent work has revealed a fascinating collection of cave creatures, many with surprising biogeographic histories. This has especially been so in the karstic regions of north-western Australia (Cape Range peninsula, Barrow Island, Pilbara), which is home not only to a diverse collection of subterranean invertebrates, but also to the continent’s only known underworld-adapted vertebrates, which includes the cave fish in the genus Milyeringa. These cave gudgeons have recently been in a state of taxonomic flux, with species being both split and lumped, but this was done with limited data (incomplete geographic sampling and no nuclear DNA sequence data). Therefore, we have revisited the systematic status of Milyeringa in a total-evidence molecular approach by integrating all existing data (mitochondrial, allozymes) with new DNA sequences from nuclear and mitochondrial loci and new multilocus allozyme data. Our conclusion, that there are two species, matches the most recent taxonomic treatment, with Milyeringa veritas present on both the eastern and western sides of the Cape Range peninsula, and Milyeringa justitia on Barrow Island. This has implications for future research in the linked fields of biogeography and conservation.

DNA barcoding of freshwater fishes and the development of a quantitative qPCR assay for the species-specific detection and quantification of fish larvae from plankton samples

The barcoding of mitochondrial cytochrome c oxidase subunit 1 (coI) gene was amplified and sequenced from 16 species of freshwater fishes found in Lake Wivenhoe (south-eastern Queensland, Australia) to support monitoring of reservoir fish populations, ecosystem function and water health. In this study, 630-650 bp sequences of the coI barcoding gene from 100 specimens representing 15 genera, 13 families and two subclasses of fishes allowed 14 of the 16 species to be identified and differentiated. The mean ± s.e. Kimura 2 parameter divergence within and between species was 0.52 ± 0.10 and 23.8 ± 2.20% respectively, indicating that barcodes can be used to discriminate most of the fish species accurately. The two terapontids, Amniataba percoides and Leiopotherapon unicolor, however, shared coI DNA sequences and could not be differentiated using this gene. A barcoding database was established and a qPCR assay was developed using coI sequences to identify and quantify proportional abundances of fish species in ichthyoplankton samples from Lake Wivenhoe. These methods provide a viable alternative to the time-consuming process of manually enumerating and identifying ichthyoplankton samples.

The past and present of an estuarine-resident fish, the "four-eyed fish" Anableps anableps (Cyprinodontiformes, Anablepidae), revealed by mtDNA sequences

Historical events, such as changes in sea level during the Pleistocene glacial cycles, had a strong impact on coastal habitats, limiting connectivity and promoting the genetic divergence of various species. In this study, we evaluated the influence of climate oscillations and the possibility of estuary function as a barrier to gene flow among populations of the four-eyed fish, Anableps anableps. This species is fully estuarine-resident, has internal fertilization, is viviparous and does not migrate across long distances. These features make the four-eyed fish an excellent model for the study of evolutionary processes related to genetic differentiation of species and populations in estuaries. The evolutionary history of A. anableps was inferred from phylogeographic and population analyses using sequences of the mitochondrial DNA Control Region of 13 populations distributed in the Amazon and Northeast Coast of Brazil from Calcoene (Amapa) to Parnaiba (Piaui). The 83 retrieved haplotypes show a pattern of four distinct mitochondrial lineages, with up to 3.4% nucleotide divergence among them. The evolutionary reconstruction suggests that these lineages diverged recently in the late Pleistocene/early Holocene after the Atlantic Ocean reaching current levels. Analysis of variability, neutrality and the genetic expansion pattern revealed that the lineages have distinct characteristics, which were shaped by the different geomorphological features of coastal regions combined with sea level oscillations over a very long period of time. Only few neighboring populations show a discreet gene flow. This study may also be helpful for designing new experiments to better understand the geomorphological evolutionary history of the estuaries of the Amazon and the Northeast Coast of Brazil using estuarine-resident species as a model.

Low interbasin connectivity in a facultatively diadromous fish: evidence from genetics and otolith chemistry

Southern smelts (Retropinna spp.) in coastal rivers of Australia are facultatively diadromous, with populations potentially containing individuals with diadromous or wholly freshwater life histories. The presence of diadromous individuals is expected to reduce genetic structuring between river basins due to larval dispersal via the sea. We use otolith chemistry to distinguish between diadromous and nondiadromous life histories and population genetics to examine interbasin connectivity resulting from diadromy. Otolith strontium isotope ((87) Sr:(86) Sr) transects identified three main life history patterns: amphidromy, freshwater residency and estuarine/marine residency. Despite the potential for interbasin connectivity via larval mixing in the marine environment, we found unprecedented levels of genetic structure for an amphidromous species. Strong hierarchical structure along putative taxonomic boundaries was detected, along with highly structured populations within groups using microsatellites (FST  = 0.046-0.181), and mtDNA (ΦST  = 0.498-0.816). The presence of strong genetic subdivision, despite the fact that many individuals reside in saline water during their early life history, appears incongruous. However, analysis of multielemental signatures in the otolith cores of diadromous fish revealed strong discrimination between river basins, suggesting that diadromous fish spend their early lives within chemically distinct estuaries rather than the more homogenous marine environment, thus avoiding dispersal and maintaining genetic structure.