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.

Modeling population effects of the Deepwater Horizon oil spill on a long-lived species.

The 2010 Deepwater Horizon (DWH) oil spill exposed common bottlenose dolphins (Tursiops truncatus) in Barataria Bay, Louisiana to heavy oiling that caused increased mortality and chronic disease and impaired reproduction in surviving dolphins. We conducted photographic surveys and veterinary assessments in the decade following the spill. We assigned a prognostic score (good, fair, guarded, poor, or grave) for each dolphin to provide a single integrated indicator of overall health, and we examined temporal trends in prognostic scores. We used expert elicitation to quantify the implications of trends for the proportion of the dolphins that would recover within their lifetime. We integrated expert elicitation, along with other new information, in a population dynamics model to predict the effects of observed health trends on demography. We compared the resulting population trajectory with that predicted under baseline (no spill) conditions. Disease conditions persisted and have recently worsened in dolphins that were presumably exposed to DWH oil: 78% of those assessed in 2018 had a guarded, poor, or grave prognosis. Dolphins born after the spill were in better health. We estimated that the population declined by 45% (95% CI 14-74) relative to baseline and will take 35 years (95% CI 18-67) to recover to 95% of baseline numbers. The sum of annual differences between baseline and injured population sizes (i.e., the lost cetacean years) was 30,993 (95% CI 6607-94,148). The population is currently at a minimum point in its recovery trajectory and is vulnerable to emerging threats, including planned ecosystem restoration efforts that are likely to be detrimental to the dolphins' survival. Our modeling framework demonstrates an approach for integrating different sources and types of data, highlights the utility of expert elicitation for indeterminable input parameters, and emphasizes the importance of considering and monitoring long-term health of long-lived species subject to environmental disasters. Article impact statement: Oil spills can have long-term consequences for the health of long-lived species; thus, effective restoration and monitoring are needed.

El derrame de petróleo Deepwater Horizon (DWH) en 2010 expuso gravemente a este hidrocarburo a los delfines (Tursiops truncatus) de la Bahía Barataria, Luisiana, causando un incremento en la mortalidad y en las enfermedades crónicas, y deteriorando la reproducción de los delfines sobrevivientes. Realizamos censos fotográficos y evaluaciones veterinarias durante la década posterior al derrame. Asignamos un puntaje pronóstico (bueno, favorable, moderado, malo, o grave) a cada delfín para proporcionar un indicador integrado único de la salud en general. También examinamos las tendencias temporales de estos puntajes. Usamos información de expertos para cuantificar las implicaciones de las tendencias para la proporción de delfines que se recuperaría dentro de su periodo de vida. Integramos esta información, junto con información nueva, a un modelo de dinámica poblacional para predecir los efectos sobre la demografía de las tendencias observadas en la salud. Comparamos la trayectoria poblacional resultante con aquella pronosticada bajo condiciones de línea base (sin derrame). Las condiciones de enfermedad persistieron y recientemente han empeorado en los delfines que supuestamente estuvieron expuestos al petróleo de DWH: 78% de aquellos evaluados en 2018 tuvieron un pronóstico moderado, malo o grave. Los delfines que nacieron después del derrame contaron con mejor salud. Estimamos que la población declinó en un 45% (95% CI 14-74) relativo a la línea base y tardará 35 años (95% CI 18-67) en recuperar el 95% de los números de línea base. La suma de las diferencias anuales entre el tamaño poblacional de línea base y el dañado (es decir, los años cetáceos perdidos) fue de 30,993 (95% CI 6,607-94,148). La población actualmente está en un punto mínimo de su trayectoria de recuperación y es vulnerable a las amenazas emergentes, incluyendo los esfuerzos de restauración ambiental planeada que probablemente sean nocivos para la supervivencia de los delfines. Nuestro marco de modelado demuestra una estrategia para la integración de diferentes fuentes y tipos de datos, resalta la utilidad de la información de expertos para los parámetros de aportación indeterminable, y enfatiza la importancia de la consideración y el monitoreo de la salud a largo plazo de las especies longevas sujetas a los desastres ambientales. Modelado de los Efectos Poblacionales del Derrame de Petróleo Deepwater Horizon sobre Especies Longevas.

Ecological divergence and speciation in common bottlenose dolphins in the western South Atlantic.

Coastal and offshore ecotypes of common bottlenose dolphins have been recognized in the western South Atlantic, and it is possible that trophic niche divergence associated with social interactions is leading them to genetic and phenotypic differentiation. The significant morphological differentiation observed between these ecotypes suggests they represent two different subspecies. However, there is still a need to investigate whether there is congruence between morphological and genetic data to rule out the possibility of ecophenotypic variation accompanied by gene flow. Mitochondrial DNA (mtDNA) control region sequence data and 10 microsatellite loci collected from stranded and biopsied dolphins sampled in coastal and offshore waters of Brazil as well as 106 skulls for morphological analyses were used to determine whether the morphological differentiation was supported by genetic differentiation. There was congruence among the data sets, reinforcing the presence of two distinct ecotypes. The divergence may be relatively recent, however, given the moderate values of mtDNA nucleotide divergence (dA = 0.008), presence of one shared mtDNA haplotype and possibly low levels of gene flow (around 1% of migrants per generation). Results suggest the ecotypes may be in the process of speciation and reinforce they are best described as two different subspecies until the degree of nuclear genetic divergence is thoroughly evaluated: Tursiops truncatus gephyreus (coastal ecotype) and T. t. truncatus (offshore ecotype). The endemic distribution of T. t. gephyreus in the western South Atlantic and number of anthropogenic threats in the area reinforces the importance of protecting this ecotype and its habitat.

Oceanographic barriers, divergence, and admixture: Phylogeography and taxonomy of two putative subspecies of short-finned pilot whale.

Genomic phylogeography plays an important role in describing evolutionary processes and their geographic, ecological, or cultural drivers. These drivers are often poorly understood in marine environments, which have fewer obvious barriers to mixing than terrestrial environments. Taxonomic uncertainty of some taxa (e.g., cetaceans), due to the difficulty in obtaining morphological data, can hamper our understanding of these processes. One such taxon, the short-finned pilot whale, is recognized as a single global species but includes at least two distinct morphological forms described from stranding and drive hunting in Japan, the "Naisa" and "Shiho" forms. Using samples (n = 735) collected throughout their global range, we examine phylogeographic patterns of divergence by comparing mitogenomes and nuclear SNP loci. Our results suggest three types within the species: an Atlantic Ocean type, a western/central Pacific and Indian Ocean (Naisa) type, and an eastern Pacific Ocean and northern Japan (Shiho) type. mtDNA control region differentiation indicates these three types form two subspecies, separated by the East Pacific Barrier: Shiho short-finned pilot whale, in the eastern Pacific Ocean and northern Japan, and Naisa short-finned pilot whale, throughout the remainder of the species' distribution. Our data further indicate two diverging populations within the Naisa subspecies, in the Atlantic Ocean and western/central Pacific and Indian Oceans, separated by the Benguela Barrier off South Africa. This study reveals a process of divergence and speciation within a globally-distributed, mobile marine predator, and indicates the importance of the East Pacific Barrier to this evolutionary process.

Demography or selection on linked cultural traits or genes? Investigating the driver of low mtDNA diversity in the sperm whale using complementary mitochondrial and nuclear genome analyses.

Mitochondrial DNA has been heavily utilized in phylogeography studies for several decades. However, underlying patterns of demography and phylogeography may be misrepresented due to coalescence stochasticity, selection, variation in mutation rates and cultural hitchhiking (linkage of genetic variation to culturally-transmitted traits affecting fitness). Cultural hitchhiking has been suggested as an explanation for low genetic diversity in species with strong social structures, counteracting even high mobility, abundance and limited barriers to dispersal. One such species is the sperm whale, which shows very limited phylogeographic structure and low mtDNA diversity despite a worldwide distribution and large population. Here, we use analyses of 175 globally distributed mitogenomes and three nuclear genomes to evaluate hypotheses of a population bottleneck/expansion vs. a selective sweep due to cultural hitchhiking or selection on mtDNA as the mechanism contributing to low worldwide mitochondrial diversity in sperm whales. In contrast to mtDNA control region (CR) data, mitogenome haplotypes are largely ocean-specific, with only one of 80 shared between the Atlantic and Pacific. Demographic analyses of nuclear genomes suggest low mtDNA diversity is consistent with a global reduction in population size that ended approximately 125,000 years ago, correlated with the Eemian interglacial. Phylogeographic analysis suggests that extant sperm whales descend from maternal lineages endemic to the Pacific during the period of reduced abundance and have subsequently colonized the Atlantic several times. Results highlight the apparent impact of past climate change, and suggest selection and hitchhiking are not the sole processes responsible for low mtDNA diversity in this highly social species.

Looking into a whale's heart: investigating a genetic basis for cardiomyopathy in a non-model species.

Understanding the pathogenesis of complex diseases can benefit from multi-species comparative studies. Yet these studies rarely include natural populations of non-model species. Here, we focused on the cause of a heart muscle disease, cardiomyopathy (CM), affecting multiple mammalian species including humans, cats, dogs, and certain species of whales. Mutations in genes coding for sarcomeric proteins have been identified as a leading cause for CM in humans, and some were also revealed to be responsible for CM in cats. We investigated whether similar mutations could be detected in the deep-diving pygmy sperm whale (Kogia breviceps), which is one of two cetacean species known to display CM. We sequenced portions of two candidate genes (MYH7: 3153 bp and MYBPC3: 3019 bp) in 55 whales including affected and unaffected individuals. Mutation screening revealed six nonsynonymous substitutions that were predicted to have an effect on protein function. However, the etiology of CM is likely complex and probably multi-factorial as three of these mutations were observed in unaffected individuals from our control group. This incomplete penetrance could be partly age-related and could also be due to the influence of environmental factors on the development of CM, as seen in humans.

From genome-wide to candidate gene: an investigation of variation at the major histocompatibility complex in common bottlenose dolphins exposed to harmful algal blooms.

The role the major histocompatibility complex (MHC) plays in response to exposure to environmental toxins is relatively poorly understood, particularly in comparison to its well-described role in pathogen immunity. We investigated associations between MHC diversity and resistance to brevetoxins in common bottlenose dolphins (Tursiops truncatus). A previous genome-wide association study investigating an apparent difference in harmful algal bloom (HAB) resistance among dolphin populations in the Gulf of Mexico identified genetic variation associated with survival in close genomic proximity to multiple MHC class II loci. Here, we characterized genetic variation at DQA, DQB, DRA, and DRB loci in dolphins from central-west Florida and the Florida Panhandle, including dolphins that died during HABs and dolphins presumed to have survived HAB exposure. We found that DRB and DQB exhibited patterns of genetic differentiation among geographic regions that differed from neutral microsatellite loci. In addition, genetic differentiation at DRB across multiple pairwise comparisons of live and dead dolphins was greater than differentiation observed at neutral loci. Our findings at these MHC loci did not approach the strength of association with survival previously described for a nearby genetic variant. However, the results provide evidence that selective pressures at the MHC vary among dolphin populations that differ in the frequency of HAB exposure and that the overall composition of DRB variants differs between dolphin survivors and non-survivors of HABs. These results may suggest a potential role of MHC diversity in variable survival of bottlenose dolphins exposed to HABs.

Reproductive outcome and survival of common bottlenose dolphins sampled in Barataria Bay, Louisiana, USA, following the Deepwater Horizon oil spill.

Common bottlenose dolphins (Tursiops truncatus) inhabit bays, sounds and estuaries across the Gulf of Mexico. Following the Deepwater Horizon oil spill, studies were initiated to assess potential effects on these ecologically important apex predators. A previous study reported disease conditions, including lung disease and impaired stress response, for 32 dolphins that were temporarily captured and given health assessments in Barataria Bay, Louisiana, USA. Ten of the sampled dolphins were determined to be pregnant, with expected due dates the following spring or summer. Here, we report findings after 47 months of follow-up monitoring of those sampled dolphins. Only 20% (95% CI: 2.50-55.6%) of the pregnant dolphins produced viable calves, as compared with a previously reported pregnancy success rate of 83% in a reference population. Fifty-seven per cent of pregnant females that did not successfully produce a calf had been previously diagnosed with moderate-severe lung disease. In addition, the estimated annual survival rate of the sampled cohort was low (86.8%, 95% CI: 80.0-92.7%) as compared with survival rates of 95.1% and 96.2% from two other previously studied bottlenose dolphin populations. Our findings confirm low reproductive success and high mortality in dolphins from a heavily oiled estuary when compared with other populations. Follow-up studies are needed to better understand the potential recovery of dolphins in Barataria Bay and, by extension, other Gulf coastal regions impacted by the spill.

Genomewide investigation of adaptation to harmful algal blooms in common bottlenose dolphins (Tursiops truncatus).

Harmful algal blooms (HABs), which can be lethal in marine species and cause illness in humans, are increasing worldwide. In the Gulf of Mexico, HABs of Karenia brevis produce neurotoxic brevetoxins that cause large-scale marine mortality events. The long history of such blooms, combined with the potentially severe effects of exposure, may have produced a strong selective pressure for evolved resistance. Advances in next-generation sequencing, in particular genotyping-by-sequencing, greatly enable the genomic study of such adaptation in natural populations. We used restriction site-associated DNA (RAD) sequencing to investigate brevetoxicosis resistance in common bottlenose dolphins (Tursiops truncatus). To improve our understanding of the epidemiology and aetiology of brevetoxicosis and the potential for evolved resistance in an upper trophic level predator, we sequenced pools of genomic DNA from dolphins sampled from both coastal and estuarine populations in Florida and during multiple HAB-associated mortality events. We sequenced 129 594 RAD loci and analysed 7431 single nucleotide polymorphisms (SNPs). The allele frequencies of many of these polymorphic loci differed significantly between live and dead dolphins. Some loci associated with survival showed patterns suggesting a common genetic-based mechanism of resistance to brevetoxins in bottlenose dolphins along the Gulf coast of Florida, but others suggested regionally specific mechanisms of resistance or reflected differences among HABs. We identified candidate genes that may be the evolutionary target for brevetoxin resistance by searching the dolphin genome for genes adjacent to survival-associated SNPs.

Hierarchical population structure and habitat differences in a highly mobile marine species: the Atlantic spotted dolphin.

Recent molecular studies have shown that highly mobile species with continuous distributions can exhibit fine-scale population structure. In this context, we assessed genetic structure within a marine species with high dispersal potential, the Atlantic spotted dolphin (Stenella frontalis). Using 19 microsatellite loci and mitochondrial control region sequences, population structure was investigated in the western North Atlantic, the Gulf of Mexico and the Azores Islands. Analyses of the microsatellite data identified four distinct genetic clusters, which were supported by the control region sequences. The highest level of divergence was seen between two clusters corresponding to previously described morphotypes that inhabit oceanic and shelf waters. The combined morphological and genetic evidence suggests these two lineages are on distinct evolutionary trajectories and could be considered distinct subspecies despite their parapatry. Further analysis of the continental shelf cluster resulted in three groups: animals inhabiting shelf waters in the western North Atlantic, the eastern Gulf of Mexico and the western Gulf of Mexico. Analyses of environmental data indicate the four genetic clusters inhabit distinct habitats in terms of depth and sea surface temperature. Contemporary dispersal rate estimates suggest all of these populations should be considered as distinct management units. Conversely, no significant genetic differentiation was observed between S. frontalis from offshore waters of the western North Atlantic and the Azores, which are separated by approximately 4500 km. Overall, the hierarchical structure observed within the Atlantic spotted dolphin shows that the biogeography of the species is complex because it is not shaped solely by geographic distance.

Review of historical unusual mortality events (UMEs) in the Gulf of Mexico (1990-2009): providing context for the multi-year northern Gulf of Mexico cetacean UME declared in 2010.

An unusual mortality event (UME) was declared for cetaceans in the northern Gulf of Mexico (GoM) for Franklin County, Florida, west through Louisiana, USA, beginning in February 2010 and was ongoing as of September 2014. The 'Deepwater Horizon' (DWH) oil spill began on 20 April 2010 in the GoM, raising questions regarding the potential role of the oil spill in the UME. The present study reviews cetacean mortality events that occurred in the GoM prior to 2010 (n = 11), including causes, durations, and some specific test results, to provide a historical context for the current event. The average duration of GoM cetacean UMEs prior to 2010 was 6 mo, and the longest was 17 mo (2005-2006). The highest number of cetacean mortalities recorded during a previous GoM event was 344 (in 1990). In most previous events, dolphin morbillivirus or brevetoxicosis was confirmed or suspected as a causal factor. In contrast, the current northern GoM UME has lasted more than 48 mo and has had more than 1000 reported mortalities within the currently defined spatial and temporal boundaries of the event. Initial results from the current UME do not support either morbillivirus or brevetoxin as primary causes of this event. This review is the first summary of cetacean UMEs in the GoM and provides evidence that the most common causes of previous UMEs are unlikely to be associated with the current UME.

Ancient dolphin genomes reveal rapid repeated adaptation to coastal waters.

Parallel evolution provides strong evidence of adaptation by natural selection due to local environmental variation. Yet, the chronology, and mode of the process of parallel evolution remains debated. Here, we harness the temporal resolution of paleogenomics to address these long-standing questions, by comparing genomes originating from the mid-Holocene (8610-5626 years before present, BP) to contemporary pairs of coastal-pelagic ecotypes of bottlenose dolphin. We find that the affinity of ancient samples to coastal populations increases as the age of the samples decreases. We assess the youngest genome (5626 years BP) at sites previously inferred to be under parallel selection to coastal habitats and find it contained coastal-associated genotypes. Thus, coastal-associated variants rose to detectable frequencies close to the emergence of coastal habitat. Admixture graph analyses reveal a reticulate evolutionary history between pelagic and coastal populations, sharing standing genetic variation that facilitated rapid adaptation to newly emerged coastal habitats.

Vulnerability to climate change of United States marine mammal stocks in the western North Atlantic, Gulf of Mexico, and Caribbean.

Climate change and climate variability are affecting marine mammal species and these impacts are projected to continue in the coming decades. Vulnerability assessments provide a framework for evaluating climate impacts over a broad range of species using currently available information. We conducted a trait-based climate vulnerability assessment using expert elicitation for 108 marine mammal stocks and stock groups in the western North Atlantic, Gulf of Mexico, and Caribbean Sea. Our approach combined the exposure (projected change in environmental conditions) and sensitivity (ability to tolerate and adapt to changing conditions) of marine mammal stocks to estimate vulnerability to climate change, and categorize stocks with a vulnerability index. The climate vulnerability score was very high for 44% (n = 47) of these stocks, high for 29% (n = 31), moderate for 20% (n = 22), and low for 7% (n = 8). The majority of stocks (n = 78; 72%) scored very high exposure, whereas 24% (n = 26) scored high, and 4% (n = 4) scored moderate. The sensitivity score was very high for 33% (n = 36) of these stocks, high for 18% (n = 19), moderate for 34% (n = 37), and low for 15% (n = 16). Vulnerability results were summarized for stocks in five taxonomic groups: pinnipeds (n = 4; 25% high, 75% moderate), mysticetes (n = 7; 29% very high, 57% high, 14% moderate), ziphiids (n = 8; 13% very high, 50% high, 38% moderate), delphinids (n = 84; 52% very high, 23% high, 15% moderate, 10% low), and other odontocetes (n = 5; 60% high, 40% moderate). Factors including temperature, ocean pH, and dissolved oxygen were the primary drivers of high climate exposure, with effects mediated through prey and habitat parameters. We quantified sources of uncertainty by bootstrapping vulnerability scores, conducting leave-one-out analyses of individual attributes and individual scorers, and through scoring data quality for each attribute. These results provide information for researchers, managers, and the public on marine mammal responses to climate change to enhance the development of more effective marine mammal management, restoration, and conservation activities that address current and future environmental variation and biological responses due to climate change.

The occurrence of mtDNA heteroplasmy in multiple cetacean species.

In population genetics and phylogenetic studies, mitochondrial DNA (mtDNA) is commonly used for examining differences both between and within groups of individuals. For these studies, correct interpretation of every nucleotide position is crucial but can be complicated by the presence of ambiguous bases resulting from heteroplasmy. Particularly for non-model taxa, the presence of heteroplasmy in mtDNA is rarely reported, therefore, it is unclear how commonly it occurs and how it can affect phylogenetic relationships among taxa and the overall understanding of evolutionary processes. We examined the occurrence of both site and length heteroplasmy within the mtDNA of ten marine mammal species, for most of which mtDNA heteroplasmy has never been reported. After sequencing a portion of the mtDNA control region for 5,062 individuals, we found heteroplasmy in at least 2% of individuals from seven species, including Stenella frontalis where 58.9% were heteroplasmic. We verified the presence of true heteroplasmy, ruling out artifacts from amplification and sequencing methods and the presence of nuclear copies of mitochondrial genes. We found no evidence that mtDNA heteroplasmy influenced phylogenetic relationships, however, its occurrence does have the potential to increase the genetic diversity for all species in which it is found. This study stresses the importance of both detecting and reporting the occurrence of heteroplasmy in wild populations in order to enhance the knowledge of both the introduction and the persistence of mutant mitochondrial haplotypes in the evolutionary process.

No leading-edge effect in North Atlantic harbor porpoises: Evolutionary and conservation implications.

Understanding species responses to past environmental changes can help forecast how they will cope with ongoing climate changes. Harbor porpoises are widely distributed in the North Atlantic and were deeply impacted by the Pleistocene changes with the split of three subspecies. Despite major impacts of fisheries on natural populations, little is known about population connectivity and dispersal, how they reacted to the Pleistocene changes, and how they will evolve in the future. Here, we used phylogenetics, population genetics, and predictive habitat modeling to investigate population structure and phylogeographic history of the North Atlantic porpoises. A total of 925 porpoises were characterized at 10 microsatellite loci and one quarter of the mitogenome (mtDNA). A highly divergent mtDNA lineage was uncovered in one porpoise off Western Greenland, suggesting that a cryptic group may occur and could belong to a recently discovered mesopelagic ecotype off Greenland. Aside from it and the southern subspecies, spatial genetic variation showed that porpoises from both sides of the North Atlantic form a continuous system belonging to the same subspecies (Phocoena phocoena phocoena). Yet, we identified important departures from random mating and restricted dispersal forming a highly significant isolation by distance (IBD) at both mtDNA and nuclear markers. A ten times stronger IBD at mtDNA compared with nuclear loci supported previous evidence of female philopatry. Together with the lack of spatial trends in genetic diversity, this IBD suggests that migration-drift equilibrium has been reached, erasing any genetic signal of a leading-edge effect that accompanied the predicted recolonization of the northern habitats freed from Pleistocene ice. These results illuminate the processes shaping porpoise population structure and provide a framework for designing conservation strategies and forecasting future population evolution.

Selection on ancestral genetic variation fuels repeated ecotype formation in bottlenose dolphins.

Studying repeated adaptation can provide insights into the mechanisms allowing species to adapt to novel environments. Here, we investigate repeated evolution driven by habitat specialization in the common bottlenose dolphin. Parapatric pelagic and coastal ecotypes of common bottlenose dolphins have repeatedly formed across the oceans. Analyzing whole genomes of 57 individuals, we find that ecotype evolution involved a complex reticulated evolutionary history. We find parallel linked selection acted upon ancient alleles in geographically distant coastal populations, which were present as standing genetic variation in the pelagic populations. Candidate loci evolving under parallel linked selection were found in ancient tracts, suggesting recurrent bouts of selection through time. Therefore, despite the constraints of small effective population size and long generation time on the efficacy of selection, repeated adaptation in long-lived social species can be driven by a combination of ecological opportunities and selection acting on ancestral standing genetic variation.

Testing mitochondrial sequences and anonymous nuclear markers for phylogeny reconstruction in a rapidly radiating group: molecular systematics of the Delphininae (Cetacea: Odontoceti: Delphinidae).

BACKGROUND: Many molecular phylogenetic analyses rely on DNA sequence data obtained from single or multiple loci, particularly mitochondrial DNA loci. However, phylogenies for taxa that have undergone recent, rapid radiation events often remain unresolved. Alternative methodologies for discerning evolutionary relationships under these conditions are desirable. The dolphin subfamily Delphininae is a group that has likely resulted from a recent and rapid radiation. Despite several efforts, the evolutionary relationships among the species in the subfamily remain unclear. RESULTS: Here, we compare a phylogeny estimated using mitochondrial DNA (mtDNA) control region sequences to a multi-locus phylogeny inferred from 418 polymorphic genomic markers obtained from amplified fragment length polymorphism (AFLP) analysis. The two sets of phylogenies are largely incongruent, primarily because the mtDNA tree provides very poor resolving power; very few species' nodes in the tree are supported by bootstrap resampling. The AFLP phylogeny is considerably better resolved and more congruent with relationships inferred from morphological data. Both phylogenies support paraphyly for the genera Stenella and Tursiops. The AFLP data indicate a close relationship between the two spotted dolphin species and recent ancestry between Stenella clymene and S. longirostris. The placement of the Lagenodelphis hosei lineage is ambiguous: phenetic analysis of the AFLP data is consistent with morphological expectations but the phylogenetic analysis is not. CONCLUSION: For closely related, recently diverged taxa, a multi-locus genome-wide survey is likely the most comprehensive approach currently available for phylogenetic inference.

Examining the relationships between blubber steroid hormones and persistent organic pollutants in common bottlenose dolphins.

Odontocete cetaceans bioaccumulate high concentrations of endocrine disrupting persistent organic pollutants (POPs), including dichlorodiphenyltrichloroethane (DDT), dichlorodiphenyltrichloroethylene (DDE), and dichlorodiphenyldichloroethane (DDD) - collectively DDTs - but few studies have explored DDTs-mediated endocrine disruption in cetaceans. Herein, we use remotely collected blubber biopsies from common bottlenose dolphins (Tursiops truncatus) inhabiting a site with high localized DDTs contamination to study the relationships between DDTs exposure and steroid hormone homeostasis in cetaceans. We quantified blubber steroid hormone concentrations by liquid chromatography-tandem mass spectrometry and blubber POP concentrations by gas chromatography-mass spectrometry. We detected six steroid hormones in blubber, including progesterone (P4), 17-hydroxyprogesterone (17OHP4), androstenedione (AE), testosterone (T), cortisol (F), and cortisone (E). Sampled dolphins (n = 62) exhibited exposure to DDT, DDE, DDD, chlordanes (CHLDs), mirex, dieldrin, hexachlorobenzene, polychlorinated biphenyls (PCBs), and brominated diphenyl ethers (BDEs). Using principal components analysis (PCA), we determined that blubber DDTs primarily loaded to the first principal component (PC1) explaining 81.6% of the total variance in POP exposure, while the remaining POPs primarily loaded to the PC2 (10.4% of variance). PC1 scores were negatively correlated with blubber T in males and blubber F in females, suggesting that exposure to DDTs impacted androgen and corticosteroid homeostasis. These conclusions were further supported by observed negative correlations between T and o,p'-DDE, o,p'-DDD, and p,p'-DDD in males sampled in the fall, and between F and the six individual DDTs and ∑6DDTs in females. Overall, these results suggest that POP-mediated endocrine disruption may have occurred in this stock of dolphins, which could negatively impact their health and fitness. However, this study relied on uncontrolled incidental exposures, making it impossible to establish a causal relationship between DDTs exposure and endocrine effects. Importantly, this study demonstrates that remotely collected blubber biopsies are a useful matrix for studying endocrine disruption in marine mammals.

Skin Transcriptomes of common bottlenose dolphins (Tursiops truncatus) from the northern Gulf of Mexico and southeastern U.S. Atlantic coasts.

Common bottlenose dolphins serve as sentinels for the health of their coastal environments as they are susceptible to health impacts from anthropogenic inputs through both direct exposure and food web magnification. Remote biopsy samples have been widely used to reveal contaminant burdens in free-ranging bottlenose dolphins, but do not address the health consequences of this exposure. To gain insight into whether remote biopsies can also identify health impacts associated with contaminant burdens, we employed RNA sequencing (RNA-seq) to interrogate the transcriptomes of remote skin biopsies from 116 bottlenose dolphins from the northern Gulf of Mexico and southeastern U.S. Atlantic coasts. Gene expression was analyzed using principal component analysis, differential expression testing, and gene co-expression networks, and the results correlated to season, location, and contaminant burden. Season had a significant impact, with over 60% of genes differentially expressed between spring/summer and winter months. Geographic location exhibited lesser effects on the transcriptome, with 23.5% of genes differentially expressed between the northern Gulf of Mexico and the southeastern U.S. Atlantic locations. Despite a large overlap between the seasonal and geographical gene sets, the pathways altered in the observed gene expression profiles were somewhat distinct. Co-regulated gene modules and differential expression analysis both identified epidermal development and cellular architecture pathways to be expressed at lower levels in animals from the northern Gulf of Mexico. Although contaminant burdens measured were not significantly different between regions, some correlation with contaminant loads in individuals was observed among co-expressed gene modules, but these did not include classical detoxification pathways. Instead, this study identified other, possibly downstream pathways, including those involved in cellular architecture, immune response, and oxidative stress, that may prove to be contaminant responsive markers in bottlenose dolphin skin.

Partially observed epidemics in wildlife hosts: modelling an outbreak of dolphin morbillivirus in the northwestern Atlantic, June 2013-2014.

Morbilliviruses cause major mortality in marine mammals, but the dynamics of transmission and persistence are ill understood compared to terrestrial counterparts such as measles; this is especially true for epidemics in cetaceans. However, the recent outbreak of dolphin morbillivirus in the northwestern Atlantic Ocean can provide new insights into the epidemiology and spatio-temporal spread of this pathogen. To deal with uncertainties surrounding the ecology of this system (only stranded animals were observed), we develop a statistical framework that can extract key information about the underlying transmission process given only sparse data. Our self-exciting Poisson process model suggests that individuals are infectious for at most 24 days and can transfer infection up to two latitude degrees (220 km) within this time. In addition, the effective reproduction number is generally below one, but reaches 2.6 during a period of heightened stranding numbers near Virginia Beach, Virginia, in summer 2013. Network analysis suggests local movements dominate spatial spread, with seasonal migration facilitating wider dissemination along the coast. Finally, a low virus transmission rate or high levels of pre-existing immunity can explain the lack of viral spread into the Gulf of Mexico. More generally, our approach illustrates novel methodologies for analysing very indirectly observed epidemics.