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.     

Revival of the genus Tropicoperdix Blyth 1859 (Phasianidae,Aves) using multilocus sequence data

Although Tropicoperdix has been considered to be either a full genus or a species complex within the Phasianid genus Arborophila (hill partridges), there is long‐standing uncertainty regarding the degree of difference that warrants generic separation, including reported anatomical cranial differences. In addition, the intra‐generic taxonomy remains under dispute. Most studies hypothesize that Tropicoperdix comprises three species, while others postulate from one to four species. However, no molecular study has been performed to clarify the systematic and taxonomic uncertainties surrounding Tropicoperdix. In the present study, we performed a series of molecular phylogenetic analyses of Tropicoperdix and Arborophila taxa based on two mitochondrial genes and five nuclear introns. All the results are consistent with the finding that Tropicoperdix and Arborophila are phylogenetically distinct and distant genera, although the precise phylogenetic position of Tropicoperdix remains undetermined. Retrospective examination of external characteristics also supports the generic separation, as well as providing evidence of remarkable multiple character convergence. We propose that Tropicoperdix comprises at least two full species based on mitochondrial data obtained from museum specimens by using a next‐generation sequencing method. © 2015 The Linnean Society of London     

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 relationships of hellbenders in the Ozark highlands of Missouri and conservation implications for the Ozark subspecies (<Emphasis Type="Italic">Cryptobranchus alleganiensis bishopi</Emphasis>)

The hellbender (Cryptobranchus alleganiensis) is an obligately aquatic salamander that is in decline due to habitat loss and disease. Two subspecies of hellbender have been described based on morphological characteristics: C. a. alleganiensis (eastern subspecies) and C. a. bishopi (Ozark hellbender). Current conservation strategies include captive propagation for restorative releases even though information regarding the current levels of genetic variability and structure within populations is not sufficient to effectively plan for conservation of the genetic diversity of the species. To investigate patterns of population structure in the hellbender, we genotyped 276 hellbenders from eight Missouri River drainages, representing both subspecies. Our results showed low levels of within-drainage diversity but strong population structure among rivers, and three distinct genetic clusters. F _ST values ranged from 0.00 to 0.61 and averaged 0.40. Our results confirmed previous reports that C. a. bishopi and C. a. alleganiensis are genetically distinct, but also revealed an equidistant relationship between two groups within C. a. bishopi and all populations of C. a. alleganiensis. Current subspecies delineations do not accurately incorporate genetic structure, and for conservation purposes, these three groups should be considered evolutionarily significant units.     

Genomic data indicate ubiquitous evolutionary distinctiveness among populations of California metalmark butterflies

Conservation geneticists have argued that evolutionarily significant units (ESUs) must be both genetically distinct and adaptively significant to be recognized for conservation protection. High-throughput DNA approaches can greatly increase the power to identify genetic distinctiveness, even if inferring adaptive significance remains a challenge. Here we present the first genomic evaluation of Lange’s metalmark, Apodemia mormo langei (Lepidoptera: Riodinidae), a U.S. federally endangered subspecies restricted to sand dune habitats in a single National Wildlife Refuge in California. Previous work based on very few genetic markers detected little genetic distinction for Lange’s metalmark. We use several thousand genome-wide single nucleotide polymorphisms to characterize the population structure of the A. mormo complex across California and determine if Lange’s metalmark qualifies as an ESU. We found that Lange’s metalmark is genetically identifiable, but is no more distinct than many other isolated populations across the study area. It remains unclear whether this genetic variation is adaptive, and so conservation efforts would benefit from more ecological characterization to determine conservation priorities.     

Molecular systematics of the critically-endangered North American spinymussels (Unionidae: <Emphasis Type="Italic">Elliptio</Emphasis> and <Emphasis Type="Italic">Pleurobema</Emphasis>) and description of <Emphasis Type="Italic">Parvaspina</Emphasis> gen. nov.

Despite being common in numerous marine bivalve lineages, lateral spines are extremely rare among freshwater bivalves (Bivalvia: Unionidae), with only three known species characterized by the presence of spines: Elliptio spinosa, Elliptio steinstansana, and Pleurobema collina. All three taxa are endemic to the Atlantic Slope of southeastern North America, critically endangered, and protected by the US Endangered Species Act. Currently, these species are recognized in two genera and remain a source of considerable taxonomic confusion. Because spines are rare in freshwater mussels and restricted to a small region of North America, we hypothesized that spinymussels represent a monophyletic group. We sequenced two mtDNA gene fragments (COI and ND1) and a fragment of the nuclear ITS-1 locus from >70 specimens. Bayesian and maximum-likelihood phylogenetic reconstructions suggest that the spinymussels do not comprise a monophyletic group. Elliptio steinstansana is sister to P. collina, forming a monophyletic clade that was estimated to have diverged from its most recent ancestor in the late Miocene and is distinct from both Elliptio and Pleurobema; we describe a new genus (Parvaspina gen. nov.) to reflect this relationship. Additionally, E. spinosa forms a monophyletic clade that diverged from members of the core Elliptio lineage in the mid-Pliocene. Furthermore, E. spinosa is genetically divergent from the other spinymussel species, suggesting that spines, while extremely rare in freshwater mussels worldwide, may have evolved independently in two bivalve lineages. Recognizing the genetic distinctiveness and inter-generic relationships of the spinymussels is an important first step towards effectively managing these imperiled species and lays the groundwork for future conservation genetics studies.     

Biogeography and designatable units of <Emphasis Type="Italic">Bombus occidentalis</Emphasis> Greene and <Emphasis Type="Italic">B. terricola</Emphasis> Kirby (Hymenoptera: Apidae) with implications for conservation status assessments

Conservation action for species of concern requires that “designatable units” (e.g., species, subspecies, geographic races, genetically distinct forms) are clearly defined, or that the species complex is treated as a whole. Several species of bumble bee are currently threatened, and some of these have cryptic colouration (resembling other species), or form complexes that vary considerably in colour patterning. Here we address the taxonomy and distribution of Bombus occidentalis Greene and B. terricola Kirby, both of which are currently of conservation concern in North America. Bombus occidentalis includes two apparently monophyletic groups of COI barcode haplotypes (recently considered as subspecies) with ranges mostly separated by that of their sister species, B. terricola. The southern B. o. occidentalis ranges throughout the western United States and into western Canada from southern Saskatchewan and Alberta, and throughout British Columbia north to ca. 55°N; the northern B. o. mckayi Ashmead, is restricted to north of this in British Columbia, westernmost Northwest Territories, Yukon Territory and Alaska. Bombus o. mckayi exists, as far as is known, only with a “banded” colour pattern. By contrast, B. o. occidentalis occurs in both banded and non-banded colour patterns, although the southern banded colour pattern is geographically isolated from the northern subspecies. Bombus o. occidentalis has declined throughout its range, perhaps due in part to exposure to novel parasites. Despite having similar levels of parasitism (ca. 40 %) as the southern subspecies, B. o. mckayi appears to have stable populations at present. There is therefore compelling evidence that the two subspecies should be distinguished for conservation and management purposes. We present the evidence for their distinction and provide tools for subspecies recognition.     

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.     

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.     

Genetic distinctiveness of the damselfly <Emphasis Type="Italic">Coenagrion puella</Emphasis> in North Africa: an overlooked and endangered taxon

North African odonates are facing conservation challenges, not only by increased degradation and loss of habitat, but also by having poorly understood taxonomy. Coenagrion puella is a widely distributed damselfly but there is debate about the taxonomic status of North African populations, where the species is very rare. We evaluate the genetic distinctiveness of North African C. puella using mitochondrial and nuclear genetic markers. We found a clear genetic differentiation between North African and European populations (3.4 % mtDNA) and a lack of shared haplotypes between individuals from the two continents. These results suggest that the damselfly C. puella comprises two genetically distinct phylogenetic lineages: one in Europe and one in North Africa, and re-invigorate the debate on the validity of the North African endemic C. puella kocheri. We propose that these two lineages of C. puella should be managed as distinct molecular operational taxonomic units. More generally, this study reinforces the important role of North Africa as centre of speciation and differentiation for odonates, and highlights the relevance of incorporating genetic data to understand the evolutionary history and taxonomy for effective biodiversity conservation.     

Characterizing range-wide divergence in an alpine-endemic bird: a comparison of genetic and genomic approaches

The delineation of intraspecific units that are evolutionarily and demographically distinct is an important step in the development of species-specific management plans. Neutral genetic variation has served as the primary data source for delineating “evolutionarily significant units,” but with recent advances in genomic technology, we now have an unprecedented ability to utilize information about neutral and adaptive variation across the entire genome. Here, we use traditional genetic markers (microsatellites) and a newer reduced-representation genomic approach (single nucleotide polymorphisms) to delineate distinct groups of white-tailed ptarmigan (Lagopus leucura), an alpine-obligate species that is distributed in naturally fragmented habitats from Alaska to New Mexico. Five subspecies of white-tailed ptarmigan are currently recognized but their distinctiveness has not been verified with molecular data. Based on analyses of 436 samples at 12 microsatellite loci and 95 samples at 14,866 single nucleotide polymorphism loci, we provide strong support for treating two subspecies as distinct intraspecific units—L. l. altipetens, found in Colorado and neighboring states; and L. l. saxatilis, found on British Columbia’s Vancouver Island—but our findings reveal more moderate patterns of divergence within the remainder of the species’ range. Results based on genetic and genomic datasets generally agreed with one another, indicating that in many cases microsatellite loci may be sufficient for describing major patterns of genetic structure across species’ ranges. This work will inform future conservation and management decisions for the white-tailed ptarmigan, a species that may be vulnerable to future changes in climate.     

Assessing phylogeographic variation in the Rosyside Dace (Teleostei,Leuciscidae), a widespread morphologically variable taxon

Clinostomus (Leuciscidae) is a wide-ranging freshwater fish genus that occurs throughout eastern North America and southern portions of Canada with two species currently recognized: C. elongatus and C. funduloides. A previous taxonomic study of C. funduloides recognized two subspecies (estor and funduloides) and one diagnosed, but undescribed subspecies based on morphological characteristics and geographic distribution. In this study, we used three molecular markers (cytochrome b, S7 intron 1 and growth hormone intron 4) to test the three lineage hypothesis and evaluate genetic variation of C. funduloides across the range using Bayesian inference. Our results indicate that C. funduloides is not monophyletic, as individuals of C. elongatus nest within C. funduloides in both the mtDNA and nDNA phylogenetic analyses, although the position of C. elongatus varies between data sets. In addition, some of the recovered clades are deeply divergent from one another, further supporting the distinctiveness of many of the populations. Overall, these results suggest that subspecies designations are not warranted and a taxonomic revision is needed as Clinostomus is likely more diverse than is currently recognized.     

Tigers of Sundarbans in India: Is the Population a Separate Conservation Unit?

The Sundarbans tiger inhabits a unique mangrove habitat and are morphologically distinct from the recognized tiger subspecies in terms of skull morphometrics and body size. Thus, there is an urgent need to assess their ecological and genetic distinctiveness and determine if Sundarbans tigers should be defined and managed as separate conservation unit. We utilized nine microsatellites and 3 kb from four mitochondrial DNA (mtDNA) genes to estimate genetic variability, population structure, demographic parameters and visualize historic and contemporary connectivity among tiger populations from Sundarbans and mainland India. We also evaluated the traits that determine exchangeability or adaptive differences among tiger populations. Data from both markers suggest that Sundarbans tiger is not a separate tiger subspecies and should be regarded as Bengal tiger (P. t. tigris) subspecies. Maximum likelihood phylogenetic analyses of the mtDNA data revealed reciprocal monophyly. Genetic differentiation was found stronger for mtDNA than nuclear DNA. Microsatellite markers indicated low genetic variation in Sundarbans tigers (He= 0.58) as compared to other mainland populations, such as northern and Peninsular (Hebetween 0.67- 0.70). Molecular data supports migration between mainland and Sundarbans populations until very recent times. We attribute this reduction in gene flow to accelerated fragmentation and habitat alteration in the landscape over the past few centuries. Demographic analyses suggest that Sundarbans tigers have diverged recently from peninsular tiger population within last 2000 years. Sundarbans tigers are the most divergent group of Bengal tigers, and ecologically non-exchangeable with other tiger populations, and thus should be managed as a separate “evolutionarily significant unit” (ESU) following the adaptive evolutionary conservation (AEC) concept.     

Genetic diversity of trout (genus Salmo) from its most eastern native range based on mitochondrial DNA and nuclear gene variation

Russia and western Asia harbour trout populations that have been classified as distinct species and subspecies, most often on the basis of morphological and ecological variation. In order to assess their origins and to verify whether traditional taxonomy reflects their evolutionary distinctiveness, we documented their genetic relationships on the basis of mitochondrial DNA (mtDNA) RFLP, mtDNA sequence analysis, and allozyme variation. Both mtDNA and nuclear gene variation defined two ancient phylogenetic assemblages of populations distributed among northern (Baltic, White, Barents), and southern (Black, Caspian, Aral) sea basins, between which gene flow has been possible but limited in postglacial times. These results supported the traditional taxonomic differentiation between populations of these two regions. They provided weak support for the taxonomic distinction of southern brown trout (Salmo trutta) populations based on their basin of origin. They also refuted the hypothesis that L. Sevan trout (Salmo ischchan) diverged from a primitive brown trout ancestor. Nevertheless, all trout populations from southern sea basins possessed private alleles or mtDNA genotypes and were genetically distinct Therefore, they represent unique gene pools that warrant individual recognition for conservation and management.     

Subspecies and units for conservation and management of the northern bobwhite in the eastern United States

The northern bobwhite (Colinus virginianus) is a small game bird with sedentary lifestyle and has experienced population declines throughout most of its native distribution in the eastern United States. We investigated intraspecific genetic relationships among 14 local populations covering four putative subspecies (C. v. marilandicus, C. v. virginianus, C. v. mexicanus, and C. v. floridanus) in the United States. Analysis of mitochondrial DNA sequences revealed a small, but significant, genetic structure of northern bobwhite populations or subspecies in the eastern US. However, our results did not support current subspecies limits as distinct evolutionarily significant units, based on the amount of population genetic divergences and insufficient lineage sorting of mtDNA haplotypes among subspecies. Instead, our results suggest that C. v. virginianus, C. v. marilandicus, and C. v. mexicanus be merged into a single management unit, and C. v. floridanus be considered as another distinct unit for conservation and management.     

Taxonomy and ecological niche modeling: Implications for the conservation of wood partridges (genus Dendrortyx)

Ecological niche models (ENMs) have a wide range of biological applications, particularly in conservation. To build these models, two sources of information are needed: occurrence records for the species of interest and environmental variables. However, taxonomic limits are often unclear, and the selection of occurrence data depends on the species concept being used. In this study we generated ENMs based on different taxonomic levels within the Dendrortyx group, which is comprised of three species and several subspecies; we analyzed the geographic and ecological distribution patterns and discuss the implications for the biogeography and conservation of this group. Our results suggest that the area with suitable climate depends on the taxonomic category used in the model, which in turn affects the interpretation of the importance of different biogeographic barriers and introduces variation into the potential differentiation of Dendrortyx. In terms of conservation, Dendrortyx macroura and Dendrortyx leucophrys are in a low risk category, that of “least concern,” although they may be amended to a higher category when their allopatric lineages are considered as the units for modeling. We suggest carrying out an a priori taxonomic analysis to facilitate the empirical identification of the units to be modeled in order to allow for a better ecological and biogeographic interpretation and more sound conservation policies.     

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.     

Taxonomic Boundaries and Craniometric Variation in the Treeshrews (Scandentia,Tupaiidae) from the Palawan Faunal Region

The taxonomy of treeshrews (Order Scandentia) has long been complicated by ambiguous morphological species boundaries, and the treeshrews of the Palawan faunal region of the Philippines are no exception. Four named forms in the genus Tupaia Raffles, 1821, have been described from four island groups based on subtle qualitative morphological characters, and as many as three distinct species have been recognized. A recent molecular phylogenetic study of relationships among Tupaia species suggests that the two currently-recognized treeshrew species from the Palawan faunal region diverged very recently relative to other sister-species divergences within the genus and may not represent species-level taxonomic entities. Here we review the taxonomic and biogeographic histories of the Tupaia taxa from this region. We also collected craniodental data from 133 skulls of all four named forms, representing five island populations, and conducted univariate and multivariate analyses on these data. Our morphometric results are consistent with molecular results, further suggesting that there is insufficient evidence to recognize T. moellendorffi Matschie, 1898, as a separate species from T. palawanensis Thomas, 1894. Our analyses also revealed a craniodentally divergent population from the island of Balabac, which has never been considered a distinct subspecies (or species) from the population on Palawan. These results have conservation implications for the island populations in our analyses, but additional surveys and molecular evidence will be required to fully assess conservation priorities for the treeshrews of the Palawan faunal region.     

Molecular survey of morphological subspecies reveals new mitochondrial lineages in Podarcis muralis (Squamata: Lacertidae) from the Tuscan Archipelago (Italy)

Recent analyses of molecular markers have significantly revised the traditional taxonomy of Podarcis species (Squamata: Lacertidae), leading to critically reconsider the taxonomic value of several subspecies described only on morphological bases. In fact, lizards often exhibit high morphological plasticity both at the intra‐specific and the intra‐population level, especially on islands, where phenotypic divergences are mainly due to local adaptation, rather than to evolutionary differentiation. The Common wall lizard Podarcis muralis exhibits high morphological variability in biometry, pholidosis values and colour pattern. Molecular analyses have confirmed the key role played by the Italian Peninsula as a multi‐glacial refuge for P. muralis, pointing out the lack of congruence between mitochondrial lineages and the four peninsular subspecies currently recognized. Here, we analyse a portion of the protein‐encoding cytochrome b gene in the seven subspecies described for the Tuscan Archipelago (Italy), in order to test whether the mitochondrial haplotypes match the morphologically based taxonomy proposed for Common wall lizard. We also compare our haplotypes with all the others from the Italian Peninsula to investigate the presence of unique genetic lineages in insular populations. Our results do not agree completely with the subspecific division based on morphology. In particular, the phylogenetic analyses show that at least four subspecies are characterized by very similar haplotypes and fall into the same monophyletic clade, whereas the other three subspecies are closer to peninsular populations from central Italy. From these results, we conclude that at least some subspecies could be better regarded as simple eco‐phenotypes; in addition, we provide an explanation for the distinctiveness of exclusive lineages found in the archipelago, which constituted a refuge for this species during last glacial periods.     

Genetic variation of complete mitochondrial genome sequences of the Sumatran rhinoceros (<Emphasis Type="Italic">Dicerorhinus sumatrensis</Emphasis>)

The Sumatran rhinoceros (Dicerorhinus sumatrensis) is the smallest and one of the most endangered rhinoceros species, with less than 100 individuals estimated to live in the wild. It was originally divided into three subspecies but only two have survived, D. sumatrensis sumatrensis (Sumatran subspecies), and D. s. harrissoni (Bornean). Questions regarding whether populations of the Sumatran rhinoceros should be treated as different management units to preserve genetic diversity have been raised, particularly in light of its severe decline in the wild and low breeding success in captivity. This work aims to characterize genetic differentiation between Sumatran rhinoceros subspecies using complete mitochondrial genomes, in order to unravel their maternal evolutionary history and evaluate their status as separate management units. We identified three major phylogenetic groups with moderate genetic differentiation: two distinct haplogroups comprising individuals from both the Malay Peninsula and Sumatra, and a third group from Borneo. Estimates of divergence time indicate that the most recent common ancestor of the Sumatran rhinoceros occurred approximately 360,000 years ago. The three mitochondrial haplogroups showed a common divergence time about 80,000 years ago corresponding with a major biogeographic event in the Sundaland region. Patterns of mitochondrial genetic differentiation may suggest considering Sumatran rhinoceros subspecies as different conservation units. However, the management of subspecies as part of a metapopulation may appear as the last resource to save this species from extinction, imposing a conservation dilemma.