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

HiFi long-read amplicon sequencing for full-spectrum variants of human mtDNA

BACKGROUND

Mitochondrial diseases (MDs) can be caused by single nucleotide variants (SNVs) and structural variants (SVs) in the mitochondrial genome (mtDNA). Presently, identifying deletions in small to medium-sized fragments and accurately detecting low-percentage variants remains challenging due to the limitations of next-generation sequencing (NGS).

METHODS

In this study, we integrated targeted long-range polymerase chain reaction (LR-PCR) and PacBio HiFi sequencing to analyze 34 participants, including 28 patients and 6 controls. Of these, 17 samples were subjected to both targeted LR-PCR and to compare the mtDNA variant detection efficacy.

RESULTS

Among the 28 patients tested by long-read sequencing (LRS), 2 patients were found positive for the m.3243 A > G hotspot variant, and 20 patients exhibited single or multiple deletion variants with a proportion exceeding 4%. Comparison between the results of LRS and NGS revealed that both methods exhibited similar efficacy in detecting SNVs exceeding 5%. However, LRS outperformed NGS in detecting SNVs with a ratio below 5%. As for SVs, LRS identified single or multiple deletions in 13 out of 17 cases, whereas NGS only detected single deletions in 8 cases. Furthermore, deletions identified by LRS were validated by Sanger sequencing and quantified in single muscle fibers using real-time PCR. Notably, LRS also effectively and accurately identified secondary mtDNA deletions in idiopathic inflammatory myopathies (IIMs).

CONCLUSIONS

LRS outperforms NGS in detecting various types of SNVs and SVs in mtDNA, including those with low frequencies. Our research is a significant advancement in medical comprehension and will provide profound insights into genetics.

Evolution of myxozoan mitochondrial genomes: insights from myxobolids

BACKGROUND

Myxozoa is a class of cnidarian parasites that encompasses over 2,400 species. Phylogenetic relationships among myxozoans remain highly debated, owing to both a lack of informative morphological characters and a shortage of molecular markers. Mitochondrial (mt) genomes are a common marker in phylogeny and biogeography. However, only five complete myxozoan mt genomes have been sequenced: four belonging to two closely related genera, Enteromyxum and Kudoa, and one from the genus Myxobolus. Interestingly, while cytochrome oxidase genes could be identified in Enteromyxum and Kudoa, no such genes were found in Myxobolus squamalis, and another member of the Myxobolidae (Henneguya salminicola) was found to have lost its entire mt genome. To evaluate the utility of mt genomes to reconstruct myxozoan relationships and to understand if the loss of cytochrome oxidase genes is a characteristic of myxobolids, we sequenced the mt genome of five myxozoans (Myxobolus wulii, M. honghuensis, M. shantungensis, Thelohanellus kitauei and, Sphaeromyxa zaharoni) using Illumina and Oxford Nanopore platforms.

RESULTS

Unlike Enteromyxum, which possesses a partitioned mt genome, the five mt genomes were encoded on single circular chromosomes. An mt plasmid was found in M. wulii, as described previously in Kudoa iwatai. In all new myxozoan genomes, five protein-coding genes (cob, cox1, cox2, nad1, and nad5) and two rRNAs (rnl and rns) were recognized, but no tRNA. We found that Myxobolus and Thelohanellus species shared unidentified reading frames, supporting the view that these mt open reading frames are functional. Our phylogenetic reconstructions based on the five conserved mt genes agree with previously published trees based on the 18S rRNA gene.

CONCLUSIONS

Our results suggest that the loss of cytochrome oxidase genes is not a characteristic of all myxobolids, the ancestral myxozoan mt genome was likely encoded on a single circular chromosome, and mt plasmids exist in a few lineages. Our findings indicate that myxozoan mt sequences are poor markers for reconstructing myxozoan phylogenetic relationships because of their fast-evolutionary rates and the abundance of repeated elements, which complicates assembly.

The invasive land flatworm Arthurdendyus triangulatus has repeated sequences in the mitogenome, extra-long cox2 gene and paralogous nuclear rRNA clusters

Using a combination of short- and long-reads sequencing, we were able to sequence the complete mitochondrial genome of the invasive 'New Zealand flatworm' Arthurdendyus triangulatus (Geoplanidae, Rhynchodeminae, Caenoplanini) and its two complete paralogous nuclear rRNA gene clusters. The mitogenome has a total length of 20,309 bp and contains repetitions that includes two types of tandem-repeats that could not be solved by short-reads sequencing. We also sequenced for the first time the mitogenomes of four species of Caenoplana (Caenoplanini). A maximum likelihood phylogeny associated A. triangulatus with the other Caenoplanini but Parakontikia ventrolineata and Australopacifica atrata were rejected from the Caenoplanini and associated instead with the Rhynchodemini, with Platydemus manokwari. It was found that the mitogenomes of all species of the subfamily Rhynchodeminae share several unusual structural features, including a very long cox2 gene. This is the first time that the complete paralogous rRNA clusters, which differ in length, sequence and seemingly number of copies, were obtained for a Geoplanidae.

Assembly and phylogeographical analysis of novel Taenia solium mitochondrial genomes suggest stratification within the African-American genotype

BACKGROUND

Taenia solium is a parasite of public health concern, causing human taeniasis and cysticercosis. Two main genotypes have been identified: Asian and African-American. Although characterizing T. solium genotypes is crucial to understanding the genetic epidemiology of its diseases, not much is known about the differences between T. solium mitochondrial genomes from different genotypes. Also, little is known about whether genotypes are further subdivided. Therefore, this study aimed to identify a set of point mutations distributed throughout the T. solium mitochondrial genome that differentiate the African-American from the Asian genotype. Another objective was to identify whether T. solium main genotypes are further stratified.

METHODS

One Mexican and two Peruvian T. solium mitochondrial genomes were assembled using reads available in the NCBI Sequence Read Archive and the reference genome from China as a template. Mutations with respect to the Chinese reference were identified by multiple genome alignment. Jensen-Shannon and Grantham scores were computed for mutations in protein-coding genes to evaluate whether they affected protein function. Phylogenies by Bayesian inference and haplotype networks were constructed using cytochrome c oxidase subunit 1 and cytochrome b from these genomes and other isolates to infer phylogeographical relationships.

RESULTS

A set of 31 novel non-synonymous point mutations present in all genomes of the African-American genotype were identified. These mutations were distributed across the mitochondrial genome, differentiating the African-American from the Asian genotype. All occurred in non-conserved protein positions. Furthermore, the analysis suggested a stratification of the African-American genotypes into an East African and a West African sublineage.

CONCLUSIONS

A novel set of 31 non-synonymous mutations differentiating the main T. solium genotypes was identified. None of these seem to be causing differences in mitochondrial protein function between parasites of the two genotypes. Furthermore, two sublineages within the African-American genotype are proposed for the first time. The presence of the East African sublineage in the Americas suggests an underestimated connection between East African and Latin American countries that might have arisen in the major slave trade between Portuguese Mozambique and the Americas. The results obtained here help to complete the molecular epidemiology of the parasite.

Revealing novel cytb and nad5 genes-based population diversity and benzimidazole resistance in Echinococcus granulosus of bovine origin

Cystic echinococcosis (CE) is a neglected zoonotic disease caused by Echinococcus granulosus (sensu stricto). The parasite affects a wide range of livestock and wild animals. In this study, the population diversity of the Echinococcus species was investigated based on mitochondrial cytochrome b (cytb) and NADH dehydrogenase subunit 5 (nad5) genes. In addition to this, β-tubulin gene isoforms of Echinococcus granulosus were amplified to determine the resistance against benzimidazoles. For this purpose, 40 cyst samples from cattle (n = 20) and buffaloes (n = 20) were collected from the main abattoir of Sialkot. DNA extraction was performed using Qiagen Blood and Tissue Kits. Amplification was performed through PCR. Each amplicon was confirmed by GelRed™ stained agarose gel (2%). Samples were sequenced in a DNA analyzer and viewed for any misread nucleotide by using MEGA (v.11). Corrections in nucleotide sequence and multiple sequence alignment were made through the same software. NCBI-BLAST was used for sample specific sequences to identify them as belonging to a particular species. Diversity indices were estimated using DnaSP (v.6) while phylogenetic analysis was inferred using the Bayesian method using MrBayes (v.1.1). β-tubulin gene isoforms sequence analysis was performed to find out the candidate gene causing benzimidazole resistance. All 40 isolates were found positive for E. granulosus. BLAST-based searches of sequences of each isolate for each gene (nad5 and cytb) confirmed their maximum similarity with the G1 genotype. Overall, high haplotype diversity (Hd nad5 = 1.00; Hd cytb = 0.833) and low nucleotide diversity (π nad5 = 0.00560; π = cytb = 0.00763) was identified based on diversity indices. For both the genes, non-significant values of Tajima's D (nad5 = -0.81734; cytb = -0.80861) and Fu's Fs (nad5 = -1.012; cytb = 0.731) indicate recent population expansion. Bayesian phylogeny-based results of nad5 and cytb sequences confirmed their genotypic status as distinct from other Echinococcus species. This study shed light on the status of benzimidazole resistance in Echinococcus granulosus for the very first time from Pakistan. The findings of this study will significantly add in the information available on genetic diversity of Echinoccous granulosus based on cytb and nad5 genes sequences.

Molecular Phylogeny of Thoracotreme Crabs Including Nine Newly Determined Mitochondrial Genomes

Mitochondrial genomes are used widely for the molecular phylogenetic analysis of animals. Although phylogenetic analyses based on the mitogenomes of brachyurans often yield well-resolved phylogenies, most interfamilial phylogenetic relationships in Thoracotremata remain unclear. We determined nine new mitogenomes of Thoracotremata, including mitogenomes of Camptandriidae (Deiratonotus japonicus), Dotillidae (Ilyoplax integra, Ilyoplax pusilla, and Tmethypocoelis choreutes), Macrophthalmidae (Ilyograpsus nodulosus), Pinnotheridae (Arcotheres sp. and Indopinnixa haematosticta), Plagusiidae (Guinusia dentipes), and Percnidae (Percnon planissimum). Interestingly, Percnon planissimum (Percnidae) was found to possess ≥ 19 repeated sequences in the control region. The gene orders of Il. nodulosus, Arcotheres sp., and In. haematosticta were revealed to be unique among thoracotreme crabs. Although the results of Bayesian and maximum likelihood (ML) phylogenetic analyses of three datasets were incongruent, highly supported clades (PP ≥ 0.99 or BS ≥ 99%) were not contradictory among the analyses. All analyses suggested the paraphyly of Grapsoidea and Ocypodoidea, corroborating the findings of previous studies based on molecular phylogenies of thoracotreme crabs. The phylogenetic positions of symbiotic thoracotreme crabs, Pinnotheridae and Cryptochiridae, were highly supported (Pinnotheridae + Ocypodidae and Cryptochiridae + Grapsidae, respectively) for the Bayesian analyses but not for the ML analyses. Analyses of more thoracotreme species' mitogenome sequences in additional studies will further strengthen the framework for thoracotreme evolution.

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.

Update on the genetic diversity and population structure of Echinococcus granulosus in Gansu Province, Tibet Autonomous Region, and Xinjiang Uygur Autonomous Region, Western China, inferred from mitochondrial cox1, nad1, and nad5 sequences

The identification of additional Echinococcus granulosus sensu lato (s.l.) complex species/genotypes in recent years raises the possibility that there might be more variation among this species in China than is currently understood. The aim of this study was to explore intra- and inter-species variation and population structure of Echinococcus species isolated from sheep in three areas of Western China. Of the isolates, 317, 322, and 326 were successfully amplified and sequenced for cox1, nad1, and nad5 genes, respectively. BLAST analysis revealed that the majority of the isolates were E. granulosus s.s., and using the cox1, nad1, and nad5 genes, respectively, 17, 14, and 11 isolates corresponded to Elodea canadensis (genotype G6/G7). In the three study areas, G1 genotypes were the most prevalent. There were 233 mutation sites along with 129 parsimony informative sites. A transition/transversion ratio of 7.5, 8, and 3.25, respectively, for cox1, nad1, and nad5 genes was obtained. Every mitochondrial gene had intraspecific variations, which were represented in a star-like network with a major haplotype with observable mutations from other distant and minor haplotypes. The Tajima's D value was significantly negative in all populations, indicating a substantial divergence from neutrality and supporting the demographic expansion of E. granulosus s.s. in the study areas. The phylogeny inferred by the maximum likelihood (ML) method using nucleotide sequences of cox1-nad1-nad5 further confirmed their identity. The nodes assigned to the G1, G3, and G6 clades as well as the reference sequences utilized had maximal posterior probability values (1.00). In conclusion, our study confirms the existence of a significant major haplotype of E. granulosus s.s. where G1 is the predominant genotype causing of CE in both livestock and humans in China.

The complete mitogenome of the Asian lung fluke Paragonimus skrjabini miyazakii and its implications for the family Paragonimidae (Trematoda: Platyhelminthes)

The complete circular mitogenome of Paragonimus skrjabini miyazakii (Platyhelminthes: Paragonimidae) from Japan, obtained by PacBio long-read sequencing, was 17 591 bp and contained 12 protein-coding genes (PCGs), 2 mitoribosomal RNA and 22 transfer RNA genes. The atp8 gene was absent, and there was a 40 bp overlap between nad4L and nad4. The long non-coding region (4.3 kb) included distinct types of long and short repeat units. The pattern of base usage for PCGs and the mtDNA coding region overall in Asian and American Paragonimus species (P. s. miyazakii, P. heterotremus, P. ohirai and P. kellicotti) and the Indian form of P. westermani was T > G > A > C. On the other hand, East-Asian P. westermani used T > G > C > A. Five Asian and American Paragonimus species and P. westermani had TTT/Phe, TTG/Leu and GTT/Val as the most frequently used codons, whereas the least-used codons were different in each species and between regional forms of P. westermani. The phylogenetic tree reconstructed from a concatenated alignment of amino acids of 12 PCGs from 36 strains/26 species/5 families of trematodes confirmed that the Paragonimidae is monophyletic, with 100% nodal support. Paragonimus skrjabini miyazakii was resolved as a sister to P. heterotremus. The P. westermani clade was clearly separate from remaining congeners. The latter clade was comprised of 2 subclades, one of the East-Asian and the other of the Indian Type 1 samples. Additional mitogenomes in the Paragonimidae are needed for genomic characterization and are useful for diagnostics, identification and genetic/ phylogenetic/ epidemiological/ evolutionary studies of the Paragonimidae.

Phylogenetic Analysis of Mitochondrial Genome of Tabanidae (Diptera: Tabanidae) Reveals the Present Status of Tabanidae Classification

Simple Summary

Tabanidae suck the blood of humans and animals, are important biological vectors for the transmission of diseases, and are of considerable economic and medical significance. However, current knowledge about the mitochondrial genome of this family is limited. Therefore, six newly completed mitochondrial genomes of four genera of Tabanidae (Haematopota turkestanica, Chrysops vanderwulpi, Chrysops dissectus, Tabanus chrysurus, Tabanus pleskei, and Hybomitra sp. species) were sequenced and analyzed. The results show that the six newly mitochondrial genomes have quite similar structures and features. Phylogeny was inferred by analyzing the 13 amino acid sequences coded by mitochondrial genes of 22 mitogenomes (all available complete mitochondrial genomes of tabanidae). Bayesian inference, maximum likelihood trees, and maximum parsimony inference analyses all showed consistent results. This study supports the concept of monophyly of all groups, ratifies the current taxonomic classification, and provides useful genetic markers for studying the molecular ecology, systematics, and population genetics of Tabanidae.

Abstract

Tabanidae suck the blood of humans and animals, are important biological vectors for the transmission of diseases, and are of considerable economic and medical significance. However, current knowledge about the mitochondrial genome of this family is limited. More complete mitochondrial genomes of Tabanidae are essential for the identification and phylogeny. Therefore, this study sequenced and analyzed six complete mitochondrial (mt) genome sequences of four genera of Tabanidae for the first time. The complete mt genomes of the six new sequences are circular molecules ranging from 15,851 to 16,107 base pairs (bp) in size, with AT content ranging from 75.64 to 77.91%. The six complete mitochondrial genomes all consist of 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (RRNA), 22 transfer RNA genes (tRNAs), and a control region, making a total of 37 functional subunits. ATT/ATG was the most common start codon, and the stop codon was TAA of all PCGS. All tRNA except tRNA Ser1 had a typical clover structure. Phylogeny was inferred by analyzing the 13 concatenated amino acid sequences of the 22 mt genomes. Bayesian inference, maximum-likelihood trees, and maximum-parsimony inference analyses all showed consistent results. This study supports the concept of monophyly of all genus, ratifies the current taxonomic classification, and provides effective genetic markers for molecular classification, systematics, and genetic studies of Tabanidae.

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.

Echinococcus spp. and genotypes infecting humans in Tibet Autonomous Region of China: a molecular investigation with near-complete/complete mitochondrial sequences

Background

Molecular markers are essential to identify Echinococcus species and genotypes in areas with multiple Echinococcus species to understand their epidemiology and pathology. Tibet Autonomous Region (TAR) is one of the areas worst hit by echinococcosis. However, molecular epidemiology is still missing among echinococcosis patients in TAR. This research explored the Echinococcus species and genotypes infecting humans in TAR and the population diversity and the possible origin of G1 in TAR.

Methods

Cyst samples were collected in one echinococcosis-designated hospital in TAR. Echinococcus species and genotypes were identified through a maximum-likelihood approach with near-complete/complete mtDNA using IQ-TREE. Phylogenetic networks were built with PopART, and the phylogeographical diffusion pattern was identified using a Bayesian discrete phylogeographic method.

Results

Using phylogenetic trees made with near-complete/complete mtDNA obtained from 92 cysts from TAR patients, the Echinococcus species and genotypes infecting humans in TAR were identified as Echinococcus granulosus (s.s.) G1 (81, 88.04%), accounting for the majority, followed by G6 of the E. canadensis cluster (6, 6.52%), E. granulosus (s.s.) G3 (3, 3.26%), and E. multilocularis (2, 2.17%). An expansion trend and a possible recent bottleneck event were confirmed among the G1 samples in TAR. Adding the other near-complete mtDNA of G1 samples globally from the literature, we identified the possible phylogeographic origin of the G1 samples in TAR as Turkey.

Conclusions

Using near-complete/complete mtDNA sequences of Echinococcus spp. obtained from echinococcosis patients, a variety of Echinococcus species and genotypes infecting humans throughout TAR were identified. As far as we know, this is the first comprehensive molecular investigation of Echinococcus species and genotypes infecting humans throughout TAR. We identified, for the first time to our knowledge, the possible origin of the G1 in TAR. We also enriched the long mtDNA database of Echinococcus spp. and added two complete E. multilocularis mtDNA sequences from human patients. These findings will improve our knowledge of echinococcosis, help to refine the targeted echinococcosis control measures, and serve as a valuable baseline for monitoring the Echinococcus species and genotypes mutations and trends of the Echinococcus spp. population in TAR.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05199-6.

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.

A high-quality genome for the parasitic tapeworm, Echinococcus granulosus, provides further insight into the EG95 vaccine target for cystic echinococcosis.

Expansion of Cyclophyllidea Biodiversity in Rodents of Qinghai-Tibet Plateau and the "Out of Qinghai-Tibet Plateau" Hypothesis of Cyclophyllideans

The Cyclophyllidea comprises the most species-rich order of tapeworms (Platyhelminthes, Cestoda) and includes species with some of the most severe health impact on wildlife, livestock, and humans. We collected seven Cyclophyllidea specimens from rodents in Qinghai-Tibet Plateau (QTP) and its surrounding mountain systems, of which four specimens in QTP were unsequenced, representing “putative new species.” Their complete mitochondrial (mt) genomes were sequenced and annotated. Phylogenetic reconstruction of partial 28S rDNA, cox1 and nad1 datasets provided high bootstrap frequency support for the categorization of three “putative new species,” assigning each, respectively, to the genera Mesocestoides, Paranoplocephala, and Mosgovoyia, and revealing that some species and families in these three datasets, which contain 291 species from nine families, may require taxonomic revision. The partial 18S rDNA phylogeny of 29 species from Taeniidae provided high bootstrap frequency support for the categorization of the “putative new species” in the genus Hydatigera. Combined with the current investigation, the other three known Taeniidae species found in this study were Taenia caixuepengi, T. crassiceps, and Versteria mustelae and may be widely distributed in western China. Estimates of divergence time based on cox1 + nad1 fragment and mt protein-coding genes (PCGs) showed that the differentiation rate of Cyclophyllidea species was strongly associated with the rate of change in the biogeographic scenarios, likely caused by the uplift of the QTP; i.e., species differentiation of Cyclophyllidea might be driven by host-parasite co-evolution caused by the uplift of QTP. We propose an “out of QTP” hypothesis for the radiation of these cyclophyllidean tapeworms.

Mitophylogenomics of the zoonotic fluke Echinostoma malayanum confirms it as a member of the genus Artyfechinostomum Lane, 1915 and illustrates the complexity of Echinostomatidae systematics

The complete mitochondrial genome (mitogenome or mtDNA) of the trematode Echinostoma malayanum Leiper, 1911 was fully determined and annotated. The circular mtDNA molecule comprised 12 protein-coding genes (PCGs) (cox1 - 3, cob, nad1 - 6, nad4L, atp6), two mitoribosomal RNAs (MRGs) (16S or rrnL and 12S or rrnS), and 22 transfer RNAs (tRNAs or trn), and a non-coding region (NCR) rich in long and short tandem repeats (5.5 LRUs/336 bp/each and 7.5 SRUs/207 bp/each). The atp8 gene is absent and the 3' end of nad4L overlaps the 5' end of nad4 by 40 bp. Special DHU-arm missing tRNAs for Serine were found for both tRNA^Ser1(AGN) and tRNA^Ser2(UCN). Codons of TTT (for phenylalanine), TTG (for leucine), and GTT (for valine) were the most, and CGC (for Arginine) was the least frequently used. A similar usage pattern was seen in base composition, AT and GC skewness for PCGs, MRGs, and mtDNA* (coding cox3 to nad5) in E. malayanum and Echinostomatidae. The nucleotide use is characterized by (T > G > A > C) for PCGs/mtDNA*, and by (T > G ≈ A > C) for MRGs. E. malayanum exhibited the lowest genetic distance (0.53%) to Artyfechinostomum sufrartyfex, relatively high to the Echinostoma congeners (13.20-13.99%), higher to Hypoderaeum conoideum (16.18%), and the highest to interfamilial Echinochasmidae (26.62%); Cyclocoelidae (30.24%); and Himasthlidae (25.36%). Topology indicated the monophyletic position between E. malayanum/A. sufrartyfex and the group of Echinostoma caproni, Echinostoma paraensei, Echinostoma miyagawai, and Echinostoma revolutum, rendering Hypoderaeum conoideum and unidentified Echinostoma species paraphyletic. The strictly closed genomic/taxonomic/phylogenetic features (including base composition, skewness, codon usage/bias, genetic distance, and topo-position) reinforced Echinostoma malayanum to retake its generic validity within the Artyfechinostomum genus.

Hammerhead flatworms (Platyhelminthes, Geoplanidae, Bipaliinae): mitochondrial genomes and description of two new species from France, Italy, and Mayotte

BACKGROUND

New records of alien land planarians are regularly reported worldwide, and some correspond to undescribed species of unknown geographic origin. The description of new species of land planarians (Geoplanidae) should classically be based on both external morphology and histology of anatomical structures, especially the copulatory organs, ideally with the addition of molecular data.

METHODS

Here, we describe the morphology and reproductive anatomy of a species previously reported as Diversibipalium "black", and the morphology of a species previously reported as Diversibipalium "blue". Based on next generation sequencing, we obtained the complete mitogenome of five species of Bipaliinae, including these two species.

RESULTS

The new species Humbertium covidum n. sp. (syn: Diversibipalium "black" of Justine et al., 2018) is formally described on the basis of morphology, histology and mitogenome, and is assigned to Humbertium on the basis of its reproductive anatomy. The type-locality is Casier, Italy, and other localities are in the Department of Pyrénées-Atlantiques, France; some published or unpublished records suggest that this species might also be present in Russia, China, and Japan. The mitogenomic polymorphism of two geographically distinct specimens (Italy vs France) is described; the cox1 gene displayed 2.25% difference. The new species Diversibipalium mayottensis n. sp. (syn: Diversibipalium "blue" of Justine et al., 2018) is formally described on the basis of external morphology and complete mitogenome and is assigned to Diversibipalium on the basis of an absence of information on its reproductive anatomy. The type- and only known locality is the island of Mayotte in the Mozambique Channel off Africa. Phylogenies of bipaliine geoplanids were constructed on the basis of SSU, LSU, mitochondrial proteins and concatenated sequences of cox1, SSU and LSU. In all four phylogenies, D. mayottensis was the sister-group to all the other bipaliines. With the exception of D. multilineatum which could not be circularised, the complete mitogenomes of B. kewense, B. vagum, B. adventitium, H. covidum and D. mayottensis were colinear. The 16S gene in all bipaliine species was problematic because usual tools were unable to locate its exact position.

CONCLUSION

Next generation sequencing, which can provide complete mitochondrial genomes as well as traditionally used genes such as SSU, LSU and cox1, is a powerful tool for delineating and describing species of Bipaliinae when the reproductive structure cannot be studied, which is sometimes the case of asexually reproducing invasive species. The unexpected position of the new species D. mayottensis as sister-group to all other Bipaliinae in all phylogenetic analyses suggests that the species could belong to a new genus, yet to be described.

Advances in the treatment, diagnosis, control and scientific understanding of taeniid cestode parasite infections over the past 50 years

In the past 50 years, enormous progress has been made in the diagnosis, treatment and control of taeniid cestode infections/diseases and in the scientific understanding thereof. Most interest in this group of parasites stems from the serious diseases that they cause in humans. It is through this lens that we summarize here the most important breakthroughs that have made a difference to the treatment of human diseases caused by these parasites, reduction in transmission of the taeniid species associated with human disease, or understanding of the parasites' biology likely to impact diagnosis or treatment in the foreseeable future. Key topics discussed are the introduction of anti-cestode drugs, including benzimidazoles and praziquantel, and the development of new imaging modalities that have transformed the diagnosis and post-treatment monitoring of human echinococcoses and neurocysticercosis. The availability of new anti-cestode drugs for use in dogs and a detailed understanding of the transmission dynamics of Echinococcus granulosus sensu lato have underpinned successful programs that have eliminated cystic echinococcosis in some areas of the world and greatly reduced the incidence of infection in others. Despite these successes, cystic and alveolar echinococcosis and neurocysticercosis continue to be prevalent in many parts of the world, requiring new or renewed efforts to prevent the associated taeniid infections. Major advances made in the development of practical vaccines against E. granulosus and Taenia solium will hopefully assist in this endeavour, as might the understanding of the parasites' biology that have come from an elucidation of the nuclear genomes of each of the most important taeniid species causing human diseases.

Cestodes in the genomic era

The first cestode genomes were obtained by an international consortium led by the Wellcome Sanger Institute that included representative institutions from countries where the sequenced parasites have been studied for decades, in part because they are etiological agents of endemic diseases (Argentina, Uruguay, Mexico, Canada, UK, Germany, Switzerland, Ireland, USA, Japan, and China). After this, several complete genomes were obtained reaching 16 species to date. Cestode genomes have smaller relative size compared to other animals including free-living flatworms. Moreover, the features genome size and repeat content seem to differ in the two analyzed orders. Cyclophyllidean species have smaller genomes and with fewer repetitive content than Diphyllobothriidean species. On average, cestode genomes have 13,753 genes with 6 exons per gene and 41% GC content. More than 5,000 shared cestode proteins were accurately annotated by the integration of gene predictions and transcriptome evidence being more than 40% of these proteins of unknown function. Several gene losses and reduction of gene families were found and could be related to the extreme parasitic lifestyle of these species. The application of cutting-edge sequencing technology allowed the characterization of the terminal sequences of chromosomes that possess unique characteristics. Here, we review the current status of knowledge of complete cestode genomes and place it within a comparative genomics perspective. Multidisciplinary work together with the implementation of new technologies will provide valuable information that can certainly improve our chances to finally eradicate or at least control diseases caused by cestodes.

Characterization of the complete mitochondrial genomes of six horseflies (Diptera: Tabanidae)

The family Tabanidae (Insecta: Diptera) is one of the economically most important group of haematophagous insects, causing millions of livestock deaths per year. However, current knowledge on the mitochondrial genomes (mitogenomes) from this family is limited. Additional tabanid mitogenomes characterization is of utmost importance for their identification, epidemiologic and phylogenetic studies. We sequenced the mt genomes of six horseflies with an Illumina platform and their phylogenetic relationship was conducted with other infraorder Tabanomorpha members with available mt genome datasets. All six newly sequenced mitogenomes were typical 37-gene circular structures retaining the gene order of Tabanomorpha. The trnQ, trnM and trnA were highly conserved among the six mitogenomes (identity = 100%). The TΨC arm and variable loop regions were relatively more variable compared to the amino acid receptor arm, anticodon arm and DHU arm of the tRNAs. Among 13 protein-coding genes (PCGs) of tabanids mitogenomes, the highest nucleotide diversity was detected in atp8, cox1, cox3, nad6 and cytb (0.1 for each). In addition, atp8 genes exhibited the highest evolutionary rate (ω = 0.24) among 13 PCGs. The interspecies K2P genetic distances among some Tabanus spp. across the mitogenome was greater (0.08) than intergeneric genetic distance between T. amaenus and Atylotus miser (0.07). Phylogenetic analyses revealed non-monophyletic relationships among horseflies of the genus Tabanus. The present study showed mt gene order is highly conserved within Tabanus species. Our mito-phylogenomic analysis supports the paraphyly of the genus Tabanus. The new data provide novel genetic markers for studies of population genetics and systematics of horseflies.

Genetic diversity of Echinococcus granulosus sensu stricto infecting humans in western Algeria

Human cystic echinococcosis is a zoonosis due to the flat worm Echinococcus granulosus sensu lato. The disease remains a major public health problem in Northern Africa. Molecular typing enables a better understanding of the parasite circulation from animals to humans. In this study, we investigated the genotypic diversity of 46 Echinococcus granulosus isolates collected from humans in the western part of Algeria by the mean of partial sequences of 4 mitochondrial loci, namely cox1a, cox1b, nd3, and atp6. Nucleotide polymorphism ranges from 0.6 (nd3) to 2.7% (cox1a). Eight alleles had not been previously reported. Multilocus analysis showed that all the isolates were from the Echinococcus granulosus sensu stricto (G1 genotype). Nineteen different haplotypes made of the concatenation of 4 sequenced loci were observed, the most common type clustering 13 isolates (36.1%). Twelve of these haplotypes had never been described previously and fifteen (41.7%) haplotypes were represented by only one isolate. Using sequences from this study and others retrieved from the GenBank database, any clustering either according to the geographic origin within Algeria or according to the human or animal origin of the isolates could be demonstrated supporting that genotype G1 population genetics has been shaped by intensive animal breeding.

The complete mitochondrial genome sequences of the cat flea Ctenocephalides felis felis (Siphonaptera: Pulicidae) support the hypothesis that C. felis isolates from China and USA were the same C. f. felis subspecies

The cat flea Ctenocephalides felis (Siphonaptera: Pulicidae) is the most important ectoparasite in cats and dogs worldwide. Over the years, there has been much dispute regarding the taxonomic and systematic status of C. felis. Mitochondrial (mt) genome sequences are useful genetic markers for the identification and differentiation of ectoparasites, but the mt genome of C. felis and its subspecies has not yet been entirely characterized. In the present study, the entire mt genome of C. f. felis from China was sequenced and compared with that of C. felis from the USA. Both contain 37 genes and a long non-coding region of >6 kbp. The molecular identity between the Chinese and American isolates was 99%, except for the non-coding region. The protein-coding genes showed differences at both the nucleotide (1.2%) and amino acid (1%) levels. Interestingly, the cox1 gene of the Chinese isolate had an unusual putative start codon (TTT). Taken together, our analyses strongly support the hypothesis that C. felis isolates from China and the USA were the same C. f. felis subspecies. The mt genome sequence of the C. f. felis China isolate presented in this study provides useful molecular markers to further address the taxonomy and systematics of C. felis.

Three species of Echinococcus granulosus sensu lato infect camels on the Arabian Peninsula

We report on the genetic identity of 36 Echinococcus cysts that were collected during a recent slaughterhouse survey of 810 locally bred camels (dromedaries) in the Eastern Province of the Kingdom of Saudi Arabia. Analysis of a partial nad1 gene sequence showed that the majority (n = 29) belonged to E. granulosus sensu stricto, four to E. canadensis G6/7, and three to E. ortleppi. Eight of the 29 E. granulosus s.s. cysts contained protoscoleces; all other cysts were calcified and non-viable. This is the first report of the presence E. ortleppi from the Arabian Peninsula, a parasite that is typically transmitted via cattle. The results indicate widespread infection of camels with CE in eastern Saudi Arabia and an active role of camels in the lifecycles of at least E. granulosus s.s.. Complete cox1 haplotype analysis of 21 E. granulosus s.s. isolates shows that the majority of variants circulating in eastern Saudi Arabia is distinct from but closely related to haplotypes from neighboring countries in the Middle East, which indicates the presence of this parasite in KSA for a longer period of time. All isolates of E. granulosus s.s. in this study belonged to the G1 cluster, although the G3 genotype has previously also been reported from the Middle East.

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.

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.

Mitochondrial genomes of two eucotylids as the first representatives from the superfamily Microphalloidea (Trematoda) and phylogenetic implications

Background

The Eucotylidae Cohn, 1904 (Superfamily: Microphalloidea), is a family of digeneans parasitic in kidneys of birds as adults. The group is characterized by the high level of morphological similarities among genera and unclear systematic value of morphological characters traditionally used for their differentiation. In the present study, we sequenced the complete or nearly complete mitogenomes (mt genome) of two eucotylids representing the genera Tamerlania (T. zarudnyi) and Tanaisia (Tanaisia sp.). They represent the first sequenced mt genomes of any member of the superfamily Microphalloidea.

Methods

A comparative mitogenomic analysis of the two newly sequenced eucotylids was conducted for the investigation of mitochondrial gene arrangement, contents and genetic distance. Phylogenetic position of the family Eucotylidae within the order Plagiorchiida was examined using nucleotide sequences of mitochondrial protein-coding genes (PCGs) plus RNAs using maximum likelihood (ML) and Bayesian inference (BI) methods. BI phylogeny based on concatenated amino acids sequences of PCGs was also conducted to determine possible effects of silent mutations.

Results

The complete mt genome of T. zarudnyi was 16,188 bp and the nearly complete mt genome of Tanaisia sp. was 13,953 bp in length. A long string of additional amino acids (about 123 aa) at the 5′ end of the cox1 gene in both studied eucotylid mt genomes has resulted in the cox1 gene of eucotylids being longer than in all previously sequenced digeneans. The rrnL gene was also longer than previously reported in any digenean mitogenome sequenced so far. The TΨC and DHU loops of the tRNAs varied greatly between the two eucotylids while the anticodon loop was highly conserved. Phylogenetic analyses based on mtDNA nucleotide and amino acids sequences (as a separate set) positioned eucotylids as a sister group to all remaining members of the order Plagiorchiida. Both ML and BI phylogenies revealed the paraphyletic nature of the superfamily Gorgoderoidea and the suborder Xiphidiata.

Conclusions

The average sequence identity, combined nucleotide diversity and Kimura-2 parameter distances between the two eucotylid mitogenomes demonstrated that atp6, nad5, nad4L and nad6 genes are better markers than the traditionally used cox1 or nad1 for the species differentiation and population-level studies of eucotylids because of their higher variability. The position of the Dicrocoeliidae and Eucotylidae outside the clade uniting other xiphidiatan trematodes strengthened the argument for the need for re-evaluation of the taxonomic content of the Xiphidiata.

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.

In the present study, we characterised the complete mitochondrial genome of Clonorchis sinensis—a carcinogenic liver fluke. To do this, we sequenced from total genomic DNA from multiple adult worms using a new method (Oxford Nanopore technology) to obtain data for long stretches of DNA, and then assembled these data to construct a mitochondrial genome of 18,304 bp, containing a > 4.5 kb-long tandem-repetitive region—not previously detected in this species. The results demonstrate that this method is effective at sequencing long and complex non-coding elements—not achievable using conventional techniques. The discovery of this long tandem-repetitive region in C. sinensis provides an opportunity to now explore its origin(s) and length/sequence diversity in populations of this species, and also to characterise its function(s). The technological approach employed here should have broad applicability to characterise previously-elusive non-coding mitochondrial genomic regions in a wide range of taxa.

Long-term (35 years) cryopreservation of Echinococcus multilocularis metacestodes

The metacestode of Echinococcus multilocularis is the etiological agent of alveolar echinococcosis. The metacestode stage used for research is maintained in rodents by serial passages. In order to determine whether cryopreservation of E. multilocularis metacestodes would be suitable for long-term maintenance and replace serial passages, isolates of different geographic origin were cryopreserved in 1984–1986. The aim of the current study was to test the viability of cryopreserved isolates following long-term cryopreservation (up to 35 years) and to determine the phylogenetic clades these isolates belonged to. Cryopreserved isolates were tested for viability in vitro and in vivo in gerbils. In vitro results of 5 isolates indicated protoscolex survival in 13 of 17 experiments (76%) and metacestode survival in 5 of 12 (42%) in vivo experiments. In vivo results showed ‘abortive lesions’ in 13 of the 36 animals, 15 were negative and 8 harboured proliferating metacestode tissue containing protoscoleces. Genetic analysis confirmed the isolates belonged to European, Asian and North-American clades. In conclusion, the results of the current study indicate that metacestodes of E. multilocularis are able to survive long-term cryopreservation. Therefore, cryopreservation is a suitable method for long-term storage of E. multilocularis metacestode isolates and reduces the number of experimental animals.

Comprehensive characterization of plasma cell-free Echinococcus spp. DNA in echinococcosis patients using ultra-high-throughput sequencing

Background

Echinococcosis is a life-threatening parasitic disease caused by Echinococcus spp. tapeworms with over one million people affected globally at any time. The Echinococcus spp. tapeworms in the human body release DNA to the circulatory system, which can be a biomarker for echinococcosis. Cell-free DNA (cfDNA) is widely used in medical research and has been applied in various clinical settings. As for echinococcosis, several PCR-based tests had been trialed to detect cell-free Echinococcus spp. DNA in plasma or serum, but the sensitivity was about 20% to 25%. Low sensitivity of PCR-based methods might be related to our limited understanding of the features of cell-free Echinococcus spp. DNA in plasma, including its concentration, fragment pattern and release source. In this study, we applied ultra-high-throughput sequencing to comprehensively investigate the characteristics of cell-free Echinococcus spp. DNA in plasma of echinococcosis patients.

Methodology/Principal findings

We collected plasma samples from 23 echinococcosis patients. Total plasma cfDNA was extracted and sequenced with a high-throughput sequencing platform. An average of 282 million read pairs were obtained for each plasma sample. Sequencing data were analyzed with bioinformatics workflow combined with Echinococcus spp. sequence database. After identification of cell-free Echinococcus spp. reads, we found that the cell-free Echinococcus spp. reads accounted for 1.8e-5 to 4.0e-9 of the total clean reads. Comparing fragment length distribution of cfDNA between Echinococcus spp. and humans showed that cell-free Echinococcus spp. DNA of cystic echinococcosis (CE) had a broad length range, while that of alveolar echinococcosis (AE) had an obvious peak at about 135 bp. We found that most of the cell-free Echinococcus spp. DNA reads were from the nuclear genome with an even distribution, which might indicate a random release pattern of cell-free Echinococcus spp. DNA.

Conclusions/Significance

With ultra-high-throughput sequencing technology, we analyzed the concentration, fragment length, release source, and other characteristics of cell-free Echinococcus spp. DNA in the plasma of echinococcosis patients. A better understanding of the characteristics of cell-free Echinococcus spp. DNA in plasma may facilitate their future application as a biomarker for diagnosis.

Echinococcosis is one of the most neglected tropical diseases caused by the metacestodes of Echinococcus spp. tapeworms, which affect both humans and livestock. Plasma cell-free DNA (cfDNA) consists of nucleic acid fragments found extracellularly and may contain DNA released from the parasites. Research shows that a variety of parasites can be detected from plasma cfDNA. Cell-free Echinococcus spp. DNA in plasma or serum had been tested with PCR-based methods, but these PCR methods had low sensitivity ranged from 20% to 25%. Low sensitivity may be due to our limited understanding of cell-free Echinococcus spp. DNA in plasma. Here, we take advantage of high-throughput sequencing to get a comprehensive characterization of cell-free Echinococcus spp. DNA. Our results showed that with high-throughput sequencing we could detect cell-free Echinococcus spp. DNA in all samples, though at a very low level. Based on the sequencing data, we found that cell-free Echinococcus spp. DNA in plasma had a different fragment length distribution to cell-free human DNA, and fragment length distribution of cell-free Echinococcus spp. DNA is also different between cystic echinococcosis (CE) and alveolar echinococcosis (AE). The sequencing data can also help trace the release source of cell-free Echinococcus spp. DNA from the genome. According to the mapping results of cell-free Echinococcus spp. DNA reads, we found that most of them were from the nuclear genome rather than the mitochondrial genome, and their release position showed an even distribution on the genome. These characteristics of cell-free Echinococcus spp. DNA in echinococcosis patients’ plasma could facilitate their future application in research or clinical settings.

Molecular phylogenetics and mitogenomics of three avian dicrocoeliids (Digenea: Dicrocoeliidae) and comparison with mammalian dicrocoeliids

Background

The Dicrocoeliidae are digenetic trematodes mostly parasitic in the bile ducts and gall bladder of various avian and mammalian hosts. Until recently their systematics was based on morphological data only. Due to the high morphological uniformity across multiple dicrocoeliid taxa and insufficient knowledge of relative systematic value of traditionally used morphological characters, their taxonomy has always been unstable. Therefore, DNA sequence data provide a critical independent source of characters for phylogenetic inference and improvement of the system.

Methods

We examined the phylogenetic affinities of three avian dicrocoeliids representing the genera Brachylecithum, Brachydistomum and Lyperosomum, using partial sequences of the nuclear large ribosomal subunit (28S) RNA gene. We also sequenced the complete or nearly complete mitogenomes of these three isolates and conducted a comparative mitogenomic analysis with the previously available mitogenomes from three mammalian dicrocoeliids (from 2 different genera) and examined the phylogenetic position of the family Dicrocoeliidae within the order Plagiorchiida based on concatenated nucleotide sequences of all mitochondrial genes (except trnG and trnE).

Results

Combined nucleotide diversity, Kimura-2-parameter distance, non-synonymous/synonymous substitutions ratio and average sequence identity analyses consistently demonstrated that cox1, cytb, nad1 and two rRNAs were the most conserved and atp6, nad5, nad3 and nad2 were the most variable genes across dicrocoeliid mitogenomes. Phylogenetic analyses based on mtDNA sequences did not support the close relatedness of the Paragonimidae and Dicrocoeliidae and suggested non-monophyly of the Gorgoderoidea as currently recognized.

Conclusions

Our results show that fast-evolving mitochondrial genes atp6, nad5 and nad3 would be better markers than slow-evolving genes cox1 and nad1 for species discrimination and population level studies in the Dicrocoeliidae. Furthermore, the Dicrocoeliidae being outside of the clade containing other xiphidiatan trematodes suggests a need for the re-evaluation of the taxonomic content of the Xiphidiata.

Comparative mitogenomics of the zoonotic parasite Echinostoma revolutum resolves taxonomic relationships within the 'E. revolutum' species group and the Echinostomata (Platyhelminthes: Digenea)

The complete mitochondrial sequence of 17,030 bp was obtained from Echinostoma revolutum and characterized with those of previously reported members of the superfamily Echinostomatoidea, i.e. six echinostomatids, one echinochasmid, five fasciolids, one himasthlid, and two cyclocoelids. Relationship within suborders and between superfamilies, such as Echinostomata, Pronocephalata, Troglotremata, Opisthorchiata, and Xiphiditata, are also considered. It contained 12 protein-coding, two ribosomal RNA, 22 transfer RNA genes and a tandem repetitive consisting non-coding region (NCR). The gene order, one way-positive transcription, the absence of atp8 and the overlapped region by 40 bp between nad4L and nad4 genes were similar as in common trematodes. The NCR located between tRNAGlu (trnE) and cox3 contained 11 long (LRUs) and short repeat units (SRUs) (seven LRUs of 317 bp, four SRUs of 207 bp each), and an internal spacer sequence between LRU7 and SRU4 specifying high-level polymorphism. Special DHU-arm missing tRNAs for Serine were found for both tRNAS1(AGN) and tRNAS2(UCN). Echinostoma revolutum indicated the lowest divergence rate to E. miyagawai and the highest to Tracheophilus cymbius and Echinochasmus japonicus. The usage of ATG/GTG start and TAG/TAA stop codons, the AT composition bias, the negative AT-skewness, and the most for Phe/Leu/Val and the least for Arg/Asn/Asp codons were noted. Topology indicated the monophyletic position of E. revolutum to E. miyagawai. Monophyly of Echinostomatidae and Fasciolidae was clearly solved with respect to Echinochasmidae, Himasthlidae, and Cyclocoelidae which were rendered paraphyletic in the suborder Echinostomata.

Genetic variation of Echinococcus spp. in yaks and sheep in the Tibet Autonomous Region of China based on mitochondrial DNA

Background

Cystic echinococcosis (CE) in humans and livestock is caused by Echinococcus granulosus (sensu lato). In China where CE is endemic, a number of studies have shown that Echinococcus granulosus (sensu stricto) is majorly responsible for CE. However, E. canadensis (G6) which is the second leading cause of CE is now being detected in most parts of the country. In this study, the species diversity and genetic variation of Echinococcus granulosus (s.l.) in four counties in Tibet Autonomous Region of China were investigated.

Methods

Infection with Echinococcus granulosus (s.s.) in yaks and sheep was identified using NADH dehydrogenase subunit 1 and 5 (nad1 and nad5) mitochondrial genes while the genotype G6 of E. canadensis initially diagnosed with NADH dehydrogenase subunit 1 (nad1) was further confirmed by analysis of the complete mitochondrial genome and a phylogenetic network constructed based on the nad2 and nad5 genes.

Results

Out of 85 hydatid cyst samples collected from slaughtered sheep (n = 54) and yaks (n = 31), 83 were identified as E. granulosus (s.s.) G1 (n = 77), G3 (n = 6) and 2 were identified as E. canadensis G6. Analysis of the nad1/nad5 genes revealed 16/17 mutations with 9/14 parsimony informative sites resulting in 15/14 haplotypes, respectively. Haplotype diversity (Hd) and nucleotide diversity (π) of E. granulosus (s.s.) population were 0.650 and 0.00127 for nad1 and 0.782 and 0.00306 for nad5, respectively, with an overall negative Tajima’s D and Fu’s Fs. A low F_ST indicated no genetic difference between isolates from sheep and yaks.

Conclusion

Pockets of infection with E. canadensis (G6, G7, G8 and G10) have been previously reported in sheep, goats, yaks and/or humans in different parts of China. While the G6 genotype has been previously reported in sheep in the Tibet Autonomous Region, the detection in a yak in the present study represents the first to the best of our knowledge. Therefore, we recommend future surveys and control efforts to comprehensively investigate other potential intermediate hosts for the prevalence and genetic diversity of the E. canadensis group (G6, G7, G8 and G10) across the country and their inclusion into the existing CE control programme.

First report of Echinococcus granulosus sensu stricto (G1) in Nigeria, West Africa

Echinococcus granulosus sensu stricto is regarded to have the highest zoonotic potential of all Echinococcus taxa. Globally, human infection due to this species constitutes over 88.44% of the total cystic echinococcosis (CE) burden. Here, we report a CE infection in a Nigerian camel caused by E. granulosus G1 genotype. To the best of our knowledge, this report is the first encounter of the G1 genotype in the West Africa sub-region where the G6 genotype is reportedly prevalent, suggesting that the epidemiology of this highly zoonotic group could have a wider host range and distribution in the sub-region, and emphasizes the need for further investigation into the genetic diversity of Echinococcus spp. in Nigeria and across the sub-region.

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.