Mem Inst Oswaldo Cruz, Rio de Janeiro, VOLUME 114 | SEPTEMBER 2019
Short communication

Genetic diversity of Echinococcus vogeli in the Western Brazilian Amazon

Daniel Daipert-Garcia1, Márcio Galvão Pavan2, Leandro Batista das Neves1, Fernanda Barbosa de Almeida1, Nilton Ghiotti Siqueira3,4, Guilherme Brzoskowski dos Santos5, Tuan Pedro Dias-Correia1, Henrique Bunselmeyer Ferreira6, Rosângela Rodrigues-Silva1,+

1Fundação Oswaldo Cruz-Fiocruz, Instituto Oswaldo Cruz, Laboratório de Helmintos Parasitos de Vertebrados, Rio de Janeiro, RJ, Brasil
2Fundação Oswaldo Cruz-Fiocruz, Instituto Oswaldo Cruz, Laboratório de Mosquitos Transmissores de Hematozoários, Rio de Janeiro, RJ, Brasil
3Fundacão Hospital Estadual do Acre, Rio Branco, AC, Brasil
4Universidade Federal do Acre, Rio Branco, AC, Brasil
5Universidade Federal do Rio Grande do Sul, Centro de Biotecnologia, Laboratório de Biologia Molecular de Cestódeos, Porto Alegre, RS, Brasil
6Universidade Federal do Rio Grande do Sul, Centro de Biotecnologia, Laboratório de Genômica Estrutural e Funcional, Porto Alegre, RS, Brasil

DOI: 10.1590/0074-02760190149
519 views 205 downloads
ABSTRACT

Human polycystic echinococcosis is a parasitic infection caused by the larval stage of Echinococcus vogeli, which occurs in rural areas of Central and South America. Until now, little information on the genetic variability of E. vogeli is available.
Here, 32 samples from human-excised E. vogeli cysts had a 396-bp sequence of the mitochondrial cytochrome oxidase I (COI) gene sequenced and compared to another 17 COI sequences representing nine Echinococcus species. A Bayesian COI tree
revealed that all E. vogeli sequences formed a monophyletic and well-supported clade with an E. vogeli reference sequence. The occurrence of geographically restricted E. vogeli COI haplotypes suggests retention of ancestral polymorphisms with little migration in Acre, Brazil.

Echinococcus vogeli is the etiological agent of PE, which is endemic in the Neotropical region, including the North region of Brazil, and has a significant impact in terms of morbidity and mortality in the affected human populations.(6) Medium-sized rodents are the intermediate hosts, and dogs (wild or domestic) are the definitive hosts of E. vogeli.(7) The total number of human cases of PE reported (220 until 2015) is probably only a small fraction of the current infections, since several countries where the disease occurs do not have cumulative reporting and there is still considerable difficulty in the diagnosis and treatment of the disease, primarily in remote places.(8, 9)

Studies on molecular genetics and evolutionary ecology combined with traditional taxonomy based on morphological characters have an essential role in understanding the biodiversity of the Echinococcus genus and the differentiation of species. It is estimated that the ancestral node of the genus Echinococcus dates from ~6 million years ago. Previous reports on the genus Echinococcus evolutionary history showed that Neotropical E. vogeli and Echinococcus oligarthra (the etiological agent of unicystic echinococcosis) were the first to diverge within the clade. They reached South America after the formation of the Isthmus of Panama, along with the immigration of wild canids and felids from North America.(10)

It has been recognised that intraspecific genetic variations can influence several factors in parasites, such as life cycle patterns, host specificity, development time, transmission dynamics, sensitivity to chemotherapeutics, antigenicity, and disease-causing ability.(5) The elucidation of intraspecific genetic variation has greatly contributed to the characterisation of local populations and has allowed a better understanding of the parasite-host relationship and clinical manifestations of the disease, serving as the basis for the identification of important antigens and the development of immunodiagnostic assays and vaccines.(11)

Thirty-two human isolates of surgically excised E. vogeli cysts [Supplementary data (Table)] were obtained from PE patients located in eight municipalities of the state of Acre, Brazil [Supplementary data (Figure, Table)]. The surgical procedure was part of the treatment and all samples (cysts) were collected with the patients’ agreement, according to the ethical standards of the Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro, RJ, Brazil and to the Helsinki Declaration of 1975 (revised in 2008).

Cysts were provisionally stored at -20ºC just after the surgical procedure and defrosted for molecular analysis. For each sample, total DNA was extracted from germinal membranes using the QIAamp DNA Mini Kit (QIAGEN, Hilden, Germany), following the manufacturer’s instructions. A 396-bp sequence of the COI sequence gene was amplified by polymerase chain reaction (PCR) from each E. vogeli isolate DNA sample. PCR was performed according to Bowles et al.,(14) with modifications in reagent concentrations according to Sánchez et al.(15) and cycling conditions as described in Santos et al.(16) The amplicons were purified using the illustra™GFX™ PCR DNA kit (GE Healthcare, Little Chalfont, United Kingdom), following the manufacturer’s instructions. Both DNA strands were sequenced using the same PCR primers and the PrimTM ABI BigDye Terminator Cycle sequencing kit (Applied Biosystems, Foster City, USA), according to the manufacturers protocol. Sanger sequencing of amplicons was performed with an automated DNA sequencer ABI 3730 analyser (Applied Biosystems, Foster City, USA). Primer sequences were removed and a consensus sequence from the forward and reverse strands was assigned with SeqMan v. 7.1 (DNASTAR, Madison, USA).

Forty-nine 396-bp COI sequences (32 E. vogeli sequences from this study, and 17COI sequences from GenBank representing nine Echinococcus species (Figure) were included in the Bayesian phylogenetic tree reconstruction under the coalescent model inferred in BEAST v. 1.8.(17) Three independent runs were performed for 5 × 107 generations, sampling every 50,000 generations. Convergence of parameters and proper mixing were confirmed through the calculation of effective sample sizes (ESS) in Tracer v. 1.6;(18) ESS estimates above 104 were considered reliable.(19) The best-fit model of nucleotide substitution was determined with jModel test v. 2.(20) Molecular diversity indices of the number of segregating sites (S), the number of haplotypes (NH), haplotype diversity (HD), and nucleotide diversity (π) were computed in Dna SP v. 5(21) for each sampling site, as well as deviations from neutrality, with Fu’s Fs(22) and Tajima’s D(23) tests. Sequence divergence between populations was calculated in Mega-X.(24) A median-joining network(25) was constructed with Network v. 4.6 (Fluxus Technology Ltd. 2008) for a better visualisation of the relationships between COI haplotypes.

The Tamura and Nei model with four gamma categories (TrN+Γ) was selected as the best evolutionary model for the data, following the Akaike and Bayesian information criteria. The Bayesian COI tree (Figure A) revealed that all Echinococcus sample sequences generated in this study formed a monophyletic and well-supported clade with an E. vogeli reference sequence (KX527916), further corroborating the species identity of the collected samples. This phylogenetic reconstruction disclosed short branches (i.e.low sequence divergence) among E. vogeli samples. Only two samples from Sena Madureira and Rio Branco clustered in a separate clade (PP = 1.0), along with another four E. vogeli sequences retrieved from GenBank (AB208546, JX315616, KX257618 and NC009462).

The molecular divergence of E. vogeli samples varied between zero and 0.8%, as expected for intraspecific comparisons. The most divergent sequence was from Capixaba (0.25-0.69%) and the least divergent sequence was from Bujari (0.10-0.27%). Inspection of the sequences revealed 12 polymorphic sites and 12 haplotypes, with five haplotypes (H1-H3, H5 and H6) shared between different localities (Figure B). Molecular diversity indices (Table) showed high haplotype diversity in three localities (HD = 0.833-0.911) and nucleotide diversity (HD = 0.0025-0.0044) comparable to E. vogeli COI data obtained in a previous study (HD = 0.0027-0.0044).(16)

Haplotypes derived from the COI sequences disclosed a network (Figure B) with weak geographic structure, since localities that are close and far apart (50-650km) shared the same haplotypes. The network had two central haplotypes (H1 and H2) that were very abundant and widespread, and to which seven other less common haplotypes were closely related (1 mutational step). The most frequent haplotype (H1, N = 13) was shared with specimens from all localities, except for Capixaba. The other common haplotype (H2, N = 7) was shared with specimens from Cruzeiro do Sul, Rio Branco and Sena Madureira. These geographically restricted haplotypes suggest sudden population expansion or retention of ancestral polymorphism with little migration. Since neutrality tests did not indicate significant departures from neutrality (p > 0.05), it is highly probable that the weak geographic structure of the network reflects retention of ancestral polymorphisms and restricted gene flow in Acre, Brazil.

Previous studies on species that are transmitted in wildlife cycles, such as E. vogeli from Brazil,(16) E. multilocularis from North America and Europe,(27) and Tibetan Echinococcus shiquicus(11) reported higher nucleotide diversity (0.0044-0.0055) when compared to species transmitted to livestock animals, as Echinococcus ortleppi collected in five African countries and Brazil (0.0001-0.0008)(28) and E. granulosus s.s. from China, Peru, Eastern Europe, and Italy (0.0002-0.0055).(11) While the species with wild-life cycles seem to retain ancestral polymorphisms, those with domestic cycles show signs of genetic homogenisation due to animal transportation and introduction of small founder populations.(10)

To date, there is little information on the genetic structure of E. vogeli. The study of parasites in wild animals is laborious and sample collection in places with difficult access and with inadequate sanitary conditions lead to a great gap of knowledge concerning both the parasite and the disease. Our results differ from a previous study, in which a more comprehensive genetic structure was verified in populations from different Amazonian Brazilian states.(16) However, it is noteworthy that our sampling strategy included a more restricted geographical area. New sampling strategies are required to provide a better picture of the E. vogeli population dynamics in the Amazon.

Overall, our results showed that E. vogeli has undergone little migration in Acre and the presence of shared haplotypes among different populations seems to reflect the retention of ancestral polymorphisms. To the best of our knowledge, this is only the second study in the literature on the genetic structure and variability of E. vogeli populations. A better understanding of the dynamics, genetic structure and diversity of this species should help the diagnosis, epidemiology, and prevention of the disease, and contribute to taxonomical and evolutionary studies.

 

ACKNOWLEDGEMENTS

 

To Samanta Cristina das Chagas Xavier (Oswaldo Cruz Institute - IOC/Fiocruz-RJ) for helping prepare the map.

 

AUTHORS' CONTRIBUTION

 

DDG, HFB and RRS - Conceived the study and designed the experiments; DDG, LBN and TPDC - carried out DNA extraction from the biological specimens and performed the PCR; NGS - provided the samples; FBA, MGP and GBS - analysed the sequences and performed phylogenetic analysis; DDG, MGP, HFB and RRS - drafted the manuscript. All authors critically revised the manuscript for intellectual content and approved the final version.

REFERENCES
01. Gazzinelli A, Correa-Oliveira R, Yang GJ, Boatin BA, Kloos H. A research agenda for helminth diseases of humans: social ecology, environmental determinants, and health systems. PLoS Negl Trop Dis. 2012; 6(4): e1603-11.
02. Parija S, Giri S. A review on diagnostic and preventive aspects of cystic echinococcosis and human cysticercosis. Trop Parasitol. 2012; 2(2): 99-108.
03. Deplazes P, Rinaldi L, Rojas CAA, Torgerson PR, Harandi MF, Romig T, et al. Global distribution of alveolar and cystic echinococcosis. Adv Parasitol. 2017; 95: 315-493.
04. Bittencourt-Oliveira F, Teixeira P, Alencar A, Menezes R, Corrêa C, Neves L, et al. First parasitological, histopathological and molecular characterization of Echinococcus vogeli Rausch and Bernstein, 1972 from Cuniculus paca Linnaeus, 1766 in the Cerrado biome (Mato Grosso do Sul, Brazil). Vet Parasitol. 2018; 250: 35-9.
05. Thompson RCA. Biology and systematics of Echinococcus. Adv Parasitol. 2017; 95: 65-109.
06. D’Alessandro A, Rausch RL. New aspects of neotropical polycystic (Echinococcus vogeli) and unicystic (Echinococcus oligarthrus) echinococcosis. Clin Microbiol Rev. 2008; 21(2): 380-401.
07. Eckert J, Thompson RCA. Historical aspects of echinococcosis. Adv Parasitol. 2017; 95: 1-64.
08. Mayor P, Baquedano LE, Sánchez E, Aramburu J, Gómez-Puerta LA, Mamani VJ, et al. Polycystic echinococcosis in Pacas, Amazon region, Peru. Emerging Infect Dis. 2015; 21(3): 456-9.
09. Stijnis K, Dijkmans AC, Bart A, Brosens LAA, Muntau B, Schoen C, et al. Echinococcus vogeli in immigrant from Suriname to the Netherlands. Emerging Infect Dis. 2015; 21(3): 528-30.
10. Lymbery AJ. Phylogenetic pattern, evolutionary processes and species delimitation in the genus Echinococcus. Adv Parasitol. 2017; 95: 111-45.
11. Nakao M, Li T, Han X, Ma X, Xiao N, Qiu J, et al. Genetic polymorphisms of Echinococcus tapeworms in China as determined by mitochondrial and nuclear DNA sequences. Int J Parasitol. 2010; 40(3): 379-85.
12. Wen H, Vuitton L, Tuxun T, Li J, Vuitton DA, Zhang W, et al. Echinococcosis: advances in the 21st Century. Clin Microbiol Rev. 2019; 32(2): 1-39.
13. Siles-Lucas M, Casulli A, Conraths FJ, Müller N. Laboratory diagnosis of Echinococcus spp. in human patients and infected animals. Adv Parasitol. 2017; 96: 159-257.
14. Bowles J, Blair D, Mcmanus D. Genetic variants within the genus Echinococcus identified by mitochondrial DNA sequencing. Mol Biochem Parasitol. 1992; 54(2): 165-73.
15. Sánchez E, Cáceres O, Náquira C, Garcia D, Patiño G, Silvia H, et al. Molecular characterization of Echinococcus granulosus from Peru by sequencing of the mitochondrial cytochrome C oxidase subunit 1 gene. Mem Inst Oswaldo Cruz. 2010; 105(6): 806-10.
16. Santos GB, Soares MCP, Brito EM, Rodrigues AL, Siqueira NG, Gomes-Gouvêa MS, et al. Mitochondrial and nuclear sequence polymorphisms reveal geographic structuring in Amazonian populations of Echinococcus vogeli (Cestoda: Taeniidae). Int J Parasitol. 2012; 42(13-14): 1115-8.
17. Drummond AJ, Suchard MA, Xie D, Rambault A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012; 29(8): 1969-73.
18. Rambaut A, Suchard MA, Xie D, Drummond AJ. Tracer v1.6. 2014. Available from: http://beast.bio.ed.ac.uk/Tracer.
19. Drummond AJ, Bouckaert RR. Bayesian evolutionary analysis with BEAST. 1st ed. Cambridge, UK: Cambridge University Press; 2015. 260 pp.
20. Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012; 9(8): 772.
21. Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009; 25(11): 1451-2.
22. Fu YX. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997; 147(2): 915-25.
23. Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989; 123(3): 585-95.
24. Kumar S, Stecher G, Li M, Knyaz C, Tamura, K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018; 35(6): 1547-9.
25. Bandelt HJ, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999; 16(1): 37-48.
26. Sharma M, Fomda BA, Mazta S, Sehgal R, Singh BB, Malla N. Genetic diversity and population genetic structure analysis of Echinococcus granulosus sensu stricto complex based on mitochondrial DNA signature. PLoS One. 2013; 8(12): e82904-8.
27. Spotin A, Boufana B, Ahmadpour E, Casulli A, Mahami-Oskouei M, Rouhani S, et al. Assessment of the global pattern of genetic diversity in Echinococcus multilocularis inferred by mitochondrial DNA sequences. Vet Parasitol. 2018; 262: 30-41.
28. Addy F, Wassermann M, Banda F, Mbaya H, Aschenborn J, Aschenborn O, et al. Genetic polymorphism and population structure of Echinococcus ortleppi. Parasitology. 2017; 144(4): 450-8.

Financial support: FIOCRUZ
+ Corresponding author: rosanrsilva@gmail.com
https://orcid.org/0000-0002-6420-3400
Received 29 April 2019
Accepted 02 September 2019

Our Location

Memórias do Instituto Oswaldo Cruz

Av. Brasil 4365, Castelo Mourisco 
sala 201, Manguinhos, 21040-900 
Rio de Janeiro, RJ, Brazil

Tel.: +55-21-2562-1222

This email address is being protected from spambots. You need JavaScript enabled to view it.

Support Program

logo fiocruz logo governo
logo faperj logo cnpq marca capes