Molecular insights into the population structure and haplotype network of Theileria annulata based on the small-subunit ribosomal RNA (18S rRNA) gene

The present study was conducted to elucidate the population genetic diversity and haplotype network of Theileria annulata based on all available nearly complete 18S rRNA gene sequences in the GenBank™. In total, 52 sequences of the nuclear 18S rRNA gene used to assess the relationship of T. annulata with their country of origin identified 34 haplotypes. Haplotype 4 was widespread, occurring in India, China, Turkey and Iran, while the remaining haplotypes were singleton and unique to one country. Haplotype 4 displayed numerous single haplotypes around it and the stellate shape of the network suggested a rapid population expansion. India exhibited the largest number of haplotypes (h = 25) followed by Turkey (h = 6), China (h = 4), and Iran and Italy (h = 1). No geographical clustering of haplotypes was recorded. Nucleotide diversity was the highest in the Turkish followed by the Indian and Chinese populations. Similarly, haplotype diversity was the highest in China followed by Turkey, and the lowest in India. The overall dataset exhibited a low nucleotide diversity (0.00253 ± 0.00035), but high haplotype diversity (0.917 ± 0.034). It suggested the presence of only minor differences (01–11 nucleotide) between haplotypes which was also evident from the haplotype network. A high level of genetic diversity was documented within the Indian, Chinese and Turkish populations of T. annulata, whereas little genetic differentiation was noticed among these populations with a very high level of gene flow. Analysis of molecular variance (AMOVA) of T. annulata sequences revealed higher genetic variation within countries (83.58%) as compared to the variation among countries (16.42%). Neutrality indices, viz., Tajima's D, Fu and Li's F, Fu's Fs, and R2, along with the unimodal mismatch distributions demonstrated a recent population expansion of T. annulata in India and the overall dataset. However, the non-significant values of Tajima's D, Fu and Li's F, and Fu's Fs for the Chinese population along with a bimodal mismatch distribution signified a constant population size. For the Turkish population, the neutrality and mismatch distribution tests either indicated a constant or a slight increase in population size. The present study provides novel insights into the population genetics and haplotype network of T. annulata based on the 18S rRNA gene for the first time.