Chromosome-scale genome assembly of Lepus oiostolus (Lepus, Leporidae).

Lepus oiostolus (L. oiostolus) is a species endemic to the Qinghai-Tibet Plateau. However, the absence of a reference genome limits genetic studies. Here, we reported a high-quality L. oiostolus genome assembly, with scaffolds anchored to 24 chromosomes and a total assembled length of 2.80 Gb (contig N50 = 64.25 Mb). Genomic annotation uncovered 22,295 protein-coding genes and identified 49.84% of the sequences as transposable elements. Long interspersed nuclear elements (LINEs) constitute a high proportion of the genome. Our study is at the first time to report the chromosome-scale genome for the species of the L. oiostolus. It provides a valuable genomic resource for future research on the evolution of the Leporidae.

HAT: de novo variant calling for highly accurate short-read and long-read sequencing data.

MOTIVATION: de novo variants (DNVs) are variants that are present in offspring but not in their parents. DNVs are both important for examining mutation rates as well as in the identification of disease-related variation. While efforts have been made to call DNVs, calling of DNVs is still challenging from parent-child sequenced trio data. We developed Hare And Tortoise (HAT) as an automated DNV detection workflow for highly accurate short-read and long-read sequencing data. Reliable detection of DNVs is important for human genomics and HAT addresses this need. RESULTS: HAT is a computational workflow that begins with aligned read data (i.e. CRAM or BAM) from a parent-child sequenced trio and outputs DNVs. HAT detects high-quality DNVs from Illumina short-read whole-exome sequencing, Illumina short-read whole-genome sequencing, and highly accurate PacBio HiFi long-read whole-genome sequencing data. The quality of these DNVs is high based on a series of quality metrics including number of DNVs per individual, percent of DNVs at CpG sites, and percent of DNVs phased to the paternal chromosome of origin. AVAILABILITY AND IMPLEMENTATION: https://github.com/TNTurnerLab/HAT.

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.

Characterization of the Tongue Worm, Linguatula serrata (Pentastomida), Identified from Hares (Lepus europaeus) in Romania.

Linguatula serrata (Frölich, 1789) is a widespread parasite known as the tongue worm belonging to the family Linguatulidae. The adult form of the parasite is usually located in the upper respiratory tract of domestic and wild carnivores while the larval forms are located in the visceral organs of intermediate hosts (various herbivorous mammals). Twenty-four European brown hares (Lepus europaeus) were examined in this study, of which two were positive with L. serrata nymphs. The collected nymphs were examined morphologically using electron-microscopic analysis and molecularly by amplification of 18S rRNA and COX1 genes. Lung tissue samples were also collected and histopathological examination was performed. Histopathological examination revealed the following lesions: generalized inflammatory oedema, granulomas with necrosis, calcification and fibrosis in the bronchial tree. The results of molecular sequencing for L. serrata specimens collected from the European brown hares are deposited in GenBank. This study presents the first report on Linguatula serrata nymphs collected from L. europaeus in Romania, using molecular and morphological characterization simultaneously.

Phylogeography of Brucella suis biovar 2 with focus on Slovenian wildlife.

Brucella suis commonly infects swine but occasionally also other animal species and humans. Wild boars are the most important reservoir of B. suis biovar 2, continually infecting susceptible hosts through close contact. Nevertheless, the genetic diversity of B. suis in wildlife remains understudied. Here, we typed 17 Slovenian B. suis biovar 2 isolates obtained in 2017-2019 from wild boars (n = 16) and a hare (n = 1) using whole-genome sequencing (WGS). To assess the global phylogenetic diversity of B. suis, we compared them to 126 publicly available B. suis genomes. All Slovenian isolates fell within the biovar 2 lineage, confirming the previous multiplex PCR typing results. According to MLST-21, the wild boar isolates were of sequence types (STs) ST16 (n = 8) and ST153 (n = 8); the maximum genetic distance between isolates of the same ST was 28 wgMLST alleles. The ST153 isolates were restricted to the Slovenian-Croatian border and clustered together with the Croatian ST153 isolates from swine, indicating cross-border transmission of B. suis ST153 strain. The hare isolate was of ST40 and was genetically distant (≥ 489 alleles) from the wild boar isolates. The genome-wide phylogeny clearly separated different B. suis biovars. The present study is the first report on the population structure of B. suis in wildlife in Slovenia and shows that the Slovenian B. suis population is genetically heterogeneous. At the species level, B. suis biovars are clearly separated in the WGS-based phylogenetic tree and can therefore be reliably predicted using WGS.

Genome-wide association study of obstructive sleep apnoea in the Million Veteran Program uncovers genetic heterogeneity by sex.

BACKGROUND: Genome-wide association studies (GWAS) for obstructive sleep apnoea (OSA) are limited due to the underdiagnosis of OSA, leading to misclassification of OSA, which consequently reduces statistical power. We performed a GWAS of OSA in the Million Veteran Program (MVP) of the U.S. Department of Veterans Affairs (VA) healthcare system, where OSA prevalence is close to its true population prevalence. METHODS: We performed GWAS of 568,576 MVP participants, stratified by biological sex and by harmonized race/ethnicity and genetic ancestry (HARE) groups of White, Black, Hispanic, and Asian individuals. We considered both BMI adjusted (BMI-adj) and unadjusted (BMI-unadj) models. We replicated associations in independent datasets, and analysed the heterogeneity of OSA genetic associations across HARE and sex groups. We finally performed a larger meta-analysis GWAS of MVP, FinnGen, and the MGB Biobank, totalling 916,696 individuals. FINDINGS: MVP participants are 91% male. OSA prevalence is 21%. In MVP there were 18 and 6 genome-wide significant loci in BMI-unadj and BMI-adj analyses, respectively, corresponding to 21 association regions. Of these, 17 were not previously reported in association with OSA, and 13 replicated in FinnGen (False Discovery Rate p-value < 0.05). There were widespread significant differences in genetic effects between men and women, but less so across HARE groups. Meta-analysis of MVP, FinnGen, and MGB biobank revealed 17 additional, previously unreported, genome-wide significant regions. INTERPRETATION: Sex differences in genetic associations with OSA are widespread, likely associated with multiple OSA risk factors. OSA shares genetic underpinnings with several sleep phenotypes, suggesting shared aetiology and causal pathways. FUNDING: Described in acknowledgements.

Ancient DNA of the Don-Hares Assumes the Existence of Two Distinct Mitochondrial Clades in Northeast Asia.

Paleoclimatic changes during the Pleistocene-Holocene transition is suggested as a main factor that led to species extinction, including the woolly mammoth (Mammuthus primigenius), Steller's sea cow (Hydrodamalis gigas) and the Don-hare (Lepus tanaiticus). These species inhabited the territory of Eurasia during the Holocene, but eventually went extinct. The Don-hare is an extinct species of the genus Lepus (Leporidae, Lagomorpha), which lived in the Late Pleistocene-Early Holocene in Eastern Europe and Northern Asia. For a long time, the Don-hare was considered a separate species, but at the same time, its species status was disputed, taking into account both morphological data and mitochondrial DNA. In this study, mitochondrial genomes of five Don-hares, whose remains were found on the territory of Northeastern Eurasia were reconstructed. Firstly, we confirm the phylogenetic proximity of the "young" specimens of Don-hare and mountain or white hare, and secondly, that samples older than 39 Kya form a completely distinct mitochondrial clade.

The evolution of white-tailed jackrabbit camouflage in response to past and future seasonal climates.

The genetic basis of adaptive traits has rarely been used to predict future vulnerability of populations to climate change. We show that light versus dark seasonal pelage in white-tailed jackrabbits (Lepus townsendii) tracks snow cover and is primarily determined by genetic variation at endothelin receptor type B (EDNRB), corin serine peptidase (CORIN), and agouti signaling protein (ASIP). Winter color variation was associated with deeply divergent alleles at these genes, reflecting selection on both ancestral and introgressed variation. Forecasted reductions in snow cover are likely to induce widespread camouflage mismatch. However, simulated populations with variation for darker winter pelage are predicted to adapt rapidly, providing a trait-based genetic framework to facilitate evolutionary rescue. These discoveries demonstrate how the genetic basis of climate change adaptation can inform conservation.

Molecular phenotyping uncovers differences in basic housekeeping functions among closely related species of hares (Lepus spp., Lagomorpha: Leporidae).

Speciation is a fundamental evolutionary process, which results in genetic differentiation of populations and manifests as discrete morphological, physiological and behavioural differences. Each species has travelled its own evolutionary trajectory, influenced by random drift and driven by various types of natural selection, making the association of genetic differences between the species with the phenotypic differences extremely complex to dissect. In the present study, we have used an in vitro model to analyse in depth the genetic and gene regulation differences between fibroblasts of two closely related mammals, the arctic/subarctic mountain hare (Lepus timidus Linnaeus) and the temperate steppe-climate adapted brown hare (Lepus europaeus Pallas). We discovered the existence of a species-specific expression pattern of 1623 genes, manifesting in differences in cell growth, cell cycle control, respiration, and metabolism. Interspecific differences in the housekeeping functions of fibroblast cells suggest that speciation acts on fundamental cellular processes, even in these two interfertile species. Our results help to understand the molecular constituents of a species difference on a cellular level, which could contribute to the maintenance of the species boundary.

Francisella tularensis PCR detection in Cape hares (Lepus capensis) and wild rabbits (Oryctolagus cuniculus) in Algeria.

Tularemia is a zoonosis caused by the bacterium Francisella tularensis. Leporids are primary sources of human infections in the northern hemisphere. Africa is classically considered free of tularemia, but recent data indicate that this dogma might be wrong. We assessed the presence of this disease in wild leporids in Algeria. Between 2014 and 2018, we collected 74 leporids carcasses from spontaneously dead or hunted animals. Francisella tularensis DNA was detected by specific real-time PCR tests in 7/36 (19.44%) Cape hares (Lepus capensis) and 5/38 (13.15%) wild rabbits (Oryctolagus cuniculus). Known tularemia arthropod vectors infested half of the PCR-positive animals. At necropsy, F. tularensis-infected animals presented with an enlarged spleen (n = 12), enlarged adrenal glands (12), liver discoloration (12), hemorrhages (11), and pneumonia (11). Immunohistological examination of liver tissue from one animal was compatible with the presence of F. tularensis. Our study demonstrates the existence of tularemia in lagomorphs in Algeria. It should encourage investigations to detect this disease among the human population of this country.

Mountain- and brown hare genetic polymorphisms to survey local adaptations and conservation status of the heath hare (Lepus timidus sylvaticus, Nilsson 1831).

We provide the first whole genome sequences from three specimens of the mountain hare subspecies the heath hare (Lepus timidus sylvaticus), along with samples from two mountain hares (Lepus timidus timidus) and two brown hares (Lepus europaeus) from Sweden. The heath hare has a unique grey winter pelage as compared to other mountain hares (white) and brown hares (mostly brown), and face regional extinction, likely due to competitive exclusion from the non-native brown hare. Whole genome resequencing from the seven hare specimens were mapped to the Lepus timidus pseudoreference genome and used for detection of 11,363,883 polymorphic nucleotide positions. The data presented here could be useful for addressing local adaptations and conservation status of mountain hares and brown hares in Sweden, including unique subspecies.

Genetic and morphological identification of filarial worm from Iberian hare in Portugal.

The Iberian hare (Lepus granatensis) is an endemic species of the Iberian Peninsula and the only hare species found in Portugal, although also being present in some areas of Spain. The reduction of wild hare populations due to several ecological and sanitary factors, has been raising growing concerns in the recent years. Despite different helminth species were already described in Iberian hares in Portugal, to this date, no filarial worms have been identified in this species. Furthermore, only a few studies on lagomorphs' onchocercid worms are available, referring to other hosts species of hares and/or rabbits. In this study, we describe the presence of filarial worms in the blood vessels of two adult Iberian hares collected in 2019 in continental Portugal. Morphology and sequencing data from the 12S rRNA, coxI, 18S rRNA, myoHC, hsp70 and rbp1 genes, showed that the filaroid species were genetically related with Micipsella numidica. However, the extension of the genetic differences found with M. numidica suggests that the filaroids specimens under study belong to a new species, that we provisionally named Micipsella iberica n. sp.. The body location of this putative new parasite species and its physiological implications indicate that it may constitute a potential menace to the already fragile Iberian hare justifying, therefore, further investigation regarding the morphological characterization, prevalence and real clinical impact of this new parasite in hares.

As fast as a hare: Did intraspecific morphological change bring the Hallands Väderö Island population of Lepus timidus close to interspecific differences in less than 150 years?

The study of insular variation has fascinated generations of biologists and has been central to evolutionary biology at least since the time of Wallace and Darwin. In this context, using 3D geometric morphometrics, I investigate whether the population of mountain hares (Lepus timidus Linnaeus, 1758) introduced in 1857 on the Swedish island of Hallands Väderö shows distinctive traits in cranial size and shape. I find that size divergence follows the island rule, but is very small. In contrast, shape differences, compared to the mainland population, are almost as large as interspecific differences among lineages separated by hundreds of thousands of years of a largely independent evolutionary history. Even if, contrary to what is documented in the scientific literature, mountain hares were present in HV before 1857, the evolutionary history of this population could not have start earlier than the end of the last glaciation (i.e., at least one order of magnitude more recently than the separation of L. timidus from other hare species in this study). My results, thus, suggest that the insular population is a significant evolutionary unit and a potentially important component of the diversity of Swedish mountain hares. This is interesting for evolutionary biologists, but even more relevant for conservationists trying to protect the disappearing population of southern Swedish L. timidus, threatened by changes in climate and the environment, as well as by disease and the introduced European hare (Lepus europaeus Pallas, 1778). Island populations of mountain hares, thus, represent a potential source for future reintroductions on the mainland and, as my research shows, an important component of variability to maximize the preservation of the evolutionary potential in a species facing huge environmental changes.

The evolutionary pathways for local adaptation in mountain hares.

Understanding the evolution of local adaptations is a central aim of evolutionary biology and key for the identification of unique populations and lineages of conservation relevance. By combining RAD sequencing and whole-genome sequencing, we identify genetic signatures of local adaptation in mountain hares (Lepus timidus) from isolated and distinctive habitats of its wide distribution: Ireland, the Alps and Fennoscandia. Demographic modelling suggested that the split of these mountain hares occurred around 20 thousand years ago, providing the opportunity to study adaptive evolution over a short timescale. Using genome-wide scans, we identified signatures of extreme differentiation among hares from distinct geographic areas that overlap with area-specific selective sweeps, suggesting targets for local adaptation. Several identified candidate genes are associated with traits related to the uniqueness of the different environments inhabited by the three groups of mountain hares, including coat colour, ability to live at high altitudes and variation in body size. In Irish mountain hares, a variant of ASIP, a gene previously implicated in introgression-driven winter coat colour variation in mountain and snowshoe hares (L. americanus), may underlie brown winter coats, reinforcing the repeated nature of evolution at ASIP moulding adaptive seasonal colouration. Comparative genomic analyses across several hare species suggested that mountain hares' adaptive variants appear predominantly species-specific. However, using coalescent simulations, we also show instances where the candidate adaptive variants have been introduced via introgressive hybridization. Our study shows that standing adaptive variation, including that introgressed from other species, was a crucial component of the post-glacial dynamics of species.

Hybridization with mountain hares increases the functional allelic repertoire in brown hares.

Brown hares (Lepus europaeus Pallas) are able to hybridize with mountain hares (L. timidus Linnaeus) and produce fertile offspring, which results in cross-species gene flow. However, not much is known about the functional significance of this genetic introgression. Using targeted sequencing of candidate loci combined with mtDNA genotyping, we found the ancestral genetic diversity in the Finnish brown hare to be small, likely due to founder effect and range expansion, while gene flow from mountain hares constitutes an important source of functional genetic variability. Some of this variability, such as the alleles of the mountain hare thermogenin (uncoupling protein 1, UCP1), might have adaptive advantage for brown hares, whereas immunity-related MHC alleles are reciprocally exchanged and maintained via balancing selection. Our study offers a rare example where an expanding species can increase its allelic variability through hybridization with a congeneric native species, offering a route to shortcut evolutionary adaptation to the local environmental conditions.

Positive selection on two mitochondrial coding genes and adaptation signals in hares (genus Lepus) from China.

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.

The Legacy of Recurrent Introgression during the Radiation of Hares.

Hybridization may often be an important source of adaptive variation, but the extent and long-term impacts of introgression have seldom been evaluated in the phylogenetic context of a radiation. Hares (Lepus) represent a widespread mammalian radiation of 32 extant species characterized by striking ecological adaptations and recurrent admixture. To understand the relevance of introgressive hybridization during the diversification of Lepus, we analyzed whole exome sequences (61.7 Mb) from 15 species of hares (1-4 individuals per species), spanning the global distribution of the genus, and two outgroups. We used a coalescent framework to infer species relationships and divergence times, despite extensive genealogical discordance. We found high levels of allele sharing among species and show that this reflects extensive incomplete lineage sorting and temporally layered hybridization. Our results revealed recurrent introgression at all stages along the Lepus radiation, including recent gene flow between extant species since the last glacial maximum but also pervasive ancient introgression occurring since near the origin of the hare lineages. We show that ancient hybridization between northern hemisphere species has resulted in shared variation of potential adaptive relevance to highly seasonal environments, including genes involved in circadian rhythm regulation, pigmentation, and thermoregulation. Our results illustrate how the genetic legacy of ancestral hybridization may persist across a radiation, leaving a long-lasting signature of shared genetic variation that may contribute to adaptation. [Adaptation; ancient introgression; hybridization; Lepus; phylogenomics.].

The Origin and Spread of Locally Adaptive Seasonal Camouflage in Snowshoe Hares.

AbstractAdaptation is central to population persistence in the face of environmental change, yet we seldom precisely understand the origin and spread of adaptive variation in natural populations. Snowshoe hares (Lepus americanus) along the Pacific Northwest coast have evolved brown winter camouflage through positive selection on recessive variation at the Agouti pigmentation gene introgressed from black-tailed jackrabbits (Lepus californicus). Here, we combine new and published whole-genome and exome sequences with targeted genotyping of Agouti to investigate the evolutionary history of local seasonal camouflage adaptation in the Pacific Northwest. We find evidence of significantly elevated inbreeding and mutational load in coastal winter-brown hares, consistent with a recent range expansion into temperate coastal environments that incurred indirect fitness costs. The genome-wide distribution of introgression tract lengths supports a pulse of hybridization near the end of the last glacial maximum, which may have facilitated range expansion via introgression of winter-brown camouflage variation. However, signatures of a selective sweep at Agouti indicate a much more recent spread of winter-brown camouflage. Through simulations, we show that the delay between the hybrid origin and subsequent selective sweep of the recessive winter-brown allele can be largely attributed to the limits of natural selection imposed by simple allelic dominance. We argue that while hybridization during periods of environmental change may provide a critical reservoir of adaptive variation at range edges, the probability and pace of local adaptation will strongly depend on population demography and the genetic architecture of introgressed variation.

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.

Purifying selection shaping the evolution of the Toll-like receptor 2 TIR domain in brown hares (Lepus europaeus) from Europe and the Middle East.

Toll-like receptors (TLRs) are transmembrane proteins of the innate immune system, composed of the ectodomain involved in pathogen recognition and the intracellular Toll/interleukin-1 receptor (TIR) domain important for downstream signal transduction. Here, we analyze the genetic variability of TIR nucleotide and amino-acid sequences of the TLR2 gene in 243 brown hares from Europe and the Middle East and tested for the presence of selection signals and spatial structuring. TLR2 TIR domain sequences were PCR amplified and sequenced, while genotyping was performed by phasing. Genetic diversity indices were calculated in DnaSP and Arlequin, while presence of selection signals was tested using MEGA and the Datamonkey web server. The presence of spatial patterns in TIR sequence distribution was tested by spatial Principal Component Analysis (sPCA) in adegenet. A total of 13 haplotypes were revealed with haplotype diversity of 0.424, and nucleotide diversity (π) of 0.00138. Two spatial clusters were revealed: "Anatolia/Middle East" and "Europe". In Anatolia the two most prevalent amino-acid variants, A and B (the latter being the most ancestral) were maintained at similar frequencies; but in Europe a shift in genotype frequencies was observed as well as a higher number of nonsynonymous substitutions giving rise to novel amino-acid protein variants originating from the evolutionarily younger protein variant. Molecular diversity (haplotype and nucleotide diversity) indices were significantly higher in the "Anatolia/Middle East" cluster. A signal of purifying selection was detected acting on the TIR sequences.