Papers by José Melo-Ferreira
Heredity, 2014
The evolution of the mitochondrial genome and its potential adaptive impact still generates vital... more 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.
Systematic Biology, 2012
Understanding recent speciation history requires merging phylogenetic and population genetics app... more 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.
Molecular Ecology, 2014
With climate warming, the ranges of many boreal species are expected to shift northward and to fr... more With climate warming, the ranges of many boreal species are expected to shift northward and to fragment in southern peripheral ranges. To understand the conservation implications of losing southern populations, we examined range-wide genetic diversity of the snowshoe hare (Lepus americanus), an important prey species that drives boreal ecosystem dynamics. We analysed microsatellite (8 loci) and mitochondrial DNA sequence (cytochrome b and control region) variation in almost 1000 snowshoe hares. A hierarchical structure analysis of the microsatellite data suggests initial subdivision in two groups, Boreal and southwestern. The southwestern group further splits into Greater Pacific Northwest and U.S. Rockies. The genealogical information retrieved from mtDNA is congruent with the three highly differentiated and divergent groups of snowshoe hares. These groups can correspond with evolutionarily significant units that might have evolved in separate refugia south and east of the Pleistocene ice sheets. Genetic diversity was highest at mid-latitudes of the species' range, and genetic uniqueness was greatest in southern populations, consistent with substructuring inferred from both mtDNA and microsatellite analyses at finer levels of analysis. Surprisingly, snowshoe hares in the Greater Pacific Northwest mtDNA lineage were more closely related to black-tailed jackrabbits (Lepus californicus) than to other snowshoe hares, which may result from secondary introgression or shared ancestral polymorphism. Given the genetic distinctiveness of southern populations and minimal gene flow with their northern neighbours, fragmentation and loss of southern boreal habitats could mean loss of many unique alleles and reduced evolutionary potential.
Biological Journal of the Linnean Society, 2014
Immunogenetics, 2014
Antigen recognition by immunoglobulins depends upon initial rearrangements of heavy chain V, D, a... more Antigen recognition by immunoglobulins depends upon initial rearrangements of heavy chain V, D, and J genes. In leporids, a unique system exists for the VH genes usage that exhibit highly divergent lineages: the VHa allotypes, the Lepus sL lineage and the VHn genes. For the European rabbit (Oryctolagus cuniculus), four VHa lineages have been described, the a1, a2, a3 and a4. For hares (Lepus sp.), one VHa lineage was described, the a2L, as well as a more ancient sL lineage. Both genera use the VHn genes in a low frequency of their VDJ rearrangements. To address the hypothesis that the VH specificities could be associated with different environments, we sequenced VDJ genes from a third leporid genus, Sylvilagus. We found a fifth and equally divergent VHa lineage, the a5, and an ancient lineage, the sS, related to the hares' sL, but failed to obtain VHn genes. These results show that the studied leporids employ different VH lineages in the generation of the antibody repertoire, suggesting that the leporid VH genes are subject to strong selective pressure likely imposed by specific pathogens.
Molecular Phylogenetics and Evolution, 2015
Molecular ecology, 2014
Hybridization drives the evolutionary trajectory of many species or local populations, and assess... more Hybridization drives the evolutionary trajectory of many species or local populations, and assessing the geographic extent and genetic impact of interspecific gene flow may provide invaluable clues to understand population divergence or the adaptive relevance of admixture. In North America, hares (Lepus spp.) are key species for ecosystem dynamics and their evolutionary history may have been affected by hybridization.
PLoS Genetics, 2013
In animals, the population genomic literature is dominated by two taxa, namely mammals and drosop... more In animals, the population genomic literature is dominated by two taxa, namely mammals and drosophilids, in which fully sequenced, well-annotated genomes have been available for years. Data from other metazoan phyla are scarce, probably because the vast majority of living species still lack a closely related reference genome. Here we achieve de novo, referencefree population genomic analysis from wild samples in five non-model animal species, based on next-generation sequencing transcriptome data. We introduce a pipe-line for cDNA assembly, read mapping, SNP/genotype calling, and data cleaning, with specific focus on the issue of hidden paralogy detection. In two species for which a reference genome is available, similar results were obtained whether the reference was used or not, demonstrating the robustness of our de novo inferences. The population genomic profile of a hare, a turtle, an oyster, a tunicate, and a termite were found to be intermediate between those of human and Drosophila, indicating that the discordant genomic diversity patterns that have been reported between these two species do not reflect a generalized vertebrate versus invertebrate gap. The genomic average diversity was generally higher in invertebrates than in vertebrates (with the notable exception of termite), in agreement with the notion that population size tends to be larger in the former than in the latter. The non-synonymous to synonymous ratio, however, did not differ significantly between vertebrates and invertebrates, even though it was negatively correlated with genetic diversity within each of the two groups. This study opens promising perspective regarding genome-wide population analyses of non-model organisms and the influence of population size on nonsynonymous versus synonymous diversity.
Virus Research, 2015
Endogenization of mouse mammary tumor virus (MMTV)-like elements in genomes of pikas (Ochotona sp... more Endogenization of mouse mammary tumor virus (MMTV)-like elements in genomes of pikas (Ochotona sp.).Virus Research http://dx.
Molecular Ecology Resources, 2012
Next-generation sequencing (NGS) technologies offer the opportunity for population genomic study ... more 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 fulllength 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.
Molecular Ecology, 2009
Extensive interspecific genetic introgression is often reported, and appraising its genomic impac... more 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 quasifixation 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.
Molecular Phylogenetics and Evolution, 2008
Global Change Biology, 2012
Molecular Ecology, 2012
Species are generally described from morphological features, but there is growing recognition of ... more 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.
PLoS ONE, 2012
The application of species distribution models (SDMs) in ecology and conservation biology is incr... more The application of species distribution models (SDMs) in ecology and conservation biology is increasing and assuming an important role, mainly because they can be used to hindcast past and predict current and future species distributions. However, the accuracy of SDMs depends on the quality of the data and on appropriate theoretical frameworks. In this study, comprehensive data on the current distribution of the Iberian hare (Lepus granatensis) were used to i) determine the species' ecogeographical constraints, ii) hindcast a climatic model for the last glacial maximum (LGM), relating it to inferences derived from molecular studies, and iii) calibrate a model to assess the species future distribution trends (up to 2080). Our results showed that the climatic factor (in its pure effect and when it is combined with the land-cover factor) is the most important descriptor of the current distribution of the Iberian hare. In addition, the model's output was a reliable index of the local probability of species occurrence, which is a valuable tool to guide species management decisions and conservation planning. Climatic potential obtained for the LGM was combined with molecular data and the results suggest that several glacial refugia may have existed for the species within the major Iberian refugium. Finally, a high probability of occurrence of the Iberian hare in the current species range and a northward expansion were predicted for future. Given its current environmental envelope and evolutionary history, we discuss the macroecology of the Iberian hare and its sensitivity to climate change.
Five species of genus Lepus occur naturally in Europe: L. europaeus, L. timidus, L. granatensis, ... more Five species of genus Lepus occur naturally in Europe: L. europaeus, L. timidus, L. granatensis, L. corsicanus, and L. castroviejoi. Of these, the latter two have restricted ranges, L. castroviejoi in the Iberian Peninsula and L. corsicanus in central and southern Italy.
European Journal of Wildlife Research, 2011
Abstract: The Italian hare, Lepus corsicanus, was first described in Corsica more than 100 years ... more Abstract: The Italian hare, Lepus corsicanus, was first described in Corsica more than 100 years ago, but the knowledge on the status of the species in this island remains scarce. Moreover, frequent introductions of thousands of individuals from other hare species, namely L. europaeus and L. granatensis, into Corsica are known to have occurred and an updated assessment of the prevalence of L. corsicanus in Corsica is therefore of the utmost importance. Here, to estimate the relative prevalence of the hare species present in Corsica we conducted a molecular analysis on 67 samples collected by hunters between 2002 and 2007 in 36 Corsican communes. Sequencing of portions of the nuclear gene transferrin and of the control region of the mitochondrial DNA allowed classifying most of the collected samples as belonging to L. corsicanus (70.1%). Of the sampled Corsican communes, 86.1% contained this species, while only in 11.1% L. europaeus was present. Three of the analysed specimens showed an inconsistent molecular assignment between markers suggesting a hybrid origin: L. corsicanus x L. europaeus, L. corsicanus x L. granatensis and L. europaeus x L. granatensis. The first two cases of hybridization had never been described in nature, even in studies focusing on hares from Italy where L. corsicanus and L. europaeus are often sympatric. These results stress the real risk of corrosion of the native gene pool of L. corsicanus via hybridization with introduced species. We highlight the need of urgently rethinking the management plan of hare populations in Corsica.
BMC Evolutionary Biology, 2011
Background: Introgression of mitochondrial DNA (mtDNA) is among the most frequently described cas... more Background: Introgression of mitochondrial DNA (mtDNA) is among the most frequently described cases of reticulate evolution. The tendency of mtDNA to cross interspecific barriers is somewhat counter-intuitive considering the key function of enzymes that it encodes in the oxidative-phosphorylation process, which could give rise to hybrid dysfunction. How mtDNA reticulation affects the evolution of metabolic functions is, however, uncertain. Here we investigated how morpho-physiological traits vary in natural populations of a common rodent (the bank vole, Myodes glareolus) and whether this variation could be associated with mtDNA introgression. First, we confirmed that M. glareolus harbour mtDNA introgressed from M. rutilus by analyzing mtDNA (cytochrome b, 954 bp) and nuclear DNA (four markers; 2333 bp in total) sequence variation and reconstructing loci phylogenies among six natural populations in Finland. We then studied geographic variation in body size and basal metabolic rate (BMR) among the populations of M. glareolus and tested its relationship with mtDNA type.
Molecular Biology and Evolution, 2012
The nearly neutral theory of molecular evolution predicts that the efficacy of both positive and ... more The nearly neutral theory of molecular evolution predicts that the efficacy of both positive and purifying selection is a function of the long-term effective population size (N e ) of a species. Under this theory, the efficacy of natural selection should increase with N e . Here, we tested this simple prediction by surveying ;1.5 to 1.8 Mb of protein coding sequence in the two subspecies of the European rabbit (Oryctolagus cuniculus algirus and O. c. cuniculus), a mammal species characterized by high levels of nucleotide diversity and N e estimates for each subspecies on the order of 1 Â 10 6 . When the segregation of slightly deleterious mutations and demographic effects were taken into account, we inferred that .60% of amino acid substitutions on the autosomes were driven to fixation by positive selection. Moreover, we inferred that a small fraction of new amino acid mutations (,4%) are effectively neutral (defined as 0 , N e s , 1) and that this fraction was negatively correlated with a gene's expression level. Consistent with models of recurrent adaptive evolution, we detected a negative correlation between levels of synonymous site polymorphism and the rate of protein evolution, although the correlation was weak and nonsignificant. No systematic X chromosome-autosome difference was found in the efficacy of selection. For example, the proportion of adaptive substitutions was significantly higher on the X chromosome compared with the autosomes in O. c. algirus but not in O. c. cuniculus. Our findings support widespread positive and purifying selection in rabbits and add to a growing list of examples suggesting that differences in N e among taxa play a substantial role in determining rates and patterns of protein evolution.
Uploads
Papers by José Melo-Ferreira