China is known for the endemic diversity of Paragonimus species wherein 38 out of the 50 globally reported species have been identified (15, 17). Hence, China is a hotspot for phylogenetic diversity among the previously reported species. The probability of phylogenetic diversification is more in naturally infected secondary hosts like crabs, wherein other ecological factors may drive diversification (18). Hence, we determined the occurrence and phylogenetic diversity of Paragonimus species in Yunnan province China. In the current study, out of the sequenced samples, 2 were P.microrchis and 1 each of P.heterotremus, P. skrjabini and P.proliferus.
Since the differential identification features of cercariae are not well marked, most often the species identification is done at the metacercariae stage from the secondary intermediate hosts (3). In this study morphological identification made at the metacercariae stage and adult stage were concordant with each other further substantiating the role of morphological identification. This was also confirmed with PCR amplification and sequencing of CO1 and ITS2 region.
Phylogenetic analysis with previously published sequences and the sequences generated in this study revealed concordance between the single locus (CO1 or ITS2) and double loci (CO1 and ITS2) phylogenetic analysis. The phylogenetic analysis revealed the presence of 3 distinct clusters of Paragonimus consisting of P.siemensis, P. westermani in cluster 1, P. skrjabini, P. miyazakii, P. proliferus, P. macrorchis and P. heterotremus in cluster 2 and P. microrchis, P. bangkokensis, P. harinasutai, P. iloktsuensis and P. sadoensis in cluster 3. Extensive branching was observed in cluster 2 wherein P.skrjabini isolated from different regions in China showed as distinct branches. Further, P.heterotremus clustered into a stable branch with P. heterotremus from south east Asia clustering into a single monophyletic branch group distinct from the P.heterotremus from India. Similarly, P.proliferus also clustered as a monophyletic group in most cases (except in NJ trees constructed by ITS2). The relationship between P. macrorchis and the first two branches is not clear, but in the MP phylogenetic tree of ITS2 data sets, COI and ITS2 combined data sets, P.heterotremus tends to be clustered into one branch.
By comparing the results of constructing adjacent NJ and MP phylogenetic trees from COI, ITS2, COI and ITS2 datasets, it is obvious that different data sets and tree-building methods have great influence on the robustness of phylogenetic trees. In phylogenetic trees constructed from single gene datasets, the number of branches with more than 75% support is less than that constructed from joint datasets, especially self-exhibition, which deserves significant improvement. The phylogenetic tree constructed by MP method using COI and ITS2 data sets has a high robustness, and P.macrorchis and P.heterotremus were clustered together, which is supported by 91%.
Among the Paragonimus spp. reported in this study, all the species were previously reported from China. A study by Lou et al., utilized similar approaches in identifying the phylogenetic diversity of Paragonimus in China and reported clustering of isolates based on geographic location which was in accordance to the current study(14). P.skrjabini was previously reported from Gansu, Shanxi, Yunnan, Guangxi, Guizhou, Sichuan and Jiangxi provinces of China. While, P heterotremus has been reported only from Yunnan and Guangxi provinces(2).
Phylogenetic analysis has varied values in different organisms. In Paragonimus, where there is no consensus on speciation with both morphological and molecular characteristics, phylogenetic analysis helps in determining the phylogeographies(7). The results of our study further substantiated the conservative phylogeographies revealed through analysis of Cox1 and ITS 2 regions.
The study is not without limitations. The study sites were localized to Yunnan province in China; hence the results may not be extended to other parts of China. Secondly, we used maximum parsimony method which is infrequently used despite the robustness mainly because of the nature of assumptions deployed in the method. Nevertheless, maximum likelihood method which was also used in other studies also has drawbacks.
In conclusion, a total of 4 morphologically different species of Paragonimus was isolated from Yunnan province, China. Phylogenetic analysis using the polymorphisms in two different loci revealed considerable variations in the species of Paragonimus identified in different geographic locations with clustering of Paragonimus based on geographic location.