Genetic diversity within Anatolian brown hares (Lepus europaeus Pallas, 1778) and differentiation among Anatolian and European populations

Patterns of evolution in MHC class II DQA and DQB exon 2 genes of Alpine mountain hares, Lepus timidus varronis, and sympatric and parapatric brown hares, L. europaeus, from Switzerland

In natural populations, hybridization is known to occur between a wide range of species. However, its evolutionary significance is less clear. Genes involved in fighting pathogens are considered excellent candidates for studying adaptive introgression, although both introgression and balancing selection can generate similar patterns of diversity and differentiation. Here, we compared DQA and DQB MHC class II and microsatellite allelic diversity of sympatric and parapatric mountain (Lepus timidus) and brown hare (L. europaeus) populations from Switzerland. We detected higher genetic diversity in brown hares compared to mountain hares at both MHC and microsatellite loci. We consider the observed patterns of microsatellite diversity both for L. europaeus and L. timidus as result of stochastic demographic processes while the pattern of MHC polymorphism of the studied hare populations can be explained by pathogen-driven selection. Rare bidirectional gene flow between both hare species seems to occur specifically for MHC alleles. However, the high number of shared alleles showing similar high frequency in both species suggests that reciprocally exchanged MHC alleles are being maintained via balancing selection. Adaptation to similar pathogen communities can also lead to parallel selection of MHC alleles. Positive selection, recombination and mutations have played different roles in shaping the patterns of MHC allelic diversity in and differentiation between both species. Results for the latter evolutionary forces do not show a better matching between the sympatric populations compared to the parapatric ones, suggesting a minor role of introgression for the observed evolutionary patterns of the studied hare species.

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.

Gut microbiota of the European Brown Hare (Lepus europaeus)

Diseases of the gastrointestinal tract due to changes in the bacterial flora have been described with increasing incidence in the European brown hare. Despite extensive demographic and phylogeographic research, little is known about the composition of its gut microbiota and how it might vary based on potential environmental or host factors. We analysed the intestinal and faecal microbiota of 3 hare populations by Illumina MiSeq 16S rRNA gene amplicon sequencing. The phyla and OTU abundance composition differed significantly between intestinal and faecal samples (PERMANOVA: P = 0.002 and P = 0.031, respectively), but in both sample types Firmicutes and Bacteroidetes dominated the microbial community composition (45.51% and 19.30% relative abundance). Intestinal samples contained an enrichment of Proteobacteria compared with faecal samples (15.71-fold change, P < 0.001). At OTU level, a significant enrichment with best BLAST hits to the Escherichia-Shigella group, Eubacterium limosum, Sphingomonas kyeonggiensis, Flintibacter butyricus and Blautia faecis were detected in intestinal samples (P < 0.05). In our statistical model, geographic location and possibly associated environmental factors had a greater impact on the microbiota composition than host factors. Population had a significant effect on the composition of abundant intestinal and faecal OTUs, and on the abundance of potential pathogenic bacteria of the family Enterobacteriaceae, regularly associated with intestinal dysbiosis in hares, in faecal samples. Our study is the first to describe the microbiota in brown hares and provides a foundation to generate hypothesis aiming to test the role of gut health in population fluctuations of the species.

Islands and hybrid zones: combining the knowledge from "Natural Laboratories" to explain phylogeographic patterns of the European brown hare

Background

The aim of the study was to use hybrid populations as well as island populations of the European brown hare (Lepus europaeus) to explore the effect of evolutionary events, such as the post-deglaciation translocations, spontaneous and human-mediated, local adaptation and the genetic drift in the shaping of the phylogeographic patterns of the species. For this purpose, we used molecular markers, both nuclear and mitochondrial, that are indicative for local adaptation as well as neutral markers to elucidate the patterns of population differentiation based on geographic isolation and the clade of origin. To broaden our analysis, we included data from our previous studies concerning mainland populations, to explore the genetic differentiation in the base of the geographic origin (mainland/island) of the populations.

Results

Our results suggest that local adaptation shapes the differentiation in both genomes, favoring specific alleles in nuclear genes (e.g. DQA) or haplotypes in mtDNA (e.g. Control Region, CR). mtDNA variation was found to be in a higher level and was able to give a phylogeographic signal for the populations. Furthermore, the degree of variation was influenced not only by the geographic origin, but also by the clade of origin, since specific island populations of Anatolian origin showed a greater degree of variation compared to specific mainland populations of the European clade. Concerning the hybrid population, we confirmed the existence of both clades in the territory and we provided a possible explanation for the lack of introgression between the clades.

Conclusion

Our results indicate that the Quaternary’s climatic oscillations played a major role in the shaping of the phylogeographic patterns of the species, by isolating populations in the distinct refugia, where they adapted and differentiate in allopatry, leading to genome incompatibilities observed nowadays.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1354-y) contains supplementary material, which is available to authorized users.

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.

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.

Major histocompatibility complex variation at class II DQA locus in the brown hare (Lepus europaeus)

The major histocompatability complex (MHC) is a multigene family of receptors that bind and present antigenic peptides to T-cells. Genes of the MHC are characterized by an outstanding genetic polymorphism, which is considered to be maintained by positive selection. Sites involved in peptide binding form binding pockets (P) that are collectively termed the peptide-binding region (PBR). In this study, we examined the level of MHC genetic diversity within and among natural populations of brown hare (Lepus europaeus) from Europe and Anatolia choosing for analysis of the second exon of the DQA locus, one of the most polymorphic class II loci. We aimed at an integrated population genetic analysis of L. europeaus by (i) correlating MHC polymorphism to genetic variability and phylogenetic status estimated previously from maternally (mtDNA) and biparentally (allozymes, microsatellites) inherited loci; and (ii) comparing full-length exon amino acid polymorphism with functional polymorphism in the PBR and the binding pockets P1, P6 and P9. A substantial level of DQA exon 2 polymorphism was detected with two completely different set of alleles between the Anatolian and European populations. However, the phylogeny of full-length exon 2 Leeu-DQA alleles did not show a strong phylogeographic signal. The presence of balancing selection was supported by a statistically significant excess of nonsynonymous substitutions over synonymous in the PBR and a trans-species pattern of evolution detected after phylogenetic reconstruction. The differentiating patterns detected between genetic and functional polymorphism, i.e. the number and the distribution of pocket variants within and among populations, indicated a hierarchical action of selection pressures.

Evolutionary history of an MHC gene in two leporid species: characterisation of Mhc-DQA in the European brown hare and comparison with the European rabbit

We surveyed the genetic diversity of the expressed major histocompatibility complex class II DQA locus in natural populations of European brown hares, Lepus europaeus, from Austria and Belgium (267 individuals in total). Based on cDNA sequences, we designed hare-specific primers to amplify the highly variable second exon of the DQA gene. Using cloning-sequencing methodology and capillary electrophoresis single-strand conformation polymorphism, we found ten alleles of the DQA exon 2 locus across these two European regions, of which eight are described for the first time. To search for signals of selection and recombination in the evolution of the DQA gene within the leporids, we augmented our sample with orthologous DQA alleles from the European rabbit, Oryctolagus cuniculus, in order to carry out a species level, species pairwise comparison. We found evidence of recombination in the history of the DQA sequences in leporids with some recombinant alleles bridging the species divide. In both species, selection on peptide binding site codons can be detected, though stronger for the rabbit. This result suggests that there may be a differential selection pressure in the deeper evolutionary history of these two species due to differences in several demographic and ecological traits likely subjecting them to differential selection by parasites. Finally, evolutionary relationships show a widespread and statistically significant intermingling of alleles from the two species. The many macroparasites shared between hares and rabbits may explain this pattern of trans-species polymorphism.

Recruitment of mitochondrial tRNA genes as auxiliary variability markers for both intra- and inter-species analysis: The paradigm of brown hare (Lepus europaeus)

We sequenced and analyzed the mitochondrial tRNA(Thr) and tRNA(Pro) genes from brown hare (Lepus europaeus) individuals of different geographic distribution and we investigated the role of various nucleotide substitutions that were detected. We compared these tRNAs with the respective available mitochondrial tRNA genes sequences within Lepus species and among mammals. The mutations that were detected represent specific and conserved polymorphisms that do not seem to affect the structural and functional features that are required for participation of tRNA molecules in mitochondrial protein synthesis. These changes however, possibly reflect on the evolutionary background of the species, which is based on the high intra-genomic variability and the evolutionary dynamic of the mitochondrial DNA. In an attempt to compare the phylogeny that is based on these specific tRNA genes with the phylogeny that is produced from sequencing data of the mitochondrial variable loop, we came up with results that indicate similar phylogeographic clusters. This observation implies that the tRNA mutations that were used for the present study have been well tolerated during evolution and they define an additional genetic and biochemical tag that can be used for such studies. Based on this notion and according to our results, we propose that mitochondrial tRNA genes can be used as valuable auxiliary molecular markers for contemporaneous and linked biochemical and genetic analyses.

Multigenic and morphometric differentiation of ground squirrels (Spermophilus, Scuiridae, Rodentia) in Turkey, with a description of a new species

This study reports on the molecular phylogenetics of ground squirrels, genus Spermophilus, in Turkey using cytochrome b (1140bp), part of the D-loop and flanking tRNAs (572bp), X chromosome (867-1051bp) and Y chromosome (983-989bp) DNA sequences. Individuals also were characterized by karyotype and with geometric morphometric analyses of mandibles and skulls. Two hundred fourteen individuals from 91 localities were studied. All the data support the recognition of a new species in SW Anatolia: the Taurus ground squirrel Spermophilus taurensis sp. nov. The new species has a small distribution in the Taurus Mountains in an area that is a hotspot for biodiversity. Molecular clock analysis suggests that the new species diverged from the European ground squirrel, Spermophilus citellus, about 2.5 million years ago and that the ancestor of these two species diverged from the widespread Anatolian ground squirrel, Spermophilus xanthoprymnus, about 5 million years ago. Morphometric differentiation in skull and mandible shape among the three species is incomplete, but statistically significant. S. xanthoprymnus is subdivided into five cytochrome b phylogroups and we use these data to infer the location of glacial refugia where the species lived during the last glacial maximum. This study illustrates the potential of combined molecular and morphometric studies to uncover new Anatolian species and to reconstruct their phylogeographic history. The new species is important for squirrel taxonomy and for understanding Eurasian mammal evolution.

Biochemical Genetic Relationships Among Tunisian Hares (Lepus sp.), South African Cape Hares (L. capensis), and European Brown Hares (L. europaeus)

Tunisian hares (n = 45), currently assigned to Lepus capensis, were assayed for allelic variation at 40 allozyme loci, and allele frequencies at 32 loci were directly compared with earlier data of South African cape hares (L. capensis, n = 9) and European brown hares (L. europaeus, n = 244) to reveal genetic relationships among them. European mountain hares (L. timidus, n = 200) were used for outgroup comparison. In the Tunisian hares 27.5% of the loci were polymorphic with 2-4 alleles. Among all alleles at polymorphic loci, 15.1% occurred exclusively in Tunisian hares, 5.7% exclusively in cape hares, and 7.5% exclusively in brown hares at low frequencies. Not a single locus showed alternately fixed alleles between the samples of the L. capensis/L. europaeus complex. Levels of absolute and relative genetic differentiation among the samples of the L. capensis/L. europaeus complex were low, relative to pairwise comparisons involving mountain hares. Diverse cluster analyses and multidimensional scaling of various pairwise genetic distance matrices concordantly grouped Tunisian hares with brown hares, and South African cape hares clustered only slightly farther apart, whereas mountain hares were distinctly separate. These results suggest regionally distinct phylogenetic units within an overall cohesive gene pool in the L. capensis/L. europaeus complex, supporting Petter's view that all North African hares belong to L. capensis except for one local population of savanna hares, and that cape hares and brown hares are conspecific.