Phylogeography of the brown hare (Lepus europaeus) in Europe: a legacy of south‐eastern Mediterranean refugia?

Aim We analysed the population genetics of the brown hare (Lepus europaeus) in order to test the hypothesis that this species migrated into central Europe from a number of late glacial refugia, including some in Asia Minor. Location Thirty‐three localities in Greece, Bulgaria, Italy, Croatia, Serbia, Poland, Switzerland, Austria, France, Germany, the Netherlands, Spain, the United Kingdom, Turkey and Israel. Methods In total, 926 brown hares were analysed for mitochondrial DNA (mtDNA) variation by restriction fragment length polymorphism (RFLP) performed on polymerase chain reaction‐amplified products spanning cytochrome b (cyt b)/control region (CR), cytochrome oxidase I (COI) and 12S–16S rRNA. In addition, sequence analysis of the mtDNA CR‐I region was performed on 69 individuals, and the data were compared with 137 mtDNA CR‐I sequences retrieved from GenBank. Results The 112 haplotypes detected were partitioned into five phylogeographically well‐defined major haplogroups, namely the ‘south‐eastern European type haplogroup’ (SEEh), ‘Anatolian/Middle Eastern type haplogroup’ (AMh), ‘European type haplogroup, subgroup A’ (EUh‐A), ‘European type haplogroup, subgroup B’ (EUh‐B) and ‘Intermediate haplogroup’ (INTERh). Sequence data retrieved from GenBank were consistent with the haplogroups determined in this study. In Bulgaria and north‐eastern Greece numerous haplotypes of all five haplogroups were present, forming a large overlap zone. Main conclusions The mtDNA results allow us to infer post‐glacial colonization of large parts of Europe from a late glacial/early Holocene source population in the central or south‐central Balkans. The presence of Anatolian/Middle Eastern haplotypes in the large overlap zone in Bulgaria and north‐eastern Greece reveals gene flow from Anatolia to Europe across the late Pleistocene Bosporus land‐bridge. Although various restocking operations could be partly responsible for the presence of unexpected haplotypes in certain areas, we nevertheless trace a strong phylogeographic signal throughout all regions under study. Throughout Europe, mtDNA results indicate that brown hares are not separated into discernable phyletic groups.     

Genetic variation in populations of the endangered fish Ladigesocypris ghigii and its implications for conservation

1. The genetic variation of the endangered freshwater fish Ladigesocypris ghigii, endemic to the island of Rhodes (Greece), was investigated for nine populations, originating from seven different stream systems and a reservoir, both at the mtDNA and nuclear level, in order to suggest conservation actions. 2. Both restriction fragment length polymorphism analysis of five segments of mitochondrial DNA (ND‐5/6, COI and 12S‐16S rRNA) amplified by polymerase chain reaction, and random amplified polymorphic DNA analysis, revealed extremely low levels of intra‐population polymorphism. It is highly likely that the low intra‐population variability is the result of successive bottleneck events evident in shrinkage and expansion of the populations year after year, which may have led to a complete loss of several genotypes and haplotypes, and an increased degree of inbreeding. 3. Inter‐population genetic structuring was high, with fixation of haplotypes within six of the nine populations and fixation of alleles within populations originating from different waterbodies. It is probable that all haplotypes and/or alleles found were initially represented in all populations. However, because of the long time of isolation coupled with successive bottleneck and subsequent genetic drift, common mtDNA haplotypes and alleles among the populations may have become rare or extinct through stochastic lineage loss. 4. Although nucleotide divergence among haplotypes was very shallow, half of the haplotypes recorded (three of six), resulted from nucleotide changes on the 12S–16S rRNA segments, which are the most conserved part of the mitochondrial genome. This fact may indicate that the observed genetic variation did not necessarily result only from the retention of ancestral polymorphism, but may have arisen through mutation and complete lineage sorting over a relatively small number of generations, once the populations had become isolated from one another. 5. Our data suggest that two of the L. ghigii populations may be on independent evolutionary trajectories. Considering that each population appears so far well adapted within each site, all populations should be managed and conserved separately.     

Biochemical genetic variability in brown hares (Lepus europaeus) from Greece

Allozyme variability of 91 brown hares (Lepus europaeus) from seven regions in Greece was compared to existing data of Bulgarian populations to test the hypothesis of the occurrence of specific alleles in Greece, likely stemming from an isolated Late Pleistocene refugial population in the southern Balkans. This hypothesis is particularly suggested by some subfossil Late Pleistocene hare remains in Greece and the reported high mtDNA diversity in Greek hares. Allozymic diversity could be higher in Greek hares than in hares from neighboring regions as a result of the accumulation of variants in a long-lasting Pleistocene refugium. Conversely, Greek hares could exhibit reduced genetic diversity because of long-lasting low effective population sizes during the Late Glacial Maximum and a lower chance of postglacial gene flow from other populations into this rather marginal part in the southern Balkans. Horizontal starch gel electrophoresis of proteins from 35~loci revealed three alleles (Es-1 ^–162, Pep-2 ¹¹⁴, Mpi ⁸⁸) at low frequencies, which were not found in Bulgarian or any other brown hare population. In contrast, some alleles from the populations from Bulgaria and other regions of Europe were absent in the Greek samples. Population genetic statistics indicated only a slight tendency of increased gene pool diversity in Greek hares, little substructuring in Greek and Bulgarian populations, respectively, as well as an only slightly lower level of gene flow between the two neighboring regions, as compared to the gene flow within each region. The results conform to the hypothesis of a Late Pleistocene refugial population in the southern Balkans, with some few specific nuclear gene pool characteristics, but little effect on the overall genetic differentiation between Greek and Bulgarian hares.     

Genetic diversity and population connectivity of the Asian green mussel Perna viridis in South China Sea,inferred from mitochondria DNA markers

The Asian green mussel Perna viridis is ecologically and economically important in the coastal regions of China. In order to characterize the genetic diversity and population connectivity of P. viridis in South China Sea, a 664 bp region of mitochondrial COI gene and a 293 bp region of 16S rRNA gene were sequenced and analyzed for 78 and 92 individuals from four populations in South China Sea, respectively. A total of 15 haplotypes were defined by 14 variable nucleotide sites in COI gene, and 7 haplotypes by 6 variable nucleotide sites in 16S rRNA gene. High haplotype diversity and low nucleotide diversity were observed in COI gene, while moderate haplotype diversity and low nucleotide diversity were observed in 16S rRNA gene. Pairwise F_ST values of COI gene were all negative and those of 16S rRNA gene ranged from −0.01409 to 0.10289. The results showed that no significant genetic divergence (or shallow genetic structure) and high levels of population connectivity among the four populations of P. viridis in South China Sea.     

Genetic characterization of brown meagre (Sciaena umbra) and Shi Drum (Umbrina cirrosa) populations

This study analyzed population structure of fish, brown meagre (Sciaena umbra) and Shi drum (Umbrina cirrosa) from the Black Sea, the Aegean Sea, and the Mediterranean Sea. The Cytb and 16S rRNA genes of Shi drum and brown meagre fish were sequenced. In brown meagre and Shi drum, 25 and 20 haplotypes, respectively, of Cytb gene were identified; while for 16S rRNA, 4 and 8 haplotypes were identified. Nucleotide diversities of 16S rRNA and Cytb gene sites were found to be 83.8% and 52.6%, respectively, for brown meagre; while that for Shi drum were of 80.5% and 73.6%, respectively. There was a significant relationship between geographic distance and the genetic distance of the fish. Since Shi drum and brown meagre are migratory species, they can migrate between the seas. The lack of barriers among different populations facilitates the gene flow among the populations belonging to different regions. Since there is no information available on Shi drum and brown meagre population genetics, this study may be useful to understand the genetic diversity of these species to assist fishery managers for the management of these resources in terms of conservation and sustainability.     

Genetic structure of the endangered species <Emphasis Type="Italic">Pinna nobilis</Emphasis> (Mollusca: Bivalvia) inferred from mtDNA sequences

This study examines the population genetic structure of the endangered bivalve Pinna nobilis (Mollusca: Bivalvia), based on novel mtDNA sequences (partial COI and 16S rDNA mtDNA genes). The analyzed nucleotide sequences of COI were 729 bp in size, coding for a 243 amino acid peptide, while the analyzed nucleotide sequences of 16S rDNA were 489 bp in size. These sequences of P. nobilis were the first DNA sequences of the species submitted to any Genetic Data Base. Population samples from four geographic regions from Greece, as well as a population sample of Atrina fragilis (as an outgroup) were used. High values of haplotypic diversity were found in the population samples of P. nobilis, based on the COI sequences. A single base in the analyzed 16S rDNA sequences was different in all analyzed individuals from a single population sample (Chios island) differentiating it from the other ones. These mtDNA sequences could be informative for further genetic analyses of the endangered species, contributing in conservation plans for its protection and/or aquaculture investigations.     

Mitochondrial DNA phylogeography of Acrossocheilus paradoxus (Cyprinidae) in Taiwan

Nucleotide sequences of 3' end of the cytochrome b gene, tRNA genes, D-loop control region, and the 5' end of the 12S rRNA of mitochondrial DNA (mtDNA) were used to assess the genetic and phylogeographic structure of Acrossocheilus paradoxus populations, a Cyprinidae fish of Taiwan. A hierarchical examination of populations in 12 major streams from three geographical regions using an analysis of molecular variance (AMOVA) indicates high genetic differentiation both among populations (PhiST = 0.511, P < 0.001) and among regions (PhiCT = 0.368, P < 0.001). Limited migration largely contributed to the genetic differentiation. High nucleotide diversity (1.13%) and haplotype diversity (0.80%) were detected among populations. The degree of genetic differentiation was correlated with geographical distance between populations, a result consistent with the one-dimensional stepping stone models. A neighbour-joining tree recovered by (DAMBE) supports the pattern of isolation by distance and reveals a closer relationship between populations of the central and southern regions. A minimum spanning network based on nucleotide substitutions reflected migration routes from populations of the central region to the northern and southern regions, respectively. Postglacial colonization and expansion can explain the phylogeographical pattern. Single and ancient migration events may have allowed the northern region to attain the monophyly of mtDNA alleles. In contrast, most populations within geographical regions are either paraphyletic or polyphyletic due to the relatively shorter time period for coalescence. Both low haplotype number and genetic variability suggest a bottleneck event in the Chingmei population of northern Taiwan. Based on coalescence theory, the monophyly of the Tungkang population of the southern region may be associated with a founder event.     

Molecular phylogeny, taxonomy and evolution of the land snail genus Iberus (Pulmonata: Helicidae)

Partial DNA sequences of two mitochondrial genes [cytochrome oxidase subunit I (COI) and 16S rRNA] from 59 specimens of Iberus were used to test the validity of the described morphospecies of this genus, and examine genetic divergences within and between main phylogenetic groups. Both gene fragments showed phylogenetic concordance. The COI gene was found to be faster evolving than the 16S gene and was fully protein-coding with no insertions or deletions. 16S rRNA was more informative than COI for resolving basal nodes. Both individual and combined analyses of the two gene fragments revealed five main phylogroups. These five groups are genetically unique lineages that are allopatrically distributed and considered to have full species status. Further subdivisions were also considered. Shell morphology was suitable for delimiting species boundaries, but several incongruences between morphology and mtDNA phylogeny were observed. These incongruences were considered consequence of hybridization between Iberus cobosi and Iberus marmoratus , and the result of shell shape polymorphism in Iberus rositai . According to spatial patterns of sequence divergence, life habits and shell morphology may be concluded that the keeled-flat shelled snails independently originated several times within Iberus and they could represent cases of similar shell adaptation to a karstic arid environment.     

Reduced genetic variation and strong genetic population structure in the freshwater killifish Valencia letourneuxi (Valenciidae) based on nuclear and mitochondrial markers

The genetic variation of the critically endangered Corfu killifish (Valencia letourneuxi), an endemic freshwater fish species of the western Balkans, was assessed for nine populations sampled in eight water systems in western continental Greece, the Peloponnese and the Ionian Island of Corfu, using mitochondrial and microsatellite markers. The analyses were based on data from three mtDNA regions (D‐loop, COI and 16S rRNA sequences) and 14 microsatellite loci. Samples from the congeneric species Valencia hispanica and the phylogenetically closely related species Aphanius fasciatus were also used in the study as outgroups. Both the mitochondrial and the microsatellite analyses revealed three distinct population groupings associated with the geographical distribution of the populations: one southern group occupying rivers draining to the Patraikos Gulf, the second one including the populations flowing into the Amvrakikos Gulf and the third, more northern group, including the other populations from rivers in Corfu Island and Epirus flowing into the Ionian Sea. Within these groupings there is limited genetic differentiation between populations; in addition, there is reduced intrapopulation genetic variation, evidenced by low heterozygosity values, number of alleles and haplotype diversity. In terms of taxonomic implications and appropriate management actions for conservation, our data suggest that the major population groups should be regarded at least as three distinct conservation units (CUs), with translocation and restocking actions to take place only within the geographical range of the CU concerned. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 334–349.     

Intraspecific mitochondrial DNA variation in the Colorado potato beetle,Leptinotarsa decemlineata (Coleoptera: Chrysomelidae)

Restriction endonuclease analyses of mitochondrial DNA (mtDNA) were used to examine genetic variability and population structure inLeptinotarsa decemlineata (Say). A group of three enzymes, EcoRI, HpaI, and PstI, was used to reveal polymorphism both within and among some of the 10 populations tested, yielding 16 haplotypes in combination. The frequencies of these 16 haplotypes differed significantly across geographic regions, indicating some partitioning of mtDNA haplotypes. Estimates of mtDNA sequence divergence (δ) between haplotypes ranged from 0.016 to 0.135%, suggesting local differentiation of mtDNA in some populations. Analysis of these data suggests that Texas was colonized by more than one mtDNA lineage, most likely originating in Mexico. We hypothesize that a larger founder size for the initial introductions or high levels of variability in the parent population at the edge of the CPB expanding range led to the initial partitioning of haplotypes observed in samples from Texas.     

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.     

A molecular phylogeny for the South African limbless lizard taxa of the subfamily Acontinae (Sauria: Scincidae) with special emphasis on relationships within Acontias

In the present study, relationships among three genera Acontias, Acontophiops, and Typhlosaurus, that comprise the South African limbless lizard subfamily Acontinae, were assessed with partial sequences of the 16S rRNA mitochondrial DNA gene. In addition, relationships within Acontias were further investigated using sequence data from the cytochrome oxidase I gene (COI). Maximum likelihood and maximum parsimony analyses of the 16S rRNA mtDNA data revealed that within this subfamily, Typhlosaurus is basal while Acontophiops and Acontias are sister taxa. Based on the 16S rRNA mtDNA data, the relationships within Acontias placed A. meleagris orientalis as the sister taxon of A. percivali tasmani, with A. m. orientalis lineacauda morph and A. m. meleagrus being the sister taxa to this group. The small-bodied skinks A. lineatus lineatus and A. l. tristis formed a monophyletic group, with the medium-bodied species A. gracilicauda gracilicauda being their sister taxon. Analyses of the COI gene for Acontias place A. m. orientalis as the sister taxon of A. p. tasmani with both A. meleagris meleagris and A. m. orientalis lineacauda being distinct. In contrast to the 16S rRNA mtDNA data, the COI data placed A. g. gracilicauda as the sister taxon to these medium-bodied species; while the subspecies status of the small-bodied taxa A. l. lineatus and A. l. tristis is reaffirmed. Combined analysis of both gene fragments for Acontias taxa recovered the same clades as found using only COI data. Systematic affinities in Acontias are discussed. These results indicate that Acontias is more species rich than previously thought.     

Mitochondrial DNA variation and population genetic structure of the northern redbelly dace (Phoxinus eos)

Two sections of the control region and the genes coding for NADH dehydrogenase sub-units 5 and 6 (ND-5/6) of mitochondrial DNA (mtDNA) were amplified from Phoxinus eos with the polymerase chain reaction. Both sections of the control region were sequenced directly while the ND-5/6 fragment was sequenced in from each end only. Additionally, the entire ND-5/6 fragment was examined for sequence variation using RFLP analysis. No sequence variation was detected in the control region among 70 individuals sampled from 18 populations across three Ontario regions (Spanish River, Madawaska R. and Cataraqui R.). To examine ND-5/6 variation, a total of 75 individuals were sampled from five populations representing two of the three regions (Madawaska River and Cataraqui R.). Six haplotypes were detected by direct sequencing and four by RFLP analysis. Estimates of population subdivision from RFLP data, sequence analysis, and the two data sets combined for the ND-5/6 fragment, suggest that gene flow is restricted within and between regions. However, estimates of sequence divergence for both sequence and RFLP analysis of this fragment suggested that populations were either founded by already differentiated populations or that populations were founded by a single stock and divergence between regions occurred prior to isolation of populations within regions. These estimates of population structure are much greater than those obtained from allozyme analysis. Additionally, high levels of heterozygosity in nuclear DNA, but low mtDNA diversity suggests that populations have experienced reductions in population size sufficient to reduce only mtDNA variation. Random lineage extinction and limited time for the accumulation of new mutations are likely responsible for low levels of mtDNA variation in ND-5/6 and the control region, while functional constraints may limit variation more than expected in the control region in dace and other fishes.     

COI is better than 16S rRNA for DNA barcoding Asiatic salamanders (Amphibia: Caudata: Hynobiidae)

The 5' region of the mitochondrial DNA (mtDNA) gene cytochrome c oxidase I (COI) is the standard marker for DNA barcoding. However, because COI tends to be highly variable in amphibians, sequencing is often challenging. Consequently, another mtDNA gene, 16S rRNA gene, is often advocated for amphibian barcoding. Herein, we directly compare the usefulness of COI and 16S in discriminating species of hynobiid salamanders using 130 individuals. Species identification and classification of these animals, which are endemic to Asia, are often based on morphology only. Analysis of Kimura 2-parameter genetic distances (K2P) documents the mean intraspecific variation for COI and 16S rRNA genes to be 1.4% and 0.3%, respectively. Whereas COI can always identify species, sometimes 16S cannot. Intra- and interspecific genetic divergences occasionally overlap in both markers, thus reducing the value of a barcoding gap to identify genera. Regardless, COI is the better DNA barcoding marker for hynobiids. In addition to the comparison of two potential markers, high levels of intraspecific divergence in COI (>5%) suggest that both Onychodactylus fischeri and Salamandrella keyserlingii might be composites of cryptic species.     

Mitochondrial DNA of the blowflyPhormia regina: restriction analysis and gene localization

A study of an invertebrate mitochondrial genome, that of the blowflyPhormia regina, has been initiated to compare its structural and functional relatedness to other metazoan mitochondrial genomes. A restriction map of mitochondrial DNA (mtDNA) isolated from sucrose gradient-purified mitochondria has been established using a combination of single and double restriction endonuclease digestions and hybridizations with isolated mtDNA fragments, revealing a genome size of 17.5 kilobases (kb). A number of mitochondrial genes including those encoding the 12 S and 16 S ribosomal RNA, the cytochromec oxidase I subunit (COI) and an unidentified open reading frame (URF2) have been located on thePhormia mtDNA by Southern blot analysis using as probes both isolated mtDNA fragments and oligonucleotides derived from the sequences of previously characterized genes from rat andDrosophila yakuba mtDNAs. These data indicate that for those regions examined, the mitochondrial genome organization of blowfly mtDNA is the same as that ofDrosophila yakuba, the order being COI-URF2-12 S-16 S. These data also report the presence of an A + T-rich region, located as a 2.5-kb region between the URF2 and the 12 S rRNA genes, and its amplification by the polymerase chain reaction is described.     

Geographic differentiation of <Emphasis Type="Italic">Cobitis shikokuensis</Emphasis> inferred from mtDNA RFLP analysis

Mitochondrial DNA divergence among populations of the Japanese spinous loach Cobitis shikokuensis, endemic to Shikoku Island, was investigated by restricted fragment length polymorphism analysis. A total of 68 restriction sites on DNA fragments from the cytochrome b to D-loop regions and from the 12S rRNA to 16S rRNA regions, amplified by PCR, were analyzed. A total of 12 haplotypes (plus 6 in outgroups) were detected in 268 specimens collected from 19 localities in seven rivers (and 41 specimens from four localities in three rivers in outgroups). Three of the seven river populations of C. shikokuensis were shown to have unique haplotypes, and four of the seven river populations were monomorphic. The nested structure of the haplotype network for populations of C. shikokuensis exhibited two large clades corresponding to (1) populations from the Shimanto River and its neighbors and (2) two genetically divergent populations in the Shigenobu and Iwamatsu Rivers. The population from the Shimanto River, the largest river inhabited by C. shikokuensis, maintains great haplotype diversity as well as the allozyme diversity previously reported. On the other hand, populations from the Hiji River, the second largest river, which exhibited the highest allozyme diversity, were monomorphic in their mtDNA. The nested clade analysis (NCA) revealed that past fragmentation between the above two clades could occur in the initial distribution process of C. shikokuensis. The large genetic divergence of two river populations from the Shigenobu and Iwamatsu Rivers was inferred to be caused by a process of long distance colonization and fragmentation. MtDNA introgression into the Hiji River population from southern river populations was suggested. Taking genetic divergence into consideration, each river population of C. shikokuensis should be conserved separately as like a distinct species, and conservation programs for the small populations showing less genetic variability should be invoked as soon as possible.     

Burmoniscus kitadaitoensis Nunomura, 2009 (Crustacea,Isopoda, Oniscidea) from southern Japan,a junior synonym of B. meeusei (Holthuis, 1947)

Re-examination of the holotype of Burmoniscus kitadaitoensis Nunomura, 2009 from Kitadaitojima Island, southern Japan reveals that this species is a junior synonym of B. meeusei (Holthuis, 1947). Partial regions of mitochondrial COI, 12S and 16S rRNA genes, and nuclear 18S and 28S rRNA genes were detected for species identification in the future.     

Reconstructing the colonisation and diversification history of the endemic freshwater crab (Seychellum alluaudi) in the granitic and volcanic Seychelles Archipelago

The endemic, monotypic freshwater crab species Seychellum alluaudi was used as a template to examine the initial colonisation and evolutionary history among the major islands in the Seychelles Archipelago. Five of the “inner” islands in the Seychelles Archipelago including Mahé, Praslin, Silhouette, La Digue and Frégate were sampled. Two partial mtDNA fragments, 16S rRNA and cytochrome oxidase subunit I (COI) was sequenced for 83 specimens of S. alluaudi. Evolutionary relationships between populations were inferred from the combined mtDNA dataset using maximum parsimony, maximum likelihood and Bayesian inferences. Analyses of molecular variance (AMOVA) were used to examine genetic variation among and within clades. A haplotype network was constructed using TCS while BEAST was employed to date the colonisation and divergence of lineages on the islands. Phylogenetic analyses of the combined mtDNA data set of 1103 base pairs retrieved a monophyletic S. alluaudi group comprised three statistically well-supported monophyletic clades. Clade one was exclusive to Silhouette; clade two included samples from Praslin sister to La Digue, while clade three comprised samples from Mahé sister to Frégate. The haplotype network corresponded to the three clades. Within Mahé, substantial phylogeographic substructure was evident. AMOVA results revealed limited genetic variation within localities with most variation occurring among localities. Divergence time estimations predated the Holocene sea level regressions and indicated a Pliocene/Pleistocene divergence between the three clades evident within S. alluaudi. The monophyly of each clade suggests that transoceanic dispersal is rare. The absence of shared haplotypes between the three clades, coupled with marked sequence divergence values suggests the presence of three allospecies within S. alluaudi.     

Genetic variation in the major histocompatibility complex of the European brown hare (Lepus europaeus) across distinct phylogeographic areas

The major histocompatibility complex is one of the best studied systems in vertebrates providing evidence for the long-term action of selection. Here, we examined the intra- and inter-population genetic diversity of the MHC class II DRB locus in European brown hare (Lepus europaeus) and correlated the results with genetic variability already estimated from the MHC DQA locus and from maternally (mitochondrial DNA (mtDNA)) and biparentally (allozymes, microsatellites) inherited loci. L. europaeus showed remarkable genetic polymorphism in both DQA and DRB1 loci. The Anatolian populations exhibited the highest genetic polymorphism for both loci. Balancing selection has established increased variability in the European populations despite the founder effects after the last glaciation. Different evolutionary rates were traced for DRB1 and DQA loci, as evidenced by the higher number of common DRB1 than DQA alleles and the greater differences between DRB1 alleles with common origin in comparison with DQA alleles. The high number of rare alleles with low frequencies detected implies that frequency-dependent selection drives MHC evolution in the brown hare through the advantage of rare alleles. Both loci were under the influence of positive selection within the peptide-binding region. The functional polymorphism, recorded as amino acid substitutions within the binding pockets, fell also within distinct geographic patterns, yet it was much narrower than the genetic polymorphism. We hypothesize that certain structural and functional characteristics of the binding pockets set limitations to the actual shape of genetic polymorphism in MHC.