Insights into the evolution of the new variant rabbit haemorrhagic disease virus (GI.2) and the identification of novel recombinant strains.

Rabbit haemorrhagic disease (RHD) is a viral disease that affects the European rabbit. RHD was detected in 1984 in China and rapidly disseminated worldwide causing a severe decline in wild rabbit populations. The aetiological agent, rabbit haemorrhagic disease virus (RHDV), is an RNA virus of the family Caliciviridae, genus Lagovirus. Pathogenic (G1-G6 or variants GI.1a-GI.1d) and non-pathogenic strains (GI.4) have been characterized. In 2010, a new variant of RHDV, RHDV2/RHDVb/GI.2, was detected in France. GI.2 arrived to the Iberian Peninsula in 2011, and several recombination events were reported. Here, we sequenced full genomes of 19 samples collected in Portugal between 2014 and 2016. New GI.2 recombinant strains were detected, including triple recombinants. These recombinants possess a non-structural protein p16 related to a non-pathogenic strain. Evolutionary analyses were conducted on GI.2 VP60 sequences. Estimated time to the most recent common ancestor (tMRCA) suggests an emergence of GI.2 in July 2008, not distant from its first detection in 2010. This is the first study on GI.2 evolution and highlights the need of continued monitoring and characterization of complete genome sequences when studying lagoviruses' evolution.

Molecular and ecological signs of mitochondrial adaptation: consequences for introgression?

The evolution of the mitochondrial genome and its potential adaptive impact still generates vital debates. Even if mitochondria have a crucial functional role, as they are the main cellular energy suppliers, mitochondrial DNA (mtDNA) introgression is common in nature, introducing variation in populations upon which selection may act. Here we evaluated whether the evolution of mtDNA in a rodent species affected by mtDNA introgression is explained by neutral expectations alone. Variation in one mitochondrial and six nuclear markers in Myodes glareolus voles was examined, including populations that show mtDNA introgression from its close relative, Myodes rutilus. In addition, we modelled protein structures of the mtDNA marker (cytochrome b) and estimated the environmental envelopes of mitotypes. We found that massive mtDNA introgression occurred without any trace of introgression in the analysed nuclear genes. The results show that the native glareolus mtDNA evolved under past positive selection, suggesting that mtDNA in this system has selective relevance. The environmental models indicate that the rutilus mitotype inhabits colder and drier habitats than the glareolus one that can result from local adaptation or from the geographic context of introgression. Finally, homology models of the cytochrome b protein revealed a substitution in rutilus mtDNA in the vicinity of the catalytic fraction, suggesting that differences between mitotypes may result in functional changes. These results suggest that the evolution of mtDNA in Myodes may have functional, ecological and adaptive significance. This work opens perspective onto future experimental tests of the role of natural selection in mtDNA introgression in this system.

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.

Cryptic speciation in the field vole: a multilocus approach confirms three highly divergent lineages in Eurasia.

Species are generally described from morphological features, but there is growing recognition of sister forms that show substantial genetic differentiation without obvious morphological variation and may therefore be considered 'cryptic species'. Here, we investigate the field vole (Microtus agrestis), a Eurasian mammal with little apparent morphological differentiation but which, on the basis of previous sex-linked nuclear and mitochondrial DNA (mtDNA) analyses, is subdivided into a Northern and a Southern lineage, sufficiently divergent that they may represent two cryptic species. These earlier studies also provided limited evidence for two major mtDNA lineages within Iberia. In our present study, we extend these findings through a multilocus approach. We sampled 163 individuals from 46 localities, mainly in Iberia, and sequenced seven loci, maternally, paternally and biparentally inherited. Our results show that the mtDNA lineage identified in Portugal is indeed a distinct third lineage on the basis of other markers as well. In fact, multilocus coalescent-based methods clearly support three separate evolutionary units that may represent cryptic species: Northern, Southern and Portuguese. Divergence among these units was inferred to have occurred during the last glacial period; the Portuguese lineage split occurred first (estimated at c. 70 000 bp), and the Northern and Southern lineages separated at around the last glacial maximum (estimated at c. 18 500 bp). Such recent formation of evolutionary units that might be considered species has repercussions in terms of understanding evolutionary processes and the diversity of small mammals in a European context.

Reference-free transcriptome assembly in non-model animals from next-generation sequencing data.

Next-generation sequencing (NGS) technologies offer the opportunity for population genomic study of non-model organisms sampled in the wild. The transcriptome is a convenient and popular target for such purposes. However, designing genetic markers from NGS transcriptome data requires assembling gene-coding sequences out of short reads. This is a complex task owing to gene duplications, genetic polymorphism, alternative splicing and transcription noise. Typical assembling programmes return thousands of predicted contigs, whose connection to the species true gene content is unclear, and from which SNP definition is uneasy. Here, the transcriptomes of five diverse non-model animal species (hare, turtle, ant, oyster and tunicate) were assembled from newly generated 454 and Illumina sequence reads. In two species for which a reference genome is available, a new procedure was introduced to annotate each predicted contig as either a full-length cDNA, fragment, chimera, allele, paralogue, genomic sequence or other, based on the number of, and overlap between, blast hits to the appropriate reference. Analyses showed that (i) the highest quality assemblies are obtained when 454 and Illumina data are combined, (ii) typical de novo assemblies include a majority of irrelevant cDNA predictions and (iii) assemblies can be appropriately cleaned by filtering contigs based on length and coverage. We conclude that robust, reference-free assembly of thousands of genes from transcriptomic NGS data is possible, opening promising perspectives for transcriptome-based population genomics in animals. A Galaxy pipeline implementing our best-performing assembling strategy is provided.

Recurrent introgression of mitochondrial DNA among hares (Lepus spp.) revealed by species-tree inference and coalescent simulations.

Understanding recent speciation history requires merging phylogenetic and population genetics approaches, taking into account the persistence of ancestral polymorphism and possible introgression. The emergence of a clear phylogeny of hares (genus Lepus) has been hampered by poor genomic sampling and possible occurrence of mitochondrial DNA (mtDNA) introgression from the arctic/boreal Lepus timidus into several European temperate and possibly American boreal species. However, no formal test of introgression, taking also incomplete lineage sorting into account, has been done. Here, to clarify the yet poorly resolved species phylogeny of hares and test hypotheses of mtDNA introgression, we sequenced 14 nuclear DNA and 2 mtDNA fragments (8205 and 1113 bp, respectively) in 50 specimens from 11 hare species from Eurasia, North America, and Africa. By applying an isolation-with-migration model to the nuclear data on subsets of species, we find evidence for very limited gene flow from L. timidus into most temperate European species, and not into the American boreal ones. Using a multilocus coalescent-based method, we infer the species phylogeny, which we find highly incongruent with mtDNA phylogeny using parametric bootstrap. Simulations of mtDNA evolution under the speciation history inferred from nuclear genes did not support the hypothesis of mtDNA introgression from L. timidus into the American L. townsendii but did suggest introgression from L. timidus into 4 temperate European species. One such event likely resulted in the complete replacement of the aboriginal mtDNA of L. castroviejoi and of its sister species L. corsicanus. It is remarkable that mtDNA introgression in hares is frequent, extensive, and always from the same donor arctic species. We discuss possible explanations for the phenomenon in relation to the dynamics of range expansions and species replacements during the climatic oscillations of the Pleistocene.

The genomic legacy from the extinct Lepus timidus to the three hare species of Iberia: contrast between mtDNA, sex chromosomes and autosomes.

Extensive interspecific genetic introgression is often reported, and appraising its genomic impact can serve to determine whether it results from selection on specific loci or from demographic processes affecting the whole genome. The three species of hares present in the Iberian Peninsula harbour high frequencies of mitochondrial DNA (mtDNA) from Lepus timidus, an arctic/boreal species now extinct in the region. This could result from the invasive replacement of L. timidus by the temperate species during deglaciation but should then have left traces in the nuclear genome. We typed single nucleotide polymorphisms (SNPs) discovered by sequencing 10 autosomal loci, two X-linked and one Y-linked in species-wide samples of the four taxa. Based on lineage-diagnostic SNPs, we detected no trace of L. timidus sex chromosomes in Iberia. From the frequencies of inferred haplotypes, autosomal introgression into L. granatensis appeared mostly sporadic but always widespread instead of restricted to the north as mtDNA. Autosomal introgression into Iberian L. europaeus, inhabiting the Pyrenean foothills, was hardly detectable, despite quasi-fixation of L. timidus mtDNA. L. castroviejoi, endemic to the Cantabrian Mountains and fixed for L. timidus mtDNA, showed little traces of autosomal introgression. The absence of sex-chromosome introgression presumably resulted from X-linked hybrid male unfitness. The contrasting patterns between the autosomes and mtDNA could reflect general gender asymmetric processes such as frequency-dependent female assortative mating, lower mtDNA migration and higher male dispersal, but adaptive mtDNA introgression cannot be dismissed. Additionally, we document reciprocal introgression between L. europaeus and both L. granatensis in Iberia and L. timidus outside Iberia.

The rise and fall of the mountain hare (Lepus timidus) during Pleistocene glaciations: expansion and retreat with hybridization in the Iberian Peninsula.

The climatic fluctuations during glaciations have affected differently arctic and temperate species. In the northern hemisphere, cooling periods induced the expansion of many arctic species to the south, while temperate species were forced to retract in southern refugia. Consequently, in some areas the alternation of these species set the conditions for competition and eventually hybridization. Hares in the Iberian Peninsula appear to illustrate this phenomenon. Populations of Iberian hare (Lepus granatensis), brown hare (Lepus europaeus) and broom hare (Lepus castroviejoi) in Northern Iberia harbour mitochondrial haplotypes from the mountain hare (Lepus timidus), a mainly boreal and arctic species presently absent from the peninsula. To understand the history of this past introgression we analysed sequence variation and geographical distribution of mitochondrial control region and cytochrome b haplotypes of L. timidus origin found in 378 specimens of these four species. Among 124 L. timidus from the Northern Palaearctic and the Alps we found substantial nucleotide diversity (2.3%) but little differentiation between populations. Based on the mismatch distribution of the L. timidus sequences, this could result from an expansion at a time of temperature decrease favourable to this arctic species. The nucleotide diversity of L. timidus mtDNA found in Iberian L. granatensis, L. europaeus and L. castroviejoi (183, 70 and 1 specimens, respectively) was of the same order as that in L. timidus over its range (1.9%), suggesting repeated introgression of multiple lineages. The structure of the coalescent of L. granatensis sequences indicates that hybridization with L. timidus was followed by expansion of the introgressed haplotypes, as expected during a replacement with competition, and occurred when temperatures started to rise, favouring the temperate species. Whether a similar scenario explains the introgression into Iberian L. europaeus remains unclear but it is possible that it hybridized with already introgressed L. granatensis.

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.