Complete mitochondrial exploration of Echinococcus multilocularis from French alveolar echinococcosis patients

Alveolar echinococcosis (AE) is a parasitosis that is expanding worldwide, including in Europe. The development of genotypic markers is essential to follow its spatiotemporal evolution. Sequencing of the commonly used mitochondrial genes cob, cox1, and nad2 shows low discriminatory power, and analysis of the microsatellite marker EmsB does not allow nucleotide sequence analysis. We aimed to develop a new method for the genotyping of Echinococcus multilocularis based on whole mitochondrial genome (mitogenome) sequencing, to determine the genetic diversity among 30 human visceral samples from French patients, and compare this method with those currently in use. Sequencing of the whole mitochondrial genome was carried out after amplification by PCR, using one uniplex and two multiplex reactions to cover the 13,738 bp of the mitogenome, combined with Illumina technology. Thirty complete mitogenome sequences were obtained from AE lesions. One showed strong identity with Asian genotypes (99.98% identity) in a patient who had travelled to China. The other 29 mitogenomes could be differentiated into 13 haplotypes, showing higher haplotype and nucleotide diversity than when using the cob, cox1, and nad2 gene sequences alone. The mitochondrial genotyping data and EmsB profiles did not overlap, probably because one method uses the mitochondrial genome and the other the nuclear genome. The pairwise fixation index (Fst) value between individuals living inside and those living outside the endemic area was high (Fst = 0.222, P = 0.002). This is consistent with the hypothesis of an expansion from historical endemic areas to peripheral regions