http://orcid.org/0000-0003-1055-5299
Figures
Abstract
In threatened wildlife populations, it is important to determine whether observed low genetic diversity may be due to recent anthropogenic pressure or the consequence of historic events. Historical size of the Irrawaddy dolphin (Orcaella brevirostris) population inhabiting the Mekong River is unknown and there is significant concern for long-term survival of the remaining population as a result of low abundance, slow reproduction rate, high neonatal mortality, and continuing anthropogenic threats. We investigated population structure and reconstructed the demographic history based on 60 Irrawaddy dolphins samples collected between 2001 and 2009. The phylogenetic analysis indicated reciprocal monophyly of Mekong River Orcaella haplotypes with respect to haplotypes from other populations, suggesting long-standing isolation of the Mekong dolphin population from other Orcaella populations. We found that at least 85% of all individuals in the two main study areas: Kratie and Stung Treng, bore the same mitochondrial haplotype. Out of the 21 microsatellite loci tested, only ten were polymorphic and exhibited very low levels of genetic diversity. Both individual and frequency-based approaches suggest very low and non-significant genetic differentiation of the Mekong dolphin population. Evidence for recent bottlenecks was equivocal. Some results suggested a recent exponential decline in the Mekong dolphin population, with the current size being only 5.2% of the ancestral population. In order for the Mekong dolphin population to have any potential for long-term survival, it is imperative that management priorities focus on preventing any further population fragmentation or genetic loss, reducing or eliminating anthropogenic threats, and promoting connectivity between all subpopulations.
Citation: Krützen M, Beasley I, Ackermann CY, Lieckfeldt D, Ludwig A, Ryan GE, et al. (2018) Demographic collapse and low genetic diversity of the Irrawaddy dolphin population inhabiting the Mekong River. PLoS ONE 13(1): e0189200. https://doi.org/10.1371/journal.pone.0189200
Editor: Tzen-Yuh Chiang, National Cheng Kung University, TAIWAN
Received: February 10, 2017; Accepted: November 21, 2017; Published: January 3, 2018
Copyright: © 2018 Krützen et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The raw data (microsatellite genotypes and mtDNA haplotypes) were deposited on DRYAD (http://www.datadryad.org). doi:10.5061/dryad.45512.
Funding: The authors would like to thank Gordon Congdon (WWF-Greater Mekong Program) for his significant contribution to Mekong dolphin research and conservation efforts, and Verne Dove for her assistance with the stranding program and collection of samples used in this analysis. Many thanks to Martin Gilbert (WCS- Cambodia Program) who managed the Cambodian stranding program from 2005-2007. IB would like to thank the funding organisations that supported the Mekong Dolphin Conservation Project from 2001-2007; Ocean Park Conservation Foundation Hong Kong, the Mekong River Commission Environment program, Rufford Foundation Small Grants Program, Society for Marine Mammalogy Small Grants Program, IUCN Mekong Wetlands Biodiversity Project, James Cook University and the Wildlife Conservation Society Research Fellowship Program. DL and AL thank Petr Obrdlik (sample organisation) and the WWF Germany (financial support). Thanks to all the Cambodian staff and volunteers that assisted to collect samples and conduct necropsies, particularly Lor Kim San, Phay Somany, Sean Kin and Department of Fisheries officers. Our deepest gratitude and thanks to all local villagers throughout the lower Mekong River for their generous hospitality, assistance with reporting and recovering numerous carcasses, and sharing their knowledge about the area and the Mekong dolphin population.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Freshwater cetaceans are among the world’s most threatened mammal species resulting from impacts of multiple human stressors, including habitat loss and degradation, pollution, entanglements in fishing nets and direct killing for bait, traditional medicine purposes and meat and oil for consumption [1–3]. Of the 90 cetacean species worldwide, three species (i.e., Amazon river dolphin Inia geoffrensis, South Asian river dolphin Platanista gangetica and Yangtze river dolphin Lipotes vexillifer) only reside in freshwater habitats (‘obligate’ or ‘true’ river dolphins), whereas three species (Yangtze finless porpoise, Neophocaena asiaorientalis ssp. asiaorientalis, Irrawaddy dolphin, Orcaella brevirostris and tucuxi Sotalia fluviatilis and S. guianensis) have populations in both coastal and riverine environments (‘facultative’ river dolphins). All obligate and facultative river dolphin populations in Asia have suffered dramatic population declines and are listed as Endangered or Critically Endangered on the International Union for Conservation of Nature (IUCN) Red List of Threatened Species. Although the Yangtze River dolphin or baiji is currently listed as Critically Endangered, the species is considered functionally extinct [4–6], which illustrates the vulnerability of these small isolated river dolphin populations. Demographic and genetic data are urgently needed to understand the vulnerability of freshwater dolphins to growing anthropogenic pressures, and contribute information required to prevent further extinctions.
The Irrawaddy dolphin is a small delphinid distributed along Asia’s coastal, lacustrine and riverine waters from the north-western Bay of Bengal southeast to the east coast of Kalimantan, Borneo, and south to the Indonesian island of Sumatra [7]. Lacustrine populations occur in Chilika Lake (India [8]) and Songkhla Lake (Thailand [9]), and riverine populations occur in the Mahakam River (Kalimantan [10]), Ayeyarwardy River (Myanmar [11) and Mekong River (Cambodia and southern Lao People’s Democratic Republic [hereafter referred to as Laos] [12]). Although some coastal Irrawaddy dolphin populations number in the thousands (e.g., Bangladesh [13–14]), most, if not all, lacustrine and riverine populations are small (less than 100 individuals) and face eminent local extinction. Apart from the Chilika Lake population, all lacustrine and riverine populations are classified as Critically Endangered by the IUCN [15–18].
The Irrawaddy dolphin population that inhabits the Mekong River of southern Laos and northern Cambodia (hereafter referred to as the Mekong dolphin population) has been the subject of research since 1973 [19–22], with extensive effort occurring over the past 15 years [23–26]. Based on photo-identification studies, recent population estimates number less than 100 individuals [25,26,12]. The remaining Mekong dolphin population is now primarily restricted to the lower stretch of the Mekong River and some associated tributaries from Khone Falls (five kilometres north of the Laos/Cambodian border) south to Kratie in Cambodia (Fig 1).
Each coloured dot represents the approximate sampling location and mitochondrial DNA haplotype.
There is immediate concern for the Mekong dolphin population’s long-term survival as a result of their low abundance, slow rate of reproduction and continued anthropogenic threats. The main anthropogenic threats are accidental catch in gillnet fisheries [27], un-regulated harassment caused by dolphin-watching vessels [28], direct deaths caused by illegal electric and dynamite fishing methods [27], contaminants and pollutants [29], and dam construction on the mainstream Mekong River in Laos and Cambodia causing significant environmental change and fragmentation of sub-populations [30–31]. An additional concern is a lack of recruitment into the population as a result of unusually high mortality of neonates and calves [25–26]. Little is understood of the genetics of this population, or how genetics studies may support its management. Some recent research suggests distinctness from nearby coastal populations, and signs of adaptation toward the freshwater habitat [32].
Based on photo-identification studies conducted from 2001–2007, there appear to be three main locations within the lower Mekong where Irrawaddy dolphins occur: ‘Cheuteal’, the transboundary pool below Khone Falls on the Laos/Cambodian border, ‘Stung Treng’, from the township of Stung Treng south to Tbong Klar pool, and ‘Koh Pidau/Kampi’, from below Tbong Klar pool south to Kampi Pool located north of the Kratie Township [26,33]. The Mekong dolphin population was estimated to be 93 individuals (95% confidence interval (CI) of 86–101) in 2007, based on research from 2004–2007 [26]. By 2010, the population was estimated to be 85 individuals (95% CI 78–91), based on data from 2007–2010 [25]. The most recent research, which extended the data gathered by [25] to April 2015, estimated 80 individuals, though the uncertainty around this estimate is larger (95% CI 64–100: [34]). Connectivity between the three locations is presumed to be limited, however potential wet season movements remain largely unknown [26].
In this paper, we aim to elucidate the conservation status of the Irrawaddy dolphin population inhabiting the Mekong River by: (a) investigating Irrawaddy dolphin genetic samples from Cambodia in a phylogenetic context; (b) analysing the potential presence and extent of population structure among Irrawaddy dolphins in the Mekong River; and (c) determining demographic parameters of the Irrawaddy dolphin population in the Mekong River.
Materials and methods
Study area
The Mekong dolphin population inhabits a 190 km section of the Cambodian Mekong River, from Kratie to Khone Falls on the Laos/Cambodian border (Fig 1). This river stretch consists of a variety of habitat types, ranging from deep pools (up to 90m) and shallow riffle habitat (i.e. short, relatively shallow and coarse-bedded length of stream over which the stream flows at slower velocity but a higher turbulence than a pool), with numerous islands that separate sections of the river into extensive channels of varying depth and width [15].
Sampling and laboratory procedures
Samples were obtained from the Cambodian cetacean carcass recovery program, which has been operational since 2001 [27,24,33,35]. From 2000–2009, tissue and teeth samples were collected from 67 carcasses, comprising 29 adults (10 females, 8 males, 11 unknown), 35 calves (15 females, 14 males, 6 unknown), two juveniles (both females) and one foetus of unknown sex. Tissue samples were stored in 20% dimethyl sulfoxide (DMSO; [36]) at room temperature in the field and at -20°C upon arrival in the laboratory. Based on sampling locality and DNA quality, individuals sampled in the carcass recovery program and used in this study were assigned to one of the following three sites: Cheuteal (n = 2), Stung Treng (n = 12) Kratie (n = 42) (Fig 1). In addition to these individuals, two individuals were recovered from 100km south of Kratie, near Kampong Cham (added to ‘Kratie’ individuals), one individual was recovered from upstream the Sekong River, Mondulkiri Province (added to ‘Stung Teng’ individuals), and one individual was recovered from Sre Ambel, coastal Cambodia (Fig 1, SI Table 1).
We performed DNA extractions using the Gentra Puregene DNA Purification Kit/DNeasy® Blood and Tissue Kit (Qiagen, Hilden) according to manufacturer’s instructions. For mitochondrial DNA, we amplified a segment of 384-bp containing the hyper-variable region I of the control region by means of the polymerase chain reaction (PCR), using primers dlp1.5 and dlp5 [37]. The PCR products were cleaned using a purification column (Quiagen) according to manufacturer’s instructions. We then sequenced PCR products with the ABI BigDye® Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems), using the forward primer dlp1.5, followed by analysis on an ABI 3730 DNA Analyser sequencer system (Applied Biosystems). Resulting sequences were edited with Sequencing Analysis, version 5.2 (Applied Biosystems). The alignment was carried out manually using the software LaserGene, version 6 (DNASTAR).
We also amplified the same samples using 21 autosomal microsatellites, which had previously shown to be highly polymorphic in different cetacean species: MK3, MK5, MK6, MK9 [38]; 66, 87,91, 98, 105, 108, 111, 117, 128, 138, 141, 142, 153, 162, F10 [39]; D22 [40]; and KWM12 [41]. PCR products were run on an ABI 3730 DNA Analyser. We measured fragment sizes using GeneMapper, version 4.0 (Applied Biosystems) software.
Given the poor quality of the samples, some of which originated from exhumed carcasses in a humid tropical environment, we were concerned about undue influences that allelic dropout could have on our analyses. Thus, we introduced rigid quality control measures to obtain high quality genotypes. First, we amplified each locus at least twice and independently. For individuals that appeared to be homozygous after this procedure, we carried out an additional three, independent PCRs to confirm homozygosity. Second, we also generated genotypes from six individuals independently to estimate genotyping error rates for each locus.
Phylogenetic analysis
To assess phylogenetic relationships among dolphins from the Mekong River and other taxa of the genus Orcaella, we carried out a phylogenetic analysis using a Bayesian framework. Mitochondrial DNA sequences were obtained from GenBank from Orcaella brevirostris from India, the Philippines, Indonesia, Thailand, and the Cambodian coast, as well as from Australian snubfin dolphin Orcealla heinsohni from Australia [42]. As outgroups, we used mtDNA sequences from Orcinus orca and Steno bredanensis. We performed the analysis using program MrBayes version 3.2.2 [43–44]. Instead of selecting a substitution model, we sampled across the general time reversible (GTR) model space in the Bayesian Markov Chain Monte Carlo (MCMC) analysis itself, obviating the need to test for models a priori. We ran 20 x 106 generations with a sample frequency of 10,000, discarding the first 25% of saved trees. Convergence of runs was determined by tools built into MrBayes, i.e. standard deviations of split frequencies below 0.01.
Population genetic analyses
Mitochondrial DNA.
Using mitochondrial DNA, we assessed genetic diversity within sampling localities by calculating haplotype and nucleotide diversity using Arlequin version 3.5.1.3 [45–46]. To test the degree of population differentiation among sampling localities, we calculated ɸST values using Arlequin and applied a Tamura and Nei [47] with a γ-correction of 0.144, as computed by JModeltest version 2.1.3 [48]. Significance was obtained by performing 10,000 permutations.
Autosomal microsatellite data.
The software package Structure version 2.3.3 [49] was used to determine the genetic structure and number of genetic clusters in our dataset. The Structure algorithm divides sampled individuals into a number of clusters (K) independent of locality information and by minimizing deviations from Hardy–Weinberg and linkage equilibrium in each cluster. The program uses a MCMC procedure to estimate P(X|K), the posterior probability that the data fit the hypothesis of K clusters. For all analyses, the length of the burn-in period was set to 105, followed by 106 MCMC steps. For each K, the analysis was run ten times. We chose the Locprior model, which improves clustering when the signal is weak without spuriously inferring structure, if absent [50].
Using data generated from the autosomal microsatellites, we estimated genetic variation within and genetic differentiation between sampling sites. Prior to this, we tested for significant deviations from Hardy-Weinberg equilibrium (HWE) expectations for each locus and population, as well as the presence of linkage disequilibrium (LD) among all loci, using the software Genepop, version 4.2.1 [51–52]. A Bonferroni correction for multiple tests was subsequently applied for each population and locus [53].
Genetic variation within sampling locations was determined by estimating the number of alleles per locus (Na), number of effective alleles per locus (Ne), and calculating unbiased levels of expected (uHe) and observed heterozygosity (Ho), using the software GenAlEx, version 6.5 [54]. To estimate genetic differentiation between sampling locations, we calculated Weir and Cockerham’s θ [55] in Genepop and Dest [56–57] using Smogd, version 1.2.5 [58]. In addition, the number of private alleles per population was calculated as described in [59] with GenAlEx.
Demographic modelling
We detected change in the effective population size using three approaches. The first two used summary statistics to detect a change in allele frequencies while the last one used a full likelihood Bayesian method that also allowed us to detect, quantify and date the change in size. The first approach was implemented in Arlequin and developed by [60]. This approach may identify bottlenecks if they had lasted several generations, the pre-bottleneck θ was large, and the population had made a demographic recovery. The approach takes advantage of the fact that reductions in population size can be used to estimate the mean ratio of the number of alleles to the range in allele size [60] using a statistic Mk. The value of Mk has been shown to decrease when a population is reduced in size over recent times (<10 generations) and as such, the statistic Mk can distinguish between populations that have been recently reduced in size to those that have been small for a long time [60]. An expected distribution for Mk under equilibrium conditions was generated using simulations. The critical value is at the lower 95th percentile of this distribution, and a bottleneck identified when a value of Mk is lower than this critical value [60].
The second approach was developed by [61], and inferred genetic bottlenecks as a result of an excess of heterozygotes under a stepwise mutation model. We employed the software Bottleneck, version 1.2 [62] to test for distortion of allele frequency distributions as a result of rare alleles being more likely to be lost during a bottleneck than common alleles. This test for a genetic bottleneck is appropriate for populations that have been reduced very recently with little severity, and a small pre-bottleneck value of θ [63]. We tested three models of microsatellite mutation (Infinite Allele Model—IAM, Two-Step Model—TPM, and the Stepwise Mutation Model—SMM). For the TPM, we assumed that single step mutations account for 90% of all mutation events, and a variance among multiple steps of 12, as suggested by Piry et al. (1999). Due to the relatively small number of loci analysed (n = 10), we used a Wilcoxon sign-rank test to estimate significance [62]. Despite its use, we believe that with 10 polymorphic loci, low heterozygosity, and a significant reduction in pre/post bottleneck size (data from our results; see [64]), there will be only very low power to accumulate a detectable excess.
We also employed a Bayesian approach developed by [65] for detecting changes in population size. This approach makes use of coalescent theory to create posterior distributions from different prior demographic parameters, given the observed allele frequency distributions at the sampled loci. The main assumptions of the approach are a population of constant size of N1, followed by an increase (or decrease) ta generations ago to the current population size of N0. Changes are assumed to be either linear or exponential at a rate of θ = 2N0μ (μ = mutation rate) under a stepwise mutation model.
First, we employed a Bayesian coalescent-based MCMC approach implemented in the software MsVar version 0.4.2 to estimate posterior probabilities of the size of population size change (r = N0/N1), the timing of when the change occurred (tf = ta/N0—a measure scaled by N0), and the scaled mutation rate θ. We ran five independent chains with different starting demographic histories, using wide uniform priors of between -5 and 5 (on a log10-scale) for θ, r and tf.
With the software MsVar version 1.3, we also used a method by [66] to quantify the effective population size of the current (N0) and the ancestral (N1) population, as well as the time (T) since the population size changed. In this model, the prior distributions of θ, N0, and T are assumed to be log-normal. As such, we used wide uninformative priors, multiple runs (at least three independent runs as a test of convergence, with differing starting seeds) and varying hyperprior parameters. MsVar was run for 180 simulations with 109 iterations each using a cluster array at http://www.ivec.org/. Current population size was set to N0 = 101 [26]. We used the Irrawaddy dolphin generation time of 20 years based on [67].
Outputs for MSVAR were interpreted using purpose-written scripts in R v2.10.1 (R Development Core Team 2009). We discarded the initial 10% of sampled points for each run and checked convergence of chains visually and formally using the Gelman-Rubin statistic [67–69], with R library: boa [70]) for Log10 (θ, r and tf). It has been shown that a point estimate (and 95% quartile) of <1.2 is indicative of satisfactory convergence [71]. Approximate posterior densities for 0.1, 0.5 and 0.05 (highest posterior density limits; HPD) were calculated and graphically represented (R library: boa, hdrcde and Locfit).
Results
The phylogenetic analysis indicates reciprocal monophyly of the Mekong dolphin population haplotypes with respect to haplotypes from other populations, suggesting a long-standing isolation of the Mekong dolphin population from other Orcaella populations (Fig 2). Posterior probability values support a clear split between O. brevirostris from the Mekong River and populations from other areas, including coastal Cambodia, approximately 700km from the Mekong Delta.
Numbers indicate Bayesian posterior probability values for each clade.
Mitochondrial haplotype and nucleotide diversity indices were estimated from 55 individuals, from which we were able to obtain mtDNA sequence data. Based on their origin, individuals for this analysis represented three sampling regions: Kratie (n = 42), Stung Treng (n = 12) and Cheuteal (n = 1). In the Kratie and Stung Treng regions, haplotype D was the most common, with at least 85% of all individuals in each region bearing the same haplotype (Table 1). There was only data for one individual from the Cheuteal region, which also showed the most common haplotype D. The low haplotype diversity was reflected in the nucleotide diversity at each site, which was extremely low (Table 1). The AMOVA revealed a complete lack of genetic structure for mitochondrial DNA within our dataset (Table 2). Only a non-significant (P = 0.32), very small portion (1.73%) of the total genetic variation was attributed to genetic differences between the sampling sites.
Given the challenging sample quality, we were able to generate microsatellite genotypes for 44 individuals. Out of the 21 loci tested, eleven were monomorphic (Table 3). Over all 21 loci, expected heterozygosity values were extremely low (Kratie: uHe21 = 0.22, Stung Treng uHe21 = 0.19, combined: uHe21 = 0.21). The number of private alleles was higher in Kratie (n = 6) than in Stung Treng (n = 1), although this was expected given the almost threefold difference in sample size between the two.
Eleven loci were monomorphic, thus, summary statistics are also presented for the ten polymorphic loci. N = mean number of samples/locus; Na = average number of alleles/locus; Ne = mean number of effective alleles/locus; Ho = mean observed heterozygosity; uHe = mean unbiased expected heterozygosity; S.E. = standard error. Numerical indices in the top refer to number of microsatellite loci on which summary statistics are based.
Genetic differentiation between all sampling localities was low. Both FST and Jost’s D were small (FST = 0.044, D = 0.004) and non-significant. The Structure analysis corroborated these findings, as we did not detect any population structure (Log-likelihood values for five independent runs for each K indicated that K = 1 had the highest probability) even when using a model with the highest chance to detect weak structure. Combining these results and the lack of differentiation found mitochondrial DNA, we considered both sampling locations as one population for subsequent analyses of demographic change.
Evidence about recent bottlenecks was equivocal. Estimates of Mk (the mean ratio over all loci of the number of alleles over the range in allele size) were low, ranging from 0.26 to 0.67 (Mk10 loci = 0.39, s.d. = 0.13). However, no values were significantly lower than those expected by chance (P>0.05), even under wide θ starting–values and current effective population size spanning 100 to 10,000 dolphins. In the Bottleneck analysis, we did not detect a deviation from a normal L-shaped distribution of alleles, and no significant bottleneck using a Wilcoxon rank sum test under a SMM (P = 0.55), TPM (P = 1.00), or IAM (P = 0.19) model of mutation, as predicted (SI Table 2).
Posterior distributions of the rate of population change in the Mekong dolphin population, as calculated with MsVar, version 0.4.2, suggested that there has been a recent exponential decline in the population (Fig 3). All models invoking a linear change of population size did not converge, as indicated by the Gelman statistic [65], while all runs with exponential changes converged and were remarkably consistent, despite a wide range of prior parameters (SI Table 3). The median log10(r) ratio of current (N0) to ancestral (N1) effective population sizes was -1.20 (95% HDR -2.85–0.18) across all runs (SI Table 4), indicating that a decline in the effective size of all populations of about one order of magnitude and no support for stability (log10(r) = 0) or expansion (log10(r) > 0). The decline occurred fairly recently, as median log10tf (2.2) was small.
Results from the MsVar, version 1.3 runs, using the method by [66] to quantify changes in the effective populations size, suggested an exponential decline of the Mekong dolphin population of about 19 times its original source size (Fig 3, SI Table 4). Similar to simulations from MsVar 0.4.2, none of the linear runs converged, as indicated by the Gelman statistic [65]. The posterior distributions for N0 and N1 showed only a moderate overlap, confirming the occurrence of a population decline. Current N0 (median log N0 = 1.89, 95% HDR = 0.11–3.15) was more than an order of magnitude smaller than N1 (median log N1 = 3.17, 95% HDR = 1.57–4.91). These data suggest that the historic population size was at least 19 times larger than the current one, or that the current effective population size is only 5.2% of that of the ancestral population. Assuming a generation time of 20 years for Mekong River Orcaella, the decline might have started ca. 145,000 years ago (95% HDR 159–1,584,893 years).
Discussion
The Irrawaddy dolphin population occupying the Mekong River is small, declining, isolated, and potentially suffering from inbreeding and loss of genetic diversity. It is apparent that this population is in serious threat of extirpation in the near future. Of the three riverine populations of Irrawaddy dolphins, the population that inhabits the Mekong River is likely the largest and certainly the most extensively studied. Thus, the capacity for this population to recover, and our ability to assist with its recovery, are potentially greater than the other riverine populations of Irrawaddy dolphins.
Our phylogenetic analysis indicated reciprocal monophyly of Mekong dolphin population haplotypes with respect to haplotypes from other populations, suggesting a long-standing isolation of the Mekong dolphin population from other Irrawaddy dolphin populations, including coastal Cambodia. However, further Irrawaddy dolphin genetic samples, particularly from coastal areas near the Vietnamese Mekong Delta and other riverine populations, are required in order to elucidate whether sub-species level designation is justified. Moreover, to increase resolution and support, data from full mitochondrial genomes should be generated. Regardless of the need for further samples, our limited dataset indicates that the Mekong dolphin population conservation efforts should be expanded and intensified as a matter of priority, to reduce the local extinction probability of this potentially unique sub-species.
Our data suggest that the different sub-populations along the Mekong River should be treated as a panmictic population. However, even at a combined level, genetic variability is extremely low. The causes for this could be manifold, and unfortunately our data are insufficient to clearly disentangle different scenarios of how the low genetic diversity originated. Given the limited occurrence of Mekong River Irrawaddy dolphins and their genetic isolation, a potential explanation is that the founding population is fairly young and had a small effective population size to start with. Moreover, given the slow generation time of 20 years [67], this population might have never accumulated a substantial amount of genetic diversity via mutation and/or immigration.
Being part of the Sunda shelf, the Mekong River system was subjected to repeated geological transformations during multiple glacial periods during the Pleistocene (between 2.58 Mya and 11.8 kya: [72,42,73–74]). Consecutive cycles of cooling and warming resulted in different sea levels, so that land-masses and river systems fluctuated greatly in size and geographic extent over time [75]. These changes in length and depth of different river stretches would certainly have had consequences for the connectivity of different river stretches to one another, as well as for the coastal area and deep water pools within the upper Cambodian Mekong River. Based on these complex geological features and the broad confidence intervals of our demographic analyses, it is unfortunately not possible to determine whether the current low genetic diversity observed in the Mekong dolphin population is an ancient effect.
In recent times, a significant loss of individuals was probably caused by the intensive bombing of Laos, Cambodia and Vietnam during wars from the 1950s to 1970s, as well as intentional killings of dolphins in Tonle Sap Lake for their oil during the same period [21–22,76]. There are no estimates on the number of dolphins killed, but the area along the eastern Cambodian and Vietnamese border, and especially the Mekong River, was considered a potential supply route of the Vietcong, and struck by more than 2.7 million tons of explosives between 1964 and 1975 [76]. However, given the slow life history of Irrawaddy dolphins, the observed extremely low current genetic diversity cannot be entirely explained by the substantial number of dolphins killed during the 1960s and 1970s.
Conservation considerations
To ensure both the long-term persistence and evolutionary potential of the Mekong dolphin population, drastic management measures are required immediately by the relevant Governments, conservation agencies and local Cambodian and Laos communities.
Since the 1970’s, dolphins have been shot for target practice, accidentally by-caught in subsistence fisheries, and have been directly killed through dynamite and electric fishing [20–22,33,77]. Recent anthropogenic threats include potential contaminants from gold-mining and agricultural run-off [78] and potential harassment from dolphin-watching vessels [28]. The high newborn mortality is an additional concern, with no apparent cause being discovered despite extensive investigations [79,24]. Since 2009, an ad hoc team of international scientists have been working with the Cambodian Government and WWF-Cambodia to co-ordinate and develop Mekong dolphin conservation efforts, with the landmark ‘Kratie Declaration on the Conservation of the Mekong River Irrawaddy dolphins’ being finalized in 2012 [79]. Although progress has been made with patrolling activities and the refinement of protocols to investigate mortality, entanglement in gillnets is reportedly one of the most critical and immediate threat to survival of dolphins in the Mekong River [80–81,26]
The construction and operation of the Don Sahong hydropower dam on the Laos/Cambodian border is also a major concern, particularly to the small sub-population of dolphins (now numbering only three individuals: [34] that inhabit the Laos/Cambodian border region of Cheuteal [81,31]. Construction and operation of the Don Sahong dam would almost certainly cause local extirpation of the already small Cheuteal sub-population as a result of increased noise (particularly from explosives used during the construction phase), increased boat traffic, and reduced prey from blockage of a major fish passage during the dry season [31,82]. The proposed Sambor Dam, which would be constructed across the mainstream Mekong River north of Kratie, is of significant concern due to its likelihood of fragmenting the remaining Mekong dolphin population and potential environmental impacts [83,84,85]. Construction of these dams increases extinction risk of the entire Mekong dolphin population [31,82]. Fragmentation of river dolphin meta-populations by dams and irrigation barrages is an increasing problem for the long-term survival of river dolphin populations world-wide [11,86–87], where loss of connectivity, reduced water and prey availability, and close proximity to anthropogenic activities accentuates the difficulties of river dolphin conservation [28].
Gene flow among small fragmented populations is critical for maintaining genetic diversity, and therefore the evolutionary potential of a species [88]. When gene flow is restricted, low abundance sub-populations, or population fragments, become increasingly vulnerable to loss of genetic variation via genetic drift, accumulation of deleterious mutations, inbreeding depression, and inability to adapt to change [89], all of which are correlated with the risk of extinction. Alleviating inbreeding depression by transferring new breeding stock into severely inbred and isolated populations (‘genetic rescue’), as carried out in panthers [90,91], prairie chickens [92], adders [93], and ibex [94], are not considered a viable options for cetacean populations.
Given that ‘genetic rescue’ is not feasible for the Mekong dolphin population, or cetaceans in general, management efforts should focus on effectively managing all locations where Mekong dolphins occur (i.e. Cheuteal, Stung Treng, Koh Pdao, and Kratie); with the Cheuteal population being a high priority, given its small size of three individuals and isolated location 70km north of the Steng Treng sub-population. Important considerations are ensuring continued connectivity between sub-populations as well as reduction of anthropogenic threats, such as accidental entanglement in gillnets and boat-based dolphin watching tourism by changing to land-based tourism. Population monitoring is an essential component of any species recovery program, where photo-identification to assess trends in abundance, and carcass recovery programs to monitor mortality rates and causes, is required to monitor the remaining population and changes in response to management initiatives [25,26].
Captive breeding and reintroduction have not played a major role in the conservation of small cetaceans [95]. At present, the prospects of establishing an Irrawaddy dolphin captive breeding and reintroduction program are not recognised as viable options because of the inherent risks associated with the capture process, a lack of reliable care and husbandry techniques for Irrawaddy dolphins, the known poor survival rate for Irrawaddy dolphins in captivity [96], and the lack of substantial Cambodian/Vietnamese coastal Irrawaddy dolphin populations from which a large number of translocatable individuals could be taken. There is also a lack of knowledge about possible local adaptions required for the translocated individuals, such as prey requirements, site fidelity or social structure [33]. Importantly, as discussed in [91], there is little justification to ‘rescue’ small isolated populations if known threats are continuing.
The potential for extinction is now a reality for the Mekong dolphin population, particularly with construction of the Don Sahong Dam, and potential construction of the Sambor Dam. This study suggests a long-standing isolation of the Mekong dolphin population from other Orcaella populations, with the remaining population now experiencing low genetic diversity. In order for the Mekong dolphin population to have any potential for long-term survival, it is imperative that any further population fragmentation or genetic loss is prevented, anthropogenic threats are reduced, and continued connectivity between all regions is a high management priority.
Supporting information
S1 Fig. Permission to publish Fig 1 under a CC BY license.
https://doi.org/10.1371/journal.pone.0189200.s001
(XLSX)
Acknowledgments
The authors would like to thank Gordon Congdon (WWF-Greater Mekong Program) for his significant contribution to Mekong dolphin research and conservation efforts, and Verne Dove for her assistance with the stranding program and collection of samples used in this analysis from 2007 onwards. Many thanks to Martin Gilbert (WCS- Cambodia Program) who managed the Cambodian stranding program from 2005–2007. Thanks to all the Cambodian staff and volunteers that assisted to collect samples and conduct necropsies, particularly Channa Phan, Lor Kim San, Phay Somany, Sean Kin and Department of Fisheries officers. Particular thanks to the Cambodian Department of Fisheries, especially His Excellency Nao Thouk (Director) for permission to research the Mekong dolphin population since 2001. Our deepest gratitude and thanks to all local villagers throughout the lower Mekong River for their generous hospitality, assistance with reporting and recovering numerous carcasses, and sharing their knowledge about the area and the Mekong dolphin population.
References
- 1. Jefferson TA, Hung SK, Würsig B. Protecting small cetaceans from coastal development: Impact assessment and mitigation experience in Hong Kong. Mar Policy. 2009;33(2):305–311.
- 2. Smith BD, Reeves RR. River cetaceans and habitat change: Generalist resilience or specialist vulnerability? J Mar Biol. 2012.;
- 3. Gomez-Salazar C, Coll M, Whitehead H. River dolphins as indicators of ecosystem degradation in large tropical rivers. Ecol. Indic. 2012;23:19–26.
- 4. Turvey ST, Pitman RL, Taylor BL, Barlow J, Akamatsu T, Barrett LA, et al. First human-caused extinction of a cetacean species? Biol. Lett. 2007;3(5):537–540. pmid:17686754
- 5. Wang D, Zhang X, Wang K, Wei Z, Wursig B, Braulik GT, et al. Conservation of the Baiji: No simple solution. Conserv. Biol. 2006;20(3):623–625. pmid:16909544
- 6. Reeves RR, Gales N. Realities of Baiji conservation. Conserv. Biol. 2006;20(3):626–628. pmid:16909545
- 7.
Stacey PJ, Arnold P. Orcaella brevirostris. Mamm.o Species. 1999;1–8.
- 8. Sutaria D, Marsh H. Abundance estimates of Irrawaddy dolphins in Chilka Lagoon, India, using photo-identification based mark-recapture methods. Mar. Mamm. Sci. 2011;27(4):338–348.
- 9. Beasley IL, Chooruk S, Piwpong N. The status of the Irrawaddy dolphin, Orcaella brevirostris, in Songkhla Lake, Southern Thailand. Raffles Bull. Zool. 2002;75–84.
- 10. Kreb D. Abundance of freshwater Irrawaddy dolphins in the Mahakam River in East Kalimantan, Indonesia, based on mark-recapture analysis of photo-identified individuals. J. Cetacean Res. Manag. 2004;6(3):269–277.
- 11. Smith BD, Hobbs L. Status of Irrawaddy dolphins Orcaella brevirostris in the upper reaches of the Ayeyarwady River, Myanmar. Raffles Bull. Zool. 2002;10:67–74.
- 12. Beasley IL, Pollock K, Jefferson TA, Arnold P, Morse L, Yim S, et al. Likely future extirpation of another Asian river dolphin: The critically endangered population of the Irrawaddy dolphin in the Mekong River is small and declining. Mar. Mamm. Sci. 2013;29(3):226–252.
- 13. Smith BD, Ahmed B, Mowgli RM, Strindberg S. Species occurrence and distributional ecology of nearshore cetaceans in the Bay of Bengal, Bangladesh, with abundance estimates for Irrawaddy dolphins Orcaella brevirostris and finless porpoises Neophocaena phocaenoides. J. Cetacean Res. Manag. 2008;10(1):45–58.
- 14. Smith BD, Braulik GT, Strindberg S, Ahmed B, Mansur R. Abundance of Irrawaddy dolphins (Orcaella brevirostris) and Ganges river dolphins (Platanista gangetica gangetica) estimated using concurrent counts made by independent teams in waterways of the Sundarbans mangrove forest in Bangladesh. Mar. Mamm. Sci. 2006;22(3):527–547.
- 15.
Smith BD, Beasley IL. Orcaella brevirostris (Mekong River subpopulation). 2004 IUCN Red List of threatened species, www.redlist.org. 2004a; Downloaded on 10 September 2016; http://www.iucnredlist.org/details/44555/0
- 16.
Smith BD, Beasley IL. Orcaella brevirostris (Songkhla Lake subpopulation). 2004 IUCN Red List of threatened species, www.redlist.org. 2004b; Downloaded on 10 September 2016; http://www.iucnredlist.org/details/44557/0
- 17.
Kreb D, Smith BD. Orcaella brevirostris (Mahakam subpopulation). 2004 Red List of Threatened Species, www.redlist.org. 2000; Downloaded on 10 September 2016; http://www.iucnredlist.org/details/39428/0
- 18.
Smith BD. Orcaella brevirostris (Ayeyarwady River subpopulation). 2004 IUCN Red List of Threatened Species, www.redlist.org. 2004; Downloaded on 10 September 2016; http://www.iucnredlist.org/details/44556/0
- 19.
Lloze R. Contribution a l’etude anatomique, histologique et biologique de l’Orcaella brevirostris (Gray 1866) (Cetacea, Delphinidae) du Mekong. Dissertation, L’Universite Paul Sabatier de Toulouse, Toulouse, France; 1973.
- 20. Baird IG, Beasley IL. Irrawaddy dolphin Orcaella brevirostris in the Cambodian Mekong River: an initial survey. Oryx 2005;39(3):301–310.
- 21. Baird IG, Mounsouphom B. Irrawaddy dolphins (Orcaella brevirostris) in southern Lao PDR and northeastern Cambodia. Nat. Hist. Bull. Siam. Soc. 1994;42(2):59–175.
- 22. Baird IG, Mounsouphom B. Distribution, mortality, diet and conservation of Irrawaddy dolphins (Orcaella brevirostris Gray) in Lao PDR. Asian Mar. Biol. 1997;14:41–48.
- 23. Stacey PJ, Hvenegaard GT. Habitat use and behavior of Irrawaddy dolphins (Orcaella brevirostris) in the Mekong River of Laos. Aquat. Mamm. 2002;28(1):1–13.
- 24.
Dove V. Mortality investigation of the Mekong Irrawaddy River dolphin (Orcaella brevirostris) in Cambodia based on necropsy sample analysis. WWF Cambodia Technical Report; 2009.
- 25.
Dove V, Truijillo F, Zanre R. Abundance estimation of the Mekong Irrawaddy dolphin Orcaella brevirostris based on mark and recapture analysis of photo-identified individuals. WWF Cambodia Technical Report; 2008.
- 26. Ryan GE, Dove V, Trujillo F, Doherty PF. Irrawaddy dolphin demography in the Mekong River: An application of mark-resight models. Ecosphere 2011;2(5):1–14.
- 27.
Gilbert M, Beasley IL. Mekong River Irrawaddy dolphin stranding and mortality summary: January 2001-December 2005. Wildlife Conservation Society-Cambodia Program, Phnom Penh; 2006.
- 28.
Beasley IL, Bejder L, Marsh H. Concerns about cetacean-watching in developing countries: a case study from the Mekong River. In: Higham J, Bejder L, Williams R, editor. Whale-watching, sustainable tourism and ecological management. Cambridge University Press, Cambridge; 2014. p. 242–262.
- 29.
Tran, T., Barzen, J., Choowaew, S., Engels, M., Duong, V.N., Nguyen, et al. Persistent organic pollutants in wetlands of the Mekong Basin: U.S. Geological Survey Scientific Investigations Report 2013–5196. 2014;140 p., http://dx.doi.org/10.3133/sir20135196.
- 30.
Ryan GE. Last chance for dolphins in Laos: A review of the history, threats and status. Technical Report from WWF-Greater Mekong Programme; 2012.
- 31.
Ryan GE. The Don Sahong Dam and the Mekong dolphin: An updated review of the potential impacts of the Don Sahong hydropower project on the Mekong River’s critically endangered Irrawaddy dolphins, Orcaella brevirostris). WWF Science Brief; 2014.
- 32. Caballero S, Duchêne S, Garavito MF, Slikas B, Baker CS. Initial evidence for adaptive selection on the NADH subunit two of freshwater dolphins by analyses of mitochondrial genomes. PLoS ONE. 2015;10(5):1–17. https://doi.org/10.1371/journal.pone.0123543
- 33.
Beasley IL. Conservation of the Irrawaddy dolphin, Orcaella brevirostris (Owen in Gray, 1866) in the Mekong River: Biological and social considerations influencing management. Dissertation, School of Earth and Environmental Science, James Cook University, Townsville, Australia; 2017.
- 34.
Phan, C., Hang, S., Tan, S. and Lor, K. Population monitoring of the critically endangered Mekong dolphin based on mark-resight models. WWF-Cambodia Technical Report; 2015.
- 35. Beasley IL, Davidson PJ. Conservation status of marine mammals in Cambodian waters, including seven new cetacean records of occurrence. Aquat. Mamm. 2007;33:368–379.
- 36.
Amos W, Hoelzel AR. Long-term preservation of whale skin for DNA analysis. In: Hoelzel AR, Donovan GP (ed) Genetic Ecology of Whales and Dolphins, Report to the International Whaling Commission, Cambridge, UK; 1991.
- 37. Baker C, Gilbert D, Weinrich MT, Lambertsen R, Calambokidis J, McArdle B, et al. Population characteristics of DNA fingerprints in humpback whales (Megaptera novaeangliae). J. Hered. 1993;84(4):281–290. pmid:8340617
- 38. Krützen M, Valsecchi E, Connor RC, Sherwin WB. Characterization of microsatellite loci in Tursiops aduncus. Mol. Ecol. Notes 2001;1(3):170–172.
- 39. Nater A, Kopps AM, Krützen M. New polymorphic tetranucleotide microsatellites improve scoring accuracy in bottlenose dolphins (Tursiops aduncus). Mol. Ecol. Resour. 2009;9:531–534.0. pmid:21564684
- 40. Shinohara M, Domingo-Roura X, Takenaka O. Microsatellites in the bottlenose dolphin Tursiops truncatus. Mol. Ecol. 1997;6(7):695–696. pmid:9226950
- 41. Hoelzel AR, Dahlheim M, Stern SJ. Low genetic variation among killer whales (Orcinus orca) in the eastern north Pacific and genetic differentiation between foraging specialists. J. Hered. 1998;89(2):121–128. pmid:9542159
- 42. Beasley IL, Robertson K, Arnold P. Description of a new dolphin, the Australian snubfin dolphin Orcaella heinsohni sp.n. (Cetacean, Delphinidae). Mar. Mamm. Sci. 2005;21(3):365–400.
- 43. Huelsenbeck JP, Ronquist F. MrBayes: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17(8):754–755. pmid:11524383
- 44. Ronquist F, Huelsenbeck J. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003;19(12):1572–1574. pmid:12912839
- 45. Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol. Bioinformatics 2005;1:47.
- 46.
Schneider S, Roessli D, Excoffier L. Arlequin: a software for population genetics data analysis, version 2.000. Genetics and Biometry Laboratory, University of Geneva, Switzerland; 2000.
- 47. Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol. Biol. Evol. 1993;10(3):512–526. pmid:8336541
- 48. Posada D. jModelTest: Phylogenetic model averaging. Mol. Biol. Evol. 2008;25(7):1253–1256. pmid:18397919
- 49. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–959. pmid:10835412
- 50. Hubisz MJ, Falush D, Stephens M, Pritchard J. Inferring weak population structure with the assistance of sample group information. Mol. Ecol. Resour. 2009;9(5):1322–1332. pmid:21564903
- 51. Raymond M, Rousset F. GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J. Hered. 1995;86(3):248–249.
- 52. Rousset F. Genepop’007: a complete re-implementation of the genepop software for Windows and Linux. Mol. Ecol. Resour. 2008;8(1):103–106. pmid:21585727
- 53. Rice WR. Analyzing tables of statistical tests. Evolution 1989;43(1):223–225. pmid:28568501
- 54. Peakall R, Smouse PE. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics 2012;28(19):2537–2539. pmid:22820204
- 55. Weir BS, Cockerham CC. Estimating F-Statistics for the analysis of population structure. Evolution 1984;38(6):1358–1370. pmid:28563791
- 56. Hedrick PW. A standardized genetic differentiation measure. Evol. 2005;59(8):1633–1638.
- 57. Jost L. GST and its relatives do not measure differentiation. Mol. Ecol. 2008;17(18):4015–4026. pmid:19238703
- 58. Crawford NG. SMOGD: Software for the measurement of genetic diversity. Mol. Ecol. Resour. 2010;10(3):556–557. pmid:21565057
- 59. Slatkin M. Rare alleles as indicators of gene flow. Evol. 1985;39(1):53–65.
- 60. Garza JC, Williamson EG. Detection of reduction in population size using data from microsatellite loci. Mol. Ecol. 2001;10(2):305–318. pmid:11298947
- 61. Luikart G. Distortion of allele frequency distributions provides a test for recent population bottlenecks. J. Hered. 1998;89:238–247. pmid:9656466
- 62.
Piry S, Luikart G, Cornuet JM. BOTTLENECK: A program for detecting recent effective population size reductions from allele data frequencies. Montpellier, France; 1999.
- 63. Williamson-Natesan EG. Comparison of methods for detecting bottlenecks from microsatellite loci. Conserv. Genet. 2005; 6(4): 551–562.
- 64. Cornuet JM, Luikart G. Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics. 1996;144(4):2001–2014. pmid:8978083
- 65. Beaumont MA. Detecting population expansion and decline using microsatellites. Genetics. 1999; 153(4):2013–2029. pmid:10581303
- 66. Storz JF, Beaumont MA. Testing for genetic evidence of population expansion and contraction: an empirical analysis of microsatellite DNA variation using a hierarchical Bayesian model. Evolution. 2002;56(1):154–166. pmid:11913661
- 67.
Moore J. Intrinsic growth (rmax) and generation time (T) estimates for odontocetes of the genus Sousa, Orcaella, and Neophocaena, in support of IUCN Red List Assessments. NOAA Technical Memorandum; 2015.
- 68. Gelman A, Rubin DB. Inference from iterative simulation using multiple sequences. Stat. Sci. 1992;457–472.
- 69. Brooks SP, Gelman A. General methods for monitoring convergence of iterative simulations. J. Comp. Graph. Stat. 1998; 7(4): 434–455.
- 70. Smith B.J. boa: An R Package for MCMC output convergence assessment and posterior inference. J. Stat. Software 2007; 21(11): 1–37.
- 71.
Hill J, Gelman A. Data analysis using regression and multilevel/hierarchical models. Cambridge University Press; 2007.
- 72.
Pillans B, Gibbard P. The Quaternary Period. In: Gradstein F, Ogg GJ, Schmitz M, Ogg GM editors The Geologic Time Scale, 1st ed. Elsevier, Amsterdam; 2012; p. 979–1004.
- 73.
Holloway JD, Hall R. Southeast Asian geology and biogeography: An introduction. In: Hall R, Holloway JD, editors. Biogeography and Geological Evolution of Southeast Asia. Backhuys, Leiden; 1998. p 1–23
- 74. Woodland DJ. Zoogeography of the Siganidae (Pisces): An interpretation of distribution and richness patterns. Bull. Mar. Sci. 1983;33(3):713–717.
- 75. Voris HK. Maps of Pleistocene sea levels in Southeast Asia: Shorelines, river systems and time durations. J. Biogeogr. 2000;27(5):1153–1167.
- 76.
Owen T, Kiernan B. Bombs over Cambodia. Third World Resurgence. 2007;201:41.
- 77. Jayasankar P, Patel A, Khan M, Das P, Panda S. Mitochondrial DNA diversity and PCR-based sex determination of Irrawaddy dolphin (Orcaella brevirostris) from Chilika Lagoon, India. Mol. Biol. Rep. 2011;38(3):1661–1668. pmid:20857220
- 78. Murphy PT, Irvine KN, Sampson M, Guo J, Parr T. Mercury contamination along the Mekong River, Cambodia. Asian Journal of Water, Environment and Pollution. 2007;6(1): 1–9.
- 79.
Reeves RR, Bownell RL, Gulland F, Smith BD, Turvey ST, Ding W. Assessment of mortality of Irrawaddy dolphins in the Mekong River and recommendations for a population recovery plan. Unpublished report submitted to WWF-Cambodia; 2009.
- 80.
McLellan WA. General observations of Cambodian river dolphin images. Unpublished report to WWF Cambodia Program; 2012.
- 81.
Thuok N, Ath CA, Reeves RR. Conclusions, recommendations and working group reports from the workshop on the conservation of Irrawaddy dolphins in the Mekong River. Unpublished report submitted to the IUCN Cetacean Specialist Group; 2014.
- 82.
Bezuijen MR, Zanre R, Goichot M. The Don Sahong Dam and the Irrawaddy dolphin. WWF Greater Mekong Programme, Vientiane; 2007.
- 83.
CSG. The Kingdom of Cambodia Feasibility Study of Sambor Hydropower Station. China Southern Power Grid Co. Ltd; 2008.
- 84.
Natural Heritage Institute. Sambor hydropower dam alternatives assessment: interim report. Unpublished report submitted to the Royal Government of Cambodia; 2017.
- 85. Brownell RL Jr, Randall R, Smith BD Ryan GE. Dams threaten rare Mekong dolphins. Science 2017;355(6327):805.
- 86. Braulik GT, Arshad M, Noureen U, Northridge SP. Habitat fragmentation and species extirpation in freshwater ecosystems; causes of range decline of the Indus River dolphin (Platanista gangetica minor): PLoS ONE 2014; e101657. pmid:25029270
- 87. Araújo CC, Wang JY. The dammed river dolphins of Brazil: impacts and conservation. Oryx 2015;49(01):17–24.
- 88. Hamner RM, Pichler FB, Heimeier D, Constantine R, Baker SC. Genetic differentiation and limited gene flow among fragmented populations of New Zealand endemic Hector’s and Maui’s dolphins. Conserv. Genet. 2012;13(4):987–1002.
- 89. Frankham R. Conservation Genetics. Annu. Rev. Genet. 1995;29(1):305–327.
- 90. Maehr DS, Lacy RC. Avoiding the lurking pitfalls in Florida panther recovery. Wildlife Society Bulletin. 2002;30:971–978.
- 91. Pimm SL, Dollar L, Bass OL. The genetic rescue of the Florida panther. Anim. Conserv. 2006;9(2):115–122.
- 92. Westemeier RL, Brawn JD, Simpson SA, Esker TL, Jansen RW, Walk JW, et al. Tracking the long-term decline and recovery of an isolated population. Science 1998;282(5394):1695–1698. pmid:9831558
- 93. Madsen T, Ujvari B, Olsson M. Novel genes continue to enhance population growth in adders (Vipera berus). Biol. Conserv. 2004;120(1):145–147.
- 94. Turcek FJ, Hickey JJ. Effect of introduction on two game populations in Czechoslovakia. J. Wildl. Manage. 1951;15(1):113–114.
- 95. Curry BE, Ralls K, Brownell RLJ. Prospects for captive breeding of poorly known small cetacean species. Endangered Species Res. 2013;19(3):223–243.
- 96. Tas’an M, Leatherwood S. Cetaceans live-captured for Jaya Ancol Oceanarium, Djakarta, 1974–1982. Report of the International Whaling Commission. 1984;34:485–489.