The complete mitochondrial genome of the causative agent of the human cercarial dermatitis, the visceral bird schistosome species Trichobilharzia szidati (platyhelminthes: Trematoda: Schistosomatidae)

A report on the complete mitochondrial genome of the trematode Azygia robusta Odhner, 1911, its new definitive host from the Russian Far East, and unexpected phylogeny of Azygiidae within Digenea, as inferred from mitogenome sequences

New data on the complete mitochondrial genome of Azygia robusta (Azygiidae) were obtained by the next-generation sequencing (NGS) approach. The mitochondrial DNA (mtDNA) of A. robusta had a length of 13 857 bp and included 12 protein-coding genes, two ribosomal genes, 22 transfer RNA genes, and two non-coding regions. The nucleotide sequences of the complete mitochondrial genomes of two A. robusta specimens differed from each other by 0.12 ± 0.03%. Six of 12 protein-coding genes demonstrated intraspecific variation. The difference between the nucleotide sequences of the complete mitochondrial genomes of A. robusta and Azygia hwangtsiyui was 26.95 ± 0.35%; the interspecific variation of protein-coding genes between A. robusta and A. hwangtsiyui ranged from 20.5 ± 0.9% (cox1) to 30.7 ± 1.2% (nad5). The observed gene arrangement in the mtDNA sequence of A. robusta was identical to that of A. hwangtsiyui. Codon usage and amino acid frequencies were highly similar between A. robusta and A. hwangtsiyui. The results of phylogenetic analyses based on mtDNA protein-coding regions showed that A. robusta is closely related to A. hwangtsiyui (belonging to the same suborder, Azygiida) that formed a distinct early-diverging branch relative to all other Digenea. A preliminary morphological analysis of paratypes of the two azygiid specimens studied showed visible morphological differences between them. The specimen extracted from Sakhalin taimen (Parahucho perryi) was most similar to A. robusta. Thus, we here provide the first record of a new definitive host, P. perryi, for A. robusta and also molecular characteristics of the trematode specimens.

Characterization of complete mitochondrial genome and ribosomal operon for Carassotrema koreanum Park, 1938 (Digenea: Haploporidae) by means of next-generation sequencing data

We obtained new data on the complete mitochondrial DNA (mtDNA) and the ribosomal operon of the trematode Carassotrema koreanum (Digenea: Haploporata: Haploporidae), an intestinal parasite of Carassius auratus, using next-generation sequencing. The mtDNA of C. koreanum contained 13,965 bp, including 12 protein-coding genes, two ribosomal genes, 22 transport RNA (tRNA) genes and a non-coding region. The ribosomal operon of C. koreanum was 10,644 bp in length, including ETS1 (1449 bp), 18S ribosomal RNA (rRNA) gene (1988 bp), ITS1 ribosomal DNA (rDNA) (558 bp), 5.8S rRNA gene (157 bp), ITS2 rDNA (274 bp), 28S rRNA gene (4152 bp) and ETS2 (2066 bp). Phylogenetic analysis based on mtDNA protein-coding regions showed that C. koreanum was closely related to Parasaccocoelium mugili, a species from the same suborder Haploporata. Bayesian phylogenetic tree topology was the most reliable and confirmed the validity of the Haploporata. The results of sequence cluster analysis based on codon usage bias demonstrated some agreement with the results of the phylogenetic analysis. In particular, Schistosoma spp. were differentiated from the other members of Digenea and the members of Pronocephalata were localized within the same cluster. Carassotrema koreanum and P. mugili fell within different clusters. The grouping of C. koreanum and P. mugili within the same cluster was obtained on the basis of frequencies of 13 specified codons, of which three codon pairs were degenerate. A similarity was found between two haploporid species and two Dicrocoelium spp. in the presence of TTG start codon of the mitochondrial nad5 gene. Our results confirmed the taxonomical status of the Haploporata identified in the previous studies and revealed some characteristic features of the codon usage in its representatives.

Non-coding Regions of Mitochondrial DNA and the cox1 Gene Reveal Genetic Variability Among Local Belarusian Populations of the Causative Agent of Cercarial Dermatitis, Bird Schistosome Trichobilharzia szidati (Digenea: Schistosomatidae)

INTRODUCTION

The cercariae of avian blood flukes Trichobilharzia szidati (Digenea, Schistosomatidae) are known to cause cercarial allergic dermatitis ("swimmer's itch") in humans. Global epidemics can have significant impacts on local tourism-related economies in recreational areas. Little is known about the genetic polymorphism of the parasite population, or about the variability of the non-coding regions of mitochondrial DNA (mtDNA) and the possibility of using this as a genetic marker.

MATERIALS AND METHODS

The T. szidati cercariae were collected over 7 years from 33 naturally infected Lymnaea stagnalis snails from five sites at two neighboring lakes in Belarus. We investigated the variability of the short (SNR) and long (LNR) non-coding regions of mt DNA and the genetic diversity within the 1125-bp sequences of the gene for subunit 1 of cytochrome c oxidase (cox1).

RESULTS

In the SNR sequences, we found only length variability caused by changes in the number of bases in the mononucleotide tracts T6-T8. LNR demonstrates high variability in nucleotide sequence length (182-260 bp) depending on the presence of two long deletions of 59 and 78 nucleotides. Both mitochondrial loci (LNR and cox1) are characterized by high haplotype diversity (H = 0.922 and H = 1.0, respectively); the nucleotide diversity is significantly higher for LNR (π = 1.926 ± 0.443) compared to cox1 (π = 0.704 ± 0.059). Phylogenetic reconstructions based on the variability of each of the loci (LNR and cox1) and their concatenated sequences revealed their shallow structure and the absence of a correlation between the distribution of single-nucleotide polymorphisms and the geographic origin of parasites from two Belarusian lakes. We identified at last four weakly sublineages in the phylogenetic pattern of T. szidati. The carriers of each deletion have specific patterns for each of the two loci and form their own phylogeographic sublineages. An association between two fixed LNR substitutions and a fixed non-synonymous substitution in cox1 was found in four representatives of one lineage that had a short deletion in the LNR.

CONCLUSIONS

This study clarified the phylogeographic structure of the Belarusian population of T. szidati. Our data provide the basis for the use two mt markers in large-scale population studies of the parasite, as well as for studying the molecular evolution of coding and non-coding mtDNA in trematodes.

First Record of Trichobilharzia physellae (Talbot, 1936) in Europe, a Possible Causative Agent of Cercarial Dermatitis

Several species of avian schistosomes are known to cause dermatitis in humans worldwide. In Europe, this applies above all to species of the genus Trichobilharzia. For Austria, a lot of data are available on cercarial dermatitis and on the occurrence of Trichobilharzia, yet species identification of trematodes in most cases is doubtful due to the challenging morphological determination of cercariae. During a survey of trematodes in freshwater snails, we were able to detect a species in the snail Physella acuta (Draparnaud, 1805) hitherto unknown for Austria, Trichobilharzia physellae; this is also the first time this species has been reported in Europe. Species identification was performed by integrative taxonomy combining morphological investigations with molecular genetic analyses. The results show a very close relationship between the parasite found in Austria and North American specimens (similarity found in CO1 ≥99.57%). Therefore, a recent introduction of T. physellae into Europe can be assumed.

First next-generation sequencing data for Haploporidae (Digenea: Haploporata): characterization of complete mitochondrial genome and ribosomal operon for Parasaccocoelium mugili Zhukov, 1971

The first data on a whole mitochondrial genome of Haploporidae, Parasaccocoelium mugili (Digenea: Haploporata: Haploporidae) was generated using the next-generation sequencing (NGS) approach. We sequenced the complete mitochondrial DNA (mtDNA) and ribosomal operon of Parasaccocoelium mugili, intestine parasite of mullet fish. The mtDNA of P. mugili contained 14,021 bp, including 12 protein-coding genes, two ribosomal genes, 22 tRNA genes, and non-coding region. The ribosomal operon of P. mugili was 8308 bp in length, including 18S rRNA gene (1981 bp), ITS1 rDNA (955 bp), 5.8S rRNA gene (157 bp), ITS2 rDNA (268 bp), 28S rRNA gene (4180 bp), and ETS (767 bp). We used the mtDNA protein-coding regions to make phylogenetic reconstructions of Haploporidae. Additionally, we performed the sequence cluster analysis based on codon usage bias of most of currently available mitochondrial genome data for trematodes. The observed gene arrangement in mtDNA sequence of P. mugili is identical to those of Plagiorchis maculosus (Rudolphi, 1802). Results of maximum likelihood (ML) phylogenetic analysis showed that P. mugili was closely related to Paragonimus species from the suborder Xiphidiata. The results of sequence cluster analysis based on codon usage bias showed that P. mugili has the highest similarity with Plagiorchis maculosus (Xiphidiata). Our results do not contradict to proposing a new suborder for Haploporoidea-Haploporata. On the basis of obtained results, the relationship between mitochondrial protein-coding gene rearrangements and synonymous nucleotide substitutions in mitochondrial genomes has been suggested.

HERMES: An improved method to test mitochondrial genome molecular synapomorphies among clades

Mitochondrial chromosomes have diversified among eukaryotes and many different architectures and features are now acknowledged for this genome. Here we present the improved HERMES index, which can measure and quantify the amount of molecular change experienced by mitochondrial genomes. We test the improved approach with ten molecular phylogenetic studies based on complete mitochondrial genomes, representing six bilaterian Phyla. In most cases, HERMES analysis spotted out clades or single species with peculiar molecular synapomorphies, allowing to identify phylogenetic and ecological patterns. The software presented herein handles linear, circular, and multi-chromosome genomes, thus widening the HERMES scope to the complete eukaryotic domain.

Characterization and comparative analysis of the complete mitochondrial genome of Azygia hwangtsiyui Tsin, 1933 (Digenea), the first for a member of the family Azygiidae

Azygia hwangtsiyui (Trematoda, Azygiidae), a neglected parasite of predatory fishes, is little-known in terms of its molecular epidemiology, population ecology and phylogenetic study. In the present study, the complete mitochondrial genome of A. hwangtsiyui was sequenced and characterized: it is a 13,973 bp circular DNA molecule and encodes 36 genes (12 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes) as well as two non-coding regions. The A+T content of the A. hwangtsiyui mitogenome is 59.6% and displays a remarkable bias in nucleotide composition with a negative AT skew (-0.437) and a positive GC skew (0.408). Phylogenetic analysis based on concatenated amino acid sequences of twelve protein-coding genes reveals that A. hwangtsiyui is placed in a separate clade, suggesting that it has no close relationship with any other trematode family. This is the first characterization of the A. hwangtsiyui mitogenome, and the first reported mitogenome of the family Azygiidae. These novel datasets of the A. hwangtsiyui mt genome represent a meaningful resource for the development of mitochondrial markers for the identification, diagnostics, taxonomy, homology and phylogenetic relationships of trematodes.

Mitochondrial genome characterization and phylogenetic analysis of bird schistosome Trichobilharzia szidati

In the present study, the complete mitochondrial genome of T. szidati was assembled by next generation sequencing (NGS). We found that the complete mitochondrial genome of T.szidati is 14, 303 bp in length and consists of 3023 (21.1%) adenine, 1153 (8.1%) cytosine, 3432 (24.0%) guanosine and 6695 (46.8%) thymine. The genome contains 12 conserved core protein-coding genes (atp6, cox1, cox2, cox3, nad1, nad2, nad3, nad4, nad4L, nad5, nad6, CYTB), 21 tRNA genes, 2 rRNA genes and 1 D-loop region. Phylogenetic analysis showed that T. szidati has a close relationship with T. regent. Knowledge of mitochondrial genome of T. szidati could provide useful information for the further studies of evolutionary biology, epidemiology and species identification.