Evaluating the use of non-invasive hair sampling and ddRAD to characterize populations of endangered species: Application to a peripheral population of the European mink.

The application of next-generation sequencing (NGS) to non-invasive samples is one of the most promising methods in conservation genomics, but these types of samples present significant challenges for NGS. The European mink (Mustela lutreola) is critically endangered throughout its range. However, important aspects such as census size and inbreeding remain still unknown in many populations, so it is crucial to develop new methods to monitor this species. In this work, we placed hair tubes along riverbanks in a border area of the Iberian population, which allowed the genetic identification of 76 European mink hair samples. We then applied a reduced representation genomic sequencing (ddRAD) technique to a subset of these samples to test whether we could extract sufficient genomic information from them. We show that several problems with the DNA, including contamination, fragmentation, oxidation, and possibly sample mixing, affected the samples. Using various bioinformatic techniques to reduce these problems, we were able to unambiguously genotype 19 hair samples belonging to six individuals. This small number of individuals showed that the demographic status of the species in this peripheral population is worse than expected. The data obtained also allowed us to perform preliminary analyses of relatedness and inbreeding. Although further improvements in sampling and analysis are needed, the application of the ddRAD technique to non-invasively obtained hairs represents a significant advance in the genomic study of endangered species.

Methodological challenges in the genomic analysis of an endangered mammal population with low genetic diversity.

Recently, populations of various species with very low genetic diversity have been discovered. Some of these persist in the long term, but others could face extinction due to accelerated loss of fitness. In this work, we characterize 45 individuals of one of these populations, belonging to the Iberian desman (Galemys pyrenaicus). For this, we used the ddRADseq technique, which generated 1421 SNPs. The heterozygosity values of the analyzed individuals were among the lowest recorded for mammals, ranging from 26 to 91 SNPs/Mb. Furthermore, the individuals from one of the localities, highly isolated due to strong barriers, presented extremely high inbreeding coefficients, with values above 0.7. Under this scenario of low genetic diversity and elevated inbreeding levels, some individuals appeared to be almost genetically identical. We used different methods and simulations to determine if genetic identification and parentage analysis were possible in this population. Only one of the methods, which does not assume population homogeneity, was able to identify all individuals correctly. Therefore, genetically impoverished populations pose a great methodological challenge for their genetic study. However, these populations are of primary scientific and conservation interest, so it is essential to characterize them genetically and improve genomic methodologies for their research.

Divergence time estimation using ddRAD data and an isolation-with-migration model applied to water vole populations of Arvicola.

Molecular dating methods of population splits are crucial in evolutionary biology, but they present important difficulties due to the complexity of the genealogical relationships of genes and past migrations between populations. Using the double digest restriction-site associated DNA (ddRAD) technique and an isolation-with-migration (IM) model, we studied the evolutionary history of water vole populations of the genus Arvicola, a group of complex evolution with fossorial and semi-aquatic ecotypes. To do this, we first estimated mutation rates of ddRAD loci using a phylogenetic approach. An IM model was then used to estimate split times and other relevant demographic parameters. A set of 300 ddRAD loci that included 85 calibrated loci resulted in good mixing and model convergence. The results showed that the two populations of A. scherman present in the Iberian Peninsula split 34 thousand years ago, during the last glaciation. In addition, the much greater divergence from its sister species, A. amphibius, may help to clarify the controversial taxonomy of the genus. We conclude that this approach, based on ddRAD data and an IM model, is highly useful for analyzing the origin of populations and species.

The genome of the Pyrenean desman and the effects of bottlenecks and inbreeding on the genomic landscape of an endangered species.

The Pyrenean desman (Galemys pyrenaicus) is a small semiaquatic mammal endemic to the Iberian Peninsula. Despite its limited range, this species presents a strong genetic structure due to past isolation in glacial refugia and subsequent bottlenecks. Additionally, some populations are highly fragmented today as a consequence of river barriers, causing substantial levels of inbreeding. These features make the Pyrenean desman a unique model in which to study the genomic footprints of differentiation, bottlenecks and extreme isolation in an endangered species. To understand these processes, the complete genome of the Pyrenean desman was sequenced and assembled using a Bloom filter-based approach. An analysis of the 1.83 Gb reference genome and the sequencing of five additional individuals from different evolutionary units allowed us to detect its main genomic characteristics. The population differentiation of the species was reflected in highly distinctive demographic trajectories. In addition, a severe population bottleneck during the postglacial recolonization of the eastern Pyrenees created one of the lowest genomic heterozygosity values recorded in a mammal. Moreover, isolation and inbreeding gave rise to a high proportion of runs of homozygosity (ROH). Despite these extremely low levels of genetic diversity, two key multigene families from an eco-evolutionary perspective, the major histocompatibility complex and olfactory receptor genes, showed heterozygosity excess in the majority of individuals, revealing that functional diversity can be maintained up to a certain extent. Furthermore, these two classes of genes were significantly less abundant than expected within ROH. In conclusion, the genomic landscape of each analysed Pyrenean desman turned out to be strikingly distinctive and was a clear reflection of its recent ancestry and current conservation conditions. These results may help characterize the genomic health of each individual, and can be crucial for the conservation and management of the species.

Publisher Correction: Size increase without genetic divergence in the Eurasian water shrew Neomys fodiens.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Size increase without genetic divergence in the Eurasian water shrew Neomys fodiens.

When a population shows a marked morphological change, it is important to know whether that population is genetically distinct; if it is not, the novel trait could correspond to an adaptation that might be of great ecological interest. Here, we studied a subspecies of water shrew, Neomys fodiens niethammeri, which is found in a narrow strip of the northern Iberian Peninsula. This subspecies presents an abrupt increase in skull size when compared to the rest of the Eurasian population, which has led to the suggestion that it is actually a different species. Skulls obtained from owl pellets collected over the last 50 years allowed us to perform a morphometric analysis in addition to an extensive multilocus analysis based on short intron fragments successfully amplified from these degraded samples. Interestingly, no genetic divergence was detected using either mitochondrial or nuclear data. Additionally, an allele frequency analysis revealed no significant genetic differentiation. The absence of genetic divergence and differentiation revealed here indicate that the large form of N. fodiens does not correspond to a different species and instead represents an extreme case of size increase, of possible adaptive value, which deserves further investigation.

Quantitative analysis of connectivity in populations of a semi-aquatic mammal using kinship categories and network assortativity.

Analysing the impact of anthropogenic and natural river barriers on the dispersal of aquatic and semi-aquatic species may be critical for their conservation. Knowledge of kinship relationships between individuals and reconstructions of pedigrees obtained using genomic data can be extremely useful, not only for studying the social organization of animals, but also inferring contemporary dispersal and quantifying the effect of specific barriers on current connectivity. In this study, we used kinship data to analyse connectivity patterns in a small semi-aquatic mammal, the Pyrenean desman (Galemys pyrenaicus), in an area comprising two river systems with close headwaters and dams of various heights and types. Using a large SNP dataset from 70 specimens, we obtained kinship categories and reconstructed pedigrees. To quantify the barrier effect of specific obstacles, we built kinship networks and devised a method based on the assortativity coefficient, which measures the proportion between observed and expected kinship relationships across a barrier. The estimation of this parameter enabled us to infer that the most important barrier in the area was the watershed divide between the rivers, followed by a dam on one of the rivers. Other barriers did not significantly reduce the expected number of kinship relationships across them. This strategy and the information obtained with it may be crucial in determining the most important connectivity problems in an area and help develop conservation plans aimed at improving genetic exchange between populations of threatened species.

Evolutionary history of the endemic water shrew Neomys anomalus: Recurrent phylogeographic patterns in semi-aquatic mammals of the Iberian Peninsula.

The Cabrera's water shrew (Neomys anomalus) is a small semi-aquatic mammal whose taxonomic status was recently elevated from subspecies to species; as a consequence of this change, this species is now endemic to the Iberian Peninsula. In this study, we looked at its evolutionary history by combining phylogeography, the spatial distribution of genetic diversity, and species distribution modeling. To perform these analyses, we used noninvasive samples collected across the species distribution range and sequenced partial mitochondrial cytochrome b and D-loop genes. Maximum-likelihood and Bayesian phylogenetic trees derived from these sequences indicated that N. anomalus is divided into two main phylogroups that correlate strongly with geography, with two contact zones between the groups that showed limited spatial mixing between them. River basins were responsible for only a small percentage of the structure of the genetic diversity of this species despite its riparian habitat. The nucleotide diversity variation map showed the highest genetic diversity to be in the north of the Iberian Peninsula. Finally, species distribution modeling allowed the inference of an optimal area during the Last Interglacial in the north of the Iberian Peninsula, and multiple glacial refugia during the Last Glacial Maximum. The phylogeographic pattern of N. anomalus is strikingly similar to that of another semi-aquatic Iberian mammal, the Pyrenean desman (Galemys pyrenaicus), revealing how Pleistocene glaciations could have had equivalent effects on species of similar ecology and distribution. This phylogeographic structure is consistent with N. anomalus having been isolated for long periods in multiple glacial refugia within the Iberian Peninsula, in agreement with the "refugia-within-refugia" hypothesis, and further supporting its status as a distinct species.

A pipeline for metabarcoding and diet analysis from fecal samples developed for a small semi-aquatic mammal.

Metabarcoding allows the genetic analysis of pooled samples of various sources. It is becoming popular in the study of animal diet, especially because it allows the analysis of the composition of feces without the need of handling animals. In this work, we studied the diet of the Pyrenean desman (Galemys pyrenaicus), a small semi-aquatic mammal endemic to the Iberian Peninsula and the Pyrenees, by sequencing COI minibarcodes from feces using next-generation sequencing techniques. For the identification of assembled sequences, we employed a tree-based identification method that used a reference tree of sequences of freshwater organisms. The comparison of freshly collected fecal samples and older samples showed that fresh samples produced significantly more sequencing reads. They also rendered more operational taxonomical units (OTUs), but not significantly. Our analyses of 41 samples identified 224 OTUs corresponding to species of the reference tree. Ephemeroptera, Diptera excl. Chironomidae, and Chironomidae were the most highly represented groups in terms of reads as well as samples. Other groups of freshwater organisms detected were Plecoptera, Trichoptera, Neuropteroida, Coleoptera, Crustacea, and Annelida. Our results are largely in line with previous morphological and genetic studies on the diet of the Pyrenean desman, but allowed the identification of a higher diversity of OTUs in each sample. Additionally, the bioinformatic pipeline we developed for deep sequencing of fecal samples will enable the quantitative analysis of the diet of this and other species, which can be highly useful to determine their ecological requirements.

Postglacial dispersal patterns and mitochondrial genetic structure of the Pyrenean desman (Galemys pyrenaicus) in the northwestern region of the Iberian Peninsula.

The genetic structure of small semiaquatic animals may be influenced by dispersal across both rivers and land. The relative importance of these two modes of dispersal may vary across different species and with ecological conditions and evolutionary periods. The Pyrenean desman (Galemys pyrenaicus) is an endemic mammal of the Iberian Peninsula with a strong phylogeographic structure and semiaquatic habits, thus making it an ideal model to study the effects of river and overland dispersal on its genetic structure. Thanks to different types of noninvasive samples, we obtained an extensive sampling of the Pyrenean desman from the northwestern region of the Iberian Peninsula and sequenced two mitochondrial DNA fragments. We then analyzed, using an isolation-by-distance approach, the correlation between phylogenetic distances and geographical distances measured along both river networks and land to infer the relative importance of river and overland dispersal. We found that the correlations in the whole area and in a large basin were consistent with an effect of overland dispersal, which may be due to the postglacial colonization of new territories using terrestrial corridors and, possibly, a more extensive fluvial network that may have been present during the Holocene. However, in a small basin, likely to be less influenced by the impact of ancient postglacial dispersal, the correlations suggested significant overall effects of both overland and river dispersal, as expected for a semiaquatic mammal. Therefore, different scales and geographical regions reflect different aspects of the evolutionary history and ecology of this semiaquatic species using this isolation-by-distance method. The results we obtained may have crucial implications for the conservation of the Pyrenean desman because they reinforce the importance of interbasin dispersal for this species in the studied area and the need to protect the whole riverine ecosystem, including rivers, upland streams and terrestrial corridors between basins.

Using relatedness networks to infer contemporary dispersal: Application to the endangered mammal Galemys pyrenaicus.

Information about the degree of contemporary dispersal is important when trying to understand how populations interchange individuals and identify the specific barriers that prevent these movements. In the case of endangered species, this can represent crucial information when designing appropriate conservation strategies. Here we analyse relatedness between individuals from different localities and use these data to infer whether dispersal occurred in recent generations. We applied this approach to the Pyrenean desman (Galemys pyrenaicus), a semiaquatic and endangered species endemic to the Iberian Peninsula. We studied this species in four primary rivers in the Iberian Range, where two ancient mitochondrial lineages are separated by a strict contact zone, suggesting the existence of complex dispersal patterns. Using next-generation sequencing, we obtained 912 SNPs from each specimen and estimated relatedness values between them. While relatedness networks were dense within each river, we found surprisingly few relationships between individuals from different rivers despite their close proximity in some cases, indicating much lower dispersal between rivers compared to dispersal within a single river. In agreement with this result, the degree of inbreeding was exceedingly high in most individuals. These data show that relatedness information can be crucial to understand the contemporary dispersal patterns and conservation status of specific populations of endangered species.

Multilocus species trees and species delimitation in a temporal context: application to the water shrews of the genus Neomys.

BACKGROUND: Multilocus data are becoming increasingly important in determining the phylogeny of closely related species and delimiting species. In species complexes where unequivocal fossil calibrations are not available, rigorous dating of the coalescence-based species trees requires accurate mutation rates of the loci under study but, generally, these rates are unknown. Here, we obtained lineage-specific mutation rates of these loci from a higher-level phylogeny with a reliable fossil record and investigated how different choices of mutation rates and species tree models affected the split time estimates. We implemented this strategy with a genus of water shrews, Neomys, whose taxonomy has been contentious over the last century. RESULTS: We sequenced 13 introns and cytochrome b from specimens of the three species currently recognized in this genus including two subspecies of N. anomalus that were originally described as species. A Bayesian multilocus species delimitation method and estimation of gene flow supported that these subspecies are distinct evolutionary lineages that should be treated as distinct species: N. anomalus (sensu stricto), limited to part of the Iberian Peninsula, and N. milleri, with a larger Eurasian range. We then estimated mutation rates from a Bayesian relaxed clock analysis of the mammalian orthologues with several fossil calibrations. Next, using the estimated Neomys-specific rates for each locus in an isolation-with-migration model, the split time for these sister taxa was dated at 0.40 Myr ago (with a 95 % confidence interval of 0.26 - 0.86 Myr), likely coinciding with one of the major glaciations of the Middle Pleistocene. We also showed that the extrapolation of non-specific rates or the use of simpler models would lead to very different split time estimates. CONCLUSIONS: We showed that the estimation of rigorous lineage-specific mutation rates for each locus allows the inference of robust split times in a species tree framework. These times, in turn, afford a better understanding of the timeframe required to achieve isolation and, eventually, speciation in sister lineages. The application of species delimitation methods and an accurate dating strategy to the genus Neomys helped to clarify its controversial taxonomy.

Development of rapidly evolving intron markers to estimate multilocus species trees of rodents.

One of the major challenges in the analysis of closely related species, speciation and phylogeography is the identification of variable sequence markers that allow the determination of genealogical relationships in multiple genomic regions using coalescent and species tree approaches. Rodent species represent nearly half of the mammalian diversity, but so far no systematic study has been carried out to detect suitable informative markers for this group. Here, we used a bioinformatic pipeline to extract intron sequences from rodent genomes available in databases and applied a series of filters that allowed the identification of 208 introns that adequately fulfilled several criteria for these studies. The main required characteristics of the introns were that they had the maximum possible mutation rates, that they were part of single-copy genes, that they had an appropriate sequence length for amplification, and that they were flanked by exons with suitable regions for primer design. In addition, in order to determine the validity of this approach, we chose ten of these introns for primer design and tested them in a panel of eleven rodent species belonging to different representative families. We show that all these introns can be amplified in the majority of species and that, overall, 79% of the amplifications worked with minimum optimization of the annealing temperature. In addition, we confirmed for a pair of sister species the relatively high level of sequence divergence of these introns. Therefore, we provide here a set of adequate intron markers that can be applied to different species of Rodentia for their use in studies that require significant sequence variability.

Phylogeography and postglacial expansion of the endangered semi-aquatic mammal Galemys pyrenaicus.

BACKGROUND: Species with strict ecological requirements may provide new insights into the forces that shaped the geographic variation of genetic diversity. The Pyrenean desman, Galemys pyrenaicus, is a small semi-aquatic mammal that inhabits clean streams of the northern half of the Iberian Peninsula and is endangered in most of its geographic range, but its genetic structure is currently unknown. While the stringent ecological demands derived from its aquatic habitat might have caused a partition of the genetic diversity among river basins, Pleistocene glaciations would have generated a genetic pattern related to glacial refugia. RESULTS: To study the relative importance of historical and ecological factors in the genetic structure of G. pyrenaicus, we used mitochondrial and intronic sequences of specimens covering most of the species range. We show, first, that the Pyrenean desman has very low levels of genetic diversity compared to other mammals. In addition, phylogenetic and dating analyses of the mitochondrial sequences reveal a strong phylogeographic structure of a Middle Pleistocene origin, suggesting that the main lineages arose during periods of glacial isolation. Furthermore, both the spatial distribution of nuclear and mitochondrial diversity and the results of species distribution modeling suggest the existence of a major glacial refugium in the northwestern part of the Iberian Peninsula. Finally, the main mitochondrial lineages show a striking parapatric distribution without any apparent exchange of mitochondrial haplotypes between the lineages that came into secondary contact (although with certain permeability to nuclear genes), indicating incomplete mixing after the post-glacial recolonization. On the other hand, when we analyzed the partition of the genetic diversity among river basins, the Pyrenean desman showed a lower than expected genetic differentiation among main rivers. CONCLUSIONS: The analysis of mitochondrial and intronic markers in G. pyrenaicus showed the predominant effects of Pleistocene glaciations on the genetic structure of this species, while the distribution of the genetic diversity was not greatly influenced by the main river systems. These results and, particularly, the discovery of a marked phylogeographic structure, may have important implications for the conservation of the Pyrenean desman.

Diversification rates and the latitudinal gradient of diversity in mammals.

The latitudinal gradient of species richness has frequently been attributed to higher diversification rates of tropical groups. In order to test this hypothesis for mammals, we used a set of 232 genera taken from a mammalian supertree and, additionally, we reconstructed dated Bayesian phylogenetic trees of 100 genera. For each genus, diversification rate was estimated taking incomplete species sampling into account and latitude was assigned considering the heterogeneity in species distribution ranges. For both datasets, we found that the average diversification rate was similar among all latitudinal bands. Furthermore, when we used phylogenetically independent contrasts, we did not find any significant correlation between latitude and diversification parameters, including different estimates of speciation and extinction rates. Thus, other factors, such as the dynamics of dispersal through time, may be required to explain the latitudinal gradient of diversity in mammals.

Impact of deep coalescence on the reliability of species tree inference from different types of DNA markers in mammals.

An important challenge for phylogenetic studies of closely related species is the existence of deep coalescence and gene tree heterogeneity. However, their effects can vary between species and they are often neglected in phylogenetic analyses. In addition, a practical problem in the reconstruction of shallow phylogenies is to determine the most efficient set of DNA markers for a reliable estimation. To address these questions, we conducted a multilocus simulation study using empirical values of nucleotide diversity and substitution rates obtained from a wide range of mammals and evaluated the performance of both gene tree and species tree approaches to recover the known speciation times and topological relationships. We first show that deep coalescence can be a serious problem, more than usually assumed, for the estimation of speciation times in mammals using traditional gene trees. Furthermore, we tested the performance of different sets of DNA markers in the determination of species trees using a coalescent approach. Although the best estimates of speciation times were obtained, as expected, with the use of an increasing number of nuclear loci, our results show that similar estimations can be obtained with a much lower number of genes and the incorporation of a mitochondrial marker, with its high information content. Thus, the use of the combined information of both nuclear and mitochondrial markers in a species tree framework is the most efficient option to estimate recent speciation times and, consequently, the underlying species tree.

Novel intron markers to study the phylogeny of closely related mammalian species.

BACKGROUND: Multilocus phylogenies can be used to infer the species tree of a group of closely related species. In species trees, the nodes represent the actual separation between species, thus providing essential information about their evolutionary history. In addition, multilocus phylogenies can help in analyses of species delimitation, gene flow and genetic differentiation within species. However, few adequate markers are available for such studies. RESULTS: In order to develop nuclear markers that can be useful in multilocus studies of mammals, we analyzed the mammalian genomes of human, chimpanzee, macaque, dog and cow. Rodents were excluded due to their unusual genomic features. Introns were extracted from the mammalian genomes because of their greater genetic variability and ease of amplification from the flanking exons. To an initial set of more than 10,000 one-to-one orthologous introns we applied several filters to select introns that belong to single-copy genes, show neutral evolutionary rates and have an adequate length for their amplification. This analysis led to a final list of 224 intron markers randomly distributed along the genome. To experimentally test their validity, we amplified twelve of these introns in a panel of six mammalian species. The result was that seven of these introns gave rise to a PCR band of the expected size in all species. In addition, we sequenced these bands and analyzed the accumulation of substitutions in these introns in five pairs of closely related species. The results showed that the estimated genetic distances in the five species pairs was quite variable among introns and that this divergence cannot be directly predicted from the overall intron divergence in mammals. CONCLUSIONS: We have designed a new set of 224 nuclear introns with optimal features for the phylogeny of closely related mammalian species. A large proportion of the introns tested experimentally showed a perfect amplification and enough variability in most species, indicating that this marker set can be very helpful in multilocus phylogenetics of mammals. Due to the lower variability and stronger stochasticity of nuclear markers with respect to mitochondrial genes, studies should be designed to make use of several markers like the ones designed here.

Paleogenomics in a temperate environment: shotgun sequencing from an extinct Mediterranean caprine.

BACKGROUND: Numerous endemic mammals, including dwarf elephants, goats, hippos and deers, evolved in isolation in the Mediterranean islands during the Pliocene and Pleistocene. Most of them subsequently became extinct during the Holocene. Recently developed high-throughput sequencing technologies could provide a unique tool for retrieving genomic data from these extinct species, making it possible to study their evolutionary history and the genetic bases underlying their particular, sometimes unique, adaptations. METHODOLOGY/PRINCIPALS FINDINGS: A DNA extraction of a approximately 6,000 year-old bone sample from an extinct caprine (Myotragus balearicus) from the Balearic Islands in the Western Mediterranean, has been subjected to shotgun sequencing with the GS FLX 454 platform. Only 0.27% of the resulting sequences, identified from alignments with the cow genome and comprising 15,832 nucleotides, with an average length of 60 nucleotides, proved to be endogenous. CONCLUSIONS: A phylogenetic tree generated with Myotragus sequences and those from other artiodactyls displays an identical topology to that generated from mitochondrial DNA data. Despite being in an unfavourable thermal environment, which explains the low yield of endogenous sequences, our study demonstrates that it is possible to obtain genomic data from extinct species from temperate regions.

Distribution and evolutionary trends of photoprotective isoprenoids (xanthophylls and tocopherols) within the plant kingdom.

The earliest land photosynthesis would have increased the risk of photo-oxidations and the demand of anti-oxidative protection. In this work, we aimed to determine the evolutionary trends in photoprotection across a wide representation of the plant kingdom and to verify whether the non-ubiquitous lutein-epoxide (Lx) cycle is a polyphyletic or an ancient character. Carotenoids and alpha-tocopherol (alpha-toc) were analysed by HPLC in 266 species. Phylogenetic analyses of the presence of photoprotective compounds and zeaxanthin-epoxidase (ZE) sequences were performed. Violaxanthin-cycle pigments (VAZ) and alpha-toc were taxonomically ubiquitous. Ancient groups showed higher contents of VAZ than vascular plants, while alpha-toc showed the opposite pattern. Lutein-epoxide was present in 45% of the species. It showed a remarkable variation across groups but with a clear increasing trend from algae to basal angiosperms. Lutein-epoxide was also related to woody trait and leaf longevity. No correlation between the presence of Lx and recurrent mutations in ZE sequences, including the duplications, was found. Thus, there is an evolutionary trend to increase the content of alpha-toc and to decrease the total amount of VAZ pigments. Absence of Lx in algae discards an ancestral origin. Present results are also inconsistent with a polyphyletic origin of Lx in angiosperms.

Estimation of phylogenetic inconsistencies in the three domains of life.

Discrepancies in phylogenetic trees of bacteria and archaea are often explained as lateral gene transfer events. However, such discrepancies may also be due to phylogenetic artifacts or orthology assignment problems. A first step that may help to resolve this dilemma is to estimate the extent of phylogenetic inconsistencies in trees of prokaryotes in comparison with those of higher eukaryotes, where no lateral gene transfer is expected. To test this, we used 21 proteomes each of eukaryotes (mainly opisthokonts), proteobacteria, and archaea that spanned equivalent levels of genetic divergence. In each domain of life, we defined a set of putative orthologous sequences using a phylogenetic-based orthology protocol and, as a reference topology, we used a tree constructed with concatenated genes of each domain. Our results show, for most of the tests performed, that the magnitude of topological inconsistencies with respect to the reference tree was very similar in the trees of proteobacteria and eukaryotes. When clade support was taken into account, prokaryotes showed some more inconsistencies, but then all values were very low. Discrepancies were only consistently higher in archaea but, as shown by simulation analysis, this is likely due to the particular tree of the archaeal species used here being more difficult to reconstruct, whereas the trees of proteobacteria and eukaryotes were of similar difficulty. Although these results are based on a relatively small number of genes, it seems that phylogenetic reconstruction problems, including orthology assignment problems, have a similar overall effect over prokaryotic and eukaryotic trees based on single genes. Consequently, lateral gene transfer between distant prokaryotic species may have been more rare than previously thought, which opens the way to obtain the tree of life of bacterial and archaeal species using genomic data and the concatenation of adequate genes, in the same way as it is usually done in eukaryotes.