Positive selection and precipitation effects on the mitochondrial NADH dehydrogenase subunit 6 gene in brown hares (Lepus europaeus) under a phylogeographic perspective.

Previous studies in hares and jackrabbits have indicated that positive selection has shaped the genetic diversity of mitochondrial genes involved in oxidative phosphorylation, which may affect cellular energy production and cause regional adaptation to different environmental (climatic) pressures. In the present study, we sequenced the NADH dehydrogenase subunit 6 (MT-ND6) gene of 267 brown hares (L. europaeus) from Europe and Asia Minor and tested for positive selection and adaptations acting on amino acid sequences (protein variants). Molecular diversity indices and spatial clustering were assessed by DnaSP, Network, and Geneland, while the presence of selection signals was tested by codeml in PAML, and by using the Datamonkey Adaptive Evolution web server. The SPSS software was used to run multinomial regression models to test for possible effects of climate parameters on the currently obtained protein variants. Fifty-eight haplotypes were revealed with a haplotype diversity of 0.817, coding for 17 different protein variants. The MT-ND6 phylogeographic pattern as determined by the nucleotide sequences followed the earlier found model based on the neutrally evolving D-loop sequences, and reflected the earlier found phylogeographic Late Pleistocene scenario. Based on several selection tests, only one codon position consistently proved to be under positive selection. It did occur exclusively in the evolutionarily younger hares from Europe and it gave rise to several protein variants from the southeastern and south-central Balkans. The occurrence of several of those variants was significantly favored under certain precipitation conditions, as proved by our multinomial regression models. Possibly, the great altitudinal variation in the Balkans may have lead to bigger changes in precipitation across that region and this may have imposed an evolutionarily novel selective pressure on the protein variants and could have led to regional adaptation.

Phylogeography and ecological niche modeling unravel the evolutionary history of the Yarkand hare, Lepus yarkandensis (Mammalia: Leporidae), through the Quaternary.

BACKGROUND: The Taklimakan Desert in China is characterized by unique geological and historical dynamics and endemic flora and fauna, but the influence of historical climate oscillations on the evolutionary history of endemic animals is poorly understood. Lepus yarkandensis is an oases-dependent Near Threatened species that lives in fragmented oasis habitats in the Taklimakan Desert, China. We investigated the geological and climatic impacts on its geographical differentiation, demographic history and influence of Pleistocene glacial-interglacial cycles on the evolutionary history of L. yarkandensis. Further, studied the impact of climatic oscillation based modification on phylogeography, distribution and diversification pattern of Yarkand hare by using Cytb (1140 bp), MGF (592 bp) and SPTBN1 (619 bp) markers. Ecological niche modeling (ENM) revealed the evolutionary history of this species in response to climate change during the Quaternary. Paleodistribution modeling was used to identify putative refugia and estimate their historical distributions. RESULTS: Both historical demographic analyses and climatic niche modeling revealed strong effects of glacial climate changes, suggesting recurrent range contractions and expansions. The EBSP results indicated clear population expansion of L. yarkandensis since the Pleistocene. In the "early Pleistocene", the demographic expansion continued from 0.83 MYA to the last glacial period. The ENM analysis supported a wide distribution of Lepus yarkandensis at high altitudes during the last interglacial (LIG) period. During the last glacial maximum (LGM), the suitable climate was reduced and restricted to the western part of the Taklimakan Desert. CONCLUSIONS: Inland aridification, oasis evolution and river flow played major roles in the population differentiation and demographic history of Yarkand hares. Historically, the large, continuous oases in the Taklimakan Desert contained a viable and unique population of L. yarkandensis. The fragmented desert environment might have caused low gene flow between individuals or groups, thus leading to predominant genetic differentiation. The Pleistocene climatic cycles triggered the diversification and expansion of this species during cold and warm periods, respectively, leading to multiple colonization events within the Taklimakan Desert. These events might be due to the expansion of the Taklimakan Desert during the Middle Pleistocene. Yarkand hare previously occupied vast areas at low and intermediate altitudes in Xinjiang, Gansu, Shanxi, Henan and Shaanxi Provinces in China. The past aridification, climate change-induced oasis modifications, changes in river volumes and flow directions, and human activities all affected the population demography and phylogeography of the Yarkand hare.

Contrasting phylogeographic histories between the continent and islands of East Asia: Massive mitochondrial introgression and long-term isolation of hares (Lagomorpha: Lepus).

Hares of the genus Lepus are distributed worldwide, and introgressive hybridization is thought to be pervasive among species, leading to reticulate evolution and taxonomic confusion. Here, we performed phylogeographic analyses of the following species of hare across East Asia: L. timidus, L. mandshuricus, L. coreanus, and L. brachyurus collected from far-eastern Russia, South Korea, and Japan. Nucleotide sequences of one mitochondrial DNA and eight nuclear gene loci were examined, adding sequences of hares in China from databases. All nuclear DNA analyses supported the clear separation of three phylogroups: L. timidus, L. brachyurus, and the L. mandshuricus complex containing L. coreanus. On the other hand, massive mitochondrial introgression from two L. timidus lineages to the L. mandshuricus complex was suggested in continental East Asia. The northern population of the L. mandshuricus complex was mainly associated with introgression from the continental lineage of L. timidus, possibly since the last glacial period, whereas the southern population of the L. mandshuricus complex experienced introgression from another L. timidus lineage related to the Hokkaido population, possibly before the last glacial period. In contrast to continental hares, no evidence of introgression was found in L. brachyurus in the Japanese Archipelago, which showed the oldest divergence amongst East Asian hare lineages. Our findings suggest that glacial-interglacial climate changes in the circum-Japan Sea region promoted distribution shifts and introgressive hybridization among continental hare species, while the geographic structure of the region contributed to long-term isolation of hares on the islands, preventing inter-species gene flow.

Widespread introgression of mountain hare genes into Fennoscandian brown hare populations.

In Fennoscandia, mountain hare (Lepus timidus) and brown hare (Lepus europaeus) hybridize and produce fertile offspring, resulting in gene flow across the species barrier. Analyses of maternally inherited mitochondrial DNA (mtDNA) show that introgression occur frequently, but unavailability of appropriate nuclear DNA markers has made it difficult to evaluate the scale- and significance for the species. The extent of introgression has become important as the brown hare is continuously expanding its range northward, at the apparent expense of the mountain hare, raising concerns about possible competition. We report here, based on analysis of 6833 SNP markers, that the introgression is highly asymmetrical in the direction of gene flow from mountain hare to brown hare, and that the levels of nuclear gene introgression are independent of mtDNA introgression. While it is possible that brown hares obtain locally adapted alleles from the resident mountain hares, the low levels of mountain hare alleles among allopatric brown hares suggest that hybridization is driven by stochastic processes. Interspecific geneflow with the brown hare is unlikely to have major impacts on mountain hare in Fennoscandia, but direct competition may.

Mitogenomic analysis in European brown hare (Lepus europaeus) proposes genetic and functional differentiation between the distinct lineages.

European brown hare is a small game species spreading across Europe to Asia Minor, with important economic traits. Population genetics studies using mitochondrial DNA markers have revealed the existence of two major phylogeographic lineages, the European and the Anatolian. European lineage is further divided in the European type halpogroup and south-eastern European type halpogroup, while Anatolian consists only by the Anatolian/Middle Eastern type halpogroup. All three haplogroups show a discrete geographical distribution, with an overlapping zone forming in North-East Greece and Bulgaria, forming a contact zone. Despite the existence of a contact zone, European haplotype was never detected in Anatolia and vice versa, proposing the presence of genetic barriers responsible for this phenomenon. In this study, we analyzed the whole mitochondrial genomes of specimens originating from both lineages, aiming to detect the genetic and functional differentiation of the oxidative phosphorylation complexes that are encoded by mtDNA that could lead gradually to the reproductive isolation of the lineages.

Mitochondrial and nuclear DNA reveals reticulate evolution in hares (Lepus spp., Lagomorpha, Mammalia) from Ethiopia.

For hares (Lepus spp., Leporidae, Lagomorpha, Mammalia) from Ethiopia no conclusive molecular phylogenetic data are available. To provide a first molecular phylogenetic model for the Abyssinian Hare (Lepus habessinicus), the Ethiopian Hare (L. fagani), and the Ethiopian Highland Hare (L. starcki) and their evolutionary relationships to hares from Africa, Eurasia, and North America, we phylogenetically analysed mitochondrial ATPase subunit 6 (ATP6; n = 153 / 416bp) and nuclear transferrin (TF; n = 155 / 434bp) sequences of phenotypically determined individuals. For the hares from Ethiopia, genotype composition at twelve microsatellite loci (n = 107) was used to explore both interspecific gene pool separation and levels of current hybridization, as has been observed in some other Lepus species. For phylogenetic analyses ATP6 and TF sequences of Lepus species from South and North Africa (L. capensis, L. saxatilis), the Anatolian peninsula and Europe (L. europaeus, L. timidus) were also produced and additional TF sequences of 18 Lepus species retrieved from GenBank were included as well. Median joining networks, neighbour joining, maximum likelihood analyses, as well as Bayesian inference resulted in similar models of evolution of the three species from Ethiopia for the ATP6 and TF sequences, respectively. The Ethiopian species are, however, not monophyletic, with signatures of contemporary uni- and bidirectional mitochondrial introgression and/ or shared ancestral polymorphism. Lepus habessinicus carries mtDNA distinct from South African L. capensis and North African L. capensis sensu lato; that finding is not in line with earlier suggestions of its conspecificity with L. capensis. Lepus starcki has mtDNA distinct from L. capensis and L. europaeus, which is not in line with earlier suggestions to include it either in L. capensis or L. europaeus. Lepus fagani shares mitochondrial haplotypes with the other two species from Ethiopia, despite its distinct phenotypic and microsatellite differences; moreover, it is not represented by a species-specific mitochondrial haplogroup, suggesting considerable mitochondrial capture by the other species from Ethiopia or species from other parts of Africa. Both mitochondrial and nuclear sequences indicate close phylogenetic relationships among all three Lepus species from Ethiopia, with L. fagani being surprisingly tightly connected to L. habessinicus. TF sequences suggest close evolutionary relationships between the three Ethiopian species and Cape hares from South and North Africa; they further suggest that hares from Ethiopia hold a position ancestral to many Eurasian and North American species.

Variation of partial transferrin sequences and phylogenetic relationships among hares (Lepus capensis, Lagomorpha) from Tunisia.

North African hares are currently included in cape hares, Lepus capensis sensu lato, a taxon that may be considered a superspecies or a complex of closely related species. The existing molecular data, however, are not unequivocal, with mtDNA control region sequences suggesting a separate species status and nuclear loci (allozymes, microsatellites) revealing conspecificity of L. capensis and L. europaeus. Here, we study sequence variation in the intron 6 (468 bp) of the transferrin nuclear gene, of 105 hares with different coat colour from different regions in Tunisia with respect to genetic diversity and differentiation, as well as their phylogenetic status. Forty-six haplotypes (alleles) were revealed and compared phylogenetically to all available TF haplotypes of various Lepus species retrieved from GenBank. Maximum Likelihood, neighbor joining and median joining network analyses concordantly grouped all currently obtained haplotypes together with haplotypes belonging to six different Chinese hare species and the African scrub hare L. saxatilis. Moreover, two Tunisian haploypes were shared with L. capensis, L timidus, L. sinensis, L. yarkandensis, and L. hainanus from China. These results indicated the evolutionary complexity of the genus Lepus with the mixing of nuclear gene haplotypes resulting from introgressive hybridization or/and shared ancestral polymorphism. We report the presence of shared ancestral polymorphism between North African and Chinese hares. This has not been detected earlier in the mtDNA sequences of the same individuals. Genetic diversity of the TF sequences from the Tunisian populations was relatively high compared to other hare populations. However, genetic differentiation and gene flow analyses (AMOVA, FST, Nm) indicated little divergence with the absence of geographically meaningful phylogroups and lack of clustering with coat colour types. These results confirm the presence of a single hare species in Tunisia, but a sound inference on its phylogenetic position would require additional nuclear markers and numerous geographically meaningful samples from Africa and Eurasia.

SNP Identification through Transcriptome Analysis of the European Brown Hare (Lepus europaeus): Cellular Energetics and Mother's Curse.

The European brown hare (Lepus europaeus, Pallas 1778) is an important small game species in Europe. Due to its size and position in the food chain, as well as its life history, phenotypic variation and the relatively recent speciation events, brown hare plays an important role in the structure of various ecosystems and has emerged as an important species for population management and evolutionary studies. In order to identify informative SNPs for such studies, heart and liver tissues of three samples from the European lineage and a three-sample pool from the Anatolian lineage were subjected to RNA-Sequencing analysis. This effort resulted in 9496 well-assembled protein-coding sequences with close homology to human. After applying very stringent filtering criteria, 66185 polymorphic sites were identified in 7665 genes/cds and 2050 of those polymorphic sites are potentially capable of distinguishing the European from the Anatolian lineage. From these distinguishing mutations we focused on those in genes that are involved in cellular energy production, namely the glycolysis, Krebs cycle and the OXPHOS machinery. A selected set of SNPs was also validated by Sanger sequencing. By simulating the three European individuals as one pool, no substantial informative-SNP identification was lost, making it a cost-efficient approach. To our knowledge this is the first attempt to correlate the differentiation in both nuclear and mitochondrial genome between the two different lineages of L. europaeus with the observed spatial partitioning of the lineages of the species, proposing a possible mechanism that is maintaining the reproductive isolation of the lineages.

The complete mitochondrial DNA sequence of the cape hare Lepus capensis pamirensis.

Lepus capensis pamirensis is a specific cape hare species that is only distributed in the Pamirs, with a cruel natural environment. The complete mitochondrial genome of the cape hare L. c. pamirensis was first determined (accession number LC073697), in a total length of 17 597 nucleotides, of which 15 310 nucleotides are coding DNA and 2287 nucleotides are non-coding DNA. It has the common feature with those of other hares with respect to genome structure and gene arrangement. It contained 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer tRNA genes and 1 control region (D-loop region). The complete mitochondrial genome of the cape hare L. c. pamirensis provides an important data set for further study on its classification.

Molecular evidence revealed Lepus hainanus and L. peguensis have a conspecific relationship.

Accurate species delimitation in Lepus was often hindered by highly conserved morphology and frequent introgression. In this study, we used rigorous molecular species delimitation methods to evaluate the taxonomic status of Hainan hare (Lepus hainanus) which has been traditionally identified as a distinct species, or a subspecies of Burmese hare (L. peguensis). The genetic distance and phylogenetic network support L. hainanus and L. peguensis are conspecific. However, the phylogenetic species concept and Bayesian species delimitation analysis based on combined mtDNA supported they are different species. The discordance between different methods can be explained by different species criterion. By taking into account our conflict results, we hold the opinion that adoption of the phylogenetic species concept and Bayesian species delimitation analysis would increase the risk of taxonomic inflation of island biota or otherwise spatially isolated population. Conservatively, we suggest that L. hainanus and L. peguensis are conspecific based on the results of our genetic divergence and phylogenetic network exclusively.

Leporid immunoglobulin G shows evidence of strong selective pressure on the hinge and CH3 domains.

Immunoglobulin G (IgG) is the predominant serum immunoglobulin and has the longest serum half-life of all the antibody classes. The European rabbit IgG has been of significant importance in immunological research, and is therefore well characterized. However, the IgG of other leporids has been disregarded. To evaluate the evolution of this gene in leporids, we sequenced the complete IGHG for six other genera: Bunolagus, Brachylagus, Lepus, Pentalagus, Romerolagus and Sylvilagus. The newly sequenced leporid IGHG gene has an organization and structure similar to that of the European rabbit IgG. A gradient in leporid IgG constant domain diversity was observed, with the CH1 being the most conserved and the CH3 the most variable domain. Positive selection was found to be acting on all constant domains, but with a greater incidence in the CH3 domain, where a cluster of three positively selected sites was identified. In the hinge region, only three polymorphic positions were observed. The same hinge length was observed for all leporids. Unlike the variation observed for the European rabbit, all 11 Lepus species studied share exactly the same hinge motif, suggesting its maintenance as a result of an advantageous structure or conformation.

The hidden history of the snowshoe hare, Lepus americanus: extensive mitochondrial DNA introgression inferred from multilocus genetic variation.

Hybridization drives the evolutionary trajectory of many species or local populations, and assessing the geographic extent and genetic impact of interspecific gene flow may provide invaluable clues to understand population divergence or the adaptive relevance of admixture. In North America, hares (Lepus spp.) are key species for ecosystem dynamics and their evolutionary history may have been affected by hybridization. Here we reconstructed the speciation history of the three most widespread hares in North America - the snowshoe hare (Lepus americanus), the white-tailed jackrabbit (L. townsendii) and the black-tailed jackrabbit (L. californicus) - by analysing sequence variation at eight nuclear markers and one mitochondrial DNA (mtDNA) locus (6240 bp; 94 specimens). A multilocus-multispecies coalescent-based phylogeny suggests that L. americanus diverged ~2.7 Ma and that L. californicus and L. townsendii split more recently (~1.2 Ma). Within L. americanus, a deep history of cryptic divergence (~2.0 Ma) was inferred, which coincides with major speciation events in other North American species. While the isolation-with-migration model suggested that nuclear gene flow was generally rare or absent among species or major genetic groups, coalescent simulations of mtDNA divergence revealed historical mtDNA introgression from L. californicus into the Pacific Northwest populations of L. americanus. This finding marks a history of past reticulation between these species, which may have affected other parts of the genome and influence the adaptive potential of hares during climate change.

Pseudogenization of the MCP-2/CCL8 chemokine gene in European rabbit (genus Oryctolagus), but not in species of Cottontail rabbit (Sylvilagus) and Hare (Lepus).

BACKGROUND: Recent studies in human have highlighted the importance of the monocyte chemotactic proteins (MCP) in leukocyte trafficking and their effects in inflammatory processes, tumor progression, and HIV-1 infection. In European rabbit (Oryctolagus cuniculus) one of the prime MCP targets, the chemokine receptor CCR5 underwent a unique structural alteration. Until now, no homologue of MCP-2/CCL8a, MCP-3/CCL7 or MCP-4/CCL13 genes have been reported for this species. This is interesting, because at least the first two genes are expressed in most, if not all, mammals studied, and appear to be implicated in a variety of important chemokine ligand-receptor interactions. By assessing the Rabbit Whole Genome Sequence (WGS) data we have searched for orthologs of the mammalian genes of the MCP-Eotaxin cluster. RESULTS: We have localized the orthologs of these chemokine genes in the genome of European rabbit and compared them to those of leporid genera which do (i.e. Oryctolagus and Bunolagus) or do not share the CCR5 alteration with European rabbit (i.e. Lepus and Sylvilagus). Of the Rabbit orthologs of the CCL8, CCL7, and CCL13 genes only the last two were potentially functional, although showing some structural anomalies at the protein level. The ortholog of MCP-2/CCL8 appeared to be pseudogenized by deleterious nucleotide substitutions affecting exon1 and exon2. By analyzing both genomic and cDNA products, these studies were extended to wild specimens of four genera of the Leporidae family: Oryctolagus, Bunolagus, Lepus, and Sylvilagus. It appeared that the anomalies of the MCP-3/CCL7 and MCP-4/CCL13 proteins are shared among the different species of leporids. In contrast, whereas MCP-2/CCL8 was pseudogenized in every studied specimen of the Oryctolagus - Bunolagus lineage, this gene was intact in species of the Lepus - Sylvilagus lineage, and was, at least in Lepus, correctly transcribed. CONCLUSION: The biological function of a gene was often revealed in situations of dysfunction or gene loss. Infections with Myxoma virus (MYXV) tend to be fatal in European rabbit (genus Oryctolagus), while being harmless in Hares (genus Lepus) and benign in Cottontail rabbit (genus Sylvilagus), the natural hosts of the virus. This communication should stimulate research on a possible role of MCP-2/CCL8 in poxvirus related pathogenicity.

Recurrent introgression of mitochondrial DNA among hares (Lepus spp.) revealed by species-tree inference and coalescent simulations.

Understanding recent speciation history requires merging phylogenetic and population genetics approaches, taking into account the persistence of ancestral polymorphism and possible introgression. The emergence of a clear phylogeny of hares (genus Lepus) has been hampered by poor genomic sampling and possible occurrence of mitochondrial DNA (mtDNA) introgression from the arctic/boreal Lepus timidus into several European temperate and possibly American boreal species. However, no formal test of introgression, taking also incomplete lineage sorting into account, has been done. Here, to clarify the yet poorly resolved species phylogeny of hares and test hypotheses of mtDNA introgression, we sequenced 14 nuclear DNA and 2 mtDNA fragments (8205 and 1113 bp, respectively) in 50 specimens from 11 hare species from Eurasia, North America, and Africa. By applying an isolation-with-migration model to the nuclear data on subsets of species, we find evidence for very limited gene flow from L. timidus into most temperate European species, and not into the American boreal ones. Using a multilocus coalescent-based method, we infer the species phylogeny, which we find highly incongruent with mtDNA phylogeny using parametric bootstrap. Simulations of mtDNA evolution under the speciation history inferred from nuclear genes did not support the hypothesis of mtDNA introgression from L. timidus into the American L. townsendii but did suggest introgression from L. timidus into 4 temperate European species. One such event likely resulted in the complete replacement of the aboriginal mtDNA of L. castroviejoi and of its sister species L. corsicanus. It is remarkable that mtDNA introgression in hares is frequent, extensive, and always from the same donor arctic species. We discuss possible explanations for the phenomenon in relation to the dynamics of range expansions and species replacements during the climatic oscillations of the Pleistocene.

Interspecific X-chromosome and mitochondrial DNA introgression in the Iberian hare: selection or allele surfing?

Introgression from a resident species into an invading one is predicted to occur through the demographic process of "allele surfing," and to particularly affect genomic regions transmitted by the lower migrating sex, such as mtDNA. This could explain that northern Iberian populations of Lepus granatensis harbor high frequencies of mtDNA from L. timidus, an arctic hare it replaced there after deglaciation. We report that variation of introgressed timidus-like mtDNA reflects several predicted effects of this process: increasing frequency and diversity in the direction of expansion, strong perpendicular phylogeographic structure and signs of postglacial demographic growth. However, demographic inferences for the granatensis and timidus-like mtDNA lineages suggest the latter may have outcompeted the former in northern Iberia. Autosomal introgression occurs at low frequencies and species-wide rather than only in the north. If this difference with mtDNA resulted from sex-biased migration, an intermediate pattern should prevail for the X-chromosome, but we report species-wide and high-frequency introgression of an X-fragment. Either selection favored this ubiquitous X-introgression, or more complex postglacial expansion patterns prevailed, with different consequences depending on the genomic and geographic region. This illustrates the difficulty of distinguishing demographic and selective effects and the need for genome and species-wide based demographic models.

Reticulate evolution: frequent introgressive hybridization among Chinese hares (genus lepus) revealed by analyses of multiple mitochondrial and nuclear DNA loci.

BACKGROUND: Interspecific hybridization may lead to the introgression of genes and genomes across species barriers and contribute to a reticulate evolutionary pattern and thus taxonomic uncertainties. Since several previous studies have demonstrated that introgressive hybridization has occurred among some species within Lepus, therefore it is possible that introgressive hybridization events also occur among Chinese Lepus species and contribute to the current taxonomic confusion. RESULTS: Data from four mtDNA genes, from 116 individuals, and one nuclear gene, from 119 individuals, provides the first evidence of frequent introgression events via historical and recent interspecific hybridizations among six Chinese Lepus species. Remarkably, the mtDNA of L. mandshuricus was completely replaced by mtDNA from L. timidus and L. sinensis. Analysis of the nuclear DNA sequence revealed a high proportion of heterozygous genotypes containing alleles from two divergent clades and that several haplotypes were shared among species, suggesting repeated and recent introgression. Furthermore, results from the present analyses suggest that Chinese hares belong to eight species. CONCLUSION: This study provides a framework for understanding the patterns of speciation and the taxonomy of this clade. The existence of morphological intermediates and atypical mitochondrial gene genealogies resulting from frequent hybridization events likely contribute to the current taxonomic confusion of Chinese hares. The present study also demonstrated that nuclear gene sequence could offer a powerful complementary data set with mtDNA in tracing a complete evolutionary history of recently diverged species.

[Genetic structure and subspecies divergence of Lepus capensis in Xinjiang].

To date, the genetic structure and genetic diversity of Lepus capensis in Xinjiang has not been systematically studied at the molecular level, and its subspecies taxonomic status has been under debate for years. According to traditional morphology, there are three subspecies of L. capensis distributed in Xinjiang: L.c. centrasiaticus, L.c. lehmanni and L.c. pamirensis. In this study, we determined 592 bp D-loop sequences of 87 cape hares from Xinjiang Province. Forty-four haplotypes were defined based on 148 polymorphic sites. Both the haplotype diversity (0.977+/-0.005) and nucleotide diversity (0.064+/-0.031) are high. FST P values are significantly high and no haplotype was shared among the four geographic populations, indicating that genetic differentiation among populations is significant. AMOVA shows that most of the genetic differentiation occurred among geographic groups, indicating that geographic isolation such as mountains and deserts might make an effective barrier against gene flow. Both the phylogenetic tree and median-joining network grouped 44 haplotypes into four distinct clades corresponding to four geographic areas, indicating an obvious phylogeographic pattern. Our data supported the subspecies status of L. c. lehmanni. The fact that haplotypes of L. c. centrasiaticus were grouped into two distinct clades suggests that this traditional subspecies should be considered as two subspecies. In addition, L. c. pamirensis shows a significantly higher sequence divergence compared to other subspecies, and the difference even reached the level of species.

Bidirectional introgressive hybridization between Lepus capensis and Lepus yarkandensis.

Lepus capensis and Lepus yarkandensis are two well-established hare species with different habitat preferences and obvious morphological distinctions. L. capensis is common and widespread while L. yarkandensis is an endemic species with a restricted distribution in the Tarim Basin of western China. During field sampling, individuals with a distinct phenotype were found in the contact zone between the two species. To understand the origin of these unclassified individuals, we analyzed morphological and genetic data from 700 individuals including L. yarkandensis, L. capensis and the unclassified individuals. Morphological analyses of the unclassified individuals revealed that they lie between L. yarkandensis and L. capensis. Genetic analyses based on species-specific molecular markers (mitochondrial control region and SRY) showed that the unclassified populations have mtDNA and SRY genes from both L. capensis and L. yarkandensis, suggesting that the unclassified populations result from hybridization between the two species. Interestingly, many identical and/or very similar mtDNA haplotypes were shared between L. capensis populations and L. yarkandensis populations close to their contact zones, which further suggests the existence of extensive bidirectional mtDNA introgression. Similarly, we found evidence for SRY introgression, though it existed at a lower level compared to mtDNA introgression.

A phantom extinction? New insights into extinction dynamics of the Don-hare Lepus tanaiticus.

The Pleistocene to Holocene transition was accompanied by a worldwide extinction event affecting numerous mammalian species. Several species such as the woolly mammoth and the giant deer survived this extinction wave, only to go extinct a few thousand years later during the Holocene. Another example for such a Holocene extinction is the Don-hare, Lepus tanaiticus, which inhabited the Russian plains during the late glacial. After being slowly replaced by the extant mountain hare (Lepus timidus), it eventually went extinct during the middle Holocene. Here, we report the phylogenetic relationship of L. tanaiticus and L. timidus based on a 339-basepair (bp) fragment of the mitochondrial D-loop. Phylogenetic tree- and network reconstructions do not support L. tanaiticus and L. timidus being different species. Rather, we suggest that the two taxa represent different morphotypes of a single species and the extinction of 'L. tanaiticus' represents the disappearance of a local morphotype rather than the extinction of a species.

The genomic legacy from the extinct Lepus timidus to the three hare species of Iberia: contrast between mtDNA, sex chromosomes and autosomes.

Extensive interspecific genetic introgression is often reported, and appraising its genomic impact can serve to determine whether it results from selection on specific loci or from demographic processes affecting the whole genome. The three species of hares present in the Iberian Peninsula harbour high frequencies of mitochondrial DNA (mtDNA) from Lepus timidus, an arctic/boreal species now extinct in the region. This could result from the invasive replacement of L. timidus by the temperate species during deglaciation but should then have left traces in the nuclear genome. We typed single nucleotide polymorphisms (SNPs) discovered by sequencing 10 autosomal loci, two X-linked and one Y-linked in species-wide samples of the four taxa. Based on lineage-diagnostic SNPs, we detected no trace of L. timidus sex chromosomes in Iberia. From the frequencies of inferred haplotypes, autosomal introgression into L. granatensis appeared mostly sporadic but always widespread instead of restricted to the north as mtDNA. Autosomal introgression into Iberian L. europaeus, inhabiting the Pyrenean foothills, was hardly detectable, despite quasi-fixation of L. timidus mtDNA. L. castroviejoi, endemic to the Cantabrian Mountains and fixed for L. timidus mtDNA, showed little traces of autosomal introgression. The absence of sex-chromosome introgression presumably resulted from X-linked hybrid male unfitness. The contrasting patterns between the autosomes and mtDNA could reflect general gender asymmetric processes such as frequency-dependent female assortative mating, lower mtDNA migration and higher male dispersal, but adaptive mtDNA introgression cannot be dismissed. Additionally, we document reciprocal introgression between L. europaeus and both L. granatensis in Iberia and L. timidus outside Iberia.