Intra- and interindividual variation in ITS1 of paragonimus westermani (Trematoda: digenea) and related species: implications for phylogenetic studies

The ribosomal transcription units of five echinostomes and their taxonomic implications for the suborder Echinostomata (Trematoda: Platyhelminthes)

Five newly obtained nuclear ribosomal transcription unit (rTU) sequences from Echinostomatidae and Echinochasmidae are presented. The inter- and intrafamilial relationships of these and other families in the suborder Echinostomata are also analyzed. The sequences obtained are the complete rTU of Artyfechinostomum malayanum (9,499 bp), the near-complete rTU of Hypoderaeum conoideum (8,076 bp), and the coding regions (from 5'-terminus of 18S to 3'-terminus of 28S rRNA gene) in Echinostoma revolutum (6,856 bp), Echinostoma miyagawai (6,854 bp), and Echinochasmus japonicus (7,150 bp). Except for the longer first internal transcribed spacer (ITS1) in Echinochasmus japonicus, all genes and spacers were almost identical in length. Comprehensive maximum-likelihood phylogenies were constructed using the PhyML software package. The datasets were either the concatenated 28S + 18S rDNA sequences (5.7-5.8 kb) from 60 complete rTUs of 19 families or complete 28S sequences only (about 3.8-3.9 kb) from 70 strains or species of 22 families. The phylogenetic trees confirmed Echinostomatoidea as monophyletic. Furthermore, a detailed phylogeny constructed from alignments of 169 28S D1-D3 rDNA sequences (1.1-1.3 kb) from 98 species of 50 genera of 10 families, including 154 echinostomatoid sequences (85 species/42 genera), clearly indicated known generic relationships within Echinostomatidae and Echinochasmidae and relationships of families within Echinostomata and several other suborders. Within Echinostomatidae, Echinostoma, Artyfechinostomum, and Hypoderaeum appeared as monophyletic, while Echinochasmus (Echinochasmidae) was polyphyletic. The Echinochasmidae are a sister group to the Psilostomidae. The datasets provided here will be useful for taxonomic reappraisal as well as studies of evolutionary and population genetics in the superfamily Echinostomatoidea, the sole superfamily in the suborder Echinostomata.

Potential intermediate hosts of Angiostrongylus cantonensis in the European Mediterranean region (Mallorca, Spain)

Angiostrongylus cantonensis is the main etiological agent of eosinophilic meningitis, a disease that often leads to severe neurological manifestations in mammals and birds. In recent years, the prevalence of this zoonotic nematode has dramatically increased as it expands into new territories beyond its native range in Southeast Asia and the Pacific Basin.Its arrival to Europe, the last continent to be invaded, has raised many questions concerning the parasite's life cycle, particularly in the Mediterranean region, where it is now endemic. This study aims to collect evidence about intermediate hosts (snail and slug species) involved in the transmission of the A. cantonensis in Mallorca. We have conducted a systematic surveillance of gastropods within 2 km radius areas, surrounding the specific locations where infected hedgehogs were found. We employed a sequence-based detection approach that included a species-specific PCR amplification followed by DNA sequencing of the internal transcribed spacer I (ITS-1). This conservative approach was essential to avoid cross reactions with the numerous metastrongylid species known to be circulating in Mallorca. Overall, we collected and identified 398 gastropods belonging to 17 species, of which 11% tested positive for A. cantonensis. These gastropods were collected from touristic settlements to agricultural lands. Five gastropod species: Theba pisana, Cornu aspersum, Deroceras reticulatum, Milax gagates and Otala lactea, had been previously incriminated as A. cantonensis intermediate hosts, while 6 species: Papillifera papillaris, Cochlicella acuta, Eobania vermiculata, Ganula lanuginosa, Milax nigricans and Rumina decollata, are newly recorded hosts. The findings of this study have important epidemiological implications, and further measures are discussed to prevent neuroangiostrongylosis cases.

Mitogenomic and nuclear ribosomal transcription unit datasets support the synonymy of Paragonimus iloktsuenensis and P. ohirai (Paragonimidae: Platyhelminthes)

The complete mitogenome (mtDNA) of nominal Paragonimus iloktsuenensis (Paragonimidae: Trematoda) and the nuclear ribosomal transcription unit (rTU) coding region (rTU*: from 5'-terminus of 18S to 3'-terminus of 28S rRNA gene, excluding the external spacer region) of this species and of P. ohirai were obtained and used to further support the previously suggested synonymy of these taxa in the P. ohirai complex. The complete mitogenome of P. iloktsuenensis was 14,827 bp long (GenBank: ON961029) and nearly identical to that of P. ohirai (14,818 bp; KX765277), with a 99.12% nucleotide identity. The rTU* was 7543 bp and 6932 bp in these two taxa, respectively. All genes and spacers in the rTU were identical in length, with exception of the first internal transcribed spacer, which contained multiple tandem repeat units (6.7 for P. iloktsuenensis and 5.7 for P. ohirai). There was near 100% identity for the rTU genes. The phylogenetic topology inferred from the mtDNA and from individual gene regions (partial cox1 of 387 bp and the ITS-2 of 282 bp - 285 bp) indicated a very close relationship consistent with synonymy of P. iloktsuenensis and P. ohirai. The datasets provided here will be useful for taxonomic reappraisal as well as studies of evolutionary and population genetics of the genus Paragonimus and family Paragonimidae.

Intrageneric Relationship of Datnioides (Lobotiformes) Inferred from the Complete Nuclear Ribosomal DNA Operon

Tiger fish (genus Datnioides) are critical ornamental and economic fish and are valuable freshwater fish worldwide, belonging to the order Lobotiformes. Currently, there are five extant species (Datnioides campbelli, D. microlepis, D. polota, D. pulcher, and D. undecimradiatus) of Datnioides in the world, usually inhabiting in south and southeast Asia. Due to the decline of wild population sizes of tiger fish and the lack of molecular research on them, in the present study, we sequenced, assembled, and characterized the complete nuclear ribosomal DNA (nrDNA) operon of all five extant tiger fish species, in order to elucidate the phylogenetic relationship among the genus Datnioides. The nrDNA sequences of five tiger fish species were 8548-9182 bp in length, encompassing complete 18S rDNA, ITS1, 5.8S rDNA, ITS2, 28S rDNA, and IGS regions. Numerous repetitive sequences were detected, substantially influencing the sequence length of different regions in each species. We employed maximum-likelihood (ML) method and Bayesian inference (BI) method to construct phylogenetic trees for Datnioides. Phylogenetic analyses indicated that each region in nrDNA operon is not sufficiently phylogenetically informative to delineate the species in Datnioides; nevertheless, the whole operon is able to delineate five tiger fish species much better, three of five species were successfully partitioned. Particularly, regardless of employed markers, it was strongly supported that D. campbelli was considerably partitioned from the other four species, possibly due to the geographical separation. In spite of the fact that discrimination of Datnioides species requires further investigation, our study provides reference genome resources for the Lobotiformes, as well as insights into the phylogenetic position of Lobotiformes and further biological conservation.

Identification and Genetic Characterization of Fasciola hepatica Isolated from Cattle in Jeddah, Saudi Arabia Based on Sequence Analysis of Mitochondrial (COI) Gene

Background

In Saudi Arabia, more than US$ 0.2 million annual losses are caused by liver condemnations due to fascioliasis. Data obtained from the genetic characterization of Fasciola population sheds light on parasite transmission which could eventually help in development of effective parasite control measures. So, the aim of this study was to investigate the genetic diversity of Fasciola spp. isolated from cattle in Saudi Arabia by sequence analyses of COI gene.

Materials and Methods

A total of 325 cows slaughtered at the central municipal abattoir in Jeddah city, Jeddah Province, Saudi Arabia were examined for fascioliasis in the period from 1^st of June to 1^st of July 2020. DNA was extracted from adult Fasciola worms and used for PCR and DNA sequence using a primer pair targeting COI gene. Analysis of the obtained sequences was done using BLAST search and phylogenetic analysis.

Results

Bovine fascioliasis was diagnosed in 18 out of 325 cattle (5.5%). Forty-eight flukes were extracted from infected animals and DNA was successfully amplified from all flukes. Overall 12 different DNA sequences were obtained. BLAST search showed that all obtained sequences were F. hepatica and had >97% similarity with F. hepatica isolates from Tanzania, Europe and Iran. Phylogenetic analysis of the obtained sequences showed that Fasciola isolates from the current study were clustered in one subclade closely related to isolates from North and South Africa and Italy.

Conclusion

Reports on the molecular characterization of Fasciola spp. in Saudi Arabia are limited. In the current study, our findings showed that F. hepatica was the only Fasciola species parasitizing cattle in Jeddah city, Saudi Arabia. Further studies using a large number of samples from different localities in Saudi Arabia are needed to provide data that will help the development of control measures against fascioliasis.

What lies behind the curtain: Cryptic diversity in helminth parasites of human and veterinary importance

Parasite cryptic species are morphologically indistinguishable but genetically distinct organisms, leading to taxa with unclear species boundaries. Speciation mechanisms such as cospeciation, host colonization, taxon pulse, and oscillation may lead to the emergence of cryptic species, influencing host-parasite interactions, parasite ecology, distribution, and biodiversity. The study of cryptic species diversity in helminth parasites of human and veterinary importance has gained relevance, since their distribution may affect clinical and epidemiological features such as pathogenicity, virulence, drug resistance and susceptibility, mortality, and morbidity, ultimately affecting patient management, course, and outcome of treatment. At the same time, the need for recognition of cryptic species diversity has implied a transition from morphological to molecular diagnostic methods, which are becoming more available and accessible in parasitology. Here, we discuss the general approaches for cryptic species delineation and summarize some examples found in nematodes, trematodes and cestodes of medical and veterinary importance, along with the clinical implications of their taxonomic status. Lastly, we highlight the need for the correct interpretation of molecular information, and the correct use of definitions when reporting or describing new cryptic species in parasitology, since molecular and morphological data should be integrated whenever possible.

•

Cryptic diversity has been described in helminths of human and animal importance.

•

Cryptic species are morphologically indistinguishable but genetically distinct organisms.

•

These entities emerge by different evolutionary and speciation mechanisms.

•

Analysis of molecular and morphological data is needed for cryptic species delimitation.

•

Cryptic diversity may affect pathogenicity, virulence and drug resistance of helminths.

Exploring the genetic structure of Parastrigea diovadena Dubois and Macko, 1972 (Digenea: Strigeidae), an endoparasite of the white ibis, Eudocimus albus, from the Neotropical region of Mexico

Parastrigea diovadena Dubois and Macko, 1972, is an allogenic trematode species that infects the intestine of white ibis. This widely distributed Neotropical species has been studied poorly, and nothing is known about its population genetic structure. In the current study, we attempt to fill this gap for the first time and to explore the genetic diversity in P. diovadena populations from three biogeographic provinces (Sierra Madre Oriental, Sierra Madre Occidental, and Sierra Madre del Sur) in the Neotropical region of Mexico. Newly generated sequences of the internal transcribed spacers (ITS) from ribosomal DNA and cytochrome c oxidase subunit 1 (cox 1) from mitochondrial DNA were compared with sequences available from the GenBank data set. Phylogenetic analyses performed with the ITS and cox 1 data sets using maximum likelihood and Bayesian inference unequivocally showed that new sequences of P. diovadena recovered from the white ibis formed a clade with other sequences of specimens previously identified as P. diovadena. The intraspecific genetic divergence among the isolates was very low, ranging from 0 to 0.38% for ITS and from 0 to 1.5% for cox 1, and in combination with the phylogenetic trees confirmed that the isolates belonged to the same species. The cox 1 haplotype network (star-shaped) inferred with 62 sequences revealed 36 haplotypes. The most frequent haplotype (H3, n = 18) corresponded to specimens from all the populations (except Tecolutla, Veracruz). In addition to the common haplotype, we identified four other shared haplotypes (H2, H9, H12, and H14) and 31 unique haplotypes (singlets). In addition, high haplotype diversity (Hd = 0.913), low nucleotide diversity (Pi = 0.0057), and null genetic differentiation or population structure (Fst = 0.0167) were found among the populations from the three biogeographic provinces. The results suggest that the biology of the definitive host has played a key role in the population genetic structure of Parastrigea diovadena in the Neotropical region of Mexico.

Comparative analyses of complete ribosomal DNA sequences of Clonorchis sinensis and Metorchis orientalis: IGS sequences may provide a novel genetic marker for intraspecific variation

Both Clonorchis sinensis and Metorchis orientalis are the fish-borne zoonotic trematodes, and have a wide distribution of southeastern Asia, especially in China. Due to the similar morphology, life cycle, and parasitic positions are difficult to differentiate between both metacercariae. In the present study, the complete rDNA sequences of five C. sinensis and five M. orientalis were obtained and compared for the first time. And the IGS rDNA sequences were tested as a genetic marker. The results showed complete rDNA lengths of C. sinensis were range from 8049 bp to 8391 bp, including 1991 bp, 1116 bp, 3854 bp, and 1088-1430 bp belonging to 18S, ITS, 28S and IGS, respectively. And the complete rDNA lengths of M. orientalis were range from 7881 bp to 9355 bp, including 1991 bp, 1077 bp, 3856 bp, and 957-2431 bp belonging to 18S, ITS, 28S and IGS, respectively. Comparative analyses reveal length difference main in IGS, which has higher intraspecific and interspecific variations than other ribosomal regions. Forty four repeat (forward and inverted) sequences were found in the complete rDNAs of C. sinensis and M. orientalis. The phylogenetic analyses showed that the sequences of ITS1, ITS2, 18S and 28S could be used as different level genetic markers. In IGS phylogenetic tree, Opisthorchiidae, Paramphistomidae, Dicrocoeliidae, and Schistosomatidae formed monophyletic groups, and the same length sequences were clustered together in the same species. These findings of the present study provide the new molecular data for studying the complete rDNA of C. sinensis and M. orientalis, and indicate IGS sequences may used as a novel genetic marker for studying intraspecific variation in trematodes.

Recent and old duplications in crustaceans "Internal Transcribed Spacer 1″: structural and phylogenetic implications

Internal Transcribed Spacer structures are important in preserving accessibility to specific enzymes for the maturation of rRNAs. ITS1 sequences reported in the literature in Crustaceans range between 182 and 820 bp and are characterized by the absence of repeats or the presence of only a limited number of microsatellites. Here, we sequenced ITS1 for a range of shrimp families (infraorder Caridea) and show that most taxa have much larger ITS1 sequences. We find a high number of microsatellites in Alpheus hebes and Crangon crangon and we report repeat units in Pandalidae, Palaemonidae and mainly in Alpheidae species. Up to four repeats were found in A. vanderbilti (1915 bp), A. rostratus (1635 bp) and A. lottini (1625 bp). In general, four helices were found in ITS1. Repeat units led to extra hairpins and loops. No conserved positions occurred except in helix 4. Three clades were defined in A. lottini for the first time. We estimated the ITS1 divergence rate for the three clades of A. lottini collected in French Polynesia using existing calibrations of substitution rates. Rates of sequence evolution are largely influenced by repeat units, which likely evolve separately. By comparison with COI marker, we estimated the divergence rate of the whole ITS1 sequence to range from 0.5 to 1.4% Pmy and between 0.12 and 0.5% for the 3' end of ITS1 located outside the repeat units. Given the degree of identity between repeats, we suggest that a duplication event recently occurred in A. floridanus (98% identity) whereas an ancient duplication happened in A. sulcatus (50% identity) early at the origination of the group Alpheidae, approximately 50 mya ago. In conclusion, our results highlight an over representation of shorter ITS1 sequences in public repositories, and underlines the importance to further understand patterns of molecular evolution of this functionally important gene.

The complete mitochondrial genome of Paragonimus ohirai (Paragonimidae: Trematoda: Platyhelminthes) and its comparison with P. westermani congeners and other trematodes

We present the complete mitochondrial genome of Paragonimus ohirai Miyazaki, 1939 and compare its features with those of previously reported mitochondrial genomes of the pathogenic lung-fluke, Paragonimus westermani, and other members of the genus. The circular mitochondrial DNA molecule of the single fully sequenced individual of P. ohirai was 14,818 bp in length, containing 12 protein-coding, two ribosomal RNA and 22 transfer RNA genes. As is common among trematodes, an atp8 gene was absent from the mitogenome of P. ohirai and the 5' end of nad4 overlapped with the 3' end of nad4L by 40 bp. Paragonimusohirai and four forms/strains of P. westermani from South Korea and India, exhibited remarkably different base compositions and hence codon usage in protein-coding genes. In the fully sequenced P. ohirai individual, the non-coding region started with two long identical repeats (292 bp each), separated by tRNAGlu . These were followed by an array of six short tandem repeats (STR), 117 bp each. Numbers of the short tandem repeats varied among P. ohirai individuals. A phylogenetic tree inferred from concatenated mitochondrial protein sequences of 50 strains encompassing 42 species of trematodes belonging to 14 families identified a monophyletic Paragonimidae in the class Trematoda. Characterization of additional mitogenomes in the genus Paragonimus will be useful for biomedical studies and development of molecular tools and mitochondrial markers for diagnostic, identification, hybridization and phylogenetic/epidemiological/evolutionary studies.

How to become a successful invasive tapeworm: a case study of abandoned sexuality and exceptional chromosome diversification in the triploid carp parasite Atractolytocestus huronensis Anthony, 1958 (Caryophyllidea: Lytocestidae)

Background

A cytogenetic analysis of the new local triploid population of the caryophyllidean tapeworm Atractolytocestus huronensis, a unique parthenogenetic species with the ability to colonise new regions, was performed to understand the inner structure of its chromosome complement.

Methods

A karyotype analysis was carried out using classical Giemsa staining and C-banding combined with fluorescent DAPI staining. A hypothesis that triplets are composed from three homologue chromosomes of approximately the same length and same centromere position was tested statistically for multiple dependent variables using a non-parametric Friedman’s ANOVA. The chromosomal location of ribosomal DNA clusters within the nucleolar organization region (NORs) and telomeric (TTAGGG)_n sequences were detected by fluorescent in situ hybridization (FISH). Chromosomes were subjected to AgNO₃ staining in order to determine whether the rDNA sites represent active NORs.

Results

The cytogenetic analysis confirmed the karyotype composed from eight chromosome triplets (3n = 24) as well as the existence of a pair of NORs located on each chromosome of the second triplet. Six NORs varied their activity from cell to cell, and it was reflected in the numbers of nucleoli (from 1 to 5). A huge morphological diversification of homologue chromosomes was originally detected in six out of eight triplets; the homologue elements differed significantly either in length and/or morphology, and some of them carried discernible interstitial telomeric sequences (ITSs), while the end telomeres were minute. The heterochromatin bands with high AT content varied irregularly, and the course of aberrant spermatogenesis was evident.

Conclusions

Diversification of homologues is a unique phenomenon very likely caused by the long-term absence of a recombination and consequential accumulation of chromosome rearrangements in the genome of A. huronensis during species evolution. Unalterable asexual reproduction of the tapeworm, along with international trade in its host (carp), is facilitating its ongoing spread.

Parasite detection in the ornamental fish trade using environmental DNA

Effective border control relies on stringent biosecurity protocols to detect and prevent introductions of exotic pests and diseases. Detection of pathogens and parasites in the live ornamental fish trade using environmental DNA (eDNA) techniques has the potential to improve current biosecurity practices. We examined water samples from 11 target consignments (cyprinids susceptible to Dactylogyrus spp. infections) and seven non-target fish consignments (non-cyprinids, not susceptible to Dactylogyrus spp. infections) imported from Southeast Asia to Australia for the presence of eDNA from five Dactylogyrus species (Monogenea: Dactylogyridae). A four-step predictive framework was used to predict putative positive and putative negative detections from quantitative PCR assays. Both target and non-target consignments were positive for Dactylogyrus spp. eDNA as confirmed by Sanger sequencing. Positive detections for Dactylogyrus spp. eDNA in non-target fish consignments demonstrates the possibility of source water contamination, limiting the applicability of eDNA screening methods at border control. This study suggests that screening for parasite eDNA within ornamental fish consignments should be tested during pre-export quarantine periods to avoid false positive detections at border control. Lastly, the proposed predictive framework has a broad utility for minimizing false positive and false negative eDNA detections of aquatic organisms.

Phylogenetic analysis of Spirocerca lupi and Spirocerca vulpis reveal high genetic diversity and intra-individual variation

Background

Spirocerca lupi is a parasitic nematode of canids that can lead to a severe and potentially fatal disease. Recently, a new species, Spirocerca vulpis, was described from red foxes in Europe, suggesting a high genetic diversity of the Spirocerca spp. infecting canids. The genetic variation and phylogenetic relationships of S. lupi collected from naturally-infected domestic dogs from Australia, Hungary, Israel, Italy, India and South Africa, and S. vulpis from red foxes from Bosnia and Herzegovina, Italy and Spain, was studied using mitochondrial and rDNA markers.

Results

A high intra-individual variation was found in the first internal transcribed spacer (ITS1) locus in all Spirocerca spp., ranging between 0.37–2.84%, with up to six haplotypes per specimen. In addition, a combination of phylogenetic and haplotype analyses revealed a large variability between S. lupi specimens collected from different geographical locations using the ITS1 (0.37–9.33%) and the cytochrome c oxidase subunit 1 (cox1) gene (1.42–6.74%). This genetic diversity led to the identification of two S. lupi genotypes circulating among dogs (PTP support > 0.829), including genotype 1 found in S. lupi from Australia, India, Israel and South Africa, and genotype 2 represented by specimens from Hungary and Italy. These genotypes presented pairwise nucleotide distances of 0.14%, 8.06% and 6.48 ± 0.28% in the small rDNA subunit (18S), ITS1 and cox1 loci, respectively. Additionally, Nei’s genetic distance in the ITS1 showed a further subdivision of genotype 1 worms into 1A (Israel and South Africa) and 1B (Australia and India). A morphological analysis of the anterior and posterior extremities of genotype 1 and genotype 2 worms using scanning electron microscopy did not show any differences between the specimens, contrary to the morphological differences between S. lupi and S. vulpis.

Conclusions

These findings demonstrate the high genetic variability among Spirocerca spp. from different geographical locations, thereby expanding our understanding of the epidemiology, evolution and phylogenetic variability within the genus.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3202-0) contains supplementary material, which is available to authorized users.

Tour around the globe: The case of invasive tapeworm Atractolytocestus huronensis (Cestoda: Caryophyllidea), a parasite of common carp

The monozoic tapeworm Atractolytocestus huronensis Anthony, 1958 (Cestoda: Caryophyllidea), an intestinal parasite of the common carp, is characterized by its invasive character and potential to colonize new territories. It was initially described from North America and has also been found in several European countries. The most recent findings of A. huronensis originated from China and South Africa; however, no data on genetic relationships of these populations were available. The current study provides the first molecular characterisation of A. huronensis from South Africa and China using a partial sequence of mitochondrial cytochrome c oxidase subunit 1 (cox1) and a complete ribosomal ITS2 spacer. Ribosomal and mitochondrial data were applied for phylogenetic analyses in order to assess the genetic interrelationships among global A. huronensis populations. Divergent intragenomic copies of ribosomal ITS2 were detected in all analysed specimens; the structure and frequency of the ITS2 variants of tapeworms from China and South Africa corresponded with the data on ITS2 paralogues observed previously in A. huronensis from Slovakia, the United States and the United Kingdom. The phylogenetic analysis of cox1 indicated that A. huronensis exist in two slightly differentiated clusters; one cluster was supported by all phylogenetic approaches (NJ, ML, BI) and was represented by samples from China, the USA and the UK. A second cluster was represented by tapeworms from continental Europe (Slovakia, Hungary, Romania, Croatia) and South Africa. Haplotype network analysis revealed that the highest population diversity occurs in China. The results provide useful pilot information about the interrelationships of A. huronensis on four continents and indicate that China, or the eastern Palaearctic, served as the original source population for the global expansion of this invasive tapeworm. Data on the origin and distribution of the common carp, the only specific host of A. huronensis, are also discussed.

Prevalence and Identity of Taenia multiceps cysts "Coenurus cerebralis" in Sheep in Egypt

Coenurosis is a parasitic disease caused by the larval stage (Coenurus cerebralis) of the canids cestode Taenia multiceps. C. cerebralis particularly infects sheep and goats, and pose a public health concerns. The present study aimed to determine the occurrence and molecular identity of C. cerebralis infecting sheep in Egypt. Infection rate was determined by postmortem inspection of heads of the cases that showed neurological manifestations. Species identification and genetic diversity were analyzed based on PCR-sequence analysis of nuclear ITS1 and mitochondrial cytochrome oxidase (COI) and nicotinamide adenine dinucleotide dehydrogenase (ND1) gene markers. Out of 3668 animals distributed in 50 herds at localities of Ashmoun and El Sadat cities, El Menoufia Province, Egypt, 420 (11.45%) sheep showed neurological disorders. Postmortem examination of these animals after slaughter at local abattoirs indicated to occurrence of C. cerebralis cysts in the brain of 111 out of 420 (26.4%), with overall infection rate 3.03% of the involved sheep population. Molecular analysis of representative samples of coenuri at ITS1 gene marker showed extensive intra- and inter-sequence diversity due to deletions/insertions in the microsatellite regions. On contrast to the nuclear gene marker, considerably low genetic diversity was seen in the analyzed mitochondrial gene markers. Phylogenetic analysis based on COI and ND1 gene sequences indicated that the generated sequences in the present study and the reference sequences in the database clustered in 4 haplogroups, with more or less similar topologies. Clustering pattern of the phylogenetic tree showed no effect for the geographic location or the host species.

Characterization of the complete nuclear ribosomal DNA sequences of Eurytrema pancreaticum

Eurytrema pancreaticum is one of the most common trematodes of cattle and sheep, and also infects humans occasionally, causing great economic losses and medical costs. In this study, the sequences of the complete nuclear ribosomal DNA (rDNA) repeat units of five E. pancreaticum individuals were determined for the first time. They were 8306–8310 bp in length, including the small subunit (18S) rDNA, internal transcribed spacer 1 (ITS1), 5.8S rDNA, internal transcribed spacer 2 (ITS2), large subunit (28S) rDNA and intergenic spacer (IGS). There were no length variations in any of the investigated 18S (1996 bp), ITS1 (1103 bp), 5.8S (160 bp), ITS2 (231 bp) or 28S (3669 bp) rDNA sequences, whereas the IGS rDNA sequences of E. pancreaticum had a 4-bp length variation, ranging from 1147 to 1151 bp. The intraspecific variations within E. pancreaticum were 0–0.2% for 18S rDNA, 0–0.5% for ITS1, 0% for 5.8S rDNA and ITS2, 0–0.2% for 28S rDNA and 2.9–20.2% for IGS. There were nine types of repeat sequences in ITS1, two types in 28S rDNA, but none in IGS. A phylogenetic analysis based on the 18S rDNA sequences classified E. pancreaticum in the family Dicrocoeliidae of Plagiorchiata, closely related to the suborder Opisthorchiata. These results provide useful information for the further study of Dicrocoeliidae trematodes.

Genetic divergence of human pathogens Nanophyetus spp. (Trematoda: Troglotrematidae) on the opposite sides of the Pacific Rim

Human and animal nanophyetiasis is caused by intestinal flukes belonging to the genus Nanophyetus distributed on both North American and Eurasian coasts of Northern Pacific. In spite of the wide geographical distribution and medical and veterinary importance of these flukes, the intra-generic taxonomy of Nanophyetus spp. remains unresolved. The two most widely distributed nominal species, Nanophyetus salmincola and Nanophyetus schikhobalowi, both parasitizing humans and carnivorous mammals, were described from North America and eastern Eurasia, respectively. However, due to their high morphological similarity their interrelationships remained unclear and taxonomic status unstable. In this study, we explored genetic diversity of Nanophyetus spp. from the Southern Russian Far East in comparison with that of samples from North America based on the sequence variation of the nuclear ribosomal gene family (18S, internal transcribed spacers, ITS1-5·8S-ITS2 and 28S). High levels of genetic divergence in each rDNA region (nucleotide substitutions, indels, alterations in the secondary structures of the ITS1 and ITS2 transcripts) as well as results of phylogenetic analysis provided strong support for the status of N. salmincola and N. schikhobalowi as independent species.

Molecular characterization of Stictodora tridactyla (Trematoda: Heterophyidae) from Kuwait Bay using rDNA ITS and mtCO1

Stictodora tridactyla is an intestinal fluke in the family Heterophyidae that parasitizes shorebirds and mammals, including humans. Its metacercarial cyst stage was reported in the Arabian killifish, Aphanius dispar, at Kuwait Bay. In the present study, Cerithidea cingulata was found to serve as the first intermediate host of S. tridactyla. In order to establish the snail-fish link in the life cycle of S. tridactyla, complete sequences of ribosomal DNA internal transcribed spacer region 1 and 2 (rDNA ITS1 and ITS2) and partial sequence of cytochrome oxidase subunit 1 were obtained for metacercarial cysts isolated from the fish A. dispar and rediae isolated from the snail C. cingulata. Sequence alignment demonstrated that these larval stages belong to the same heterophyid species, S. tridactyla. Phylogenetic analysis based on rDNA ITS1, ITS2, and mtCO1 confirmed the position of S. tridactyla within the Heterophyidae and found it to cluster with Haplorchis spp. The present study represents the first molecular study correlating the larval stages of S. tridactyla using rDNA ITS1, ITS2, and mtCO1 and examining the phylogenetic relationships of S. tridactyla with different heterophyid species.

Characterization of the complete nuclear ribosomal DNA sequences of Paramphistomum cervi

Sequences of the complete nuclear ribosomal DNA (rDNA) gene from five individual Paramphistomum cervi were determined for the first time. The five complete rDNA sequences, which included the 18S rDNA, the internal transcribed spacer 1 (ITS1), the 5.8S rDNA, the internal transcribed spacer 2 (ITS2), the 28S rDNA, and the intergenic spacer (IGS) regions, had a length range of 8,493–10,221 bp. The lengths of the investigated 18S, ITS1, 5.8S, ITS2, and 28S rDNA sequences, which were 1,994 bp, 1,293 bp, 157 bp, 286 bp, and 4,186 bp, respectively, did not vary. However, the IGS rDNA sequences had a length range of 577–2,305 bp. The 5.8S and ITS-2 rDNA sequences had 100% identity among the five investigated samples, while the identities among the IGS had a range of 53.7–99.8%. A comparative analysis revealed that different types and numbers of repeats were found within each ITS1 and IGS region, which may be related to the length polymorphism of IGS. The phylogenetic position of P. cervi in Paramphistomatidae was analyzed based on the 18S rDNA sequences. These results will aid in studying the intra- and interspecific variation of the Paramphistomatidae and the systematics and phylogenetics of Digenea.

Integrative taxonomy identifies a new species of Phyllodistomum (Digenea: Gorgoderidae) from the twospot livebearer, Heterandria bimaculata (Teleostei: Poeciliidae), in Central Veracruz, Mexico

Phyllodistomum inecoli n. sp. is described from the twospot livebearer, Heterandria bimaculata (Teleostei: Poeciliidae), collected in the Río La Antigua basin, Veracruz, Mexico. The new species is described and characterised by using a combination of morphology, scanning electron microscopy, and sequences of nuclear and mitochondrial genes. Diagnostic characters of the new species of Phyllodistomum include a genital pore opening at the level of the caecal bifurcation; oval vitellarium, situated just posterior to the ventral sucker and not extended laterally and anterior extracaecal uterine loops variable in extension (reaching the anterior, median or posterior margin of the ventral sucker). P. inecoli n. sp. most closely resembles P. brevicecum, a species described as a parasite of the central mudminnow, Umbra limi, in other parts of North America; however, the genital pore in P. brevicecum is situated between the caecal bifurcation and the ventral sucker, the ovary is larger, the vitellarium is lobed and extended laterally and the anterior portion of the uterus extends to the posterior margin of the ventral sucker. Comparison of about 1,500–2,200 nucleotides of cox1 and 28S rDNA and ITS1 strongly supports the status of P. inecoli as a new species. Bayesian inference analysis of combined datasets of 28S rDNA and cox1 sequences showed that P. inecoli n. sp. and the other species found in freshwater fishes of Mexico, including the species complex of P. lacustri, are not sister species. Phylogenetic analysis based on 28S rDNA sequences of several gorgoderid taxa revealed the close relationship of P. inecoli n. sp. with several species of Phyllodistomum, Gorgodera and Gorgoderina with cystocercous cercariae developing in sphaeriid bivalves. Dot-plot analysis of ITS1 sequences of P. inecoli n. sp. revealed the presence of eight repetitive elements with different length, which together represent almost half the length of ITS1.

A review of molecular approaches for investigating patterns of coevolution in marine host-parasite relationships

Parasites and their relationships with hosts play a crucial role in the evolutionary pathways of every living organism. One method of investigating host-parasite systems is using a molecular approach. This is particularly important as analyses based solely on morphology or laboratory studies of parasites and their hosts do not take into account historical evolutionary interactions that can shape the distribution, abundance and population structure of parasites and their hosts. However, the predominant host-parasite coevolution literature has focused on terrestrial hosts and their parasites, and there still is a lack of studies in marine environments. Given that marine systems are generally more open than terrestrial ones, they provide fascinating opportunities for large-scale (as well as small-scale) geographic studies. Further, patterns and processes of genetic structuring and systematics are becoming more available across many different taxa (but especially fishes) in many marine systems, providing an excellent basis for examining whether parasites follow host population/species structure. In this chapter, we first highlight the factors and processes that challenge our ability to interpret evolutionary patterns of coevolution of hosts and their parasites in marine systems at different spatial, temporal and taxonomic scales. We then review the use of the most commonly utilized genetic markers in studying marine host-parasite systems. We give an overview and discuss which molecular methodologies resolve evolutionary relationships best and also discuss the applicability of new approaches, such as next-generation sequencing and studies utilizing functional markers to gain insights into more contemporary processes shaping host-parasite relationships.

A case of canine Angiostrongylus vasorum in Scotland confirmed by PCR and sequence analysis

This report describes the first recognised case of Angiostrongylus vasorum infection in a native Scottish dog and its diagnostic confirmation using polymerase chain reaction. The use of this technique to definitively diagnose canine angiostrongylosis confirms the presence of this epidemiologically significant case. The implications of the changing distribution of this parasite are discussed.

North American paragonimiasis (Caused by Paragonimus kellicotti) in the context of global paragonimiasis

Paragonimus species are highly evolved parasites with a complex life cycle that involves at least three different hosts, i.e., snails, crustaceans, and mammals. The adult forms of Paragonimus species reside and mate in the lungs of a variety of permissive mammalian hosts, including humans. Although human paragonimiasis is uncommonly encountered in North America, both autochthonous and imported disease may be encountered. Paragonimus kellicotti, the species endemic to North America, is a well-known pathogen in wild and domestic animals. Five patients with North American paragonimiasis have been reported in the recent medical literature. The biologic, clinical, radiologic, and laboratory features of paragonimiasis are reviewed, with emphasis on North American paragonimiasis whenever possible.

Internal repetition and intraindividual variation in the rDNA ITS1 of the anopheles punctulatus group (Diptera: Culicidae): multiple units and rates of turnover

The rapid divergence of repetitive sequences makes them desirable markers for phylogenetic studies of closely related groups, provided that a high level of sequence homogeneity has been maintained within species. Intraspecific polymorphisms are found in an increasing number of studies now, and this highlights the need to determine why these occur. In this study we examined intraindividual variation present in the first ribosomal internal transcribed spacer (ITS1) from a group of cryptic mosquito species. Individuals of the Anopheles punctulatus group contained multiple ITS1 length variants that ranged from 1.2 to 8.0 kb. Nucleotide and copy number variation for several homologous internal repeats is common, yet the intraspecific sequence divergence of cloned PCR isolates is comparable to that of other mosquito species (~0.2-1.5%). Most of the length variation is comprised of a 5'-ITS1 repeat that was identified as a duplication of a conserved ITS2 region. Secondary structure conservation for this repeat is pronounced and several repeat types that are highly homogenized have formed. Significant interspecific divergence indicates a high rate of evolutionary change for this spacer. A maximum likelihood tree constructed here was congruent with previous phylogenetic hypotheses and suggests that concerted evolution is also accompanied by interpopulation divergence. The lack of interindividual differences and the presence of homogenized internal repeats suggest that a high rate of turnover has reduced the overall level of variation. However, the intraindividual variation also appears to be maintained by the absence of a single turnover rate and the complex dynamics of ongoing recombination within the spacer.

Description of a new species of Crassicutis Manter, 1936, parasite of Cichlasoma beani Jordan (Osteichthyes: Cichlidae) in Mexico, based on morphology and sequences of the ITS1 and 28S ribosomal RNA genes

A new species of Crassicutis Manter, 1936 is described from the Sinaloan cichlid Cichlasoma beani (Jordan) (Osteichthyes: Cichlidae) in the upper Río Santiago basin. Crassicutis choudhuryi n. sp. differs from most of the other nominal species by having testes located in a symmetrical position. The only other species of the genus that includes some specimens exhibiting this trait is Crassicutis intermedius (Szidat 1954), a species found in 5 species of siluriforms and 1 species of characiform in South America. However, this species differs from Cr. choudhuryi n. sp. by having testes almost half of the size, and vitelline follicles extending anteriorly to the region between the acetabulum and the intestinal bifurcation. The new species is morphologically very similar to Crassicutis cichlasomae Manter, 1936, but clearly differs from this species because of the constantly symmetrical position of the testes. Additionally, Cr. choudhuryi n. sp. is found in the Santiago River basin on the Pacific slope of Mexico, parasitizing specifically the endemic Ci. beani that does not co-occur with any other cichlid. Cr. cichlasomae exhibits more hosts (about 25 species of cichlids only in Mexico) and a wider distribution range that extends from northeastern Mexico southward to Central America, Cuba, and Brazil. To corroborate that our specimens were not conspecific with Cr. cichlasomae, sequences of the internal transcribed spacer (ITS1) and the 28S ribosomal RNA genes of individuals from several populations (recently collected in southeastern Mexico) were obtained and compared to the species described herein. Sequence divergence (1.3% for the 28S and 4.0% for the ITS1) gives further support to the erection of a new species.

Intraspecific concerted evolution of the rDNA ITS1 in Anopheles farauti sensu stricto (Diptera: Culicidae) reveals recent patterns of population structure

We examined the intraindividual variation present in the first ribosomal internal transcribed spacer (ITS1) of Anopheles farauti to determine the level of divergence among populations for this important malarial vector. We isolated 187 clones from 70 individuals and found regional variation among four internal tandem repeats. The data were partitioned prior to analysis given the presence of a paralogous ITS2 sequence, called the 5'-subrepeat, inserted in the ITS1 of most clones. A high level of homogenization and population differentiation was observed for this repeat, which indicates a higher rate of turnover relative to the adjacent 'core' region. Bayesian analysis was performed using several substitutional models on both a combined and a partitioned data set. On the whole, the ITS1 phylogeny and geographic origin of the samples appear to be congruent. Some interesting exceptions indicate the spread of variant repeats between populations and the retention of ancestral polymorphism. Our data clearly demonstrate concerted evolution at the intraspecific level despite intraindividual variation and a complex internal repeat structure from a species that occupies a continuous coastal distribution. A high rate of genomic turnover in combination with a high level of sequence divergence appears to be a major factor leading to its concerted evolution within these populations.

Structural evolution of nrDNA ITS in Pinaceae and its phylogenetic implications

Nuclear ribosomal DNA (nrDNA) has been considered as an important tool for inferring phylogenetic relationships at many taxonomic levels. In comparison with its fast concerted evolution in angiosperms, nrDNA is symbolized by slow concerted evolution and substantial ITS region length variation in gymnosperms, particularly in Pinaceae. Here we studied structure characteristics, including subrepeat composition, size, GC content and secondary structure, of nrDNA ITS regions of all Pinaceae genera. The results showed that the ITS regions of all taxa studied contained subrepeat units, ranging from 2 to 9 in number, and these units could be divided into two types, longer subrepeat (LSR) without the motif (5'-GGCCACCCTAGTC) and shorter subrepeat (SSR) with the motif. Phylogenetic analyses indicate that the homology of some SSRs still can be recognized, providing important informations for the evolutionary history of nrDNA ITS and phylogeny of Pinaceae. In particular, the adjacent tandem SSRs are not more closely related to one another than they are to remote SSRs in some genera, which may imply that multiple structure variations such as recombination have occurred in the ITS1 region of these groups. This study also found that GC content in the ITS1 region is relevant to its sequence length and subrepeat number, and could provide some phylogenetic information, especially supporting the close relationships among Picea, Pinus, and Cathaya. Moreover, several characteristics of the secondary structure of Pinaceae ITS1 were found as follows: (1) the structure is dominated by several extended hairpins; (2) the configuration complexity is positively correlated with subrepeat number; (3) paired subrepeats often partially overlap at the conserved motif (5'-GGCCACCCTAGTC), and form a long stem, while other subrepeats fold onto itself, leaving part of the conserved motif exposed in hairpin loops.

Coevolution of Retrovarium n. gen. (Digenea: Cryptogonimidae) in Lutjanidae and Haemulidae (Perciformes) in the Indo-West Pacific

We describe 11 new species of cryptogonimids belonging to Retrovarium n. gen., from eight species of Lutjanidae and one species of Haemulidae, from the Great Barrier Reef, French Polynesia and the Maldives. We also transfer Neoparacryptogonimus saccatus (Manter, 1963) and Neoparacryptogonimus sphericus Nahhas, Sey & Nishimoto, 1998 to Retrovarium. The morphologically based taxonomic approach was augmented with DNA sequence data from three nuclear ribosomal DNA regions (28S, ITS1 and ITS2) to explore the species integrity, biogeographic distribution and evolution of the species recognised here. Sequencing included multiple replicates and revealed 11 distinct genotypes which corroborated our morphologically based hypotheses of putative species present in the system. There was no intraspecific variation and all three rDNA regions differed between every combination of species. Two species exhibited wide geographic ranges, having identical rDNA sequences between the Great Barrier Reef and the Maldives, localities separated by over 9600 km. One host species, Symphorus nematophorus, proved to be exceptionally rich, harbouring six species. Minimum evolution analyses were conducted on each of the rDNA datasets independently; minimum evolution, maximum likelihood and Bayesian inference analyses were conducted on a combined (28S, ITS1 and ITS2) dataset for sequence comparison purposes and to explore the evolutionary history of these parasites. To examine the coevolutionary history of this complex, assessment of phylogenetic relationships between the 23 species of Lutjanidae and two species of Haemulidae collected during this survey was performed with data from 16S and cytochrome b mtDNA using Bayesian inference analysis. Despite the high host specificity observed in most of the species, mapping of the parasites on the host phylogeny revealed an absence of strict coevolution or co-descent within this complex. Overall, Retrovarium appears to have had an exceptionally patchy radiation, failing to infect many taxa, infecting species with no readily discernible pattern, and radiating dramatically within one species.

Cryptic long internal repeat sequences in the ribosomal DNA ITS1 gene of the dinoflagellate Cochlodinium polykrikoides (dinophyceae): a 101 nucleotide six-repeat track with a palindrome-like structure

Extremely long PCR fragments were generated by PCR amplification of ITS and 5.8S rDNA from Cochlodinium polykrikoides against other dinoflagellates. These patterns were consistent among geographically different isolates of C. polykrikoies. DNA sequencing reactions revealed that the PCR products were 1,166 bp in length and consisted of 813 bp of ITS1, 160 bp of 5.8S rDNA and 193 bp of ITS2. Thus, the long length was caused mainly by the long ITS1 sequence. Cryptically, the ITS1 contained a tract of 101 bp that occurs six times in tandem. The six repeated elements had identical nucleotide sequences. ITS1, therefore, separated three distinct regions: the 5' end (122 bp), the six parallel repeats (606 bp), and the 3' region (85 bp). Interestingly, both the single and six-repeat sequences should be palindrome-like sequences. In inferred secondary structures, both repeat sequences formed a long helical structure. This is the first reported discovery of comparatively long internal repeats in the ITS1 of dinoflagellates.