# Shaking the Tree: Multi-locus Sequence Typing Usurps Current Onchocercid (Filarial Nematode) Phylogeny

Emilie Lefoulon, Odile Bain, Jérôme Bourret, Kerstin Junker, Ricardo

Guerrero, Israel Cañizales, Yuriy Kuzmin, Tri Baskoro T. Satoto, Jorge

Manuel Cardenas-Callirgos, Sueli de Souza Lima, et al.

► To cite this version:

Emilie Lefoulon, Odile Bain, Jérôme Bourret, Kerstin Junker, Ricardo Guerrero, et al. Shaking the Tree: Multi-locus Sequence Typing Usurps Current Onchocercid (Filarial Nematode) Phylogeny. PLoS Neglected Tropical Diseases, 2015, 9 (11), pp.e0004233. 10.1371/journal.pntd.0004233. mnhn-02505542

HAL Id: mnhn-02505542

https://mnhn.hal.science/mnhn-02505542v1

Submitted on 4 Dec 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
RESEARCH ARTICLE

# Shaking the Tree: Multi-locus Sequence Typing Usurps Current Onchocercid (Filarial Nematode) Phylogeny

Emilie Lefoulon¹, Odile Bain¹†, Jérôme Bourret¹, Kerstin Junker², Ricardo Guerrero³, Israel Cañizales³, Yuriy Kuzmin⁴, Tri Baskoro T. Satoto⁵, Jorge Manuel Cardenas-Callirgos⁶, Sueli de Souza Lima⁷, Christian Raccurt⁸, Yasen Mutafchiev⁹, Laurent Gavotte¹⁰, Coralie Martin¹*

<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

1. Unité Molécules de Communication et Adaptation des Microorganismes, Sorbonne Universités, Muséum national d'Histoire naturelle, Paris, France,
2. ARC-Onderstepoort Veterinary Institute, Onderstepoort, South Africa,
3. Laboratorio de Biologia de Vectores y Parásitos, Instituto de Zoología y Ecología Tropical, Facultad de Ciencias, Universidad Central de Venezuela, Caracas, Venezuela,
4. Institute of Zoology, Ukrainian National Academy of Sciences, Kiev, Ukraine,
5. Center For Tropical Medicine, Universitas Gadjah Mada, Jl Teknika Utara Yogyakarta, Indonesia,
6. Asociación Peruana de Helmintología e Invertebrados A fines (APHIA), Peru,
7. LTEH ODILE BAIN, Departamento de Zoologia, Instituto de Ciências Biológicas Universidade Federal de Juiz de Fora, Minas Gerais, Brasil,
8. Laboratoire National de Santé Publique, Port-au-Prince, Haiti,
9. Institute of Biodiversity and Ecosystem Research, Bulgarian Academy of Sciences, Sofia, Bulgaria,
10. Université Montpellier, place Eugène Bataillon, Montpellier, France

† Deceased.

* cmartin@mnhn.fr

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Lefoulon E, Bain O, Bourret J, Junker K, Guerrero R, Cañizales I, et al. (2015) Shaking the Tree: Multi-locus Sequence Typing Usurps Current Onchocercid (Filarial Nematode) Phylogeny. PLoS Negl Trop Dis 9(11): e0004233. doi:10.1371/journal.pntd.0004233

**Editor:** John Pius Dalton, McGill University, CANADA

**Received:** May 5, 2015

**Accepted:** October 23, 2015

**Published:** November 20, 2015

**Copyright:** © 2015 Lefoulon et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All 345 sequences were deposited in the GenBank Data Library: accession numbers KP760116 to KP760460

**Funding:** This work was supported by the Interdisciplinary Programs of the MNHN (ATM Biodiversité actuelle et fossile. Crises, stress, restaurations et panchronisme: le message systématique and ATM Taxonomie moléculaire: DNA Barcode et gestion durable des collections) and by an European Community grant FP7-HEALTH-EPIAF-243121. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## Abstract

During the past twenty years, a number of molecular analyses have been performed to determine the evolutionary relationships of Onchocercidae, a family of filarial nematodes encompassing several species of medical or veterinary importance. However, opportunities for broad taxonomic sampling have been scarce, and analyses were based mainly on 12S rDNA and *cox1* gene sequences. While being suitable for species differentiation, these mitochondrial genes cannot be used to infer phylogenetic hypotheses at higher taxonomic levels. In the present study, 48 species, representing seven of eight subfamilies within the Onchocercidae, were sampled and sequences of seven gene loci (nuclear and mitochondrial) analysed, resulting in the hitherto largest molecular phylogenetic investigation into this family. Although our data support the current hypothesis that the Oswaldofilariinae, Waltohellinae and Icosiellinae subfamilies separated early from the remaining onchocercids, Setariinae was recovered as a well separated clade. *Dirofilaria*, *Loxodontofilaria* and *Onchocerca* constituted a strongly supported clade despite belonging to different subfamilies (Onchocercinae and Dirofilariinae). Finally, the separation between Splendidofilariinae, Dirofilariinae and Onchocercinae will have to be reconsidered.
**Competing Interests:** The authors have declared that no competing interests exist.

## Author Summary

Filariae are predominantly tissue-dwelling nematodes of the Onchocercidae (Spirurida). They are parasites of terrestrial vertebrates and some of them are particularly well known as agents of human diseases in tropical environments (e.g. onchocercosis, lymphatic filariosis, loaosis). Because of their predilection for host tissue sites, filariae are difficult to collect and their diversity is far from being well-investigated. Given the lack of fossilized material, it is not possible to formulate a comprehensive evolutionary hypothesis for this group based exclusively on morphological studies. Molecular analyses on the other hand have so far been based on a limited number of filarial species only or have used genetic markers convenient for species differentiation but not for retracing filarial history. Thus a consistent evolutionary framework for this family is still not in place. The originality and strength of the present study rest within its ability to access an unusually large range of biological samples, comprising 48 species representing seven of the eight onchocercid subfamilies. Using a new multi-gene dataset analysis, five major clades within the family were defined, including i) a group of ancestrally-derived subfamilies, ii) Setariinae as a sister group to all remaining subfamilies and iii) one large clade comprising genera of the Dirofilariinae, Onchocercinae and Splendidofilariinae.

## Introduction

The Onchocercidae (Spirurida), commonly referred to as filariae, are a family of parasitic nematodes characterised by a wide host-range in squamates, crocodilians, amphibians, mammals and birds [1]. As several filarial species are agents of human and veterinary diseases [2–4], the Onchocercidae have been the subject of numerous studies. Presently, the family is divided into eight subfamilies, including 88 genera [5–8]. They are nematodes with an evolved life cycle, involving blood- or skin-inhabiting first-stage larvae transmitted by haematophagous arthropods vectors [1,9].

A pronounced regression of morphological features and a number of convergences resulting from a parasitic life style, as well as the absence of fossilised material, make it difficult to produce phylogenetic hypotheses for this family [10]. Thus, major questions regarding the classification of the Onchocercidae, their origin and evolution remain as yet unresolved. In numerous morphological and biological studies an attempt has been made to elucidate the evolution of the Onchocercidae [1,10]. First drawn up by Wehr in 1935 [11], the classification of the family has been reconstructed several times [12,13], with the most comprehensive analysis being that of Anderson and Bain in 1976 [6].

In 1994 [14], the development of molecular techniques made the construction of the first molecular phylogeny of the Onchocercidae, based on the 5S rDNA sequences, possible. This pioneering study focused on ten species from six genera, representing two subfamilies (Dirofilariinae and Onchocercinae). With the rise of the DNA barcoding approach to distinguish between species, studies on filarial nematode phylogeny have been predominantly based on 12S rDNA and *coxI* gene analyses [15–18], and phylogenies for selected onchocercid species have been proposed. However, biological material has been scarce, impeding broad taxonomic sampling, and the markers used were not suitable to resolve the internal nodes which describe the evolution of the Onchocercidae [19]. Finally, a few studies have been conducted using alternative genes, but were either focused on nematodes in general [20–22], a specific genus [23,24] or based on the mitochondrial genome [25], all including a limited number of species.
The aim of the present study was to propose a robust phylogenetic hypothesis of the relationships within the Onchocercidae based on the concatenation of seven loci (two mitochondrial and five nuclear genes) of 48 belonging to seven subfamilies. This hitherto most comprehensive sampling of the Onchocercidae casts doubt on the biological validity of some of the classically defined subfamilies.

## Materials and Methods

### Material

Sixty-two specimens comprising 48 species belonging to 25 genera and representing seven of eight onchocercid subfamilies were analysed (Table 1 and S1 Table) [8,9,26–59]. Due to a lack of material, the Lemdaninae could not be included in this study. All procedures were conducted in compliance with the rules and regulations of the respective competent national ethical bodies (S2 Table). Nematode specimens and DNA samples were deposited in the National Nematode collection of the Muséum National d'Histoire Naturelle (MNHN), Paris, France; accession numbers are recorded in Tables 1 and S1. Below we briefly list the specimens used in this study, followed by the source/collector in parentheses. For the author(s) and year of parasite and host (authors only) species, the reader is referred to Table 1 and S1 Appendix.

DNA, L3, microfilariae.

From the following species, DNA had been extracted from adult worms for previous [18,19,41] and deposited in the MNHN collection (Table 1): *Madathamugadia hiepei*, *Rumenfilaria andersoni*, *Cercopithifilaria rugosicauda*, *Aproctella alessandroi*, *Loxodontofilaria caprini*, *Mansonella (Cutifilaria) perforata*, *Onchocerca dewittei japonica*, *Onchocerca eberhardi*, *Onchocerca skrjabini* and *Ochoterenella* sp. 1. DNA of adult *Setaria labiatopapillosa* (Dr Ben Makepeace, UK) was given to us. DNA of *Cercopithifilaria bainae* had been extracted from infective larvae from *Ixodes ricinus* (Prof Domenico Otranto, Italy). From the following species, DNA had been extracted from microfilaria: *Brugia timori* (Dr Tri Baskoro, Indonesia); *Mansonella (Mansonella) ozzardi* (Prof. Christian Raccurt, Haiti); *Dirofilaria (Dirofilaria) immitis* (Bayer Animal Health GmbH, Germany).

### Adults

From the following species, DNA had been extracted from adult females: *Acanthocheilonema viteae*, *Brugia malayi* and *Brugia pahangi* (NIAID/NIH Filariasis Research Reagent Resource Center (MTA University of Wisconsin Oshkosh—SJ 770–12); *Loa loa* (Prof Jean Dupouy Camet, Hopital Cochin, France); *Litomosoides sigmodontis* (experimentally infected jirds, *Meriones unguiculatus*; MNHN).

Additional specimens of adult filariae were recovered during dissections of their vertebrate hosts captured in the wild from different geographic areas (Table 1). They were mainly recovered from body cavities and lymphatic vessels or extracted from the dermis, subcutaneous or connective tissue, or tendons of limbs. Identification were made by several experts (OB, KJ, RG, YM, YK, SL, and AR), based on morphological studies. A few species have not yet been named, as their description is still in progress. Samples were fixed and stored in 70% ethanol. The anterior and posterior parts of worms were used for morphological studies, whereas the median part was processed for molecular analysis.

### Ethical statement

Jirds are maintained in the animal facilities of UMR7245, MNHN, as hosts of L. sigmodontis; they are inoculated intraperitoneally with 70 infective larvae and sacrificed 60 days post infection,
<table><caption>**Table 1. List of filarial species used in phylogenetic analyses in this study.** The onchocercid subfamilies are indicated in the first column. The name of the species as well as the author(s) and date are given in the second column. Columns 3 to 5 contain sample information, i.e. host species, MNHN accession number and collection locality.</caption><thead><tr><th>Subfamilies</th><th>Species</th><th>Host</th><th>N° MNHN</th><th>Collection locality</th></tr></thead><tbody><tr><td rowspan="2">Outgroups</td><td>Filaria latala Chabaud and Mohammad, 1989</td><td>Panthera leo</td><td>62YT</td><td>South Africa</td></tr><tr><td>Protospirura muricola* Gedoelst, 1916</td><td>Gorilla gorilla</td><td>97YU</td><td>Central African Republic</td></tr><tr><td rowspan="2">Oswaldofilariinae</td><td>Oswaldofilaria chabaudi Pereira, Souza and Bain, 2010</td><td>Tropidurus torquatus</td><td>191YU</td><td>Brasil</td></tr><tr><td>Oswaldofilaria petersi Bain and Sulahian 1974</td><td>Crocodilurus amazonicus</td><td>34PF</td><td>Peru</td></tr><tr><td rowspan="3">Waltonellinae</td><td>Ochoterenella sp.1 Guerrero et al. work in progress</td><td>Rhinella granulosa</td><td>200YU</td><td>Venezuela</td></tr><tr><td>Ochoterenella sp.2 Guerrero et al. work in progress</td><td>Rhinella marina</td><td>47YT</td><td>Venezuela</td></tr><tr><td>Ochoterenella sp.3 Kuzmin et al. work in progress</td><td>Phyllomedusa bicolor</td><td>194JW</td><td>French Guyana</td></tr><tr><td rowspan="2">Icosiellinae</td><td>Icosiella neglecta (Diesing, 1851)</td><td>Pelophylax ridibunda</td><td>44YT</td><td>Ukraine</td></tr><tr><td></td><td>Pelophylax kl. esculeta</td><td>45YT</td><td>France</td></tr><tr><td rowspan="2">Setariinae</td><td>Setaria labiatopapillosa (Alessandrini, 1848)</td><td>Bos taurus</td><td>413YU</td><td>Cameroon</td></tr><tr><td>Setaria tundra Bain, 1974</td><td>Rangifer tarandus</td><td>71YT</td><td>Finlande</td></tr><tr><td rowspan="6">Dirofilariinae</td><td>Dirofilaria (Dirofilaria) immitis (Leidy, 1856)</td><td>Canis familiaris</td><td>79YT</td><td>Bayer strain</td></tr><tr><td>Dirofilaria (Nochtiella) repens Railliet and Henry, 1911</td><td>Canis familiaris</td><td>297YU</td><td>Italy</td></tr><tr><td>Foleyella candeezei (Fraipont, 1882)</td><td>Agama agama</td><td>68CE</td><td>Togo</td></tr><tr><td>Loa loa (Cobbold, 1864)</td><td>Homo sapiens</td><td>80YT</td><td>France</td></tr><tr><td>Pelecitus fulicaeatrae (Diesing, 1861)</td><td>Podiceps nigricollis</td><td>49YT</td><td>Spain</td></tr><tr><td>Aproctella alessandroi Bain, Petit, Kosek and Chabaud, 1981</td><td>Saltator similis</td><td>117YU</td><td>Brasil</td></tr><tr><td rowspan="2">Splendidofilariinae</td><td>Cardiofilaria pavlovskyi Storm, 1937</td><td>Oriolus oriolus</td><td>180YU</td><td>Bulgaria</td></tr><tr><td>Madathamugadia hiepei Hering-Hagenbeck, Boomker, Petit, Killick-Kendrick and Bain, 2000</td><td>Pachycactylus turneri</td><td>81YU</td><td>South Africa</td></tr><tr><td rowspan="10">Onchocercinea</td><td>Rumenfilaria andersoni Lankester and Snider, 1982</td><td>Rangifer tarandus</td><td>94YU</td><td>Finlande</td></tr><tr><td>Acanthocheilonema odendhali (Perry, 1967)</td><td>Callorhinus ursinus</td><td>401YU</td><td>Alaska</td></tr><tr><td>Acanthocheilonema viteae (Krepkogorskaya, 1933)</td><td>Meriones unguiculatus</td><td>7YT</td><td>FR3 strain</td></tr><tr><td>Breinlia (Breinlia) jittapalapongi Veciana, Bain, Morand, Chaisiri, Douanghoupha, Miquel and Ribas, 2015</td><td>Rattus tanezumi</td><td>78YT</td><td>Laos</td></tr><tr><td>Brugia malayi (Brug, 1927)</td><td>Meriones unguiculatus</td><td>8YT</td><td>FR3 strain</td></tr><tr><td>Brugia pahangi (Buckley and Edeson, 1956)</td><td>Meriones unguiculatus</td><td>46YT</td><td>FR3 strain</td></tr><tr><td>Brugia timori Partono, 1977</td><td>Homo sapiens</td><td>6YT</td><td>Indonesia</td></tr><tr><td>Cercopithifilaria bainae Almeida and Vicente, 1984</td><td>Canis familiaris</td><td>9YT</td><td>experimental</td></tr><tr><td>Cercopithifilaria rugosicauda (Böhm and Supperer, 1953)</td><td>Capreolus capreolus</td><td>350YU</td><td>France</td></tr><tr><td>Cruorifilaria tuberocauda Eberhard, Morales and Orihel, 1976</td><td>Hydrochoerus hydrochaeris</td><td>55YT</td><td>Venezuela</td></tr><tr><td rowspan="2">Dipetalonema caudispina (Molin, 1858)</td><td>Ateles paniscus</td><td>362YU</td><td>Guyana</td></tr><tr><td>Ateles sp.</td><td>64YT</td><td>Guyana</td></tr><tr><td rowspan="3">Dipetalonema gracile (Rudolphi, 1809)</td><td>Cebus olivaceus</td><td>124CV</td><td>Venezuela</td></tr><tr><td>Cebus apella</td><td>215YU</td><td>Peru</td></tr><tr><td>Ateles sp.</td><td>63YT</td><td>Guyana</td></tr><tr><td>Dipetalonema graciliformis (Freitas, 1964)</td><td>Saimiri scuireus</td><td>220YU</td><td>Peru</td></tr><tr><td>Dipetalonema robini Petit, Bain and Roussilhon, 1985</td><td>Lagothrix poeppigii</td><td>216YU</td><td>Peru</td></tr><tr><td>Litomosoides brasiliensis Lins de Almeida, 1936</td><td>Carollia perspicillata</td><td>35/37PF</td><td>Peru</td></tr><tr><td>Litomosoides hamletti Sandground, 1934</td><td>glossophaga soricina</td><td>36PF1</td><td>Peru</td></tr><tr><td>Litomosoides sigmodontis Chandler, 1931</td><td>Meriones unguiculatus</td><td>186MS</td><td>MNHN strain</td></tr><tr><td>Litomosoides solarii Guerrero, Martin, Gardner and Bain, 2002</td><td>Trachops cirrhosus</td><td>213YU</td><td>Venezuela</td></tr><tr><td>Loxodontofilaria caprini Uni and Bain, 2006</td><td>Naemorhedus crispus</td><td>YG2-58</td><td>Japan</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1. (Continued)**</caption><thead><tr><th>Subfamilies</th><th>Species</th><th>Host</th><th>N° MNHN</th><th>Collection locality</th></tr></thead><tbody><tr><td></td><td>Mansonella (Cutifilaria) perforata Uni, Bain and Takaoka, 2004</td><td>Cervus nippon</td><td>216JW</td><td>Japan</td></tr><tr><td></td><td>Mansonella (Mansonella) ozzardi (Manson, 1897)</td><td>Homo sapiens</td><td>77YT</td><td>Haiti</td></tr><tr><td></td><td>Monanema martini Bain, Bartlett and Petit, 1986</td><td>Arvicanthis niloticus</td><td>31NC</td><td>Senegal</td></tr><tr><td></td><td>Onchocerca armillata Railliet and Henry, 1909</td><td>Bos taurus</td><td>54FK</td><td>Cameroon</td></tr><tr><td></td><td>Onchocerca dewittei japonica Uni, Bain and Takaoka, 2001</td><td>Sus scrofa<br>leucomystax</td><td>OB9</td><td>Japan</td></tr><tr><td></td><td>Onchocerca eberhardi Uni and Bain, 2007</td><td>Cervus nippon</td><td>S63-5</td><td>Japan</td></tr><tr><td></td><td>Onchocerca gutturosa Neumann, 1910</td><td>Bos taurus</td><td>54FK</td><td>Cameroon</td></tr><tr><td></td><td>Onchocerca ochengi Bwangamoi, 1969</td><td>Bos taurus</td><td>54FK</td><td>Cameroon</td></tr><tr><td></td><td>Onchocerca skrjabini Ruklyadev, 1964</td><td>Cervus nippon</td><td>S63-6</td><td>Japan</td></tr><tr><td></td><td>Yatesia hydrochoerus (Yates, 1980)</td><td>Hydrochoerus<br>hydrochaeris</td><td>52YT</td><td>Venezuela</td></tr></tbody></table>

doi:10.1371/journal.pntd.0004233.t001

upon which adult worms are recovered from the peritoneal cavity. All experimental procedures were carried out in strict accordance with the EU Directive 2010/63/UE and the relevant national legislation, namely the French "Décret no 2013-118, 1er février 2013, Ministère de l'Agriculture, de l'Agroalimentaire et de la Forêt", National licence number 75-1415 approved animal experiments: protocols were approved by the ethical committee of the Museum National d'Histoire Naturelle (Comité Cuvier, Licence: 68-002) and by the "Direction départementale de la cohésion sociale et de la protection des populations" (DDCSPP) (No. C75-05-15).

Some non-human vertebrates were captured for experimental procedures, subject to the ethics approval of the relevant national bodies (S2 Table), while others were obtained at abattoirs or donated to the MNHN by hunters or veterinarians (S2 Table). The MNHN does neither solicit nor compensate for these donations.

Human samples were not collected specifically for this study. They were provided by third parties. Human blood samples were collected with the ethics approval of the relevant national bodies (S2 Table).

The adult female of *Loa loa* had been surgically removed from a male patient in Hospital Cochin, France. The adult subject had given oral consent to Prof J. Dupouy Camet, head of parasitology department at Hospital Cochin, to donate the worm to the MNHN (S2 Table). The MNHN did not and will not compensate for the donation.

## Molecular screening

DNA was extracted using a commercial kit, following manufacturer's instructions (QIAamp micro kit, Qiagen, Germany). A preliminary step of disruption for two cycles of 30 seconds at 30 Hz using a TissueLyser II (Qiagen, Germany) was added. PCR screening of filarial nematodes was based on seven partial sequences of seven different genes: two mitochondrial genes, 12S rDNA (approximately 450 base-pair (bp) sequence) and cytochrome oxidase subunit I (*coxI*; approximately 600 bp); five nuclear genes, 18S rDNA (approximately 740 bp), 28S rDNA (approximately 900 bp), the myosin heavy chain (*MyoHC*; approximately 785 bp), RNA polymerase II large subunit (*rbp1*; approximately 640 bp), 70 kilodalton heat shock proteins (*hsp70*; approximately 610 bp). Amplification of the 12S rDNA and *coxI* sequences was conducted according to Casiraghi et al. [15]. For the remaining five genes, primer pairs were designed (S3

<footnote>*Protospirura muricola is typically a parasite of rodents, but has been reported from primates [78].</footnote>
Table) based on regions that had been found conserved among nine species for which drafts of or complete genome were accessible: from NCBI database *Brugia malayi* (PRJNA27801), *Loa loa* (PRJNA37757), *Onchocerca flexuosa* (Weld, 1856) (PRJEB512), *Dirofilaria immitis* (PRJEB1797) *Onchocerca ochengi* (PRJEB1809); from Nematode Genomes from the Blaxter lab, University of Edinburg (www.nematode.org) *Acanthocheilonema viteae* (nAv.1.0), *Litomosoides sigmodontis* (nLs.2.1); and from Filarial worms Sequencing Project, Broad Institute of Harvard and MIT (http://www.broadinstitute.org/) *Onchocerca volvulus* (Leuckart, 1893) and *Wuchereria bancrofti* (Cobbold, 1877). For each gene, primer pairs adapted to nested PCR were designed (S3 Table). PCRs were processed in a final volume of 20 µl under the conditions summarised in S3 Table. Obtained PCR products were purified using the SV Wizard PCR Purification Kit (Promega, USA) and sequenced directly. A total of 345 sequences were deposited in the GenBank Data Library: KP760116 to KP760460 (S4 Table)

## Phylogenetic analyses

Sequences generated during the current study and previously published sequences from draft/ complete genomes were aligned using MAFFT [60]. To check for the absence of stop codons, the alignment of coding genes was translated using EMBOSS Transeq [61], and a comparison with available transcript sequences made. Using the corrected version of the Akaike Information Criterion (AICc), JModelTest analysis [62] was performed to establish the evolution model best adapted to the sequences alignment for each individual gene and for the concatenation of all genes. The General time-reversible plus Invariant sites plus Gamma distributed model (GTR+I+Γ) offered the best fit for *coxI*, *hsp70*, 12S rRNA and 28S rRNA sequences as well as the concatenated alignment; the General Time-Reversible plus Gamma distributed model (GTR+Γ) for *MyoHC* and *rbp1* sequences alignment; and the Kimura 2-parameter Invariant sites plus Gamma distributed model (K80+I+Γ) for 18S rRNA sequences alignment. To root the trees, two species were included as outgroups: *Filaria latala* (Spirurida: Filariidae) and *Protospirura muricola* (Spirurida: Spiruridae). Phylogenetic relationships between onchocercid taxa, based on the concatenated dataset, were performed by Bayesian inference using MrBayes [63]. A partitioned model was implemented to estimate evolution parameters separately for each gene. Two runs were performed using five millions steps with four chains, with tree sampling every 1,000 generations; the first 1,250 points were discarded as burn-in and Posterior Probabilities were calculated from these post-burning trees. In addition, Maximum Likelihood (ML) was used to infer a phylogenetic tree based on the partitioned concatenated dataset, and was executed with 1000 slow bootstrap replicates using RaxML [64]; presented in Supplementary data (S1 Fig).

## Accession numbers

All new sequences generated in this study were deposited in GenBank (http://www.ncbi.nlm.nih.gov/genbank) under the accession numbers: KP760168 to KP760211 for *coxI*, KP760314 to KP760357 for 12S rDNA, KP760116 to KP760167 for 18S rRNA, KP760212 to KP760262 for *MyoHC*, KP760410 to KP760460 for *hsp70*, KP760263 to KP760313 for *rbp1*, KP760358 to KP760409 for 28S rRNA.

## Results/Discussion

### Choosing an outgroup

An important consideration when proposing a comprehensive phylogeny for any taxonomic group is the suitable rooting of the tree and, thus, the choice of appropriate outgroups. Previous studies of the Onchocercidae, based on *coxI* and 12S rRNA gene analyses, have used spirurid
representatives of the Thelaziidae [15,16] and Filariidae [16] as outgroups. Furthermore, some phylogenetic analyses of nematodes based on SSU rDNA [21] or complete mitochondrial genome sequences [22] suggested other spirurids as the Physalopteroidea or Diplotriaenoidea to be appropriate outgroups for the Onchocercidae. In the present study, *P. muricola* (Spiruridae) and *F. latala* (Filariidae) were investigated as potential outgroups. Our multi-gene analyses indicate that this choice is reasonable (Fig 1 and S1 Fig). To support this result, the average genetic divergence between the outgroup and ingroup taxa was about 20% (21% for *P. muricola*; 25% for *F. latala*), and only about 14% between onchocercid specimens themselves. It

<figure><img src="image_4.png" /></figure>

Subfamilies :

<table><caption>**Fig 1. Phylogeny of Onchocercidae based on partitioned concatenated datasets of 12S rDNA, cox1, rbp1, hsp70, myoHC, 18S rDNA, and 28S rDNA sequences using Bayesian Inference.** The total length of datasets is approximately 4950 bp. Sixty onchocercid specimens (representing 48 species) were analysed. *Filaria latala* and *Protospirura muricola* were used as outgroups. The topology was inferred using Bayesian (B) inference. Nodes are associated with Bayesian posterior probabilities based on one run of 5 million generations. The onchocercid subfamilies are indicated by colour: blue for Onchocercinae, dark green for Dirofilariinae, purple for Splendidofilariinae, pale green for Setariinae, yellow for Waltonellinae, orange for Icosiellinae and red for Oswaldofilariinae.</caption><tr><td>Oswaldofilariinae</td><td>Dirofilariinae</td></tr><tr><td>Icosiellinae</td><td>Onchocercinae</td></tr><tr><td>Waltonellinae</td><td>Splendidofilariinae</td></tr><tr><td>Setariinae</td><td></td></tr></table>

doi:10.1371/journal.pntd.0004233.g001
was, therefore, concluded that both *P. muricola* and *F. latala* are appropriate outgroups for the present multi-gene dataset. In congruence with previous systematic analyses [5,6], the chosen rooting clearly supported the monophyly of the Onchocercidae.

## Monophyly of the Onchocercidae

All onchocercid taxa examined in the present study formed a monophyletic group and five strongly supported clades within the family were identified (Fig 1). Based on their monophyly and to facilitate the description of tree topology, these clades will be referred to as ONC1 to ONC5 (Fig 2).

<figure><img src="image_5.png" /></figure>

Subfamilies :

<table><tr><td>Oswaldofilariinae</td><td>Dirofilariinae</td></tr><tr><td>Icosiellinae</td><td>Onchocercinae</td></tr><tr><td>Waltonellinae</td><td>Splendidofilariinae</td></tr><tr><td>Setariinae</td><td></td></tr></table>

**Fig 2.** Onchocercid clades based on partitioned concatenated datasets of 12S rDNA, *cox1*, *rbp1*, *hsp70*, *myoHC*, 18S rDNA, and 28S rDNA sequences using Bayesian Inference. The total length of datasets is approximately 4950 bp. Sixty onchocercid specimens (representing 48 species) were analysed. *Filaria lata*la and *Protospirura muricola* were used as outgroups. Bayesian posterior probabilities support five clades named ONC1 to ONC5. Dotted areas represented species of the same genus. The topology was inferred using Bayesian (B) inference on one run of 5 million generations. The onchocercid subfamilies are indicated by colour: blue for Onchocercinae, dark green for Dirofilariinae, purple for Splendidofilariinae, pale green for Setariinae, yellow for Waltonellinae, orange for Icosiellinae and red for Oswaldofilariinae.

doi:10.1371/journal.pntd.0004233.g002
## The ancestrally derived subfamilies are comforted

In the present phylogenetic analysis, the Oswaldofilariinae, Icosiellinae and Waltonellinae occupied a commun position as a highly supported clade (ONC1) and sister group of all other taxa studied within the Onchocercidae (Fig 1.). The three subfamilies appeared to be closely related, with ONC1 having evolved independently from the remaining onchocercid representatives (Fig 2., S1 Fig). Previous morphological studies revealed a number of ancestral characters in the Icosiellinae and Oswaldofilariinae, such as infective larvae possessing cephalic spines (Icosiellinae), and a vulva that is positioned posterior to the oesophagus (Oswaldofilariinae) [10,65]. Therefore, these subfamilies have traditionally been considered “ancient” [1,10,65]. In addition, it has been suggested that their Gondwanian distribution is an indication of an ancient evolution which occurred prior to the supercontinent’s break-up [1,10]. Interestingly, the Oswaldofilariinae appear to have evolved independently from the remaining two subfamilies of clade ONC1 and were placed as sister group to a clade that consists of *Ochoterenella* spp. (Waltonellinae) and *Icosiella neglecta* (Icosiellinae) (S1 Fig). Indeed, previous studies emphasised the presence of numerous plesiomorphic characters (long oesophagus, large buccal capsule and the presence of deirids) within *Oswaldofilaria* spp. [10]. Moreover, the host spectrum of these subfamilies differs: oswaldofilariines are parasites of squamates, whereas members of Icosiellinae and Waltonellinae are exclusively parasites of anurans [5].

## Setariinae are unusual Onchocercidae with an intermediate position

The systematic position of the Setariinae among the Filarioidea superfamily has long been questioned. Initially the Setariinae had been placed into Filariidae because they share certain morphological similarities with representatives of this family (spines on anterior end of first stage larva and adult cuticular elevation) [11,13]. Nevertheless, based on other biological and morphological traits (e.g. viviparous, long microfilaria with pointed spineless tail), it would have been equally possible to consign them to the Onchocercidae. At a later stage, the Setariidae was proposed as a family with the purpose of removing the Setariinae from the Filariidae [12], and current classifications list the Setariinae as a subfamily within the Onchocercidae [6]. Such an affiliation between the Setariinae and Onchocercidae has since been supported by many molecular phylogenies [16,19,25,65]. The same was true for the present study, where the two *Setaria* species grouped in one clade (ONC2) with good support (Fig 2). Excepting ONC1, the clade consisting of *Setaria* spp, was sister to all remaining onchocercid taxa sampled in this study. In the past, the Setariinae, Oswaldofilariinae and Waltonellinae were grouped together on the previous phylogenetic trees despite an unresolved topology [16,65]. For some time, morphological studies of infective stages from *Setaria* spp. (long larvae, presence of deirids) have supported the hypothesis that Setariinae may have emerged early in the evolution of Onchocercidae [10,59]. In contrast, the study of larval development of *S. labiotopapillosa* suggested that Setariinae, like the remaining Onchocercidae, could have been derived from habronematid ancestors, but that speciation leading to the Setariinae occurred later and independently from the remaining Onchocercidae [1,66]. For example, the development of the glandular cells of the oesophagus in the Setariinae was shown to be closed to the development observed in some potential ancestors of the Onchocercidae, namely the Habronematidae [66,67]. The current phylogeny strongly supported the position of Setariinae as a group derived from an independent speciation to the examined Onchocercidae.

## Subfamilies within the Onchocercidae: A need for systematic revision

The remaining taxa of the Onchocercidae analysed in this study grouped in three distinct and well supported clades (ONC3, ONC4 and ONC5) (Fig 2), including representatives of the
Dirofilariinae, Onchocercinae and Splendidofilariinae. It is noteworthy that the tree topology
demonstrated here is not congruent with classic systematic delineations:

I. The first clade, ONC3 (Fig 2), contained representatives of the genera *Dirofilaria* (Dirofilariinae), *Loxodontofilaria* and *Onchocerca* (Onchocercinae) as suggested by several earlier molecular analyses [14,16,18,25]. Although these genera belong to two different subfamilies when looking at established classifications [13,68], a possible need to merge these subfamilies is reinforced by several biological traits [65]. Firstly, the infective stages of *Dirofilaria* spp. and *Onchocerca* spp. possess morphological similarities (normally developed buccal capsule; cylindrical tail with tiny caudal lappets) [59,65]. Secondly, there are shared life history traits: an early first moult in the vertebrate host takes place early, *i.e.* 2–3 days postinfection, in *Dirofilaria* spp. and *Onchocerca* spp. [1], as opposed to a late moult, *i.e.* 6–10 days postinfection, in all other onchocercid species studied thus far [68]. No information is currently available on the infective larvae and development in the final host of representatives of *Loxodontofilaria*. This clade (ONC3) (Fig 2) constitutes the most ancestrally derived group belonging to modern Onchocercidae. Clade ONC3 was sister group of both clades ONC4 and ONC5 (Fig 2), suggesting an independent evolution.

II. The second clade, ONC4 (Fig 2), was composed of the genera *Dipetalonema*, *Acanthocheilonema*, *Monanema*, *Cruorifilaria*, *Cercopithifilaria*, *Yatesia* and *Litomosoides*. For a long time, taxonomists used *Dipetalonema* as a “catch all” genus [69]. In years to follow, numerous subgenera formerly included in this genus were elevated to generic rank, whereas other species were transferred in newly erected genera; together, they were referred to as the “*Dipetalonema* lineage” [69,70]. This lineage combined those species of the Onchocercinae that present adults with a long tail, a buccal capsule divided into two segments (or three in *Skrjabinofilaria* spp.), and a caudal extremity with two subterminal lappets [70]. These structures might, however, be more or less atrophied in the more specialised forms, such as seen in the buccal capsule of *Monanema* and *Cercopithifilaria* species [69,71]. The present phylogenetic tree supports the accuracy of the *Dipetalonema* lineage, as all genera included in the clade ONC4, with the exception of *Litomosoides*, have been assigned to this lineage [69,70]. Prior to this study, molecular analyses had only supported the clustering of *Dipetalonema* with *Acanthocheilonema* [16]. A closer look at clade ONC4 revealed its division into two groups (Fig 2):

i) *Dipetalonema* spp. were sister to all remaining taxa within clade ONC4 (Fig 2). This position agrees with earlier hypotheses that, based on ancestral morphological particularities (divided oesophagus, well developed and highly cuticularized buccal capsule) the genus *Dipetalonema* should be considered distinct from the remaining genera in the *Dipetalonema* lineage [70]. Moreover, it had been proposed that the genus may have been derived through geographic isolation in the Neotropical region after the Gondwana fragmentation and so would be a paleoendemic South American genus [10,70].

ii) The second subclade group comprised members of *Acanthocheilonema*, *Monanema*, *Cercopithifilaria*, *Yatesia*, *Cruorifilaria* and *Litomosoides* seems more derived. However, species of *Acanthocheilonema* and *Monanema* seemed to be more closely related and were sister to the remaining taxa in this group. Earlier taxonomic studies described *Acanthocheilonema* as an ancestral genus due to several plesiomorph characters of its representatives (strongly divided oesophagus, highly cuticularized buccal capsule, shape of right spicule of male) [10,65,70]. In contrast, the morphology of *Monanema* includes both ancestral (numerous caudal papillae distributed along length of tail) as well as derived characters (undivided oesophagus) [69,71]. Bain *et al*. (70) considered *Monanema* as
derived from Ethiopian *Acanthocheilonema* species [70]. The remaining genera formed two sister groups: the one uniting taxa of *Cruorifilaria*, *Yatesia* and *Cercopithifilaria*, the other containing species of *Litomosoides*. Earlier studies had suggested *Cercopithifilaria* to be derived from *Acanthocheilonema* and, as indicated by the derived morphological characters of its members (reduced buccal capsule, undivided oesophagus, shape of right spicule of male, number and repartition of caudal papillae), to represent the most evolved forms of the *Dipetalonema* lineage [10,70]. *Litomosoides* had not been included in the *Dipetalonema* lineage in the 1976 classification by Chabaud and Bain [69]. The latter authors did, however, point out that the reason for this omission was somewhat arbitrary. In fact, *Litomosoides* species do share morphological similarities with other representatives of the *Dipetalonema* lineage (long tail, long buccal capsule) [6,14,69].

### III. The third clade (ONC5) united genera belonging to the Dirofilariinae (*Loa, Foleyella, Pelecitus*), Onchocercinae (*Wuchereria, Brugia, Breinlia, Mansonella*) and Splendidofilariinae (*Aproctella, Cardiofilaria, Madathamugadia, Rumenfilaria*) (Fig 2). Most human filariae were included in this clade, with the exception of *O. volvulus*, which groups in clade ONC3. Three groups could be pinpointed in clade ONC5:

i) The first group, including *Brugia* and *Wuchereria*, was sister to all remaining representatives of clade ONC5. Previous molecular studies highlighted a close relation between these two genera [14,16,25,65], which is supported by shared morphological similarities (reduced buccal capsule, glandular part of oesophagus, vagina with bends and chamber in adults; infective larvae with a long tail and three lappets) [1,72]. Moreover, their mode of transmission and development is similar (culicid vector, lymphatic location in the host) [9].

ii) The second group comprised species exclusively assigned to the Dirofilariinae: *L. loa*, *P. fulicaeatrae* and *F. candeezei*. These parasites did not group together in previous molecular analyses [16,19]. However, they share morphological similarities, based on which this subfamily had originally been introduced (highly developed caudal alae of adult males) [9,12]. Furthermore, *L. loa* and *Pelecitus scapiceps* Leidy, 1886 present similar biological traits, such as a late third moult (respectively on day 8 and day 6 postinfection) and a fourth moultprior to day 20 postinfection, both of which are close to the developmental pattern seen in many species of the Onchocercinae and Splendidofilariinae [31,68].

iii) The last group in clade ONC5 was composed of representatives of Splendidofilariinae and Onchocercinae. Up to now, Splendidofilariinae have rarely been studied [19,25]. In the present analysis, species of *Madathamugadia* and *Rumenfilaria* (Splendidofilariinae) appeared more closely related to the two species of *Mansonella* (Onchocercinae), of which they were a sister group, than they were to the remaining splendidofilariine genera, *Cardiofilaria* and *Aproctella* species. The close relationship between *Cardiofilaria* and *Aproctella* has been underlined by morphological similarities in infective larvae (body length exceeding 1300µm, thick and short tail, ending in two small lateral terminal papillae) [33]. *Mansonella* is a complex and taxonomically challenging genus and numerous attempts to reorganise its classification, including description of several subgenera, have been made throughout the years [1,69,73]. Its relationship with representatives of the Splendidofilariinae has not been resolved, but it does share similarities with some splendidofilariine genera such as the *Madathamugadia* (e.g short and thin infective larvae with long tail) [35,59]. Species of *Breinlia* have been described as Australian representatives of the *Dipetalonema* lineage which have more recently migrated into Asia, and are morphologically distinct from the remaining species of this lineage [70]. Their somewhat isolated position
was upheld by the current phylogeny. It is interesting to note that these various species (e. g. *Madathamugadia* or *Foleyella*) have been described as an evolved group [74].

One of the main difficulties in comparing the results of the present study with those of earlier molecular analyses was either the lack of broad taxonomic sampling or incompletely resolved deeper phylogenetic relationships in preceding studies [15,16,18,19]. However, an equivalent to clade ONC5 clade had been identified before, in the first phylogeny including *L. loa*, *B. malayi* and *W. bancrofti* and *Mansonella (M.) perstans* (Manson, 1891) [14]. A similar clade, including *B. malayi*, *W. bancrofti*, *L. loa* and *Chanderella quiscali* Linstow, 1904, a Splendidofilariinae, was confirmed by a recent analysis of mitochondrial genomes [25].

To summarise, the present phylogeny did neither support the monophyly of the Dirofilariinae, Onchocercinae nor Splendidofilariinae (see ONC3 and ONC5, Fig 2). Interestingly, however, the relationships exhibited between the various taxa in clades ONC3 and ONC5, although traditionally included in different subfamilies, are supported by morphological as well as biological features (see discussion above). Based on the combined analysis of molecular, structural and developmental characters, we, thus, conclude that a reassessment of the boundaries between these three subfamilies is called for.

## Evolution of filariae, their vector and host: Is there a relationship?

The parasitic life style of the Onchocercidae has led to the reduction of certain morphological traits and promoted convergent evolution in others [10], reducing the usefulness of structural characters in developing a phylogenetic framework for this family. One might thus look at the role of vectors and hosts for further clues as to the speciation of the Onchocercidae. Three (ONC1, ONC2 and ONC3) of the five clades exhibited on the current phylogenetic tree use exclusively dipterans as vectors. In addition, dipterans are involved in the transmission of members of clades ONC4 and ONC5, supporting a hypothesis that old evolutionary bonds exist between this order and the Onchocercidae (Fig 3, S1 Table). The diversity of vectors used in clades ONC4 and ONC5 (Fig 3), on the other hand, suggests that the evolution of filarial nematodes is not strictly constrained by their vector but does permit changes to other, more appropriate groups of vectors, to ensure parasite transmission, colonisation of new ecological niches and diversification. Clade ONC4 shows an adaptation to a new group of vectors-ticks and mites. It should, however, be emphasised that some *Acanthocheilonema* spp. (ONC4) are known to use fleas, hippoboscid flies or sucking lice as vectors [9], and *Dipetalonema* spp. are transmitted by ceratopogonid biting midges (Fig 3). Considering relationships of the Onchocercidae with their definitive hosts, the present analysis revealed a group of parasites adapted exclusively to cold blooded tetrapods, namely squamates, crocodilians and amphibians (ONC1), while their sister group diversified in various groups of mammals (ONC2-5), as well as birds and squamates (ONC5) (Fig 4). It has been suggested that the parasite-host relationships of *M. hiepei* (parasite of the gecko *Chondrodactylus turneri*) and *F. candeezi* (parasite of lizards) may be the result of secondary capture with the narrow geographical localisation of *Foleyella* and *Madathamugadia*, i.e. Africa and Madagascar, supporting the hypothesis of capture phenomena [74]. The present phylogenetic analysis does not support a close relationship between parasites of birds, but representation of avian parasites was limited. All three parasites of birds formed part of clade ONC5, which presented the most heterogeneous host range within the Onchocercidae, comprising parasites of three classes of vertebrates, reptiles, birds and mammals. The diverse mammalian hosts represented in clades ONC3-5 indicate multiple events of host-switching and radiation within the evolutionary history of the group. In summary, the overall picture revealed by the current analysis does not suggest a broad pattern of onchocercid parasite-host/vector coevolution and is consistent with earlier studies suggesting
<figure><img src="image_6.png" /><figcaption>**Fig 3.** Phylogeny of Onchocercidae based on partitioned concatenated datasets of 12S rDNA, *cox1*, *rbp1*, *hsp70*, *myoHC*, 18S rDNA, and 28S rDNA sequences using Bayesian Inference with indication of vector range. The total length of datasets is approximately 4950 bp. Sixty onchocercid specimens (representing 48 species) were analysed. *Filaria latala* and *Protospirura muricola* were used as outgroups. The topology was inferred using Bayesian (B) inference. Nodes are associated with Bayesian posterior probabilities based on one run of 5 million generations. The vector family for each filariae species is indicated to the right of the tree using the specified symbols.</figcaption></figure>

doi:10.1371/journal.pntd.0004233.g003

that parasitic nematode speciation could be related to a variety of events of host switching and host acquisition and/or geographical and ecological drift [75–77].

## Concluding remarks

The present study casts doubt on the validity of the currently eight traditional onchocercid subfamilies. Traditional classifications of the Onchocercidae have been established mainly on
<figure><img src="image_7.png" /><figcaption>**Fig 4.** Phylogeny of Onchocercidae based on partitioned concatenated datasets of 12S rDNA, *cox1*, *rbp1*, *hsp70*, *myoHC*, 18S rDNA, and 28S rDNA sequences using Bayesian Inference with indication of host vertebrate range. The total length of datasets is approximately 4950 bp. Sixty onchocercid specimens (representing 48 species) were analysed. *Filaria latala* and *Protospirura muricola* were used as outgroups. The topology was inferred using Bayesian (B) inference. Nodes are associated with Bayesian posterior probabilities based on one run of 5 million generations. The host vertebrate order for each filariae species is indicated to the right of the tree using the specified symbols.</figcaption></figure>

doi:10.1371/journal.pntd.0004233.g004

the basis of morphological characters of adult worms. However, as suggested by previous authors [10,68], biological features, such as larval biology or the morphology of the third-stage larva, may offer a better phylogenetic resolution. It is hoped that an approach combining broad
molecular taxonomic sampling and traditional morphological as well as life history studies, will eventually lead to the development of more comprehensive phylogenetic hypotheses for this fascinating group of nematodes.

## Supporting Information

S1 Fig. Phylogeny of Onchocercidae based on partitioned concatenated datasets of 12S rDNA, *coxI*, *rbp1*, *hsp70*, *myoHC*, 18S rDNA, and 28S rDNA sequences using Maximum Likelihood Inference. The total length of datasets is approximately 4950 bp. Sixty onchocercid specimens (representing 48 species) were analysed. *Filaria latala* and *Protospirura muricola* were used as outgroups. The topology was inferred using Maximum Likelihood (ML) inference. Nodes are associated with Bootstrap values based on 1000 replicates. The onchocercid subfamilies are indicated by colour: blue for Onchocercinae, dark green for Dirofilariinae, purple for Splendidofilariinae, pale green for Setariinae, yellow for Waltonellinae, orange for Icosiellinae and red for Oswaldofilariinae.
(TIF)

S1 Table. Supplementary information of filarial species used in phylogenetic analyses in this study. The species were classified by subfamilies indicated in the first column. The name of the species as well as the author(s) and date are given in the second column. Columns 3 to 8 present species information, i.e. family host, order host, localisation in the host, genera vector, family vector and geographical distribution. The references of morphological studies on which species identifications were based are listed in the last column. For the most part of the sampling some supplementary information were based on Anderson *et al.* [9].
(XLSX)

S2 Table. Ethical statement on vertebrate hosts and human samples. Vertebrate hosts are indicated in the first column and their collection place is given in the second column. The name of the collected filarial species is given in the third column. Columns 4 and 5 present ethical information on the host collection, i.e. the capture permit number and the ethical committee or legal entity for scientific procedures; Columns 6 to 9 describe non-scientific procedures in which filariae are recovered from vertebrates after hunting, necropsies in slaughter houses, or post-veterinary or post-medical procedures. The number of MNHN collection registration is listed in the last column.
(XLS)

S3 Table. Primers and PCR programs used in this study. Abbreviations: Step 1: denaturation; Step 2: annealing; Step 3: elongation; T: temperature (°C); D: duration (sec); N: number of cycles. * Indicates the primers designed for nested PCR. References: PS means newly design for the present study.
(XLS)

S4 Table. Accession number list of filarial species and outgroups used for phylogenetic analyses. The name of the species are given in the first column. The identification number of samples are given in the second column. Columns 3 to 9 present accession number for each analysed genes, i.e. *coxI* (cytochrome oxidase subunit I), 12S rRNA, 18S rRNA, *MyoHC* (myosin heavy chain), *hsp70* (70 kilodalton heat shock proteins), *rbp1* (RNA polymerase II large subunit) and 28S rRNA. Abbreviations: ext. indicated external samples.
(XLSX)
## S1 Appendix. List of author(s) and date associated with name of the vertebrate hosts or vector species.
(DOCX)

## Acknowledgments

The authors wish to thank the following collaborators for the material they provided: Dr. M. Casiraghi (DNA of *Aproctella* sp., R. andersoni, D. gracile, M. martini, F. candezei), Dr. S. Uni (DNA of O. eberhardi, O. d. japonica, O. skrjabini, L. caprini), Dr. S. Nikander (S. tundra), Dr. D. Otranto (C. bainae), Dr. M. Diagne (M. martini), Dr. S. Wanji (O. armillata, O. gutturosa, O. ochengi), Dr. B. Makepeace (S. labiatopapillosa), Dr. Alexis Ribas (B. (B.) jittapalapongi), Dr. J. Dupouy-Camet (*Loa loa*), Dr. S. Michalski from FR3 center (B. malayi, B. pahangi, A. viteae), S. Landi (B. Timori), Bayer Animal Health GmbH (*Dirofilaria immitis*).

## Author Contributions

Conceived and designed the experiments: CM OB LG EL. Performed the experiments: EL JB. Analyzed the data: EL LG. Contributed reagents/materials/analysis tools: EL CM KJ RG YK YM CR IC TBTS SdSL JMCC. Wrote the paper: EL CM LG.

## References

1. Bain O (2002) Evolutionary relationships among filarial nematodes. In: Rajan K, editor. World Class Parasites: Kluwer Acad. pp. 181.

2. Bouchery T, Lefoulon E, Karadjian G, Nieguitsila A, Martin C (2013) The symbiotic role of *Wolbachia* in Onchocercidae and its impact on filariasis. Clin Microbiol Infect 19: 131–140. doi: 10.1111/1469-0691.12069 PMID: 23398406

3. Hoerauf A, Pfarr K, Mand S, Debrah AY, Specht S (2011) Filariasis in Africa-treatment challenges and prospects. Clin Microbiol Infect 17: 977–985. doi: 10.1111/j.1469-0691.2011.03586.x PMID: 21722251

4. Taylor MJ, Hoerauf A, Bockarie M (2010) Lymphatic filariasis and onchocerciasis. Lancet 376: 1175–1185. doi: 10.1016/S0140-6736(10)60586-7 PMID: 20739055

5. Bain O, Mutafchiev Y, Junker K (2013) Order Spirurida. In: Schmidt-Rhaesa A, editor. Handbook of Zoology Vol 2 Nematoda. Berlin: De Gruyter. pp. 661–732.

6. Anderson RC, Bain O (1976) Keys to genera of the order Spirurida. Part 3. Diplotriaenoidea, Aproctoidea and Filarioidea. In: Anderson RC, Chabaud AG, Willmott S, editors. CIH Keys to the nematode parasites of vertebrates. England: CAB International. pp. 59–116.

7. Gibbons LM (2010) Keys to the nematode parasites of vertebrates. Supplementary volume. Wallingford: CAB International. i–viii, pp 416.

8. Bain O, Mutafchiev Y, Junker K, Guerrero R, Martin C, et al. (2015) Review of the genus *Mansonella* Faust, 1929 sensu lato (Nematoda: Onchocercidae), with descriptions of a new subgenus and a new subspecies. Zootaxa 3918: 151–193. doi: 10.11646/zootaxa.3918.2.1 PMID: 25781087

9. Anderson CR (2000) Nematode Parasites of Vertebrates: Their Development and Transmission 2nd Edition. Wallingford, Oxon (UK): C.A.B. International. pp 672.

10. Chabaud AG, Bain O (1994) The evolutionary expansion of the Spirurida. Int J Parasitol 24: 1179–1201. PMID: 7729976

11. Wehr EE (1935) A revised classification of the nematode superfamily Filarioidea. Proc Helminthol Soc Wash 2: 84–88.

12. Chabaud AG, Anderson CR (1959) Nouvel essai de classification des filaires (superfamille des Filarioidea) II. Ann Parasitol Hum Comp 34: 64–87.

13. Chabaud AG, Choquet MT (1953) Nouvel essai de classification des filaires (super-famille des Filarioidea). Ann Parasitol Hum Comp 28: 172–192. PMID: 13092687

14. Xie H, Bain O, Williams SA (1994) Molecular phylogenetic studies on filarial parasites based on 5S ribosomal spacer sequences. Parasite 1: 141–151. PMID: 9140481

15. Casiraghi M, Anderson TJ, Bandi C, Bazzocchi C, Genchi C (2001) A phylogenetic analysis of filarial nematodes: comparison with the phylogeny of *Wolbachia* endosymbionts. Parasitology 122 Pt 1: 93–103. PMID: 11197770
16. Casiraghi M, Bain O, Guerrero R, Martin C, Pocacqua V, et al. (2004) Mapping the presence of *Wolbachia pipientis* on the phylogeny of filarial nematodes: evidence for symbiont loss during evolution. Int J Parasitol 34: 191–203. PMID: 15037105

17. Ferri E, Barbuto M, Bain O, Galimberti A, Uni S, et al. (2009) Integrated taxonomy: traditional approach and DNA barcoding for the identification of filarioid worms and related parasites (Nematoda). Front Zool 6: 1. doi: 10.1186/1742-9994-6-1 PMID: 19128479

18. Ferri E, Bain O, Barbuto M, Martin C, Lo N, et al. (2011) New insights into the evolution of *Wolbachia* infections in filarial nematodes inferred from a large range of screened species. PLoS One 6: e20843. doi: 10.1371/journal.pone.0020843 PMID: 21731626

19. Lefoulon E, Gavotte L, Junker K, Barbuto M, Uni S, et al. (2012) A new type *F. Wolbachia* from *Splendidofilariinae* (Onchocercidae) supports the recent emergence of this supergroup. Int J Parasitol 42: 1025–1036. doi: 10.1016/j.ijpara.2012.09.004 PMID: 23041355

20. Blaxter ML, De Ley P, Garey JR, Liu LX, Scheldeman P, et al. (1998) A molecular evolutionary framework for the phylum Nematoda. Nature 392: 71–75. PMID: 9510248

21. Nadler SA, Carreno RA, Mejia-Madrid H, Ullberg J, Pagan C, et al. (2007) Molecular phylogeny of clade III nematodes reveals multiple origins of tissue parasitism. Parasitology 134: 1421–1442. PMID: 17506928

22. Park JK, Sultana T, Lee SH, Kang S, Kim HK, et al. (2011) Monophyly of clade III nematodes is not supported by phylogenetic analysis of complete mitochondrial genome sequences. BMC Genomics 12: 392. doi: 10.1186/1471-2164-12-392 PMID: 21813000

23. Krueger A, Fischer P, Morales-Hojas R (2007) Molecular phylogeny of the filaria genus *Onchocerca* with special emphasis on Afrotropical human and bovine parasites. Acta Trop 101: 1–14. PMID: 17174932

24. Xie H, Bain O, Williams SA (1994) Molecular phylogenetic studies on *Brugia* filariae using Hha I repeat sequences. Parasite 1: 255–260. PMID: 9140492

25. McNulty SN, Mullin AS, Vaughan JA, Tkach VV, Weil GJ, et al. (2012) Comparing the mitochondrial genomes of *Wolbachia*-dependent and independent filarial nematode species. BMC Genomics 13: 145. doi: 10.1186/1471-2164-13-145 PMID: 22530989

26. Pereira FB, Souza SL, Bain O (2010) *Oswaldofilaria chabaudi* n. sp. (Nematoda: Onchocercidae) from a South American tropidurid lizard (Squamata: Iguania) with an update on *Oswaldofilariinae*. Parasite 17: 307–318. PMID: 21275236

27. Bain O, Sulahian A (1974) Trois nouvelles filaires du genre *Oswaldofilaria* chez des Lézards sud-américains; essai de classification des Oswaldofilariinae. Bull Mus Natl Hist Nat 156: 827–841.

28. Bain O, Chabaud AG (1975) Développement chez des moustiques de trois filaires de lézards Sud-Americains du genre *Oswaldofilaria*. Ann Parasitol Hum Comp 50: 209–222. PMID: 240307

29. Desportes C (1941) Nouvelles recherches sur la morphologie et sur l'évolution d'Icosiella neglecta (Diesing) filaire commune de la grenouille verte. Ann Parasitol Hum Comp 18: 46–67.

30. Nikander S, Laaksonen S, Saari S, Oksanen A (2007) The morphology of the filaroid nematode *Setaria tundra*, the cause of peritonitis in reindeer *Rangifer tarandus*. J Helminthol 81: 49–55. PMID: 17381867

31. Bartlett CM, Greiner EC (1986) A revision of *Pelecitus Railliet* and Henry, 1910 (Filarioidea, Dirofilariinae) and evidence for the 'capture' by mammals of filarioids from birds. Bull Mus Natl Hist Nat Ser. 4, Section A: 47–99.

32. Sonin MD (1975) Filariata of animals and man and disease caused by them. Onchocercidae.; Skrjabin KI, editor. Moscow: Nauka Publisher.

33. Bain O, Petit G, Kozek WJ, Chabaud AG (1981) Sur les Filaires Splendidofilariinae du genre *Aproctella*. Ann Parasitol Hum Comp 56: 95–105. PMID: 7259001

34. Bartlett CM, Anderson RC (1980) Filarioid nematodes (Filarioidae: Onchocercidae) of corvus *Brachyrhynchos brachyrhynchos* brehm in southern Ontario, Canada, and a consideration of epizootiology of avian Filariasis. Syst Parasitol 2: 77–102.

35. Hering-Hagenbeck S, Boomker J, Petit G, Killick-Kendrick M, Bain O (2000) Description of Madathumugadia hiepei n. sp. (nematoda: Splendidofilariinae), a parasite of a South African gecko, and its development in laboratory bred Phlebotomus dubosqi (Diptera: Psychodidae). Syst Parasitol 47: 207–213. PMID: 11071156

36. Laaksonen S, Saari S, Nikander S, Oksanen A, Bain O (2010) Lymphatic dwelling filarioid nematodes in reindeer *Rangifer tarandus tarandus*, (cervidae) in Finland, identified as *Rumenfilaria andersoni* Lankester & Snider, 1982 (nematoda: Onchocercidae: Splendidofilariinae). Parasite 17: 23–31. PMID: 20387735

37. Kuzmina TA, Kuzmin YI, Tkach VV, Spraker TR, Lyons ET (2013) Ecological, morphological, and molecular studies of *Acanthocheilonema odendhali* (Nematoda: Filarioidea) in northern fur seals
(Callorhinus ursinus) on St. Paul Island, Alaska. Parasitol Res 112: 3091-3100. doi: 10.1007/s00436-013-3483-3 PMID: 23760875

38. Veciana M, Bain O, Morand S, Chaisiri K, Douangboupha B, et al. (2015) *Breinlia* (*Breinlia*) *jittapalapongoi* n. sp. (Nematoda: Filarioidea) from the Asian house rat *Rattus tanezumi* Temminck in Lao PDR. Syst Parasitol 90: 237-245. doi: 10.1007/s11230-014-9544-x PMID: 25693458

39. Partono F, Oemijati S, Dennis DT, Purnomo, Atmosoedjono S, et al. (1976) The Timor filaria on Flores and experimental transmission of the parasite. Trans R Soc Trop Med Hyg 70: 354-355. PMID: 12599

40. Otranto D, Varcasia A, Solinas C, Scala A, Brianti E, et al. (2013) Redescription of *Cercopithifilaria bai-nae* Almeida & Vicente, 1984 (Spirurida, Onchocercidae) from a dog in Sardinia, Italy. Parasit Vectors 6: 132. doi: 10.1186/1756-3305-6-132 PMID: 23642161

41. Lefoulon E, Kuzmin Y, Plantard O, Mutafchiev Y, Otranto D, et al. (2014) Redescription of *Cercopithifilaria rugosicauda* (Bohm & Supperer, 1953) (Spirurida: Onchocercidae) of roe deer, with an emended diagnosis of the genus *Cercopithifilaria* and a genetic characterisation. Parasitol Int 63: 808-816. doi: 10.1016/j.parint.2014.07.011 PMID: 25108130

42. Eberhard ML, Morales GA, Orihel TC (1976) *Cruorifilaria tuberocauda* gen. et sp. n. (Nematoda: Filar-ioidea) from the capybara, *Hydrochoerus hydrochaeris* in colombia. J Parasitol 62: 604-607. PMID: 957038

43. Eberhard ML, Lowrie RC Jr., Orihel TC (1979) Development of *Dipetalonema gracile* and *D. caudispina* to the infective stage in *Culicoides hollensis*. J Parasitol 65: 89-95. PMID: 571910

44. Bain O, Petit G, Rosales-Loesener L (1986) Filaires de Singes sud-américains. Bull Mus Natl Hist. Nat 4e ser., 8: 513-542.

45. Bain O, Diagne M, Muller R (1987) Une cinquième filaire du genre *Dipetalonema*, parasite de singes Sud-Americains. Ann Parasitol Hum Comp 62: 262-270. PMID: 3662331

46. Guerrero R, Martin C, Gardner SL, Bain O (2002) New and Known Species of *Litomosoides* (Nema-toda: Filarioidae): Important Adult and Larval Characters and Taxonomic Changes. Comp Parasitol 69: 177-195.

47. Uni S, Bain O, Agatsuma T, Katsumi A, Baba M, et al. (2006) New filarial nematode from Japanese serows (*Naemorhedus crispus*: Bovidae) close to parasites from elephants. Parasite 13: 193-200. PMID: 17007210

48. Uni S, Bain O, Takaoka H (2004) Affinities between *Cutifilaria* (Nematoda: Filarioidea), parasites of deer, and *Mansonella* as seen in a new onchocercid, *M*. (C.) *perforata* n. sp., from Japan. Parasite 11: 131-140. PMID: 15224573

49. Bain O, Petit G, Gueye A (1985) Transmission experimentale de *Monanema nilotica* El Bihari et coll., 1977, filaire a microfilaires dermiques parasite de Murides africains. Ann Parasitol Hum Comp 60: 83-89. PMID: 3985535

50. Fukuda M, Takaoka H, Uni S, Bain O (2008) Infective larvae of five *Onchocerca* species from experi-mentally infected *Simulium* species in an area of zoonotic onchocerciasis in Japan. Parasite 15: 111-119. PMID: 18642503

51. Uni S, Bain O, Takaoka H, Miyashita M, Suzuki Y (2001) *Onchocerca dewittei japonica* n. subsp., a common parasite from wild boar in Kyushu Island, Japan. Parasite 8: 215-222. PMID: 11584751

52. Bain O (1975) Redescription de cinq especes d'onchocerques. Ann Parasitol Hum Comp 50: 763-788. PMID: 1227381

53. Uni S, Bain O, Agatsuma T, Harada M, Torii H, et al. (2007) *Onchocerca eberhardi* n. sp. (Nematoda: Filarioidea) from sika deer in Japan; relationships between species parasitic in cervids and bovids in the Holarctic region. Parasite 14: 199-211. PMID: 17933297

54. Bain O, Schulz-Key H (1974) Les onchocerques du cerf europeen: redescription d'O. *flexuosa* (Wedl, 1856) et description d'O. *tubingensis* n.sp. et O. *tarsicola* n.sp. Tropenmed Parasitol 25: 437-449. PMID: 4446067

55. Bain O (1979) Transmission de l'Onchocerque bovine, *Onchocerca gutturosa*, par *Culicoides*. Ann Parasitol Hum Comp 54: 483-488. PMID: 533114

56. Yagi K, Bain O, Shoho C (1994) *Onchocerca suzukii* n. sp. and O. *skrjabini* (= O. *tarsicola*) from a relict bovid, *Capricornis crispus*, in Japan. Parasite 1: 349-356. PMID: 9140501

57. Voelker J, Garms R (1972) Morphological studies of unknown filarial larvae from *Simulium damnosum*, the vector of onchocerciasis, and from S. *kenyae* in Liberia with reference to their possible final hosts. Z Tropenmed Parasitol 23: 285-301. PMID: 4634262

58. Yates JA, Lowrie RC (1984) Development of *Yatesia hydrochoerus* (Nematoda: Filarioidae) to the infective stage in ixodid ticks. Proceeding of Helminthological Society of Washington 51: 187-190.
59. Bain O, Chabaud AG (1986) Atlas des larves infestantes de Filaires. Tropical Medicine and Parasitol- ogy 37: 237-340.

60. Katoh K, Toh H (2008) Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform 9: 286-298. doi: 10.1093/bib/bbn013 PMID: 18372315

61. Li W, Cowley A, Uludag M, Gur T, McWilliam H, et al. (2015) The EMBL-EBI bioinformatics web and programmatic tools framework. Nucleic Acids Res 43: W580–584. doi: 10.1093/nar/gkv279 PMID: 25845596

62. Posada D (2008) jModelTest: phylogenetic model averaging. Mol Biol Evol 25: 1253-1256. doi: 10.1093/molbev/msn083 PMID: 18397919

63. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572-1574. PMID: 12912839

64. Stamatakis A (2014) RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30: 1312-1313. doi: 10.1093/bioinformatics/btu033 PMID: 24451623

65. Bain O, Casiraghi M, Martin C, Uni S (2008) The Nematode Filarioidea: critical analysis linking molecular and traditional approaches. Parasite 15: 342-348. PMID: 18814705

66. Bain O (1981) Filariids and their evolution. Parasitology 82: 161-174.

67. Bain O (1970) Morphologie larvaire de *Setaria labiatopapillosa* (Nematoda, Filarioidea) chez *Aedes aegypti*. Ann Parasitol Hum Comp 45: 431-439. PMID: 5511124

68. Bain O, Wanji S, Enyong P, Petit G, Noireau F, et al. (1998) New features on the moults and morphogenesis of the human filaria *Loa loa* by using rodent hosts consequences. Parasite 5: 37-46. PMID: 9754295

69. Chabaud AG, Bain O (1976) La lignée *Dipetalonema*. Nouvel essai de classification. Ann Parasitol Hum Comp 51: 365-397. PMID: 988773

70. Bain O, Baker M, Chabaud AG (1982) Nouvelles données sur la lignée *Dipetalonema* (Filarioidea, Nematoda). Ann Parasitol Hum Comp 57: 593-620. PMID: 6891995

71. Junker K, Medger K, Lutermann H, Bain O (2012) *Monanema joopi* n. sp. (Nematoda, Onchocercidae) from *Acomys* (*Acomys* spinosissimus Peters, 1852 (Muridae) in South Africa, with comments on the filarial genus. Parasite 19: 331–340. PMID: 23193517

72. Buckley JJ, Edeson JF (1956) On the adult morphology of *Wuchereria* sp. (malayi?) from a monkey (*Macaca irus*) and from cats in Malaya, and on *Wuchereria pahangi* n.sp. from a dog and a cat. J Helminthol 30: 1-20. PMID: 13319678

73. Eberhard ML, Orihel TC (1984) The genus *Mansonella* (syn. *Tetrapetalonema*): a new classification. Ann Parasitol Hum Comp 59: 483-496. PMID: 6508144

74. Bain O (1974) Description de nouvelles Filaires Oswaldofilariinae de Lézards sud-américains; hypothèse sur l'évolution des Filaires de Reptiles. Bull Mus Natl Hist Nat 3º ser n°208, 138: 169-200.

75. Chilton NB, Jabbar A, Huby-Chilton F, Jex A, Gasser RB, et al. (2012) Genetic variation within the *Hypodontus macropi* (Nematoda: Strongyloidea) complex from macropodid marsupial hosts in Australia. Electrophoresis 33: 3544-3554. doi: 10.1002/elp.s.201200364 PMID: 23161349

76. Susoy V, Herrmann M (2014) Preferential host switching and codivergence shaped radiation of bark beetle symbionts, nematodes of *Micoletzkya* (Nematoda: Diplogastridae). J Evol Biol 27: 889-898. doi: 10.1111/jeb.12367 PMID: 24731122

77. Quentin JC (1971) Sur les modalités d'évolution chez quelques lignées d'Helminthes de Rongeurs Murioidea. Cahiers ORSTOM (Série entomologie médicale et parasitologie) IX: 103-176.

78. Masi S, Chauffour S, Bain O, Todd A, Guillot J, Krief S (2012) Seasonal Effects on Great Ape Health: A Case Study of Wild Chimpanzees and Western Gorillas. PLoS ONE 7(12): e49805. doi: 10.1371/journal.pone.0049805 PMID: 23227152