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.

Phylogenotyping of infectious bursal disease virus in Vietnam according to the newly unified genotypic classification scheme.

Since 1987, infectious bursal disease virus (IBDV) has circulated and evolved in Vietnam, but little is known about the genotypes present. IBDV samples were collected in 1987, 2001-2006, 2008, 2011, 2015-2019, and 2021 in 18 provinces. We conducted phylogenotyping analysis based on an alignment of 143 VP2-HVR (hypervariable region) sequences from 64 Vietnamese isolates (26 previous and 38 additional sequences and two vaccines, and alignment of 82 VP1 B-marker sequences, including one vaccine and four Vietnamese field strains. The analysis identified three A-genotypes, A1, A3, and A7, and two B-genotypes, B1 and B3, among the Vietnamese IBDV isolates. The lowest average evolutionary distance (8.6%) was seen between the A1 and A3 genotypes, and the highest (21.7%) was between A5 and A7, while there was a distance of 14% between B1 and B3 and 17% between B3 and B2. Unique signature residues were observed for genotypes A2, A3, A5, A6, and A8, which could be used for genotypic discrimination. A timeline statistical summary revealed that the A3-genotype predominated (79.8% presence) in Vietnam from 1987 to 2021 and that it remained the dominant IBDV genotype over the last five years (2016-2021). The current study contributes to a better understanding of the circulating genotypes and evolution of IBDV in Vietnam and worldwide.

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.

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.

Molecular genotypic analysis of porcine circovirus type 2 reveals the predominance of PCV2d in Vietnam (2018-2020) and the association between PCV2h, the recombinant forms, and Vietnamese vaccines.

We conducted nucleotide and amino acid sequence alignment and phylogenetic analysis of porcine circovirus ORF2 (Cap protein) from 17 PCV2-positive clinical samples from nine different northern Vietnamese provinces (Mar 2018-Nov 2020), four local vaccines, and 77 reference strains. We identified one PCV2a (1/17 = 5.9%), five PCV2b (5/17 = 29.9%), and 11 PCV2d (11/17 = 64.7%) isolates, while only PCV2d was detected in 2020. Timeline analysis indicated an increasing predominance of PCV2d nationwide (2018-2020). With strong nodal support (98% for nucleotides and 74% for amino acids), the phylogenetic tree topology revealed a distinct PCV2h clade including recombinant/intermediate strains and local vaccines. The Cap protein sequences from 11 PCV2d field strains had the 2d-genotype-typical motif 86SNPLSV91 in loop CD, the motif TGID in loop GH-HI, and the motif 230PLNPK234 in loop CT. The PCV2h isolates (and vaccines) had the 86SNPLSV91, SAID, and 230L(N/H)PK234 motifs. Selection pressure analysis indicated positive selection at seven sites: A68N in immunoreactive region (IRR)-A; 119G and 130V in IRR-B; and 167L, T190(A/S), 194D and 202F in IRR-C. We identified PCV2h as the genotype of the recombinant strains, which resulted from intergenotype recombination of PCV2a, PCV2b, and PCV2d. The current data provide new information about the diversity, distribution, and dominance of the PCV2 genotype in Vietnam.

Canine parvovirus type 2c in Vietnam continues to produce distinct descendants with new mutations restricted to Vietnamese variants.

Viral protein 2 (VP2) of canine parvovirus (CPV) exhibits a high degree of genetic and antigenic diversity. We analyzed 88 Vietnamese CPV-VP2 sequences (1755 bp), 34 from this study and 54 from previous studies, and discovered a new sublineage, "new var.", within the lineage CPV-2c-"new", characterized by the mutation 5G/447M, which is restricted to the Vietnamese isolates. These new mutants appear to have emerged in recent years, accounting for 65.5% of the total. With strong nodal support (98%), the distinct Vietnamese 2c-"new-var." sublineage (5G/426E/447M) was found to be separate from the 2c-"new" sublineage (5G/426E/447I) within the 2c-(Asia)/Asia-2c lineage. Amino acid changes in epitopes of VP2 might have led to the generation of subvariants and affected the antigenicity, immunogenicity, or virulence of the virus, resulting in vaccine failure worldwide.

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.

The ribosomal transcription units of Haplorchis pumilio and H. taichui and the use of 28S rDNA sequences for phylogenetic identification of common heterophyids in Vietnam.

BACKGROUND: Heterophyidiasis is now a major public health threat in many tropical countries. Species in the trematode family Heterophyidae infecting humans include Centrocestus formosanus, Haplorchis pumilio, H. taichui, H. yokogawai, Procerovum varium and Stellantchasmus falcatus. For molecular phylogenetic and systematic studies on trematodes, we need more prospective markers for taxonomic identification and classification. This study provides near-complete ribosomal transcription units (rTU) from Haplorchis pumilio and H. taichui and demonstrates the use of 28S rDNA sequences for identification and phylogenetic analysis. RESULTS: The near-complete ribosomal transcription units (rTU), consisting of 18S, ITS1, 5.8S, ITS2 and 28S rRNA genes and spacers, from H. pumilio and H. taichui from human hosts in Vietnam, were determined and annotated. Sequence analysis revealed tandem repetitive elements in ITS1 in H. pumilio and in ITS2 in H. taichui. A phylogenetic tree inferred from 28S rDNA sequences of 40 trematode strains/species, including 14 Vietnamese heterophyid individuals, clearly confirmed the status of each of the Vietnamese species: Centrocestus formosanus, Haplorchis pumilio, H. taichui, H. yokogawai, Procerovum varium and Stellantchasmus falcatus. However, the family Heterophyidae was clearly not monophyletic, with some genera apparently allied with other families within the superfamily Opisthorchioidea (i.e. Cryptogonimidae and Opisthorchiidae). These families and their constituent genera require substantial re-evaluation using a combination of morphological and molecular data. Our new molecular data will assist in such studies. CONCLUSIONS: The 28S rDNA sequences are conserved among individuals within a species but varied between genera. Based on analysis of 40 28S rDNA sequences representing 19 species in the superfamily Opisthorchioidea and an outgroup taxon (Alaria alata, family Diplostomidae), six common human pathogenic heterophyids were identified and clearly resolved. The phylogenetic tree inferred from these sequences again confirmed anomalies in molecular placement of some members of the family Heterophyidae and demonstrates the need for reappraisal of the entire superfamily Opisthorchioidea. The new sequences provided here supplement those already available in public databases and add to the array of molecular tools that can be used for the diagnosis of heterophyid species in human and animal infections.

Molecular genotyping of duck hepatitis A viruses (DHAV) in Vietnam.

INTRODUCTION: The aim of this study was to identify the genetic characteristics and molecular genotyping of duck hepatitis A virus (DHAV) isolated in Vietnam during 2009-2013. METHODOLOGY: Thirty duckling livers from outbreaks between 2009 and 2013 in seven provinces were collected and identified by polymerase chain reaction (PCR). Then, VP1 genes of eleven positive samples and two attenuated vaccine strains were sequenced and analyzed. RESULTS: Genotypic and phylogenetic analyses indicated that the 13 Vietnamese isolates were classified into two genotypes, DHAV-1 and DHAV-3. The rate of identity and homology was 91%-100% between the 10 Vietnamese and 26 global strains of DHAV-3, and 92%-100% between 3 Vietnamese and 16 strains of DHAV-1. Between the DHAV-3 and DHAV-1 strains, the divergence reached 30%. At the C-terminal of VP1 for the different strains, a hypervariable region was observed, and notably, six of the Vietnamese DHAV-3 strains in this study showed four consistent differences (at positions T184M, Q200H, K207N, and K214R) within this group that were distinct from all other DHAV-3 strains. CONCLUSIONS: This is the first report of molecular characterization of DHAVs in Vietnam. At least two genotypes were identified, DHAV-1 and DHAV-3, with diversified clades within and between genotypes. DHAV-3 seemed to be dominant in Vietnam.

A complete mitochondrial genome from Echinochasmus japonicus supports the elevation of Echinochasminae Odhner, 1910 to family rank (Trematoda: Platyhelminthes).

The complete mitochondrial (mt) genome of the trematode Echinochasmus japonicus Tanabe, 1926 was fully determined and annotated. The circular mt molecule of this species is 15,865bp in length, containing 12 protein-coding genes (arranged in the following order: cox3-cob-nad4L-nad4-atp6-nad2-nad1-nad3-cox1-cox2-nad6-nad5), two ribosomal RNA genes (rrnL and rrnS) and 22 transfer RNA genes (trnH; trnQ; trnF; trnM; trnV; trnA; trnD; trnN; trnP; trnI; trnK; trnS1(AGN); trnW; trnT; trnC; trnL1(CUN); trnS2(UCN); trnL2(UUN); trnG; and trnE). The atp8 gene is absent. The 3' end of nad4L overlaps the 5' end of nad4 by 40bp. An array of eight identical tandem repeats (240bp each) was found between trnG and trnE in the long non-coding region of the individual worm sequenced. Numbers of these repeats varied among E. japonicus samples. Phylogenetic analysis of concatenated mt protein sequences of 40 trematode species/strains supports the elevation of Echinochasminae Odhner, 1910 to family rank, close to the families Echinostomatidae and Fasciolidae. As echinochasmid and echinostomatid species can parasitize humans, the future characterization of additional mt genomes is needed for development of mt markers for identification and phylogenetic, population, epidemiological and hybridization studies.

Toxocara malaysiensis infection in domestic cats in Vietnam--An emerging zoonotic issue?

Toxocara canis of canids is a parasitic nematode (ascaridoid) that infects humans and other hosts, causing different forms of toxocariasis. This species of Toxocara appears to be the most important cause of human disease, likely followed by Toxocara cati from felids. Although some studies from Malaysia and China have shown that cats can harbor another congener, T. malaysiensis, no information is available about this parasite for other countries. Moreover, the zoonotic potential of this parasite is unknown at this point. In the present study, we conducted the first investigation of domestic dogs and cats for Toxocara in Vietnam using molecular tools. Toxocara malaysiensis was identified as a common ascaridoid of domestic cats (in the absence of T. cati), and T. canis was commonly found in dogs. Together with findings from previous studies, the present results emphasize the need to explore the significance and zoonotic potential of T. malaysiensis in Vietnam and other countries where this parasite is endemic and prevalent in cats.

Development and evaluation of a single-step duplex PCR for simultaneous detection of Fasciola hepatica and Fasciola gigantica (family Fasciolidae, class Trematoda, phylum Platyhelminthes).

A single-step multiplex PCR (here referred to as a duplex PCR) has been developed for simultaneous detection and diagnosis of Fasciola hepatica and F. gigantica. These species overlap in distribution in many countries of North and East Africa and Central and Southeast Asia and are similar in egg morphology, making identification from fecal samples difficult. Based on a comparative alignment of mitochondrial DNA (mtDNA) spanning the region of cox1-trnT-rrnL, two species-specific forward primers were designed, FHF (for F. hepatica) and FGF (for F. gigantica), and a single reverse primer, FHGR (common for both species). Conventional PCR followed by sequencing was applied using species-specific primer pairs to verify the specificity of primers and the identity of Fasciola DNA templates. Duplex PCR (using three primers) was used for testing with the DNA extracted from adult worms, miracidia, and eggs, producing amplicons of 1,031 bp for F. hepatica and 615 bp for F. gigantica. The duplex PCR failed to amplify from DNA of other common liver and intestinal trematodes, including two opisthorchiids, three heterophyids, an echinostomid, another fasciolid, and a taeniid cestode. The sensitivity assay showed that the duplex PCR limit of detection for each Fasciola species was between 0.012 ng and 0.006 ng DNA. Evaluation using DNA templates from 32 Fasciola samples (28 adults and 4 eggs) and from 25 field-collected stools of ruminants and humans revealed specific bands of the correct size and the presence of Fasciola species. This novel mtDNA duplex PCR is a sensitive and fast tool for accurate identification of Fasciola species in areas of distributional and zonal overlap.