Assessment of conservation units for the Sumatran rhinoceros (Dicerorhinus sumatrensis)

An assessment of conservation units for the Sumatran rhinoceros (Dicerorhinus sumatrensis) was conducted using a population aggregation analysis (PAA) of mitochondrial DNA site substitutions. Populations were defined as the three geographically separated regions of West Malaysia, Sumatra, and Borneo. The intent of this assessment was to explore management options for this highly endangered lineage rather than conduct a traditional taxonomic revision. Individual DNA positions were not diagnostic for any population. A single haplotype provided a character as support for diagnosing the West Malaysian and Bornean population. The haplotypes on West Malaysia and Sumatra were more similar to each other than either was to the one on Borneo. These data, and a review of the morphological characters, support the option of treating Sumatran rhinos as a single conservation unit, providing managers with greater flexibility in managing the unique Dicerorhine lineage. © 1995 Wiley-Liss, Inc.     

Molecular and morphological divergence in the butterfly genus Lycaeides (Lepidoptera: Lycaenidae) in North America: evidence of recent speciation

Male genital morphology, allozyme allele frequencies and mtDNA sequence variation were surveyed in the butterfly species Lycaeides idas and L. melissa from across much of their range in North America. Despite clear differences in male genital morphology, wing colour patterns and habitat characteristics, genetic variation was not taxonomically or geographically structured and the species were not identifiable by either genetic data set. Genetic distances (Nei's D=0.002–0.078, calculated from allozyme data) between all populations of both species were within the range commonly observed for conspecific populations of other butterflies. The most frequent mtDNA haplotype was present in individuals of both species in populations from southern California to Wisconsin. We conclude that speciation has probably happened recently and the lack of genetic differentiation between the species is the product of either (1) recent or ongoing gene flow at neutral loci, and/or (2) an insufficiency of time for lineage sorting. The evolution of male genital morphology, wing colour patterns and ecological characteristics has proceeded more rapidly than allozyme or mtDNA evolution.     

Evaluating the utility of cox1 for cetacean species identification

The mitochondrial cytochrome c oxidase I (cox1) gene has been promoted as a universal reference gene, or barcode, to identify organisms to the species level. We evaluated whether cox1 would be appropriate to diagnose cetacean species. The 5′ end of cox1 (686 base pairs, bp) was sequenced for 46 of 86 recognized species of cetaceans. In addition, we included 105 sequences from GenBank, increasing our taxonomic coverage to 61 species. Particular focus was placed on sampling two subfamilies that contain closely related taxa: the Delphininae and the Globicephalinae. Species‐specific sequences were observed for all but three taxa (Delphinus delphis, D. capensis, and Stenella coeruleoalba). Although correct assignment was seen for most species, significant overlap between intra‐ and interspecific variation makes cox1 an imperfect barcode for cetaceans. The efficacy of cox1 was compared to the 5′ end of the cytochrome b (cytb) gene, a mitochondrial region routinely used for cetacean species identification. Although cytb performed better than cox1 for some species, this marker could not differentiate other closely related taxa (Eubalaena spp.). Species identification for taxa not reliably identified using cox1 or cytb might be best addressed through use of multiple mitochondrial DNA fragments or other newly developed markers.     

Assessment of genetic diversity in the critically endangered Australian corroboree frogs, Pseudophryne corroboree and Pseudophryne pengilleyi, identifies four evolutionarily significant units for conservation

The iconic and brightly coloured Australian northern corroboree frog, Pseudophryne pengilleyi , and the southern corroboree frog, Pseudophryne corroboree are critically endangered and may be extinct in the wild within 3 years. We have assembled samples that cover the current range of both species and applied hypervariable microsatellite markers and mitochondrial DNA sequences to assess the levels and patterns of genetic variation. The four loci used in the study were highly variable, the total number of alleles observed ranged from 13 to 30 and the average number of alleles per locus was 19. Expected heterozygosity of the four microsatellite loci across all populations was high and varied between 0.830 and 0.935. Bayesian clustering analyses in structure strongly supported four genetically distinct populations, which correspond exactly to the four main allopatric geographical regions in which the frogs are currently found. Individual analyses performed on the separate regions showed that breeding sites within these four regions could not be separated into distinct populations. Twelve mtND2 haplotypes were identified from 66 individuals from throughout the four geographical regions. A statistical parsimony network of mtDNA haplotypes shows two distinct groups, which correspond to the two species of corroboree frog, but with most of the haplotype diversity distributed in P. pengilleyi . These results demonstrate an unexpectedly high level of genetic diversity in both species. Our data have important implications for how the genetic diversity is managed in the future. The four evolutionarily significant units must be protected and maintained in captive breeding programmes for as long as it is possible to do.     

Phylogeny,phylogeography, and systematics of the American pea crab genus Calyptraeotheres Campos, 1990, inferred from molecular markers

We used mitochondrial cytochrome oxidase I (COI) and the large ribosomal subunit (16S) genes to establish evolutionary relationships amongst species of Calyptraeotheres, evaluate their usefulness as DNA‐barcoding genes, and assess molecular diversity at the population level within Calyptraeotheres garthi. Bayesian, maximum likelihood, and maximum parsimony phylogenies confirmed the monophyly of Calyptraeotheres, showing that the ancestor of C. garthi, Calyptraeotheres hernandezi, and Calyptraeotheres granti radiated after the formation of the Panamanian isthmus. This finding contradicts the austral/tropical hypothesis previously proposed based on morphological data. The COI and 16S distance matrices supported separation of species as well as the genera, and corroborated that DNA barcoding is a useful tool and complements the classical taxonomy in Pinnotheridae. Phylogenetic and genetic distance analyses suggested that C. hernandezi is a junior synonym of C. garthi. Finally, C. garthi did not show a population structure across its distribution range, and showed a pattern consistent with a recent population expansion event that began 230–300 Kya. © 2013 The Linnean Society of London     

Phylogeography and conservation genetics of the Iowa pleistocene snail

The Iowa Pleistocene snail, Discus macclintocki, is an endangered species that survives only in relictual populations on algific (cold-air) talus slopes in northeast Iowa and northwest Illinois in the central region of the USA. These populations are believed to have been isolated since the temperatures began to warm at the end of the last glacial period around 16 500 years ago. DNA sequencing of the 16s rRNA gene of the mitochondria was used to determine the genetic relationship among 10 populations and the genetic diversity within these populations. Genetic diversity is extremely high within this species with 40 haplotypes spread across the 10 populations sampled within a 4000 km2 region. Phylogenetic analyses showed that haplotypes formed monophyletic groups by the watershed on which they were found, suggesting that watersheds were important historical avenues of gene flow. Genetic distances were strongly related to the geographical distance among all populations, but this relationship was dependent on the scale being considered.     

A survey of equid mitochondrial DNA: Implications for the evolution, genetic diversity and conservation of Equus

The evolution, taxonomy and conservation of the genus Equuswere investigated by examining the mitochondrial DNA sequences of thecontrol region and 12S rRNA gene. The phylogenetic analysis of thesesequences provides further evidence that the deepest node in thephylogeny of the extant species is a divergence between twolineages; one leading to the ancestor of modern horses (E.ferus, domestic and przewalskii) and the other to thezebra and ass ancestor, with the later speciation events of the zebrasand asses occurring either as one or more rapid radiations, or withextensive secondary contact after speciation. Examination of the geneticdiversity within species suggested that two of the E. hemionussubspecies (E. h. onager and E. h. kulan) onlyrecently diverged, and perhaps, are insufficiently different to beclassified as separate subspecies. The genetic divergence betweendomestic and wild forms of E. ferus (horse) and E.africanus (African ass) was no greater than expected within anequid species. In E. burchelli (plains zebra) there was anindication of mtDNA divergence between populations increasing withdistance. The implications of these results for equid conservation arediscussed and recommendations are made for conservation action.     

The effect of canopy cover and seasonal change on host plant quality for the endangered Karner blue butterfly (Lycaeidesmelissasamuelis)

Larvae of the Karner blue butterfly, Lycaeidesmelissasamuelis, feed solely on wild lupine, Lupinusperennis, from the emergence to summer senescence of the plant. Wild lupine is most abundant in open areas but Karner blue females oviposit more frequently on lupines growing in moderate shade. Can differences in lupine quality between open and shaded areas help explain this disparity in resource use? Furthermore, many lupines are senescent before the second larval brood completes development. How does lupine senescence affect larval growth? We addressed these questions by measuring growth rates of larvae fed lupines of different phenological stages and lupines growing under different shade conditions. The habitat conditions under which lupines grew and plant phenological stage did not generally affect final larval or pupal weight but did significantly affect duration of the larval period. Duration was shortest for larvae fed leaves from flowering lupines and was negatively correlated with leaf nitrogen concentration. Ovipositing in areas of moderate shade should increase?second-brood larval exposure to flowering lupines. In addition, larval growth was significantly faster on shade-grown lupines that were in seed than on similar sun-grown lupines. These are possible advantages of the higher-than-expected oviposition rate on shade-grown lupines. Given the canopy-related trade-off between lupine?abundance and quality, maintenance of canopy heterogeneity is an important conservation management goal. Larvae were also fed leaves growing in poor soil conditions and leaves with mildew infection. These and other feeding treatments that we anticipated would inhibit larval growth often did not. In particular, ant-tended larvae exhibited the highest weight gain per amount of lupine eaten and a relatively fast growth rate. This represents an advantage of ant tending to Karner blue larvae.     

On the comparison of population-level estimates of haplotype and nucleotide diversity: a case study using the gene cox1 in animals

Estimates of genetic diversity represent a valuable resource for biodiversity assessments and are increasingly used to guide conservation and management programs. The most commonly reported estimates of DNA sequence diversity in animal populations are haplotype diversity (h) and nucleotide diversity (π) for the mitochondrial gene cytochrome c oxidase subunit I (cox1). However, several issues relevant to the comparison of h and π within and between studies remain to be assessed. We used population-level cox1 data from peer-reviewed publications to quantify the extent to which data sets can be re-assembled, to provide a standardized summary of h and π estimates, to explore the relationship between these metrics and to assess their sensitivity to under-sampling. Only 19 out of 42 selected publications had archived data that could be unambiguously re-assembled; this comprised 127 population-level data sets (n ≥ 15) from 23 animal species. Estimates of h and π were calculated using a 456-base region of cox1 that was common to all the data sets (median h=0.70130, median π=0.00356). Non-linear regression methods and Bayesian information criterion analysis revealed that the most parsimonious model describing the relationship between the estimates of h and π was π=0.0081 h(2). Deviations from this model can be used to detect outliers due to biological processes or methodological issues. Subsampling analyses indicated that samples of n>5 were sufficient to discriminate extremes of high from low population-level cox1 diversity, but samples of n ≥ 25 are recommended for greater accuracy.     

Identifying evolutionarily significant units for conservation of the endangered Malus sieversii using genome‐wide RADseq data

Malus sieversii, a wild progenitor of the domesticated apple, is an endangered species and is assigned second conservation priority by the China Plant Red Data Book. It is urgent to carry out in situ conservation of this species, but previous studies have not identified evolutionarily significant units (ESUs) for conservation management. In this study, we investigated the genetic diversity and relationships of six M. sieversii populations from China using integrated analysis of microsatellite (nSSR) data, genome‐wide SNPs and previous results in order to propose a reasonable conservation management. The results showed that levels of genetic diversity were inconsistently reflected by our nSSR and previous studies, suggesting that indices of genetic diversity are not effective to identify priority conservation areas for M. sieversii. Based on the selection criteria of ESUs for endangered species conservation, ESUs should reflect lineage divergence, geographical separation and different adaptive variation. Our phylogenetic tree based on genome‐wide SNPs yielded a clear relationship of divergent lineages among M. sieversii populations, leading to new different from those of previous studies. Three independent lineages, including the pairs of populations Huocheng‐Yining, Gongliu‐Xinyuan and Tuoli‐Emin, were identified. The geographic distances between populations among the different phylogenetic lineages were much greater than those within the same phylogenetic lineage. A cluster analysis on environmental variables showed that the three independent lineages inhabit different environmental conditions, suggesting that they may have adapted to different environments. Based on the results, we propose that three independent ESUs should be recognized as conservation units for M. sieversii in China.     

Assessment of the stocked or wild origin of anadromous brown trout (Salmo trutta L.) in a Danish river system, using mitochondrial DNA RFLP analysis

Declines in the number of anadromous brown trout in the Karup River in Denmark, due to environmental degradation, led to the stocking of large numbers of hatchery trout during the 1980s. This practice was gradually replaced by stocking with the offspring of electrofished local trout The genetic contribution of the hatchery fish to the current population of anadromous trout in the river was estimated by restriction fragment length polymorphism analysis of mitochondrial DNA, using seven restriction endonucleases. Fish from the hatchery strain as well as from five locations in the river system, and from a further unstocked river were screened. Eight haplotypes were observed. The distribution and frequencies of the observed haplotypes revealed little genetic differentiation among stocked populations. The hatchery strain differed significantly from the stocked populations. One haplotype which was found at a high frequency in the hatchery strain was almost absent from the stocked populations. This suggests that the genetic contribution of the hatchery trout to the current population is much less than would be expected from the number of stocked fish. The possible reasons for the failure of the hatchery trout to contribute to the gene pool, and also the implications for conservation biology, are discussed.     

High genetic diversity in a small population: the case of Chilean blue whales

It is generally assumed that species with low population sizes have lower genetic diversities than larger populations and vice versa. However, this would not be the case for long‐lived species with long generation times, and which populations have declined due to anthropogenic effects, such as the blue whale (Balaenoptera musculus). This species was intensively decimated globally to near extinction during the 20th century. Along the Chilean coast, it is estimated that at least 4288 blue whales were hunted from an apparently pre‐exploitation population size (k) of a maximum of 6200 individuals (Southeastern Pacific). Thus, here, we describe the mtDNA (control region) and nDNA (microsatellites) diversities of the Chilean blue whale aggregation site in order to verify the expectation of low genetic diversity in small populations. We then compare our findings with other blue whale aggregations in the Southern Hemisphere. Interestingly, although the estimated population size is small compared with the pre‐whaling era, there is still considerable genetic diversity, even after the population crash, both in mitochondrial (N = 46) and nuclear (N = 52) markers (H_d = 0.890 and H_o = 0.692, respectively). Our results suggest that this diversity could be a consequence of the long generation times and the relatively short period of time elapsed since the end of whaling, which has been observed in other heavily‐exploited whale populations. The genetic variability of blue whales on their southern Chile feeding grounds was similar to that found in other Southern Hemisphere blue whale feeding grounds. Our phylogenetic analysis of mtDNA haplotypes does not show extensive differentiation of populations among Southern Hemisphere blue whale feeding grounds. The present study suggests that although levels of genetic diversity are frequently used as estimators of population health, these parameters depend on the biology of the species and should be taken into account in a monitoring framework study to obtain a more complete picture of the conservation status of a population.     

The complete mitochondrial genome of Rhynchocypris oxycephalus (Teleostei: Cyprinidae) and its phylogenetic implications

Rhynchocypris oxycephalus (Teleostei: Cyprinidae) is a typical small cold water fish, which is distributed widely and mainly inhabits in East Asia. Here, we sequenced and determined the complete mitochondrial genome of R. oxycephalus and studied its phylogenetic implication. R. oxycephalus mitogenome is 16,609 bp in length (GenBank accession no.: MH885043), and it contains 13 protein‐coding genes (PCGs), two rRNA genes, 22 tRNA genes, and two noncoding regions (the control region and the putative origin of light‐strand replication). 12 PCGs started with ATG, while COI used GTG as the start codon. The secondary structure of tRNA‐Ser (AGN) lacks the dihydrouracil (DHU) arm. The control region is 943bp in length, with a termination‐associated sequence, six conserved sequence blocks (CSB‐1, CSB‐2, CSB‐3, CSB‐D, CSB‐E, CSB‐F), and a repetitive sequence. Phylogenetic analysis was performed with maximum likelihood and Bayesian methods based on the concatenated nucleotide sequence of 13 PCGs and the complete sequence without control region, and the result revealed that the relationship between R. oxycephalus and R. percnurus is closest, while the relationship with R. kumgangensis is farthest. The genus Rhynchocypris is revealed as a polyphyletic group, and R. kumgangensis had distant relationship with other Rhynchocypris species. In addition, COI and ND2 genes are considered as the fittest DNA barcoding gene in genus Rhynchocypris. This work provides additional molecular information for studying R. oxycephalus conservation genetics and evolutionary relationships.     

Identification of molecular markers for DNA barcoding in the Aphidiinae (Hym. Braconidae)

Reliable identification of Aphidiinae species (Braconidae) is a prerequisite for conducting studies on aphid–parasitoid interactions at the community level. However, morphological identification of Aphidiinae species remains problematic even for specialists and is almost impossible with larval stages. Here, we compared the efficiency of two molecular markers [mitochondrial cytochrome c oxydase I (COI) and nuclear long wavelength rhodopsin (LWRh)] that could be used to accurately identify about 50 species of Aphidiinae that commonly occur in aphid–parasitoid networks in northwestern Europe. We first identified species on a morphological basis and then assessed the consistency of genetic and morphological data. Probably because of mitochondrial introgression, Aphidius ervi and A. microlophii were indistinguishable on the basis of their COI sequences, whereas LWRh sequences discriminated these species. Conversely, because of its lower variability, LWRh failed to discriminate two pairs of species (Aphidius aquilus, Aphidius salicis, Lysiphlebus confusus and Lysiphlebus fabarum). Our study showed that no unique locus but a combination of two genes should be used to accurately identify members of Aphidiinae.     

The genetics of amphibian declines: population substructure and molecular differentiation in the yosemite toad, Bufo canorus (Anura, bufonidae) based on single-strand conformation polymorphism analysis (SSCP) and mitochondrial DNA sequence data

We present a comprehensive survey of genetic variation across the range of the narrowly distributed endemic Yosemite toad Bufo canorus, a declining amphibian restricted to the Sierra Nevada of California. Based on 322 bp of mitochondrial cytochrome b sequence data, we found limited support for the monophyly of B. canorus and its closely related congener B. exsul to the exclusion of the widespread western toad B. boreas. However, B. exsul was always phylogenetically nested within B. canorus, suggesting that the latter may not be monophyletic. SSCP (single-strand conformation polymorphism) analysis of 372 individual B. canorus from 28 localities in Yosemite and Kings Canyon National Parks revealed no shared haplotypes among these two regions and lead us to interpret these two parks as distinct management units for B. canorus. Within Yosemite, we found significant genetic substructure both at the level of major drainages and among breeding ponds. Kings Canyon samples show a different pattern, with substantial variation among breeding sites, but no substructure among drainages. Across the range of B. canorus as well as among Yosemite ponds, we found an isolation-by-distance pattern suggestive of a stepping stone model of migration. However, in Kings Canyon we found no hint of such a pattern, suggesting that movement patterns of toads may be quite different in these nearby parklands. Our data imply that management for B. canorus should focus at the individual pond level, and effective management may necessitate reintroductions if local extirpations occur. A brief review of other pond-breeding anurans suggests that highly structured populations are often the case, and thus that our results for B. canorus may be general for other species of frogs and toads.     

Islands within an island: phylogeography and conservation genetics of the endangered Hawaiian tree snail Achatinella mustelina

Mitochondrial DNA (mtDNA) sequences were used to evaluate phylogeographic structure within and among populations of three endangered Hawaiian tree snail species (n = 86). The primary focus of this investigation was on setting conservation priorities for Achatinella mustelina. Limited data sets for two additional endangered Hawaiian tree snails, A. livida and A. sowerbyana, were also developed for comparative purposes. Pairwise genetic distance matrices and phylogenetic trees were generated, and an analysis of molecular variance was performed on 675-base pair cytochrome oxidase I gene sequences from multiple populations of Hawaiian tree snails. Sequence data were analysed under distance-based maximum-likelihood, and maximum-parsimony optimality criteria. Within the focal species, A. mustelina, numbers of variable and parsimony informative sites were 90 and 69, respectively. Pairwise intraspecific mtDNA sequence divergence ranged from 0 to 5.3% in A. mustelina, from 0 to 1.0% in A. livida and from 0 to 1.9% in A. sowerbyana. For A. mustelina, population genetic structure and mountain topography were strongly correlated. Maximum genetic distances were observed across deep, largely deforested valleys, and steep mountain peaks, independent of geographical distance. However, in certain areas where forest cover is presently fragmented, little mtDNA sequence divergence exists despite large geographical scales (8 km). Genetic data were used to define evolutionarily significant units for conservation purposes including decisions regarding placement of predator exclusion fences, captive propagation, re-introduction and translocation.