Matthew S. Leslie

Targeted multiplex next-generation sequencing: advances in techniques of mitochondrial and nuclear DNA sequencing for population genomics

Next‐generation sequencing (NGS) is emerging as an efficient and cost‐effective tool in population genomic analyses of nonmodel organisms, allowing simultaneous resequencing of many regions of multi‐genomic DNA from multiplexed samples. Here, we detail our synthesis of protocols for targeted resequencing of mitochondrial and nuclear loci by generating indexed genomic libraries for multiplexing up to 100 individuals in a single sequencing pool, and then enriching the pooled library using custom DNA capture arrays. Our use of DNA sequence from one species to capture and enrich the sequencing libraries of another species (i.e. cross‐species DNA capture) indicates that efficient enrichment occurs when sequences are up to about 12% divergent, allowing us to take advantage of genomic information in one species to sequence orthologous regions in related species. In addition to a complete mitochondrial genome on each array, we have included between 43 and 118 nuclear loci for low‐coverage sequencing of between 18 kb and 87 kb of DNA sequence per individual for single nucleotide polymorphisms discovery from 50 to 100 individuals in a single sequencing lane. Using this method, we have generated a total of over 500 whole mitochondrial genomes from seven cetacean species and green sea turtles. The greater variation detected in mitogenomes relative to short mtDNA sequences is helping to resolve genetic structure ranging from geographic to species‐level differences. These NGS and analysis techniques have allowed for simultaneous population genomic studies of mtDNA and nDNA with greater genomic coverage and phylogeographic resolution than has previously been possible in marine mammals and turtles.

Mitogenomic phylogenetics of fin whales (Balaenoptera physalus spp.): genetic evidence for revision of subspecies.

There are three described subspecies of fin whales (Balaenoptera physalus): B. p. physalus Linnaeus, 1758 in the Northern Hemisphere, B. p. quoyi Fischer, 1829 in the Southern Hemisphere, and a recently described pygmy form, B. p. patachonica Burmeister, 1865. The discrete distribution in the North Pacific and North Atlantic raises the question of whether a single Northern Hemisphere subspecies is valid. We assess phylogenetic patterns using ∼16 K base pairs of the complete mitogenome for 154 fin whales from the North Pacific, North Atlantic - including the Mediterranean Sea - and Southern Hemisphere. A Bayesian tree of the resulting 136 haplotypes revealed several well-supported clades representing each ocean basin, with no haplotypes shared among ocean basins. The North Atlantic haplotypes (n = 12) form a sister clade to those from the Southern Hemisphere (n = 42). The estimated time to most recent common ancestor (TMRCA) for this Atlantic/Southern Hemisphere clade and 81 of the 97 samples from the North Pacific was approximately 2 Ma. 14 of the remaining North Pacific samples formed a well-supported clade within the Southern Hemisphere. The TMRCA for this node suggests that at least one female from the Southern Hemisphere immigrated to the North Pacific approximately 0.37 Ma. These results provide strong evidence that North Pacific and North Atlantic fin whales should not be considered the same subspecies, and suggest the need for revision of the global taxonomy of the species.

Using Genome-Wide SNPs to Detect Structure in High-Diversity and Low-Divergence Populations of Severely Impacted Eastern Tropical Pacific Spinner (Stenella longirostris) and Pantropical Spotted Dolphins (S. attenuata)

Millions of spinner (Stenella longirostris) and pantropical spotted dolphins (Stenella attenuata) died since the 1960’s as bycatch in tuna nets in the eastern tropical Pacific Ocean. Despite three decades of protection, they show little-to-no sign of recovery (although recent fisheries-independent abundance estimates are not available). In efforts to establish biologically meaningful management boundaries for recovery, endemic subspecies and multiple stocks have been proposed. However, genetic differentiation among most of these units has been difficult to identify, possibly due to low statistical power stemming from large historical abundances, ongoing gene flow, and recent divergence. We tested for genetic structure at multiple hierarchical levels by analyzing the largest dataset to date brought to bear on these questions. Single nucleotide polymorphisms (SNPs) were collected from nuclear DNA regions associated with the restriction enzyme site PstI from 72 spinner dolphins and 58 pantropical spotted dolphins using genotype-by-sequencing (GBS). Our results support the current subspecies for both species and indicate stock-level separation for Tres Marias spinner dolphins and the two offshore pantropical spotted dolphin stocks in this area. Although bycatch has been reduced a small fraction of pre-protection levels, incidental mortality continues to impact these populations. Our results are important for the ongoing management and recovery of these highly-impacted pelagic dolphins in the eastern tropical Pacific Ocean.

Structure and phylogeography of two tropical predators, spinner (Stenella longirostris) and pantropical spotted (S. attenuata) dolphins, from SNP data

Little is known about global patterns of genetic connectivity in pelagic dolphins, including how circumtropical pelagic dolphins spread globally following the rapid and recent radiation of the subfamily delphininae. In this study, we tested phylogeographic hypotheses for two circumtropical species, the spinner dolphin (Stenella longirostris) and the pantropical spotted dolphin (Stenella attenuata), using more than 3000 nuclear DNA single nucleotide polymorphisms (SNPs) in each species. Analyses for population structure indicated significant genetic differentiation between almost all subspecies and populations in both species. Bayesian phylogeographic analyses of spinner dolphins showed deep divergence between Indo-Pacific, Atlantic and eastern tropical Pacific Ocean (ETP) lineages. Despite high morphological variation, our results show very close relationships between endemic ETP spinner subspecies in relation to global diversity. The dwarf spinner dolphin is a monophyletic subspecies nested within a major clade of pantropical spinner dolphins from the Indian and western Pacific Ocean populations. Population-level division among the dwarf spinner dolphins was detected—with the northern Australia population being very different from that in Indonesia. In contrast to spinner dolphins, the major boundary for spotted dolphins is between offshore and coastal habitats in the ETP, supporting the current subspecies-level taxonomy. Comparing these species underscores the different scale at which population structure can arise, even in species that are similar in habitat (i.e. pelagic) and distribution.

Mitogenomic differentiation in spinner (Stenella longirostris) and pantropical spotted dolphins (S. attenuata) from the eastern tropical Pacific Ocean

Endemic subspecies and multiple stocks have been proposed for spinner (Stenella longirostris) and spotted dolphins (Stenella attenuata), two species with historically high mortality due to bycatch in the eastern tropical Pacific yellowfin tuna fishery. However, there has been low statistical power in tests for genetic differentiation among most of these stocks, due to large historical abundances, ongoing gene flow, and recent divergence. We tested for structure at multiple hierarchical levels by collecting whole mitochondrial genome sequences (mtDNA) and nuclear SNPs (nuDNA) from 104 spinner and 76 spotted dolphins using capture array library enrichment and highly paralleled DNA sequencing. MtDNA showed weak but significant differences between subspecies of spotted (FST: 0.0125; P = 0.0402) and spinner dolphins (FST: 0.0133; P = 0.034). NuDNA supported subspecies of spotted but not spinner dolphins. Relatively strong and significant differentiation was detected between whitebelly and eastern spinner stocks using nuDNA (FST: 0.0297; P = 0.0059). Neither mtDNA nor nuDNA supported the division of existing offshore stocks of spotted dolphins or Tres Marias spinner dolphins. This work identifies a genetic basis for biologically meaningful management units of these two species, a critical component in understanding their response in the face of historical and continued fishery interactions.

Impacts of phylogenetic nomenclature on the efficacy of the U.S. Endangered Species Act

Cataloging biodiversity is critical to conservation efforts because accurate taxonomy is often a precondition for protection under laws designed for species conservation, such as the U.S. Endangered Species Act (ESA). Traditional nomenclatural codes governing the taxonomic process have recently come under scrutiny because taxon names are more closely linked to hierarchical ranks than to the taxa themselves. A new approach to naming biological groups, called phylogenetic nomenclature (PN), explicitly names taxa by defining their names in terms of ancestry and descent. PN has the potential to increase nomenclatural stability and decrease confusion induced by the rank-based codes. But proponents of PN have struggled with whether species and infraspecific taxa should be governed by the same rules as other taxa or should have special rules. Some proponents advocate the wholesale abandonment of rank labels (including species); this could have consequences for the implementation of taxon-based conservation legislation. I examined the principles of PN as embodied in the PhyloCode (an alternative to traditional rank-based nomenclature that names biological groups based on the results of phylogenetic analyses and does not associate taxa with ranks) and assessed how this novel approach to naming taxa might affect the implementation of species-based legislation by providing a case study of the ESA. The latest version of the PhyloCode relies on the traditional rank-based codes to name species and infraspecific taxa; thus, little will change regarding the main targets of the ESA because they will retain rank labels. For this reason, and because knowledge of evolutionary relationships is of greater importance than nomenclatural procedures for initial protection of endangered taxa under the ESA, I conclude that PN under the PhyloCode will have little impact on implementation of the ESA.