Revised taxonomy of eastern North Pacific killer whales (Orcinus orca): Bigg's and resident ecotypes deserve species status.

Killer whales (Orcinus orca) are currently recognized as a single ecologically and morphologically diverse, globally distributed species. Multiple morphotypes or ecotypes have been described, often associated with feeding specialization, and several studies have suggested taxonomic revision to include multiple subspecies or species in the genus. We review the ecological, morphological and genetic data for the well-studied 'resident' and Bigg's (aka 'transient') ecotypes in the eastern North Pacific and use quantitative taxonomic guidelines and standards to determine whether the taxonomic status of these killer whale ecotypes should be revised. Our review and new analyses indicate that species-level status is justified in both cases, and we conclude that eastern North Pacific Bigg's killer whales should be recognized as Orcinus rectipinnus (Cope in Scammon, 1869) and resident killer whales should be recognized as Orcinus ater (Cope in Scammon, 1869).

The genomic footprint of whaling and isolation in fin whale populations.

Twentieth century industrial whaling pushed several species to the brink of extinction, with fin whales being the most impacted. However, a small, resident population in the Gulf of California was not targeted by whaling. Here, we analyzed 50 whole-genomes from the Eastern North Pacific (ENP) and Gulf of California (GOC) fin whale populations to investigate their demographic history and the genomic effects of natural and human-induced bottlenecks. We show that the two populations diverged ~16,000 years ago, after which the ENP population expanded and then suffered a 99% reduction in effective size during the whaling period. In contrast, the GOC population remained small and isolated, receiving less than one migrant per generation. However, this low level of migration has been crucial for maintaining its viability. Our study exposes the severity of whaling, emphasizes the importance of migration, and demonstrates the use of genome-based analyses and simulations to inform conservation strategies.

Diversity of mitochondrial DNA in 3 species of great whales before and after modern whaling.

The 20th century commercial whaling industry severely reduced populations of great whales throughout the Southern Hemisphere. The effect of this exploitation on genetic diversity and population structure remains largely undescribed. Here, we compare pre- and post-whaling diversity of mitochondrial DNA (mtDNA) control region sequences for 3 great whales in the South Atlantic, such as the blue, humpback, and fin whale. Pre-whaling diversity is described from mtDNA extracted from bones collected near abandoned whaling stations, primarily from the South Atlantic island of South Georgia. These bones are known to represent the first stage of 20th century whaling and thus pre-whaling diversity of these populations. Post-whaling diversity is described from previously published studies reporting large-scale sampling of living whales in the Southern Hemisphere. Despite relatively high levels of surviving genetic diversity in the post-whaling populations, we found evidence of a probable loss of mtDNA lineages in all 3 species. This is evidenced by the detection of a large number of haplotypes found in the pre-whaling samples that are not present in the post-whaling samples. A rarefaction analysis further supports a loss of haplotypes in the South Atlantic humpback and Antarctic blue whale populations. The bones from former whaling stations in the South Atlantic represent a remarkable molecular archive for further investigation of the decline and ongoing recovery in the great whales of the Southern Hemisphere.

Guidelines and quantitative standards for improved cetacean taxonomy using full mitochondrial genomes.

In many organisms, especially those of conservation concern, traditional lines of evidence for taxonomic delineation, such as morphological data, are often difficult to obtain. In these cases, genetic data are often the only source of information available for taxonomic studies. In particular, population surveys of mitochondrial genomes offer increased resolution and precision in support of taxonomic decisions relative to conventional use of the control region or other gene fragments of the mitochondrial genome. To improve quantitative guidelines for taxonomic decisions in cetaceans, we build on a previous effort targeting the control region and evaluate, for whole mitogenome sequences, a suite of divergence and diagnosability estimates for pairs of recognized cetacean populations, subspecies, and species. From this overview, we recommend new guidelines based on complete mitogenomes, combined with other types of evidence for isolation and divergence, which will improve resolution for taxonomic decisions, especially in the face of small sample sizes or low levels of genetic diversity. We further use simulated data to assist interpretations of divergence in the context of varying forms of historical demography, culture, and ecology.

Possible niche compression and individual specialization in Pacific Arctic beluga (Delphinapterus leucas) from the 19th to 20th century.

Cetaceans have shown a potential to be used as sentinel species for tracking environmental change in marine ecosystems, yet our assessment of change is typically limited to recent decades and lacks ecological baselines. Using historical museum specimens, we compared community niche metrics and degree of individual dietary specialization in groups of Pacific Arctic beluga (Delphinapterus leucas) from the 1800s (n = 5) to 1900s (n = 10) using stable carbon and nitrogen isotopes drilled from teeth. Beluga occupied a broader trophic niche and demonstrated a higher degree of individual specialization in the 1800s than in the 1900s. The cause of this shift is difficult to confirm given long timescales and constraints of specimen-based research but could indicate changes in the prey base or competition. The scale and nature of this detected shift provide perspective for continued research on these climate-vulnerable species.

Historical Mitogenomic Diversity and Population Structuring of Southern Hemisphere Fin Whales.

Fin whales Balaenoptera physalus were hunted unsustainably across the globe in the 19th and 20th centuries, leading to vast reductions in population size. Whaling catch records indicate the importance of the Southern Ocean for this species; approximately 730,000 fin whales were harvested during the 20th century in the Southern Hemisphere (SH) alone, 94% of which were at high latitudes. Genetic samples from contemporary whales can provide a window to past population size changes, but the challenges of sampling in remote Antarctic waters limit the availability of data. Here, we take advantage of historical samples in the form of bones and baleen available from ex-whaling stations and museums to assess the pre-whaling diversity of this once abundant species. We sequenced 27 historical mitogenomes and 50 historical mitochondrial control region sequences of fin whales to gain insight into the population structure and genetic diversity of Southern Hemisphere fin whales (SHFWs) before and after the whaling. Our data, both independently and when combined with mitogenomes from the literature, suggest SHFWs are highly diverse and may represent a single panmictic population that is genetically differentiated from Northern Hemisphere populations. These are the first historic mitogenomes available for SHFWs, providing a unique time series of genetic data for this species.

Global genetic diversity status and trends: towards a suite of Essential Biodiversity Variables (EBVs) for genetic composition.

Biodiversity underlies ecosystem resilience, ecosystem function, sustainable economies, and human well-being. Understanding how biodiversity sustains ecosystems under anthropogenic stressors and global environmental change will require new ways of deriving and applying biodiversity data. A major challenge is that biodiversity data and knowledge are scattered, biased, collected with numerous methods, and stored in inconsistent ways. The Group on Earth Observations Biodiversity Observation Network (GEO BON) has developed the Essential Biodiversity Variables (EBVs) as fundamental metrics to help aggregate, harmonize, and interpret biodiversity observation data from diverse sources. Mapping and analyzing EBVs can help to evaluate how aspects of biodiversity are distributed geographically and how they change over time. EBVs are also intended to serve as inputs and validation to forecast the status and trends of biodiversity, and to support policy and decision making. Here, we assess the feasibility of implementing Genetic Composition EBVs (Genetic EBVs), which are metrics of within-species genetic variation. We review and bring together numerous areas of the field of genetics and evaluate how each contributes to global and regional genetic biodiversity monitoring with respect to theory, sampling logistics, metadata, archiving, data aggregation, modeling, and technological advances. We propose four Genetic EBVs: (i) Genetic Diversity; (ii) Genetic Differentiation; (iii) Inbreeding; and (iv) Effective Population Size (Ne ). We rank Genetic EBVs according to their relevance, sensitivity to change, generalizability, scalability, feasibility and data availability. We outline the workflow for generating genetic data underlying the Genetic EBVs, and review advances and needs in archiving genetic composition data and metadata. We discuss how Genetic EBVs can be operationalized by visualizing EBVs in space and time across species and by forecasting Genetic EBVs beyond current observations using various modeling approaches. Our review then explores challenges of aggregation, standardization, and costs of operationalizing the Genetic EBVs, as well as future directions and opportunities to maximize their uptake globally in research and policy. The collection, annotation, and availability of genetic data has made major advances in the past decade, each of which contributes to the practical and standardized framework for large-scale genetic observation reporting. Rapid advances in DNA sequencing technology present new opportunities, but also challenges for operationalizing Genetic EBVs for biodiversity monitoring regionally and globally. With these advances, genetic composition monitoring is starting to be integrated into global conservation policy, which can help support the foundation of all biodiversity and species' long-term persistence in the face of environmental change. We conclude with a summary of concrete steps for researchers and policy makers for advancing operationalization of Genetic EBVs. The technical and analytical foundations of Genetic EBVs are well developed, and conservation practitioners should anticipate their increasing application as efforts emerge to scale up genetic biodiversity monitoring regionally and globally.

Charting a course for genetic diversity in the UN Decade of Ocean Science.

The health of the world's oceans is intrinsically linked to the biodiversity of the ecosystems they sustain. The importance of protecting and maintaining ocean biodiversity has been affirmed through the setting of the UN Sustainable Development Goal 14 to conserve and sustainably use the ocean for society's continuing needs. The decade beginning 2021-2030 has additionally been declared as the UN Decade of Ocean Science for Sustainable Development. This program aims to maximize the benefits of ocean science to the management, conservation, and sustainable development of the marine environment by facilitating communication and cooperation at the science-policy interface. A central principle of the program is the conservation of species and ecosystem components of biodiversity. However, a significant omission from the draft version of the Decade of Ocean Science Implementation Plan is the acknowledgment of the importance of monitoring and maintaining genetic biodiversity within species. In this paper, we emphasize the importance of genetic diversity to adaptive capacity, evolutionary potential, community function, and resilience within populations, as well as highlighting some of the major threats to genetic diversity in the marine environment from direct human impacts and the effects of global climate change. We then highlight the significance of ocean genetic diversity to a diverse range of socioeconomic factors in the marine environment, including marine industries, welfare and leisure pursuits, coastal communities, and wider society. Genetic biodiversity in the ocean, and its monitoring and maintenance, is then discussed with respect to its integral role in the successful realization of the 2030 vision for the Decade of Ocean Science. Finally, we suggest how ocean genetic diversity might be better integrated into biodiversity management practices through the continued interaction between environmental managers and scientists, as well as through key leverage points in industry requirements for Blue Capital financing and social responsibility.

Reference genome and demographic history of the most endangered marine mammal, the vaquita.

The vaquita is the most critically endangered marine mammal, with fewer than 19 remaining in the wild. First described in 1958, the vaquita has been in rapid decline for more than 20 years resulting from inadvertent deaths due to the increasing use of large-mesh gillnets. To understand the evolutionary and demographic history of the vaquita, we used combined long-read sequencing and long-range scaffolding methods with long- and short-read RNA sequencing to generate a near error-free annotated reference genome assembly from cell lines derived from a female individual. The genome assembly consists of 99.92% of the assembled sequence contained in 21 nearly gapless chromosome-length autosome scaffolds and the X-chromosome scaffold, with a scaffold N50 of 115 Mb. Genome-wide heterozygosity is the lowest (0.01%) of any mammalian species analysed to date, but heterozygosity is evenly distributed across the chromosomes, consistent with long-term small population size at genetic equilibrium, rather than low diversity resulting from a recent population bottleneck or inbreeding. Historical demography of the vaquita indicates long-term population stability at less than 5,000 (Ne) for over 200,000 years. Together, these analyses indicate that the vaquita genome has had ample opportunity to purge highly deleterious alleles and potentially maintain diversity necessary for population health.

Acquisition of obligate mutualist symbionts during the larval stage is not beneficial for a coral host.

Theory suggests that the direct transmission of beneficial endosymbionts (mutualists) from parents to offspring (vertical transmission) in animal hosts is advantageous and evolutionarily stable, yet many host species instead acquire their symbionts from the environment (horizontal acquisition). An outstanding question in marine biology is why some scleractinian corals do not provision their eggs and larvae with the endosymbiotic dinoflagellates that are necessary for a juvenile's ultimate survival. We tested whether the acquisition of photosynthetic endosymbionts (family Symbiodiniaceae) during the planktonic larval stage was advantageous, as is widely assumed, in the ecologically important and threatened Caribbean reef-building coral Orbicella faveolata. Following larval acquisition, similar changes occurred in host energetic lipid use and gene expression regardless of whether their symbionts were photosynthesizing, suggesting the symbionts did not provide the energetic benefit characteristic of the mutualism in adults. Larvae that acquired photosymbionts isolated from conspecific adults on their natal reef exhibited a reduction in swimming, which may interfere with their ability to find suitable settlement substrate, and also a decrease in survival. Larvae exposed to two cultured algal species did not exhibit differences in survival, but decreased their swimming activity in response to one species. We conclude that acquiring photosymbionts during the larval stage confers no advantages and can in fact be disadvantageous to this coral host. The timing of symbiont acquisition appears to be a critical component of a host's life history strategy and overall reproductive fitness, and this timing itself appears to be under selective pressure.

Disentangling genetic structure for genetic monitoring of complex populations.

Genetic monitoring estimates temporal changes in population parameters from molecular marker information. Most populations are complex in structure and change through time by expanding or contracting their geographic range, becoming fragmented or coalescing, or increasing or decreasing density. Traditional approaches to genetic monitoring rely on quantifying temporal shifts of specific population metrics-heterozygosity, numbers of alleles, effective population size-or measures of geographic differentiation such as FST. However, the accuracy and precision of the results can be heavily influenced by the type of genetic marker used and how closely they adhere to analytical assumptions. Care must be taken to ensure that inferences reflect actual population processes rather than changing molecular techniques or incorrect assumptions of an underlying model of population structure. In many species of conservation concern, true population structure is unknown, or structure might shift over time. In these cases, metrics based on inappropriate assumptions of population structure may not provide quality information regarding the monitored population. Thus, we need an inference model that decouples the complex elements that define population structure from estimation of population parameters of interest and reveals, rather than assumes, fine details of population structure. Encompassing a broad range of possible population structures would enable comparable inferences across biological systems, even in the face of range expansion or contraction, fragmentation, or changes in density. Currently, the best candidate is the spatial Λ-Fleming-Viot (SLFV) model, a spatially explicit individually based coalescent model that allows independent inference of two of the most important elements of population structure: local population density and local dispersal. We support increased use of the SLFV model for genetic monitoring by highlighting its benefits over traditional approaches. We also discuss necessary future directions for model development to support large genomic datasets informing real-world management and conservation issues.

Comparative growth rates of cultured marine dinoflagellates in the genus Symbiodinium and the effects of temperature and light.

Many dinoflagellate microalgae of the genus Symbiodinium form successful symbioses with a large group of metazoans and selected protists. Yet knowledge of growth kinetics of these endosymbionts and their ecological and evolutionary implications is limited. We used a Bayesian biphasic generalized logistic model to estimate key parameters of the growth of five strains of cultured Symbiodinium, S. microadriaticum (cp-type A194; strain 04-503), S. microadriaticum (cp-type A194; strain CassKB8), S. minutum (cp-type B184; strain Mf 1.05b.01.SCI.01), S. psygmophilum (cp-type B224; strain Mf 11.05b.01) and S. trenchii (cp-type D206; strain Mf 2.2b), grown in four different combinations of temperature and light. Growth kinetics varied among Symbiodinium strains and across treatments. Biphasic growth was especially evident for S. minutum and S. psygmophilum across all treatments. Monophasic growth was more common when final asymptotic densities were relatively low (~ 200 million cells ml-1). All species tended to grow faster and / or reached a higher asymptote at 26°C than at 18°C. The fastest growth was exhibited by S. minutum, with an approximate four-fold increase in estimated cell density after 60 days. The strongest effect of light was seen in S. trenchii, in which increasing light levels resulted in a decrease in initial growth rate, and an increase in asymptotic density, time when growth rate was at its maximum, final growth rate, and maximum growth rate. Results suggest that Symbiodinium species have different photokinetic and thermal optima, which may affect their growth-related nutritional physiology and allow them to modify their response to environmental changes.

Meta-mass shift chemical profiling of metabolomes from coral reefs.

Untargeted metabolomics of environmental samples routinely detects thousands of small molecules, the vast majority of which cannot be identified. Meta-mass shift chemical (MeMSChem) profiling was developed to identify mass differences between related molecules using molecular networks. This approach illuminates metabolome-wide relationships between molecules and the putative chemical groups that differentiate them (e.g., H2, CH2, COCH2). MeMSChem profiling was used to analyze a publicly available metabolomic dataset of coral, algal, and fungal mat holobionts (i.e., the host and its associated microbes and viruses) sampled from some of Earth's most remote and pristine coral reefs. Each type of holobiont had distinct mass shift profiles, even when the analysis was restricted to molecules found in all samples. This result suggests that holobionts modify the same molecules in different ways and offers insights into the generation of molecular diversity. Three genera of stony corals had distinct patterns of molecular relatedness despite their high degree of taxonomic relatedness. MeMSChem profiles also partially differentiated between individuals, suggesting that every coral reef holobiont is a potential source of novel chemical diversity.

stratag: An r package for manipulating, summarizing and analysing population genetic data.

We introduce the r package stratag as a user-friendly population genetics toolkit. stratag provides easy access to a suite of standard genetic summaries as well as the ability to rapidly manipulate stratified genetic data for custom analyses. Tests of population subdivision with most common measures of population subdivision (e.g., FST , GST , ΦST , Χ2 ) can be conducted within a single function. The package also provides wrapper functions that allow users to configure and run popular external programs such as genepop, structure, and fastsimcoal from within r, and smoothly interface with popular r packages adegenet and pegas. stratag is intended to be an open-source dynamic package that will grow with future needs and user input.

skelesim: an extensible, general framework for population genetic simulation in R.

Simulations are a key tool in molecular ecology for inference and forecasting, as well as for evaluating new methods. Due to growing computational power and a diversity of software with different capabilities, simulations are becoming increasingly powerful and useful. However, the widespread use of simulations by geneticists and ecologists is hindered by difficulties in understanding these softwares' complex capabilities, composing code and input files, a daunting bioinformatics barrier and a steep conceptual learning curve. skelesim (an R package) guides users in choosing appropriate simulations, setting parameters, calculating genetic summary statistics and organizing data output, in a reproducible pipeline within the R environment. skelesim is designed to be an extensible framework that can 'wrap' around any simulation software (inside or outside the R environment) and be extended to calculate and graph any genetic summary statistics. Currently, skelesim implements coalescent and forward-time models available in the fastsimcoal2 and rmetasim simulation engines to produce null distributions for multiple population genetic statistics and marker types, under a variety of demographic conditions. skelesim is intended to make simulations easier while still allowing full model complexity to ensure that simulations play a fundamental role in molecular ecology investigations. skelesim can also serve as a teaching tool: demonstrating the outcomes of stochastic population genetic processes; teaching general concepts of simulations; and providing an introduction to the R environment with a user-friendly graphical user interface (using shiny).

Geographic and temporal dynamics of a global radiation and diversification in the killer whale.

Global climate change during the Late Pleistocene periodically encroached and then released habitat during the glacial cycles, causing range expansions and contractions in some species. These dynamics have played a major role in geographic radiations, diversification and speciation. We investigate these dynamics in the most widely distributed of marine mammals, the killer whale (Orcinus orca), using a global data set of over 450 samples. This marine top predator inhabits coastal and pelagic ecosystems ranging from the ice edge to the tropics, often exhibiting ecological, behavioural and morphological variation suggestive of local adaptation accompanied by reproductive isolation. Results suggest a rapid global radiation occurred over the last 350 000 years. Based on habitat models, we estimated there was only a 15% global contraction of core suitable habitat during the last glacial maximum, and the resources appeared to sustain a constant global effective female population size throughout the Late Pleistocene. Reconstruction of the ancestral phylogeography highlighted the high mobility of this species, identifying 22 strongly supported long-range dispersal events including interoceanic and interhemispheric movement. Despite this propensity for geographic dispersal, the increased sampling of this study uncovered very few potential examples of ancestral dispersal among ecotypes. Concordance of nuclear and mitochondrial data further confirms genetic cohesiveness, with little or no current gene flow among sympatric ecotypes. Taken as a whole, our data suggest that the glacial cycles influenced local populations in different ways, with no clear global pattern, but with secondary contact among lineages following long-range dispersal as a potential mechanism driving ecological diversification.

Taxonomic and environmental variation of metabolite profiles in marine dinoflagellates of the genus symbiodinium.

Microorganisms in terrestrial and marine ecosystems are essential to environmental sustainability. In the marine environment, invertebrates often depend on metabolic cooperation with their endosymbionts. Coral reefs, one of the most important marine ecosystems, are based on the symbiosis between a broad diversity of dinoflagellates of the genus Symbiodinium and a wide phyletic diversity of hosts (i.e., cnidarian, molluscan, poriferan). This diversity is reflected in the ecology and physiology of the symbionts, yet the underlying biochemical mechanisms are still poorly understood. We examined metabolite profiles of four cultured species of Symbiodinium known to form viable symbioses with reef-building corals, S. microadriaticum (cp-type A194), S. minutum (cp-type B184), S. psygmophilum (cp-type B224) and S. trenchii (cp-type D206). Metabolite profiles were shown to differ among Symbiodinium species and were found to be affected by their physiological response to growth in different temperatures and light regimes. A combined Random Forests and Bayesian analysis revealed that the four Symbiodinium species examined primarily differed in their production of sterols and sugars, including a C29 stanol and the two sterols C28Δ5 and C28Δ5,22, as well as differences in metabolite abundances of a hexose and inositol. Inositol levels were also strongly affected by changes in temperature across all Symbiodinium species. Our results offer a detailed view of the metabolite profile characteristic of marine symbiotic dinoflagellates of the genus Symbiodinium, and identify patterns of metabolites related to several growth conditions.

Nuclear and mitochondrial patterns of population structure in North Pacific false killer whales (Pseudorca crassidens).

False killer whales (Pseudorca crassidens) are large delphinids typically found in deep water far offshore. However, in the Hawaiian Archipelago, there are 2 resident island-associated populations of false killer whales, one in the waters around the main Hawaiian Islands (MHI) and one in the waters around the Northwestern Hawaiian Islands (NWHI). We use mitochondrial DNA (mtDNA) control region sequences and genotypes from 16 nuclear DNA (nucDNA) microsatellite loci from 206 individuals to examine levels of differentiation among the 2 island-associated populations and offshore animals from the central and eastern North Pacific. Both mtDNA and nucDNA exhibit highly significant differentiation between populations, confirming limited gene flow in both sexes. The mtDNA haplotypes exhibit a strong pattern of phylogeographic concordance, with island-associated populations sharing 3 closely related haplotypes not found elsewhere in the Pacific. However, nucDNA data suggest that NWHI animals are at least as differentiated from MHI animals as they are from offshore animals. The patterns of differentiation revealed by the 2 marker types suggest that the island-associated false killer whale populations likely share a common colonization history, but have limited contemporary gene flow.

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.

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.