# Emergent and Neglected Equine Filariosis in Egypt:
## Species Diversity and Host Immune Response

Faten a M Abo-Aziza, Seham H M Hendawy, Hend H a M Abdullah, Amira

El Namaky, Younes Laidoudi, Oleg Mediannikov

► To cite this version:

Faten a M Abo-Aziza, Seham H M Hendawy, Hend H a M Abdullah, Amira El Namaky, Younes Laidoudi, et al.. Emergent and Neglected Equine Filariosis in Egypt: Species Diversity and Host Immune Response. Pathogens, 2022, 11 (9), pp.979. 10.3390/pathogens11090979. hal-03948730

HAL Id: hal-03948730

https://amu.hal.science/hal-03948730

Submitted on 20 Jan 2023

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 Attribution 4.0 International License
Article

# Emergent and Neglected Equine Filariosis in Egypt: Species Diversity and Host Immune Response

Faten A. M. Abo-Aziza<fnref n="1" />, Seham H. M. Hendawy <fnref n="1" />,<fnref n="2" />, Hend H. A. M. Abdullah <fnref n="1" />,<fnref n="3" /> ID, Amira El Namaky <fnref n="1" />, Younes Laidoudi <fnref n="3" />,<fnref n="4" /> ID and Oleg Mediannikov <fnref n="3" />,* ID

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

check for
updates

**Citation:** Abo-Aziza, F.A.M.; Hendawy, S.H.M.; Abdullah, H.H.A.M.; El Namaky, A.; Laidoudi, Y.; Mediannikov, O. Emergent and Neglected Equine Filariosis in Egypt: Species Diversity and Host Immune Response. *Pathogens* **2022**, *11*, 979.
https://doi.org/10.3390/pathogens11090979

Academic Editors:
Francesca Mancianti and Valentina
Virginia Ebani

Received: 24 July 2022

Accepted: 17 August 2022

Published: 27 August 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** Equine filariosis (EF) is a neglected vector-borne disease caused by nematode species belonging to the Onchocercidae and Setariidae families. Aside from their zoonotic potential, some species are responsible for serious health problems in equids worldwide, leading to significant economic difficulties. Here, we molecularly investigated equine blood samples (320 horses and 109 donkeys from Egypt) and four adult worms isolated from the peritoneal cavity of 5 out of the 94 slaughtered donkeys. In addition, quantitative enzyme-linked immunoassays (ELISAs) targeting circulating cytokines were used to identify whether the immunological profile of the infected animals is a Th1 (i.e., INF-gamma as indicator) or Th2 (i.e., IL-5 and IL-10 as indicators) response type. Overall, 13.8% and 0.3% of the donkeys and horses, respectively, were scored as positive for filaroid DNA. The 18S phylogeny revealed the occurrence of three different filaroid species, identified here as *Mansonella* (*Tetrapetalonema*) sp., *Setaria digitata* and *Dirofilaria repens*. Th1 (INF-gamma and IL-5) and Th2 (IL-10) immune response types were identified in equines infected with *S. digitata* and *Mansonella* (T.) sp., respectively. These results provide new data on the species diversity of EF in Egypt and extend knowledge of the downregulation of the protective immune response by the potentially zoonotic *Mansonella* (T) sp. There is an urgent need to implement control measures to preserve equine health and limit the propagation of these vector-borne filaroids in Egypt.

**Keywords:** equids; zoonotic filariosis; cytokines; phylogeny; *Mansonella* sp.; *Setaria digitata*; *Dirofilaria repens*

## 1. Introduction

Equine filariosis (EF) is a neglected vector-borne disease caused by heteroxenous parasitic nematodes belonging to the Onchocercidae and Setariidae families. Adult nematode parasitise tissues and the body cavities of equids produce skin- or blood-dwelling microfilariae. The larvae are ingested by blood-feeding arthropods and undergo two-stage development followed by the third stage (L3), where they are transmitted to a new receptive host [1–5]. The most known filarioid species infecting equids belong to the *Onchocerca*, *Dirofilaria* and *Setaria* genera. Infection by these parasites leads to a myriad of problems ranging from fibrinous myocarditis and/or peritonitis to tendinitis when the worms infect the muscles and ligaments [6]. In all cases, the infection is often responsible for considerably limited movement and lameness [7–9]. In addition, the unusual localisations of the worms in the eye and the central nervous system may occur and may lead to blindness and neurological disorders [10,11]. Moreover, the deep localisation of the adult worms inside the body cavities and tissues and the heterogeneous localisation of their microfilariae require invasive methods, such as skin biopsies or the surgical removal of the parasite,

<footnote n="1">¹ Department of Parasitology and Animal Diseases, Veterinary Research Institute, National Research Centre, Dokki, Giza 12622, Egypt</footnote>

<footnote n="2">² Tick and Tick-Borne Diseases Research Unit, Veterinary Research Institute, National Research Centre, Dokki, Giza 12622, Egypt</footnote>

<footnote n="3">³ Aix Marseille Université, IRD, AP-HM, MEPHI, IHU-Méditerranée Infection, 13005 Marseille, France</footnote>

<footnote n="4">⁴ PADESCA Laboratory, Veterinary Science Institute, University of Constantine 1, El-Khroub 25100, Algeria</footnote>

<footnote>* Correspondence: olegusss1@gmail.com; Tel.: +33-04-1373-24-01</footnote>
which represents a major limitation for the routine diagnosis of such infection. Therefore, in the absence of specific serological tests, the use of molecular identification remains the most effective and accurate method for identifying these parasites from blood or tissue fragments of infected animals [12–14]. This significant lack of knowledge of their epidemiology and biology makes their control and management difficult.

Host–parasite interaction relies on the expressed cytokine profile which is modulated by the parasite itself [15–17]. The predominant immunological response during filarial infection is antigen-specific Th2 stimulation, leading to an expansion of IL-10, along with the inhibition of Th1 response, Babu and Nutman [18]. This explains the parasite longevity within the infected hosts [19,20]. It is admitted that filarial parasites produce and secrete a range of immunomodulatory molecules, such as the ES-62, that interfere directly with the natural immune response by suppressing the inflammatory phenomenon and countering the pathology associated with the disease [15]. However, host–parasite interaction differs between species, resulting in the specific modulation of the immune response according to the filaroid species. For example, the production of the cystatins by some filarial parasites such as *Onchocerca volvulus* and *Acanthocheilonema viteae* leads to the inhibition of the cysteine protease group (i.e., legumains and cathepsin L and S) and induces the production of anti-inflammatory interleukin (i.e., IL-10) [21,22]. The precise immune mechanisms that govern *S. equina* infection in donkeys have been partially investigated and depend on the TNF-α and IL-4 as Th1/Th2 balance with regard to the existence of the adult stages [6]. For example, *S. equina* asymptomatic infections are often characterised by the down- and upregulation of TNF-α, and IL-4, respectively [6], while chronic infections are associated with an increase in IL-10 expression, leading to the down- and upregulation of IFN-γ (Th1) and IL-5 (Th2), respectively [15,23].

In Egypt, equids play a significant role in numerous sectors, including agriculture, the police services, tourism and the pharmaceutical and transport sectors, as well as in recreational activities (i.e., sport, gaming and entertainment) [24,25]. Thus, understanding the impact of infectious pathogens on equines should not be overlooked as a vital part of Egypt's culture and economy. Studies on neglected filaroid parasites (i.e., *Setaria equina*, *Onchocerca cervicalis*, *Onchocerca reticulata* and *Parafilaria multipapillosa*) remain scant and outdated from many Egyptian provinces [7,12,26–32]. Therefore, this study aimed to investigate the species diversity of equine filariod and their impact on the host immune response of both horses and donkeys from Egypt.

**2. Results**

**2.1. Prevalence and Phylogeny of Filaroids**

Overall, 22 (5.1%, 95% CI) (1 horse and 21 donkeys) out of 429 equids were investigated and tested positive using pan-filarial qPCR. Filaroid DNA was successfully amplified from 16 out of the 22 qPCR-positive blood samples, and all (n = 4) retrieved worms using the 18S primers. Infection rates of 3.7% (0.3% in horse (1/320) and 13.8% in donkeys (15/109)) were recorded. Likewise, all DNA amplicons were sequenced and grouped into three sequence groups according to BLAST analysis (Table 1): group i) amplified from five donkeys and one worm retrieved from a donkey's peritoneal cavity and showed 100% (739/739) and 99% (736/737) identity with *S. digitata* from the UK (GenBank DQ094175) and Egypt (GenBank MN728217), respectively. Group ii) amplified from a microfilaremic donkey and a worm from a donkey's peritoneal cavity and were 100% (739/739) identical with those of *D. repens* previously detected in dog from Egypt (GenBank MN728215). Group iii) amplified from 10 blood samples (one horse and nine donkeys) and one worm retrieved from a donkey's peritoneal cavity and were 100% (1057/1057) identical to *Mansonella* sp. (genotype OM- 2015, GenBank MT786947-49), which is detected in donkeys from Algeria (Table 1). Accordingly, the maximum likelihood phylogeny confirmed the BLAST results for the three filaroid species (*Mansonella* sp., *S. digitata* and *D. repens*) with high bootstrap values (Figure 1). In addition, the ML phylogeny resolved the species identification of the *Mansonella* sp. herein isolated at the subgenus *Tetrapetalonema* (Figure 1).
<table><caption>**Table 1.** Results of molecular identification of equine filaroids from the analysed samples (i.e., blood or filaroid worms). Species identification, prevalence and sequence accession number are shown.</caption><thead><tr><th>Animal Hosts (no.)</th><th>No. of Filaroid-Positive Animals (%)</th><th>Species Name (Acc. No., Size bps)</th><th>100% Identical Sequences from GenBank</th></tr></thead><tbody><tr><td>Horses (320)</td><td>1 (0.3)</td><td rowspan="2">Mansonella (T.) sp. (MN728216)</td><td rowspan="2">Mansonella sp. (MT786947-49)</td></tr><tr><td>Total</td><td>1 (0.3)</td></tr><tr><td rowspan="3">Donkeys (109)</td><td>9 (8.3)</td><td>Mansonella (T.) sp. (MN728182, MN728181 *)</td><td>Mansonella sp. (MT786947-49)</td></tr><tr><td>5 (4.6)</td><td>Setaria digitata (MN728184, MN728183 **)</td><td>S. digitata (DQ094175)</td></tr><tr><td>1 (0.9)</td><td>Dirofilaria repens (MN728180, MN728180 *)</td><td>D. repens (MN728215)</td></tr><tr><td>Total</td><td>15 (13.8)</td><td></td><td></td></tr></tbody></table>

## 2.2. Immunological Studies

The slopes of regression (b) among the graded log-doses and their corresponding binding were 0.467, 0.529 and 0.402 for INF-gamma, IL-5 and IL-10, respectively. *Setaria digitata* infection induced a higher inflammatory response, as INF-gamma was significantly higher in *S. digitata*-infected equines (*p* < 0.05) compared to controls. Both *Mansonella* (T) sp. and *S. digitata* infection resulted in an increase in Th2 cytokines, but the elevation in *S. digitata*-infected animals was more pronounced. It was found that the IL-5 level in *S. digitata* was significantly higher (*p* < 0.01) than in the control group. Similarly, equines infected with *Mansonella* (T) sp. showed a significant elevation in the IL-5 level relative to the control group (*p* < 0.05). It was observed that equines infected with *Mansonella* (T) sp. showed a significant elevation in IL-10 level (*p* < 0.01) compared to the control group (Table 2). Moreover, comparing the effect of the infection with the two filarial species in terms of cytokine levels, it was found that INF-gamma and IL-5 levels were significantly higher in equines infected with *S. digitata* than those of equines infected with *Mansonella* (T) sp. (*p* < 0.05). However, *Mansonella* (T) sp. infections resulted in a significant elevation in IL-10 level (*p* < 0.01) compared to *S. digitata* infections. In addition, Th1/Th2 in *S. digitata* was 1.99, while in *Mansonella* (T) sp., Th1/Th2 was 1.71 (Table 2). In the control group, a moderate positive correlation was observed between INF-gamma and IL-5, though a weak correlation between INF-gamma and IL-10 serum levels was recorded. A strong positive correlation was recorded between INF-gamma, IL-5 and IL-10 serum levels in animals infected with *S. digitata*. In contrast, a moderate positive correlation between INF-gamma and IL-5 serum levels and a negative moderate correlation between INF-gamma and IL-10 serum levels were observed in *Mansonella* sp. (T)-infected groups (Figure 2).

<table><caption>**Table 2.** Serum cytokines profile of *Mansonella* (T) sp.- and *S. digitata*-infected equids.</caption><thead><tr><th></th><th>INF-gamma<br>(pg/mL)</th><th>IL-5 (pg/mL)</th><th>Th1/Th2<br>(INF-<br>gamma/IL-5)</th><th>IL-10 (pg/mL)</th></tr></thead><tbody><tr><td>Control</td><td>23.634 ± 1.66</td><td>6.160 ± 0.34</td><td>3.84</td><td>8.726 ± 2.11</td></tr><tr><td>Mansonella (T) sp.</td><td>21.992 ± 3.18</td><td>11.054 ± 2.61 <sup>a</sup></td><td>1.99</td><td>17.868 ± 1.02 <sup>bB</sup></td></tr><tr><td>S. digitata</td><td>28.006 ± 1.09 <sup>aA</sup></td><td>16.362 ± 1.35 <sup>bA</sup></td><td>1.71</td><td>9.254 ± 2.99</td></tr></tbody></table>

<footnote>* and ** indicate accession numbers of sequence of the filaroid worms from one infected donkey or two infected donkeys, respectively.</footnote>

<footnote>a and b indicate significant difference in cytokines concentration at *p*-value of 0.01 and 0.05, respectively, between infected and control animals. A and B indicate significant difference in cytokine profile at *p*-value of 0.01 and 0.05, respectively, between animals of the same host species infected with different filaroid species (*Mansonella* (T) sp. and *S. digitata*).</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** The 18S-rRNA-based phylogeny of filaroid species identified in this study and the representative members of the five Onchocercidae Clades (ONC1-5) available in the GenBank database. The tree corresponds to the IQTREE inferred from 39 partial (713 bps) sequences with 9.7% (i.e., 52 and 71 of parsimony informative and distinct sites, respectively) of informative sites using the K2P (+R2) model under 1000 bootstrap replicates and the ML method. GenBank accession number, species name, host, and geographical origin are indicated at the tip of each branch when available. Blue bolded label indicates sequence types are amplified in the present study. Bootstraps values higher than 50% are printed at the branch nodes.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Correlation between INF-gamma serum level and IL-5 and IL-10 levels. The level of INF-gamma correlated with the levels of IL-5 and IL-10 in control (a) and infected animals (b,c). In the control group, a moderate positive correlation was observed between INF-gamma and IL-5; however, a weak correlation between INF-gamma and IL-10 serum levels was recorded. A strong positive correlation was recorded between INF-gamma and IL-5 and IL-10 serum levels in animals infected with *S. digitata*. In contrast, a moderate positive correlation between INF-gamma and IL-5 serum levels and a negative moderate correlation between INF-gamma and IL-10 serum levels were observed in *Mansonella* sp. (T)-infected groups.</figcaption></figure>
### 3. Discussion

Equine filariosis is a neglected vector-borne disease causing a reduction in the working capacity of equids. Consequently, health threat and economic losses are the main observed problems in endemic areas [32,33]. Therefore, their accurate detection coupled with a better understanding of cytokine profiles and the immunological signalling associated with the disease should not be overlooked for disease management and control [6,34]. The present study provides new data on the species diversity and disease-associated cytokines of EF from Egypt.

In the current study, the overall prevalence of EF was 3.7%, with a higher prevalence among donkeys (13.8%) than among horses (0.3%). Based on our designed 18S rRNA primers, three distinct filaroid species were detected in equines: *Mansonella* sp., *S. digitata* and *D. repens*.

The infection rate of *Mansonella* was 0.3% (1/320) in horses from Cairo and 8.3% (9/109) in donkeys from Beni Suef and Al-Faiyum. Phylogenetically, this species clustered together with *Mansonella* sp. OM-2015 genotype and belongs to the *Tetrapetalonema* subgenus within *Mansonella* (*T.*) *atelensis* (GenBank: MT336173). According to GenBank entries, the *Mansonella* sp. OM-2015 genotype was identified either in donkeys from Algeria using 18S, ITS1, *cox1* and 5S sequencing (genotype OM-2015, MT786947-49; MT786950-51, 55–58; MT791217 and MT795713, 15–17, respectively), and from Senegal using ITS1 sequencing (genotype OM-2015, MT786952), in horses from Senegal using ITS1 and 5S sequencing (genotype OM-2015, MT786953-54 and MT795714, respectively) or in biting midges *Culicoid esenderleini* from Senegal using generic ITS1 sequencing (genotype OM-2015, KR080175). Despite the scant data, this species seems to be widely distributed and was only identified in equids (i.e., donkeys and horses) and in biting midges *C. enderleini*, suggesting the reservoir and vector role of these hosts, respectively (Figure 3).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Geographical mapping of the available molecular data (i.e., 18S, ITS1, *cox1* and 5S) on *Mansonella* (*T.*) sp., genotype OM-2015 based on host source. The map was generated using Microreact server (available at: https://microreact.org accessed on 25 May 2022). Maps © Mapbox (www.mapbox.com/about/maps).</figcaption></figure>

This phylogenetic analysis resolved the identification of the *Mansonella* sp. genotype OM-2015 at the subgenus level (*Tetrapetalonema*). However, in the absence of morphological data on this species isolated in this study and given the absence of DNA sequences from the other morphologically described *Mansonella* (*T*) species, the possibility of the *Mansonella* (*T*) sp. described in this study being a new species cannot yet be ruled out. The subgenus *Tetrapetalonema* encompass 14 valid species all isolated from South American monkeys, with the exception of *Mansonella* (*T*) zakii [35,36] (syns. *Parlitomosazakii* [35]; *Tetrapetalonema*
zakii [35,37]; Dipetalonema zakii [35]) in Leontopithecus (= Leontocebus) rosalia (Linnaeus), which was isolated in Egypt (in captivity, originating from Brazil) [35] and is considered to be a species studied by Eberhard and Orihel [36]. More studies are needed to confirm or describe this potential new species of the *Mansonella* group and determine its importance in pathology.

We detected *S. digitata* and its microfilariae in donkeys derived from Al-Faiyum province, with a total prevalence rate of 4.6% (5/109). Setiosis caused by different *Setaria* spp. is a filarial disease affecting bovines and can accidentally be transmitted into unusual hosts such as horses, donkeys, sheep, goats, and camels [38–41]. *Setaria digitata* has been detected in horses in Korea [42,43], Malaysia [8], India [13], Iran [5], and China [9]. *S. equina* has been detected in donkeys in Egypt [7,12,32]. The identification of adult *S. digitata* in donkeys and microfilariae in the blood indicate that donkeys are competent hosts of *S. digitata* in Egypt. The relatively high prevalence of *S. digitata* microfilaraemia in donkeys in Egypt (4.6%) may lead to a significant impact on animal health.

We detected adult *D. repens* and its microfilariae in one donkey from Al-Faiyum province, with a total prevalence rate of 0.9% (1/109). To the best of our knowledge, *D. repens* has never been reported in donkeys anywhere in the world. *D. repens* affects dogs and other carnivores, causing subcutaneous dirofilariasis [44,45]. Other mammals, including humans, may also be affected [46], but they very rarely become the final hosts or develop microfilaraemia [47]. A generalisation stating that donkeys play a significant role as a reservoir of subcutaneous dirofilariasis cannot be made from the single animal we describe here; however, it is evidence that this animal may also be a competent host of *D. repens*.

An assay of cytokine responses (INF-gamma as Th1 and IL-5 and IL-10 as Th2 indicators) to *Mansonella* (*T*) sp. and *S. digitata* infection was performed to elucidate immunological differences, motivating an immune response. In this study, it was noted that equines infected with *Mansonella* (*T*) sp. showed a significant elevation in IL-10 levels compared to the control group. Moreover, comparing the effect of the infection with the two filarioid species on cytokine levels, it was found that INF-gamma and IL-5 levels were significantly higher in equines infected with *S. digitata* than those infected with *Mansonella* (*T*) sp. It is well known that polymorphous cells are essential for killing adult worms through the evoked inflammatory nodule formation and production of IL-5 in response to helminths [48]. IL-5 is also essential for containing parasitaemia [49]. In IL-5-deficient mice, microfilaraemia was greatly enhanced and pronounced in addition to increased adult worm burden [50,51]. In contrast, many studies have shown that Th1 cytokines (INF-gamma) can be involved in the protection against helminths through proper neutrophil migration [52]. However, *Mansonella* (*T*) sp. infection was correlated to a significant elevation of IL-10 level compared to *S. digitata* infection. In addition, Th1/Th2 in *S. digitata* was 1.99, but in *Mansonella* (*T*) sp., Th1/Th2 was 1.71. In the control group, a moderate positive correlation was observed between INF-gamma and IL-5, while a weak correlation between INF-gamma and IL-10 serum levels was recorded. A strong positive correlation was recorded between INF-gamma and IL-5 and IL-10 serum levels in animals infected with *S. digitata*. In contrast, a moderate positive correlation between INF-gamma and IL-5 serum levels and a negative moderate correlation between INF-gamma and IL-10 serum levels were observed in *Mansonella* (*T*) sp. infected groups. Helminths are able to downregulate host immune answers in order to facilitate their long-term existence. This ability to modulate the host immunological situation is a key part of the evolutionary achievement of helminths and could be attributed to immunomodulatory mediators such as ES-62, cystatins, legumains and cathepsin L and S, which are secreted by the nematode parasite as part of their immune-evasion strategy. They have an immunomodulatory effect to enhance the production of anti-inflammatory IL10 [17,22]. Such a strategy was more pronounced in *Mansonella* (*T*) sp. than in *S. digitata*. Inflammatory procedures are blocked and pathogenic damage to the host is reduced or even subclinical. This ultimately lets the parasite continue within the host for a long time [53]. Mechanisms that downregulate the immune responses include the induction of regulatory T cells [54]. T helper cells constitute
common T lymphocyte activation followed by differentiation into Th1 and Th2 phenotypes associated with the progress of type-2 cytokines and the impairment of type-1 cytokine production [55]. This shifting plays a key role in regulating the balance between infection and disease. In the case of Th1 and Th2 phenotypes, they are predominantly related to susceptibility and protection, respectively [56]. Subash and Nutman [57] reported that early filarial infection was accompanied by the elevation of Th1 over Th2 cytokines, which is essential to understanding the pathogenesis of infection and the host-parasite relationships. This response could be the beginning of acute filariasis and the formation of host resistance to the helminth infection [57]. Th1 and Th2 cytokines orchestrate different immune pathways to fight Strongylus; Th1 cytokines coordinate cellular immune responses, and Th2 cytokines coordinate humoral immune responses [58]. Helminths infections are usually associated with polarised Th2/Th1 immune responses [16,57]. In this study, Th1/Th2 was 1.99 in *S. digitata* and 1.71 in *Mansonella* (T) sp. These results indicated that both infections directed the effort of cytokines to the humoral immune response, as reflected by the elevation of IL-5. Conversely, this effort was reflected in cellular immune response, as indicated by the measured balance.

**4. Materials and Methods**

**4.1. Study Design and Animal Sampling**

Between 2016 and 2018, a cross-sectional study was conducted on equids from three Egyptian provinces (Cairo, Beni-Suef and Al-Faiyum; Figure 4). Using the convenience sampling strategy, a total of 429 equids, including 320 live horses and 109 donkeys (94 donkeys slaughtered at the zoo and 15 live donkeys), were randomly assigned to the study (Table 3). EDTA-blood and sera samples were, respectively, collected either from the jugular vein or the carotid arteries/jugular vein of live equids and the donkeys' corpses immediately after slaughter. The collected blood and serum samples were stored at -20 °C until analysis. Adult filaroids were searched for in all donkey corpses, and when found were stored in 70% alcohol until further investigation.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Map of Egypt shows the locations of provinces where *Equidae* blood samples were collected. https://en.wikipedia.org/wiki/Governorates_of_Egypt (accessed on 25 May 2022) and the picture has CC BY-SA 3.0.</figcaption></figure>
<table><caption>**Table 3.** The data of collected samples.</caption><thead><tr><th>Provinces</th><th>Lat./Long.</th><th>Hosts</th><th>Locations</th><th>Total</th></tr></thead><tbody><tr><td rowspan="3">Cairo</td><td rowspan="3">30°03'45.47" N,<br>31°14'58.81" E</td><td rowspan="3">Horses</td><td>Police Academy (El-Abbasia)</td><td>94</td></tr><tr><td>Police Academy (El-Tagamoa)</td><td>70</td></tr><tr><td>Police Academy (El-Pasateen)</td><td>147</td></tr><tr><td rowspan="2">Beni-Suef</td><td rowspan="2">29°03'60.00" N,<br>31°04'60.00" E</td><td>Horses</td><td>households</td><td>9</td></tr><tr><td>Donkeys</td><td>households</td><td>22</td></tr><tr><td>Al-Faiyum</td><td>29°18'35.82" N,<br>30°50'30.48" E</td><td>Donkeys</td><td>households</td><td>87</td></tr></tbody></table>

## 4.2. Molecular Studies

### 4.2.1. Polymerase Chain Reactions and Sequencing

First, all EDTA-blood samples (200 µL) and a piece of tissue (~25 mg) from each retrieved worm were subjected individually to mechanical and enzymatic lysis steps prior to DNA extraction as previously described [4]. DNA extraction was performed by the EZ1 biorobot Qiagen (Hilden, Germany), using the EZ1 DNA Tissue Kit Qiagen (Hilden, Germany) according to the manufacturer's instructions. DNA was eluted in 200 µL and stored at -20 °C until molecular analysis.

All DNA samples were molecularly screened for filaroid DNA using the pan-filarial real time qPCR (Pan-Fil 28S qPCR), targeting the 28S rRNA gene as described by Laidoudi and colleagues [45]. All samples which were filaroid-positive by the Pan-Fil 28S qPCR were subjected to PCR amplification using the newly designed primers (Fwd-18S-Nem.58: AATGGTGAAACCGCGAAC and Rwd-18S-Nem.998: AACACCGCTTGTCCCTCTAA). Briefly, primers were designed according to PCR design protocol described by Laidoudi and colleagues [45] and targeting a partial (1325-35 bps) small subunit (SSU) 18S rRNA gene of nematodes. PCR reactions were carried out in a total volume of 50 µL, consisting of 1 µL of each primer (50 nm), 25 µL of AmpliTaq Gold® 360 Master Mix (Thermo Fisher Scientific, Applied Biosystems, Foster City, CA, USA), 18 µL of DNAse-RNAse free water (Eurogentec, Liège, Belgium) and 5 µL of DNA template. PCR reactions were performed using a thermocycler (Applied Biosystems, Paris, France) with initial denaturation at 95 °C for 15 min, followed by 40 cycles including denaturation at 95 °C for 1 min, annealing at 55 °C for 30 s, and elongation at 72 °C for 1.5 min, followed by a final extension at 72 °C for 10 min. Positive and negative controls were included in each amplification. PCR amplification was confirmed in 1.5% agarose electrophoresis.

Finally, NucleoFast 96 PCR plates (Macherey Nagel, EURL, Hoerdt, France) were used for the purification of PCR products, in accordance with the manufacturer's recommendations. The purified PCR products were sequenced using the Big Dye Terminator Cycle Sequencing Kit (3130X1 Genetic Analyzer, ABI PRISM) within the ABI automated sequencer (Applied Biosystems). The obtained sequences were assembled and edited using ChromasPro software (ChromasPro 1.7, Technelysium Pty Ltd., Tewantin, Australia). The corrected sequences were subjected to preliminary analysis using the BLAST server [59].

### 4.2.2. Phylogenetic Analyses

Molecular sequences obtained via 18S PCR sequencing from worms and filaroid-positive samples were aligned against the representative members of the five Onchocercidae clades retrieved from the GenBank database using MAFFT v7.490 [60], with adjustment to the direction of the first sequence. Bioedit software (Bioedit version 7.2.5) was used to manually refine the multisequence alignment [61]. The maximum likelihood phylogeny was performed using IQTREE (IQ-TREE multicore version 1.6.12 for Mac OS X 64-bit) software [62], under 1000 Ultra-Fast bootstrap replications [63]. The K2P + R2 model was selected by ModelFinder before computing the tree [64]. DNA sequences of *Ascaris* sp. (GenBank accession number: JN256985), *Toxascaris leonine* (GenBank accession number: JN256984), *Abbreviata caucasica* (GenBank accession number: MN956825) and *Gongylonema*
*nepalensis* (GenBank accession number: LC278392) were used as out-groups to root the tree.
Finally, the tree was annotated within iTOL v5 software [65].

4.3. *Immunological Study*

Commercially available sandwich-type enzyme immunoassay kits using biotin-streptavidin chemistry (Genorise Scientific, INC., Glen Mills, PA, USA) were used to quantify the level of serum INF-gamma, IL-5 and IL-10 according to the manufacturer's instructions. Optical density was measured at a wavelength of 450 nm. To compare the amount of circulating interleukins, the slope of regression "b" of the log dose response curves (n = 5) was obtained from a respective standard preparation, calculated according to the method described by Spiegel [66]. In addition, the correlation of the level of circulating INF-gamma to IL-5 and IL-10 in control and filaroid-positive equids was statistically performed using the Minitap programme version 21. Finally, data from each cytokine concentration were presented as mean ± and standard errors (SE). Descriptive statistics and simple one-way analysis of variance (ANOVA) were performed using the SPSS programme version 19 (IBM Corp., Armonk, NY, USA) to characterise the immunological profile associated with each infection.

**5. Conclusions**

Our study provided an epidemiological picture of species diversity and immune
response modulation in equids that were naturally infected with filaroid parasites (*Mans-
onella, Setaria* and *Dirofilaria*) from Egypt. In addition to *S. digitata*, which is a well-known
filaroid infecting equids, the presence of both adult and larval forms of *Mansonella* (*T*)
sp. and *D. repens* extends knowledge on the host suitability and reservoir role of equids
for these zoonotic/potentially zoonotic species. In contrast to *S. digitata*, *Mansonella* (*T*)
sp. seems to downregulate the protective immune response. However, the pathogeny
of EF in general and equine mansonellosis caused by *Mansonella* (*T*) sp. remains unclear,
and further studies on the biology of these equine vector-borne nematodes are needed.
Likewise, further studies are needed to assist in the quest to identify the elusive vector of
the *Mansonella* (*T*) sp. as well as to promote its full morphological description.

**Author Contributions:** Conceptualisation, F.A.M.A.-A., S.H.M.H., H.H.A.M.A., Y.L. and O.M.; methodology, H.H.A.M.A., Y.L., S.H.M.H., F.A.M.A.-A. and A.E.N.; software, F.A.M.A.-A., H.H.A.M.A. and Y.L.; validation, F.A.M.A.-A., S.H.M.H. and H.H.A.M.A.; formal analysis and investigation, H.H.A.M.A., F.A.M.A.-A. and S.H.M.H.; resources, H.H.A.M.A., F.A.M.A.-A., S.H.M.H. and A.E.N.; data curation, H.H.A.M.A. and F.A.M.A.-A.; writing—original draft preparation, H.H.A.M.A., F.A.M.A.-A. and S.H.M.H.; writing—review and editing, Y.L. and O.M.; visualisation, F.A.M.A.-A., S.H.M.H., H.H.A.M.A., Y.L. and O.M.; supervision: O.M. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Science and Technology Development Fund (STDF), InstitutFrancaisd'Egypte (IFE) (ID: 30652) and the Institut Hospitalo-Universitaire (IHU) Méditerranée Infection, the National Research Agency under the program "Investissements d'avenir", reference ANR-10-IAHU-03, the Région Provence-Alpes-Côte d'Azur and European funding FEDER PRIMI.

**Institutional Review Board Statement:** This study was approved by the Medical and Veterinary Research Ethics Committee at the National Research Centre, Egypt under No. 20147, dated 1 October 2020).

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Acknowledgments:** The authors would like to acknowledge the National Research Centre, Egypt.

**Conflicts of Interest:** The authors declared no conflict of interest with respect to the research, author-
ship and publication of this article.
## References

1. Taylor, L.H.; Latham, S.M.; Woolhouse, M.E.J. Risk factors for human disease emergence. *Philos. Trans. R. Soc. Lond. B* **2001**, *356*, 983-989. [CrossRef] [PubMed]

2. Perumal, A.N.I.; Gunawardene, Y.I.N.S.; Dassanayake, R.S. *Setaria digitata* in advancing our knowledge of human lymphatic filariasis. *J. Helminthol.* **2016**, *90*, 129-138. [CrossRef]

3. Nabie, R.; Spotin, A.; Rouhani, S. Subconjunctival setariasis due to *Setaria equina* infection; A case report and a literature review. *Parasitol. Int.* **2016**, *66*, 930-932. [CrossRef] [PubMed]

4. Laidoudi, Y.; Medkour, H.; Levasseur, A.; Davoust, B.; Mediannikov, O. New Molecular Data on *Filaria* and its *Wolbachia* from Red Howler Monkeys (*Alouattamacconnelli*) in French Guiana-A Preliminary Study. *Pathogens* **2020**, *9*, 626. [CrossRef] [PubMed]

5. Khamesipour, F.; Taktaz-Hafshejani, T.; Tebit, K.E.; Razavi, S.M.; Hosseini, S.R. Prevalence of endo- and ecto-parasites of equines in Iran: A systematic review. *Vet. Med.* **2021**, *7*, 25-34. [CrossRef] [PubMed]

6. El Namaky, A.H.; Hendawy, S.H.M.; Abo-Aziza, F.A.M.; Ashry, H.M. Cytokines and immunoglobulin G response in donkeys with spontaneous *Setaria equina* infection. *Bulg. J. Vet. Med.* **2019**, *22*, 180-189. [CrossRef]

7. Radwan, A.M.; Ahmed, N.E.; Elakabawy, L.M.; Ramadan, M.Y.; Elmadawy, R.S. Prevalence and pathogenesis of some filarial nematodes infecting donkeys in Egypt. *Vet. World* **2016**, *9*, 888-892. [CrossRef]

8. Peng, T.L.; Armiladiana, M.M.; Ruhil, H.H.; Maizan, M.; Choong, S.S. First report of equine *Setaria digitata* (von Linstow 1906) infestation in Malaysia. *Vet. Parasitol. Reg. Stud. Rep.* **2019**, *17*, 100310. [CrossRef]

9. Yu, F.; Liu, B.; Chen, S.; Yi, Z.; Liu, X.; Zhu, Y.; Li, J. First Molecular Confirmation of Equine Ocular *Setaria digitata* in China. *Vet. Sci.* **2021**, *8*, 55. [CrossRef]

10. Tamilmahan, P.; Zama, M.M.S.; Pathak, R.; Muneeswaran, N.S.; Karthik, K. A retrospective study of ocular occurrence in domestic animals: 799 cases. *Vet. World.* **2013**, *6*, 274-276. [CrossRef]

11. Hillyer, L.; Coles, G.; Randle, T. *Setaria equina* in the UK. *Vet. Rec.* **2001**, *149*, 464. [PubMed]

12. Abbas, I.; Al-Araby, M.; Al-Kappany, Y. Molecular characterization of *Setaria equina* infecting donkeys (*Equusasinus*) from Egypt. *Res. J. Parasitol.* **2016**, *11*, 73-78. [CrossRef]

13. Maharana, B.R.; Potliya, S.; Ganguly, A.; Bisla, R.S.; Mishra, C.; Ganguly, I. First report of the isolation and phylogenetic characterization of equine *Setaria digitata* from India based on mitochondrial COI, 12S rDNA, and nuclear ITS2 sequence data. *Parasitol. Res.* **2020**, *119*, 473-481. [CrossRef] [PubMed]

14. Laidoudi, Y.; Bedjaoui, S.; Medkour, H.; Latrofa, M.S.; Mekroud, A.; Bitam, I.; Davoust, B.; Otranto, D.; Mediannikov, O. Molecular Approach for the Diagnosis of Blood and Skin Canine Filarioids. *Microorganisms* **2020**, *8*, 1671. [CrossRef] [PubMed]

15. Harnett, W.; Harnett, M.M. Filarial nematode secreted product ES-62 is an anti-inflammatory agent: Therapeutic potential of small molecule derivatives and ES-62 peptide mimetics. *Clin. Exp. Pharmacol. Physiol.* **2006**, *33*, 511-518. [CrossRef] [PubMed]

16. Abo-Aziza, F.A.M.; Hendawy, S.H.M.; El Namaky, A.H.; Ashry, H.M. Th1/Th2 balance and humoral immune response to potential antigens as early diagnostic method of equine Strongylus nematode infection. *Vet. World* **2017**, *10*, 679-687. [CrossRef]

17. Hendawy, S.H.M. Immunity to gastrointestinal nematodes in ruminants: Effector cell mechanisms and cytokines. *J. Parasit. Dis.* **2018**, *42*, 471-482. [CrossRef]

18. Babu, S.; Nutman, T. Immunology of lymphatic filariasis. *Parasite Immunol.* **2013**, *36*, 2. [CrossRef]

19. Specht, S.; Volkmann, L.; Wynn, T.; Hoerauf, A. Interleukin-10 (IL-10) counter regulates IL-4- dependent effector mechanisms in murine filariasis. *Infect. Immun.* **2004**, *72*, 6287-6293. [CrossRef]

20. Metenou, S.; Dembe, B.; Konate, S.; Dolo, H.; Coulibaly, S.Y.; Coulibaly, Y.I.; Diallo, A.A.; Soumaoro, L.; Coulibaly, M.E.; Sanogo, D.; et al. Patent filarial infection modulates malaria-specific type 1 cytokine responses in an IL-10-dependent manner in a filaria/malaria-coinfected population. *J. Immunol.* **2009**, *183*, 916-924. [CrossRef]

21. Schierack, P.; Lucius, R.; Sonnenburg, B.; Schilling, K.; Hartmann, S. Parasite-specific immunomodulatory functions of filarial cystatin. *Infect. Immun.* **2003**, *71*, 2422-2429. [CrossRef] [PubMed]

22. Cooper, D.; Eleftherianos, I. Parasitic nematode immunomodulatory strategies: Recent advances and perspectives. *Pathogens* **2016**, *5*, 58. [CrossRef] [PubMed]

23. Anuradha, R.; George, P.J.; Hanna, L.E.; Kumaran, P.; Chandrasekaran, V.; Nutman, T.B.; Babu, S. Expansion of parasite-specific CD4+ and CD8+ T cells expressing IL-10 superfamily cytokine members and their regulation in human lymphatic filariasis. *PLoS Negl. Trop. Dis.* **2014**, *8*, e0002762. [CrossRef] [PubMed]

24. Pritchard, J.C.; Lindberg, A.C.; Main, D.C.; Whay, H.R. Assessment of the welfare of working horses, mules and donkeys, using health and behaviour parameters. *Prev. Vet. Med.* **2005**, *69*, 265-283. [CrossRef]

25. Valette, D. *Invisible Workers. The Economic Contributions of Working Donkeys, Horses and Mules to Livelihoods*; The Brooke: Louisville, KY, USA, 2015; pp. 1-23.

26. Abu El-Magd, A.; Ahmed, Z.G. The occurrence of *Setaria equina* in donkey's eyes and their treatment. *Assiut Vet. Med. J.* **1994**, *31*, 86-90.

27. Abd El-Wahab, T.M.; Ashour, A.A. Scanning electron microscopy of the two filariid nematodes *Setaria equina* and *Onchocerca cervicalis* from Kafr El-Sheikh area, Egypt. *Alex. J. Vet. Sci.* **1999**, *15*, 541-547. [CrossRef]

28. Mahmoud, A.E. Laboratory Diagnosis of Filariasis in Assiut Governorate. Master's Thesis, Faculty of Veterinary Medicine, University of Assiut, Assiut, Egypt, 1998.
29. Mahmoud, A.E. A study of *Setaria equina* (Abildgaard, 1789) by light and scanning electron microscopy. *EL-Minia Med. Bullet*. **2006**, *17*, 157–173.

30. Marzok, M.A.; Desouky, A.Y. Ocular infection of donkeys (*Equusasinus*) with *Setaria equina*. *Trop. Anim. Health Prod.* **2009**, *41*, 859–863. [CrossRef]

31. Ahmed, N.E.; El-Akabawy, L.M.; Ramadan, M.Y.; Radwan, A.M.M. Studies on helminth parasites in necropsied donkeys in Egypt. *Benha Vet. Med. J.* **2011**, *1*, 153–162.

32. Abdel Rahman, M.M.I. Morphological and molecular characterization of *Setaria equina* in donkeys. *Beni-Suef Univ. J. Basic Appl. Sci.* **2020**, *9*, 17. [CrossRef]

33. World Health Organization. Global programme to eliminate lymphatic filariasis: Annual report on lymphatic filariasis. *Wkly. Epidemiol. Record*. **2006**, *82*, 361–380.

34. Devi, A.; Sudan, V.; Shanker, D. Phylogenetic characterization of *Setaria equina* and its association with other filarids. *Parasitol. Res.* **2020**, *119*, 4267–4270. [CrossRef] [PubMed]

35. Nagaty, H.F. *Parlitomosazakii* (Filariinae), a new genus and species and its microfilaria from *Leontocebusrosalia*. *J. Egypt. Med. Assoc.* **1935**, *18*, 483–496.

36. Eberhard, M.L.; Orihel, T.C. The genus *Mansonella* (syn. *Tetrapetalonema*): A new classification. *Ann. Parasitol. Hum. Comp.* **1984**, *59*, 483–496. [CrossRef] [PubMed]

37. Sandground, J.H. Report on the nematode parasites collected by the Kelley-Roosevelts expedition to Indo-China with descriptions of several new species. Part 1. Parasites of birds. Part 2. Parasites of mammals. *Z. Parasitenkd.* **1933**, *5*, 542–583. [CrossRef]

38. Wijesundera, W.S.; Chandrasekharan, N.V.; Karunanayake, E.H. A sensitive polymerase chain reaction-based assay for the detection of *Setaria digitata*: The causative organism of cerebrospinal nematodiasis in goats, sheep and horses. *Vet. Parasitol.* **1999**, *81*, 225–233. [CrossRef]

39. BinoSundar, S.T.; D'Souza, P.E. Morphological characterization of *Setaria* worms collected from cattle. *J. Parasit. Dis.* **2015**, *39*, 572–576. [CrossRef]

40. Kaur, D.; Ganai, A.; Parveen, S.; Borkataki, S.; Yadav, A.; Katoch, R.; Godara, R. Occurrence of *Setaria digitata* in a cow. *J. Parasit. Dis.* **2015**, *39*, 477–478. [CrossRef]

41. Abdullah, H.H.A.M.; Amanzougaghene, N.; Dahmana, H.; Louni, M.; Raoult, D.; Mediannikov, O. Multiple vector-borne pathogens of domestic animals in Egypt. *PLoS Negl. Trop. Dis.* **2021**, *15*, e0009767. [CrossRef]

42. Shin, J.; Ahn, K.S.; Suh, G.H.; Kim, H.J.; Jeong, H.S.; Kim, B.S.; Choi, E.; Shin, S.S. First Blindness Cases of Horses Infected with *Setaria Digitata* (Nematoda: Filarioidea) in the Republic of Korea. *Korean J. Parasitol.* **2017**, *55*, 667–671. [CrossRef]

43. Lee, H.; Hwang, H.; Ro, Y.; Kim, J.; Lee, K.; Choi, E.; Bae, Y.; So, B.; Lee, I. *Setaria digitata* was the main cause of equine neurological ataxia in Korea: 50 cases (2015–2016). *J. Vet. Med. Sci.* **2021**, *83*, 869–875. [CrossRef]

44. Genchi, C.; Kramer, L. Subcutaneous dirofilariosis (*Dirofilaria repens*): An infection spreading throughout the old world. *Parasit. Vectors* **2017**, *10*, 517. [CrossRef]

45. Laidoudi, Y.; Davoust, B.; Varloud, M.; Niang, E.H.A.; Fenollar, F.; Mediannikov, O. Development of a multiplex qPCR-based approach for the diagnosis of *Dirofilaria immitis*, *D. repens* and *Acanthocheilonema reconditum*. *Parasit. Vectors* **2020**, *13*, 319. [CrossRef] [PubMed]

46. Laidoudi, Y.; Otranto, D.; Stolowy, N.; Amrane, S.; SanthakumariManoj, R.R.; Polette, L.; Watier-Grillot, S.; Mediannikov, O.; Davoust, B.; L'Ollivier, C. Human and Animal Dirofilariasis in Southeast of France. *Microorganisms* **2021**, *9*, 1544. [CrossRef] [PubMed]

47. Pupić-Bakrač, A.; Pupić-Bakrač, J.; Beck, A.; Jurković, D.; Polkinghorne, A.; Beck, R. *Dirofilaria repens* microfilaremia in humans: Case description and literature review. *One Health* **2021**, *13*, 100306. [CrossRef] [PubMed]

48. Al-Qaoud, K.M.; Pearlman, E.; Hartung, T.; Klukowski, J.; Fleischer, B.; Hoerauf, A. A new mechanism for IL-5-dependent helminth control: Neutrophil accumulation and neutrophilmediated worm encapsulation in murine filariasis are abolished in the absence of IL-5. *Int. Immunol.* **2000**, *12*, 899–908. [CrossRef] [PubMed]

49. Allen, J.E.; Maizels, R.M. Diversity and dialogue in immunity to helminths. *Nat. Rev. Immunol.* **2011**, *11*, 375–388. [CrossRef]

50. Volkmann, L.; Saeftel, M.; Bain, O.; Fischer, K.; Fleischer, B.; Hoerauf, A. Interleukin-4 is essential for the control of microfilariae in murine infection with the filaria Litomosoides sigmodontis. *Infect. Immun.* **2001**, *69*, 2950–2956. [CrossRef]

51. Volkmann, L.; Bain, O.; Saeftel, M.; Specht, S.; Fischer, K.; Brombacher, F.; Matthaei, K.I.; Hoerauf, A. Murine filariasis: Interleukin 4 and interleukin 5 lead to containment of different worm developmental stages. *Med. Microbial. Immunol.* **2003**, *192*, 23–31. [CrossRef]

52. Saeftel, M.; Volkmann, L.; Korten, S.; Brattig, N.; Al-Qaoud, K.; Fleischer, B.; Hoerauf, A. Lack of interferon-gamma confers impaired neutrophil granulocyte function and imparts prolonged 111 survival of adult filarial worms in murine filariasis. *Microbes Infect.* **2001**, *3*, 203–213. [CrossRef]

53. Maizels, R.M.; Balic, A.; Gomez-Escobar, N.; Nair, M.; Taylor, M.D.; Allen, J.E. Helminth parasites-masters of regulation. *Immunol. Rev.* **2004**, *201*, 89–116. [CrossRef]

54. Taylor, M.D.; LeGoff, L.; Harris, A.; Malone, E.; Allen, J.E.; Maizels, R.M. Removal of regulatory T cell activity reverses hyporesponsiveness and leads to filarial parasite clearance in vivo. *J. Immunol.* **2005**, *174*, 4924–4933. [CrossRef] [PubMed]

55. Kwarteng, A.S.T.; Ahuno, S.T.; Akoto, F.O. Killing filarial nematode parasites: Role of treatment options and host immune response. *Infect. Dis. Poverty* **2016**, *5*, 86. [CrossRef] [PubMed]
56. Carvalho, L.; Sun, J.; Kane, C.; Marshall, F.; Krawczyk, C.; Pearce, E.J. Review series on helminthes, immune modulation and the hygiene hypothesis: Mechanisms underlying helminthes modulation of dendritic cell function. *J. Immunol.* **2009**, *126*, 28–34. [CrossRef] [PubMed]

57. Subash, B.; Nutman, T. Proinflammatory cytokines dominate the early immune response to filarial parasites. *J. Immunol.* **2003**, *171*, 6723–6732. [CrossRef]

58. Seidl, A.; Panzer, M.; Voehringer, D. Protective immunity against the gastrointestinal nematode Nippostrongylus brasiliensis requires a broad T-cell receptor repertoire. *J. Immunol.* **2011**, *134*, 214–223. [CrossRef]

59. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. *J. Mol. Biol.* **1990**, *215*, 403–410. [CrossRef]

60. Nakamura, T.; Yamada, K.D.; Tomii, K.; Katoh, K. Parallelization of MAFFT for large-scale multiple sequence alignments. *Bioinformatics* **2018**, *34*, 2490–2492. [CrossRef]

61. Hall, T.; Biosciences, I.; Carlsbad, C. BioEdit: An important software for molecular biology. *GERF Bull. Biosci.* **2011**, *2*, 60–61.

62. Minh, B.Q.; Schmidt, H.A.; Chernomor, O.; Schrempf, D.; Woodhams, M.D.; von Haeseler, A.; Lanfear, R. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. *Mol. Biol. Evol.* **2020**, *37*, 1530–1534. [CrossRef]

63. Hoang, D.T.; Chernomor, O.; von Haeseler, A.; Minh, B.Q.; Vinh, L.S. UFBoot2: Improving the Ultrafast Bootstrap Approximation. *Mol. Biol. Evol.* **2018**, *35*, 518–522. [CrossRef] [PubMed]

64. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.F.; von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phylogenetic estimates. *Nat. Methods* **2017**, *14*, 587–589. [CrossRef] [PubMed]

65. Letunic, I.; Bork, P. Interactive Tree of Life (iTOL) v5: An online tool for phylogenetic tree display and annotation. *Nucleic Acids Res.* **2021**, *49*, 293–296. [CrossRef]

66. Spiegel, M.R. *Statistics*; McGraw-Hill Natl. Book Co.: Singapore, 1981.