Revisiting the genotypes of Theileria equi based on the V4 hypervariable region of the 18S rRNA gene

Introduction

Equine theileriosis, an economically important disease that affects horses and other equids worldwide, is caused by a tick-borne intracellular apicomplexan protozoa Theileria equi. Genotyping of T. equi based on the 18S rRNA gene revealed the presence of two, three, four or five genotypes. In previous published reports, these genotypes have been labelled either alphabetically or numerically, and there is no uniformity in naming of these genotypes. The present study was aimed to revisit the phylogeny, genetic diversity and geographical distribution of T. equi based on the nucleotide sequences of the V4 hypervariable region of the 18S rRNA gene available in the nucleotide databases.

Methods

Out of 14792 nucleotide sequences of T. equi available in the GenBank™, only 736 sequences of T. equi containing the complete V4 hypervariable region of the 18S rRNA gene (>207 bp) were used in multiple sequence alignment. Subsequently, a maximum likelihood phylogenetic tree was constructed based on the Kimura 2-parameter model (K2+I).

Results

The phylogenetic tree placed all the sequences into four distinct clades with high bootstrap values which were designated as T. equi clades/ genotypes A, B, C and D. Our results indicated that the genotype B of Nagore et al. and genotype E of Qablan et al. together formed the clade B with a high bootstrap value (95%). Furthermore, all the genotypes probably originated from clade B, which was the most dominant genotype (52.85%) followed by clades A (27.58%), and C (9.78%) and D (9.78%). Genotype C manifested a comparatively higher genetic diversity (91.0-100% identity) followed by genotypes A (93.2-99.5%), and B and D (95.7-100%). The alignment report of the consensus nucleotide sequences of the V4 hypervariable region of the 18S rRNA gene of four T. equi genotypes (A-D) revealed significant variations in one region, between nucleotide positions 113-183, and 41 molecular signatures were recognized. As far as geographical distribution is concerned, genotypes A and C exhibited far-extending geographical distribution involving 31 and 13 countries of the Asian, African, European, North American and South American continents, respectively. On the contrary, the genotypes B and D exemplified limited distribution with confinement to 21 and 12 countries of Asian, African and European continents, respectively. Interestingly, genotypes A and C have been reported from only two continents, viz., North and South America. It was observed that genotypes A and C, and B and D exhibit similar geographical distribution.

Discussion

The present study indicated the presence of only four previously described T. equi genotypes (A, B, C and D) after performing the molecular analyses of all available sequences of the complete V4 hypervariable region of the 18S rRNA gene of T. equi isolates in the GenBank™.

Population genetic characterization of Theileria annulata based on the cytochrome b gene, with genetic insights into buparvaquone susceptibility in Haryana (India)

The present investigation was aimed at population genetic characterization of Theileria annulata on the basis of the cytochrome b (cyt b) gene along with the evaluation of status of buparvaquone resistance in Haryana (India). The sequences originating from China, Egypt, India, Iran, Iraq, Tunisia, Turkey and Sudan were included in the analysis. The maximum likelihood tree based on the Tamura-Nei (TN93+G) model placed all the sequences of T. annulata into a single clade. The median-joining haplotype network exemplified geographical clustering between T. annulata haplotypes originating from each country. Only five haplotypes (7.81 %) were shared between any two countries, while the remaining 59 haplotypes (92.19 %) were singleton and unique to one country. The values of pairwise genetic distance (FST) between all the populations indicated huge genetic differentiation (> 0.25) between different T. annulata populations, barring the FST value between Iraq and Turkey (0.14454) which suggested a moderate differentiation. Contrary to the FST index, the values of gene flow (Nm) between T. annulata populations were very low. The neutrality indices and mismatch distributions indicated a population expansion in the Indian T. annulata population. Furthermore, the secondary structure and homology modeling of the partial cyt b protein is also reported. The molecular analysis of newly generated sequences for buparvaquone resistance revealed that all the isolates were susceptible to buparvaquone treatment. However, two novel mutations at positions V203I and V219I in between the Q01 and Q02 drug-binding regions of the cyt b gene were observed for the first time.

Phylogenetic analysis and haplotype networking of Hepatozoon felis infecting wild animals in Gir National Park, Gujarat, India

The present study was aimed to determine the phylogenetic relationship, haplotype network, and demographic dynamics of H. felis infecting the endangered Asiatic lions in Gir National Park, Gujarat, India, on the basis of partial 18S rRNA gene. The phylogenetic analysis based on the partial 18S rRNA gene sequences of H. felis exhibited the presence of two distinct genotypes of H. felis (HfG1 and HfG2) infecting the Indian wild felids, viz., the Asiatic lion, Royal Bengal tiger, and Indian leopard. The HfG1 and HfG2 genotypes exhibited 97.6-100% and 99.7-100%, and 96.9-98.7% nucleotide identity within and between themselves, respectively. The HfG1 genotype exhibited a higher genetic diversity as compared to HfG2. A total of 22 molecular signatures were identified in the 18S rRNA gene between these genotypes. Further, analysis of a total 67 sequences of H. felis (13 different host species from 13 countries of Africa, South America, Asia, and Europe) that were downloaded from GenBankTM, generated 30 haplotypes. Among all the haplotypes, Hap_17 (h=12) was the most frequent followed by Hap_12 (h=09) and Hap_4 (h=05). Out of 13 location-wise populations, India (h=12) shared the highest number of haplotypes followed by Japan (h=08), and the least number of haplotypes were found in Hungary (h=02). Population dynamics study involving neutrality tests and mismatch distribution, and genetic differentiation indices, revealed the presence of phylogeographic population structure and a constant population size indicating a uniform gene flow among the populations worldwide. In conclusion, a high genetic diversity along with the presence of two distinct genotypes of H. felis were observed on the basis of 18S rRNA gene sequence analysis.

An insight into misidentification of the small-subunit ribosomal RNA (18S rRNA) gene sequences of Theileria spp. as Theileria annulata

BACKGROUND

There had been isolated reports of the presence of novel Theileria annulata genotypes based on the 18S rRNA gene sequence data from India, Pakistan and Saudi Arabia; but, these studies were restricted to limited field samples. Additionally, no comparative study has been conducted on all the isolates of this parasite from different countries whose sequences are available in the nucleotide databases. Therefore, we aimed to study the genetic diversity of T. annulata based on all available nearly complete 18S rRNA gene sequences in the GenBank™. Out of a total of 312 gene sequences of T. annulata available in the NCBI database, only 70 nearly complete sequences (> 1527 bp) were used for multiple sequence alignment.

RESULTS

The maximum likelihood tree obtained using TN93 + G + I model manifested two major clades. All the valid host-cell transforming Theileria species clustered in one clade. The T. annulata designated sequences occupying this clade clustered together, excluding two isolates (DQ287944 and EU083799), and represented the true T. annulata sequences (n = 54). DQ287944 and EU083799 exhibited close association with Theileria lestoquardi. In addition, 14 Indian sequences formed a large monophyletic group with published Theileria orientalis sequences. The broad range of sequence identity (95.8-100%) of T. annulata designated sequences indicated the presence of different Theileria spp. A closer analysis revealed the presence of three Theileria spp., namely, T. annulata, T. orientalis, and two isolates (DQ287944 and EU083799) closely related to T. lestoquardi. The true T. annulata sequences manifested 98.8-100% nucleotide identity within them. EU083799 and 14 misidentified Indian T. annulata sequences exhibited the highest similarity with T. lestoquardi (98.6-98.8%) and T. orientalis (98.0-99.9%) in comparison with the other Theileria spp. of domestic and wild ruminants.

CONCLUSION

In the course of analyzing the genetic diversity of T. annulata, we identified the nearly complete 18S rRNA gene sequences of other Theileria spp. that have not only been misidentified as T. annulata in the GenBank™, but are also published as T. annulata. Moreover, a high level of sequence conservation was noticed in the 18S rRNA gene of true T. annulata and T. orientalis sequences.

Genetic and population diversity of Toxocara cati (Schrank, 1788) Brumpt, 1927, on the basis of the internal transcribed spacer (ITS) region

The present investigation was aimed to study the sequence, phylogenetic and haplotype analyses of Toxocara cati based on the ITS region, along with the genetic diversity, demographic history and population-genetic structure. The maximum likelihood tree based on Kimura 2-parameter model was constructed using the complete ITS region of all the nucleotide sequences (n = 57) of Toxocara spp. and other related ascarid worms available in the GenBank™. It placed all the sequences of T. cati into four major clades designated as T. cati genotypes 1-4 (TcG1-G4). A total of 66 signature nucleotides were identified in the ITS region between genotypes. The median-joining haplotype network displayed a total of 24 haplotypes, with China exhibiting the highest number of haplotypes (h = 20) followed by India (h = 4), and Japan and Russia (h = 1). It indicated a clear distinction between all the four genotypes. The pairwise F_ST values between all the genotypes indicated huge genetic differentiation (> 0.25) between different T. cati genotypes. Moreover, the gene flow (Nm) between T. cati genotypes was very low. Results of AMOVA revealed higher genetic variation between genotypes (92.82%) as compared to the variation within genotypes (7.18%). The neutrality indices and mismatch distributions for the G1-G4 genotypes, Indian isolates and the overall dataset of T. cati indicated either a constant population size or a slight population increase. The geographical distribution of all the genotypes of T. cati is also reported. This is the first report of genotyping of T. cati on the basis of the ITS region.

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 F_s, and R_2, 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 F_s 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.

Genetic diversity in the Tams1 gene of Theileria annulata (Duschunkowsky and Luhs, 1904) infecting cattle

The present study describes the genetic diversity in the Tams1 gene (733 bp) of Theileria annulata along with the sequence, phylogenetic and haplotype analyses of the Indian isolates. The phylogenetic analyses displayed distinct clustering of the Indian isolates into three groups suggesting the presence of three genotypes, hitherto designated as T. annulata genotypes 1-3 (G1-G3). Genotype 3 seems to be novel containing only two newly generated sequences. Indian isolates displayed 88.4-100% and 82.2-100% similarity with each other at nucleotide (nt) and amino acid (aa) levels, respectively. However, the newly generated sequences (n = 36) showed 90.5-100% and 84.3-100% identity between them at nt and aa levels, respectively. The most diverse and heterogeneous genotype, G1, exhibited the highest number of polymorphic sites (S = 148), haplotypes (h = 16) and nucleotide differences (k = 43.23) besides haplotype (Hd = 0.903 ± 0.031) and nucleotide (π = 0.059 ± 0.005) diversities. Neutrality indices suggested a respective decrease and increase in population sizes of G1 and G2 genotypes in India. The nucleotide sequence analyses indicated the presence of extensive sequence variations between nucleotide positions 1-124, 194-257 and 396-494. The N-terminus of Tams1 protein displayed a considerable sequence variability with extensive variations in two regions, between amino acid positions 1-39 and 127-172, as compared to the conserved carboxyl terminus.

Phylogenetic analysis, genetic diversity and geographical distribution of Babesia caballi based on 18S rRNA gene

The present investigation was aimed to study the presence of Babesia caballi clades upon phylogenetic analysis of all available V4 hypervariable 18S rRNA gene sequences in GenBank in addition to the intra- and interclade genetic diversity in B. caballi and the distribution of parasite clades in different countries. Out of altogether 155 small-subunit ribosomal RNA gene sequences of B. caballi available in the database, only 92 sequences with a complete V4 hypervariable region (>293 bp) were used in multiple sequence alignment. The phylogenetic tree placed all the sequences into two distinct clades with high bootstrap values which are designated as B. caballi clades A and B. Clade A was further divided into two subclades A1 and A2 with 98% bootstrap support. On the contrary, clade B contained multiple small subclades which either lacked bootstrap support or did not have enough bootstrap support to further group them into subclades. All the sequences of B. caballi were 91.5-100% identical with each other. Clade B manifested a comparatively higher genetic diversity (95.2-100% identity) amongst sequences as opposed to clade A (97.3-100% identity). Moreover, it indicated 91.5-93.5%, 92.9-94.6% and 91.5-94.6% nucleotide identity with B. caballi subclades A1, A2, and clade A, respectively. Significant nucleotide variations were observed in one region, between nucleotide positions 126-178, in some of the sequences. A total of 21 molecular signature residues were identified in the V4 hypervariable region. The alignment report of the V4 hypervariable region of 18S rRNA gene of clades A and B exhibited nucleotide variation at nine and 24 places, respectively. The distribution map of all the clades of B. caballi is also reported. The number of 18S rRNA gene sequences employed in the study is relatively high compared to previous studies. Therefore, a fair comparison of definite genetic variations between isolates/sequences from different countries was carried out.

Climate change and heat stress: Impact on production, reproduction and growth performance of poultry and its mitigation using genetic strategies

Heat stress is an important environmental determinant which adversely affects the performance of poultry worldwide. The present communication reviews the impact of heat stress on production, reproduction and growth performance of poultry, and its alleviation using genetic strategies. The adverse effects of high environmental temperature on poultry include decrease in growth rate, body weight, egg production, egg weight, egg quality, meat quality, semen quality, fertility and hatchability, which cause vast financial losses to the poultry industry. High ambient temperature has an antagonistic effect on performance traits of the poultry. Thus, selection of birds for high performance has increased their susceptibility to heat stress. Additionally, heat burden during transportation of birds from one place to another leads to reduced meat quality, increased mortality and welfare issues. Molecular markers are being explored nowadays to recognize the potential candidate genes related to production, reproduction and growth traits for selecting poultry birds to enhance thermo-tolerance and resistance against diseases. In conclusion, there is a critical need of formulating selection strategies based on genetic markers and exploring more genes in addition to HSP25, 70, 90, H1, RB1CC, BAG3, PDK, ID1, Na, F, dw and K responsible for thermoregulation, to improve the overall performance of poultry along with their ability to tolerate heat stress conditions.

Genetic characterization and phylogenetic analysis of Sarcocystis suihominis infecting domestic pigs (Sus scrofa) in India

A total of 57 tissue samples of domestic pigs (Sus scrofa) were collected from the meat outlets of five north Indian states and examined for sarcocystosis by histological and molecular methods. The genomic DNA extracted from five representative positive isolates was subjected to PCR amplification of the partial 18S rRNA gene followed by cloning and sequencing. Sequence analysis of the newly generated Indian isolates recorded 96.9-100.0% identity with published sequences of Sarcocystis suihominis. Two new haplotypes that have not been previously described manifested 99.5-100.0% nucleotide homology within themselves. In the phylogenetic analysis, Indian isolates of S. suihominis grouped together with S. suihominis originating from Italy, and they collectively formed a sister clade with Sarcocystis miescheriana within a clade containing various Sarcocystis spp. of ruminants having felids as final hosts. At the same time, this clade separated from a sister clade containing Sarcocystis spp. of bovid or cervid ruminants using canids as known or surmised definitive host. The current study established the phylogenetic relationship of Indian isolates of S. suihominis with various Sarcocystis spp. as well as with other taxa of Sarcocystidae family based on 18S rRNA gene for the first time.

Immune response against subclinical haemonchosis in Himalayan hill goats

Haemonchosis commonly occurs as chronic and subclinical infection in small ruminants, and understanding of immunological response against subclinical haemonchosis is of paramount importance for designing and implementing effective control strategies. The present study was designed to evaluate immunological response during subclinical haemonchosis, experimentally established in goats. Sixteen 5-6 month-old helminth naive kids were randomly allocated into one of two groups, infected and uninfected; the infected group being infected per os with 250 Haemonchus contortus larvae per kg body weight. Faecal, blood and serum samples were collected every third day up to 30 days post-infection (DPI), thereafter weekly up to 58 DPI to record changes in faecal egg count (FEC), haemoglobin (Hb), packed cell volume (PCV), peripheral eosinophil percentage and immunological parameters, such as macrophage cytokine interleukin-12 (IL-12), Th1 cytokine (IFN-γ), Th2 cytokines (IL-4, 13, 25, 33) and immunoglobulins (IgG and IgE). Pre-patent period of H. contortus in the present study was 18 days and eggs per gram (EPG) peaked on 30 DPI. The total reduction in body weight gain in the infected group was 26 g per day when compared with uninfected animals. Hb (7.35 ± 0.34 g/dL in infected animals compared with 9.76 ± 0.67 in control animals) and PCV levels (22 ± 1.54 g/dL in infected animals compared with 29.2 ± 1.27 in control animals) decreased significantly up to 44 DPI in infected group (P = 0.000). IL-4, IL-13, IL-33, IgG and IgE showed significant increase in infected animals at different periods. IFN-γ, IL-12 and IL-25 did not show any significant changes barring a steep rise of IFN-γ on 27 DPI. A positive correlation was observed between IgE and IL-4 in subclinical haemonchosis. Of particular note was that all the major cytokines, such as IFN-γ (P = 0.000), IL-4 (P = 0.000), IL-13 (P = 0.009), and both IgG (P = 0.000) and IgE (P = 0.003), were observed at the lowest concentration on 24 DPI. The effect of infection was found to be significant on cytokines with a strong interaction with time. Taken together, the data suggest that Th2 immune response is predominating in subclinical haemonchosis. The economic loss in term of body weight gain due to subclinical haemonchosis was considerable.