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

Species-specific variation in mitochondrial genome tandem repeat polymorphisms in hares (Lepus spp., Lagomorpha, Leporidae) provides insight into their evolution

The non-coding regions of the mitochondrial DNAs (mtDNAs) of hares, rabbits, and pikas (Lagomorpha) contain short (∼20 bp) and long (130-160 bp) tandem repeats, absent in related mammalian orders. In the presented study, we provide in-depth analysis for mountain hare (Lepus timidus) and brown hare (L. europaeus) mtDNA non-coding regions, together with a species- and population-level analysis of tandem repeat variation. Mountain hare short tandem repeats (SRs) as well as other analyzed hare species consist of two conserved 10 bp motifs, with only brown hares exhibiting a single, more variable motif. Long tandem repeats (LRs) also differ in sequence and copy number between species. Mountain hares have four to seven LRs, median value five, while brown hares exhibit five to nine LRs, median value six. Interestingly, introgressed mountain hare mtDNA in brown hares obtained an intermediate LR length distribution, with median copy number being the same as with conspecific brown hare mtDNA. In contrast, transfer of brown hare mtDNA into cultured mtDNA-less mountain hare cells maintained the original LR number, whereas the reciprocal transfer caused copy number instability, suggesting that cellular environment rather than the nuclear genomic background plays a role in the LR maintenance. Due to their dynamic nature and separation from other known conserved sequence elements on the non-coding region of hare mitochondrial genomes, the tandem repeat elements likely to represent signatures of ancient genetic rearrangements. clarifying the nature and dynamics of these rearrangements may shed light on the possible role of NCR repeated elements in mitochondria and in species evolution.

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.

Single-cell RNA sequencing and UPHLC-MS/MS targeted metabolomics offer new insights into the etiological basis for male cattle-yak sterility

The variety of species can be efficiently increased by interspecific hybridization. However, because the males in the hybrid progeny are usually sterile, this heterosis cannot be employed when other cattle and yaks are hybridized. While some system-level studies have sought to explore the etiological basis for male cattle-yak sterility, no systematic cellular analyses of this phenomenon have yet been performed. Here, single-cell RNA sequencing and UPHLC-MS/MS targeted metabolomics methods were used to study the differences in testicular tissue between 4-year-old male yak and 4-year-old male cattle-yak, providing new and comprehensive insights into the causes of male cattle-yak sterility. Cattle-yak testes samples detected 6 somatic cell types and one mixed germ cell type. Comparisons of these cell types revealed the more significant differences in Sertoli cells (SCs) and [Leydig cells and myoid cells (LCs_MCs)] between yak and cattle-yak samples compared to other somatic cell clusters. Even though the LCs and MCs from yaks and cattle-yaks were derived from the differentiation of the same progenitor cells, a high degree of overlap between LCs and MCs was observed in yak samples. Still, only a small overlap between LCs and MCs was observed in cattle-yak samples. Functional enrichment analyses revealed that genes down-regulated in cattle-yak SCs were primarily enriched in biological activity, whereas up-regulated genes in these cells were enriched for apoptotic activity. Furthermore, the genes of up-regulated in LCs_MCs of cattle-yak were significantly enriched in enzyme inhibitor and molecular function inhibitor activity. On the other hand, the genes of down-regulated in these cells were enriched for signal receptor binding, molecular function regulation, positive regulation of biological processes, and regulation of cell communication activity. The most significant annotated differences between yak and cattle-yak LCs_MCs were associated with cell-to-cell communication. While yak LCs_MCs regulated spermatogenic cells at spermatogonia, spermatocyte, and spermatid levels, no such relationships were found between cattle-yak LCs_MCs and germ cells. This may suggest that the somatic niche in male cattle-yak testes is a microenvironment that is ultimately not favorable for spermatogenesis.

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.

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.

Measurement(s)	whole genome sequencing
Technology Type(s)	Illumina HiSeq X
Factor Type(s)	Species
Sample Characteristic - Organism	Lepus timidus • Lepus europaeus • Lepus timidus sylvaticus
Sample Characteristic - Location	Sweden