Frederick I. Archer

Mitogenome Phylogenetics: The Impact of Using Single Regions and Partitioning Schemes on Topology, Substitution Rate and Divergence Time Estimation

The availability of mitochondrial genome sequences is growing as a result of recent technological advances in molecular biology. In phylogenetic analyses, the complete mitogenome is increasingly becoming the marker of choice, usually providing better phylogenetic resolution and precision relative to traditional markers such as cytochrome b (CYTB) and the control region (CR). In some cases, the differences in phylogenetic estimates between mitogenomic and single-gene markers have yielded incongruent conclusions. By comparing phylogenetic estimates made from different genes, we identified the most informative mitochondrial regions and evaluated the minimum amount of data necessary to reproduce the same results as the mitogenome. We compared results among individual genes and the mitogenome for recently published complete mitogenome datasets of selected delphinids (Delphinidae) and killer whales (genus Orcinus). Using Bayesian phylogenetic methods, we investigated differences in estimation of topologies, divergence dates, and clock-like behavior among genes for both datasets. Although the most informative regions were not the same for each taxonomic group (COX1, CYTB, ND3 and ATP6 for Orcinus, and ND1, COX1 and ND4 for Delphinidae), in both cases they were equivalent to less than a quarter of the complete mitogenome. This suggests that gene information content can vary among groups, but can be adequately represented by a portion of the complete sequence. Although our results indicate that complete mitogenomes provide the highest phylogenetic resolution and most precise date estimates, a minimum amount of data can be selected using our approach when the complete sequence is unavailable. Studies based on single genes can benefit from the addition of a few more mitochondrial markers, producing topologies and date estimates similar to those obtained using the entire mitogenome.

Applied Conservation Genetics and the Need for Quality Control and Reporting of Genetic Data Used in Fisheries and Wildlife Management

Genetic data are often critical for defining populations for management purposes (e.g., identifying geographic boundaries or diagnostic characters for genetically discrete subunits) but can be called into question by both scientific and legal review. This can result in reversed or delayed implementation of management actions. We discuss methods for data quality control and quality analysis and describe examples of steps applied to 2 of the most common types of genetic data, mitochondrial DNA sequences, and microsatellite genotypes. These steps can serve both as guides to conservation geneticists and as an initial protocol for managers to determine whether genetic data will hold up against legal and scientific challenges. In addition, we suggest types of data and quality measures that should be reported as supplementary materials to published reports. These supplementary data serve to reduce the occurrence of legal and conservation controversies and improve reproducibility over time in population genetics studies where genetic monitoring is likely to play an increasing role.

Significant deviations from Hardy-Weinberg equilibrium caused by low levels of microsatellite genotyping errors

Microsatellite genotyping from samples with varying quality can result in an uneven distribution of errors. Previous studies reporting error rates have focused on estimating the effects of both randomly distributed and locus‐specific errors. Sample‐specific errors, however, can also significantly affect results in population studies despite a large sample size. From two studies including six microsatellite markers genotyped from 272 sperm whale DNA samples, and 33 microsatellites genotyped from 213 bowhead whales, we investigated the effects of sample‐ and locus‐specific errors on calculations of Hardy–Weinberg equilibrium. The results of a jackknife analysis in these two studies identified seven individuals that were highly influential on estimates of Hardy–Weinberg equilibrium for six different markers. In each case, the influential individual was homozygous for a rare allele. Our results demonstrate that Hardy–Weinberg P values are very sensitive to homozygosity in rare alleles for single individuals, and that > 50% of these cases involved genotype errors likely due to low sample quality. This raises the possibility that even small, normal levels of laboratory errors can result in an overestimate of the degree to which markers are out of Hardy–Weinberg equilibrium and hence overestimate population structure. To avoid such bias, we recommend routine identification of influential individuals and multiple replications of those samples.

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.

Sperm whale population structure in the eastern and central North Pacific inferred by the use of single-nucleotide polymorphisms, microsatellites and mitochondrial DNA.

We use mitochondrial DNA (mtDNA) (400 bp), six microsatellites and 36 single‐nucleotide polymorphisms (SNPs), 20 of which were linked, to investigate population structure of sperm whales (Physeter macrocephalus) in the eastern and central North Pacific. SNP markers, reproducible across technologies and laboratories, are ideal for long‐term studies of globally distributed species such as sperm whales, a species of conservation concern because of both historical and contemporary impacts. We estimate genetic differentiation among three strata in the temperate to tropical waters where females are found: California Current, Hawai`i and the eastern tropical Pacific. We then consider how males on sub‐Arctic foraging grounds assign to these strata. The California Current stratum was differentiated from both the other strata (P < 0.05) for mtDNA, microsatellites and SNPs, suggesting that the region supports a demographically independent population and providing the first indication that males may exhibit reproductive philopatry. Comparisons between the Hawai`i stratum and the eastern tropical Pacific stratum are not conclusive at this time. Comparisons with Alaska males were statistically significant, or nearly so, from all three strata and individuals showed mixed assignment to, and few exclusions from, the three potential source strata, suggesting widespread origin of males on sub‐Arctic feeding grounds. We show that SNPs have sufficient power to detect population structure even when genetic differentiation is low. There is a need for better analytical methods for SNPs, especially when linked SNPs are used, but SNPs appear to be a valuable marker for long‐term studies of globally dispersed and highly mobile species.

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.

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.

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.

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.

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.