Cystic echinococcosis in Iceland: a brief history and genetic analysis of a 46-year-old Echinococcus isolate collected prior to the eradication of this zoonotic disease.

Cystic echinococcosis (CE) is considered the most severe parasitic disease that ever affected the human population in Iceland. Before the start of eradication campaign in the 1860s, Iceland was a country with very high prevalence of human CE, with approximately every fifth person infected. Eradication of CE from Iceland by 1979 was a huge success story and served as a leading example for other countries on how to combat such a severe One Health problem. However, there is no genetic information on Echinococcus parasites before eradication. Here, we reveal the genetic identity for one of the last Echinococcus isolates in Iceland, obtained from a sheep 46 years ago (1977). We sequenced a large portion of the mitochondrial genome (8141 bp) and identified the isolate as Echinococcus granulosus sensu stricto genotype G1. As G1 is known to be highly infective genotype to humans, it may partly explain why such a large proportion of human population in Iceland was infected at a time . The study demonstrates that decades-old samples hold significant potential to uncover genetic identities of parasites in the past.

Exploring the genetic diversity of genotypes G8 and G10 of the Echinococcus canadensis cluster in Europe based on complete mitochondrial genomes (13 550-13 552 bp).

Echinococcus granulosus sensu lato is a group of tapeworm species known to cause cystic echinococcosis. Within this group, the Echinococcus canadensis cluster includes genotypes G8 and G10 that have a predominantly sylvatic life cycle ­ transmission occurs between wild cervids and wolves. Relatively few studies have explored the genetic variation of the elusive G8 and G10, and their extent of genetic variation is yet to be investigated at the complete mitochondrial (mt) genome level. The aim was to explore the genetic variation of these 2 genotypes in Europe using complete mtDNA sequences and provide a high-quality reference dataset for future studies. Sequences of complete mt genomes were produced for 29 samples of genotype G8 and G10 from wolves, moose, reindeer and roe deer, originating from Finland, Sweden, Russia, Poland, Latvia and Estonia. Genetic variation was explored based on phylogenetic network analysis, revealing marked differences between G8 and G10 (over 400 mutations), and more detailed patterns of variability within the 2 genotypes than previously observed. Understanding the mt genetic composition of a species provides a baseline for future studies aiming to understand whether this mt distinctiveness is mirrored in the nuclear genome and whether it has any impact on any phenotypic traits or parasite transmission.

Genome-wide exploration reveals distinctive northern and southern variants of Clonorchis sinensis in the Far East.

Clonorchis sinensis is a carcinogenic liver fluke that causes clonorchiasis-a neglected tropical disease (NTD) affecting ~35 million people worldwide. No vaccine is available, and chemotherapy relies on one anthelmintic, praziquantel. This parasite has a complex life history and is known to infect a range of species of intermediate (freshwater snails and fish) and definitive (piscivorous) hosts. Despite this biological complexity and the impact of this biocarcinogenic pathogen, there has been no previous study of molecular variation in this parasite on a genome-wide scale. Here, we conducted the first extensive nuclear genomic exploration of C. sinensis individuals (n = 152) representing five distinct populations from mainland China, and one from Far East Russia, and revealed marked genetic variation within this species between "northern" and "southern" geographical regions. The discovery of this variation indicates the existence of biologically distinct variants within C. sinensis, which may have distinct epidemiology, pathogenicity and/or chemotherapic responsiveness. The detection of high heterozygosity within C. sinensis specimens suggests that this parasite has developed mechanisms to readily adapt to changing environments and/or host species during its life history/evolution. From an applied perspective, the identification of invariable genes could assist in finding new intervention targets in this parasite, given the major clinical relevance of clonorchiasis. From a technical perspective, the genomic-informatic workflow established herein will be readily applicable to a wide range of other parasites that cause NTDs.

Chromosome-scale Echinococcus granulosus (genotype G1) genome reveals the Eg95 gene family and conservation of the EG95-vaccine molecule.

Cystic echinococcosis is a socioeconomically important parasitic disease caused by the larval stage of the canid tapeworm Echinococcus granulosus, afflicting millions of humans and animals worldwide. The development of a vaccine (called EG95) has been the most notable translational advance in the fight against this disease in animals. However, almost nothing is known about the genomic organisation/location of the family of genes encoding EG95 and related molecules, the extent of their conservation or their functions. The lack of a complete reference genome for E. granulosus genotype G1 has been a major obstacle to addressing these areas. Here, we assembled a chromosomal-scale genome for this genotype by scaffolding to a high quality genome for the congener E. multilocularis, localised Eg95 gene family members in this genome, and evaluated the conservation of the EG95 vaccine molecule. These results have marked implications for future explorations of aspects such as developmentally-regulated gene transcription/expression (using replicate samples) for all E. granulosus stages; structural and functional roles of non-coding genome regions; molecular 'cross-talk' between oncosphere and the immune system; and defining the precise function(s) of EG95. Applied aspects should include developing improved tools for the diagnosis and chemotherapy of cystic echinococcosis of humans.

"Escalibur"-A practical pipeline for the de novo analysis of nucleotide variation in nonmodel eukaryotes.

The revolution in genomics has enabled large-scale population genetic investigations of a wide range of organisms, but there has been a relatively limited focus on improving analytical pipelines. To efficiently analyse large data sets, highly integrated and automated software pipelines, which are easy to use, efficient, reliable, reproducible and run in multiple computational environments, are required. A number of software workflows have been developed to handle and process such data sets for population genetic analyses, but effective, specialized pipelines for genetic and statistical analyses of nonmodel organisms are lacking. For most species, resources for variomes (sets of genetic variations found in populations of species) are not available, and/or genome assemblies are often incomplete and fragmented, complicating the selection of the most suitable reference genome when multiple assemblies are available. Additionally, the biological samples used often contain extraneous DNA from sources other than the species under investigation (e.g., microbial contamination), which needs to be removed prior to genetic analyses. For these reasons, we established a new pipeline, called Escalibur, which includes: functionalities, such as data trimming and mapping; selection of a suitable reference genome; removal of contaminating read data; recalibration of base calls; and variant-calling. Escalibur uses a proven gatk variant caller and workflow description language (WDL), and is, therefore, a highly efficient and scalable pipeline for the genome-wide identification of nucleotide variation in eukaryotes. This pipeline is available at https://gitlab.unimelb.edu.au/bioscience/escalibur (version 0.3-beta) and is essentially applicable to any prokaryote or eukaryote.

Nanopore Sequencing Resolves Elusive Long Tandem-Repeat Regions in Mitochondrial Genomes.

Long non-coding, tandem-repetitive regions in mitochondrial (mt) genomes of many metazoans have been notoriously difficult to characterise accurately using conventional sequencing methods. Here, we show how the use of a third-generation (long-read) sequencing and informatic approach can overcome this problem. We employed Oxford Nanopore technology to sequence genomic DNAs from a pool of adult worms of the carcinogenic parasite, Schistosoma haematobium, and used an informatic workflow to define the complete mt non-coding region(s). Using long-read data of high coverage, we defined six dominant mt genomes of 33.4 kb to 22.6 kb. Although no variation was detected in the order or lengths of the protein-coding genes, there was marked length (18.5 kb to 7.6 kb) and structural variation in the non-coding region, raising questions about the evolution and function of what might be a control region that regulates mt transcription and/or replication. The discovery here of the largest tandem-repetitive, non-coding region (18.5 kb) in a metazoan organism also raises a question about the completeness of some of the mt genomes of animals reported to date, and stimulates further explorations using a Nanopore-informatic workflow.

High-quality reference genome for Clonorchis sinensis.

The Chinese liver fluke, Clonorchis sinensis, causes the disease clonorchiasis, affecting ~35 million people in regions of China, Vietnam, Korea and the Russian Far East. Chronic clonorchiasis causes cholangitis and can induce a malignant cancer, called cholangiocarcinoma, in the biliary system. Control in endemic regions is challenging, and often relies largely on chemotherapy with one anthelmintic, called praziquantel. Routine treatment carries a significant risk of inducing resistance to this anthelmintic in the fluke, such that the discovery of new interventions is considered important. It is hoped that the use of molecular technologies will assist this endeavour by enabling the identification of drug or vaccine targets involved in crucial biological processes and/or pathways in the parasite. Although draft genomes of C. sinensis have been published, their assemblies are fragmented. In the present study, we tackle this genome fragmentation issue by utilising, in an integrated way, advanced (second- and third-generation) DNA sequencing and informatic approaches to build a high-quality reference genome for C. sinensis, with chromosome-level contiguity and curated gene models. This substantially-enhanced genome provides a resource that could accelerate fundamental and applied molecular investigations of C. sinensis, clonorchiasis and/or cholangiocarcinoma, and assist in the discovery of new interventions against what is a highly significant, but neglected disease-complex.

Mitochondrial genome of Bulinus truncatus (Gastropoda: Lymnaeoidea): Implications for snail systematics and schistosome epidemiology.

Many freshwater snails of the genus Bulinus act as intermediate hosts in the life-cycles of schistosomes in Africa and adjacent regions. Currently, 37 species of Bulinus representing four groups are recognised. The mitochondrial cytochrome c oxidase subunit 1 (cox1) gene has shown utility for identifying and differentiating Bulinus species and groups, but taxonomic relationships based on genetic data are not entirely consistent with those inferred using morphological and biological features. To underpin future systematic studies of members of the genus, we characterised here the mitochondrial genome of Bulinus truncatus (from a defined laboratory strain) using a combined second- and third-generation sequencing and informatics approach, enabling taxonomic comparisons with other planorbid snails for which mitochondrial (mt) genomes were available. Analyses showed consistency in gene order and length among mitochondrial genomes of representative planorbid snails, with the lowest and highest nucleotide diversities being in the cytochrome c oxidase and nicotinamide dehydrogenase subunit genes, respectively. This first mt genome for a representative of the genus Bulinus should provide a useful resource for future investigations of the systematics, population genetics, epidemiology and/or ecology of Bulinus and related snails. The sequencing and informatic workflow employed here should find broad applicability to a range of other snail intermediate hosts of parasitic trematodes.

Marked mitochondrial genetic variation in individuals and populations of the carcinogenic liver fluke Clonorchis sinensis.

Clonorchiasis is a neglected tropical disease caused by the Chinese liver fluke, Clonorchis sinensis, and is often associated with a malignant form of bile duct cancer (cholangiocarcinoma). Although some aspects of the epidemiology of clonorchiasis are understood, little is known about the genetics of C. sinensis populations. Here, we conducted a comprehensive genetic exploration of C. sinensis from endemic geographic regions using complete mitochondrial protein gene sets. Genomic DNA samples from C. sinensis individuals (n = 183) collected from cats and dogs in China (provinces of Guangdong, Guangxi, Hunan, Heilongjiang and Jilin) as well as from rats infected with metacercariae from cyprinid fish from the Russian Far East (Primorsky Krai region) were deep sequenced using the BGISEQ-500 platform. Informatic analyses of mitochondrial protein gene data sets revealed marked genetic variation within C. sinensis; significant variation was identified within and among individual worms from distinct geographical locations. No clear affiliation with a particular location or host species was evident, suggesting a high rate of dispersal of the parasite across endemic regions. The present work provides a foundation for future biological, epidemiological and ecological studies using mitochondrial protein gene data sets, which could aid in elucidating associations between particular C. sinensis genotypes/haplotypes and the pathogenesis or severity of clonorchiasis and its complications (including cholangiocarcinoma) in humans.

First record of a tandem-repeat region within the mitochondrial genome of Clonorchis sinensis using a long-read sequencing approach.

BACKGROUND: Mitochondrial genomes provide useful genetic markers for systematic and population genetic studies of parasitic helminths. Although many such genome sequences have been published and deposited in public databases, there is evidence that some of them are incomplete relating to an inability of conventional techniques to reliably sequence non-coding (repetitive) regions. In the present study, we characterise the complete mitochondrial genome-including the long, non-coding region-of the carcinogenic Chinese liver fluke, Clonorchis sinensis, using long-read sequencing. METHODS: The mitochondrial genome was sequenced from total high molecular-weight genomic DNA isolated from a pool of 100 adult worms of C. sinensis using the MinION sequencing platform (Oxford Nanopore Technologies), and assembled and annotated using an informatic approach. RESULTS: From > 93,500 long-reads, we assembled a 18,304 bp-mitochondrial genome for C. sinensis. Within this genome we identified a novel non-coding region of 4,549 bp containing six tandem-repetitive units of 719-809 bp each. Given that genomic DNA from pooled worms was used for sequencing, some variability in length/sequence in this tandem-repetitive region was detectable, reflecting population variation. CONCLUSIONS: For C. sinensis, we report the complete mitochondrial genome, which includes a long (> 4.5 kb) tandem-repetitive region. The discovery of this non-coding region using a nanopore-sequencing/informatic approach now paves the way to investigating the nature and extent of length/sequence variation in this region within and among individual worms, both within and among C. sinensis populations, and to exploring whether this region has a functional role in the regulation of replication and transcription, akin to the mitochondrial control region in mammals. Although applied to C. sinensis, the technological approach established here should be broadly applicable to characterise complex tandem-repetitive or homo-polymeric regions in the mitochondrial genomes of a wide range of taxa.

Cystic echinococcosis in wild boars (Sus scrofa) from southern Italy: Epidemiological survey and molecular characterization.

Cystic Echinococcosis (CE) caused by Echinococcus granulosus sensu lato (s.l.) is one of the most important parasitic zoonotic diseases in the world and it represents an important public health and socio-economic concern. In the Mediterranean basin, CE is widespread and it is endemic in Italy, with major prevalence in southern areas. Several studies have investigated CE in domestic pigs, however, such data in wild boars are scant. In the last decades the wild boar population in Italy has increased and this ungulate could play an important role in the spreading of CE in the wild. Here we report on the prevalence and fertility rate of hydatid cysts in wild boars that were shot during two hunting seasons (2016-2017) in the Campania region of southern Italy. For each animal, a detailed inspection of the carcass and organs (lungs, liver and spleen) was performed and when cysts were found, their number, morphology and fertility were determined by visual and microscopic examination. Cysts were classified morphologically as fertile, sterile, caseous and calcified. Protoscoleces and germinal layers were collected from individual cysts and DNA was extracted to identify different strains/genotypes of E. granulosus s.l. Out of a total of 2108 wild boars 93 (4.4%) were found positive for CE. Infected animals were 45 males and 48 females, aged between 1 and 8 years. The average number of cysts per wild boar was 1.3 (min 1 - max 13). The total number of cysts collected was 123, of which 118 (95.9%) in the liver, 4 (3.3%) in the lungs and 1 (0.8%) in the spleen. Of all analyzed cysts, 70 (56.9%) were fertile and 53 (43.1%) sterile/acephalous. The presence of fertile cysts in 19.4% of CE-positive animals is noteworthy. Overall, molecular diagnosis showed 19 wild boars infected with the pig strain (G7).

Analysis of nad2 and nad5 enables reliable identification of genotypes G6 and G7 within the species complex Echinococcus granulosus sensu lato.

The larval stages of tapeworms in the species complex Echinococcus granulosus sensu lato cause a zoonotic disease known as cystic echinococcosis (CE). Within this species complex, genotypes G6 and G7 are among the most common genotypes associated with human CE cases worldwide. However, our understanding of ecology, biology and epidemiology of G6 and G7 is still limited. An essential first step towards this goal is correct genotype identification, but distinguishing genotypes G6 and G7 has been challenging. A recent analysis based on complete mitogenome data revealed that the conventional sequencing of the cox1 (366 bp) gene fragment mistakenly classified a subset of G7 samples as G6. On the other hand, sequencing complete mitogenomes is not practical if only genotype or haplogroup identification is needed. Therefore, a simpler and less costly method is required to distinguish genotypes G6 and G7. We compared 93 complete mitogenomes of G6 and G7 from a wide geographical range and demonstrate that a combination of nad2 (714 bp) and nad5 (680 bp) gene fragments would be the best option to distinguish G6 and G7. Moreover, this method allows assignment of G7 samples into haplogroups G7a and G7b. However, due to very high genetic variability of G6 and G7, we suggest to construct a phylogenetic network based on the nad2 and nad5 sequences in order to be absolutely sure in genotype assignment. For this we provide a reference dataset of 93 concatenated nad2 and nad5 sequences (1394 bp in total) containing representatives of G6 and G7 (and haplogroups G7a and G7b), which can be used for the reconstruction of phylogenetic networks.

Long-read sequencing reveals a 4.4 kb tandem repeat region in the mitogenome of Echinococcus granulosus (sensu stricto) genotype G1.

BACKGROUND: Echinococcus tapeworms cause a severe helminthic zoonosis called echinococcosis. The genus comprises various species and genotypes, of which E. granulosus (sensu stricto) represents a significant global public health and socioeconomic burden. Mitochondrial (mt) genomes have provided useful genetic markers to explore the nature and extent of genetic diversity within Echinococcus and have underpinned phylogenetic and population structure analyses of this genus. Our recent work indicated a sequence gap (> 1 kb) in the mt genomes of E. granulosus genotype G1, which could not be determined by PCR-based Sanger sequencing. The aim of the present study was to define the complete mt genome, irrespective of structural complexities, using a long-read sequencing method. METHODS: We extracted high molecular weight genomic DNA from protoscoleces from a single cyst of E. granulosus genotype G1 from a sheep from Australia using a conventional method and sequenced it using PacBio Sequel (long-read) technology, complemented by BGISEQ-500 short-read sequencing. Sequence data obtained were assembled using a recently-developed workflow. RESULTS: We assembled a complete mt genome sequence of 17,675 bp, which is > 4 kb larger than the complete mt genomes known for E. granulosus genotype G1. This assembly includes a previously-elusive tandem repeat region, which is 4417 bp long and consists of ten near-identical 441-445 bp repeat units, each harbouring a 184 bp non-coding region and adjacent regions. We also identified a short non-coding region of 183 bp, which includes an inverted repeat. CONCLUSIONS: We report what we consider to be the first complete mt genome of E. granulosus genotype G1 and characterise all repeat regions in this genome. The numbers, sizes, sequences and functions of tandem repeat regions remain to be studied in different isolates of genotype G1 and in other genotypes and species. The discovery of such 'new' repeat elements in the mt genome of genotype G1 by PacBio sequencing raises a question about the completeness of some published genomes of taeniid cestodes assembled from conventional or short-read sequence datasets. This study shows that long-read sequencing readily overcomes the challenges of assembling repeat elements to achieve improved genomes.

First detection of zoonotic tapeworm Echinococcus granulosus sensu lato genotype G7 in continental Italy.

The larval stage of the species complex Echinococcus granulosus sensu lato (s.l.) is the cause of a widespread zoonotic disease known as cystic echinococcosis (CE). The disease is highly prevalent in southern Italy and represents a serious public health issue. The main aim of this study was to characterize E. granulosus s.l. genotypes from wild boar on a continental area of Italy (Campania region), using recently developed mtDNA markers of nad2 and nad5 for reliable identification of different genotypes. Here, nad5 (680 bp) allowed for a clear identification of G1 and G3, whereas a combination of nad2 (714 bp) and nad5 (1394 bp in total) did the same for genotype G7 and its haplogroups G7a and G7b. The results of this study revealed for the first time the presence of genotype G7 in continental Italy. While haplogroup G7b was previously shown to be restricted to the islands of Corsica and Sardinia, here we demonstrate that haplogroup G7b is also present on the mainland of Italy. This work has implications in designing future strategies to reduce CE in Italy.

Distinguishing Echinococcus granulosus sensu stricto genotypes G1 and G3 with confidence: A practical guide.

Cystic echinococcosis (CE), a zoonotic disease caused by tapeworms of the species complex Echinococcus granulosus sensu lato, represents a substantial global health and economic burden. Within this complex, E. granulosus sensu stricto (genotypes G1 and G3) is the most frequent causative agent of human CE. Currently, there is no fully reliable method for assigning samples to genotypes G1 and G3, as the commonly used mitochondrial cox1 and nad1 genes are not sufficiently consistent for the identification and differentiation of these genotypes. Thus, a new genetic assay is required for the accurate assignment of G1 and G3. Here we use a large dataset of near-complete mtDNA sequences (n = 303) to reveal the extent of genetic variation of G1 and G3 on a broad geographical scale and to identify reliable informative positions for G1 and G3. Based on extensive sampling and sequencing data, we developed a new method, that is simple and cost-effective, to designate samples to genotypes G1 and G3. We found that the nad5 is the best gene in mtDNA to differentiate between G1 and G3, and developed new primers for the analysis. Our results also highlight problems related to the commonly used cox1 and nad1. To guarantee consistent identification of G1 and G3, we suggest using the sequencing of the nad5 gene region (680 bp). This region contains six informative positions within a relatively short fragment of the mtDNA, allowing the differentiation of G1 and G3 with confidence. Our method offers clear advantages over the previous ones, providing a significantly more consistent means to distinguish G1 and G3 than the commonly used cox1 and nad1.

The benefits of analysing complete mitochondrial genomes: Deep insights into the phylogeny and population structure of Echinococcus granulosus sensu lato genotypes G6 and G7.

Cystic echinococcosis (CE) is a zoonotic disease caused by the larval stage of the species complex Echinococcus granulosus sensu lato. Within this complex, genotypes G6 and G7 have been frequently associated with human CE worldwide. Previous studies exploring the genetic variability and phylogeography of genotypes G6 and G7 have been based on relatively short mtDNA sequences, and the resolution of these studies has often been low. Moreover, using short sequences, the distinction between G6 and G7 has in some cases remained challenging. The aim here was to sequence complete mitochondrial genomes (mitogenomes) to obtain deeper insight into the genetic diversity, phylogeny and population structure of genotypes G6 and G7. We sequenced complete mitogenomes of 94 samples collected from 15 different countries worldwide. The results demonstrated that (i) genotypes G6 and G7 can be clearly distinguished when mitogenome sequences are used; (ii) G7 is represented by two major haplogroups, G7a and G7b, the latter being specific to islands of Corsica and Sardinia; (iii) intensive animal trade, but also geographical isolation, have likely had the largest impact on shaping the genetic structure and distribution of genotypes G6 and G7. In addition, we found phylogenetically highly divergent haplotype from Mongolia (Gmon), which had a higher affinity to G6.

Molecular phylogeny based on six nuclear genes suggests that Echinococcus granulosus sensu lato genotypes G6/G7 and G8/G10 can be regarded as two distinct species.

Tapeworms of the species complex of Echinococcus granulosus sensu lato (s. l.) are the cause of a severe zoonotic disease - cystic echinococcosis, which is listed among the most severe parasitic diseases in humans and is prioritized by the World Health Organization. A stable taxonomy of E. granulosus s. l. is essential to the medical and veterinary communities for accurate and effective communication of the role of different species in this complex on human and animal health. E. granulosus s. l. displays high genetic diversity and has been divided into different species and genotypes. Despite several decades of research, the taxonomy of E. granulosus s. l. has remained controversial, especially the species status of genotypes G6-G10. Here the Bayesian phylogeny based on six nuclear loci (7387 bp in total) demonstrated, with very high support, the clustering of G6/G7 and G8/G10 into two separate clades. According to the evolutionary species concept, G6/G7 and G8/G10 can be regarded as two distinct species. Species differentiation can be attributed to the association with distinct host species, largely separate geographical distribution and low level of cross-fertilization. These factors have limited the gene flow between genotypic groups G6/G7 and G8/G10, resulting in the formation of distinct species. We discuss ecological and epidemiological differences that support the validity of these species.

Global phylogeography and genetic diversity of the zoonotic tapeworm Echinococcus granulosus sensu stricto genotype G1.

Echinococcus granulosus sensu stricto (s.s.) is the major cause of human cystic echinococcosis worldwide and is listed among the most severe parasitic diseases of humans. To date, numerous studies have investigated the genetic diversity and population structure of E. granulosus s.s. in various geographic regions. However, there has been no global study. Recently, using mitochondrial DNA, it was shown that E. granulosus s.s. G1 and G3 are distinct genotypes, but a larger dataset is required to confirm the distinction of these genotypes. The objectives of this study were to: (i) investigate the distinction of genotypes G1 and G3 using a large global dataset; and (ii) analyse the genetic diversity and phylogeography of genotype G1 on a global scale using near-complete mitogenome sequences. For this study, 222 globally distributed E. granulosus s.s. samples were used, of which 212 belonged to genotype G1 and 10 to G3. Using a total sequence length of 11,682 bp, we inferred phylogenetic networks for three datasets: E. granulosus s.s. (n = 222), G1 (n = 212) and human G1 samples (n = 41). In addition, the Bayesian phylogenetic and phylogeographic analyses were performed. The latter yielded several strongly supported diffusion routes of genotype G1 originating from Turkey, Tunisia and Argentina. We conclude that: (i) using a considerably larger dataset than employed previously, E. granulosus s.s. G1 and G3 are indeed distinct mitochondrial genotypes; (ii) the genetic diversity of E. granulosus s.s. G1 is high globally, with lower values in South America; and (iii) the complex phylogeographic patterns emerging from the phylogenetic and geographic analyses suggest that the current distribution of genotype G1 has been shaped by intensive animal trade.

Genetic diversity and phylogeography of the elusive, but epidemiologically important Echinococcus granulosus sensu stricto genotype G3.

Cystic echinococcosis (CE) is a severe parasitic disease caused by the species complex Echinococcus granulosus sensu lato. Human infections are most commonly associated with E. granulosus sensu stricto (s.s.), comprising genotypes G1 and G3. The objective of the current study was to provide first insight into the genetic diversity and phylogeography of genotype G3. Despite the epidemiological importance of the genotype, it has remained poorly explored due to the ambiguity in the definition of the genotype. However, it was recently demonstrated that long sequences of mitochondrial DNA (mtDNA) provide a reliable method to discriminate G1 and G3 from each other. Therefore, we sequenced near-complete mtDNA of 39 G3 samples, covering most of the known distribution range and host spectra of the genotype. The phylogenetic network revealed high genetic variation within E. granulosus s.s. G3 and while G3 is significantly less prevalent worldwide than G1, the genetic diversity of both of the genotypes is equally high. We also present the results of the Bayesian phylogeographic analysis, which yielded several well-supported diffusion routes of genotype G3 originating from Turkey and Iran, suggesting the Middle East as the origin of the genotype.

New mitogenome and nuclear evidence on the phylogeny and taxonomy of the highly zoonotic tapeworm Echinococcus granulosus sensu stricto.

Cystic echinococcosis, a zoonotic disease caused by Echinococcus granulosus sensu lato (s. l.), is a significant global public health concern. Echinococcus granulosus s. l. is currently divided into numerous genotypes (G1-G8 and G10) of which G1-G3 are the most frequently implicated genotypes in human infections. Although it has been suggested that G1-G3 could be regarded as a distinct species E. granulosus sensu stricto (s. s.), the evidence to support this is inconclusive. Most importantly, data from nuclear DNA that provide means to investigate the exchange of genetic material between G1-G3 is lacking as none of the published nuclear DNA studies have explicitly included G2 or G3. Moreover, the commonly used relatively short mtDNA sequences, including the complete cox1 gene, have not allowed unequivocal differentiation of genotypes G1-G3. Therefore, significantly longer mtDNA sequences are required to distinguish these genotypes with confidence. The main aim of this study was to evaluate the phylogenetic relations and taxonomy of genotypes G1-G3 using sequences of nearly complete mitogenomes (11,443bp) and three nuclear loci (2984bp). A total of 23 G1-G3 samples were analysed, originating from 5 intermediate host species in 10 countries. The mtDNA data demonstrate that genotypes G1 and G3 are distinct mitochondrial genotypes (separated by 37 mutations), whereas G2 is not a separate genotype or even a monophyletic cluster, but belongs to G3. Nuclear data revealed no genetic separation of G1 and G3, suggesting that these genotypes form a single species due to ongoing gene flow. We conclude that: (a) in the taxonomic sense, genotypes G1 and G3 can be treated as a single species E. granulosus s. s.; (b) genotypes G1 and G3 should be regarded as distinct genotypes only in the context of mitochondrial data; (c) we recommend excluding G2 from the genotype list.