The subspecies concept in butterflies: has its application in taxonomy and conservation biology outlived its usefulness?

Subspecies lie at the interface between systematics and population genetics, and represent a unit of biological organization in zoology that is widely used in the disciplines of taxonomy and conservation biology. In this review, we explore the utility of subspecies in relation to their application in systematics and biodiversity conservation, and briefly summarize species concepts and criteria for their diagnosis, particularly from an invertebrate perspective. The subspecies concept was originally conceived as a formal means of documenting geographical variation within species based on morphological characters; however, the utility of subspecies is hampered by inconsistencies by which they are defined conceptually, a lack of objective criteria or properties that serve to delimit their boundaries, and their frequent failure to reflect distinct evolutionary units according to population genetic structure. Moreover, the concept has been applied to populations largely comprising different components of genetic diversity reflecting contrasting evolutionary processes. We recommend that, under the general lineage (unified) species concept, the definition of subspecies be restricted to extant animal groups that comprise evolving populations representing partially isolated lineages of a species that are allopatric, phenotypically distinct, and have at least one fixed diagnosable character state, and that these character differences are (or are assumed to be) correlated with evolutionary independence according to population genetic structure. Phenotypic character types include colour pattern, morphology, and behaviour or ecology. Under these criteria, allopatric subspecies are a type of evolutionarily significant unit within species in that they show both neutral divergence through the effects of genetic drift and adaptive divergence under natural selection, and provide an historical context for identifying biodiversity units for conservation. Conservation of the adaptedness and adaptability of gene pools, however, may require additional approaches. Recent studies of Australian butterflies exemplify these points. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

Genetic Analysis of the <Emphasis Type="Italic">Indri</Emphasis> Reveals No Evidence of Distinct Subspecies

Subspecies were traditionally defined by identifying gaps between phenotypes across the geographic range of a species, and may represent important units in the development of conservation strategies focused on preserving genetic diversity. Previous taxonomic research proposed that phenotypic variation between scattered Indri indri populations warranted the naming of two distinct subspecies, I. i. indri and I. i. variegatus. We tested these subspecific designations using mitochondrial sequence data generated from the control region or D-loop (569 bp) and a large section (2362 bp) of multiple genes and tRNAs known as Pastorini’s fragment and nuclear microsatellite markers. This study used 114 samples of I. indri from 12 rainforest sites in eastern Madagascar, encompassing the entire range of the species. These genetic samples represent multiple populations from low- and high-elevation forests from both putative subspecies. Molecular analyses of the mitochondrial sequence data did not support the two proposed subspecies. Furthermore, the microsatellite analyses showed no significant differences across the range beyond population level differentiation. This study demonstrates the utility of incorporating multiple lines of evidence in addition to phenotypic traits to define species or subspecies.     

Investigating Concordance among Genetic Data, Subspecies Circumscriptions and Hostplant Use in the Nymphalid Butterfly Polygonia faunus

Subspecies are commonly used taxonomic units to formally describe intraspecific geographic variation in morphological traits. However, the concept of subspecies is not clearly defined, and there is little agreement about what they represent in terms of evolutionary units, and whether they can be used as reliably useful units in conservation, evolutionary theory and taxonomy. We here investigate whether the morphologically well-characterized subspecies in the North American butterfly Polygonia faunus are supported by genetic data from mitochondrial sequences and eight microsatellite loci. We also investigate the phylogeographic structure of P. faunus and test whether similarities in host-plant use among populations are related to genetic similarity. Neither the nuclear nor the mitochondrial data corroborated subspecies groupings. We found three well defined genetic clusters corresponding to California, Arizona and (New Mexico+Colorado). There was little structuring among the remaining populations, probably due to gene flow across populations. We found no support for the hypothesis that similarities in host use are related to genetic proximity. The results indicate that the species underwent a recent rapid expansion, probably from two glacial refugia in western North America. The mitochondrial haplotype network indicates at least two independent expansion phases into eastern North America. Our results clearly demonstrate that subspecies in P. faunus do not conform to the structuring of genetic variation. More studies on insects and other invertebrates are needed to better understand the scope of this phenomenon. The results of this study will be crucial in designing further experiments to understand the evolution of hostplant utilization in this species.     

Subspecies of the Central American Squirrel Monkey (Saimiri oerstedii) as Units for Conservation

The accurate diagnosis of conservation units now typically includes recognition of genetic diversity and unique evolutionary lineages and is necessary to inform the conservation management of endangered species. We evaluated whether the two currently recognized subspecies of the endangered Central American squirrel monkey (Saimiri oerstedii) in Costa Rica are evolutionarily significant units (ESUs) that should be managed separately in conservation efforts. We used previously published sequences of 50 individuals of Saimiri oerstedii for 880 bp of the mtDNA d-loop and genotypes of 244 individuals for 16 microsatellites and conducted novel analyses to characterize genetic differentiation between subspecies of Saimiri oerstedii. We measured sequence differentiation and inferred an intraspecific molecular phylogeny and a haplotype network, and found consistent results supporting statistically significant divergence and reciprocal monophyly between subspecies. A population aggregation analysis also supported Saimiri oerstedii citrinellus and S. o. oerstedii as diagnosably distinct units. These results confirm previous genetic studies with smaller sample sizes and are consistent with other factors including differences in pelage and morphology and divergence at nuclear markers. Conservation managers should manage these subspecies separately to prevent the loss of genetic diversity via artificially induced outbreeding. High levels of genetic diversity may buffer populations against outside extinction pressures, to which Saimiri oerstedii are vulnerable because of their dwindling habitat and small population size.     

Conservation genetics of the yellow-bellied toad (Bombina variegata): population structure,genetic diversity and landscape effects in an endangered amphibian

Revealing patterns of genetic diversity and barriers for gene flow are key points for successful conservation in endangered species. Methods based on molecular markers are also often used to delineate conservation units such as evolutionary significant units and management units. Here we combine phylo-geographic analyses (based on mtDNA) with population and landscape genetic analyses (based on microsatellites) for the endangered yellow-bellied toad Bombina variegata over a wide distribution range in Germany. Our analyses show that two genetic clusters are present in the study area, a northern and a southern/central one, but that these clusters are not deeply divergent. The genetic data suggest high fragmentation among toad occurrences and consequently low genetic diversity. Genetic diversity and genetic connectivity showed a negative relationship with road densities and urban areas surrounding toad occurrences, indicating that these landscape features act as barriers to gene flow. To preserve a maximum of genetic diversity, we recommend considering both genetic clusters as management units, and to increase gene flow among toad occurrences with the aim of restoring and protecting functional meta-populations within each of the clusters. Several isolated populations with especially low genetic diversity and signs of inbreeding need particular short-term conservation attention to avoid extinction. We also recommend to allow natural gene flow between both clusters but not to use individuals from one cluster for translocation or reintroduction into the other. Our results underscore the utility of molecular tools for species conservation, highlight outcomes of habitat fragmentation onto the genetic structure of an endangered amphibian and reveal particularly threatened populations in need for urgent conservation efforts.

     

Limited structure and widespread diversity suggest potential buffers to genetic erosion in a threatened grassland shrub Pimelea spinescens (Thymelaeaceae)

Pimelea spinescens is a critically endangered species of the temperate grasslands of southeastern Australia. Two subspecies are recognised. Subspecies, P. spinescens subsp. spinescens, formerly common and widespread, is found in isolated remnants of previously extensive grasslands. The second subspecies, subsp. pubiflora, is thought to have been historically rare with only two geographically-isolated extant populations. The grassland communities exist now as fragmented remnants representing <1 % of their extent prior to European settlement in the early 1800s. Conservation management strategies for species in these critically endangered ecosystems rely on an understanding of genetic diversity and population structure to ensure long-term evolutionary potential. We used chloroplast DNA (cpDNA) and microsatellite markers to examine the population genetic structure of both subspecies. Analysis of cpDNA revealed 14 haplotypes with high divergence between the single haplotype found in subsp. pubiflora and most remaining haplotypes restricted to subsp. spinescens. Microsatellites also indicated high genetic differentiation between subspecies but little evidence of sub-structuring within either subspecies. Results suggest that seed dispersal has not been as limited as previously thought. In this fragmented habitat, a lack of genetic structure suggests buffers to genetic erosion, with current patterns reflecting genetic diversity prior to fragmentation. Plant longevity and the presence of seed banks may contribute to the maintenance of these patterns, resulting in a lag between fragmentation and genetic erosion. Whilst factors such as longevity and seed banks may be preserving historic genetic diversity, management is required to ensure the maintenance of this diversity into the future.     

Unraveling evolutionary lineages in the limbless fossorial skink genus Acontias (Sauria: Scincidae): are subspecies equivalent systematic units?

Subspecies in the limbless, endemic African fossorial skink genus Acontias constitute ill-defined operational taxonomic units, consequently considerable systematic debate has lingered on the systematic diversity within Acontias. In the present study, the systematic affinities among acontine taxa are explored with the utility of partial sequence data from two mitochondrial gene loci (16S rRNA and cytochrome oxidase subunit 1 (COI)) for all taxa, while two additional loci (12S rRNA, cytochrome b) were used to investigate relationships within the Acontias meleagris complex. Phylogenetic results, derived from the combined analysis, revealed two monophyletic clades. Clade 1 is comprised of small-bodied skinks while clade 2 comprised the medium bodied skinks. Within clade 2 none of the traditionally recognized subspecies formed reciprocally monophyletic groups. Furthermore, constraining the topology and enforcing sister taxa relationships between the assumed subspecies, consistently recovered a topology that was statistically significant worse, indicating that the traditionally designated subspecies groupings probably represent invalid taxonomic units, thus clearly reflecting considerable discord with current taxonomy. The burrowing life style of these lizards has probably led to marked convergent evolution and constrained the development of diagnostic morphological characters among these species. Morphological similarities in color as well as scale architecture within Acontias are labile and highly homoplaseous and do not reflect the evolutionary history of the group. Taxonomic implications of these results are discussed.     

A role for molecular genetics in the recognition and conservation of endangered species

Taxonomies based on morphological traits alone sometimes provide inadequate or misleading guides to phylogenetic distinctions at the subspecies and species levels. Yet taxonomic assignments inevitably shape perceptions of biotic diversity, including recognition of endangered species. Case histories are discussed in which the data of molecular genetics revealed prior systematic errors of the two possible kinds: taxonomic recognition of groups showing little evolutionary differentiation, and lack of taxonomic recognition of phylogenetically distinct forms. In such cases, conservation efforts for 'endangered species' can be misdirected with respect to the goal of protecting biological diversity.     

Application of molecular genetics and geometric morphometrics to taxonomy and conservation of cave beetles in central Italy

Integration of molecular genetic techniques and geometric morphometrics represent a valuable tool in the resolution of taxonomic uncertainty and the identification of significant units for conservation. We combined mitochondrial DNA cytochrome c oxidase subunit II gene sequence data and geometric morphometric analysis to examine taxonomic status and identify units for conservation in four species of the hypogean beetle Duvalius (Coleoptera, Trechinae) using mainly museum specimens collected in central Italy. Previous taxonomic studies based on morphological traits described several subspecies often inhabiting geographically distinct caves. Phylogenetic analysis identified two well supported monophyletic lineages and a number of different clades with relatively small genetic differences, suggesting a short divergence time in line with known geological history of the study area. Geometric morphometrics, on the other hand, recovered a high level of distinctiveness among specimens. Both genetic and morphometric analyses did not entirely corroborate former taxonomic nomenclature, suggesting possible rearrangements and the definition of evolutionary significant units. Beetles of the genus Duvalius are protected by regional laws and the majority of taxa considered in this study inhabit caves located outside protected areas. Our study advocates the importance of devoting protection efforts to networks of cave ecosystems rather than single locations or species.     

Taxonomy and conservation of Vietnam's primates: a review

Vietnam has the highest number of primate taxa overall (24-27) and the highest number of globally threatened primate taxa (minimum 20) in Mainland Southeast Asia. Conservation management of these species depends in part on resolving taxonomic uncertainties, which remain numerous among the Asian primates. Recent research on genetic, morphological, and acoustic diversity in Vietnam's primates has clarified some of these uncertainties, although a number of significant classification issues still remain. Herein, we summarize and compare the major current taxonomic classifications of Vietnam's primates, discuss recent advances in the context of these taxonomies, and suggest key areas for additional research to best inform conservation efforts in a region crucial to global primate diversity. Among the most important next steps for the conservation of Vietnam's primates is a new consensus list of Asian primates that resolves current differences between major taxonomies, incorporates recent research advances, and recognizes units of diversity at scales below the species-level, whether termed populations, morphs, or subspecies. Priority should be placed on recognizing distinct populations, regardless of the species concept in use, in order to foster the evolutionary processes necessary for primate populations to cope with inevitable environmental changes. The long-term conservation of Vietnam's primates depends not only on an accepted and accurate taxonomy but also on funding for on-the-ground conservation activities, including training, and the continued dedication and leadership of Vietnamese researchers and managers.     

Drawbacks with the use of microsatellites in phylogeography: the song sparrow Melospiza melodia as a case study

Microsatellites are commonly used molecular markers in phylogeography, and many view them as superior to mitochondrial DNA (mtDNA) gene trees. Being based on frequencies of alleles, and not gene trees, microsatellites exhibit the same analytical drawbacks that resulted in the abandonment of allozymes in genetic studies of population history. I illustrate some these familiar drawbacks by reanalyzing microsatellite data on the song sparrow. Subspecies were previously evaluated with hierarchical analyses of molecular variance, suggesting that subspecies explain 8% of the total variance in microsatellite frequencies. However, this useful heuristic technique only evaluates a priori groupings, and the objective of the study ought to be to discover such groupings, not assume them. In fact, other arbitrary groupings of samples explained the same or greater amounts of variance, and I suggest that for testing subspecies limits, a gene tree is preferable. Grouping population samples by subspecies in the San Francisco Bay area accounts for 1.2% of the microsatellite variation, and despite claims that this informs conservation planning, the data do not support any particular population or subspecies as being genetically or evolutionarily significant. A distance phenogram was used to infer a sequential colonization of the Aleutian Islands, but because individuals were pooled into a priori groups and the phenogram was arbitrarily rooted, this conclusion is tenuous. A plot of heterozygosity vs number of alleles per sample showed that an equally parsimonious interpretation is that current genetic diversity tracks effective population size. Microsatellites should be replaced in nuclear-gene phylogeography by analyses of sequences, which will benefit the study of phylogeography, comparison of nuclear and mtDNA results, and aid in interpreting the results in a conservation context.     

Fuzzy Boundaries: Color and Gene Flow Patterns among Parapatric Lineages of the Western Shovel-Nosed Snake and Taxonomic Implication

Accurate delineation of lineage diversity is increasingly important, as species distributions are becoming more reduced and threatened. During the last century, the subspecies category was often used to denote phenotypic variation within a species range and to provide a framework for understanding lineage differentiation, often considered incipient speciation. While this category has largely fallen into disuse, previously recognized subspecies often serve as important units for conservation policy and management when other information is lacking. In this study, we evaluated phenotypic subspecies hypotheses within shovel-nosed snakes on the basis of genetic data and considered how evolutionary processes such as gene flow influenced possible incongruence between phenotypic and genetic patterns. We used both traditional phylogenetic and Bayesian clustering analyses to infer range-wide genetic structure and spatially explicit analyses to detect possible boundary locations of lineage contact. Multilocus analyses supported three historically isolated groups with low to moderate levels of contemporary gene exchange. Genetic data did not support phenotypic subspecies as exclusive groups, and we detected patterns of discordance in areas where three subspecies are presumed to be in contact. Based on genetic and phenotypic evidence, we suggested that species-level diversity is underestimated in this group and we proposed that two species be recognized, Chionactis occipitalis and C. annulata. In addition, we recommend retention of two subspecific designations within C. annulata (C. a. annulata and C. a. klauberi) that reflect regional shifts in both genetic and phenotypic variation within the species. Our results highlight the difficultly in validating taxonomic boundaries within lineages that are evolving under a time-dependent, continuous process.     

Ancestral polymorphisms in genetic markers obscure detection of evolutionarily distinct populations in the endangered Florida grasshopper sparrow (Ammodramus savannarum floridanus)

Genetic analyses of bird subspecies designated as conservation units can address whether they represent units with independent evolutionary histories and provide insights into the evolutionary processes that determine the degree to which they are genetically distinct. Here we use mitochondrial DNA control region sequence and six microsatellite DNA loci to examine phylogeographical structure and genetic differentiation among five North American grasshopper sparrow (Ammodramus savannarum) populations representing three subspecies, including a population of the endangered Florida subspecies (A. s. floridanus). This federally listed taxon is of particular interest because it differs phenotypically from other subspecies in plumage and behaviour and has also undergone a drastic decline in population size over the past century. Despite this designation, we observed no phylogeographical structure among populations in either marker: mtDNA haplotypes and microsatellite genotypes from floridanus samples did not form clades that were phylogenetically distinct from variants found in other subspecies. However, there was low but significant differentiation between Florida and all other populations combined in both mtDNA (FST = 0.069) and in one measure of microsatellite differentiation (theta = 0.016), while the non-Florida populations were not different from each other. Based on analyses of mtDNA variation using a coalescent-based model, the effective sizes of these populations are large (approximately 80,000 females) and they have only recently diverged from each other (< 26,000 ybp). These populations are probably far from genetic equilibrium and therefore the lack of phylogenetic distinctiveness of the floridanus subspecies and minimal genetic differentiation is due most probably to retained ancestral polymorphism. Finally, levels of variation in Florida were similar to other populations supporting the idea that the drastic reduction in population size which has occurred within the last 100 years has not yet had an impact on levels of variation in floridanus. We argue that despite the lack of phylogenetic distinctiveness of floridanus genotypes the observed genetic differentiation and previously documented phenotypic differences justify continued designation of this subspecies as a protected population segment.     

Phylogeography and the conservation of coral reef fishes

Here we present a review of how the study of the geographic distribution of genetic lineages (phylogeography) has helped identify management units, evolutionary significant units, cryptic species, and areas of endemism, and how this information can help efforts to achieve effective conservation of coral reefs. These studies have confirmed the major biogeographic barriers that were originally identified by tropical species distributions. Ancient separations, identified primarily with mtDNA sequence comparisons, became apparent between populations on each side of the barriers. The general lack of correlation between pelagic larval duration and genetic connectivity across barriers indicates that life history and ecology can be as influential as oceanography and geography in shaping evolutionary partitions within ocean basins. Hence, conservation strategies require a recognition of ecological hotspots, those areas where habitat heterogeneity promotes speciation, in addition to more traditional approaches based on biogeography. Finally, the emerging field of genomics will add a new dimension to phylogeography, allowing the study of genes that are pertinent to recent and ongoing differentiation, and ultimately providing higher resolution to detect evolutionary significant units that have diverged in an ecological time scale.     

Reconstruction of caribou evolutionary history in Western North America and its implications for conservation

The role of Beringia as a refugium and route for trans-continental exchange of fauna during glacial cycles of the past 2million years are well documented; less apparent is its contribution as a significant reservoir of genetic diversity. Using mitochondrial DNA sequences and 14 microsatellite loci, we investigate the phylogeographic history of caribou (Rangifer tarandus) in western North America. Patterns of genetic diversity reveal two distinct groups of caribou. Caribou classified as a Northern group, of Beringian origin, exhibited greater number and variability in mtDNA haplotypes compared to a Southern group originating from refugia south of glacial ice. Results indicate that subspecies R. t. granti of Alaska and R. t. groenlandicus of northern Canada do not constitute distinguishable units at mtDNA or microsatellites, belying their current status as separate subspecies. Additionally, the Northern Mountain ecotype of woodland caribou (presently R. t. caribou) has closer kinship to caribou classified as granti or groenlandicus. Comparisons of mtDNA and microsatellite data suggest that behavioural and ecological specialization is a more recently derived life history characteristic. Notably, microsatellite differentiation among Southern herds is significantly greater, most likely as a result of human-induced landscape fragmentation and genetic drift due to smaller population sizes. These results not only provide important insight into the evolutionary history of northern species such as caribou, but also are important indicators for managers evaluating conservation measures for this threatened species.     

Historical range contraction,and not taxonomy,explains the contemporary genetic structure of the Australian tree <Emphasis Type="Italic">Acacia dealbata</Emphasis> Link

Irrespective of its causes, strong population genetic structure indicates a lack of gene flow. Understanding the processes that underlie such structure, and the spatial patterns it causes, is valuable for conservation efforts such as restoration. On the other hand, when a species is invasive outside its native range, such information can aid management in the non-native range. Here we explored the genetic characteristics of the Australian tree Acacia dealbata in its native range. Two subspecies of A. dealbata have previously been described based on morphology and environmental requirements, but recent phylogeographic data raised questions regarding the validity of this taxonomic subdivision. The species has been widely planted within and outside its native Australian range and is also a highly successful invasive species in many parts of the world. We employed microsatellite markers to investigate the population genetic diversity and structure among 42 A. dealbata populations from across the species’ native range. We also tested whether environmental variables purportedly relevant for the putative separation of subspecies are linked with population genetic differentiation. We found no relationship between population genetic structure of A. dealbata in Australia and these environmental features. Rather, we identified two geographically distinct genetic clusters that corresponded with populations in the northeastern part of mainland Australia, and the southern mainland and Tasmanian range of the species. Our results do not support the taxonomic subdivision of the species into two distinct subspecies based on environmental features. We therefore assume that the observed morphological differences between the putative subspecies are plastic phenotypic responses. This study provides population genetic information that will be useful for the conservation of the species within Australia as well as to better understand the invasion dynamics of A. dealbata.     

Genetic diversity and population structure in the Barrens Topminnow (<Emphasis Type="Italic">Fundulus julisia</Emphasis>): implications for conservation and management of a critically endangered species

The Barrens Topminnow (Fundulus julisia) has undergone a rapid and dramatic decline. In the 1980s, at least twenty localities with Barrens Topminnows were known to exist in the Barrens Plateau region of middle Tennessee; currently only three areas with natural (not stocked) populations remain. The long-term survival of the Barrens Topminnow will depend entirely on effective management and conservation efforts. Captive propagation and stocking of captive-reared juveniles to suitable habitats have successfully established a handful of self-sustaining populations. However, very little is known about the genetic composition of source and introduced populations including levels of genetic diversity and structuring of genetic variation. Here we use both mitochondrial sequence data and genotypes from 14 microsatellite loci to examine patterns of genetic variation among ten sites, including all sites with natural populations and a subset of sites with introduced (stocked) populations of this species. Mitochondrial sequence analysis reveals extremely low levels of variation within populations and fixed differences between drainages. Microsatellite genotype data shows higher levels of genetic variability and a molecular signature consistent with a recent history of population bottlenecks. Measures of genetic diversity at microsatellite loci including allelic richness are similar within source and introduced populations. Bayesian assignment tests and analysis of molecular variation (AMOVA) support two distinct populations, consistent with drainage boundaries. Results from AMOVA analysis also suggest low levels of genetic connectivity between isolated populations within the same drainage. Here we propose two distinct evolutionary significant units (ESUs) and two management units that reflect this population substructure and warrant consideration in future management efforts.     

Colonization and divergence: phylogeography and population genetics of the Atlantic coast beach mice

Barrier island taxa provide an opportunity to investigate recent evolutionary processes, such as colonization and isolation of recently diverged taxa, and provide important insights into understanding contemporary diversity and the assessment of conservation units. Using rapidly evolving genetic markers (mitochondrial DNA and microsatellites), we studied the Atlantic coast beach mouse subspecies (Peromyscus polionotus decoloratus, P. p. niveiventris, and P. p. phasma). Our data indicate that each of the extant coastal subspecies (P. p. niveiventris and P. p. phasma) is comprised of unique haplotypes indicative of their isolation, while the extinct subspecies, P. p. decoloratus, contain a single haplotype, which was shared with P. p. phasma. Moreover, all the coastal haplotypes originate from a single mainland haplotype found in central Florida, USA. The microsatellite data indicated high levels of genetic structure among our sampled populations. Additionally, these data group the populations into three distinct genetic clusters, with each of the extant coastal subspecies belonging to their own cluster and the mainland individuals forming a separate cluster. The extant Atlantic coast beach mice are on separate evolutionary trajectories, thus representative of separate taxonomic units. Therefore, the data support that two extant subspecies on the Florida Atlantic coast fit the Distinct Population Segment designation and should be managed and conserved as two separate independent units.     

Low genetic diversity and evidence of population structure among subspecies of Nerodia harteri, a threatened water snake endemic to Texas

Nerodia harteri is a threatened small-bodied water snake that occupies one of the most restricted ranges of any snake within the continental United States. It is found closely associated with rivers and tributaries in the Colorado and Brazos river basins, which flow through north-central Texas. Nerodia harteri has been at the center of debate owing to conflicts between conservation efforts and the construction of dams that change or destroy its preferred habitat. Additionally, its taxonomic status has also been under contention with some authors recognizing two subspecies, the Brazos water snake (N. h. harteri) and the Concho water snake (N. h. paucimaculata), whereas other authors consider each separate species. Despite its relatively recent discovery during the 1940s, N. harteri has been the subject of several ecological studies, yet no population genetic assessment of either subspecies has been performed to date. We first evaluated the phylogenetic placement of both subspecies among other North American Nerodia using partial sequence data from the mitochondrial gene cytochrome-b. We then tested for population subdivision among four rivers encompassing the range of N. harteri and tested for the presence of admixture between river basins using mitochondrial sequence data (920?bp of cyt-b) and five cross-species amplified microsatellite loci. We found low mitochondrial haplotype diversity represented by two unique haplotypes in each river basin, which were separated by no more than four nucleotide changes. Nuclear loci showed low genetic diversity and population structuring within and among river basins. We did not find conclusive evidence of admixture between basins, and we support the presence of two separate evolutionarily significant units and two separate management units corresponding to each major river basin. Given increasing natural and anthropogenic threats, we recommend continued ecological and genetic monitoring of both subspecies.     

Conservation Genetics of the Philippine Tarsier: Cryptic Genetic Variation Restructures Conservation Priorities for an Island Archipelago Primate

Establishment of conservation priorities for primates is a particular concern in the island archipelagos of Southeast Asia, where rates of habitat destruction are among the highest in the world. Conservation programs require knowledge of taxonomic diversity to ensure success. The Philippine tarsier is a flagship species that promotes environmental awareness and a thriving ecotourism economy in the Philippines. However, assessment of its conservation status has been impeded by taxonomic uncertainty, a paucity of field studies, and a lack of vouchered specimens and genetic samples available for study in biodiversity repositories. Consequently, conservation priorities are unclear. In this study we use mitochondrial and nuclear DNA to empirically infer geographic partitioning of genetic variation and to identify evolutionarily distinct lineages for conservation action. The distribution of Philippine tarsier genetic diversity is neither congruent with expectations based on biogeographical patterns documented in other Philippine vertebrates, nor does it agree with the most recent Philippine tarsier taxonomic arrangement. We identify three principal evolutionary lineages that do not correspond to the currently recognized subspecies, highlight the discovery of a novel cryptic and range-restricted subcenter of genetic variation in an unanticipated part of the archipelago, and identify additional geographically structured genetic variation that should be the focus of future studies and conservation action. Conservation of this flagship species necessitates establishment of protected areas and targeted conservation programs within the range of each genetically distinct variant of the Philippine tarsier.