The Scientific World Journal

Volume 2014 (2014), Article ID 867839, 8 pages

http://dx.doi.org/10.1155/2014/867839

Research Article

Genetic Diversity of Echinococcus granulosus in Southwest China Determined by the Mitochondrial NADH Dehydrogenase Subunit 2 Gene

Jiahai Wang,¹ Ning Wang,¹ Dandan Hu,¹ Xiuqin Zhong,¹ Shuxian Wang,¹ Xiaobin Gu,¹ Xuerong Peng,² and Guangyou Yang¹

¹Department of Parasitology, College of Veterinary Medicine, Sichuan Agricultural University, Ya'an 625014, China

²Department of Chemistry, College of Life and Science, Sichuan Agricultural University, Ya'an 625014, China

Received 26 August 2013; Accepted 24 November 2013; Published 30 January 2014

Academic Editors: A. Gonzalez and K. Y. Mumcuoglu

Copyright © 2014 Jiahai Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We evaluated genetic diversity and structure of Echinococcus granulosus by analyzing the complete mitochondrial NADH dehydrogenase subunit 2 (ND2) gene in 51 isolates of E. granulosus sensu stricto metacestodes collected at three locations in this region. We detected 19 haplotypes, which formed a distinct clade with the standard sheep strain (G1). Hence, all 51 isolates were identified as E. granulosus sensu stricto (G1–G3). Genetic relationships among haplotypes were not associated with geographical divisions, and fixation indices (Fst) among sampling localities were low. Hence, regional populations of E. granulosus in the southwest China are not differentiated, as gene flow among them remains high. This information is important for formulating unified region-wide prevention and control measures. We found large negative Fu’s Fs and Tajima’s D values and a unimodal mismatch distribution, indicating that the population has undergone a demographic expansion. We observed high genetic diversity among the E. granulosus s. s. isolates, indicating that the parasite population in this important bioregion is genetically robust and likely to survive and spread. The data from this study will prove valuable for future studies focusing on improving diagnosis and prevention methods and developing robust control strategies.

1. Introduction

Cystic echinococcosis (hydatid disease) is an important and globally distributed parasitic zoonosis caused by the larval stage of the cestode parasite Echinococcus granulosus complex [1]. Intermediate hosts, which include humans, sheep, goats, cattle, yak, camels, and other wild mammals, become infected by ingesting the parasite’s eggs from infected carnivores (the definitive hosts). Subsequently, a larval stage (metacestode) develops as a cyst in the internal organs (mainly in liver and lungs) of the intermediate host.

The causative agent of cystic echinococcosis was traditionally regarded to be a single species, E. granulosus. However, recent research has shown that E. granulosus is a species complex consisting of several taxa that differ in adult morphology, their preferences for intermediate hosts, and their pathogenicity to animals and humans [2]. But now, this species complex is differentiated into ten genotypes (G1–G10) [3–8]. Moreover, some researchers have suggested that E. granulosus should be classified as four species based on the substantial molecular differences in both mitochondrial and nuclear DNA genes: E. granulosus sensu stricto (genotypes G1–G3), E. equinus (genotype G4), E. ortleppi (genotype G5), and E. canadensis (genotypes G6–G10) [9], though the status of E. canadensis is still disputed [9–12]. Meanwhile, a new independent taxon named E. felidis (lion strain) was isolated from South Africa [13].

In China, cystic echinococcosis has been reported in more than twenty provinces and is particularly prevalent [14, 15]. However, to date, infections have been ascribed to just two E. granulosus genotypes; G1 (a sheep strain) and G6 (a camel strain) [16] Southwest China is one of the most serious areas of E. granulosus infections in China. The past geologic events and climate fluctuations lead to a high biodiversity of species in this area [17, 18]. In addition, E. shiquicus, a new species of Echinococcus, has been recently discovered in this region [19]. Recently, the first human CE case infected with G5 genotype (cattle strain) in Asia has been reported [20]. For these reasons, it is critical to understand the genetic composition and structure of the E. granulosus complex in this region. In this study, we provide the first investigation of the molecular diagnostics of cystic echinococcosis infections in Southwest China.

Mitochondrial DNA has been widely used in population genetics to elucidate phylogenies, as it experiences high mutation and low recombination rates and thus best reflects population genetic structure, population differentiation, and species relationships [21]. The NADH dehydrogenase subunit 2 gene (ND2 gene) evolves faster than other mitochondrial genes and is widely applied in molecular systematics and population genetics studies [22–25]. We used the ND2 gene as a genetic marker to investigate the genetic diversity and structure of Echinococcus granulosus within Southwest China. This information will be essential for further studies investigating the biology and transmission dynamics of these parasites, especially to humans, and will underpin research on the diagnosis, control, and prevention of this disease [2, 26–29].

2. Materials and Methods

2.1. Sample Collection

Larval isolates of E. granulosus () were collected from three intermediate host species (sheep, yak, and humans) in three locations in Southwest China (Qinghai, Sichuan, and Tibet; Table 1).

Table 1: Number of Echinococcus granulosus isolates obtained from the three intermediate host species in the three localities.

Locality
Host
Total

Sheep
Humans
Yaks

Qinghai
35

35

Sichuan

7
2
9

Tibet
5

2
7

2.2. PCR Amplification, Purification, and Sequencing

Total DNA was extracted using standard phenol-chloroform techniques [30] and then stored at −20°C. The complete ND2 gene was amplified using primers (P1: 5′-ATTGGACATTGTGTCTAGG-3′ and P2: 5′-GTTACTCCCATCAATGAGA-3′) that were designed based on the G1 genotype of E. granulosus (AF297617). The PCR mixture was prepared in a final volume of 25 μL containing 1 μL of template DNA, 1 μL of each primer, 12.5 μL of 2 × Taq PCR Master Mix, and 9.5 μL of the reaction buffer supplied by the CoWin Company (Beijing). Thermal cycling was performed with initial denaturation for 4 min at 94°C followed by 35 cycles of 50 s at 94°C, 45 s at 48°C, 50 s at 72°C, and a final extension of 10 min at 72°C. PCR products were sequenced three times by the Invitrogen Trading Company (Shanghai).

2.3. Data Analysis

The DNAMAN program was used to align and compare the reference sequences to the nucleotide sequences identified in our study (Table 2). The percentage divergence of nucleotide sequences was determined using MEGA 5.0 [31] applying Kimura’s two-parameter model with a -shaped parameter (alpha = 0.05) [32]. The maximum likelihood tree was also constructed using MEGA 5.0 [31] applying Kimura’s two-parameter model with 1000 bootstrap replications. Bayesian phylogenetic analyses were performed and tested using MrBayes version 3.1.2 [33]. We used TCS version 1.21 to construct networks based on the criterion of statistical parsimony [34]. The Arlequin package (version 3.5.1.2) was employed to calculated genetic diversity indices (number of haplotypes, haplotype diversity, and nucleotide diversity) [35]. Pairwise fixation indices (Fst), which estimate the degree of gene flow between two populations, were calculated with the Arlequin package. Mismatch distributions were used to test for demographic signatures of population expansions within mtDNA lineages [36]. The neutrality indices of Tajima’s D and Fu’s Fs were calculated using the population genetics package Arlequin [37, 38].

Table 2: ND2 gene sequence information for Echinococcus species.

Species

(genotype)
Location
Host
Database

accession number

E. granulosus (G1)
Britain
Sheep
AF297617

E. equinus (G4)
Britain
Horse
AF346403

E. ortleppi (G5)
Argentina
Buffalo
AB235846

E. canadensis (G6)
Kazakhstan
Camel
AB208063

E. canadensis (G7)
Poland
Pig
AB235847

E. canadensis (G8)
America
Moose
AB235848

E. multilocularis
Japan
Hamster
AB018440

3. Results

3.1. Variations in Nucleotide Sequences

All 51 isolates were successfully amplified for the complete ND2 mitochondrial gene (882 sites; similarity with E. granulosus G1 [AF297617] was 96.71% to 99.89%). Base frequencies were unequal (A = 0.1634, T = 0.5159, C = 0.0745, and G = 0.2.462), and an anti-G bias was detected, which is characteristic of the mitochondrial genome (the frequency of A + T [67.93%] was higher than the frequency of G + C [32.07%]) [39]. Overall, there were 52 point mutations, of which 34 were parsimony informative and 17 single-point mutations. We obtained 39 transition mutations and 16 transversion mutations (Table 3). No deletions and indel sites were obtained.

Table 3: Multiple alignments of the 19 mtDNA haplotypes found in this study. Haplotype identification (Hp), GenBank accession number, number of isolates (), percentage (%), and nucleotide positions (Np) of point mutations are reported.

Haplotypes
h1
h2
h3
h4
h5
h6
h7
h8
h9
h10
h11
h12
h13
h14
h15
h16
h17
h18
h19

6
14
2
5
3
3
3
1
1
1
1
1
1
2
1
1
1
2
2

%
11.8
27.5
3.9
9.8
5.9
5.9
5.9
2.0
2.0
2.0
2.0
2.0
2.0
3.9
2.0
2.0
2.0
3.9
3.9

NP

80b

C¹

90c

G¹

152b

T²

153c

G²

155b

T²

156c

C¹

157a

T¹

159c

G²

162c

T²

164b

T²

167a

T²

227b

C¹

262a

C¹

309c

G¹

339c

G¹

357c

C¹

362b
C¹

C¹

385a

414c
A²

434b

C¹

457a

A¹

460a

G¹

470b

A¹

471c

T¹

T¹

T¹

475a

C¹

487a

G²

C¹

546c

A¹

552a

C¹

578b

G¹

601a

C¹

613a

G²

618c

A¹

621c

A¹

A¹

625a

G¹

653b

G¹

700b

C¹

715a
C¹

794b

G²

801c

C¹

808a

C¹

809b

T¹

812b

G¹

813c

A¹

822c

C¹

823a

G²

834c

A¹

836b

T²

837c

G²

840c

G²

842b

G¹

843c

A²

First codon position.

bSecond codon position.

cThird codon position.

1Transition mutation.

2Transversion mutation.

3.2. Phylogenetic Analyses and Genotyping

We detected 19 mtDNA haplotypes (labeled H1 to H19; GenBank ID: KC897670-KC897688) within the 51 isolates (the localities of the haplotypes are shown in Figure 1). The maximum-likelihood phylogram clearly showed that ND2 formed one clade and was not divided into regional clades according to the allopatric distributions of the isolates (Figure 2(a)). All 19 haplotypes grouped with the G1–G3 genotype in a large clade, which was distinct from the other strains. Accordingly, all 51 isolates were classified as E. granulosus sensu stricto (genotypes G1–G3). Results from the Bayesian tree analysis depicted a similar topology to the maximum-likelihood phylogram (Figure 2(b)).

Figure 1: The localities of E. granulosus mitochondrial NADH dehydrogenase subunit 2 (ND2) haplotypes. The gray part of map is Southwest China.

(a)

(b)

Figure 2: (a) Phylogenetic relationships among the 19 haplotypes inferred by maximum-likelihood methods. (b) Phylogenetic relationships among the three populations of the E. granulosus haplotypes based on the ND2 gene inferred by Bayesian methods.

3.3. Genetic Polymorphism Analysis and Population Expansion

The overall haplotype diversity of E. granulosus in Southwest China was high (hd = 0.898 > 0.5%), although nucleotide diversity was low (Pi = 0.005 < 5%). The genetic distance between haplotypes (Kimura 2-parameters) ranged from 0.001 to 0.028 (average genetic distance = 0.007; Table 4). Results of AMOVA showed that the majority of the variation existed within regions/populations, because there was 90.73% variation within the Southwest China population and 9.27% variation of E. granulosus s. s. among Southwest China subpopulations (Table 5).

Table 4: Diversity and neutrality indices of E. granulosus sensu stricto in Southwest China calculated from nucleotide sequences derived from the ND2 gene.

Locality
Number of samples
Haplotypes
Haplotype diversity
Nucleotide diversity
Pairwise differences
Tajima’s
Fu's FS

Qinghai
35
10

−1.695*
−0.212

Sichuan
9
6

−1.422
−1.549

Tibet
7
6

−0.142
−0.550

Southwest China
51
19

−2.187***
−5.472*

, .

Table 5: Results of the analysis of molecular variance to examine genetic differences among E. granulosus populations.

Source of variation
Degrees of freedom
Sum of squares
Variance components
Percentage of variance
Fst

Interpopulation
2
8.870
0.202
9.270
0.093

Intrapopulation
48
94.688
1.972
90.730

Total variance
50
103.549
2.174

: variance component of interpopulation; : variance component of intrapopulation.

Assuming that the ancestral haplotype is still present in the population, statistical parsimony networks were constructed, in order to discern the genealogical relationship among the haplotypes. Haplotypes of different regions were mixed together and there was no relationship between haplotype affinities and their locations (Figure 3). However, the resulting network showed a star-like expansion, with one common ancestral haplotype (H2) occupying the center of the network (Figure 3). There were one to ten mutational steps between the ancestor and the other haplotypes.

Figure 3: Haplotype network for the three sampling localities of E. granulosus. Each haplotype is represented by a circle. Circle sizes are proportional to the corresponding haplotype frequencies. Small black dots represent hypothetical missing or unsampled ancestral haplotypes.

Fst values among the three sampling regions of Southwest China were low and ranged from 0.037 to 0.143 (Table 6), indicating that the geographical populations were not genetically differentiated from one another. However, based on the significantly large negative Fu’s Fs and Tajima’s values (Table 4), we can infer that the population of E. granulosus in Southwest China has undergone a demographic expansion. The unimodal distribution of the mismatch distribution supports this hypothesis of a sudden-expansion model (Figure 4).

Table 6: Pairwise fixation indices of the E. granulosus isolates in the three localities calculated from nucleotide sequences derived from the ND2 gene.

Region
Qinghai
Sichuan

Qinghai

Sichuan
0.037

Tibet
0.142*
0.116

.

Figure 4: Mismatch distributions of the ND2 haplotypes of E. granulosus in Southwest China.

4. Discussion

Southwest China is becoming a model region for biodiversity research and a major area for echinococcosis within China. This study used ND2 gene to investigate the molecular systematics of E. granulosus in this region for the first time and revealed there to be considerable genetic diversity within species complex within this region but no evidence of complete population differentiation.

We observed a distinct anti-G bias in the nucleotide sequences of the ND2 gene of E. granulosus. A similar result was found for other Echinococcus genes by Nakao et al. [40]. Variable sites occurred mainly in the third codon position, while the second codon position exhibited the least variation, supporting observations that the third codon position of mitochondrial protein genes evolves fastest, while the second codon position evolves slowest [41]. We obtained 39 transitions and 16 transversions (a transitions/transversions ratio of 2.4 > 2.0), indicating that the mutations of mitochondria ND2 gene of E. granulosus are not saturated and are suitable for the analysis of genetic variation [42].

Of the 19 haplotypes we defined (from the 51 E. granulosus isolates), one (H2) was shared between the Qinghai and Sichuan populations and two (H5 and H7) were shared between the Qinghai and Tibet populations, while no haplotypes were shared between the Sichuan and Tibet populations. The haplotype phylogenetic tree showed that E. multilocularis was genetically distinct from E. granulosus s. s. (G1, G4, G5, G6, G7, and G8) and that all 19 haplotypes (and thus all isolates) grouped with the standard sheep strain (G1) and were thus classified as E. granulosus s. s. (genotypes G1–G3). This finding is similar to previous studies that found that E. granulosus s. s. (G1–G3) is the major strains of E. granulosus throughout China [16, 43–45]. However, in this study, the G1–G3 isolates displayed distinct nucleotide differences to the reference sequences of G1–G3. Further studies are needed to ascertain the reasons for these differences.

Both the phylogenetic tree and the parsimony network showed that genetic relationships among haplotypes were not associated with geographical divisions, as haplotypes from all regions grouped together. This finding indicates that the regional populations of Southwest China are not fully differentiated from each other. The low pairwise Fst values we found are consistent with this observation of low genetic variation among the three sampling locations. Interestingly, an analogous genetic structure was found for E. granulosus in Tibet plateau by Yan et al. [45] using 28 ND1+ATP6 gene haplotypes and in eastern Tibet and Xinjiang by Nakao et al. [46] using 43 CO1 gene haplotypes. These results suggest that the population structure of E. granulosus s. s. may be highly uniform throughout China. However, a larger number of samples from more regions would need to be collected and analyzed to confirm this idea.

Overall, the haplotypes showed low nucleotide (Pi < 5%) and high haplotype diversities (Hd > 0.5) consistent with the hypothesis of sudden demographic expansion [47]. This hypothesis was supported by the finding of large negative Fu’s Fs and Tajima’s values and the single peak observed in the mismatch distribution [48]. This proposed demographic expansion may be caused by the migration of large numbers of host species (sheep, yaks, and dogs) to new areas in response to environmental change or by the artificial introduction of new hosts.

The ND2 gene similarity of the 51 E. granulosus isolates was found to be high (96.71% to 100%), though no sequence displayed 100% homology with the G1 genotype. The genetic diversity observed among the isolates in this study was significantly higher than that reported in previous studies, which used the CO1 gene to investigate genetic polymorphisms in E. granulosus in China [46]. However, the number of haplotypes identified in the previous study was much higher than that identified in the current study (43 versus 19). These differences are likely to be a consequence of the faster evolution rate of the ND2 gene than the CO1 gene of E. granulosus. Nevertheless, the high degree of genetic diversity in E. granulosus in this region revealed by all the studies indicates that the populations are genetically robust and likely to survive and spread.

In summary, our study provides the basic information about genetic diversity analysis of a wider range of isolates from different regions of Southwest China in order to understand in detail the genetic structure of E. granulosus populations and transmission dynamics of echinococcosis in these regions. Having the useful information, focus should now be directed to strengthen disease surveillance in these regions and improve diagnostic and prevention methods, in addition to developing a robust control strategy.

Conflict of Interests

The authors declare that they have no conflict of interests.

Acknowledgments

The authors wish to thank many colleagues and field workers (Zhongrong Jiang, Aiguo Yang, Shijin Deng, Li Guo, and Dawa Tsering, Weishu Kong) for their helpful assistance in collecting samples. This work was supported by Grants from the CHN National Key Technology R&D Program (no. 2006BAI06B09).

References

P. S. Craig, D. P. McManus, M. W. Lightowlers et al., “Prevention and control of cystic echinococcosis,” The Lancet Infectious Diseases, vol. 7, no. 6, pp. 385–394, 2007.

R. C. A. Thompson and D. P. McManus, “Towards a taxonomic revision of the genus Echinococcus,” Trends in Parasitology, vol. 18, no. 10, pp. 452–457, 2002.

J. Bowles, D. Blair, and D. P. McManus, “Genetic variants within the genus Echinococcus identified by mitochondrial DNA sequencing,” Molecular and Biochemical Parasitology, vol. 54, no. 2, pp. 165–173, 1992.

J. Bowles and D. P. McManus, “Molecular variation in Echinococcus,” Acta Tropica, vol. 53, no. 3-4, pp. 291–305, 1993.

R. C. A. Thompson, A. J. Lymbery, and C. C. Constantine, “Variation in Echinococcus: towards a taxonomic revision of the genus,” Advances in Parasitology, vol. 35, pp. 145–176, 1995.

J. Eckert and R. C. A. Thompson, “Intraspecific variation of Echinococcus granulosus and related species with emphasis on their infectivity to humans,” Acta Tropica, vol. 64, no. 1-2, pp. 19–34, 1997.

J. C. Scott, J. Stefaniak, Z. S. Pawlowski, and D. P. McManus, “Molecular genetic analysis of human cystic hydatid cases from Poland: identification of a new genotypic group (G9) of Echinococcus granulosus,” Parasitology, vol. 114, no. 1, pp. 37–43, 1997.

A. Lavikainen, M. J. Lehtinen, T. Meri, V. Hirvelä-Koski, and S. Meri, “Molecular genetic characterization of the Fennoscandian cervid strain, a new genotypic group (G10) of Echinococcus granulosus,” Parasitology, vol. 127, no. 3, pp. 207–215, 2003.

M. Nakao, T. Yanagida, M. Okamoto et al., “State-of-the-art Echinococcus and Taenia: phylogenetic taxonomy of human-pathogenic tapeworms and its application to molecular diagnosis,” Infection, Genetics and Evolution, vol. 10, no. 4, pp. 444–452, 2010.

R. C. A. Thompson, “The taxonomy, phylogeny and transmission of Echinococcus,” Experimental Parasitology, vol. 119, no. 4, pp. 439–446, 2008.

J. Knapp, M. Nakao, T. Yanagida et al., “Phylogenetic systematics of the genus Echinococcus (Cestoda: Taeniidae),” International journal for parasitology, vol. 43, no. 12-13, pp. 1017–1029, 2013.

D. McManus, “Current status of the genetics and molecular taxonomy of Echinococcus species,” Parasitology, vol. 140, no. 13, pp. 1617–1623, 2013.

M. Hüttner, M. Nakao, T. Wassermann et al., “Genetic characterization and phylogenetic position of Echinococcus felidis Ortlepp, 1937 (Cestoda: Taeniidae) from the African lion,” International Journal for Parasitology, vol. 38, no. 7, pp. 861–868, 2008.

P. S. Craig, “Epidemiology of human alveolar echinococcosis in China,” Parasitology International, vol. 55, pp. S221–S225, 2006.

P. Moro and P. M. Schantz, “Echinococcosis: a review,” International Journal of Infectious Diseases, vol. 13, no. 2, pp. 125–133, 2009.

L.-H. Zhang, J.-J. Chai, W. Jiao, Y. Osman, and D. P. McManus, “Mitochondrial genomic markers confirm the presence of the camel strain (G6 genotype) of Echinococcus granulosus in north-western China,” Parasitology, vol. 116, no. 1, pp. 29–33, 1998.

Q.-J. Yuan, Z.-Y. Zhang, H. Peng, and S. Ge, “Chloroplast phylogeography of Dipentodon (Dipentodontaceae) in southwest China and northern Vietnam,” Molecular Ecology, vol. 17, no. 4, pp. 1054–1065, 2008.

Z. Liu, B. Ren, F. Wei, Y. Long, Y. Hao, and M. Li, “Phylogeography and population structure of the Yunnan snub-nosed monkey (Rhinopithecus bieti) inferred from mitochondrial control region DNA sequence analysis,” Molecular Ecology, vol. 16, no. 16, pp. 3334–3349, 2007.

N. Xiao, J. Qiu, M. Nakao et al., “Echinococcus shiquicus n. sp., a taeniid cestode from Tibetan fox and plateau pika in China,” International Journal for Parasitology, vol. 35, no. 6, pp. 693–701, 2005.

M. Sharma, R. Sehgal, B. A. Fomda, A. Malhotra, and N. Malla, “Molecular characterization of Echinococcus granulosus cysts in North Indian patients: identification of G1, G3, G5 and G6 Genotypes,” PLoS Neglected Tropical Diseases, vol. 7, no. 6, Article ID e2262, 2013.

J. C. Avise, Phylogeography: the History and Formation of Species, Harvard University Press, 2000.

C. Dai, N. Zhao, W. Wang et al., “Profound climatic effects on two east Asian black-throated tits (Ave: Aegithalidae), revealed by Ecological niche models and phylogeographic analysis,” PLoS ONE, vol. 6, no. 12, Article ID e29329, 2011.

W. R. Ardren, J. Baumsteiger, and C. S. Allen, “Genetic analysis and uncertain taxonomic status of threatened Foskett Spring speckled dace,” Conservation Genetics, vol. 11, no. 4, pp. 1299–1315, 2010.

L.-F. Cao, H.-Q. Zheng, C.-Y. Hu, S.-Y. He, H.-O. Kuang, and F.-L. Hu, “Phylogeography of apis dorsata (Hymenoptera: Apidae) from China and neighboring asian areas,” Annals of the Entomological Society of America, vol. 105, no. 2, pp. 298–304, 2012.

Y.-T. Jin, R. P. Brown, and N.-F. Liu, “Cladogenesis and phylogeography of the lizard Phrynocephalus vlangalii (Agamidae) on the Tibetan plateau,” Molecular Ecology, vol. 17, no. 8, pp. 1971–1982, 2008.

A. J. Lymbery and R. C. A. Thompson, “Electrophoretic analysis of genetic variation in Echinococcus granulosus from domestic hosts in Australia,” International Journal for Parasitology, vol. 18, no. 6, pp. 803–811, 1988.

D. P. McManus, Z. Ding, and J. Bowles, “A molecular genetic survey indicates the presence of a single, homogeneous strain of Echinococcus granulosus in north-western China,” Acta Tropica, vol. 56, no. 1, pp. 7–14, 1994.

R. B. Gasser, “PCR-based technology in veterinary parasitology,” Veterinary Parasitology, vol. 84, no. 3-4, pp. 229–258, 1999.

D. P. McManus, “Molecular discrimination of taeniid cestodes,” Parasitology International, vol. 55, pp. S31–S37, 2006.

T. Maniatis, Molecular Cloning: A Laboratory Manual, Cold Spring Harbor Laboratory Press, New York, NY, USA, 1989.

K. Tamura, D. Peterson, N. Peterson, G. Stecher, M. Nei, and S. Kumar, “MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods,” Molecular Biology and Evolution, vol. 28, no. 10, pp. 2731–2739, 2011.

M. Kimura, “A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences,” Journal of Molecular Evolution, vol. 16, no. 2, pp. 111–120, 1980.

J. P. Huelsenbeck and F. Ronquist, “MRBAYES: bayesian inference of phylogenetic trees,” Bioinformatics, vol. 17, no. 8, pp. 754–755, 2001.

M. Clement, D. Posada, and K. A. Crandall, “TCS: a computer program to estimate gene genealogies,” Molecular Ecology, vol. 9, no. 10, pp. 1657–1659, 2000.

L. Excoffier and H. E. L. Lischer, “Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows,” Molecular Ecology Resources, vol. 10, no. 3, pp. 564–567, 2010.

A. R. Rogers, “Genetic evidence for a Pleistocene population explosion,” Evolution, vol. 5, no. 4, pp. 608–615, 1995.

F. Tajima, “Statistical method for testing the neutral mutation hypothesis by DNA polymorphism,” Genetics, vol. 123, no. 3, pp. 585–595, 1989.

Y.-X. Fu, “Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection,” Genetics, vol. 147, no. 2, pp. 915–925, 1997.

P. Cantatore, M. Roberti, G. Pesole et al., “Evolutionary analysis of cytochrome b sequences in some Perciformes: evidence for a slower rate of evolution than in mammals,” Journal of Molecular Evolution, vol. 39, no. 6, pp. 589–597, 1994.

M. Nakao, D. P. McManus, P. M. Schantz, P. S. Craig, and A. Ito, “A molecular phylogeny of the genus Echinococcus inferred from complete mitochondrial genomes,” Parasitology, vol. 134, no. 5, pp. 713–722, 2007.

A. Meyer, “Evolution of mitochondrial DNA in fishes,” Biochemistry and Molecular Biology of Fishes, vol. 2, no. 2, pp. 1–38, 1993.

A. Knight and D. P. Mindell, “Substitution bias, weighting of DNA sequence evolution, and the phylogenetic position of FEA's viper,” Systematic Biology, vol. 42, no. 1, pp. 18–31, 1993.

Y. R. Yang, M. C. Rosenzvit, L. H. Zhang, J. Z. Zhang, and D. P. McManus, “Molecular study of Echinococcus in west-central China,” Parasitology, vol. 131, no. 4, pp. 547–555, 2005.

S. M. Ma, S. Maillard, H. L. Zhao et al., “Assessment of Echinococcus granulosus polymorphism in Qinghai Province, People's Republic of China,” Parasitology Research, vol. 102, no. 6, pp. 1201–1206, 2008.

N. Yan, H. M. Nie, Jiang, Z. R, et al., “Genetic variability of Echinococcus granulosus from the Tibetan,” Veterinary Parasitology, vol. 196, no. 1-2, pp. 179–183, 2013.

M. Nakao, T. Li, X. Han et al., “Genetic polymorphisms of Echinococcus tapeworms in China as determined by mitochondrial and nuclear DNA sequences,” International Journal for Parasitology, vol. 40, no. 3, pp. 379–385, 2010.

B. W. Bowen and W. S. Grant, “Phylogeography of the sardines (Sardinops spp.): assessing biogeographic models and population histories in temperate upwelling zones,” Evolution, vol. 51, no. 5, pp. 1601–1610, 1997.

B. Su, Y. Fu, Y. Wang, L. Jin, and R. Chakraborty, “Genetic diversity and population history of the red panda (Ailurus fulgens) as inferred from mitochondrial DNA sequence variations,” Molecular Biology and Evolution, vol. 18, no. 6, pp. 1070–1076, 2001.