Search

Background Animal mitochondria play a central role in energy production in the cells through the oxidative phosphorylation (OXPHOS) pathway. Recent studies of selection on different mitochondrial OXPHOS genes have revealed the adaptive implications of amino acid changes in these subunits. In hares, climatic variation and/or introgression were suggested to be at the origin of such adaptation. Here we looked for evidence of positive selection in three mitochondrial OXPHOS genes, using tests of selection, protein structure modelling and effects of amino acid substitutions on the protein function and stability. We also used statistical models to test for climate and introgression effects on sites under positive selection. Results Our results revealed seven sites under positive selection in ND4 and three sites in Cytb. However, no sites under positive selection were observed in the COX1 gene. All three subunits presented a high number of codons under negative selection. Sites under positive selection were mapped on the tridimensional structure of the predicted models for the respective mitochondrial subunit. Of the ten amino acid replacements inferred to have evolved under positive selection for both subunits, six were located in the transmembrane domain. On the other hand, three codons were identified as sites lining proton translocation channels. Furthermore, four codons were identified as destabilizing with a significant variation of Δ vibrational entropy energy between wild and mutant type. Moreover, our PROVEAN analysis suggested that among all positively selected sites two fixed amino acid replacements altered the protein functioning. Our statistical models indicated significant effects of climate on the presence of ND4 and Cytb protein variants, but no effect by trans-specific mitochondrial DNA introgression, which is not uncommon in a number of hare species. Conclusions Positive selection was observed in several codons in two OXPHOS genes. We found that substitutions in the positively selected codons have structural and functional impacts on the encoded proteins. Our results are concordantly suggesting that adaptations have strongly affected the evolution of mtDNA of hare species with potential effects on the protein function. Environmental/climatic changes appear to be a major trigger of this adaptation, whereas trans-specific introgressive hybridization seems to play no major role for the occurrence of protein variants.

In natural populations, allelic diversity of the major histocompatibility complex (MHC) is commonly interpreted as resulting from positive selection in varying spatiotemporal pathogenic landscapes. Composite pathogenic landscape data are, however, rarely available. We studied the spatial distribution of allelic diversity at two MHC class II loci (DQA, DQB) in hares, Lepus capensis, along a steep ecological gradient in North Africa and tested the role of climatic parameters for the spatial distribution of DQA and DQB proteins. Climatic parameters were considered to reflect to some extent pathogenic landscape variation. We investigated historical and contemporary forces that have shaped the variability at both genes, and tested for differential selective pressure across the ecological gradient by comparing allelic variation at MHC and neutral loci. We found positive selection on both MHC loci and significantly decreasing diversity from North to South Tunisia. Our multinomial linear models revealed significant effects of geographical positions that were correlated with mean annual temperature and precipitation on the occurrence of protein variants, but no effects of co-occurring DQA or DQB proteins, respectively. Diversifying selection, recombination, adaptation to local pathogenic landscapes (supposedly reflected by climate parameters) and neutral demographic processes have shaped the observed MHC diversity and differentiation patterns.

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.

Toll-like receptors (TLRs) are a major group of proteins that recognize molecular components of infectious agents, known as pathogen associated molecular patterns (PAMPs). The structure of these genes is similar and characterized by the presence of an ectodomain, a signal transmembrane segment and a highly conserved cytoplasmic domain. The latter domain is homologous to the human interleukin-1 receptor (IL1R) and human IL-18 receptor (IL-18R) and designated TIR domain. The latter domain of the TLR genes was suggested to be very conservative and its evolution is driven by purifying selection. Variability and evolution of the TIR sequences of TLR2 gene were studied in three hare populations from Tunisia with different ecological characteristics (NT-North Tunisia with Mediterranean, CT-Central Tunisia with semi-arid, and ST-South Tunisia with arid climate). Sequencing of a 372bp fragment of TIR2 revealed 25 alleles among 110 hares. Twenty variable nucleotide positions were detected, of which 7 were non-synonymous. The highest variability was observed in CT, with 16 polymorphic positions. In ST, only 4 polymorphic nucleotide positions were detected with all diversity values lower than those recorded for the other two populations. By using several approaches, no positive selection was detected. However, evidence of purifying selection was found at two positions. The logistic models of the most common TIR2 protein variant that we run to examine whether its occurrence was affected by climatic variation independent of the geographic sample location suggested only a longitudinal effect. Finally, the mapping of the non-synonymous mutations to the inferred tertiary protein structure showed that they were all localized in the different loop regions. Among all non-synonymous substitutions, three were suggested to be deleterious as evidenced by PROVEAN analysis. The observed patterns of variability characterized by low genetic diversity in ST might suggest that the TIR region was more affected, than other markers, by genetic drift or/and that these patterns were shaped by different selective pressures under different ecological conditions. Notably, this low diversity was not detected by other (putatively neutral) microsatellite markers analysed in the course of other studies. But low diversity was also found for two MHC class II adaptive immune genes. As expected from functionally important regions, the evolution of the TIR2 domain is mainly driven by purifying selection. However, the occurrence of deleterious non-synonymous substitutions might highlight the flexible evolution of the TIR genes and/or their interactions with other proteins.

In Bulgaria, chamois (Rupicapra rupicapra balcanica) are protected by law and it has been recommended to assess their genetic diversity, level of inbreeding, and extent of introgression by Alpine chamois (R. r. rupicapra) that were released in the last century in the Rhodope Mountains. Chamois (n = 81) collected between 2009 and 2013 at the “Izvora State Hunting Reserve” (ISHR, western Rhodopes, a major source for chamois translocations in Bulgaria; census size: 200–250) were screened for allelic variability at 12 microsatellite loci and compared to chamois (R. r. rupicapra, n = 135) from six populations in the Austrian Alps. Indices of genetic variability were significantly lower for the ISHR population than for the Austrian populations. A significant bottleneck signal was observed for the ISHR population, but no distinct inbreeding signal due to non-random mating. F-statistics, absolute genetic distances, Bayesian STRUCTURE and assignment analyses, and a factorial correspondence analysis concordantly indicated distinct differentiation between the ISHR and Austrian chamois. STRUCTURE identified some few ISHR chamois as possibly introgressed by Alpine chamois. ONeSAMP revealed effective population size estimates of Ne ≈ 50 for the ISHR population. A reduction of shooting quotas, a reinforcement of poaching control, and exchanging chamois with other Bulgarian populations should increase its Ne to prevent further loss of genetic variability. Screening for introgression should be extended to other Bulgarian populations, particularly to areas where Alpine chamois were released in the past; and it should accompany all translocations, to prevent possible dissemination of introgressed chamois.

Background Recent studies of selection on mitochondrial (mt) OXPHOS genes suggest adaptation due mainly to environmental variation. In this context, Tunisian hares that display several external phenotypes with phylogenetically rather homogenous gene pool and shallow population structure provide a good precondition to detect positive selection on mt genes related to environmental/climatic variation, specifically ambient temperature and precipitation. Results We used codon-based methods along with population genetic data to test for positive selection on ATPase synthase 6 (ATP6) and NADH dehydrogenase 2 (ND2) of cape hares (Lepus capensis) collected along a steep ecological gradient in Tunisia. We found significantly higher differentiation at the ATP6 locus across Tunisia, with sub-humid Mediterranean, semi-arid, and arid Sahara climate than for fourteen unlinked supposedly neutrally evolving nuclear microsatellites and mt control region sequences. This suggested positive selection on ATP6 sequences, which was confirmed by several codon-based tests for one sequence site that together with a second site translated into four different amino acids. Positive selection on ND2 sequences was also confirmed by several codon-based tests. The corresponding frequencies of the two most prevalent variants at each locus varied significantly across climate regions, and our logistic general linear models of occurrence of those proteins indicated significant effects of mean annual temperature for ATP6 and mean minimum temperature of the coldest month of the year for ND2, independent of geographical location, annual precipitation, and the respective co-occurring protein at the second locus. Moreover, presence of the ancestral ATP6 protein, as inferred from phylogenetic networks, was positively affected by the simultaneous presence of the derived ND2 protein and vice versa, independent of temperature, precipitation, or geographic location. Finally, we obtained a significant coevolution signal for the ancestral ATP6 and derived ND2 sequences and vice versa. Conclusions positive selection was strongly suggested by the population genetic approach and the codon-based tests in both mtDNA genes. Moreover, the two most prevalent proteins at the ATP6 locus were distributed at significantly varying frequencies across the study area with a significant effect of mean annual temperature on the occurrence of the ATP6 proteins independent of geographical coordinates and the co-occuring ND2 protein variant. For ND2, occurrence of the two most frequent protein variants was significantly influenced by the mean minimum temperature of the coldest month, independent of the co-occurring ATP6 protein variant and geographical coordinates. This strongly suggests direct involvement of ambient temperature in the adaptation of the studied mtOXPHOS genes. Electronic supplementary material The online version of this article (doi:10.1186/s12862-017-0896-0) contains supplementary material, which is available to authorized users.

The European ground squirrel (Spermophilus citellus) is endangered and in decline. Populations are increasingly fragmented, and only a coordinated conservation effort at the European level may guarantee its long-term survival. To obtain a general population genetic picture on a larger geographic scale, we screened 117 individuals from seven local populations in Hungary, Romania, and Austria for allelic variation at eleven microsatellite loci. We found a high (23.4%) proportion of private alleles, and a moderate to somewhat elevated level (15.27%) of partitioning of genetic diversity among populations, compared to that found in many other terrestrial mammals. Genetic variability was significantly higher than in earlier studied Czech populations that are considered genetically depleted, but significantly lower than in undisturbed populations of S. suslicus and S. brunneus, that are similar to the European ground squirrel in their ecological requirements, reproductive biology, and social organization. Genetic diversity was also lower than in most presumably “undisturbed” populations of other Sciurid species. This, together with the observed level and pattern of genetic differentiation among populations, such as no significant increase of genetic differentiation with geographic distance and similar variance of genetic differentiation between populations independent of geographic distance, indicated the prevalence of relatively strong drift effects for all populations. A Bayesian STRUCTURE analysis and a factorial correspondence analysis concordantly revealed a fairly complex genetic composition of local populations, but no major geographic trend in the pattern of the genetic composition. Overall, the results suggest disintegration of local colonies that might earlier have been more connected genetically. The STRUCTURE analysis also suggested anthropogenic translocations among single Hungarian populations. Our data on genetic diversity and its distribution do not object to such conservation measures. Translocation of individuals particularly from nearby populations may increase the chances of survival of small and isolated populations and counteract inbreeding at low densities.

Conventionally, Lepus capensis is considered to range across large parts of Africa, the Middle East, Central and Far East Asia. However, a recent morphological study restricts cape hares tentatively to a small range in the Western Cape Region of South Africa and groups all other L. capensis-type hares from South Africa into a new species: L. centralis. Here, we studied molecular relationships among L. capensis-type hares from South Africa. Phenotypically and morphologically the individuals matched either the newly described L. capensis or L. centralis. We examined 66 hares for allelic variation at 13 microsatellite loci and for sequence variation of the hypervariable domain 1 of the mitochondrial control region. All tree-generating analyses of the currently obtained sequences and all South African cape hare sequences downloaded from GenBank revealed monophyly when compared to sequences of various other Lepus species. A network analysis indicated close evolutionary relationships between hares of the “L. capensis-phenotype” and the “L. centralis-phenotype” (according to Palacios et al. 2008) from the southwest of the Western Cape, relative to their pronounced evolutionary divergence from all other more central, northern, and north-eastern L. capensis-type hares. F-statistics, a Bayesian admixture STRUCTURE model, as well as a principal coordinate analysis of microsatellite data indicated close genetic relationships among all South African L. capensis-type hares studied presently. A coalescence model-based migration analysis for microsatellite alleles indicated gene flow between most of the considered subspecies of cape hare, including L. capensis capensis and L. capensis centralis, theoretically sufficient to balance stochastic drift effects. Concordantly, AMOVA models revealed only little effects of partitioning microsatellite variation into the two suggested morpho-species “L. capensis” and “L. centralis”. Under an “Interbreeding Species Concept” (e.g. a strict or relaxed Biological Species Concept), the current molecular data demonstrate conspecificity of the two proposed morpho-species “L. capensis” and “L. centralis”. Based on the present molecular data the differentiation of subspecies of cape hares from southern Africa is discussed.

Variability and phylogenetic relationships of sequences of the hypervariable domain I (HVI) of the mitochondrial DNA was studied in 46 brown hares (Lepus europaeus) from Anatolia, to test the hypotheses that (i) hares from several islands off the Anatolian coast and from Cyprus are phylogenetically close to mainland Anatolian hares, (ii) Anatolian hare sequence variability is higher than that of typical European brown hares, and to iii) infer possible Anatolian source populations of hares from some islands in the eastern Mediterranean. Neighbor joining and Maximum Parsimony analyses revealed reciprocal monophyly for sequences from Anatolia, the considered eastern Mediterranean islands off the Anatolian coast, Cyprus, and those sequences published earlier form NE Greece that were supposed to originate from earlier immigration via the late-Pleistocene/early-Holocene land bridge that connected SE Europe and W Anatolia (Kasapidis et al., 2005. Mol. Phylogenet. Evol. 34, 55–66). A high sequence idiosyncrasy was found among the Anatolian samples. Almost all approaches to compare variability between Anatolian and the downloaded European sequence data indicated higher sequence diversity in Anatolia, in accordance with earlier findings for allozyme loci. Network and principal coordinate analyses of the Anatolian sequences and those from the islands off the Anatolian coast as well as the Anatolian-type NE Greek sequences suggested high mitochondrial gene exchange among local populations in Anatolia with little effect of possible geographic barriers, and did not provide clues for tracing possible origins of island populations.

Phylogenetic relationships and systematics of hares and jackrabbits (genus Lepus, Lagomorpha) are notoriously difficult. One intriguing issue is the relationship between cape hares (Lepus capensis L., 1758) and brown hares (Lepus europaeus Pallas, 1778). According to Petter, F., [1959. Elements d'une revision des lievres africains du sous-genre Lepus. Mammalia 23, 41–67; 1961. Elements d'une revision des lievres europeens et asiatique du sous-genre Lepus. Z. Saugetierkd. 26, 30–40], these two phenotypically different taxa are conspecific. On the contrary, recent arguments were put forward to split African capensis into several species. Eleven variable microsatellite loci were used to determine the spatial genetic structure of African and South Israeli hares, currently considered cape hares, and European, Anatolian, and North Israeli hares, currently considered brown hares. Patterns of genetic differentiation of 294 hares collected at 18 locations from both species ranges were examined for accordance to the hypothesis of conspecificity or to the concept of two distinct species, the cape hare and the brown hare. Genetic variation was high (14.7 alleles/locus) and the distribution of alleles and genotypes was not homogeneous across the locations. Results from FST, RST, and GST′-statistics including various AMOVA models, comparisons of private alleles with corresponding relative genetic differentiation between neighbouring populations, as well as absolute genetic differentiation patterns were by and large congruent. They all indicated that the degree of differentiation between populations belonging to the two different taxa was not higher than within the taxa. Partitioning of genetic diversity was congruent with a isolation-by-distance model, and this also suggested that the differentiation was due to geographic distances rather than the occurrence of two species-specific gene pools. Homoplasy likely had no or only a marginal effect on our results, as indicated by the standardized genetic distance GST′ and the presence of differently fixed alleles between the South African cape hare and all European brown hare populations at the highly polymorphic Sol08 locus. The results are concordant with earlier multilocus allozyme and mtDNA RFLP data, but contrast with mtDNA sequence data published earlier, that suggest distinctly separate phylogenetic units in Africa, the Near East and Europe. Nevertheless, as concluded from the microsatellite data brown hares and African cape hares might be connected via the Near East by gene flow along a network of intergrading populations that have regionally distinct gene pools.

Backround The melon fruit fly Bactrocera cucurbitae and the peach fruit fly Bactrocera zonata are serious pests, native to Asia, that have recently invaded Egypt. In Tunisia, no report of these species has yet been made, but pest risk data suggest that both Bactrocera species are likely to establish in other countries of the Mediterranean region. Results A quick method, based on a PCR-RFLP of the mitochondrial COI gene, has been developed to enhance species identification for quarantine purposes. The restriction enzymes DdeI and XmnI were used to digest COI PCR products, enabling a rapid diagnosis of B. cucurbitae and B. zonata and their distinction from Tunisia's most devastating fruit pest, Ceratitis capitata. Conclusion The simplicity and relatively low cost of this molecular approach will replace the need to rear immature stages through to adults for identification and will facilitate rapid quarantine decisions providing greater plant protection.

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.

Background In mammals, males typically have shorter lives than females. This difference is thought to be due to behavioural traits which enhance competitive abilities, and hence male reproductive success, but impair survival. Furthermore, in many species males usually show higher parasite burden than females. Consequently, the intensity of selection for genetic factors which reduce susceptibility to pathogens may differ between sexes. High variability at the major histocompatibility complex (MHC) genes is believed to be advantageous for detecting and combating the range of infectious agents present in the environment. Increased heterozygosity at these immune genes is expected to be important for individual longevity. However, whether males in natural populations benefit more from MHC heterozygosity than females has rarely been investigated. We investigated this question in a long-term study of free-living Alpine chamois (Rupicapra rupicapra), a polygynous mountain ungulate. Results Here we show that male chamois survive significantly (P = 0.022) longer if heterozygous at the MHC class II DRB locus, whereas females do not. Improved survival of males was not a result of heterozygote advantage per se, as background heterozygosity (estimated across twelve microsatellite loci) did not change significantly with age. Furthermore, reproductively active males depleted their body fat reserves earlier than females leading to significantly impaired survival rates in this sex (P < 0.008). This sex-difference was even more pronounced in areas affected by scabies, a severe parasitosis, as reproductively active males were less likely to survive than females. However, we did not find evidence for a survival advantage associated with specific MHC alleles in areas affected by scabies. Conclusions Increased MHC class II DRB heterozygosity with age in males, suggests that MHC heterozygous males survive longer than homozygotes. Reproductively active males appear to be less likely to survive than females most likely because of the energetic challenge of the winter rut, accompanied by earlier depletion of their body fat stores, and a generally higher parasite burden. This scenario renders the MHC-mediated immune response more important for males than for females, which implies a relatively stronger selection pressure on MHC genes in males than in females.

Among the European fauna, the Sardinian hare (Lepus sp.) is peculiar in that it differs from all other hares inhabiting the continent. Here, we report on the variation of a 461 bp sequence of hypervariable domain 1 of the mitochondrial control region, examined in 42 hares collected throughout Sardinia and compared to the corresponding sequences of different Lepus taxa. Seventeen novel haplotypes were found in the Sardinian population, resulting in a haplotype diversity of 0.840 and a nucleotide diversity of 0.012. As a result of Bayesian and principal coordinates analyses, Sardinian hares were grouped with North African hares, constituting a monophyletic clade that diverges from all other Old World hares, including Cape hares from South Africa and East Asia. Hence, our data agree that populations inhabiting North Africa and Sardinia form a distinct taxon, which could possibly be included in the L. capensis superspecies. Moreover, two corresponding lineages can be found in Sardinia and Tunisia, providing evidence of a common origin of the two populations and thus supporting the hypothesis that North African hares were introduced into the island in historical times. Our data show that the two lineages differ in their geographic distribution throughout the island and that the wild Sardinian population also shows the signature of a postintroduction demographic expansion.

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.