SNP Identification through Transcriptome Analysis of the European Brown Hare (Lepus europaeus): Cellular Energetics and Mother's Curse

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.

The association between mitochondrial genetic variation and reduced colony fitness in an invasive wasp

Despite the mitochondrion's long-recognized role in energy production, mitochondrial DNA (mtDNA) variation commonly found in natural populations was assumed to be effectively neutral. However, variation in mtDNA has now been increasingly linked to phenotypic variation in life history traits and fitness. We examined whether the relative fitness in native and invasive common wasp (Vespula vulgaris) populations in Belgium and New Zealand (NZ), respectively, can be linked to mtDNA variation. Social wasp colonies in NZ were smaller with comparatively fewer queen cells, indicating a reduced relative fitness in the invaded range. Interestingly, queen cells in this population were significantly larger leading to larger queen offspring. By sequencing 1,872 bp of the mitochondrial genome, we determined mitochondrial haplotypes and detected reduced genetic diversity in NZ. Three common haplotypes in NZ frequently produced many queens, whereas the four rare haplotypes produced significantly fewer or no queens. The entire mitochondrial genome for each of these haplotypes was sequenced to identify polymorphisms associated with fitness reduction. We found 16 variable sites; however, no nonsynonymous mutation that was clearly causing impaired mitochondrial function was detected. We discuss how detected variants may alter secondary structures, gene expression or mito-nuclear interactions, or could be associated with nuclear-encoded variation. Whatever the ultimate mechanism, we show reduced fitness and mtDNA variation in an invasive wasp population as well as specific mtDNA variants associated with fitness variation within this population. Ours is one of only a few studies that confirm fitness impacts of mtDNA variation in wild nonmodel populations.

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

Background

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

Results

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

Conclusion

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

Electronic supplementary material

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

Mountain hare transcriptome and diagnostic markers as resources to monitor hybridization with European hares

We report the first mountain hare (Lepus timidus) transcriptome, produced by de novo assembly of RNA-sequencing reads. Data were obtained from eight specimens sampled in two localities, Alps and Ireland. The mountain hare tends to be replaced by the invading European hare (Lepus europaeus) in their numerous contact zones where the species hybridize, which affects their gene pool to a yet unquantified degree. We characterize and annotate the mountain hare transcriptome, detect polymorphism in the two analysed populations and use previously published data on the European hare (three specimens, representing the European lineage of the species) to identify 4 672 putative diagnostic sites between the species. A subset of 85 random independent SNPs was successfully validated using PCR and Sanger sequencing. These valuable genomic resources can be used to design tools to assess population status and monitor hybridization between species.

Mitogenomic analysis in European brown hare (Lepus europaeus) proposes genetic and functional differentiation between the distinct lineages

European brown hare is a small game species spreading across Europe to Asia Minor, with important economic traits. Population genetics studies using mitochondrial DNA markers have revealed the existence of two major phylogeographic lineages, the European and the Anatolian. European lineage is further divided in the European type halpogroup and south-eastern European type halpogroup, while Anatolian consists only by the Anatolian/Middle Eastern type halpogroup. All three haplogroups show a discrete geographical distribution, with an overlapping zone forming in North-East Greece and Bulgaria, forming a contact zone. Despite the existence of a contact zone, European haplotype was never detected in Anatolia and vice versa, proposing the presence of genetic barriers responsible for this phenomenon. In this study, we analyzed the whole mitochondrial genomes of specimens originating from both lineages, aiming to detect the genetic and functional differentiation of the oxidative phosphorylation complexes that are encoded by mtDNA that could lead gradually to the reproductive isolation of the lineages.