Genetic variation and population structure of Fasciola hepatica: an in silico analysis

Fasciola hepatica is a trematode leading to heavy economic setbacks to the livestock sector globally. The population's genetic information and intimate kinship level are frequently assessed using analysis of mitochondrial DNA. In this analysis, we retrieved cox1 (n = 247) and nad1 (n = 357) sequences of F. hepatica from the NCBI GenBank database and aligned the sequences with the respective reference sequences using MEGA software. The median joining network was drawn using PopArt software while neutrality and diversity indices were estimated with the help of DnaSp software. Neighbor-joining phylogenetic tree was constructed using the MEGA software package. A total of 46 and 98 distinctive haplotypes were observed for cox1 and nad1 genes, respectively. Diversity indices indicated high haplotype and nucleotide diversities in both genes. Positive Tajima's D and Fu's Fs values were found for the entire population of both the genes under study. The cox1 and nad1 gene segments in this study showed high Tajima's D values, suggesting a low likelihood of future population growth. The Tajima's D value of the nad1 gene sequence is lower (2.14910) than that of the cox1 gene sequence (3.40314), which suggests that the former is growing at a slower rate. However, the region-wise analysis revealed that both the cox1 and nad1 genes showed deviation from neutrality suggesting a recent population expansion as a result of an excess of low-frequency polymorphism. Furthermore, the overall host-wise analysis showed positive and significant Tajima's D values for the cox1 and nad1 gene sequences. To the best of our knowledge, this is the first attempt to provide insights into genetic variations and population structure of F. hepatica at a global scale using cox1 and nad1 genes. Our findings suggest the existence of specific variants of F. hepatica in different parts of the world and provide information on the molecular ecology of F. hepatica. The results of this study also mark a critical development in upcoming epidemiological investigations on F. hepatica and will also contribute to understanding the global molecular epidemiology and population structure of F. hepatica.