Chromosome painting refines the history of genome evolution in hares and rabbits (order Lagomorpha)

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

An Annotated Draft Genome of the Mountain Hare (Lepus timidus)

Hares (genus Lepus) provide clear examples of repeated and often massive introgressive hybridization and striking local adaptations. Genomic studies on this group have so far relied on comparisons to the European rabbit (Oryctolagus cuniculus) reference genome. Here, we report the first de novo draft reference genome for a hare species, the mountain hare (Lepus timidus), and evaluate the efficacy of whole-genome re-sequencing analyses using the new reference versus using the rabbit reference genome. The genome was assembled using the ALLPATHS-LG protocol with a combination of overlapping pair and mate-pair Illumina sequencing (77x coverage). The assembly contained 32,294 scaffolds with a total length of 2.7 Gb and a scaffold N50 of 3.4 Mb. Re-scaffolding based on the rabbit reference reduced the total number of scaffolds to 4,205 with a scaffold N50 of 194 Mb. A correspondence was found between 22 of these hare scaffolds and the rabbit chromosomes, based on gene content and direct alignment. We annotated 24,578 protein coding genes by combining ab-initio predictions, homology search, and transcriptome data, of which 683 were solely derived from hare-specific transcriptome data. The hare reference genome is therefore a new resource to discover and investigate hare-specific variation. Similar estimates of heterozygosity and inferred demographic history profiles were obtained when mapping hare whole-genome re-sequencing data to the new hare draft genome or to alternative references based on the rabbit genome. Our results validate previous reference-based strategies and suggest that the chromosome-scale hare draft genome should enable chromosome-wide analyses and genome scans on hares.

The genomic impact of historical hybridization with massive mitochondrial DNA introgression

BACKGROUND

The extent to which selection determines interspecific patterns of genetic exchange enlightens the role of adaptation in evolution and speciation. Often reported extensive interspecific introgression could be selection-driven, but also result from demographic processes, especially in cases of invasive species replacements, which can promote introgression at their invasion front. Because invasion and selective sweeps similarly mold variation, population genetics evidence for selection can only be gathered in an explicit demographic framework. The Iberian hare, Lepus granatensis, displays in its northern range extensive mitochondrial DNA introgression from L. timidus, an arctic/boreal species that it replaced locally after the last glacial maximum. We use whole-genome sequencing to infer geographic and genomic patterns of nuclear introgression and fit a neutral model of species replacement with hybridization, allowing us to evaluate how selection influenced introgression genome-wide, including for mtDNA.

RESULTS

Although the average nuclear and mtDNA introgression patterns contrast strongly, they fit a single demographic model of post-glacial invasive replacement of timidus by granatensis. Outliers of elevated introgression include several genes related to immunity, spermatogenesis, and mitochondrial metabolism. Introgression is reduced on the X chromosome and in low recombining regions.

CONCLUSIONS

General nuclear and mtDNA patterns of introgression can be explained by purely demographic processes. Hybrid incompatibilities and interplay between selection and recombination locally modulate levels of nuclear introgression. Selection promoted introgression of some genes involved in conflicts, either interspecific (parasites) or possibly cytonuclear. In the latter case, nuclear introgression could mitigate the potential negative effects of alien mtDNA on mitochondrial metabolism and male-specific traits.

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.

The telomeric sync model of speciation: species-wide telomere erosion triggers cycles of transposon-mediated genomic rearrangements, which underlie the saltatory appearance of nonadaptive characters

Charles Darwin knew that the fossil record is not overwhelmingly supportive of genetic and phenotypic gradualism; therefore, he developed the core of his theory on the basis of breeding experiments. Here, I present evidence for the existence of a cell biological mechanism that strongly points to the almost forgotten European concept of saltatory evolution of nonadaptive characters, which is in perfect agreement with the gaps in the fossil record. The standard model of chromosomal evolution has always been handicapped by a paradox, namely, how speciation can occur by spontaneous chromosomal rearrangements that are known to decrease the fertility of heterozygotes in a population. However, the hallmark of almost all closely related species is a differing chromosome complement and therefore chromosomal rearrangements seem to be crucial for speciation. Telomeres, the caps of eukaryotic chromosomes, erode in somatic tissues during life, but have been thought to remain stable in the germline of a species. Recently, a large human study spanning three healthy generations clearly found a cumulative telomere effect, which is indicative of transgenerational telomere erosion in the human species. The telomeric sync model of speciation presented here is based on telomere erosion between generations, which leads to identical fusions of chromosomes and triggers a transposon-mediated genomic repatterning in the germline of many individuals of a species. The phenotypic outcome of the telomere-triggered transposon activity is the saltatory appearance of nonadaptive characters simultaneously in many individuals. Transgenerational telomere erosion is therefore the material basis of aging at the species level.

Home-loving boreal hare mitochondria survived several invasions in Iberia: the relative roles of recurrent hybridisation and allele surfing

Genetic introgression from a resident species into an invading close relative can result from repeated hybridisation along the invasion front and/or allele surfing on the expansion wave. Cases where the phenomenon is massive and systematic, such as for hares (genus Lepus) in Iberia, would be best explained by recurrent hybridisation but this is difficult to prove because the donor populations are generally extinct. In the Pyrenean foothills, Lepus europaeus presumably replaced Lepus granatensis recently and the present species border is parallel to the direction of invasion, so that populations of L. granatensis in the contact zone represent proxies of existing variation before the invasion. Among three pairs of populations sampled across this border, we find less differentiation of mitochondrial DNA (mtDNA) across than along it, as predicted under recurrent hybridisation at the invasion front. Using autosomal microsatellite loci and X- and Y-linked diagnostic loci, we show that admixture across the border is quasi-absent, making it unlikely that lack of interspecific mtDNA differentiation results from ongoing gene flow. Furthermore, we find that the local species ranges are climatically contrasted, making it also unlikely that ongoing ecology-driven movement of the contact account for mtDNA introgression. The lack of mtDNA differentiation across the boundary is mostly due to sharing of mtDNA from a boreal species currently extinct in Iberia (Lepus timidus) whose mitochondria have thus remained in place since the last deglaciation despite successive invasions by two other species. Home-loving mitochondria thus witness past species distribution rather than ongoing exchanges across stabilised contact zones.

Low, complex and probably reticulated chromosome evolution of Sciuromorpha (Rodentia) and Lagomorpha

Lagomorpha (rabbits and pikas) and Sciuromorpha (squirrels) are grouped in the Glires superorder. Their chromosome diversification, since their separation from the eutherian mammalian common ancestor, was characterized by a low rate of chromosome rearrangements. Consequently, the structure of some chromosomes was either conserved or only slightly modified, making their comparison easy at the genus, family and even order level. Interspecific in situ hybridization (Zoo-FISH) largely corroborates classical cytogenetic data but provides much more reliability in comparisons, especially for distant species. We reconstructed common ancestral karyotypes for Glires, Lagomorpha, Sciuromorpha, and Sciuridae species, and then, determined the chromosome changes separating these ancestors from their common eutherian ancestor. We propose that reticulated evolution occurred during the diversification of Glires, which implies that several pericentric inversions and Robertsonian translocations were conserved in the heterozygous status for an extensive period. Finally, among Lagomorpha and Sciuromorpha, we focused on Leporidae and Sciuridae chromosome evolution. In the various attempts to establish dichotomic evolutionary schemes, it was necessary to admit that multiple homoplasies (convergent and reverse rearrangements) occurred in Sciuridae and in a lesser degree, in Leporidae. In Leporidae, additional rearrangements were sufficient to propose a resolved phylogeny. However, a resolved phylogeny was not possible for Sciuridae because most of the rearrangements occurred in terminal branches. We conclude that a reticulated evolution took place early during the evolution of both families and lasted longer in Sciuridae than in Leporidae. In Sciuridae, most chromosome rearrangements were pericentric inversions involving short fragments. Such rearrangements have only mild meiotic consequences, which may explain the long persistence of the heterozygous status characterizing reticulated evolution.

Reconstruction of karyotype evolution in core Glires. I. The genome homology revealed by comparative chromosome painting

Glires represent a eutherian clade consisting of rodents and lagomorphs (hares, rabbits, and pikas). Chromosome evolution of Glires is known to have variable rates in different groups: from slowly evolving lagomorphs and squirrels to extremely rapidly evolving muroids. Previous interordinal homology maps between slowly evolving Glires were based on comparison with humans. Here, we used sets of chromosome-specific probes from Tamias sibiricus (Sciuridae), Castor fiber (Castoridae) and humans to study karyotypes of six ground squirrels (genera Marmota and Spermophilus) and one tree squirrel (genus Sciurus), mountain hare (genus Lepus), and rabbit (genus Oryctolagus). These data supplemented with GTG banding comparisons allowed us to build comparative chromosome maps. Our data showed the absence of previously found squirrel associations HSA 1/8 and 2/17 in the Eurasian ground squirrels--sousliks and woodchucks, and disruptions of squirrel HSA 10/13 and HSA 8/4/8/12/22 syntenies in the four Spermophilus species studied here. We found that the karyotypes of Sciuridae and Leporidae are highly conserved and close to the Rodentia ancestral karyotype, while Castoridae chromosomes underwent many more changes. We suggest that Lagomorpha and Sciuridae (in contrast to all other rodent families) should be considered as core Glires lineages, characterized by cytogenetically conserved karyotypes which contain chromosomal elements inherent to karyotype of common Glires ancestor. Our data allowed us to further refine the putative ancestral karyotypes of Rodentia. We also describe here the putative ancestral karyotypes of Glires and lagomorphs.

Interspecific X-chromosome and mitochondrial DNA introgression in the Iberian hare: selection or allele surfing?

Introgression from a resident species into an invading one is predicted to occur through the demographic process of "allele surfing," and to particularly affect genomic regions transmitted by the lower migrating sex, such as mtDNA. This could explain that northern Iberian populations of Lepus granatensis harbor high frequencies of mtDNA from L. timidus, an arctic hare it replaced there after deglaciation. We report that variation of introgressed timidus-like mtDNA reflects several predicted effects of this process: increasing frequency and diversity in the direction of expansion, strong perpendicular phylogeographic structure and signs of postglacial demographic growth. However, demographic inferences for the granatensis and timidus-like mtDNA lineages suggest the latter may have outcompeted the former in northern Iberia. Autosomal introgression occurs at low frequencies and species-wide rather than only in the north. If this difference with mtDNA resulted from sex-biased migration, an intermediate pattern should prevail for the X-chromosome, but we report species-wide and high-frequency introgression of an X-fragment. Either selection favored this ubiquitous X-introgression, or more complex postglacial expansion patterns prevailed, with different consequences depending on the genomic and geographic region. This illustrates the difficulty of distinguishing demographic and selective effects and the need for genome and species-wide based demographic models.

Retroposon insertions provide insights into deep lagomorph evolution

The homogenous mammalian order Lagomorpha comprises about 80 species in two families, Ochotonidae (pikas) and Leporidae (rabbits and hares). However, the phylogenetic relationships among leporids are controversial. Molecular data, particularly from mitochondrial sequences, give highly homoplasious signals. To resolve the controversy between mitochondrial and nuclear data, we analyzed genomic orthologous retroposon insertion sites, a virtually homoplasy-free marker system. From a differential screen of rabbit genomic data for intronic retroposon insertions of CSINE elements, we polymerase chain reaction-amplified and sequenced 11 retroposons in eight representative lagomorphs. We found three retroposons shared among all lagomorphs but absent in outgroups, four confirmed the monophyly of leporids, and three significantly supported Pronolagus as the sister group to all other leporids. One retroposon supported the monophyly of Lepus. The position of Pronolagus outside of the remaining leporids supports the sequence-based signals of nuclear genes and clearly refutes the misleading signals of mitochondrial genes.

Phylogenomics of African guenons

The karyotypes of 28 specimens belonging to 26 species of Cercopithecinae have been compared with each other and with human karyotype by chromosome banding and, for some of them, by Zoo-FISH (human painting probes) techniques. The study includes the first description of the karyotypes of four species and a synonym of Cercopithecus nictitans. The chromosomal homologies obtained provide us with new data on a large number of rearrangements. This allows us to code chromosomal characters to draw Cercopithecini phylogenetic trees, which are compared to phylogenetic data based on DNA sequences. Our findings show that some of the superspecies proposed by Kingdon (1997 The Kingdon Field Guide to African Mammals, Academic Press.) and Groves (2001 Primates Taxonomy, Smithsonian Institution Press) do not form homogeneous groups and that the genus Cercopithecus is paraphyletic, in agreement with previous molecular analyses. The evolution of Cercopithecini karyotypes is mainly due to non-centromeric chromosome fissions and centromeric shifts or inversions. Non-Robertsonian translocations occurred in C. hamlyni and C. neglectus. The position of chromosomal rearrangements in the phylogenetic tree leads us to propose that the Cercopithecini evolution proceeded by either repeated fission events facilitated by peculiar genomic structures or successive reticulate phases, in which heterozygous populations for few rearranged chromosomes were present, allowing the spreading of chromosomal forms in various combinations, before the speciation process.

Isolation, cloning and characterization of two major satellite DNA families of rabbit (Oryctolagus cuniculus)

We report here the isolation, cloning and characterization of two abundant centromeric satellite sequences (Rsat I and Rsat II) what are not related to each other, and that of a divergent subfamily (Rsat IIE) of rabbit (Oryctolagus cuniculus). The Rsat I monomers had a 375 base pair (bp) average length, while repeat units Rsat II and Rsat IIE were approximately 585 bp long. Variable amounts of Rsat I were detected by FISH at the centromeric region of 11 chromosome pairs of the complement. Rsat II hybridized to the centromere of 12 different chromosomes, and two of these were labeled also with the Rsat IIE probe. Two-color in situ hybridizations with the satellite probes and rDNA revealed that the NOR chromosomes carried different satellites. Rsat I was abundant on chromosome 20 and 21, but it was undetectable on chromosomes 13 and 16. Large Rsat II arrays were found on chromosomes 16, 20 and 21, but reduced amount was detected on chromosome 13. The variant Rsat IIE was prominent on chromosome 16, but was absent from the other rDNA-bearing chromosomes. The rDNA signal on chromosome 21 was localized to the 21q(ter) region, what can be a useful cytological marker in comparative cytological studies. These data show that rabbit chromosomes form at least four distinct groups based on the satellite composition of their centromeres. The differences in the chromosomal distribution of satellite families will help easy FISH identification of individual chromosomes, as well as to unveil the evolutionary history of the Leporidae karyotype.

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.

Cytogenetics and cladistics

Chromosomal data have been underutilized in phylogenetic investigations despite the obvious potential that cytogenetic studies have to reveal both structural and functional homologies among taxa. In large part this is associated with difficulties in scoring conventional and molecular cytogenetic information for phylogenetic analysis. The manner in which chromosomal data have been used by most authors in the past was often conceptionally flawed in terms of the methods and principles underpinning modern cladistics. We present herein a review of the different methods employed, examine their relative strengths, and then outline a simple approach that considers the chromosomal change as the character, and its presence or absence the character state. We test this using one simulated and several empirical data sets. Features that are unique to cytogenetic investigations, including B-chromosomes, heterochromatic additions/deletions, and the location and number of nucleolar organizer regions (NORs), as well as the weighting of chromosomal characters, are critically discussed with regard to their suitability for phylogenetic reconstruction. We conclude that each of these classes of data have inherent problems that limit their usefulness in phylogenetic analyses and in most of these instances, inclusion should be subject to rigorous appraisal that addresses the criterion of unequivocal homology.