Cytochrome b phylogeny of North American hares and jackrabbits (Lepus, lagomorpha) and the effects of saturation in outgroup taxa

Seventy-eight entire mitochondrial genomes and nuclear rRNA genes provide insight into the phylogeny of the hard ticks, particularly the Haemaphysalis species, Africaniella transversale and Robertsicus elaphensis

Hoogstraal and Kim (1985) proposed from morphology, three groups of Haemaphysalis subgenera: (i) the "structurally advanced"; (ii) the "structurally intermediate"; and (iii) the "structurally primitive" subgenera. Nuclear gene phylogenies, however, did not indicate monophyly of these morphological groups but alas, only two mitochondrial (mt) genomes from the "structurally intermediate" subgenera had been sequenced. The phylogeny of Haemaphysalis has not yet been resolved. We aimed to resolve the phylogeny of the genus Haemaphysalis, with respect to the subgenus Alloceraea. We presented 15 newly sequenced and annotated mt genomes from 15 species of ticks, five species of which have not been sequenced before, and four new 18S rRNA and 28S rRNA nuclear gene sequences. Our datasets were constructed from 10 mt protein-coding genes, cox1, and the 18S and 28S nuclear rRNA genes. We found a 132-bp insertion between tRNA-Glu (E) gene and the nad1 gene in the mt genome of Haemaphysalis (Alloceraea) inermis that resembles insertions in H. (Alloceraea) kitaokai and Rhipicephalus (Boophilus) geigyi. Our mt phylogenies had the three species of Amblyomma (Aponomma) we sequenced embedded in the main clade of Amblyomma: Am. (Aponomma) fimbriatum, Am. (Aponomma) gervaisi and Am. (Aponomma) latum. This is further support for the hypothesis that the evolution of eyes appears to have occurred in the most-recent-common-ancestor of Amblyocephalus (i.e. Amblyomminae plus Rhipicephalinae) and that eyes were subsequently lost in the most-recent-common-ancestor of the subgenus Am. (Aponomma). Either Africaniella transversale or Robertsicus elaphensis, or perhaps Af. transversale plus Ro. elaphensis, appear to be the sister-group to the rest of the metastriate Ixodida. Our cox1 phylogenies did not indicate monophyly of the "structurally primitive", "structurally intermediate" nor the "structurally advanced" groups of Haemaphysalis subgenera. Indeed, the subgenus Alloceraea may be the only monophyletic subgenus of the genus Haemaphysalis sequenced thus far. All of our mt genome and cox1 phylogenies had the subgenus Alloceraea in a clade that was separate from the rest of the Haemaphysalis ticks. If Alloceraea is indeed the sister to the rest of the Haemaphysalis subgenera this would resonate with the argument of Hoogstraal and Kim (1985), that Alloceraea was a subgenus of "primitive" Haemaphysalis. Alectorobius capensis from Japan had a higher genetic-identity to A. sawaii, which was also from Japan, than to the A. capensis from South Africa. This indicates that A. capensis from Japan may be a cryptic species with respect to the A. capensis from South Africa.

Lagomorpha as a Model Morphological System

Due to their global distribution, invasive history, and unique characteristics, European rabbits are recognizable almost anywhere on our planet. Although they are members of a much larger group of living and extinct mammals [Mammalia, Lagomorpha (rabbits, hares, and pikas)], the group is often characterized by several well-known genera (e.g., Oryctolagus, Sylvilagus, Lepus, and Ochotona). This representation does not capture the extraordinary diversity of behavior and form found throughout the order. Model organisms are commonly used as exemplars for biological research, but there are a limited number of model clades or lineages that have been used to study evolutionary morphology in a more explicitly comparative way. We present this review paper to show that lagomorphs are a strong system in which to study macro- and micro-scale patterns of morphological change within a clade that offers underappreciated levels of diversity. To this end, we offer a summary of the status of relevant aspects of lagomorph biology.

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

The evolutionary history of the Cape hare (Lepus capensis sensu lato): insights for systematics and biogeography

Inferring the phylogeography of species with large distributions helps deciphering major diversification patterns that may occur in parallel across taxa. Here, we infer the evolutionary history of the Cape hare, Lepus capensis sensu lato, a species distributed from southern Africa to Asia, by analyzing variation at 18 microsatellites and 9 DNA (1 mitochondrial and 8 nuclear) sequenced loci, from field and museum-collected samples. Using a combination of assignment and coalescent-based methods, we show that the Cape hare is composed of five evolutionary lineages, distributed in distinct biogeographic regions-north-western Africa, eastern Africa, southern Africa, the Near East and the Arabian Peninsula. A deep phylogenetic break possibly dating to the Early Pleistocene was inferred between the African and Asian L. capensis groups, and the latter appear more closely related to other Eurasian hare species than to African Cape hares. The inferred phylogeographic structure is shared by numerous taxa distributed across the studied range, suggesting that environmental changes, such as the progressive aridification of the Saharo-Arabian desert and the fluctuations of savannah habitats in Sub-Saharan Africa, had comparable impacts across species. Fine-scale analyses of the western Sahara-Sahel populations showed rich fragmentation patterns for mitochondrial DNA but not for microsatellites, compatible with the environmental heterogeneity of the region and female philopatry. The complex evolutionary history of L. capensis sensu lato, which possibly includes interspecific gene flow, is not reflected by taxonomy. Integrating evolutionary inference contributes to an improved characterization of biodiversity, which is fundamental to foster the conservation of relevant evolutionary units.

Mitochondrial and nuclear DNA reveals reticulate evolution in hares (Lepus spp., Lagomorpha, Mammalia) from Ethiopia

For hares (Lepus spp., Leporidae, Lagomorpha, Mammalia) from Ethiopia no conclusive molecular phylogenetic data are available. To provide a first molecular phylogenetic model for the Abyssinian Hare (Lepus habessinicus), the Ethiopian Hare (L. fagani), and the Ethiopian Highland Hare (L. starcki) and their evolutionary relationships to hares from Africa, Eurasia, and North America, we phylogenetically analysed mitochondrial ATPase subunit 6 (ATP6; n = 153 / 416bp) and nuclear transferrin (TF; n = 155 / 434bp) sequences of phenotypically determined individuals. For the hares from Ethiopia, genotype composition at twelve microsatellite loci (n = 107) was used to explore both interspecific gene pool separation and levels of current hybridization, as has been observed in some other Lepus species. For phylogenetic analyses ATP6 and TF sequences of Lepus species from South and North Africa (L. capensis, L. saxatilis), the Anatolian peninsula and Europe (L. europaeus, L. timidus) were also produced and additional TF sequences of 18 Lepus species retrieved from GenBank were included as well. Median joining networks, neighbour joining, maximum likelihood analyses, as well as Bayesian inference resulted in similar models of evolution of the three species from Ethiopia for the ATP6 and TF sequences, respectively. The Ethiopian species are, however, not monophyletic, with signatures of contemporary uni- and bidirectional mitochondrial introgression and/ or shared ancestral polymorphism. Lepus habessinicus carries mtDNA distinct from South African L. capensis and North African L. capensis sensu lato; that finding is not in line with earlier suggestions of its conspecificity with L. capensis. Lepus starcki has mtDNA distinct from L. capensis and L. europaeus, which is not in line with earlier suggestions to include it either in L. capensis or L. europaeus. Lepus fagani shares mitochondrial haplotypes with the other two species from Ethiopia, despite its distinct phenotypic and microsatellite differences; moreover, it is not represented by a species-specific mitochondrial haplogroup, suggesting considerable mitochondrial capture by the other species from Ethiopia or species from other parts of Africa. Both mitochondrial and nuclear sequences indicate close phylogenetic relationships among all three Lepus species from Ethiopia, with L. fagani being surprisingly tightly connected to L. habessinicus. TF sequences suggest close evolutionary relationships between the three Ethiopian species and Cape hares from South and North Africa; they further suggest that hares from Ethiopia hold a position ancestral to many Eurasian and North American species.

Comparative anatomy and systematic implications of the turbinal skeleton in Lagomorpha (Mammalia)

In order to elucidate the systematic relevance of the turbinal skeleton in Lagomorpha the ethmoidal regions of 6 ochotonid, 21 leporid, and 2 outgroup species (Sciurus vulgaris, Tupaia sp.) species were investigated by high-resolution computed tomography (μCT). Number and shape of turbinals correspond to major clades and to several genera. All Lagomorpha under study have a deeply excavated nasoturbinal that is continuous with the lamina semicircularis; a feature likely to be an autapomorphy of lagomorphs. In particular, the olfactory turbinals (frontoturbinals, ethmoturbinals, and interturbinals) provide new systematic information. The plesiomorphic lagomorph pattern comprises two frontoturbinals, three ethmoturbinals, and one interturbinal between ethmoturbinal I and II. Ochotonidae are derived from the lagomorph goundplan by loss of ethmoturbinal III; an interturbinal between the two frontoturbinals is an autapomorphy of Leporidae. Pronolagus is apomorphic in having a very slender first ethmoturbinal, but shows a puzzling pattern in decreasing the number of turbinals. Pronolagus rupestris and Romerolagus diazi have independently reduced their turbinals to just two fronto- and two ethmoturbinals, which is the lowest number among the sampled lagomorphs. In contrast, the more derived leporid genera under study (Oryctolagus, Caprolagus, Sylvilagus, and Lepus) show a tendency to increase the number of turbinals, either by developing an ethmoturbinal IV (Caprolagus hispidus, Lepus arcticus) or by additional interturbinals. Intraspecific variation was investigated in Ochotona alpina, Oryctolagus cuniculus, and Lepus europaeus and is restricted to additional interturbinals in the frontoturbinal recess of the two leporids.

Northern host-parasite assemblages: history and biogeography on the borderlands of episodic climate and environmental transition

Diversity among assemblages of mammalian hosts and parasites in northern terrestrial ecosystems was structured by a deep history of biotic and abiotic change that overlies a complex geographic arena. Since the Pliocene, Holarctic ecosystems assembled in response to shifting climates (glacial and interglacial stages). Cycles of episodic dispersal/isolation and diversification defined northern diversity on landscape to regional scales. Episodes of geographic expansion and colonisation linked Eurasia and North America across Beringia and drove macroevolutionary structure of host and parasite associations. Asynchronous dispersal from centres of origin in Eurasia into the Nearctic resulted in gradients in parasite diversity in the carnivoran, lagomorph, rodent and artiodactyl assemblages we reviewed. Recurrent faunal interchange and isolation in conjunction with episodes of host colonisation have produced a mosaic structure for parasite faunas and considerable cryptic diversity among nematodes and cestodes. Mechanisms of invasion and geographic colonisation leading to the establishment of complex faunal assemblages are equivalent in evolutionary and ecological time, as demonstrated by various explorations of diversity in these high-latitude systems. Our ability to determine historical responses to episodic shifts in global climate may provide a framework for predicting the cascading effects of contemporary environmental change.

Reticulate evolution: frequent introgressive hybridization among Chinese hares (genus lepus) revealed by analyses of multiple mitochondrial and nuclear DNA loci

BACKGROUND

Interspecific hybridization may lead to the introgression of genes and genomes across species barriers and contribute to a reticulate evolutionary pattern and thus taxonomic uncertainties. Since several previous studies have demonstrated that introgressive hybridization has occurred among some species within Lepus, therefore it is possible that introgressive hybridization events also occur among Chinese Lepus species and contribute to the current taxonomic confusion.

RESULTS

Data from four mtDNA genes, from 116 individuals, and one nuclear gene, from 119 individuals, provides the first evidence of frequent introgression events via historical and recent interspecific hybridizations among six Chinese Lepus species. Remarkably, the mtDNA of L. mandshuricus was completely replaced by mtDNA from L. timidus and L. sinensis. Analysis of the nuclear DNA sequence revealed a high proportion of heterozygous genotypes containing alleles from two divergent clades and that several haplotypes were shared among species, suggesting repeated and recent introgression. Furthermore, results from the present analyses suggest that Chinese hares belong to eight species.

CONCLUSION

This study provides a framework for understanding the patterns of speciation and the taxonomy of this clade. The existence of morphological intermediates and atypical mitochondrial gene genealogies resulting from frequent hybridization events likely contribute to the current taxonomic confusion of Chinese hares. The present study also demonstrated that nuclear gene sequence could offer a powerful complementary data set with mtDNA in tracing a complete evolutionary history of recently diverged species.

The ubiquitous mountain hare mitochondria: multiple introgressive hybridization in hares, genus Lepus

Climatic oscillations during the glaciations forced dramatic changes in species distributions, such that some presently temperate regions were alternately occupied by temperate and arctic species. These species could have met and hybridized during climatic transitions. This phenomenon happened for three hare species present in Iberia (Lepus granatensis, Lepus europaeus and Lepus castroviejoi), which display high frequencies of mitochondrial DNA (mtDNA) from Lepus timidus, an arctic/boreal species presently extinct in Iberia. Here, we extend our previous geographical survey to determine whether the distribution of this mtDNA lineage extends beyond the northern half of the Iberian Peninsula, where it is found at high frequencies. We also review the taxonomy, distribution and molecular phylogeny of the genus Lepus. The phylogenetic inference reveals the presence of L. timidus-like mtDNA in several other hare species in Asia and North America, suggesting that the mitochondrial introgression observed in Iberia might be generalized. Comparison with the available nuclear gene phylogenies suggests that introgression could have happened repeatedly, possibly during different climatic transitions. We discuss demographic and adaptive scenarios that could account for the repetition in time and space of this spectacular phenomenon and suggest ways to improve our understanding of its determinants and consequences. Such high levels of introgressive hybridization should discourage attempts to revise hare taxonomy based solely on mtDNA.

Functional specialisation of the thoracic limb of the hare (Lepus europeus)

We provide quantitative anatomical data on the muscle-tendon architecture of the hare thoracic limb (specifically muscle mass, fascicle length, pennation angle, tendon mass and length). In addition, moment arms of major thoracic limb muscles were measured. Maximum isometric force and power of muscles, the moment of force about a joint, and tendon stress and strain were estimated. Data are compared with those from other cursorial mammals. The thoracic limb of the hare consists predominantly of extrinsic musculature with long parallel fascicles, specialised for generating force over a large range. A large shoulder flexor/elbow extensor muscle mass is present, in particular Triceps brachii. The pennate nature of the long head of this muscle suggests it has an important role in stabilising the elbow joint during stance, whilst moment arm curves suggest that it may also play a role in initiating shoulder flexion. In addition, Supraspinatus and Infraspinatus are capable of generating high forces, potentially to stabilise the shoulder joint during the stance phase of locomotion. Supraspinatus may in addition play an important role in forelimb protraction. The Subscapularis muscle was capable of generating surprisingly high forces, suggesting that the hare must be able to withstand/produce high forces during activities that need medio-lateral stability, such as turning. Distally, tendons were relatively short, showing little potential for elastic energy storage when compared with both their pelvic limb counterparts and their equivalents in the horse thoracic limb. Thus, a 'stiffer' thoracic limb may be beneficial in terms of behaving like a strut, simply supporting and deflecting the body during high-speed running. This more distal/less proximal distribution of limb mass is also likely to be important in retaining the manipulative/adaptive/non-locomotor capabilities of the limb.

Resolving deep phylogenetic relationships in salamanders: analyses of mitochondrial and nuclear genomic data

Phylogenetic relationships among salamander families illustrate analytical challenges inherent to inferring phylogenies in which terminal branches are temporally very long relative to internal branches. We present new mitochondrial DNA sequences, approximately 2,100 base pairs from the genes encoding ND1, ND2, COI, and the intervening tRNA genes for 34 species representing all 10 salamander families, to examine these relationships. Parsimony analysis of these mtDNA sequences supports monophyly of all families except Proteidae, but yields a tree largely unresolved with respect to interfamilial relationships and the phylogenetic positions of the proteid genera Necturus and Proteus. In contrast, Bayesian and maximum-likelihood analyses of the mtDNA data produce a topology concordant with phylogenetic results from nuclear-encoded rRNA sequences, and they statistically reject monophyly of the internally fertilizing salamanders, suborder Salamandroidea. Phylogenetic simulations based on our mitochondrial DNA sequences reveal that Bayesian analyses outperform parsimony in reconstructing short branches located deep in the phylogenetic history of a taxon. However, phylogenetic conflicts between our results and a recent analysis of nuclear RAG-1 gene sequences suggest that statistical rejection of a monophyletic Salamandroidea by Bayesian analyses of our mitochondrial genomic data is probably erroneous. Bayesian and likelihood-based analyses may overestimate phylogenetic precision when estimating short branches located deep in a phylogeny from data showing substitutional saturation; an analysis of nucleotide substitutions indicates that these methods may be overly sensitive to a relatively small number of sites that show substitutions judged uncommon by the favored evolutionary model.

Biochemical Genetic Relationships Among Tunisian Hares (Lepus sp.), South African Cape Hares (L. capensis), and European Brown Hares (L. europaeus)

Tunisian hares (n = 45), currently assigned to Lepus capensis, were assayed for allelic variation at 40 allozyme loci, and allele frequencies at 32 loci were directly compared with earlier data of South African cape hares (L. capensis, n = 9) and European brown hares (L. europaeus, n = 244) to reveal genetic relationships among them. European mountain hares (L. timidus, n = 200) were used for outgroup comparison. In the Tunisian hares 27.5% of the loci were polymorphic with 2-4 alleles. Among all alleles at polymorphic loci, 15.1% occurred exclusively in Tunisian hares, 5.7% exclusively in cape hares, and 7.5% exclusively in brown hares at low frequencies. Not a single locus showed alternately fixed alleles between the samples of the L. capensis/L. europaeus complex. Levels of absolute and relative genetic differentiation among the samples of the L. capensis/L. europaeus complex were low, relative to pairwise comparisons involving mountain hares. Diverse cluster analyses and multidimensional scaling of various pairwise genetic distance matrices concordantly grouped Tunisian hares with brown hares, and South African cape hares clustered only slightly farther apart, whereas mountain hares were distinctly separate. These results suggest regionally distinct phylogenetic units within an overall cohesive gene pool in the L. capensis/L. europaeus complex, supporting Petter's view that all North African hares belong to L. capensis except for one local population of savanna hares, and that cape hares and brown hares are conspecific.

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

Hares on ice: phylogeography and historical demographics of Lepus arcticus, L. othus, and L. timidus (Mammalia: Lagomorpha)

Phylogeographical investigations of arctic organisms provide spatial and temporal frameworks for interpreting the role of climate change on biotic diversity in high-latitude ecosystems. Phylogenetic analyses were conducted on 473 base pairs of the mitochondrial control region in 192 arctic hares (Lepus arcticus, Lepus othus, Lepus timidus) and two individual Lepus townsendii. The three arctic hare species are closely related. All L. othus individuals form one well-supported clade, L. arcticus individuals form two well-supported clades, and L. timidus individuals are scattered throughout the phylogeny. Arctic hare distribution was altered dramatically following post-Pleistocene recession of continental ice sheets. We tested for genetic signatures of population expansion for hare populations now found in deglaciated areas. Historical demographic estimates for 12 arctic hare populations from throughout their range indicate that L. arcticus and L. othus persisted in two separate North American arctic refugia (Beringia and High Canadian Arctic) during glacial advances of the Pleistocene, while the high genetic diversity in L. timidus likely reflects multiple Eurasian refugia.

Invasion from the cold past: extensive introgression of mountain hare (Lepus timidus) mitochondrial DNA into three other hare species in northern Iberia

Mitochondrial DNA introgression from Lepus timidus into Lepus granatensis and Lepus europaeus was recently reported in Iberia, although L. timidus presumably retreated from this region at the end of the last ice age. Here we assess the extent of this ancient mtDNA introgression by RFLP analysis of 695 specimens representing the three hare species present in Iberia. The introgressed L. timidus lineage was found in 23 of the 37 populations sampled. It is almost fixed in L. europaeus across its Iberian range in the Pyrenean foothills, and in L. granatensis, which occupies the rest of the peninsula, it is predominant in the north and gradually disappears further south. We also found it in Lepus castroviejoi, a species endemic to Cantabria. Multiple hybridizations and, potentially, a selective advantage for the L. timidus lineage can explain the remarkable taxonomic and geographical range of this mitochondrial introgression.

Coevolution and biogeography among Nematodirinae (Nematoda: Trichostrongylina) Lagomorpha and Artiodactyla (Mammalia): exploring determinants of history and structure for the northern fauna across the holarctic

Nematodes of the subfamily Nematodirinae are characteristic components of a Holarctic fauna. The topology of a generic-level phylogenetic hypothesis, patterns of diversity, and geographic distributions for respective nematode taxa in conjunction with data for host occurrence are consistent with primary distributions determined across Beringia for species of Murielus, Rauschia, Nematodirus, and Nematodirella. Ancestral hosts are represented by Lagomorpha, with evidence for a minimum of 1 host-switching-event and subsequent radiation in the Artiodactyla. Diversification may reflect vicariance of respective faunas along with episodic or cyclical range expansion and isolation across Beringia during the late Tertiary and Quaternary. Secondarily, species of Nematodirus attained a distribution in the Neotropical region with minimal diversification of an endemic fauna represented by Nematodirus molini among tayassuids, Nematodirus lamae among camelids and Nematodirus urichi in cervids during the Pleistocene. Nematodirines are a core component of an Arctic-Boreal fauna of zooparasitic nematodes (defined by latitude and altitude) adapted to transmission in extreme environments characterized by seasonally low temperatures and varying degrees of desiccation. The history and distribution of this fauna is examined in the context of biotic and abiotic determinants for geographic colonization and host switching with an exploration of predicted responses of complex host-parasite systems to ecological perturbation under a regime of global climate change.

Utility of nuclear DNA intron markers at lower taxonomic levels: Phylogenetic resolution among nine Tragelaphus spp

Phylogenetic relationships among the nine spiral-horn antelope species of the African bovid tribe Tragelaphini are controversial. In particular, mitochondrial DNA sequencing studies are not congruent with previous morphological investigations. To test the utility of nuclear DNA intron markers at lower taxonomic levels and to provide additional data pertinent to tragelaphid evolution, we sequenced four nuclear DNA segments (MGF, PRKCI, SPTBN, and THY) and combined these data with mitochondrial DNA sequences from three genes (cytochrome b, 12S rRNA, and 16S rRNA). Our molecular supermatrix comprised 4682 characters which were analyzed independently and in combination. Parsimony and model based phylogenetic analyses of the combined nuclear DNA data are congruent with those derived from the analysis of mitochondrial gene sequences. The corroboration between nuclear and mtDNA gene trees reject the possibility that genetic processes such as lineage sorting, gene duplication/deletion and hybrid speciation account for the conflict evident in the previously published phylogenies. It suggests rather that the morphological characters used to delimit the Tragelaphid species are subject to convergent evolution. Divergence times among species, calculated using a relaxed Bayesian molecular clock, are consistent with hypotheses proposing that climatic oscillations and their impact on habitats were the major forces driving speciation in the tribe Tragelaphini.

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.

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.

Historical biogeography at the crossroads of the northern continents: molecular phylogenetics of red-backed voles (Rodentia: Arvicolinae)

Evolutionary relationships of red-backed voles and their relatives were examined and used to test biogeographic hypotheses. Sequences of the mitochondrial cytochrome b gene were obtained for 25 individuals representing Alticola macrotis, Clethrionomys californicus, C. gapperi, C. glareolus, C. rutilus, and C. rufocanus. These were combined with 21 partial sequences from GenBank for C. regulus, C. rex, C. rufocanus, C. rutilus, Eothenomys imaizumii, E. melanogaster, Phaulomys andersoni, and P. smithii. Complete sequences of three species of Microtus (M. montanus, M. oeconomus, and M. pennsylvanicus), representative species of other arvicoline genera (Myopus, Synaptomys, Arvicola, Ellobius, Ondatra, Lemmus, Dicrostonyx, and Phenacomys), and a sigmodontine representative (Peromyscus) were included as outgroups. We used maximum parsimony, maximum likelihood, distance, and Bayesian based methods and conducted statistical tests on proposed hypotheses of phylogenetic relationships and biogeographic histories. A close relationship of species representing the genera Alticola, Clethrionomys, and Eothenomys was supported (Clethrionomyini); however, the genus Clethrionomys was paraphyletic with respect to both Alticola and Eothenomys. Three major clades were identified as Asian (Eothenomys andersoni, E. smithii, C. rex, C. regulus, and C. rufocanus), Trans-beringian (Alticola macrotis, C. californicus, C. gapperi, C. glarelolus, and C. rutilus), and Taiwanese (E. melanogaster). These results are consistent with the fossil record which indicates an initial diversification in Asia followed by colonization of the Nearctic on at least two occasions. The holarctic species, C. rutilus, appears to have either reinvaded Asia from North America or colonized North America more recently (late Pleistocene) than the two species of Clethrionomys (C. gapperi and C. californicus) that are endemic to North America (early to mid-Pleistocene). Finally, C. gapperi, appears to be comprised of an eastern and a western species, the former with affinities to the Asian C. glareolus and the latter more closely related to C. californicus.

The origin and evolution of seahorses (genus Hippocampus): a phylogenetic study using the cytochrome b gene of mitochondrial DNA

Phylogenetic relationships among 93 specimens of 22 species of seahorses (genus Hippocampus) from the Atlantic and Indo-Pacific Oceans were analysed using cytochrome b gene sequence data. A maximum sequence divergence of 23.2% (Kimura 2-parameter model) suggests a pre-Tethyan origin for the genus. Despite a greater number of seahorse species in the Indo-Pacific than in the Atlantic Ocean, there was no compelling genetic evidence to support an Indo-Pacific origin for the genus Hippocampus. The phylogenetic data suggest that high diversity in the Indo-Pacific results from speciation events dating from the Pleistocene to the Miocene, or earlier. Both vicariance and dispersal events in structuring the current global distribution of seahorses. The results suggested that several species designations need re-evaluating, and further phylogeographic studies are required to determine patterns and processes of seahorse dispersal.

Molecular genetics of Rhabdomys pumilio subspecies boundaries: mtDNA phylogeography and karyotypic analysis by fluorescence in situ hybridization

The phylogeography of the African four-striped mouse, Rhabdomys pumilio, was investigated using complete sequences of the mtDNA cytochrome b gene (1140 bp) and a combination of fluorescence in situ hybridization (FISH) and conventional cytogenetic banding techniques (G- and C-banding). Two cytotypes (2n=46 and 2n=48) were identified by cytogenetic analysis. There is no evidence of diploid number variation within populations, difference in gross chromosome morphology or of subtle interchromosomal rearrangements at levels detected by ZOO-FISH. Analysis of the mtDNA cytochrome b resulted in two major lineages that correspond roughly to the xeric and mesic biotic zones of southern Africa. One mtDNA clade comprises specimens with 2n=48 and the other representatives of two cytotypes (2n=48 and 2n=46). The mean sequence divergence (12%, range 8.3-15.6%) separating the two mtDNA clades is comparable to among-species variation within murid genera suggesting their recognition as distinct species, the prior names for which would be R. dilectus and R. pumilio. Low sequence divergences and the diploid number dichotomy within the mesic lineage support the recognition of two subspecies corresponding to R. d. dilectus (2n=46) and R. d. chakae (2n=48). Our data do not support subspecific delimitation within the nominate, R. pumilio. Molecular dating places cladogenesis of the two putative species at less than five million years, a period characterised by extensive climatic oscillations which are thought to have resulted in habitat fragmentation throughout much of the species range.

Ancient introgression of Lepus timidus mtDNA into L. granatensis and L. europaeus in the Iberian Peninsula

A 587 bp fragment of cytochrome b sequences from 90 individuals of 15 hare (Lepus) species and two outgroups were phylogenetically analysed and compared to an analysis derived from 474 bp sequences of the nuclear transferrin gene. Mountain hare (Lepus timidus) type mtDNA was observed in L. granatensis and L. europaeus from the Iberian Peninsula, far away from the extant distributional range of L. timidus. In addition to these two hare species, other hare species may also contain mtDNA from L. timidus. This species may have introgressed with other species of Lepus that occur within its present range, or where fossils indicate its historical presence during glacial periods. L. timidus mtDNA is common in the northern part of the L. granatensis range. Finally, we reassessed the phylogenetic relationships of the five European hare species based on both mitochondrial and nuclear DNA sequences.

The influence of climate on the basal metabolic rate of small mammals: a slow-fast metabolic continuum

The influence of climate (mean annual rainfall, rainfall variability, ambient temperature, T(a)) on the basal metabolic rate (BMR) of 267 small mammals (<1 kg) from six zoogeographical zones was investigated using conventional and phylogenetically independent data (linear contrasts). All climate variables varied between zones, as did BMR and body temperature ( T(b)), but not thermal conductance. Holarctic zones were more seasonal and colder, but rainfall was less variable, than non-Holarctic zones. In general, the BMR was most strongly influenced by body mass, followed by T(a) and the rainfall variables. However, there was significant variation in the strength of these relationships between zones. BMR and T(b) increased with latitude, and mass-independent BMR and T(b) were positively correlated. The latter relationship offers evidence of a slow-fast metabolic continuum in small mammals. The fast end of the continuum (high BMR) is associated with the highest latitudes where BMR is most strongly influenced by T(a) and mean annual rainfall (i.e. mean productivity). The slow end of the continuum (low BMR) is associated with the semi-tropics, low productivity zones, and climatically unpredictable zones, such as deserts. Here rainfall variability has the strongest influence on BMR after body size. The implications of a slow-fast metabolic continuum are discussed in terms of various models associated with the evolution of BMR, such as the aerobic capacity models and the "energetic definition of fitness" models.

A phylogeny of the megapodes (Aves: Megapodiidae) based on nuclear and mitochondrial DNA sequences

DNA sequences from the first intron of the nuclear gene rhodopsin (RDP1) and from the mitochondrial gene ND2 were used to construct a phylogeny of the avian family Megapodiidae. RDP1 sequences evolved about six times more slowly than ND2 and showed less homoplasy, substitution bias, and rate heterogeneity across sites. Analysis of RDP1 produced a phylogeny that was well resolved at the genus level, but RDP1 did not evolve rapidly enough for intrageneric comparisons. The ND2 phylogeny resolved intrageneric relationships and was congruent with the RDP1 phylogeny except for a single node: this node was the only aspect of tree topology sensitive to weighting in parsimony analyses. Despite differences in sequence evolution, RDP1 and ND2 contained congruent phylogenetic signal and were combined to produce a phylogeny that reflects the resolving power of both genes. This phylogeny shows an early split within the megapodes, leading to two major clades: (1) Macrocephalon and the mound-building genera Talegalla, Leipoa, Aepypodius, and Alectura, and (2) Eulipoa and Megapodius. It differs significantly from previous hypotheses based on morphology but is consistent with affiliations suggested by a recent study of parasitic chewing lice.

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

We determined mitochondrial 12S ribosomal RNA (rRNA) and cytochrome b (cyt b) gene sequences in three leporid species of Japan, the Amami rabbit Pentalagus furnessi from the Ryukyu Islands, the Japanese hare Lepus brachyurus from Honshu, and a Japanese form of the mountain hare Lepus timidus ainu from Hokkaido. We compared the sequences with those of other taxa of leporids available in databases. Phylogenetic trees of the 12S rRNA gene sequences indicated that the lineage of P. furnessi diversified during the generic radiation of the leporids at an ancient time, which was estimated to have been the middle Miocene. Cyt-b gene trees revealed that the lineage of L. brachyurus branched off at an early stage in the speciation of Lepus, probably at the beginning of the Pliocene. The cyt b sequences of L. t. ainu were somewhat distinct from those of continental conspecific populations; this lineage divergence is likely to have occurred during the middle or late Pleistocene. The results show that the three regions of the Japanese archipelago, Ryukyu, Honshu-Shikoku-Kyushu, and Hokkaido, now preserve their own leporid taxa, each with a different extent of genetic endemicity. It is possible that the zoogeographic traits of the Japanese leporids are a consequence of the evolutionary dynamics of leporids in East Asia, in that the radiation centers of leporids are likely to have shifted from tropical, through temperate, to arctic zones.

Molecular systematics and evolution of Regina and the thamnophiine snakes

Snakes of the tribe Thamnophiini represent an ecologically important component of the herpetofauna in a range of habitats across North America. Thamnophiines are the best-studied colubrids, yet little is known of their systematic relationships. A molecular phylogenetic study of 32 thamnophiine species using three complete mitochondrial genes (cytochrome b, NADH dehydrogenase subunit 2, and 12S ribosomal DNA) recovered a well-supported phylogeny with three major clades: a garter snake group, a water snake group, and a novel semifossorial group. The historically contentious genus Regina, which contains the crayfish-eating snakes, is polyphyletic. The phylogeographic pattern of Thamnophis is consistent with an hypothesis of at least one invasion of northern North America from Mexico.

Retrieval of four adaptive lineages in duiker antelope: evidence from mitochondrial DNA sequences and fluorescence in situ hybridization

Independent molecular markers (mitochondrial DNA sequences from two genes and fluorescence in situ hybridization with satellite DNA sequences as hybridization probes) were employed to investigate phylogenetic relationships among duiker antelope. When analyzed singly or taken together, the molecular and cytogenetic data allowed for the delimitation of four adaptive groups: the conservative dwarfs which are basal, a savanna specialist which groups apart from the forest duikers, the giant duikers, and the red duikers. Within the latter, a further subdivision comprising an east African and a west African red duiker clade is evident. The placement of the endangered zebra duiker and Aders' duiker remains problematic. Several of the nomenclatural divisions in current use are questioned by our results. These include the recognition of Philantomba as genus name for the blue and Maxwell's duiker and that Harvey's duiker be relegated to a subspecies of the Natal red duiker. We place our results in a biogeographic context and argue that duiker speciation has been driven predominantly by habitat fragmentation which probably led to the disruption of gene flow between geographic populations.