First assessment of the biodiversity of praying mantises (Insecta: Mantodea) in Cameroon with DNA barcoding.

Praying mantises are the apex insect predators in many ecosystems, nevertheless they receive relatively less recognition in biodiversity reviews. We report a first survey of diversity of praying mantises in Cameroon, which is situated in the Congo Basin region, one of the richest biodiversity hotspots. Combination of light trapping with manual collecting resulted in 495 specimens representing 62 species. A total of eight species are novel for the country, at least five species are likely undescribed. DNA barcodes of 72 specimens representing every collected species were obtained, curated, and submitted to NCBI database. For eight species, barcodes are published for the first time. A maximum likelihood phylogenetic tree was created using all available barcodes of Mantodea of Central African subregion. The results obtained during this study stress the importance of combining traditional and molecular approaches during biodiversity assessments of often neglected taxa, the latter aiding in uncovering new species, resolving unknown morphological divergencies and assigning conspecifics.

Evaluation of the Impact of Population Management on the Genetic Parameters of Selected Spiral-Horned Antelopes.

The rapid loss of biodiversity and the associated reduction and fragmentation of habitats means that ex situ populations have become an important part of species conservation. These populations, which are often established from a small number of founders, require careful management to avoid the negative effects of genetic drift and inbreeding. Although the inclusion of molecular data is recommended, their availability for captive breeding management remains limited. The aim of this study was to evaluate the relationship between the levels of genetic diversity in six spiral-horned antelope taxa bred under human care and their respective management strategies, conservation status, demography, and geographic origin, using 10 nuclear DNA microsatellite loci and mitochondrial control region DNA sequences. Our findings include associations between genetic diversity and management intensity but also with the diversity and contribution of wild populations to captive founders, with some populations apparently composed of animals from divergent wild lineages elevating captive genetic diversity. When population sizes are large, the potential advantages of maximizing genetic diversity in widely outcrossed populations may need careful consideration with respect to the potential disruption of adaptive diversity. Genetic data serve as a robust tool for managing captive populations, yet their interpretation necessitates a comprehensive understanding of species biology and history.

Genomic analyses reveal poaching hotspots and illegal trade in pangolins from Africa to Asia.

The white-bellied pangolin (Phataginus tricuspis) is the world's most trafficked mammal and is at risk of extinction. Reducing the illegal wildlife trade requires an understanding of its origins. Using a genomic approach for tracing confiscations and analyzing 111 samples collected from known geographic localities in Africa and 643 seized scales from Asia between 2012 and 2018, we found that poaching pressures shifted over time from West to Central Africa. Recently, Cameroon's southern border has emerged as a site of intense poaching. Using data from seizures representing nearly 1 million African pangolins, we identified Nigeria as one important hub for trafficking, where scales are amassed and transshipped to markets in Asia. This origin-to-destination approach offers new opportunities to disrupt the illegal wildlife trade and to guide anti-trafficking measures.

The extinct Sicilian wolf shows a complex history of isolation and admixture with ancient dogs.

The Sicilian wolf remained isolated in Sicily from the end of the Pleistocene until its extermination in the 1930s-1960s. Given its long-term isolation on the island and distinctive morphology, the genetic origin of the Sicilian wolf remains debated. We sequenced four nuclear genomes and five mitogenomes from the seven existing museum specimens to investigate the Sicilian wolf ancestry, relationships with extant and extinct wolves and dogs, and diversity. Our results show that the Sicilian wolf is most closely related to the Italian wolf but carries ancestry from a lineage related to European Eneolithic and Bronze Age dogs. The average nucleotide diversity of the Sicilian wolf was half of the Italian wolf, with 37-50% of its genome contained in runs of homozygosity. Overall, we show that, by the time it went extinct, the Sicilian wolf had high inbreeding and low-genetic diversity, consistent with a population in an insular environment.

Contrasting levels of hybridization across the two contact zones between two hedgehog species revealed by genome-wide SNP data.

Hybridization and introgression have played important roles in the history of various species, including lineage diversification and the evolution of adaptive traits. Hybridization can accelerate the development of reproductive isolation between diverging species, and thus valuable insight into the evolution of reproductive barrier formation may be gained by studying secondary contact zones. Hedgehogs of the genus Erinaceus, which are insectivores sensitive to changes in climate, are a pioneer model in Pleistocene phylogeography. The present study provides the first genome-wide SNP data regarding the Erinaceus hedgehogs species complex, offering a unique comparison of two secondary contact zones between Erinaceus europaeus and E. roumanicus. Results confirmed diversification of the genus during the Pleistocene period, and detected a new refugial lineage of E. roumanicus outside the Mediterranean basin, most likely in the Ponto-Caspian region. In the Central European zone, the level of hybridization was low, whereas in the Russian-Baltic zone, both species hybridise extensively. Asymmetrical gene flow from E. europaeus to E. roumanicus suggests that reproductive isolation varies according to the direction of the crosses in the hybrid zones. However, no loci with significantly different patterns of introgression were detected. Markedly different pre- and post-zygotic barriers, and thus diverse modes of species boundary maintenance in the two contact zones, likely exist. This pattern is probably a consequence of the different age and thus of the different stage of evolution of reproductive isolating mechanisms in each hybrid zone.

Genomic and Transcriptomic Characterization of Atypical Recurrent Flank Alopecia in the Cesky Fousek.

Non-inflammatory alopecia is a frequent skin problem in dogs, causing damaged coat integrity and compromised appearance of affected individuals. In this study, we examined the Cesky Fousek breed, which displays atypical recurrent flank alopecia (aRFA) at a high frequency. This type of alopecia can be quite severe and is characterized by seasonal episodes of well demarcated alopecic areas without hyperpigmentation. The genetic component responsible for aRFA remains unknown. Thus, here we aimed to identify variants involved in aRFA using a combination of histological, genomic, and transcriptomic data. We showed that aRFA is histologically similar to recurrent flank alopecia, characterized by a lack of anagen hair follicles and the presence of severely shortened telogen or kenogen hair follicles. We performed a genome-wide association study (GWAS) using 216 dogs phenotyped for aRFA and identified associations on chromosomes 19, 8, 30, 36, and 21, highlighting 144 candidate genes, which suggests a polygenic basis for aRFA. By comparing the skin cell transcription pattern of six aRFA and five control dogs, we identified 236 strongly differentially expressed genes (DEGs). We showed that the GWAS genes associated with aRFA are often predicted to interact with DEGs, suggesting their joint contribution to the development of the disease. Together, these genes affect four major metabolic pathways connected to aRFA: collagen formation, muscle structure/contraction, lipid metabolism, and the immune system.

Emergence of methicillin resistance predates the clinical use of antibiotics.

The discovery of antibiotics more than 80 years ago has led to considerable improvements in human and animal health. Although antibiotic resistance in environmental bacteria is ancient, resistance in human pathogens is thought to be a modern phenomenon that is driven by the clinical use of antibiotics1. Here we show that particular lineages of methicillin-resistant Staphylococcus aureus-a notorious human pathogen-appeared in European hedgehogs in the pre-antibiotic era. Subsequently, these lineages spread within the local hedgehog populations and between hedgehogs and secondary hosts, including livestock and humans. We also demonstrate that the hedgehog dermatophyte Trichophyton erinacei produces two ß-lactam antibiotics that provide a natural selective environment in which methicillin-resistant S. aureus isolates have an advantage over susceptible isolates. Together, these results suggest that methicillin resistance emerged in the pre-antibiotic era as a co-evolutionary adaptation of S. aureus to the colonization of dermatophyte-infected hedgehogs. The evolution of clinically relevant antibiotic-resistance genes in wild animals and the connectivity of natural, agricultural and human ecosystems demonstrate that the use of a One Health approach is critical for our understanding and management of antibiotic resistance, which is one of the biggest threats to global health, food security and development.

Identification of Tapeworm Species in Genetically Characterised Grey Wolves Recolonising Central Europe.

PURPOSE: Restored role of the grey wolf in ecological networks of newly recolonized areas can be studied via surveys of parasite communities of this predator. As helminths circulating in multi-host systems, the tapeworms directly reflect wolves' diet, while some species are also important from the One Health perspective. The Czech experienced centuries of wolves' absence, however, now it is situated on the crossroad of recolonising wolves' populations, which is opening questions of their role in ecological networks in this area and thus in sylvatic cycles of heteroxenous parasites. METHODS: Five wolf carcasses from this area were obtained and genetic affinity to a particular population was inspected. Tapeworms isolated from wolves' intestines during necropsies were molecularly identified based on sequences of COI marker. RESULTS: Three wolf haplotypes (w1, w2, w14) correspond with the dominance of haplogroup 1 (w1, w2) within Central European lowland population and haplogroup 2 (w14) within the Carpathian population. Two Taenia spp. were revealed: T. krabbei in Central European population wolves and T. hydatigena in an individual from Carpathian population. CONCLUSIONS: The results serve as a base for future monitoring and studies of the recolonising wolf population and its impact on ecosystems in the studied area to contribute to the hypothesis about differentiation of parasite communities in particular wolf population and higher parasite diversity and richness in established populations in comparison to newly settled ones.

Hedgehogs, Squirrels, and Blackbirds as Sentinel Hosts for Active Surveillance of Borrelia miyamotoi and Borrelia burgdorferi Complex in Urban and Rural Environments.

Lyme borreliosis (LB), caused by spirochetes of the Borrelia burgdorferi sensu lato (s.l.) complex, is one of the most common vector-borne zoonotic diseases in Europe. Knowledge about the enzootic circulation of Borrelia pathogens between ticks and their vertebrate hosts is epidemiologically important and enables assessment of the health risk for the human population. In our project, we focused on the following vertebrate species: European hedgehog (Erinaceus europaeus), Northern white-breasted hedgehog (E. roumanicus), Eurasian red squirrel (Sciurus vulgaris), and Common blackbird (Turdus merula). The cadavers of accidentally killed animals used in this study constitute an available source of biological material, and we have confirmed its potential for wide monitoring of B. burgdorferi s.l. presence and genospecies diversity in the urban environment. High infection rates (90% for E. erinaceus, 73% for E. roumanicus, 91% for S. vulgaris, and 68% for T. merula) were observed in all four target host species; mixed infections by several genospecies were detected on the level of individuals, as well as in particular tissue samples. These findings show the usefulness of multiple tissue sampling as tool for revealing the occurrence of several genospecies within one animal and the risk of missing particular B. burgdorferi s.l. genospecies when looking in one organ alone.

3D Geometric Morphometrics Reveals Convergent Character Displacement in the Central European Contact Zone between Two Species of Hedgehogs (Genus Erinaceus).

Hedgehogs, as medium-sized plantigrade insectivores with low basal metabolic rates and related defensive anti-predator strategies, are quite sensitive to temperature and ecosystem productivity. Their ranges therefore changed dramatically due to Pleistocene climate oscillations, resulting in allopatric speciation and the subsequent formation of secondary contact zones. Such interactions between closely related species are known to generate strong evolutionary forces responsible for niche differentiation. In this connection, here, we detail the results of research on the phenotypic evolution in the two species of hedgehog present in central Europe, as based on genetics and geometric morphometrics in samples along a longitudinal transect that includes the contact zone between the species. While in allopatry, Erinaceus europaeus is found to have a larger skull than E. roumanicus and distinct cranial and mandibular shapes; the members of the two species in sympatry are smaller and more similar to each other, with a convergent shape of the mandible. The relevant data fail to reveal any major role for either hybridisation or clinal variation. We, therefore, hypothesise that competitive pressure exerted on the studied species does not generate divergent selection sufficient for divergent character displacement to evolve, instead giving rise to convergent selection in the face of resource limitation in the direction of smaller skull size. Considering the multi-factorial constraints present in the relevant adaptive landscape, reduction in size could also be facilitated by predator pressure in ecosystems characterised by mesopredator release and other anthropogenic factors. As the function of the animals' lower jaw is mainly connected with feeding (in contrast to the cranium whose functions are obviously more complex), we interpret the similarity in shape as reflecting local adaptations to overlapping dietary resources in the two species and hence as convergent character displacement.

Evaluating the Impact of Wildlife Shelter Management on the Genetic Diversity of Erinaceus europaeus and E. roumanicus in Their Contact Zone.

Hedgehogs are among the most abundant species to be found within wildlife shelters and after successful rehabilitation they are frequently translocated. The effects and potential impact of these translocations on gene flow within wild populations are largely unknown. In this study, different wild hedgehog populations were compared with artificially created "shelter populations", with regard to their genetic diversity, in order to establish basic data for future inferences on the genetic impact of hedgehog translocations. Observed populations are located within central Europe, including the species Erinaceus europaeus and E. roumanicus. Shelters were mainly hosting one species; in one case, both species were present syntopically. Apart from one exception, the results did not show a higher genetic diversity within shelter populations, indicating that individuals did not originate from a wider geographical area than individuals grouped into one of the wild populations. Two shelters from Innsbruck hosted individuals that belonged to two potential clusters, as indicated in a distance analysis. When such a structure stems from the effects of landscape elements like large rivers, the shelter management-related translocations might lead to homogenization across the dispersal barrier.

Comparison of divergent breeding management strategies in two species of semi-captive eland in Senegal.

Breeding management of small populations may have a critical influence on the development of population characteristics in terms of levels of genetic diversity and inbreeding. Two populations of antelope sister species - Critically Endangered Western Derby eland (Tauroragus derbianus derbianus) and non-native Least Concern Cape eland (Taurotragus oryx oryx) bred under different management strategies were studied in Senegal, Western Africa. The aims of the study were to compare the population genetic parameters of the two species and to test for the presence of interspecific hybrids. In total, blood and tissue samples from 76 Western Derby elands and 26 Cape elands were investigated, using 12 microsatellite markers. No hybrid individuals were detected in the sampled animals within the multispecies enclosure in Bandia Reserve, Senegal. The parameters of genetic polymorphism indicated much lower genetic diversity in Western Derby elands compared to Cape elands. On the other hand, the coefficient of inbreeding was low in both species. It is hypothesized that this could be a positive effect of strict population management of Western Derby elands, which, despite the loss of genetic diversity, minimizes inbreeding.

Dynamic range expansion leads to establishment of a new, genetically distinct wolf population in Central Europe.

Local extinction and recolonization events can shape genetic structure of subdivided animal populations. The gray wolf (Canis lupus) was extirpated from most of Europe, but recently recolonized big part of its historical range. An exceptionally dynamic expansion of wolf population is observed in the western part of the Great European Plain. Nonetheless, genetic consequences of this process have not yet been fully understood. We aimed to assess genetic diversity of this recently established wolf population in Western Poland (WPL), determine its origin and provide novel data regarding the population genetic structure of the grey wolf in Central Europe. We utilized both spatially explicit and non-explicit Bayesian clustering approaches, as well as a model-independent, multivariate method DAPC, to infer genetic structure in large dataset (881 identified individuals) of wolf microsatellite genotypes. To put the patterns observed in studied population into a broader biogeographic context we also analyzed a mtDNA control region fragment widely used in previous studies. In comparison to a source population, we found slightly reduced allelic richness and heterozygosity in the newly recolonized areas west of the Vistula river. We discovered relatively strong west-east structuring in lowland wolves, probably reflecting founder-flush and allele surfing during range expansion, resulting in clear distinction of WPL, eastern lowland and Carpathian genetic groups. Interestingly, wolves from recently recolonized mountainous areas (Sudetes Mts, SW Poland) clustered together with lowland, but not Carpathian wolf populations. We also identified an area in Central Poland that seems to be a melting pot of western, lowland eastern and Carpathian wolves. We conclude that the process of dynamic recolonization of Central European lowlands lead to the formation of a new, genetically distinct wolf population. Together with the settlement and establishment of packs in mountains by lowland wolves and vice versa, it suggests that demographic dynamics and possibly anthropogenic barriers rather than ecological factors (e.g. natal habitat-biased dispersal patterns) shape the current wolf genetic structure in Central Europe.

Tracing genetic resurrection of pointing dog breeds: Cesky Fousek as both survivor and rescuer.

Cesky Fousek is considered to be one of the oldest pointing dog breeds in Europe and has been appreciated for its versatile working skills. Because it faced extinction in the past, the Cesky Fousek was restored from German Wirehaired and Shorthaired Pointers. Additionally, the breed was recently used in the USA with the initial intent of improvement of the Wirehaired Pointing Griffon (synonymous with Korthals Griffon) by the Bohemian Wirehaired Pointing Griffon Club of America. This study evaluates genetic diversity parameters of Cesky Fousek and compares them to the other continental pointing dogs that played a role in the formation of its gene pool. DNA from buccal swab and blood samples (n = 405) were analyzed using 18 microsatellite markers. Parameters of genetic polymorphism show that the Cesky Fousek breed has a comparable rate of variation as other hunting breeds despite the low population size and severe historical bottlenecks. Clustering analyses reveal a unique genetic status as a distinct pointing dog breed and the relatedness of the breeds is in good concordance with historical data. The present study demonstrates that despite historical admixture among lineages, separate pointing breeds constitute genetically differentiated units, mirroring unique breeding stocks and pedigree isolation among specific breed clubs, reflecting differences in breeding programs under each association.

Wolf outside, dog inside? The genomic make-up of the Czechoslovakian Wolfdog.

BACKGROUND: Genomic methods can provide extraordinary tools to explore the genetic background of wild species and domestic breeds, optimize breeding practices, monitor and limit the spread of recessive diseases, and discourage illegal crossings. In this study we analysed a panel of 170k Single Nucleotide Polymorphisms with a combination of multivariate, Bayesian and outlier gene approaches to examine the genome-wide diversity and inbreeding levels in a recent wolf x dog cross-breed, the Czechoslovakian Wolfdog, which is becoming increasingly popular across Europe. RESULTS: Pairwise FST values, multivariate and assignment procedures indicated that the Czechoslovakian Wolfdog was significantly differentiated from all the other analysed breeds and also well-distinguished from both parental populations (Carpathian wolves and German Shepherds). Coherently with the low number of founders involved in the breed selection, the individual inbreeding levels calculated from homozygosity regions were relatively high and comparable with those derived from the pedigree data. In contrast, the coefficient of relatedness between individuals estimated from the pedigrees often underestimated the identity-by-descent scores determined using genetic profiles. The timing of the admixture and the effective population size trends estimated from the LD patterns reflected the documented history of the breed. Ancestry reconstruction methods identified more than 300 genes with excess of wolf ancestry compared to random expectations, mainly related to key morphological features, and more than 2000 genes with excess of dog ancestry, playing important roles in lipid metabolism, in the regulation of circadian rhythms, in learning and memory processes, and in sociability, such as the COMT gene, which has been described as a candidate gene for the latter trait in dogs. CONCLUSIONS: In this study we successfully applied genome-wide procedures to reconstruct the history of the Czechoslovakian Wolfdog, assess individual wolf ancestry proportions and, thanks to the availability of a well-annotated reference genome, identify possible candidate genes for wolf-like and dog-like phenotypic traits typical of this breed, including commonly inherited disorders. Moreover, through the identification of ancestry-informative markers, these genomic approaches could provide tools for forensic applications to unmask illegal crossings with wolves and uncontrolled trades of recent and undeclared wolfdog hybrids.

Glacial allopatry vs. postglacial parapatry and peripatry: the case of hedgehogs.

Although hedgehogs are well-known examples of postglacial recolonisation, the specific processes that shape their population structures have not been examined by detailed sampling and fast-evolving genetic markers in combination with model based clustering methods. This study aims to analyse the impacts of isolation within glacial refugia and of postglacial expansion on the population structure of the Northern White-breasted hedgehog (Erinaceus roumanicus). It also discusses the role of the processes at edges of species distribution in its evolutionary history. The maternally inherited mitochondrial control region and the bi-parentally inherited nuclear microsatellites were used to examine samples within the Central Europe, Balkan Peninsula and adjacent islands. Bayesian coalescent inference and neutrality tests proposed a recent increase in the population size. The most pronounced pattern of population structure involved differentiation of the insular populations in the Mediterranean Sea and the population within the contact zone with E. europaeus in Central Europe. An interspecies hybrid was detected for the first time in Central Europe. A low genetic diversity was observed in Crete, while the highest genetic distances among individuals were found in Romania. The recent population in the post-refugial area related to the Balkan Peninsula shows a complex pattern with pronounced subpopulations located mainly in the Pannonian Basin and at the Adriatic and Pontic coasts. Detailed analyses indicate that parapatry and peripatry may not be the only factors that limit range expansion, but also strong microevolutionary forces that may change the genetic structure of the species. Here we present evidence showing that population differentiation may occur not only during the glacial restriction of the range into the refugia, but also during the interglacial range expansion. Population differentiation at the Balkan Peninsula and adjacent regions could be ascribed to diversification in steppe/forest biomes and complicated geomorphology, including pronounced geographic barriers as Carpathians.

Combining phylogenetic and demographic inferences to assess the origin of the genetic diversity in an isolated wolf population.

The survival of isolated small populations is threatened by both demographic and genetic factors. Large carnivores declined for centuries in most of Europe due to habitat changes, overhunting of their natural prey and direct persecution. However, the current rewilding trends are driving many carnivore populations to expand again, possibly reverting the erosion of their genetic diversity. In this study we reassessed the extent and origin of the genetic variation of the Italian wolf population, which is expanding after centuries of decline and isolation. We genotyped wolves from Italy and other nine populations at four mtDNA regions (control-region, ATP6, COIII and ND4) and 39 autosomal microsatellites. Results of phylogenetic analyses and assignment procedures confirmed in the Italian wolves a second private mtDNA haplotype, which belongs to a haplogroup distributed mostly in southern Europe. Coalescent analyses showed that the unique mtDNA haplotypes in the Italian wolves likely originated during the late Pleistocene. ABC simulations concordantly showed that the extant wolf populations in Italy and in south-western Europe started to be isolated and declined right after the last glacial maximum. Thus, the standing genetic variation in the Italian wolves principally results from the historical isolation south of the Alps.

Vocal activity of lesser galagos (Galago spp.) at zoos.

Almost nothing is known about the natural vocal behavior of lesser galagos living in zoos. This is perhaps because they are usually kept in nocturnal exhibits separated from the visitors by a transparent and acoustically insulating glass barrier. The aim of the present study was therefore to fill this gap in knowledge of the vocal behavior of lesser galagos from zoos. This knowledge might be beneficial because the vocalizations of these small primates can be used for species determination. We performed a 10-day-long acoustic monitoring of vocal activity in each of seven various groups of Galago senegalensis and G. moholi living at four zoos. We quantitatively evaluated the occurrence of four loud vocalization types present in both species, including the most species-specific advertisement call. We found that qualitative as well as quantitative differences exist in the vocal behavior of the studied groups. We confirmed that the observed vocalization types can be collected from lesser galagos living at zoos, and the success can be increased by selecting larger and more diverse groups. We found two distinct patterns of diel vocal activity in the most vocally active groups. G. senegalensis groups were most vocally active at the beginning and at the end of their activity period, whereas one G. moholi group showed an opposite pattern. The latter is surprising, as it is generally accepted that lesser galagos emit advertisement calls especially at dawn and dusk, i.e., at the beginning and at the end of their diel activity.

From Wolves to Dogs, and Back: Genetic Composition of the Czechoslovakian Wolfdog.

The Czechoslovakian Wolfdog is a unique dog breed that originated from hybridization between German Shepherds and wild Carpathian wolves in the 1950s as a military experiment. This breed was used for guarding the Czechoslovakian borders during the cold war and is currently kept by civilian breeders all round the world. The aim of our study was to characterize, for the first time, the genetic composition of this breed in relation to its known source populations. We sequenced the hypervariable part of the mtDNA control region and genotyped the Amelogenin gene, four sex-linked microsatellites and 39 autosomal microsatellites in 79 Czechoslovakian Wolfdogs, 20 German Shepherds and 28 Carpathian wolves. We performed a range of population genetic analyses based on both empirical and simulated data. Only two mtDNA and two Y-linked haplotypes were found in Czechoslovakian Wolfdogs. Both mtDNA haplotypes were of domestic origin, while only one of the Y-haplotypes was shared with German Shepherds and the other was unique to Czechoslovakian Wolfdogs. The observed inbreeding coefficient was low despite the small effective population size of the breed, possibly due to heterozygote advantages determined by introgression of wolf alleles. Moreover, Czechoslovakian Wolfdog genotypes were distinct from both parental populations, indicating the role of founder effect, drift and/or genetic hitchhiking. The results revealed the peculiar genetic composition of the Czechoslovakian Wolfdog, showing a limited introgression of wolf alleles within a higher proportion of the dog genome, consistent with the reiterated backcrossing used in the pedigree. Artificial selection aiming to keep wolf-like phenotypes but dog-like behavior resulted in a distinctive genetic composition of Czechoslovakian Wolfdogs, which provides a unique example to study the interactions between dog and wolf genomes.

A deep phylogeny of viral and cellular right-hand polymerases.

Right-hand polymerases are important players in genome replication and repair in cellular organisms as well as in viruses. All right-hand polymerases are grouped into seven related protein families: viral RNA-dependent RNA polymerases, reverse transcriptases, single-subunit RNA polymerases, and DNA polymerase families A, B, D, and Y. Although the evolutionary relationships of right-hand polymerases within each family have been proposed, evolutionary relationships between families remain elusive because their sequence similarity is too low to allow classical phylogenetic analyses. The structure of viral RNA-dependent RNA polymerases recently was shown to be useful in inferring their evolution. Here, we address evolutionary relationships between right-hand polymerase families by combining sequence and structure information. We used a set of 22 viral and cellular polymerases representing all right-hand polymerase families with known protein structure. In contrast to previous studies, which focused only on the evolution of particular families, the current approach allowed us to present the first robust phylogenetic analysis unifying evolution of all right-hand polymerase families. All polymerase families branched into discrete lineages, following a fairly robust adjacency pattern. Only single-subunit RNA polymerases formed an inner group within DNA polymerase family A. RNA-dependent RNA polymerases of RNA viruses and reverse transcriptases of retroviruses formed two sister groups and were distinguishable from all other polymerases. DNA polymerases of DNA bacteriophages did not form a monophyletic group and are phylogenetically mixed with cellular DNA polymerase families A and B. Based on the highest genetic variability and structural simplicity, we assume that RNA-dependent RNA polymerases are the most ancient group of right-hand polymerases, in agreement with the RNA World hypothesis, because RNA-dependent RNA polymerases are enzymes that could serve in replication of RNA genomes. Moreover, our results show that protein structure can be used in phylogenetic analyses of distantly related proteins that share only limited sequence similarity.