Understanding the phylogeographic patterns of European hedgehogs, Erinaceus concolor and E. europaeus using the MHC

Patterns of genetic differentiation imply distinct phylogeographic history of the mosquito species Anopheles messeae and Anopheles daciae in Eurasia

Detailed knowledge of phylogeography is important for control of mosquito species involved in the transmission of human infectious diseases. Anopheles messeae is a geographically widespread and genetically diverse dominant vector of malaria in Eurasia. A closely related species, An. daciae, was originally distinguished from An. messeae based on five nucleotide substitutions in its ribosomal DNA (rDNA). However, the patterns of phylogeographic history of these species in Eurasia remain poorly understood. Here, using internal transcribed spacer 2 (ITS2) of rDNA and karyotyping for the species identification we determined the composition of five Anopheles species in 28 locations in Eurasia. Based on the frequencies of 11 polymorphic chromosomal inversions used as genetic markers, a large-scale population genetics analysis was performed of 1932 mosquitoes identified as An. messeae, An. daciae and their hybrids. The largest genetic differences between the species were detected in the X sex chromosome suggesting a potential involvement of this chromosome in speciation. The frequencies of autosomal inversions in the same locations differed by 13%-45% between the species demonstrating a restricted gene flow between the species. Overall, An. messeae was identified as a diverse species with a more complex population structure than An. daciae. The clinal gradients in frequencies of chromosomal inversions were determined in both species implicating their possible involvement in climate adaptations. The frequencies of hybrids were low ~1% in northern Europe but high up to 50% in south-eastern populations. Thus, our study revealed critical differences in patterns of phylogeographic history between An. messeae and An. daciae in Eurasia. This knowledge will help to predict the potential of the malaria transmission in the northern territories of the continent.

Both selection and drift drive the spatial pattern of adaptive genetic variation in a wild mammal

The major histocompatibility complex (MHC) has been intensively studied for the relative effects of different evolutionary forces in recent decades. Pathogen-mediated balancing selection is generally thought to explain the high polymorphism observed in MHC genes, but it is still unclear to what extent MHC diversity is shaped by selection relative to neutral drift. In this study, we genotyped MHC class II DRB genes and 15 neutral microsatellite loci across 26 geographic populations of European badgers (Meles meles) covering most of their geographic range. By comparing variation of microsatellite and diversity of MHC at different levels, we demonstrate that both balancing selection and drift have shaped the evolution of MHC genes. When only MHC allelic identity was investigated, the spatial pattern of MHC variation was similar to that of microsatellites. By contrast, when functional aspects of the MHC diversity (e.g., immunological supertypes) were considered, balancing selection appears to decrease genetic structuring across populations. Our comprehensive sampling and analytical approach enable us to conclude that the likely mechanisms of selection are heterozygote advantage and/or rare-allele advantage. This study is a clear demonstration of how both balancing selection and genetic drift simultaneously affect the evolution of MHC genes in a widely distributed wild mammal.

Genetic diversity and structure of the hedgehogs Erinaceus europaeus and Erinaceus roumanicus: evidence for ongoing hybridization in Eastern Europe

Secondary contact zones between related species are key to understanding speciation mechanisms. The Central European sympatry zone of West European (Erinaceus europaeus) and northern white-breasted (Erinaceus roumanicus) hedgehogs is well studied, whereas data on the Eastern European sympatry zone are scarce. We examined the genetic variation in Russian populations using the mitochondrial Cytb gene, TTR intron 1 and 11 microsatellites to assess genetic variability and distribution patterns. In contrast to the Central European sympatry zone, we found evidence of ongoing hybridization between the two species in the sympatry zone of European Russia, where the proportion of individuals with mixed ancestry was c. 20%. Our data indicate bi-directional mtDNA introgression, but with a higher frequency of E. europaeus haplotypes in hybrids. The proportion of pure specimens with introgressed mitotypes is higher in E. roumanicus than in E. europaeus. Nuclear data showed the prevalence of the genetic contribution from E. roumanicus in admixed individuals. Demographic analyses indicated recent population growth in E. europaeus and little change in E. roumanicus, suggesting that E. europaeus colonized East Europe later than E. roumanicus.

Evaluating the Impact of Wildlife Shelter Management on the Genetic Diversity of Erinaceus europaeus and E. roumanicus in Their Contact Zone

Simple Summary

Hedgehogs are regularly brought to wildlife shelters. Depending on the area from where animals are accepted, translocation can occur between different regions or populations. In this study, the genetic diversity of wild hedgehog populations was compared with “shelter populations” within central Europe focusing on the western contact zone between both European hedgehog species. Some shelters were hosting both species at the same time, in one this could be shown genetically. Generally, no difference in genetic diversity between shelter individuals and wild populations was found. Two shelters from Innsbruck hosted individuals that probably belong to two subpopulations. This indicates that shelter management-related translocations could facilitate gene flow across a dispersal barrier.

Abstract

Hedgehogs are among the most abundant species to be found within wildlife shelters and after successful rehabilitation they are frequently translocated. The effects and potential impact of these translocations on gene flow within wild populations are largely unknown. In this study, different wild hedgehog populations were compared with artificially created “shelter populations”, with regard to their genetic diversity, in order to establish basic data for future inferences on the genetic impact of hedgehog translocations. Observed populations are located within central Europe, including the species Erinaceus europaeus and E. roumanicus. Shelters were mainly hosting one species; in one case, both species were present syntopically. Apart from one exception, the results did not show a higher genetic diversity within shelter populations, indicating that individuals did not originate from a wider geographical area than individuals grouped into one of the wild populations. Two shelters from Innsbruck hosted individuals that belonged to two potential clusters, as indicated in a distance analysis. When such a structure stems from the effects of landscape elements like large rivers, the shelter management-related translocations might lead to homogenization across the dispersal barrier.

Genetic structure of the European hedgehog (Erinaceus europaeus) in Denmark

Objectives

Low genetic diversity can lead to reduced average fitness in a population or even extinction. Preserving genetic connectivity across fragmented landscapes is therefore vital to counteract the negative consequences of genetic drift and inbreeding. This study aimed to assess the genetic composition and consequently the conservation status of a nationwide sample of European hedgehogs (Erinaceus europaeus) in Denmark.

Methods

We applied an adaptation of the genotyping by sequencing (GBS) technique to 178 individuals from six geographically distinct populations. We used a Bayesian clustering method to subdivide individuals into genetically distinct populations. We estimated individual observed (iH_O), observed (H_O), and unbiased expected (uH_E) heterozygosity, inbreeding coefficient (F_IS), percentage of polymorphic loci (P%) and tested for deviations from Hardy-Weinberg equilibrium (HWE). We used linear models to test for potential anthropogenic effects on the genetic variability of hedgehogs with iH_O, uH_E, P% and F_IS as response variables, and assessed the demographic history of the population.

Results

The Danish hedgehog population is composed of three genetic clusters. We found a mean P% of 54.44–94.71, a mean uH_E of 0.126–0.318 and a mean H_O of 0.124–0.293 in the six populations. The F_IS was found to be significantly positive for three of the six populations. We detected a large heterogeneity of iH_O values within populations, which can be due to inbreeding and/or fragmentation. F_IS values decreased with increasing farmland density, but there was no significant association with human population or road density.

Conclusions

We found a low level of genetic variability and evidence for genetic substructure and low effective population size, which are all consequences of habitat fragmentation. We failed to detect signs of a recent population bottleneck or population increase or decline. However, because the test only identifies recent changes in population size, we cannot reject the possibility of a longer-term decline in the Danish hedgehog population.

Comparative phylogeography of Aedes mosquitoes and the role of past climatic change for evolution within Africa

The study of demographic processes involved in species diversification and evolution ultimately provides explanations for the complex distribution of biodiversity on earth, indicates regions important for the maintenance and generation of biodiversity, and identifies biological units important for conservation or medical consequence. African and forest biota have both received relatively little attention with regard to understanding their diversification, although one possible mechanism is that this has been driven by historical climate change. To investigate this, we implemented a standard population genetics approach along with Approximate Bayesian Computation, using sequence data from two exon-primed intron-crossing (EPIC) nuclear loci and mitochondrial cytochrome oxidase subunit I, to investigate the evolutionary history of five medically important and inherently forest dependent mosquito species of the genus Aedes. By testing different demographic hypotheses, we show that Aedes bromeliae and Aedes lilii fit the same model of lineage diversification, admixture, expansion, and recent population structure previously inferred for Aedes aegypti. In addition, analyses of population structure show that Aedes africanus has undergone lineage diversification and expansion while Aedes hansfordi has been impacted by population expansion within Uganda. This congruence in evolutionary history is likely to relate to historical climate-driven habitat change within Africa during the late Pleistocene and Holocene epoch. We find differences in the population structure of mosquitoes from Tanzania and Uganda compared to Benin and Uganda which could relate to differences in the historical connectivity of forests across the continent. Our findings emphasize the importance of recent climate change in the evolution of African forest biota.

Selection on MHC in a Context of Historical Demographic Change in 2 Closely Distributed Species of Tuco-tucos (Ctenomys australis and C. talarum)

Selection necessarily acts within the same current and historical demographic framework as neutral evolutionary processes, and the outcome of the interplay between these forces may vary according to their relative strength. In this study, we compare the variation at a major histocompatibility complex (MHC) locus (DRB exon 2), typically subject to strong diversifying selection, and mitochondrial diversity (D-loop) across populations encompassing the entire distribution of 2 species of South American subterranean rodents: Ctenomys australis and C. talarum (tuco-tucos). Although these species are parapatric along most of their distribution, historically they have followed distinct demographic trajectories associated with sea level changes during the Quaternary. We surveyed 8 populations of C. australis and 15 of C. talarum, from which we analyzed 70 and 212 D-loop haplotypes and 91 and 346 DRB genotypes, respectively. Both species have gone through a recent demographic expansion; however, the signal of this process only encompasses the entire distribution of one of the species: C. australis. While balancing selection on MHC in C. talarum-enhanced DRB diversity at the local level compared to D-loop, although not promoting divergence among populations, in C. australis local diversifying selection may have driven higher population differentiation at DRB than at D-loop. Our findings reinforce the idea that the relative strength of selection acting on MHC genes varies spatially and temporally within and among species, even between species using the same macrohabitat and exposed to similar immune challenges.

Microevolution of sympatry: landscape genetics of hedgehogs Erinaceus europaeus and E. roumanicus in Central Europe

We used the mitochondrial control region and nuclear microsatellites to assess the distribution patterns, population structure, demography and landscape genetics for the hedgehogs Erinaceus europaeus and Erinaceus roumanicus in a transect of the mid-European zone of sympatry. E. roumanicus was less frequent and restricted to regions with lower altitudes. Demographic analyses suggested recent population growth in this species. A comparison of patterns in the spatial variability of mitochondrial and nuclear DNA indicated less sex-biased dispersal and higher levels of gene flow in E. roumanicus. No evidence of recent hybridisation or introgression was detected. We interpreted these results by comparing with phylogeographic and palaeontological studies as well as with the occurrence of hybridisation in the Russian contact zone. We propose that Central Europe was colonised by E. roumanicus by the beginning of the Neolithic period and that there was a subsequent reinforcement stage as well as the formation of a zone of sympatry after the complete reproductive isolation of both species.

Multiple Asian pig origins revealed through genomic analyses

Previous mitochondrial DNA (mtDNA) studies have suggested that European and Asian pig populations were derived through multiple domestication events. We investigated whether domestic pig populations were derived from distinct ancestors within their respective regions, using eight domestic breeds (five European and three Asian), and also European and Asian wild boar populations. Genomic analyses utilized 21 microsatellite markers (MS) selected for their distribution across the pig genome in addition to the mtDNA D-loop region. The number of alleles per MS loci ranged from 8 (Sw2008) to 16 (S0097 and S0218). Few significant departures from Hardy-Weinberg equilibrium were detected, suggesting the absence of heterozygote deficiencies. Analyses within populations revealed observed mean heterozygosity from 0.48 (Erhualian) to 0.68 (Dutch WB) and an expected mean heterozygosity from 0.53 (Hampshire) to 0.80 (Japanese WB) with effective alleles ranging from 2.28 (Hampshire) to 3.74 (French WB). Wild boar populations demonstrated a higher level of heterozygosity than domestic breeds. Genetic differentiation estimated by fixation indices (F(ST)) ranged from 0.021 (Yorkshire and Duroc) to 0.410 (Meishan and Hampshire) and was consistent with previous mtDNA analysis. Both phylogenetic and principal component analyses revealed a distinct separation of European and Asian derived populations with tight clustering of the European domestic breeds. Conversely, the use of both MS and mtDNA clarified that the Asian populations were comprised of three groups, one represented by Erhualian and Meishan breed, the second represented by Lanyu pigs and the third represented by the Asian wild boars. The current findings support the hypothesis that Asian domestic populations were derived from multiple Asian ancestral origins whereas the European domestic populations represent a single ancestral European lineage.

Diversity and evolutionary history of the MHC DQA gene in leporids

The European rabbit (Oryctolagus cuniculus) is used as a model for many human diseases, yet comparatively little is known of its genetics, particularly at important loci such as the major histocompatibility complex (MHC). This study investigated genetic diversity and evolutionary history of the DQA gene in a range of leporid species by analysing coding sequence diversity of exon 2 and intron 2 in 53 individuals of 16 different species. Fifty leporid DQA alleles were detected, including 13 novel European rabbit alleles. In the rabbit, the highest levels of diversity were observed in wild rabbits from Portugal, with wild rabbits from England and domestic rabbits showing less diversity. Within the sample, several recombination events were detected and trans-specific evolution of alleles was evidenced, both being general characteristics of mammalian MHC genes. Positive selection is implicated as operating on six codons within exon 2, which are also subject to positive selection in other mammals. Some of these positions are putative antigen recognition sites and underline the importance of pathogen-driven selection on these MHC genes.

Did glacial advances during the Pleistocene influence differently the demographic histories of benthic and pelagic Antarctic shelf fishes?--Inferences from intraspecific mitochondrial and nuclear DNA sequence diversity

BACKGROUND

Circum-Antarctic waters harbour a rare example of a marine species flock - the Notothenioid fish, most species of which are restricted to the continental shelf. It remains an open question as to how they survived Pleistocene climatic fluctuations characterised by repeated advances of continental glaciers as far as the shelf break that probably resulted in a loss of habitat for benthic organisms. Pelagic ecosystems, on the other hand, might have flourished during glacial maxima due to the northward expansion of Antarctic polar waters. In order to better understand the role of ecological traits in Quaternary climatic fluctuations, we performed demographic analyses of populations of four fish species from the tribe Trematominae, including both fully benthic and pelagic species using the mitochondrial cytochrome b gene and an intron from the nuclear S7 gene.

RESULTS

Nuclear and cytoplasmic markers showed differences in the rate and time of population expansions as well as the likely population structure. Neutrality tests suggest that such discordance comes from different coalescence dynamics of each marker, rather than from selective pressure. Demographic analyses based on intraspecific DNA diversity suggest a recent population expansion in both benthic species, dated by the cyt b locus to the last glacial cycle, whereas the population structure of pelagic feeders either did not deviate from a constant-size model or indicated that the onset of the major population expansion of these species by far predated those of the benthic species. Similar patterns were apparent even when comparing previously published data on other Southern Ocean organisms, but we observed considerable heterogeneity within both groups with regard to the onset of major demographic events and rates.

CONCLUSION

Our data suggest benthic and pelagic species reacted differently to the Pleistocene ice-sheet expansions that probably significantly reduced the suitable habitat for benthic species. However, the asynchronous timing of major demographic events observed in different species within both "ecological guilds", imply that the species examined here may have different population and evolutionary histories, and that more species should be analysed in order to more precisely assess the role of life history in the response of organisms to climatic changes.