Ancient introgression of Lepus timidus mtDNA into L. granatensis and L. europaeus in the Iberian Peninsula

Evolutionary characteristics of the mitochondrial NADH dehydrogenase subunit 6 gene in some populations of four sympatric Mustela species (Mustelidae, Mammalia) from central Europe

BACKGROUND

Selection on or reticulate evolution of mtDNA is documented in various mammalian taxa and could lead to misleading phylogenetic conclusions if not recognized. We sequenced the MT-ND6 gene of four sympatric Mustelid species of the genus Mustela from some central European populations. We hypothesised positive selection on MT-ND6, given its functional importance and the different body sizes and life histories of the species, even though climatic differences may be unimportant for adaptation in sympatry.

METHODS AND RESULTS

MT-ND6 genes were sequenced in 187 sympatric specimens of weasels, Mustela nivalis, stoats, M. erminea, polecats, M. putorius, and steppe polecats, M. eversmannii, from eastern Austria and of fourteen allopatric polecats from eastern-central Germany. Median joining networks, neighbour joining and maximum likelihood analyses as well as Bayesian inference grouped all species according to earlier published phylogenetic models. However, polecats and steppe polecats, two very closely related species, shared the same two haplotypes. We found only negative selection within the Mustela sequences, including 131 downloaded ones covering thirteen species. Positive selection was observed on three MT-ND6 codons of other mustelid genera retrieved from GenBank.

CONCLUSIONS

Negative selection for MT-ND6 within the genus Mustela suggests absence of both environmental and species-specific effects of cellular energy metabolism despite large species-specific differences in body size. The presently found shared polymorphism in European polecats and steppe polecats may result from ancestral polymorphism before speciation and historical or recent introgressive hybridization; it may indicate mtDNA capture of steppe polecats by M. putorius in Europe.

Bidirectional hybridisation and introgression between introduced European brown hare, Lepus europaeus and the endemic Irish hare, L. timidus hibernicus

Introduced non-native species can threaten native species through interspecific hybridisation and genetic introgression. We assessed the prevalence of hybridisation and introgression between introduced European brown hare, Lepus europaeus, and the endemic Irish hare, L. timidus hibernicus. Roadkill hares (n = 56) were sequenced for a 379bp section of the mitochondrial DNA D-loop and a 474bp segment of the nuclear transferrin (Tf) gene. A species-specific indel in the transferrin gene was present in L.t. hibernicus and absent in L. europaeus. Excluding three hares from which molecular data could not be recovered, 28 hares (53%) were native L.t. hibernicus, 7 (13%) were non-native L. europaeus and 18 (34%) were hybrids; of which 5 (28%) were first generation (F1) involving bidirectional crosses with mismatched nuclear and mtDNA (3 ♂ europaeus x ♀ hibernicus and 2 ♂ hibernicus x ♀ europaeus). Mixed nuclear transferrin sequences suggested 13 (72%) of hybrids were at least 2nd generation (F2) with 9 (69%) possessing L.t. hibernicus and 4 (31%) L. europaeus mtDNA (the latter indicative of hybrid backcrossing with the non-native). The prevalence of hybridisation at similar mountain-brown hare contact zones throughout Europe is notably lower (4–16%) and typically unidirectional (♂ europaeus x ♀ timidus). A high prevalence of bidirectional hybridisation and introgression (in association with projected climate change) may favour the introduced species over the native. Genetic surveillance and population monitoring are needed to further explore the potential conservation implications of European brown hare in Ireland.

Lagomorpha as a Model Morphological System

Due to their global distribution, invasive history, and unique characteristics, European rabbits are recognizable almost anywhere on our planet. Although they are members of a much larger group of living and extinct mammals [Mammalia, Lagomorpha (rabbits, hares, and pikas)], the group is often characterized by several well-known genera (e.g., Oryctolagus, Sylvilagus, Lepus, and Ochotona). This representation does not capture the extraordinary diversity of behavior and form found throughout the order. Model organisms are commonly used as exemplars for biological research, but there are a limited number of model clades or lineages that have been used to study evolutionary morphology in a more explicitly comparative way. We present this review paper to show that lagomorphs are a strong system in which to study macro- and micro-scale patterns of morphological change within a clade that offers underappreciated levels of diversity. To this end, we offer a summary of the status of relevant aspects of lagomorph biology.

Contribution to the ecology of the Italian hare (Lepus corsicanus)

The Italian hare (Lepus corsicanus) is endemic to Central-Southern Italy and Sicily, classified as vulnerable due to habitat alterations, low density and fragmented populations and ecological competition with the sympatric European hare (Lepus europaeus). Despite this status, only few and local studies have explored its ecological features. We provided some key traits of the ecological niche of the Italian hare as well as its potential distribution in the Italian peninsula. All data derived from genetically validated presences. We generated a habitat suitability model using maximum entropy distribution model for the Italian hare and its main competitor, the European hare. The dietary habits were obtained for the Italian hare with DNA metabarcoding and High-Throughput Sequencing on faecal pellets. The most relevant environmental variables affecting the potential distribution of the Italian hare are shared with the European hare, suggesting a potential competition. The variation in the observed altitudinal distribution is statistically significant between the two species.The diet of the Italian hare all year around includes 344 plant taxa accounted by 62 families. The Fagaceae, Fabaceae, Poaceae, Rosaceae and Solanaceae (counts > 20,000) represented the 90.22% of the total diet. Fabaceae (60.70%) and Fagaceae (67.47%) were the most abundant plant items occurring in the Spring/Summer and Autumn/Winter diets, respectively. The Spring/Summer diet showed richness (N = 266) and diversity index values (Shannon: 2.329, Evenness: 0.03858, Equitability: 0.4169) higher than the Autumn/Winter diet (N = 199, Shannon: 1.818, Evenness: 0.03096, Equitability: 0.3435). Our contribution adds important information to broaden the knowledge on the environmental (spatial and trophic) requirements of the Italian hare, representing effective support for fitting management actions in conservation planning.

An Annotated Draft Genome of the Mountain Hare (Lepus timidus)

Hares (genus Lepus) provide clear examples of repeated and often massive introgressive hybridization and striking local adaptations. Genomic studies on this group have so far relied on comparisons to the European rabbit (Oryctolagus cuniculus) reference genome. Here, we report the first de novo draft reference genome for a hare species, the mountain hare (Lepus timidus), and evaluate the efficacy of whole-genome re-sequencing analyses using the new reference versus using the rabbit reference genome. The genome was assembled using the ALLPATHS-LG protocol with a combination of overlapping pair and mate-pair Illumina sequencing (77x coverage). The assembly contained 32,294 scaffolds with a total length of 2.7 Gb and a scaffold N50 of 3.4 Mb. Re-scaffolding based on the rabbit reference reduced the total number of scaffolds to 4,205 with a scaffold N50 of 194 Mb. A correspondence was found between 22 of these hare scaffolds and the rabbit chromosomes, based on gene content and direct alignment. We annotated 24,578 protein coding genes by combining ab-initio predictions, homology search, and transcriptome data, of which 683 were solely derived from hare-specific transcriptome data. The hare reference genome is therefore a new resource to discover and investigate hare-specific variation. Similar estimates of heterozygosity and inferred demographic history profiles were obtained when mapping hare whole-genome re-sequencing data to the new hare draft genome or to alternative references based on the rabbit genome. Our results validate previous reference-based strategies and suggest that the chromosome-scale hare draft genome should enable chromosome-wide analyses and genome scans on hares.

Population Genetics and Speciation of Yellow-Bellied, Red-Naped, and Red-Breasted Sapsuckers (Sphyrapicus varius, S. nuchalis, and S. ruber)

The root of understanding speciation lies in determining the forces which drive it. In many closely-related species, including Sphyrapicus varius, S. nuchalis, and S. ruber, it is assumed that speciation occurred due to isolation in multiple Pleistocene refugia. We used genetic data from 457 samples at the control region (CR), cytochrome oxidase I (COI), and chromo-helicase DNA binding protein (CHD1Z) to examine range-wide population genetic structure and differentiation amongst these 3 species across each species' breeding range. In addition, we modeled these species' ecological niches for the Holocene (~6000 years ago), Last Glacial Maximum (~22000 years ago), and Last Interglacial (~120000-140000 years ago) to determine if Pleistocene glaciations could have contributed to allopatric distributions, therefore allowing these groups to differentiate. Population genetic data show a potential Pleistocene refugium in Haida Gwaii, an east-west split among S. varius, and low genetic differentiation within each species. Our CR data show some polyphyly, while COI and CHD1Z data show differentiation among species using composite genotypes. Ecological Niche Modeling shows a large amount of niche overlap at each time period suggesting that S. varius, S. nuchalis, and S. ruber may not have been completely allopatric, and these species likely had repeated intermittent contact. Our data support the growing body of research that suggests differentiation despite gene flow.

The evolutionary history of the Cape hare (Lepus capensis sensu lato): insights for systematics and biogeography

Inferring the phylogeography of species with large distributions helps deciphering major diversification patterns that may occur in parallel across taxa. Here, we infer the evolutionary history of the Cape hare, Lepus capensis sensu lato, a species distributed from southern Africa to Asia, by analyzing variation at 18 microsatellites and 9 DNA (1 mitochondrial and 8 nuclear) sequenced loci, from field and museum-collected samples. Using a combination of assignment and coalescent-based methods, we show that the Cape hare is composed of five evolutionary lineages, distributed in distinct biogeographic regions-north-western Africa, eastern Africa, southern Africa, the Near East and the Arabian Peninsula. A deep phylogenetic break possibly dating to the Early Pleistocene was inferred between the African and Asian L. capensis groups, and the latter appear more closely related to other Eurasian hare species than to African Cape hares. The inferred phylogeographic structure is shared by numerous taxa distributed across the studied range, suggesting that environmental changes, such as the progressive aridification of the Saharo-Arabian desert and the fluctuations of savannah habitats in Sub-Saharan Africa, had comparable impacts across species. Fine-scale analyses of the western Sahara-Sahel populations showed rich fragmentation patterns for mitochondrial DNA but not for microsatellites, compatible with the environmental heterogeneity of the region and female philopatry. The complex evolutionary history of L. capensis sensu lato, which possibly includes interspecific gene flow, is not reflected by taxonomy. Integrating evolutionary inference contributes to an improved characterization of biodiversity, which is fundamental to foster the conservation of relevant evolutionary units.

Large-scale mitochondrial DNA analysis reveals new light on the phylogeography of Central and Eastern-European Brown hare (Lepus europaeus Pallas, 1778)

European brown hare, Lepus europaeus, from Central and Eastern European countries (Hungary, Poland, Serbia, Lithuania, Romania, Georgia and Italy) were sampled, and phylogenetic analyses were carried out on two datasets: 1.) 137 sequences (358 bp) of control region mtDNA; and 2.) 105 sequences of a concatenated fragment (916 bp), including the cytochrome b, tRNA-Thr, tRNA-Pro and control region mitochondrial DNA. Our sequences were aligned with additional brown hare sequences from GenBank. A total of 52 and 51 haplotypes were detected within the two datasets, respectively, and assigned to two previously described major lineages: Anatolian/Middle Eastern (AME) and European (EUR). Furthermore, the European lineage was divided into two subclades including South Eastern European (SEE) and Central European (CE). Sympatric distribution of the lineages of the brown hare in South-Eastern and Eastern Europe revealed contact zones there. BAPS analysis assigned sequences from L. europaeus to five genetic clusters, whereas CE individuals were assigned to only one cluster, and AME and SEE sequences were each assigned to two clusters. Our findings uncover numerous novel haplotypes of Anatolian/Middle Eastern brown hare outside their main range, as evidence for the combined influence of Late Pleistocene climatic fluctuations and anthropogenic activities in shaping the phylogeographic structure of the species. Our results support the hypothesis of a postglacial brown hare expansion from Anatolia and the Balkan Peninsula to Central and Eastern Europe, and suggest some slight introgression of individual haplotypes from L. timidus to L. europaeus.

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.

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.

Multilocus approaches reveal underestimated species diversity and inter-specific gene flow in pikas (Ochotona) from southwestern China

The phylogeny of living pikas (Ochotonidae, Ochotona) remains obscure, and pika species diversity in southwestern China has never been well explored. In this study, 96 tissue samples from 11 valid species in three classified subgenera (Pika, Ochotona and Conothoa) from 23 locations were characterized using multilocus sequences of 7031bp. Two mitochondrial (CYT B and COI) and five nuclear gene segments (RAG1, RAG2, TTN, OXAIL and IL1RAPL1) were sequenced. We analysed evolutionary histories using maximum likelihood (RAxML) and Bayesian analyses (BEAST), and we also used molecular species delimitation analyses (BPP) to explore species diversity. Our study supported O. syrinx (O. huangensis) as a distinct clade from all named subgenera. Relationships among subgenera were not fully resolved, which may be due to a rapid diversification in the middle Miocene (∼13.90Ma). Conflicting gene trees implied mitochondrial introgression from O. cansus to O. curzoniae. We uncovered three cryptic species from Shaanxi, Sichuan and Yunnan with strong support, suggesting an underestimation of species diversity in the "sky-island" mountains of southwest China.

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, F_ST, 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.

Identification of blood meals in field captured sand flies by a PCR-RFLP approach based on cytochrome b gene

Leishmaniasis is a vector-borne disease transmitted by phlebotomine sand flies. Information about blood meal preferences in sand flies is essential to understand the epidemiology of the disease to adopt control measures. In previous studies, a polymerase chain reaction (PCR) of 359bp fragment of the conserved gene cytochrome b (cyt b) and further sequencing were applied in the study of blood meal sources in sand flies collected in the area of a leishmaniasis outbreak in southwest Madrid, Spain, providing significant information about blood meal preferences in the focus. In this work, a PCR-restriction fragment length polymorphism (RFLP) targeting a fragment of 359bp of vertebrate cyt b gene was developed. Restriction endonucleases HaeIII and HinfI generated specific patterns consistent with the blood meal sources found in sand flies. The protocol has been validated with twenty six engorged females collected in the field with CDC traps. Blood meals from nine vertebrates were identified based on PCR-cyt b and sequencing-human, dog, cat, horse, hare, rabbit, sheep, goat and chicken - and mixed blood meals (sheep/human; sheep/goat) - and successfully distinguished by PCR-RFLP. Therefore, this approach is an efficient and reliable alternative method to be applied in entomological surveys.

Conservation implications of the evolutionary history and genetic diversity hotspots of the snowshoe hare

With climate warming, the ranges of many boreal species are expected to shift northward and to fragment in southern peripheral ranges. To understand the conservation implications of losing southern populations, we examined range-wide genetic diversity of the snowshoe hare (Lepus americanus), an important prey species that drives boreal ecosystem dynamics. We analysed microsatellite (8 loci) and mitochondrial DNA sequence (cytochrome b and control region) variation in almost 1000 snowshoe hares. A hierarchical structure analysis of the microsatellite data suggests initial subdivision in two groups, Boreal and southwestern. The southwestern group further splits into Greater Pacific Northwest and U.S. Rockies. The genealogical information retrieved from mtDNA is congruent with the three highly differentiated and divergent groups of snowshoe hares. These groups can correspond with evolutionarily significant units that might have evolved in separate refugia south and east of the Pleistocene ice sheets. Genetic diversity was highest at mid-latitudes of the species' range, and genetic uniqueness was greatest in southern populations, consistent with substructuring inferred from both mtDNA and microsatellite analyses at finer levels of analysis. Surprisingly, snowshoe hares in the Greater Pacific Northwest mtDNA lineage were more closely related to black-tailed jackrabbits (Lepus californicus) than to other snowshoe hares, which may result from secondary introgression or shared ancestral polymorphism. Given the genetic distinctiveness of southern populations and minimal gene flow with their northern neighbours, fragmentation and loss of southern boreal habitats could mean loss of many unique alleles and reduced evolutionary potential.

Home-loving boreal hare mitochondria survived several invasions in Iberia: the relative roles of recurrent hybridisation and allele surfing

Genetic introgression from a resident species into an invading close relative can result from repeated hybridisation along the invasion front and/or allele surfing on the expansion wave. Cases where the phenomenon is massive and systematic, such as for hares (genus Lepus) in Iberia, would be best explained by recurrent hybridisation but this is difficult to prove because the donor populations are generally extinct. In the Pyrenean foothills, Lepus europaeus presumably replaced Lepus granatensis recently and the present species border is parallel to the direction of invasion, so that populations of L. granatensis in the contact zone represent proxies of existing variation before the invasion. Among three pairs of populations sampled across this border, we find less differentiation of mitochondrial DNA (mtDNA) across than along it, as predicted under recurrent hybridisation at the invasion front. Using autosomal microsatellite loci and X- and Y-linked diagnostic loci, we show that admixture across the border is quasi-absent, making it unlikely that lack of interspecific mtDNA differentiation results from ongoing gene flow. Furthermore, we find that the local species ranges are climatically contrasted, making it also unlikely that ongoing ecology-driven movement of the contact account for mtDNA introgression. The lack of mtDNA differentiation across the boundary is mostly due to sharing of mtDNA from a boreal species currently extinct in Iberia (Lepus timidus) whose mitochondria have thus remained in place since the last deglaciation despite successive invasions by two other species. Home-loving mitochondria thus witness past species distribution rather than ongoing exchanges across stabilised contact zones.

The new French 2010 Rabbit Hemorrhagic Disease Virus causes an RHD-like disease in the Sardinian Cape hare (Lepus capensis mediterraneus)

Lagovirus is an emerging genus of Caliciviridae, which includes the Rabbit Hemorrhagic Disease Virus (RHDV) of rabbits and the European brown hare syndrome virus (EBHSV) of hares that cause lethal hepatitis. In 2010, a new RHDV related virus (RHDV2) with a unique genetic and antigenic profile and lower virulence was identified in France in rabbits. Here we report the identification of RHDV2 as the cause in Sardinia of several outbreaks of acute hepatitis in rabbits and Cape hare (Lepus capensis mediterraneus). This is the first account of a lagovirus that causes fatal hepatitis in both rabbits and hares.

Polarisation of major histocompatibility complex II host genotype with pathogenesis of European Brown Hare syndrome virus

A study was conducted in order to determine the occurrence of European Brown Hare Syndrome virus (EBHSV) in Denmark and possible relation between disease pathogenesis and Major Histocompatibility Complex (MHC) host genotype. Liver samples were examined from 170 brown hares (hunted, found sick or dead), collected between 2004 and 2009. Macroscopical and histopathological findings consistent with EBHS were detected in 24 (14.1%) hares; 35 (20.6%) had liver lesions not typical of the syndrome, 50 (29.4%) had lesions in other tissues and 61 (35.9%) had no lesions. Sixty five (38.2%) of 170 samples were found to be EBHSV-positive (RT-PCR, VP60 gene). In order to investigate associations between viral pathogenesis and host genotype, variation within the exon 2 DQA gene of MHC was assessed. DQA exon 2 analysis revealed the occurrence of seven different alleles in Denmark. Consistent with other populations examined so far in Europe, observed heterozygosity of DQA (H_o = 0.1180) was lower than expected (H_e = 0.5835). The overall variation for both nucleotide and amino acid differences (2.9% and 14.9%, respectively) were lower in Denmark than those assessed in other European countries (8.3% and 16.9%, respectively). Within the peptide binding region codons the number of nonsynonymous substitutions (dN) was much higher than synonymous substitutions (dS), which would be expected for MHC alleles under balancing selection. Allele frequencies did not significantly differ between EBHSV-positive and -negative hares. However, allele Leeu-DQA*30 was detected in significantly higher (P = 0.000006) frequency among the positive hares found dead with severe histopathological lesions than among those found sick or apparently healthy. In contrast, the latter group was characterized by a higher frequency of the allele Leeu-DQA*14 as well as the proportion of heterozygous individuals (P = 0.000006 and P = 0.027). These data reveal a polarisation between EBHSV pathogenesis and MHC class II genotype within the European brown hare in Denmark.

Detection of Leishmania infantum and identification of blood meals in Phlebotomus perniciosus from a focus of human leishmaniasis in Madrid, Spain

Since 2010, the number of cases of both human visceral leishmaniasis and cutaneous leishmaniasis in southwestern Madrid region (Spain) and more specifically in the town of Fuenlabrada has increased. Direct xenodiagnosis of leishmaniasis proved that hares (Lepus granatensis) from this focus are able to infect with Leishmania infantum colonized Phlebotomus perniciosus. To a better understanding of this focus of leishmaniasis, we conducted an entomological survey using CDC light traps, at the end of the seasonal transmission period of 2011 before the beginning of control measures of the disease, to study the phlebotomine sand flies species involved. Detection of Leishmania DNA in the sand flies captured was studied by kDNA-PCR and cpb-PCR. In addition, blood fed and gravid female P. perniciosus were analysed by a PCR based in vertebrate cytochrome b (cyt b) gene. Taxonomic identification of captured sand flies (n = 174) as P. perniciosus (n = 171) and Sergentomyia minuta (n = 3) together with the analysis of blood feeding in ten sand flies that shows a high preference for hares (n = 6), followed by humans (n = 3), and cats (n = 1) confirm a strong association between P. perniciosus hares and humans in the focus. Moreover, 79 out of 135 (58.5 %) P. perniciosus were positive to L. infantum by PCR approaches. These data support the increase of human leishmaniasis cases in the area and the existence of an unusual sylvatic cycle alternative to the classical domestic one, where the dog is the main reservoir of L. infantum.

Ancient colonization and within-island vicariance revealed by mitochondrial DNA phylogeography of the mountain hare (Lepus timidus) in Hokkaido, Japan

We examined the phylogenetic status and history of the mountain hare Lepus timidus in and around Hokkaido using mitochondrial cytochrome b (cyt b) sequences from 158 samples from Hokkaido and 14 from Sakhalin, as well as four samples from the Korean hare, L. coreanus. The phylogenetic analysis of the cyt b sequences generated in this study and obtained from DNA databases showed the clear genetic specificity of the Hokkaido lineage as a clade. The Hokkaido lineage was estimated to have diverged from the other conspecific and L. coreanus lineages 0.46 and 0.30 million years ago (Mya), respectively. These results suggest that the common ancestor of the mitochondrial lineage in Hokkaido and Korea inhabited Far East Asia before colonization by the present continental lineages of L. timidus, including the Sakhalin population. We estimated the time of the most recent common ancestor of the Hokkaido population to be 0.17 Mya, and found two distinct haplogroups within the island. One group had greater genetic diversity (mean number of pairwise differences: π = 0.0188 ± 0.0108) and appears to have expanded from the west to the entire island of Hokkaido. The other had lower genetic diversity (π = 0.0038 ± 0.0037) and its distribution was concentrated in the east. These contrasting west/east trends indicate that the Hokkaido population was fragmented in the past, and then subsequently expanded. Our study suggests that Hokkaido was an important refugium for boreal species in the far eastern region, and allowed the formation of various population genetic structures within the island.

Assessing the genetic landscape of a contact zone: the case of European hare in northeastern Greece

The European hare populations of the Balkan Peninsula comprise two divergent phylogenetic lineages with discrete geographical distribution slightly overlapping in the area of northeastern Greece and Bulgaria. Here we elucidate their contact zone, by defining the spatial distributional pattern of the two highly divergent groups, detecting individuals of hybrid origin, and identifying genetic barriers present in the area of their co-existence. Specimens from northeastern Greece were assayed for lineage assignment and population genetic inference based on a 511 bp fragment of mitochondrial DNA control region and allelic data from 10 microsatellite loci. Bayesian analyses on original and simulated genotypes were performed allowing for the contact zone delineation. Our results indicate high genetic diversity in both nuclear and mitochondrial DNA, strong population structure and non random spatial distribution of the differentiated gene pools. The information provided by the two types of molecular markers yielded consistent results. This study comprises a fine scale analysis of the contact zone between the two evolutionary lineages of European brown hares in northeastern Greece. Specific questions on the spatial patterns where addressed for the first time. Furthermore, hypotheses regarding the presence of hybrids were also tested. As a result, interpretive power to the diversity patterns observed today in the Balkans was added and previously overlooked aspects of the species biology were highlighted.

Past, present and future distributions of an Iberian Endemic, Lepus granatensis: ecological and evolutionary clues from species distribution models

The application of species distribution models (SDMs) in ecology and conservation biology is increasing and assuming an important role, mainly because they can be used to hindcast past and predict current and future species distributions. However, the accuracy of SDMs depends on the quality of the data and on appropriate theoretical frameworks. In this study, comprehensive data on the current distribution of the Iberian hare (Lepus granatensis) were used to i) determine the species' ecogeographical constraints, ii) hindcast a climatic model for the last glacial maximum (LGM), relating it to inferences derived from molecular studies, and iii) calibrate a model to assess the species future distribution trends (up to 2080). Our results showed that the climatic factor (in its pure effect and when it is combined with the land-cover factor) is the most important descriptor of the current distribution of the Iberian hare. In addition, the model's output was a reliable index of the local probability of species occurrence, which is a valuable tool to guide species management decisions and conservation planning. Climatic potential obtained for the LGM was combined with molecular data and the results suggest that several glacial refugia may have existed for the species within the major Iberian refugium. Finally, a high probability of occurrence of the Iberian hare in the current species range and a northward expansion were predicted for future. Given its current environmental envelope and evolutionary history, we discuss the macroecology of the Iberian hare and its sensitivity to climate change.

Identification of hare meat by a species-specific marker of mitochondrial origin

Meat species identification in food has gained increasing interest in recent years due to public health, economic and legal concerns. Following the consumer trend towards high quality products, game meat has earned much attention. The aim of the present work was to develop a DNA-based technique able to identify hare meat. Mitochondrial cytochrome b gene was used to design species-specific primers for hare detection. The new primers proved to be highly specific to Lepus species, allowing the detection of 0.01% of hare meat in pork meat by polymerase chain reaction (PCR). A real-time PCR assay with the new intercalating EvaGreen dye was further proposed as a specific and fast tool for hare identification with increased sensitivity (1pg) compared to end-point PCR (10pg). It can be concluded that the proposed new primers can be used by both species-specific end-point PCR or real-time PCR to accurately authenticate hare meat.

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.

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.