Molecular phylogeny of Japanese Leporidae, the Amami rabbit Pentalagus furnessi, the Japanese hare Lepus brachyurus, and the mountain hare Lepus timidus, inferred from mitochondrial DNA sequences

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 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.].

Primary and immortalized cell lines derived from the Amami rabbit (Pentalagus furnessi) and evolutionally conserved cell cycle control with CDK4 and Cyclin D1

The Amami rabbit (Pentagulus furnessi) is a dark brown-furred rabbit classified as an endangered species and only found in the Amami Islands of Japan. They are often called living fossils because they retain primitive characteristics of ancient rabbits that lived approximately 1 million years ago, such as short feet and hind legs and small ears. Although the ancient rabbit has disappeared due to the competition with European rabbit (Oryctolagus cuniculus) in the most of the Asian area, Amami rabbit survived since Amami Islands has isolated from Japan and Taiwan. Although Amari rabbit is one of the protected animals, their population decreases each year due to human activities, such as deforestation and roadkill. In this study, we collected roadkill samples of Amami rabbits and established primary and immortalized fibroblast cell lines. Combined expression of human-derived mutant Cyclin-dependent kinase 4, Cyclin D1, and hTERT allowed us to immortalize fibroblasts successfully in three individuals of Amami rabbits. The immortalized fibroblasts dramatically extended the cell culture period, when it was compared with the cell culture period of wild type cells. Furthermore, the immortalized cells maintained their normal chromosomal pattern (2n = 46). Our results suggest that cellular senescence which mainly regulated by p16-RB signaling pathway is conserved in animal evolution at least from 1 million years ago.

Contrasting phylogeographic histories between the continent and islands of East Asia: Massive mitochondrial introgression and long-term isolation of hares (Lagomorpha: Lepus)

Hares of the genus Lepus are distributed worldwide, and introgressive hybridization is thought to be pervasive among species, leading to reticulate evolution and taxonomic confusion. Here, we performed phylogeographic analyses of the following species of hare across East Asia: L. timidus, L. mandshuricus, L. coreanus, and L. brachyurus collected from far-eastern Russia, South Korea, and Japan. Nucleotide sequences of one mitochondrial DNA and eight nuclear gene loci were examined, adding sequences of hares in China from databases. All nuclear DNA analyses supported the clear separation of three phylogroups: L. timidus, L. brachyurus, and the L. mandshuricus complex containing L. coreanus. On the other hand, massive mitochondrial introgression from two L. timidus lineages to the L. mandshuricus complex was suggested in continental East Asia. The northern population of the L. mandshuricus complex was mainly associated with introgression from the continental lineage of L. timidus, possibly since the last glacial period, whereas the southern population of the L. mandshuricus complex experienced introgression from another L. timidus lineage related to the Hokkaido population, possibly before the last glacial period. In contrast to continental hares, no evidence of introgression was found in L. brachyurus in the Japanese Archipelago, which showed the oldest divergence amongst East Asian hare lineages. Our findings suggest that glacial-interglacial climate changes in the circum-Japan Sea region promoted distribution shifts and introgressive hybridization among continental hare species, while the geographic structure of the region contributed to long-term isolation of hares on the islands, preventing inter-species gene flow.

Variation of partial transferrin sequences and phylogenetic relationships among hares (Lepus capensis, Lagomorpha) from Tunisia

North African hares are currently included in cape hares, Lepus capensis sensu lato, a taxon that may be considered a superspecies or a complex of closely related species. The existing molecular data, however, are not unequivocal, with mtDNA control region sequences suggesting a separate species status and nuclear loci (allozymes, microsatellites) revealing conspecificity of L. capensis and L. europaeus. Here, we study sequence variation in the intron 6 (468 bp) of the transferrin nuclear gene, of 105 hares with different coat colour from different regions in Tunisia with respect to genetic diversity and differentiation, as well as their phylogenetic status. Forty-six haplotypes (alleles) were revealed and compared phylogenetically to all available TF haplotypes of various Lepus species retrieved from GenBank. Maximum Likelihood, neighbor joining and median joining network analyses concordantly grouped all currently obtained haplotypes together with haplotypes belonging to six different Chinese hare species and the African scrub hare L. saxatilis. Moreover, two Tunisian haploypes were shared with L. capensis, L timidus, L. sinensis, L. yarkandensis, and L. hainanus from China. These results indicated the evolutionary complexity of the genus Lepus with the mixing of nuclear gene haplotypes resulting from introgressive hybridization or/and shared ancestral polymorphism. We report the presence of shared ancestral polymorphism between North African and Chinese hares. This has not been detected earlier in the mtDNA sequences of the same individuals. Genetic diversity of the TF sequences from the Tunisian populations was relatively high compared to other hare populations. However, genetic differentiation and gene flow analyses (AMOVA, F_ST, Nm) indicated little divergence with the absence of geographically meaningful phylogroups and lack of clustering with coat colour types. These results confirm the presence of a single hare species in Tunisia, but a sound inference on its phylogenetic position would require additional nuclear markers and numerous geographically meaningful samples from Africa and Eurasia.

Adaptive Gene Loss? Tracing Back the Pseudogenization of the Rabbit CCL8 Chemokine

Studies of the process of pseudogenization have widened our understanding of adaptive evolutionary change. In Rabbit, an alteration at the second extra-cellular loop of the CCR5 chemokine receptor was found to be associated with the pseudogenization of one of its prime ligands, the chemokine CCL8. This relationship has raised questions about the existence of a causal link between both events, which would imply adaptive gene loss. This hypothesis is evaluated here by tracing back the history of the genetic modifications underlying the chemokine pseudogenization. The obtained data indicate that mutations at receptor and ligand genes occurred after the lineage split of New World Leporids versus Old World Leporids and prior to the generic split of the of Old World species studied, which occurred an estimated 8-9 million years ago. More important, they revealed the emergence, before this zoographical split, of a "slippery" nucleotide motif (CCCCGGG) at the 3' region of CCL8-exon2. Such motives are liable of generating +1G or -1G frameshifts, which could, however, be overcome by "translesion" synthesis or somatic reversion. The CCL8 pseudogenization in the Old World lineage was apparently initiated by three synapomorphic point mutations at the exon2-intron2 boundary which provide at short range premature terminating codons, independently of the reading frame imposed by the slippery motif. The presence of this motif in New World Leporids might allow verifying this scenario. The importance of CCL8-CCR5 signaling in parasite-host interaction would suggest that the CCL8 knock-out in Old World populations might be related to changes in pathogenic environment.

Colorful patterns indicate common ancestry in diverged tiger beetle taxa: Molecular phylogeny, biogeography, and evolution of elytral coloration of the genus Cicindela subgenus Sophiodela and its allies

We investigated the phylogenetic relationships among tiger beetles of the subtribe Cicindelina (=Cicindela s. lat.; Coleoptera: Cicindelidae) mainly from the Oriental and Sino-Japanese zoogeographic regions using one mitochondrial and three nuclear gene sequences to examine the position of the subgenus Sophiodela, currently classified in the genus Cicindela s. str., their biogeography, and the evolution of their brilliant coloration. The subgenus Sophiodela was not related to the other subgenera of Cicindela s. str. but was closely related to the genus Cosmodela. In addition, the Oriental genus Calochroa was polyphyletic with three lineages, one of which was closely related to Sophiodela and Cosmodela. The clade comprising Sophiodela, Cosmodela and two Calochroa species, referred to here as the Sophiodela group, was strongly supported, and most species in this clade had similar brilliant coloration. The Sophiodela group was related to the genera Calomera, Cicindela (excluding Sophiodela) and Cicindelidia, and these were related to Lophyra, Hipparidium and Calochroa, except species in the Sophiodela group. Divergence time estimation suggested that these worldwide Cicindelina groups diverged in the early Oligocene, and the Sophiodela group, which is found in the Oriental and Sino-Japanese zoogeographic regions, in the mid Miocene. Some components of the elytral pattern related to maculation and coloration in the Cicindelina taxa studied contained weak, but significant, phylogenetic signals and were partly associated with habitat types. Therefore, the brilliant coloration of the Sophiodela was related to both phylogeny and habitat adaptation, although the function of coloration needs to be studied.

Evaluating the "recovery level" of endangered species without prior information before alien invasion

For maintaining social and financial support for eradication programs of invasive species, quantitative assessment of recovery of native species or ecosystems is important because it provides a measurable parameter of success. However, setting a concrete goal for recovery is often difficult owing to lack of information prior to the introduction of invaders. Here, we present a novel approach to evaluate the achievement level of invasive predator management based on the carrying capacity of endangered species estimated using long-term monitoring data. In Amami-Oshima Island, Japan, where the eradication project of introduced small Indian mongoose is ongoing since 2000, we surveyed the population densities of four endangered species threatened by the mongoose (Amami rabbit, the Otton frog, Amami tip-nosed frog, and Amami Ishikawa's frog) at four time points ranging from 2003 to 2011. We estimated the carrying capacities of these species using the logistic growth model combined with the effects of mongoose predation and environmental heterogeneity. All species showed clear tendencies toward increasing their density in line with decreased mongoose density, and they exhibited density-dependent population growth. The estimated carrying capacities of three endangered species had small confidence intervals enough to measure recovery levels by the mongoose management. The population density of each endangered species has recovered to the level of the carrying capacity at about 20–40% of all sites, whereas no individuals were observed at more than 25% of all sites. We propose that the present approach involving appropriate monitoring data of native organism populations will be widely applicable to various eradication projects and provide unambiguous goals for management of invasive species.

Ancient colonization and within-island vicariance revealed by mitochondrial DNA phylogeography of the mountain hare (Lepus timidus) in Hokkaido, Japan

We examined the phylogenetic status and history of the mountain hare Lepus timidus in and around Hokkaido using mitochondrial cytochrome b (cyt b) sequences from 158 samples from Hokkaido and 14 from Sakhalin, as well as four samples from the Korean hare, L. coreanus. The phylogenetic analysis of the cyt b sequences generated in this study and obtained from DNA databases showed the clear genetic specificity of the Hokkaido lineage as a clade. The Hokkaido lineage was estimated to have diverged from the other conspecific and L. coreanus lineages 0.46 and 0.30 million years ago (Mya), respectively. These results suggest that the common ancestor of the mitochondrial lineage in Hokkaido and Korea inhabited Far East Asia before colonization by the present continental lineages of L. timidus, including the Sakhalin population. We estimated the time of the most recent common ancestor of the Hokkaido population to be 0.17 Mya, and found two distinct haplogroups within the island. One group had greater genetic diversity (mean number of pairwise differences: π = 0.0188 ± 0.0108) and appears to have expanded from the west to the entire island of Hokkaido. The other had lower genetic diversity (π = 0.0038 ± 0.0037) and its distribution was concentrated in the east. These contrasting west/east trends indicate that the Hokkaido population was fragmented in the past, and then subsequently expanded. Our study suggests that Hokkaido was an important refugium for boreal species in the far eastern region, and allowed the formation of various population genetic structures within the island.

Genetic diversity at the hinge region of the unique immunoglobulin heavy gamma (IGHG) gene in leporids (Oryctolagus, Sylvilagus and Lepus)

Unlike other species, European rabbit (Oryctolagus cuniculus) possesses only one immunoglobulin gamma class. Allelic diversity at the Ig (immunoglobulin) gamma constant region encoded by the unique IGHG (immunoglobulin heavy gamma) gene is moreover much reduced. In the European rabbit, the genetic variation at IGGH hinge region is limited to a single nucleotide substitution, which causes a Met-Thr interchange at amino acid position 9 (IMGT hinge numbering). We have analysed the diversity at this region more in-depth by, (1) analysing the allelic variation in 11 breeds of domestic European rabbit (Oryctolagus cuniculus cuniculus), and (2) sequencing the gamma hinge exon in wild specimens of six species of rabbit (Oryctolagus and Sylvilagus) and hares (Lepus), including the two Oryctolagus subspecies (O. cuniculus cuniculus and O. cuniculus algirus). It appeared that among leporid species, amino acid changes occur exclusively at positions 8 and 9. However, while position 8 is occupied by either Pro or Ser residues, four different residues can occur at position 9 (Met, Thr, Pro and Leu). This variation concerns sites of potential O-glycosylation and/or proteolytic cleavage, suggesting that the underlying genetic diversity could be the outcome of selection. Preservation of the gamma hinge polymorphism in domestic stocks could therefore be important. We report here a polymerase chain reaction/restriction fragment length polymorphism protocol that has allowed the monitoring of the heterozygosity levels at the gamma hinge in 11 breeds of domestic European rabbit.

Genetic variation at chemokine receptor CCR5 in leporids: alteration at the 2nd extracellular domain by gene conversion with CCR2 in Oryctolagus, but not in Sylvilagus and Lepus species

Whereas in its natural host (Sylvilagus sps.) the effects of myxoma virus infections are benign, in European rabbit (Oryctolagus cuniculus), it causes a highly infectious disease with very high mortality rate, known as myxomatosis. There is evidence that, as with HIV-1 virus in human, myxoma virus may use chemokine receptors such as CCR5 of the host target cell for entry and activation of pathways of immune avoidance. We have characterized and compared CCR5 genes of leporid species with different susceptibility levels to myxomatosis. The CCR5 protein of O. cuniculus differs markedly from all those known from other species. The most striking was the replacement of a specific peptide motif of the second extracellular loop (ECL2) by a motif, which in other species characterizes the CCR2 molecules. While absent in Sylvilagus and Lepus species, this CCR2 imposed CCR5-ECL2 alteration was observed in all genomes of 25 European rabbits, representing the subspecies O. cuniculus algirus and O. cuniculus cuniculus. Allelic variation at the rabbit CCR5 locus confirmed that the gene conversion predates the subspecies split (1-2 Ma).

The evolution of the immunoglobulin heavy chain variable region (IgV H ) in Leporids: an unusual case of transspecies polymorphism

In domestic rabbit (Oryctolagus cuniculus), three serological types have been distinguished at the variable domain of the antibody H chain, the so-called V(H) a allotypes a1, a2, and a3. They correspond to highly divergent allelic lineages of the V(H) 1 gene, which is the gene rabbit utilizes in more than 80% of VDJ rearrangements. The sharing of serological V(H) a markers between rabbit and snowshoe hare (Lepus americanus) has suggested that the large genetic distances between rabbit V(H) 1 alleles (9-14% nucleotide differences) can be explained by unusually long lineage persistence times (transspecies polymorphism). Because this interpretation of the serological data is uncertain, we have determined the nucleotide sequences of V(H) genes expressed in specimens of Lepus species. Two sequence groups were distinguished, one of which occurred only in hare specimen displaying serological motifs of the rabbit V(H) a-a2 allotype. Sequences of this group are part of a monophyletic cluster containing the V(H) 1 sequences of the rabbit a2 allotype. The fact that this "transspecies a2 cluster" did not include genes of other rabbit V(H) a allotypes (a1, a3, and a4) is incompatible with the existence of a common V(H) a ancestor gene within the species, and suggests that the divergence of the V(H) a lineages preceded the Lepus vs Oryctolagus split. The sequence data are furthermore compatible with the hypothesis that the V(H)a polymorphism can be two times older than the divergence time between the Lepus and Oryctolagus lineages, which was estimated at 16-24 million years.

Molecular phylogenetics and biogeography of Lepus in Eastern Asia based on mitochondrial DNA sequences

In spite of several classification attempts among taxa of the genus Lepus, phylogenetic relationships still remain poorly understood. Here, we present molecular genetic evidence that may resolve some of the current incongruities in the phylogeny of the leporids. The complete mitochondrial cytb, 12S genes, and parts of ND4 and control region fragments were sequenced to examine phylogenetic relationships among Chinese hare taxa and other leporids throughout the World using maximum parsimony, maximum likelihood, and Bayesian phylogenetic reconstruction approaches. Using reconstructed phylogenies, we observed that the Chinese hare is not a single monophyletic group as originally thought. Instead, the data infers that the genus Lepus is monophyletic with three unique species groups: North American, Eurasian, and African. Ancestral area analysis indicated that ancestral Lepus arose in North America and then dispersed into Eurasia via the Bering Land Bridge eventually extending to Africa. Brooks Parsimony analysis showed that dispersal events followed by subsequent speciation have occurred in other geographic areas as well and resulted in the rapid radiation and speciation of Lepus. A Bayesian relaxed molecular clock approach based on the continuous autocorrelation of evolutionary rates along branches estimated the divergence time between the three major groups within Lepus. The genus appears to have arisen approximately 10.76 MYA (+/-0.86 MYA), with most speciation events occurring during the Pliocene epoch (5.65+/-1.15 MYA approximately 1.12 +/- 0.47 MYA).

Evolution and biogeography of talpid moles from continental East Asia and the Japanese islands inferred from mitochondrial and nuclear gene sequences

We sequenced the cytochrome b gene from two little-studied mammal species from the highlands of Southwest China, the long-tailed mole Scaptonyx fusicaudus and the gracile shrew-like mole Uropsilus gracilis. This data was used to examine the phylogenetic relationships among 19 talpid species within the family Talpidae (Mammalia: Eulipotyphla). Cytochrome b gene trees supported a basal placement of shrew-like moles (Uropsilus) within the Talpidae, and suggested that fossorial specializations arose twice during talpid evolution. To assess the evolutionary relationships of moles endemic to this region, we additionally sequenced the 12S rRNA gene and the nuclear recombination-activating gene-1 from eight and ten East Asian taxa, respectively. Analyses of these single and concatenated data sets suggested that East Asian shrew moles diverged prior to the evolution of fossorial Eurasian moles. However, we were unable to determine whether semi-fossorial shrew moles are monophyletic. In contrast, fossorial Eurasian genera (Talpa, Mogera and Euroscaptor) were consistently found to form a monophyletic clade, with Mogera and Euroscaptor representing sister taxa. Furthermore, this fossorial clade grouped with the semi-aquatic Desmana, although with fairly low (35-62%) bootstrap support. Mogera imaizumii was found to be more closely related to M. wogura than to M. tokudae. This implies that the ancestors of these three species entered Japan from the Asian continent in this order via a series of migration events, suggesting that the Japanese Islands have played an important role in preserving mole lineages from ancient to recent times.

A molecular supermatrix of the rabbits and hares (Leporidae) allows for the identification of five intercontinental exchanges during the Miocene

The hares and rabbits belonging to the family Leporidae have a nearly worldwide distribution and approximately 72% of the genera have geographically restricted distributions. Despite several attempts using morphological, cytogenetic, and mitochondrial DNA evidence, a robust phylogeny for the Leporidae remains elusive. To provide phylogenetic resolution within this group, a molecular supermatrix was constructed for 27 taxa representing all 11 leporid genera. Five nuclear (SPTBN1, PRKCI, THY, TG, and MGF) and two mitochondrial (cytochrome b and 12S rRNA) gene fragments were analyzed singly and in combination using parsimony, maximum likelihood, and Bayesian inference. The analysis of each gene fragment separately as well as the combined mtDNA data almost invariably failed to provide strong statistical support for intergeneric relationships. In contrast, the combined nuclear DNA topology based on 3601 characters greatly increased phylogenetic resolution among leporid genera, as was evidenced by the number of topologies in the 95% confidence interval and the number of significantly supported nodes. The final molecular supermatrix contained 5483 genetic characters and analysis thereof consistently recovered the same topology across a range of six arbitrarily chosen model specifications. Twelve unique insertion-deletions were scored and all could be mapped to the tree to provide additional support without introducing any homoplasy. Dispersal-vicariance analyses suggest that the most parsimonious solution explaining the current geographic distribution of the group involves an Asian or North American origin for the Leporids followed by at least nine dispersals and five vicariance events. Of these dispersals, at least three intercontinental exchanges occurred between North America and Asia via the Bering Strait and an additional three independent dispersals into Africa could be identified. A relaxed Bayesian molecular clock applied to the seven loci used in this study indicated that most of the intercontinental exchanges occurred between 14 and 9 million years ago and this period is broadly coincidental with the onset of major Antarctic expansions causing land bridges to be exposed.

Temporal, spatial, and ecological modes of evolution of Eurasian Mus based on mitochondrial and nuclear gene sequences

We sequenced mitochondrial (cytochrome b, 12S rRNA) and nuclear (IRBP, RAG1) genes for 17 species of the Old World murine genus Mus, drawn primarily from the Eurasian subgenus Mus. Phylogenetic analysis of the newly and previously available sequences support recognition of four subgenera within Mus (Mus, Coelomys, Nannomys, and Pyromys), with an unresolved basal polytomy. Our data further indicate that the subgenus Mus contains three distinct 'species groups': (1) a Mus booduga Species Group, also including Mus terricolor and Mus fragilicauda (probably also Mus famulus); (2) a Mus cervicolor Species Group, also including Mus caroli and Mus cookii; and (3) a Mus musculus Species Group, also including Mus macedonicus, Mus spicilegus, and Mus spretus. Species diversity in Eurasian Mus is probably explicable in terms of several phases of range expansion and vicariance, and by a propensity within the group to undergo biotope transitions. IRBP and RAG1 molecular clocks for Mus date the origin of subgenera to around 5-6 mya and the origin of Species Groups within subgenus Mus to around 2-3 mya. The temporal pattern of evolution among Eurasian Mus is more complex than that within the Eurasian temperate genus Apodemus.

Locomotor Mode, Maximum Running Speed, and Basal Metabolic Rate in Placental Mammals

The locomotor performance (absolute maximum running speed [MRS]) of 120 mammals was analyzed for four different locomotor modes (plantigrade, digitigrade, unguligrade, and lagomorph-like) in terms of body size and basal metabolic rate (BMR). Analyses of conventional species data showed that the MRS of plantigrade and digitigrade mammals and lagomorphs increases with body mass, whereas that of unguligrade mammals decreases with body mass. These trends were confirmed in plantigrade mammals and lagomorphs using phylogenetically independent contrasts. Multiple regression analyses of MRS contrasts (dependent variable) as a function of body mass and BMR contrasts (predictor variables) revealed that BMR was a significant predictor of MRS in the complete data set, as well as in plantigrade and nonplantigrade mammals. However, there was severe multicollinearity in the nonplantigrade model that may influence the interpretation of these models. Although these data show mass-independent correlation between BMR and MRS, they are not necessarily indicative of a cause-effect relationship. However, the analyses do identify a negligible role of body size associated with MRS once phylogenetic and BMR effects are controlled, suggesting that the body size increase in large mammals over time (i.e., Cope's rule) can probably rule out MRS as a driving variable.