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.

Mitophylogeny of Pangasiid Catfishes and its Taxonomic Implications for Pangasiidae and the Suborder Siluroidei.

Pangasiidae (catfish order: Siluriformes) comprises 30 valid catfish species in four genera: Pangasius, Pangasianodon, Helicophagus, and Pseudolais. Their systematics are frequently revised due to the addition of newly described species. Although Pangasiidae is known to be a monophyletic family, the generic and phylogenetic relationships among the taxa are poorly resolved. This study characterized three newly obtained complete mitogenomes of Mekong River catfishes from Vietnam (Pangasius mekongensis, Pangasius krempfi, and Pangasianodon hypophthalmus), as well as the inter-and intrafamilial relationships of the Pangasiidae and catfish families in Siluroidei. The genomic features of their mitogenomes were similar to those of previously reported pangasiids, including all regulatory elements, extended terminal associated sequences (ETAS), and conserved sequence blocks (CSBs) (CSB-1, CSB-2, CSB-3, and CSBs, A to F) in the control region. A comprehensive phylogeny constructed from datasets of multiple 13 PCG sequences from 117 complete mitogenomes of 32 recognized siluriform families established Pangasiidae as monophyletic and a sister group of Austroglanididae. The [Pangasiidae + Austroglanididae] + (Ictaluridae + Cranoglanididae) + Ariidae] clade is a sister to the "Big Africa" major clade of Siluriformes. Furthermore, both phylogenies constructed from the single barcodes (83 partial cox1 and 80 partial cytB, respectively) clearly indicate genus relationships within Pangasiidae. Pangasianodon was monophyletic and a sister to the (Pangasius + Helicophagus + Pseudolais) group. Within the genus Pangasius, P. mekongensis was placed as a sister taxon to P. pangasius. Pangasius sanitwongsei was found to be related to and grouped with Pangasianodon, but in single-gene phylogenies, it was assigned to the Pangasius + Helicophagus + Pseudolais group. The datasets in this study are useful for studying pangasiid systematics, taxonomy and evolution.

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.

Multiple genetic lineages of anadromous migratory Mekong catfish Pangasius krempfi revealed by mtDNA control region and cytochrome b.

Population genetic structure of migratory fishes can reflect ecological and evolutionary processes. Pangasius krempfi is a critically important anadromous catfish in the Mekong River, and its migration pathways and genetic structure have attracted much interest. To investigate, we quantified the genetic diversity of this species using the control region (D-loop) and Cytochrome b (Cytb) of the mitochondrial genome. Fish were sampled (n = 91) along the Mekong tributaries from upstream to estuaries and coastal areas in the Mekong Delta and compared to three samples from Pakse (Laos). The D-loop haplotype (0.941 ± 0.014) and nucleotide diversity (0.0083 ± 0.0005) were high in all populations, but that of Cytb was low (0.331 ± 0.059 and 0.00063 ± 0.00011, respectively). No genetic difference was detected between populations, indicating strong gene flow and confirming a long migration distance for this species. Pangasius krempfi was not genetically structured according to geographical populations but was delineated into three haplogroups, suggesting multiple genetic lineages. The presence of haplogroups in each sampling location implies that migration downstream is random but parallel when the fish enter two river tributaries bifurcating from the main Mekong River. Individuals can also migrate along the coast, far from the estuaries, suggesting a longer migration path than previously reported, which is crucial for maintaining diverse genetic origin and migration pathways for P. krempfi.

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 analysis reveals a distinct subgenogroup of porcine epidemic diarrhea virus in northern Vietnam in 2018-2019.

The spike protein (S) of porcine epidemic diarrhea virus (PEDV), in particular, the C-terminal domain of the S1 subunit (S1-CTD), which contains the conserved CO-26K-equivalent (COE) region (aa 499-638), which is recognized by neutralizing antibodies, exhibits a high degree of genetic and antigenic diversity. We analyzed 61 PEDV S1-CTD sequences (630 nt), including 26 from samples collected from seven provinces in northern Vietnam from 2018 to 2019 and 35 other sequences, representing the G1a and 1b, G2a and 2b, and recombinant (G1c) genotypes and vaccines. The majority (73.1%) of the strains (19/26) belonged to subgroup G2b. In a phylogenetic analysis, seven strains were clustered into an independent, distinct subgenogroup named dsG with strong nodal support (98%), separate from both G1a and G1b as well as G2a, 2b, and G1c. Sequence analysis revealed distinct changes (513T>S, 520G>D, 527V>(L/M), 591L>F, 669A>(S/P), and 691V>I) in the COE and S1D regions that were only identified in these Vietnamese strains. This cluster is a new antigenic variant subgroup, and further studies are required to investigate the antigenicity of these variants. The results of this study demonstrated the continuous evolution in the S1 region of Vietnamese PEDV strains, which emphasizes the need for frequent updates of vaccines for effective protection.

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.

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 tRNASer1(AGN) and tRNASer2(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.

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.

Comparative genomics and transcriptomics of 4 Paragonimus species provide insights into lung fluke parasitism and pathogenesis.

BACKGROUND: Paragonimus spp. (lung flukes) are among the most injurious foodborne helminths, infecting ∼23 million people and subjecting ∼292 million to infection risk. Paragonimiasis is acquired from infected undercooked crustaceans and primarily affects the lungs but often causes lesions elsewhere including the brain. The disease is easily mistaken for tuberculosis owing to similar pulmonary symptoms, and accordingly, diagnostics are in demand. RESULTS: We assembled, annotated, and compared draft genomes of 4 prevalent and distinct Paragonimus species: Paragonimus miyazakii, Paragonimus westermani, Paragonimus kellicotti, and Paragonimus heterotremus. Genomes ranged from 697 to 923 Mb, included 12,072-12,853 genes, and were 71.6-90.1% complete according to BUSCO. Orthologous group analysis spanning 21 species (lung, liver, and blood flukes, additional platyhelminths, and hosts) provided insights into lung fluke biology. We identified 256 lung fluke-specific and conserved orthologous groups with consistent transcriptional adult-stage Paragonimus expression profiles and enriched for iron acquisition, immune modulation, and other parasite functions. Previously identified Paragonimus diagnostic antigens were matched to genes, providing an opportunity to optimize and ensure pan-Paragonimus reactivity for diagnostic assays. CONCLUSIONS: This report provides advances in molecular understanding of Paragonimus and underpins future studies into the biology, evolution, and pathogenesis of Paragonimus and related foodborne flukes. We anticipate that these novel genomic and transcriptomic resources will be invaluable for future lung fluke research.

Early Postnatal and Preschool-Age Infection by Fasciola spp.: Report of Five Cases from Vietnam and Worldwide Review.

Fascioliasis is reported in five Vietnamese children aged 4 years or younger. A 10-month-old girl child and a 12-month-old boy child are the youngest patients ever diagnosed. Eggs in stools suggested an infection occurred at 5-6 months and 7-8 months of age, respectively. DNA sequencing and egg size indicated this to be the first report of a verified Fasciola gigantica infection in so small children. No specific diagnosis could be obtained in two 3-year-old children detected in the acute phase. A big and gravid ectopic F. gigantica-like worm was surgically found in a 4-year-old boy presenting with peritonitis. A worldwide review showed only 38 past cases in preschool children. They included 3, 7, 12, and 16 cases of 1, 2, 3, and 4 years, respectively, with a faster infection increase in males from 2 years onward. Reports were from all continents, except Oceania, including severe complications and death. The causal agent, when specifically diagnosed, was always Fasciola hepatica. Analyses include detection in hospital, surveys, and family outbreaks; infection sources; disease phases; parasite burden; ectopic cases; symptom onset; eosinophilia; biochemical markers; and clinical complications. C-reactive protein, creatinine, and γ-glutamyl transferase are the most useful biomarkers. A serological test and a coprological analysis are recommended for so small children, in which typical symptoms may be overlooked. Treatment problems were described with many drugs, except triclabendazole. Triclabendazole should be considered the drug of choice for such small children. The possibility of a very early infection by Fasciola spp. should be henceforth considered.

High Prevalence of Clonorchis sinensis and Other Zoonotic Trematode Metacercariae in Fish from a Local Market in Yen Bai Province, Northern Vietnam.

A small survey was performed to investigate the recent infection status of Clonorchis sinensis and other zoonotic trematode metacercariae in freshwater fish from a local market of Yen Bai city, Yen Bai province, northern Vietnam. A total of 118 fish in 7 species were examined by the artificial digestion method on March 2016. The metacercariae of 4 species of zoonotic trematodes, i.e., C. sinensis, Haplorchis pumilio, Haplorchis taichui, and Centrocestus formosanus, were detected. The metacercariae of C. sinensis were found in 62 (69.7%) out of 89 fish (5 species), and their intensity of infection was very high, 81.2 per fish infected. Prevalences of 3 intestinal flukes, H. pumilio, H. taichui and C. formosanus, were 75.0%, 47.6%, and 31.7% in positive fish species, respectively, with the metacercarial intensities of 15.5, 10.3, and 2.2 per fish infected. From the above results, it has been confirmed that various species of freshwater fish continue to play the role of the infection source of C. sinensis and other zoonotic trematodes in Yen Bai city, Yen Bai province, northern Vietnam. It is of particular note that the prevalence and intensity of C. sinensis metacercariae are much higher than those reported in previous studies in fish in northern Vietnam.

Immunization with the H5N1 Recombinant Vaccine Candidate Induces High Protection in Chickens against Vietnamese Highly Pathogenic Avian Influenza Virus Strains.

Abstract: Vietnam is one of the countries most affected worldwide by the highly pathogenic avian influenza (HPAI) virus, which caused enormous economic loss and posed threats to public health. Over nearly two decades, with the antigenic changes in the diversified H5Ny viruses, the limited protective efficacy of the available vaccines was encountered. Therefore, it is necessary to approach a technology platform for the country to accelerate vaccine production that enables quick response to new influenza subtypes. This study utilized a powerful reverse genetics technique to successfully generate a recombinant H5N1 vaccine strain (designated as IBT-RG02) containing two surface proteins (haemagglutinin (HA) and neuraminidase (NA)) from the HPAI H5N1 (A/duck/Vietnam/HT2/2014(H5N1)) of the dominant clade 2.3.2.1c in Vietnam during 2012-2014. Importantly, the IBT-RG02 vaccine candidate has elicited high antibody titres in chickens (geometric mean titre (GMT) of 6.42 and 6.92, log2 on day 14 and day 28 p.i., respectively). To test the efficacy, immunized chickens were challenged with the circulating virulent strains. As results, there was a high protection rate of 91.6% chickens against the virulent A/DK/VN/Bacninh/NCVD-17A384/2017 of the same clade and a cross-protection of 83.3% against A/duck/TG/NAVET(3)/2013 virus of clade 1.1. Our promising results showed that we can independently master the reverse genetics technology for generation of highly immunogenic vaccine candidates, and henceforth, it is a timely manner to reformulate avian influenza virus vaccines against variable H5 clade HPAI viruses in Vietnam.

Characterization and phylogenetic properties of the complete mitochondrial genome of Fascioloides jacksoni (syn. Fasciola jacksoni) support the suggested intergeneric change from Fasciola to Fascioloides (Platyhelminthes: Trematoda: Plagiorchiida).

Fascioloides jacksoni (syn. Fasciola jacksoni, Cobbold, 1869) (Platyhelminthes: Echinostomatoidea), is a liver fluke that causes severe morbidity and mortality of Asian elephants (Elephas maximus maximus). Understandings on molecular diagnosis, epidemiology, genetics and evolution of this flatworm are limited. In this study, we present the complete mitochondrial DNA (mt) sequence of 14,952 bp obtained from an individual fluke and comparative characterization of mitogenomic features with fasciolids, primarily, Fascioloides magna and other taxa in the superfamily Echinostomatoidea. Taxonomic relationship within and between Echinostomatoidea, Opisthorchioidea and Paramphistomoidea in the order Plagiorchiida, are also taxonomically considered. The complete circular mt molecule of Fas. jacksoni contained 12 protein-coding, two ribosomal RNA, 22 transfer RNA genes, and a non-coding region (NCR) rich in tandem repeat units. As common in digenean trematodes, Fas. jacksoni has the usual gene order, the absence of atp8 and the overlapped region by 40 bp between nad4L and nad4 genes. The NCR located between tRNAGlu (trnE) and cox3 contained nine nearly identical tandem repeat units (TRs of 113 bp each). Special DHU-arm missing tRNAs for Serine were found for both, tRNAS1(AGN) and tRNAS2(UCN). Base composition indicated that cox1 of Fas. jacksoni showed the lowest (11.8% to Fas. magna, 12.9 - 13.6% to Fasciola spp. and 18.1% to Fasciolopsis buski) and nad6 the highest divergence rate (19.2%, 23.8-26.5% and 27.2% to each fasciolid group), respectively. A clear bias in nucleotide composition, as of 61.68%, 62.88% and 61.54%, with a negative AT-skew of the corresponding values (-0.523, -0.225 and - 0.426) for PCGs, MRGs and mtDNA for Fas. jacksoni and likewise data for the fasciolids. Phylogenetic analysis confirmed the sister branch of Fas. jacksoni and Fas. magna with the nodal support of 100%, clearly separated from the taxonomically recognized Fasciola spp. With the previous studies, mitogenomic data presented in this study are strongly supportive for Fasciola jacksoni reappraisal as Fascioloides jacksoni in the Fascioloides genus.

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.

Adaptive Radiation of the Flukes of the Family Fasciolidae Inferred from Genome-Wide Comparisons of Key Species.

Liver and intestinal flukes of the family Fasciolidae cause zoonotic food-borne infections that impact both agriculture and human health throughout the world. Their evolutionary history and the genetic basis underlying their phenotypic and ecological diversity are not well understood. To close that knowledge gap, we compared the whole genomes of Fasciola hepatica, Fasciola gigantica, and Fasciolopsis buski and determined that the split between Fasciolopsis and Fasciola took place ∼90 Ma in the late Cretaceous period, and that between 65 and 50 Ma an intermediate host switch and a shift from intestinal to hepatic habitats occurred in the Fasciola lineage. The rapid climatic and ecological changes occurring during this period may have contributed to the adaptive radiation of these flukes. Expansion of cathepsins, fatty-acid-binding proteins, protein disulfide-isomerases, and molecular chaperones in the genus Fasciola highlights the significance of excretory-secretory proteins in these liver-dwelling flukes. Fasciola hepatica and Fasciola gigantica diverged ∼5 Ma near the Miocene-Pliocene boundary that coincides with reduced faunal exchange between Africa and Eurasia. Severe decrease in the effective population size ∼10 ka in Fasciola is consistent with a founder effect associated with its recent global spread through ruminant domestication. G-protein-coupled receptors may have key roles in adaptation of physiology and behavior to new ecological niches. This study has provided novel insights about the genome evolution of these important pathogens, has generated genomic resources to enable development of improved interventions and diagnosis, and has laid a solid foundation for genomic epidemiology to trace drug resistance and to aid surveillance.

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.

Molecular characterization of field isolates of infectious bursal disease virus from three decades, 1987-2018, reveals a distinct genotypic subgroup in Vietnam.

The complete nucleotide sequence of the viral protein 2 (VP2) ORF was determined for 26 Vietnamese infectious bursal disease isolates collected from clinical outbreaks in vaccinated flocks from 1987 to 2018 and two commercial vaccine specimens. These sequences were compared for molecular classification with 42 reference strains representing all four main classes of serotype 1, including very virulent (vvIBDV), classical (cvIBDV), antigenic variant (avIBDV) and attenuated (atIBDV) strains, and serotype 2 strains. Amino acids at nine key positions in the VP2-HVR in 20 Vietnamese isolates, A222, I242, Q253, I256, D279, A284, I294, S299, A329, which are typical of the vvIBDV class, were found to be identical in all of the isolates. Eighteen of these isolates had a unique change at residue 212 (D212N) located in the PAB loop. Phylogenetic analysis revealed a distinct lineage/subclade with strong nodal support (96%) that included recent Chinese IBDV strains that were distinct from typical vvIBDVs. Six isolates contained the amino acid substitutions P222, V242, Q253, V256, D279, A284, I294, N299, A329, which are present in two vaccine strains derived from strain 2512 and these isolates were also closely related to the classical virulent STC strain. Data from this study show that there is considerable genetic diversity among vvIBDVs, which vary according to geographic region. Antigenic drift and differences in genetic characteristics between virulent strains and IBDV vaccine strains may be the cause of vaccine failure. Better antigenic matching of vaccines to the strains circulating in Vietnam is therefore required.