Bioinformatic comparison of Kunitz protease inhibitors in Echinococcus granulosus sensu stricto and E. multilocularis and the genes expressed in different developmental stages of E. granulosus s.s.

Background Cystic and alveolar echinococcosis caused by the tapeworms Echinococcus granulosus sensu stricto ( s.s. ) and E. multilocularis , respectively, are important zoonotic diseases. Protease inhibitors are crucial for the survival of both Echinococcus spp . Kunitz-type inhibitors play a regulatory role in the control of protease activity. In this study,we identified Kunitz-type domain protease inhibitors(KDPIs) present in the genomes of these two tapeworms and analyzed the gene sequences using computational, structural bioinformatics and phylogenetic approaches to evaluate the evolutionary relationships of these genes. Hi-seq transcriptome analysis showed that E. granulosus s.s. KDPIs were differentially expressed in the different developmental stages. We validated some of the genes expressed in adult worm, protoscolex and cyst germinal membrane of E. granulosus s.s. and E. multilocularis by quantitative PCR. Results A total of 19 genes from E. multilocularis and 23 genes from E. granulosus s.s. were predicted to be KDPIs with the most containing a single Kunitz-domain. A maximum likelihood method phylogenetic tree indicated that the E. granulosus s.s. and E. multilocularis Kunitz domain peptides were divided into three branches containing 9 clusters. The ratio of positively charged residues and neutral residues are different between E. multilocularis and E. granulosus s.s. KDPIs. We also found that E. multilocularis had higher percentage of sequences containing signal peptides (17/19, 89.47%) than that of E. granulosus s.s. (14/23, 60.87%). Transcript analysis showed all the E. granulosus s.s. KDPI genes were expressed differentially in four developmental stages of the worm. Transcription analysis showed that 9 KDPIs (including EG_07244,EGR_08716 and EGR_10096) were highly upregulated in adult worm, and 2 KDPIs (EG_09268 and EG_09490) were highly expressed in the cyst germinal membrane. Quantitative gene expression analysis(qPCR) of four genes confirmed the expression of these genes. EGR_08716 and its homologous gene (EmuJ_001137000) were highly and specifically expressed in adult worms of the two worms. Conclusions A total 19 and 23 KDPIs were identified in the genomes of E. multilocularis and E. granulosus s.s. , respectively. The differential expression of these KDPIs in different stages may indicate their different roles in the different hosts. The difference in characterization of KDPIs may be associated with the different pathology of metacestode stage of these two parasites.