Sequence analysis of Paragonimus internal transcribed spacer 2 and cyclooxygenase 1 genes in freshwater crabs in Henan Province

Abstract

To investigate the sequences of internal transcribed spacer 2 (ITS2) and cyclooxygenase 1 (COX1) genes of Paragonimus metacercariae in freshwater crabs in Henan Province, identify the species of Paragonimus and evaluate its genetic relationships with Paragonimus isolates from other provinces in China. Freshwater crabs were collected from 8 survey sites in Zhengzhou, Luoyang, Pingdingshan, Nanyang and Jiyuan cities of Henan Province from 2016 to 2021, and Paragonimus metacercariae were detected in freshwater crabs. Genomic DNA was extracted from Paragonimus metacercariae, and the ITS2 and COX1 genes were amplified using PCR assay, followed by sequencing of PCR amplification products. The gene sequences were spliced and aligned using the software DNASTAR, and aligned with the sequences of Paragonimus genes in the GenBank. Phylogenetic trees were created using the MEGA6 software with the Neighbor-Joining method based on ITS2 and COX1 gene sequences, with Fasciola hepatica as the outgroup. The detection rates of Paragonimus metacercariae were 6.83% (11/161), 50.82% (31/61), 18.52% (5/26), 8.76% (12/137), 14.29% (9/63), 17.76% (19/105), 18.50% (32/173) and 42.71% (41/96) in freshwater crabs from 8 survey sites in Zhengzhou, Luoyang, Pingdingshan, Nanyang and Jiyuan cities of Henan Province, with a mean detection rate of 19.46% (160/822), and a mean infection intensity of 0.57 metacercariae/g. The amplified ITS2 and COX1 gene fragments of Paragonimus were approximately 500 bp and 450 bp in lengths, respectively. The ITS2 gene sequences of Paragonimus metacercariae from 8 survey sites of Henan Province showed the highest homology (99.8% to 100.0%) with the gene sequence of P. skrjabini (GenBank accession number: MW960209.1), and phylogenetic analysis showed that the Paragonimus in this study was clustered into the same clade with P. skrjabini from Sichuan Province (GenBank accession number: AY618747.1), Guangxi Zhuang Autonomous Region (GenBank accession number: AY618729.1) and Hubei Province (GenBank accession number: AY618751.1), and P. miyazaki from Fujian Province (GenBank accession number: AY618741.1) and Japan (GenBank accession number: AB713405.1). The COX1 gene sequences of Paragonimus metacercariae from 8 survey sites of Henan Province showed the highest homology (90.0% to 100.0%) with the gene sequence of P. skrjabini (GenBank accession number: AY618798.1), and phylogenetic analysis showed that the Paragonimus in this study was clustered into the same clade with all P. skrjabini and clustered into the same sub-clade with P. skrjabini from Hubei Province (GenBank accession numbers: AY618782.1 and AY618764.1). Paragonimus species from freshwater crabs in Henan Province were all characterized as P. skrjabini, and the ITS2 and COX1 gene sequences had the highest homology to those of P. skrjabini from Hubei Province. The results provide insights into study of Paragonimus in Henan Province and China.