Theileria annulata is a tick-associated parasite that causes tropical theileriosis in livestock and is responsible for huge economic losses. Studies have been neglected on the effect of Theileria spp. on cattle in Khyber Pakhtunkhwa (KP), Pakistan. The present study was designed to determine the genetic diversity and assess the risk factors associated with Theileria infection in selected districts of KP. Information on the risk factors related to the Theileria infection was collected through a questionnaire. Blood samples were collected from symptomatic cattle from January 2019 to February 2020, identified morphologically through microscopic examination, and processed for molecular characterization using the 18S rRNA gene as a genetic marker. Of the 555 cattle examined (136, 24.5%) and (294, 53%) were found positive for Theileria spp. by microscopic examination and a PCR test, respectively. Based on the PCR test, the highest prevalence of infection was found in district Upper Dir (46/75, 61.3%), followed by Lower Dir (54/90, 60%); Malakand (51/88, 57.9%); Peshawar (40/80, 50%); and Charsadda (52/112, 46.4%), with the lowest in Bajaur (51/110, 46.34%). A BLAST analysis of the 18S rDNA sequences showed 99.5% identity with T. annulata. In a phylogenetic tree, the 18S rDNA sequence of T. annulata clustered with sequences from Pakistan, China, and Italy. A significant association was observed between the prevalence of infection and different host characteristics. The highest infection was found in adult cattle (216/360, 60%); females (218/377, 57.8%); and Holstein Friesian (120/180, 66.6%). Theileria infection was significantly associated with management practices. Higher infection rates were observed in free-grazing cattle (190/412, 42.2%); those kept in unhygienic conditions (246/405, 60.7%); cattle in combined farming systems (165/255, 64.8%); and those in congested stall systems (150/218, 68.8%). Seasonal patterns were found to be significantly associated with infection, and a higher infection rate was observed in summer (215/350, 61.4%) than in winter (79/205, 38.5%). Identified risk factors should be considered in designing practical control approaches to reduce the burden of Theileria infection. Large scale studies are required to explore the diversity of Theileria species in KP, Pakistan.

Despite high diversity in the Oriental region, ticks of the genus Haemaphysalis have been neglected regarding their genetic data and vector potential. This study aimed to genetically characterize three species of the genus Haemaphysalis: Haemaphysalis cornupunctata, Haemaphysalis kashmirensis, and Haemaphysalis montgomeryi infesting goats and sheep, and Rickettsia spp. associated with these tick species in the Hindu Kush Himalayan range of Pakistan. Altogether, 834 ticks were collected by examining 120 hosts including goats (64/120, 53.3%) and sheep (56/120, 46.6%), in which 86 (71.6%) hosts were found to be tick-infested. The morphologically identified ticks were subjected to DNA extraction and PCR for the amplification of partial 16S rDNA and cox fragments. Rickettsia spp. associated with the collected ticks were detected through the amplification of gltA, ompA and ompB partial fragments. The 16S rDNA of H. cornupunctata and H. montgomeryi showed a maximum identity of 100% with the sequences of the same species, whereas the 16S rDNA of H. kashmirensis showed the highest identity of 93–95% with Haemaphysalis sulcata. The cox sequence of H. montgomeryi displayed 100% identity with the same species. In comparison, the cox sequences of H. cornupunctata and H. kashmirensis showed maximum identities of 87.65–89.22% with Haemaphysalis punctata and 89.34% with H. sulcata, respectively. The gltA sequence of Rickettsia sp. from H. kashmirensis showed the highest identity of 97.89% with Rickettsia conorii subsp. raoultii, while the ompA and ompB fragments from the same DNA samples revealed 100% and 98.16% identity with Rickettsia sp. and “Candidatus Rickettsia longicornii”, respectively. Another gltA sequence amplified from H. montgomeryi ticks showed 100% identity with Rickettsia hoogstraalii, while the attempts to amplify ompA and ompB for R. hoogstraalii were unsuccessful. In the phylogenetic tree, the 16S rDNA of H. cornupunctata clustered with the corresponding species while its cox clustered with H. punctata. Both 16S rDNA and cox sequences of H. kashmirensis clustered with H. sulcata. The gltA sequence of Rickettsia sp. was clustered individually in the spotted fever (SF) group of Rickettsia, while the gltA sequence of R. hoogstraalii was clustered with the same species in the transition group of Rickettsia. In the SF group, the rickettsial ompA and ompB sequence clustered with undetermined Rickettsia sp. and “Candidatus Rickettsia longicornii”, respectively. This is the earliest study regarding the genetic characterization of H. kashmirensis. This study indicated that ticks belong to the genus Haemaphysalis have the potential of harboring and/or transmitting Rickettsia spp. in the region.

Single-nucleotide polymorphisms (SNPs) are reported to be associated with many diseases, including autoimmune diseases. In rheumatoid arthritis (RA), about 152 SNPs are reported to account for ~15% of its heritability. These SNPs may result in the alteration of gene expression and may also affect the stability of mRNA, resulting in diseased protein. Therefore, in order to predict the underlying mechanism of these SNPs and identify novel therapeutic sites for the treatment of RA, several bioinformatics tools were used. The damaging effect of 23 non-synonymous SNPs on proteins using different tools suggested four SNPs, including rs2476601 in PTPN22, rs5029941 and rs2230926 in TNFAIP3, and rs34536443 in TYK2, to be the most damaging. In total, 42 of 76 RA-associated intronic SNPs were predicted to create or abolish potential splice sites. Moreover, the analysis of 11 RA-associated UTR SNPs indicated that only one SNP, rs1128334, located in 3′UTR of ETS1, caused functional pattern changes in BRD-BOX. For the identification of novel therapeutics sites to treat RA, extensive gene–gene interaction network interactive pathways were established, with the identification of 13 potential target sites for the development of RA drugs, including three novel target genes. The anticipated effect of these findings on RA pathogenesis may be further validated in both in vivo and in vitro studies.

As a vector of wide range of pathogenic agents, ticks pose health threats to wild and domestic animals, and humans. Information is unavailable about the prevalence and spatial survey of Hyalomma kumari ticks and associated Rickettsia spp. in Pakistan. Concerning this knowledge gap, the present study aimed to molecularly detect Rickettsia species associated with H. kumari infesting small ruminants in Khyber Pakhtunkhwa (KP), Pakistan. A total of 409 H. kumari ticks were collected from 163/295 infested hosts with an infestation rate of 55.25%. A total of 204 females, 158 males, and 47 nymphs were collected. Goats were heavily infested by 224 ticks having an infestation rate of 58.33% (98/168), whereas sheep were infested by 185 ticks having a lesser infestation rate of 51.18% (65/127). Genomic DNA extracted from ticks was used for the amplification of tick (cox I, 16S rRNA, ITS-2) species and Rickettsia (gltA, ompA, and ompB) partial genes. Eighty-three ticks were subjected to PCR, and 8/83 (9.6%) were found positive for rickettsial agents. The cox I and 16S rRNA sequences of H. kumari showed 98.90–99.74% identity with H. kumari sequences reported from Pakistan, and phylogenetically clustered to the corresponding species reported from Pakistan and India. The obtained rickettsial gltA, ompA, and ompB sequences showed 100% identity with Rickettsia sp. of the Rickettsia conorii reported from Pakistan. In the phylogenetic trees, rickettsial sequences clustered with uncharacterized Rickettsia sp. from Pakistan and R. conorii from Israel, Russia, South Africa, and India. The present molecular based detection of H. kumari-associated R. conorii will facilitate effective surveillance in the region.

Bartonella henselae is a Gram-negative bacterium causing a variety of clinical symptoms, ranging from cat-scratch disease to severe systemic infections, and it is primarily transmitted by infected fleas. Its status as an emerging zoonotic pathogen and its capacity to persist within host erythrocytes and endothelial cells emphasize its clinical significance. Despite progress in understanding its pathogenesis, limited knowledge exists about the virulence factors and regulatory mechanisms specific to the B. henselae strain Houston-1. Exploring these aspects is crucial for targeted therapeutic strategies against this versatile pathogen. Using reverse-vaccinology-based subtractive proteomics, this research aimed to identify the most antigenic proteins for formulating a multi-epitope vaccine against the B. henselae strain Houston-1. One crucial virulent and antigenic protein, the PAS domain-containing sensor histidine kinase protein, was identified. Subsequently, the identification of B-cell and T-cell epitopes for the specified protein was carried out and the evaluated epitopes were checked for their antigenicity, allergenicity, solubility, MHC binding capability, and toxicity. The filtered epitopes were merged using linkers and an adjuvant to create a multi-epitope vaccine construct. The structure was then refined, with 92.3% of amino acids falling within the allowed regions. Docking of the human receptor (TLR4) with the vaccine construct was performed and demonstrated a binding energy of −1047.2 Kcal/mol with more interactions. Molecular dynamic simulations confirmed the stability of this docked complex, emphasizing the conformation and interactions between the molecules. Further experimental validation is necessary to evaluate its effectiveness against B. henselae.

Rheumatoid arthritis (RA) is one of the complex diseases with the involvement of the genetic as well as environmental factors in its onset and severity. Different genome-wide association and candidate gene studies have shown the role of several genetic variants in multiple loci/genes with ethnical and geographical variations. This study was designed to detect the association of a single-nucleotide polymorphism (SNP) rs10865035 in the AFF3 gene with the genetic background of rheumatoid arthritis (RA) in the Pakistani cohort. A total of 703 individuals, including 409 RA patients and 294 healthy controls, were genotyped using TaqMan assay and Tri primer ARMS-PCR (amplification-refractory mutation system-polymerase chain reaction) methods. The association of rs10865035 with the RA was statistically determined using different models. Interestingly, besides the homozygous recessive model (G/G vs. A/G + A/A) (OR = 1.693(1.06–2.648); = 0.025), all other models, which included the codominant (χ² = 5.169; = 0.075), homozygous dominant (A/A vs. G/G + A/G) (OR = 0.867 (0.636–1.187); = 0.41), heterozygous (A/G vs. A/A + GG) (OR = 0.491 (0.667–1.215); = 0.49), and additive model (OR = 0.826 (0.665–1.027); = 0.08) showed insignificant distribution of the genotypes among the cases and controls. These findings suggest that the AFF3 gene (rs10865035) has no significant role in the onset of RA in the Pakistani population.