Genetic variation of mitochondrial genes among Echinococcus multilocularis isolates collected in western China
Parasites & Vectors volume 10, Article number: 265 (2017)
Alveolar echinococcosis (AE) is a life-threatening human disease caused by Echinococcus multilocularis transmitted between rodents and dogs/foxes in the Northern Hemisphere. The study aims to identify the genetic variation of the parasite in AE patients from China.
E. multilocularis isolates were collected from wild small mammals (n = 6) and AE patients (n = 56) from western China. Genomic DNA was extracted from different tissue samples including paraffin tissue blocks, ethanol fixed tissues and frozen tissues surgically removed. Two mitochondrial gene fragments (526 bp for cob and 474 bp for nad2) of E. multilocularis were amplified and sequenced.
The parasite fragment sequences of cob fragments from AE patients showed two haplotypes, and nad2 gene fragment sequences had four haplotypes. The gene sequences from Microtus sp. were 100% identical to the sequences of some isolates from AE patients. These haplotypes were distributed in both Qinghai and Xinjiang provinces. Alignment analysis with the sequences from the GenBank databases showed five genotypes including three Asian genotypes, one from Europe and one from North America.
Most AE patients harbored the Asian genotype 1 which may be an indication of its relative frequency in the definitive hosts and the environment or of its pathogenicity to humans, which calls for further research.
Echinococcus multilocularis is a fox/dog tapeworm which causes human alveolar echinococcosis (AE). This zoonosis is one of the most lethal infections in humans worldwide [1, 2]. The tapeworm needs two hosts to complete its life-cycle: a definitive host such as a fox or a dog, and a small mammal intermediate host (vole, pika, zokor, hamster, jird etc.) . The transmission of the disease is sustained in wildlife by the predator–prey relationship between the fox and small mammals. Domestic dogs play also an important role in the transmission of the disease to humans in central Asia including in western China [4,5,6,7,8]. Foxes or dogs are infected by feeding on small mammals infected with fertile E. multilocularis vesicles containing larval protoscoleces (PSCs). Once ingested, PSCs are released from their vesicles and develop into mature adult tapeworms producing eggs in the intestine of the carnivorous animal 28–35 days post-infection. The eggs are released through passing of dog/fox feces. Humans are infected as an intermediate host by accidently swallowing parasite eggs. However, humans are an epidemiological deadlock and are not involved in the life-cycle.
AE is more severe in humans than cystic echinococcosis which is caused by Echinococcus granulosus. The high mortality is due to the larval stage causing severe damages to the liver. About 95% of AE patients die within ten years in the absence of treatment . Recently, many human AE cases were reported on the Qinghai-Tibet plateau [3, 6, 7, 9, 10]. It is estimated that western China including Xinjiang, Sichuan, Qinghai, Ningxia and Inner Mongolia, has 91% of total AE cases in the world .
Earlier studies have clearly shown that E. multilocularis is widely distributed in the Northern Hemisphere, and the large diversity of life-cycle patterns and host communities throughout its distribution range are currently described. Large genetic variations within the parasite species might be expected from such range and diversity as a response to various hosts and environmental adaptation. Paradoxically, little genetic variation was detected within E. multilocularis which was initially classified into two genotypes only, M1 (Europe) and M2 (China, Alaska and North America) based on four nucleotide substitutions out of 471 nt of mitochondrial NADH dehydrogenase subunit 1 gene . However, based on sequence data of mitochondrial and nuclear DNA, further geographical distinctions were found with European, Asian and North American clades, and one unrelated haplotype from Inner Mongolia . In fact, there are very few mitochondrial sequences of E. multilocularis from AE patients deposited in the GenBank database. In the present study, we sequenced two mitochondrial (mt) genes (fragments) from parasitic liver tissues of 56 human AE cases and six isolates of E. multilocularis collected from rodent intermediate hosts.
Echinococcus multilocularis samples
Three isolates of E. multilocularis were maintained in Mongolian jirds in our laboratory by transplantation of protoscoleces intraperitoneally for more than 10 years. Among them, Em-NX and Em-XJ were collected from Microtus spp. from southern Ningxia and Xinjiang, respectively. Em-A, an Alaskan strain was obtained from Professor Philip Craig, University of Salford, UK. Three E. multilocularis specimens were isolated recently from the livers of voles (Microtus sp.) from Yili, Xinjiang, China. The livers containing parasite tissues were fixed in 75% ethanol until use.
To identify the variation of E. multilocularis in humans, 56 AE tissue samples were collected from AE patients including 20 paraffin blocks for pathological tests, 12 liver tissues fixed in 70% ethanol and 24 liver tissues frozen at -80 °C. The patients were mainly from Xinjiang, Qinghai, Sichuan and Gansu provinces.
Genomic DNA was extracted from fresh protoscoleces of the three laboratory maintained isolates of E. multilocularis. The protoscoleces were aspirated from the parasite mass from Mongolian jirds. After 10 washes with PBS, about 100 protoscoleces were soaked in the genomic DNA extraction buffer-PrepMan TM Ultra Sample Preparation Reagent (Appied Biosystems, Foster City, USA) and then homogenized using a microtube mortar and pestle (Sigma, St. Louis, USA). The homogenate was heated at 95 °C for 10 min. The solution was directly used as a DNA template for PCR amplification after centrifugation at 12,000× g for 3 min.
Three sources of liver tissues from AE patients were used for molecular analyses including paraffin tissue blocks, ethanol fixed tissues and frozen tissues. To extract DNA from the blocks, the QIAamp DNA FFPE Tissue Kit (Qiagen, Hilden, Germany) was used. Briefly, the paraffin blocks containing liver leisure tissues of AE patients were cut into 8 μm thick sections. About 8–10 sectioned films were placed into a microtube containing 0.5 ml xylene. The extraction procedure was performed as recommended in the manufacturer protocols. The DNA sample was eluted from the column with 50 μl of water, and then stored at -80 °C for further use.
To extract DNA from ethanol fixed tissues (both human and rodent) or frozen tissues, about 0.1 g of tissue was cut into small pieces and then transferred into a microtube containing 387 μl lysis buffer (50 mmol/l Tris, 0.1 mol/l EDTA, 0.1 mol/l NaCl, pH 8.0), 10 μl 20% SDS, and 3 μl proteinase K (20 mg/ml). The tube was incubated at 56 °C overnight. After centrifugation at 13,000× rpm for 10 min, the supernatant was transferred into a new microtube containing 400 μl ethanol. The DNA was pelleted by centrifugation at maximum speed for 15 min then the pellet was washed with 600 μl 75% ethanol and centrifuged at maximum speed for 5 min. The pellet was resuspended with 50 μl of water and stored at -20 °C until used.
PCR amplification of mitochondrial gene fragments
PCR was carried out in a 50 μl reaction mixture including 1 μl template DNA solution, 8 μl of dNTPs, 1 μl of each primer at 10 μM for cob and nad2 gene fragments (Table 1), 1 μl of Ex-Taq Polymerase (TAKARA BIO INC, Shiga, Japan), 5 μl of reaction buffer and added water to 50 μl. PCR reactions were performed for 35 cycles of denaturation (95 °C for 30 s), annealing (57 °C for 30 s) and extension (72 °C for 60 s). Two pairs of primers used were designed for the present study, in order to amplify the partial sequences of the mitochondrial genes cytochrome b (cob), and NADH dehydrogenase subunit 2 (nad2), respectively (Table 1). The size of the two gene fragments was 526 bp and 474 bp, respectively.
DNA sequence alignment and analysis
To sequence the gene fragments, we isolated PCR products from agarose gel (0.8%, w/v) by using the QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany), and then inserted the DNA into an AT vector (pMD19-T Vector, TAKARA BIO INC, Shiga,. Japan)), according to these manufacture’s instructions. Clones containing inserts were sent to sequencing company for sequencing. The sequences obtained were aligned using Bioedit version 5 (http://www.mbio.ncsu.edu/bioedit/bioedit.html/) and compared to the reference sequences from the GenBank database.
Phylogenetic analysis was performed using MEGA 6  and phylogenetic trees were constructed using neighbor-joining (NJ) method [Kimura's two-parameter distance analysis with a gamma shape parameter (a = 0.5)]. The robustness of phylogenetic trees was tested by bootstrapping with 1,000 replicates. A network of mtDNA haplotypes was illustrated by NETWORK 220.127.116.11 using statistical parsimony (http://www.fluxus-engineering.com/ sharenet.htm). We calculated the genetic distance among the subpopulations using MegAlign (www.dnastar.com/t-megalign.aspx) pairwise fixation index (Fst). Fst values close to 1 indicate extreme genetic differentiation between two subpopulations.
Sequence variation among the isolates maintained in laboratory
The size of the amplified DNA fragments obtained from E. multilocularis isolates was 526 bp and 474 bp for cob and nad2, respectively (GenBank accession numbers: KT965443–KT965494; KY290762-KY290791). Figure 1 shows the substitution nucleotides of the two gene fragments amplified from the three isolates maintained in our laboratory and the three isolates we recently collected from wild rodents in Yili, Xinjiang. The comparison with the different haplotypes of the sequences deposited in the GenBank database showed that E. multilocularis is variable at the genomic level. It showed that nad2 sequences isolated from China were variable with a single base substitution compared with the sequences from France, Japan and Slovakia, two or more nucleotide substitutions from those from Canada and USA (Alaska). One of the cob haplotypes from China was identical to the sequences isolated from France and Slovakia, whereas, there are 5 substitutions compared with the sequences from Alaska, USA. The haplotype Em-A is identical to the sequence of the isolate collected from the St. Lawrence Island in Alaska (GenBank: AB461409), which is different from the two laboratory maintained isolates collected from China. The genetic divergence of cob and nad2 gene fragments between Em-NX/Em-XJ and Em-A ranged from 0.85 to 0.95%.
Sequence analysis of E. multilocularis isolated from patients infected with alveolar echinococcosis (AE)
Liver tissue samples were collected from 56 AE patients by surgical removal. In total, 68 sequences were obtained from these samples including 32 nad2 gene sequences (including 24 from Xinjiang, 4 from Qinghai, 3 from Sichuan and 1 from Gansu); 36 cob gene sequences (23 from Xinjiang, 10 from Qinghai, one from Sichuan, Gansu and Beijing, respectively).
There are few E. multilocularis sequences from AE patients deposited in the DNA sequence databases, which may indicate the difficulty to obtain material or to amplify DNA fragments from surgical samples from AE patients. We extracted parasite DNA from three different sources of tissues. From the frozen tissues, 21 cob and 21 nad2 sequences were obtained from the 24 samples available (success rate: 87.5%). From paraffin blocks, we were successful in obtaining nad2 and cob sequences from 60% and 45% of the samples, respectively. From 12 ethanol-fixed samples, we obtained nad2 sequences only from 3 samples and cob sequences from 2 samples (Table 2).
The distance-based NJ analysis of nad2 gene sequences showed three clades referred to as NA1, NA2 and NA3 (Fig. 2a andb). The genetic divergence between NA1 and NA2 was 0.2%, that between NA1 and NA3 was 0.3%, The gap between NA1 and NA2 was not obvious with only one single nucleotide polymorphism (SNP) at position 55 (A/G). Compared to NA2 sequences, AE27 XJWQ had two substitutions at A55 to G55 and G155 to A155 (Fig. 2a) with genetic divergence of 0.4%.
We found some substitutions resulting in amino acid changes, nad2 155G in samples AE1, AE7, AE21, AE22 and AE27 is replaced by 155A in other sequences (Fig. 2a), which changes amino acid glycine (G) to serine (S). Whether the change is associated with the disease severity and pathological difference, and whether those differences are associated with different life-cycle patterns needs to be further elucidated.
The phylogenetic analysis based on the nad2 gene fragments showed that all the sequences isolated from AE patients were close to those isolates collected from rodents. In fact, some, such as sequences of AE2, AE12, AE24 and AE30, were identical to Em-XJ and Em-NX nad2 gene sequences. Alignment of these sequences combined with the sequences (n = 20) in the GenBank database shows that the amplified nad2 gene fragments (474 bp) contained 20 (4.2%) SNPs. Phylogenic analysis in a previous study showed that these sequences are divided into four clades, termed as Asian strain, European strain, North American strain and Alaska strain . All of the sequences from AE patients belonged to the Asian strain (Fig. 2b).
The divergence value of the nad2 fragments of E. multilocularis collected from the three continents ranged from 0.4 to 1.3%. In addition, parsimony network analysis also supports the clade classification of E. multilocularis divided by their geographic locations (Fig. 2c).
The alignment of the 36 cob fragments showed two haplotypes with one nucleotide substitution at position 28(T/C) (Fig. 3a), which matches with the distance-based NJ phylogram analysis showing two genotypes, named CA1 and CA2 (Fig. 3b). The divergence between CA1 and CA2 haplotypes was 0.2%. However, the substitution causes no amino acid sequence change. The two haplotypes coexist in the two major Chinese AE endemic areas, Qinghai and Xinjiang.
The statistical parsimony networks for the nad2 and cob gene fragments illustrated in Figs. 2c and 3c showed that only the Asian genotype/clade is located in the center, the others are scattered. Fig. 2c shows that the nad2 sequences of the European genotype are among the American and Asian genotypes. The European clade contains four substitutions and the Asian clade has three substitutions. The strains from Alaska and Inner Mongolia have five substituted nucleotides compared to the Asian genotype sequences.
The distribution of the genotypes in the world
Including our sequences, a total of 125 sequences from the GenBank database were employed to map the distribution of these genotypes of E. multilocularis worldwide (Fig. 4). The parasite exhibits different genotype patterns in different continents. In Europe, nad2 and cob are conserved, showing only one clade. In North America, there are three clades, the North American genotype (64–72%), the European genotype (21–29%) and the Alaskan genotype (7–8%). In Asia, cob gene sequences allowed us to distinguish two major genotypes.
The variations in the mitochondrial genes cox1 and nad1 have been used to determine the genotypes of E. granulosus (sensu lato) [12, 15, 16], which was classified into nine strains, now five species mainly according to their intermediate hosts, indicating that the host selection may play an important role in driving the genetic evolution. However, there is no evidence so far that such a selection is associated with the genomic variation of Echinococcus. In the study, we obtained 68 mitochondrial sequences for E. multilocularis from AE patients, including 36 cob and 32 nad2 gene sequences, offering basic information for identification of the relationship of host specificity and genetic variation. To our knowledge, this is the largest sequence dataset obtained from humans AE cases. We noted that cox1 of E. multilocularis has been used for identification of the variation within the genus Echinococcus . However, in our study, only one cox1 sequence was generated from the human AE cases (data not shown).
Alignment and phylogenetic analysis of nad2 sequences from the patients showed four haplotypes in three clades, whereas the cob gene sequence analysis showed two haplotypes. It will be interesting to know whether the haplotypes are associated with different pathological patterns of AE which coexisted in two major Chinese AE endemic areas, Qinghai-Tibet Plateau and Xinjiang. However, in both areas, one genotype (Asian genotype 1) was predominant, which may indicate that this genotype of E. multilocularis is suitable to the rodent host species, environment and landscape, which impact the transmission of AE .
The phylogenetic tree and genetic divergence analysis showed that the European and Asian isolates are very close, indicating they had close evolutionary relationships, whereas, the North American and Alaskan genotypes are relatively far from the Asian clade (Figs. 2 and 3).
Sequence analysis of E. multilocularis cob and nad2 gene fragments showed relatively low genetic variation among the isolates from AE patients in western China. Most of the isolates are Asian genotype 1 which may be an indication of its relative high frequency as a source of infection to humans compared to other genotypes, either because it is much more frequent in definitive hosts close to humans (e.g. dogs) or because of a stronger pathogenicity to humans. These hypotheses need to be further elucidated.
We collected E. multilocularis isolates from patients infected with alveolar echinococcosis from western China. The analysis of mitochondrial cob and nad2 gene fragments showed that the parasite exhibits variation at the genomic level with at least three genotypes in China.
- cob :
cytochrome b gene
NADH dehydrogenase subunit 2 gene
Cadavid Restrepo AM, Yang YR, McManus DP, Gray DJ, Giraudoux P, Barnes TS, et al. The landscape epidemiology of echinococcoses. Infect Dis Poverty. 2016;5:13.
McManus DP, Gray DJ, Zhang W, Yang Y. Diagnosis, treatment, and management of echinococcosis. BMJ. 2012;344:e3866:e3866.
Giraudoux P, Raoul F, Afonso E, Ziadinov I, Yang Y, Li L, et al. Transmission ecosystems of Echinococcus multilocularis in China and Central Asia. Parasitology. 2013;140(13):1655–66.
Mastin A, van Kesteren F, Torgerson PR, Ziadinov I, Mytynova B, Rogan MT, et al. Risk factors for Echinococcus coproantigen positivity in dogs from the Alay valley, Kyrgyzstan. J Helminthol. 2015;89(6):655–63.
Afonso E, Knapp J, Tete N, Umhang G, Rieffel D, van Kesteren F, et al. Echinococcus multilocularis in Kyrgyzstan: similarity in the Asian EmsB genotypic profiles from village populations of Eastern mole voles (Ellobius tancrei) and dogs in the Alay valley. J Helminthol. 2015;89(6):664–70.
Moss JE, Chen X, Li T, Qiu J, Wang Q, Giraudoux P, et al. Reinfection studies of canine echinococcosis and role of dogs in transmission of Echinococcus multilocularis in Tibetan communities, Sichuan, China. Parasitology. 2013;140(13):1685–92.
Craig PS, Li T, Qiu J, Zhen R, Wang Q, Giraudoux P, et al. Echinococcosis and Tibetan communities. Emerg Infect Dis. 2008;14(10):1674–5.
Budke CM, Campos-Ponce M, Qian W, Torgerson PR. A canine purgation study and risk factor analysis for echinococcosis in a high endemic region of the Tibetan plateau. Vet Parasitol. 2005;127(1):43–9.
Craig PS. Epidemiology of human alveolar echinococcosis in China. Parasitol Int. 2006;55(Suppl):S221–225.
Li T, Chen X, Zhen R, Qiu J, Qiu D, Xiao N, et al. Widespread co-endemicity of human cystic and alveolar echinococcosis on the eastern Tibetan Plateau, northwest Sichuan/southeast Qinghai, China. Acta Trop. 2010;113(3):248–56.
Torgerson PR, Keller K, Magnotta M, Ragland N. The global burden of alveolar echinococcosis. PLoS Negl Trop Dis. 2010;4(6):e722.
Bowles J, Blair D, McManus DP. Genetic variants within the genus Echinococcus identified by mitochondrial DNA sequencing. Mol Biochem Parasitol. 1992;54(2):165–73.
Nakao M, Xiao N, Okamoto M, Yanagida T, Sako Y, Ito A. Geographic pattern of genetic variation in the fox tapeworm Echinococcus multilocularis. Parasitol Int. 2009;58(4):384–9.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.
McManus DP. Molecular discrimination of taeniid cestodes. Parasitol Int. 2006;55(Suppl):S31–37.
Bowles J, McManus DP. Molecular variation in Echinococcus. Acta Trop. 1993;53(3–4):291–305.
Nakao M, McManus DP, Schantz PM, Craig PS, Ito A. A molecular phylogeny of the genus Echinococcus inferred from complete mitochondrial genomes. Parasitology. 2007;134(Pt 5):713–22.
We thank Mr Tao Jiang and Mr Chun Zhang for maintaining jirds for transferring Echinococcus multilocularis in the laboratory; MS Hui Liu for providing technical support for cutting slide sections and HE staining of sections.
The work was financially supported by a grant from the National Natural Science Foundation of China (31260272).
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the GenBank database under accession numbers: KT965443–KT965494, KY290762-KY290791.
JL, WBZ and LK conceived and designed the experiments; CCW and JL performed the experiments; RB, HNF, GBP, CLZ, YMS and WZ prepared materials; JL, WBZ and JK wrote the manuscript; PG and HW critically revised the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study was conducted in accordance with the “Guidelines for the Care of Laboratory Animals” by the Ministry of Science and Technology of the People’s Republic of China (2006) for using Mongolian jirds (Meriones unguiculatus) for maintaining E. multilocularis. The jirds were maintained at the Animal Center of the First Affiliated Hospital of Xinjiang Medical University by providing water and standard chow pellets. The protocols for using parasite samples from laboratory infected jirds (Approval No IACUC-20120625003) and surgically removed liver sections of AE patients (approval No. 20141201-07) were approved by the Ethic Committee of the First Affiliated Hospital of Xinjiang Medical University. All patients participating in the project were adult and have provided written informed consents.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wu, C., Zhang, W., Ran, B. et al. Genetic variation of mitochondrial genes among Echinococcus multilocularis isolates collected in western China. Parasites Vectors 10, 265 (2017). https://doi.org/10.1186/s13071-017-2172-y
- Echinococcus multilocularis
- Mitochondrial genes
- Genetic variation
- Alveolar echinococcosis