Bioinformatics analysis and construction of phylogenetic tree of aquaporins from Echinococcus granulosus

Prediction and Identification of Epitopes in the Emy162 Antigen of Echinococcus multilocularis

INTRODUCTION

Alveolar echinococcosis (AE) is a zoonotic disease caused by the parasitism of Echinococcus multilocularis larvae in the intermediate host or the final host. This study aims to identify and analyze the B-cell and T-cell (Th1, Th2 and Th17) epitopes of E. multilocularis antigen Emy162.

METHODS

(1) The secondary structural characteristics of the Emy162 protein were predicted by bioinformatics software to further predict the potential T- and B-cell epitopes. (2) The dominant antigen epitopes were detected by ELISA through the reaction of patient serum with small B-cell antigen peptide and assessing the proliferation of splenic lymphocytes of mice immunized with Emy162. (3) The expression of cytokines in splenic lymphocytes of mice stimulated by small T-cell antigen peptides was detected by ELISA, ELISpot and flow cytometry to enable the identification of the T-cell epitopes.

RESULTS

(1) The high-scored T-cell epitopes were located at positions E7-13, E36-41, E80-89, E87-96, E97-106 and E129-139, while B-cell epitopes were located at positions E7-13, E19-27, E28-36, E37-48, E78-83, E101-109, E112-121 and E129-139. (2) The three advanced antigen epitopes of Emy162 were E19-27, E112-121 and E129-139. (3) The four Th1 advanced antigen epitopes of Emy162 were E7-13, E36-41, E80-89 and E129-139. The three Th2 advanced antigen epitopes were E36-41, E87-96 and E97-106. The three Th17 advanced antigen epitopes were E36-41, E87-96 and E97-106.

CONCLUSION

(1) The Emy162 protein has advanced antigenicity and numerous potential epitopes. Six T-cell and eight B-cell dominant epitopes were revealed using bioinformatics methods. (2) There are three dominant B-cell epitopes, four dominant Th1 epitopes, three dominant Th2 epitopes, and three dominant Th17 epitopes in the Emy162 antigen.

Comparison between Echinococcus granulosus sensu stricto (G1) and E. canadensis (G6) mitochondrial genes (cox1 and nad1) and their related protein models using experimental and bioinformatics analysis

BACKGROUND

Cystic echinococcosis (CE) as a zoonotic parasitic disease, remains a health challenge in many parts of the world. There are different species of Echinococcus granulosus sensu lato with different pathogenicity and host preferences.Different procedures have been applied for characterization of Echinococcus taxa in which two mitochondrial genes, cox1 and nad1 have been used more common. They have been able to differentiate E. granulosus sensu stricto and E. canadensis species in different hosts. The affinity of E. granulosus sensu stricto and E. canadensis species for localizing different organs seems to be different. To what such affinity and related pathogenicity could be related, is not known, so far. Bioinformatics analysis may be helpful to interpret such difference by investigating the genes and their related protein models between different species infecting human and animals. The current work was designed to study the differences between E. granulosus s.s. and E. canadensis species mitochondrial genes (cox1 and nad1) and related protein models of CE cysts by experimental and bioinformatics analysis.

MATERIALS AND METHODS

Different human and animal CE cysts were collected and their DNA was extracted and sequenced based on their cox1 and nad1 genes. In order to determine the E. granulosus s.s. and E. canadensis species of the samples, BLAST analysis was performed on sequenced genes. Three sequences were selected for analysis and were deposited in GenBank. Moreover, the sequence number of KT988116.1 which belonged to E. canadensis from our already deposited in GenBank was also selected. Alignment and phylogenetic analysis were performed on the sequences using BioEdit and MEGA7 software. The raw sequences of translated proteins belonged to the mentioned genes were obtained from Protein database in NCBI. The secondary structure was determined by PSIPRED Protein Sequence Analysis Workbench. The tertiary models of COX1 and NAD1 proteins in both genotypes were constructed using Modeler 9.12 software and their physicochemical features were computed using ProtParam tool in ExPASY server.

RESULTS

BLAST analysis on sequenced genes showed that the samples belonged to E. granulosus s.s. and E. canadensis species. These sequences were deposited in GenBank with accession numbers: JN579173.1, KF437811.1, and KY924632.1. The results showed that proteins of COX1 of E. granulosus s.s., COX1of E. canadensis, NAD1of E. granulosus s.s. and NAD1of E. canadensis species, consisted of 135, 122, 120 and 124 amino acids, respectively. The aligned sequences of translated proteins belonged to COX1 and NAD1 enzymes in E. granulosus s.s. and E. canadensis species were different; such that alignment COX1 sequence between E. granulosus s.s. and E. canadensis species showed that amino acids were different in 6 positions. This difference for NAD1 sequences were different in 19 positions. The secondary structure determined by PSIPRED showed differences in coil, strand and helix chains in COX1 and NAD1 proteins in E. granulosus s.s. and E. canadensis species. Comparison between three-dimensional structures (3D) of COX1 protein model in E. granulosus s.s. and E. canadensis species demonstrated an additional helix with two conserved iron binding sites in the COX1 protein of E. granulosus s.s. species.

CONCLUSION

E. granulosus s.s. and E. canadensis species differences are reflected in two important proteins: COX1 and NAD1. These differences are demonstrable in the 3D structure of proteins of both strains. So, the present study is adding to our understanding of the difference in molecular sequences between the E. granulosus s.s. (G1) and E. canadensis (G6) which may be used for interpreting the difference between the pathogenicity and localization affinity in these two important helminthic zoonosis.