# Early Transcriptome Differences Between Pre-Infected and Naïve Kid Goats Infected With Haemonchus contortus

Hadeer Aboshady, Yoann Félicité, Julien Hira, Claude Barbier, Jean-Christophe Bambou

► To cite this version:

Hadeer Aboshady, Yoann Félicité, Julien Hira, Claude Barbier, Jean-Christophe Bambou. Early Transcriptome Differences Between Pre-Infected and Naïve Kid Goats Infected With Haemonchus contortus. Frontiers in Veterinary Science, 2022, 9, pp.873467. 10.3389/fvets.2022.873467. hal-04019718

HAL Id: hal-04019718

https://hal.inrae.fr/hal-04019718

Submitted on 8 Mar 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
<figure><img src="image_2.png" /></figure>

Check for
updates

# Early Transcriptome Differences Between Pre-Infected and Naïve Kid Goats Infected With *Haemonchus contortus*

Hadeer M. Aboshady¹,², Yoann Félicité¹, Julien Hira¹, Claude Barbier³ and Jean-Christophe Bambou¹,²*

¹ INRAE, ASSET, Petit-Bourg, France, ² Department of Animal Production, Faculty of Agriculture, Cairo University, Cairo, Egypt, ³ INRAE, Plateforme Tropicale d'Expérimentation sur l'Animal, Petit-Bourg, France

## OPEN ACCESS

**Edited by:**
Josué Martínez-de la Puente,
University of Granada, Spain

**Reviewed by:**
Amalia Naranjo-Lucena,
Department of Agriculture Food and
the Marine, Ireland
Peter Geldhof,
Ghent University, Belgium
David Bruce Conn,
Berry College, United States

*Correspondence:
Jean-Christophe Bambou
jean-christophe.bambou@inrae.fr

**Specialty section:**
This article was submitted to
Parasitology,
a section of the journal
Frontiers in Veterinary Science

**Received:** 10 February 2022
**Accepted:** 30 May 2022
**Published:** 08 July 2022

**Citation:**
Aboshady HM, Félicité Y, Hira J,
Barbier C and Bambou J-C (2022)
Early Transcriptome Differences
Between Pre-Infected and Naïve Kid
Goats Infected With Haemonchus
contortus. Front. Vet. Sci. 9:873467.
doi: 10.3389/fvets.2022.873467

In small ruminant production, gastrointestinal nematode (GIN) infection is one of the major causes of economic losses. The aim of this study was to compare the abomasal mucosa transcriptome of naïve and pre-infected goats at early time points after *Haemonchus contortus* infection, in order to identify different pathways and upstream regulators involved in the host immune response. Naïve and pre-infected Creole kids were orally infected with 10,000 *H. contortus* infective larvae (L3), and abomasal mucosa was sampled at 0, 4, and 6 days post-infection (dpi). At 6 dpi, all the animals were slaughtered to perform parasite burden counts. The mean number of L4 recovered in naïve kids was more than twice as high as that recovered in the pre-infected ones (5,860 and 2,474 respectively, *p* < 0.001). RNA-seq analysis showed a number of differentially expressed genes (DEGs) very low for both naïve and pre-infected animals when comparing day 0 vs. day 4 post-infection. A total of 2,237 and 3,206 DEGs were identified comparing 0 vs. 6 dpi in naïve and pre-infected animals, respectively. Interestingly, only 18 DEGs were found for the comparison of pre-infected vs. naïve animals at 6 dpi. Ingenuity pathway analysis (IPA) showed that several immune responses were activated in pre-infected compared with naïve animals at 0 and 4 dpi such as Th2 and Th1 pathways, natural killer cell, B cell receptor, IL-2, and IL-15 signaling. On the other hand, both naïve and pre-infected animals showed activation for those pathways comparing 6 vs. 0 dpi, with no difference between them. A similar pattern was recorded for upstream regulator genes which were related to immunity like TNF, IL-1β, IL-2, IL-5, TGFβ1, IFNγ, TCR, IL-18, IL-6, and IL-4. Our results showed that at 0 and 4 dpi the immune response was activated toward Th1 and Th2 pathways in pre-infected kids compared to the naïve ones, however, the same immune response was developed in naïve kids as earlier as 6 dpi. We conclude that repeated *H. contortus* infection in kid goats induced a concomitant early activation of a Th1 and Th2 immune response resulting in the regulation of worm establishment.

**Keywords:** gastrointestinal parasite, goats, mucosal immunity, transcriptome, *Haemonchus contortus*
# INTRODUCTION

Gastrointestinal nematodes (GIN) remain one of the major threats in sheep and goat production worldwide, both in temperate and tropical conditions. The negative impacts of these infections are both, direct and indirect on animal health, welfare, and performance (i.e., losses in milk, meat, and wool), with notably the impact of anthelmintic residues on soil biodiversity (1, 2). Furthermore, anthelmintic resistance is now reported worldwide, thus underlying the need for the development of sustainable and eco-friendly control strategies (3, 4). Host variation in resistance to GIN has been shown to be due to underlying genetic diversity which provides a promising control strategy through genetic selection (5, 6). Understanding molecular mechanisms underlying host response to GIN is the main step for developing appropriate genetic selection programs.

Numerous sheep studies showed that a strong Th2 response against GIN infections, characterized by parasite-specific local and systemic antibody response, is a key element for the development of a protective immune response and resistance to infection (7, 8). The fewer studies that have investigated the goat immune response to GIN, showed marked differences with the response in sheep (9, 10). RNA-sequencing (RNA-seq) technique is a powerful high-throughput next-generation sequencing technology for transcriptome analysis to identify and quantify genes expressed for complex quantitative traits like a response to infection (11). Indeed, numerous studies have been undertaken in sheep to identify genes and pathways, associated with the host response to GIN, from infected tissues and lymph node transcriptome (12–15). Results from these studies confirmed the role of the immune response, in particular the Th2 response for the development of host resistance to GIN. Meanwhile, the mechanisms identified in goats in response to GIN infections suggest a role of Th1, Th2, and TGF-β1 genes (16).

The results from sheep studies have also shown that the response of naïve animals to GIN infection compared with pre-infected animals is qualitatively and quantitatively different (17–19). Thus, it has been suggested that the host response developed over time to infections, and the main immune mechanism involved is acquired and not innate (15). In goats, the difference between naïve and pre-infected animals in response to GIN infection and the underlying mechanisms has been poorly investigated.

Recently, it has been shown that the response during the early time of infection is a key component of the protective host-immune response in sheep (15, 20) and goats (16). Thus, studies focused on this early time after infection and comparing naïve and pre-infected animals need to be conducted. The objective of the present study was to characterize the mucosal host response in naïve and pre-infected goats to *Haemonchus contortus* infection at very early times after infection using RNA-seq technology.

# MATERIALS AND METHODS

## Ethics Approval and Consent to Participate

All measurements and observations on animals were carried out in accordance with the current law on animal experimentation and ethics and approved by the French Ministry of Agriculture (authorization number: HC-69-2014-1) after evaluation by the Animal Care and Use Committee of French West Indies and Guyana (Comité d'Ethique en Matière d'Expérimentation Animale des Antilles et de la Guyane, C2EA-69).

## Animals, Management, and Experimental Design

This experiment was conducted at the INRAE PTEA (Plateforme Tropicale d'Expérimentation sur l'Animal) experimental farm in Guadeloupe (16° 20' North latitude, 61° 30' West longitude). All the animals owned by the INRAE PTEA were reared in this experimental farm since 1980. For this experiment, a total of 8 unrelated male kids, from 8 distinct sire families, were reared indoors and fed with parasite-free hay to avoid GIN infection. At 4- and 6-months-old, kids of the pre-infected group (n = 4) were experimentally infected with a single dose of 10,000 *H. contortus* L3 with a 10 ml syringe containing 10 ml of a suspension of L3 at 1,000 L3/ml in tap water. The naïve group (n = 4 animals) received 10 ml of tap water. After 5 weeks of infection for each challenge, the kids were drenched with moxidectine (Cydectine®, Fort Dodge Veterinaria S.A., Tours, France, 300 µg/kg). To avoid any bias due to moxidectine administration, naïve kids were drenched at the same time. For each challenge, Faecal Egg counts (FEC) was measured at week 4 and 5 post-infection. Approximately 10 g of feces were collected in plastic tubes directly from the rectum of each animal and transported from the experimental facility to the laboratory in refrigerated vials. The samples were individually analyzed using a modified McMaster method for rapid determination and FEC was expressed as the number of eggs/g feces. The mean FEC for the first and the second challenge of the pre-infected kids, at 4 and 6-month-old were 3567 ± 834 and 2765 ± 843, respectively. The FEC remained at zero for the naïve kids. At 8 months old, a fistula was surgically implanted in the abomasum of each animal to allow abomasal mucosa sampling at 0, 3-, and 6-days post-infection (dpi). At 9 months old, all the animals (pre-infected and naïves) were experimentally infected with a single dose of 10,000 *H. contortus* L3 as described above. Mucosal biopsies were collected at 0, 3-, and 6-days post-infection, then at 6 dpi, all the animals were humanly euthanized with a captive bolt stunning gun against their head, followed by exsanguination, in accordance with the procedures authorized in animal experiments. The abomasum was recovered for parasite counts.

## Surgical Procedure

The custom-designed abomasal cannula consisted of a flexible plastic tube with a length of 7 cm and a diameter of 2 cm with a rounded base of 4 cm in diameter. This flexible plastic was chosen to limit the possibility of mechanical abrasion of the mucosal
surface of the abomasum. The animals were fasted 16 h before cannula insertion surgery. The animals were premedicated with ketamine (2 mg/kg Intravenous (IV), Le Vet Pharma, Wilgenweg, Netherlands), xylazine (0.2 mg mg/kg Intramuscular (IM), Le Vet Pharma, Wilgenweg, Neitherlands), and oxytetracycline (20 mg/kg IM, Eurovet Animal Health, Handelsweg, Neitherlands). The animals were positioned in left lateral recumbency. Skin over the surgical site was shaved and prepped with povidone iodine (Vétédine, Laboratoire Vetoquinol S.A., Lure, France). A ventral midline incision was made to locate and externalize the abomasum. A 3 cm purse-string suture (Silk 2-0) was placed midway between the lesser and greater curvature and a stab incision was made in the center to insert the cannula. Then, the purse-string suture was tightened and tied off. To maintain the abomasum in an anatomically correct position, another stab incision was made in the abdominal wall at 10 cm from the laparotomy incision on the right paramedian area to enable the cannula to be passed freely through. An external flange was placed over the external part of the cannula and fixed with an adhesive fabric plaster strip. A sterile compress was inserted into the cannula as a stopper. After the surgical procedure, all the animals were housed individually with free access to fresh water and hay.

## Biopsy Sampling Procedure

Biopsy specimens were taken from the abomasal mucosa using a flexible endoscope (FG-24V, Pentax, France). The biopsies samples of 2 mm × 2 mm × 2 mm taken with the endoscopic forceps with window (model KW1815S) were quickly snap frozen into liquid nitrogen and stored at -80°C until RNA extraction. The animals were restrained in a harness made with a surgical drape allowing animal legs to protrude and which exposed the cannula. No sedation was used since no signs of discomfort or pain were observed during or after the procedure. The sterile compress inserted into the cannula was removed and the abomasal contents were collected. The endoscope was introduced into the abomasal lumen and 3 biopsies per animal and per time points were taken from the abomasal folds of the fundic mucosa. At each time point the whole fundic mucosa was observed and no sign of mucosal injury due to the previous sampling was observed.

## RNA Extraction and Sequencing

Total RNA was extracted using the NucleoSpin® RNA isolation kit (Macherey-Nagel, Hoerdt, France) following the manufacturer's instructions, except that DNase digestion was performed with twice the indicated amount of enzyme according to our adjustments within our laboratory. The total RNA concentration was measured with NanoDrop 2000 (ThermoScientific TM, France). The RNA integrity was verified using an Agilent Bioanalyzer 2100 (Agilent Technologies, France) with an RNA integrity number > 8.5. The extracted total RNA was stored at -80°C until sequencing.

High-quality RNA from all samples was processed for the preparation of cDNA libraries using an Illumina TruSeq RNA sample prep kit for mRNA analysis following Illumina's

protocols. After quality control and quantification, cDNA libraries were pooled in groups of 6 and sequenced on 5 lanes on the HiSeqTM 2000 (Illumina® NEB, USA) to obtain approximately 50 million reads (100 bp paired-end) for each sample with insert sizes ranging from 200 to 400 base pairs.

## Bioinformatics and Data Analysis

FASTQC was used to check the quality control of raw reads in FASTQ format and the Q20, Q30, and GC contents of the clean data were calculated. Transcript quantifications were made using the Salmon software (version 0.9.1) (21). The index was built within Salmon using NCBI RefSeq reference transcript of the *Capra hircus* genome (assembly ARS1). The reads from each sample were mapped to the same index and quantified. Unix commands were used to obtain corresponding gene and transcript identifiers from the NCBI RefSeq annotation of the *Capra hircus* (ARS1). Using these identifiers, the tximport (version 1.8.0) package was used to import data into the R software (v3.5.1) and summarize the TPM estimates obtained from the Salmon tool of all samples at the gene level (22). This process produced a global count file on which the statistical analyses were performed. Genes with low expression levels were removed using a threshold of greater than or equal to 5 counts across samples. Read counts table were used to identify differentially expressed genes (DEGs) using the Bioconductor package DESeq2 within R (23).

Seven comparisons were performed; three comparing pre-infected vs. naïve at day 0, 4, and 6 dpi, two comparing day 4 and 6 vs. 0 dpi in naïve kids, and another two comparing day 4 and 6 vs. 0 dpi in pre-infected kids. Genes were filtered using a Benjamini and Hochberg false discovery rate (FDR) of < 0.001 to account for multiple testing. Final DEGs were determined on the basis of their fold change values to be log2 ≥ 1.0 for upregulated genes and ≤ -1.0 for downregulated genes. DEGs in each comparison were analyzed for canonical pathways and regulator effects using Ingenuity pathway analysis (IPA) software (Ingenuity Systems, Redwood City, CA). DEGs list of each comparison was uploaded to IPA and compared with humans orthologs data set generated by IPA to calculate the matching percentage.

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Mean number of *Haemonchus contortus* fourth-stage (L4) larvae recovered from naïve and pre-infected kids at 6 days post-infection. Means are significantly different (*p* < 0.001).</figcaption></figure>
## Parasite Counts

At slaughter (6 dpi, after biopsies sampling), the abomasum of each kid was isolated with contents for parasite counts. The contents of the abomasum were collected individually and the abomasum was washed with warm distilled water and scraped with a microscope slide to recover all established larvae. The contents and the wash water were stored at 4°C until counting. The parasites were collected, counted, and sorted according to the method previously described (24).

## RESULTS

The number of L4 recovered from naïve and pre-infected kids at 6 dpi is presented in **Figure 1**. The mean number of L4 recovered in naïve kids was more than twice as high as that recovered in the pre-infected kid (5,860 and 2,474, respectively, *p* < 0.001). No other stage of larval development was found in either group.

RNA sequencing for each sample aligned to the reference *Capra hircus* genome (assembly ARSI) resulted in an average of 48.4 ± 1.8 million reads per sample. These reads correspond to 23,258 genes of the goat genome. A total of 17,311 out of the 23,258 annotated genes showed at least 5 read counts per raw and were used in the subsequent analysis.

**Table 1** shows the number of DEGs for each comparison. Comparing pre-infected versus naïve kids, the number of DEG were 2,793 and 1,231 at 0 and 4 dpi, respectively. Meanwhile, at 6 dpi there were only 18 DEGs between pre-infected and naïve. In naïve kids, a low number of DEG was observed when comparing 4 and 0 dpi (12 DEG), in contrast, a high number of DEG was observed when comparing 6 and 0 dpi (2237 DEG). The same pattern was observed for pre-infected kids; a low number of DEG when comparing 4 and 0 dpi (9 DEG) and a high number of DEG when comparing 6 and 0 dpi (3206 DEG). Whatever the comparison, the number of upregulated and downregulated genes were similar, except for the comparison between pre-infected and naïve kids at 4 dpi, the number of downregulated genes was 218 compared with 1,013 for the upregulated ones. The range for log₂-fold change was wider for the pre-infected kids comparing 6 and 0 dpi (-25.93 and 16.61)

than the naïve kids for the same comparison (-9.60, 13.24). DEG from different comparisons showed human orthologs between 75 and 91%.

**Table 2** presents the gene symbol and log₂-fold change for comparison showing a low number of DEG (pre-infected vs. naïve at 6 dpi, naïve 0 and 4 dpi, and pre-infected 0 and 4 dpi). No functional analysis was possible for these comparisons due to the low number of DEGs found. Nonetheless, the comparison of 4 and 0 dpi in naïve kids showed activation of the expression of different immune-related genes (e.g., TNFRSF12A, IL1B, CXCL8, and OLFM4; **Table 2**). Moreover, four genes involved in the IL-17 signaling pathways were found in this comparison: CXCL8, IL1B, MMP1, and MMP3. Meanwhile, pre-infected kids did not show the same pattern, and no pathways were suggested by the genes found. Besides, when comparing pre-infected and naïve kids at 6 dpi, most DEGs were not related to the immune response except MUC5AC, CCL20, KLRD1, and FAM3B.

The top 20 canonical pathways of DEG for different comparisons are presented in **Figure 2**. The majority of top canonical pathways were related to the immune response. Th1 and Th2 pathways were found in the top 10 pathways activated by 6 compared with 0 dpi, in both naïve and pre-infected kids. In general, all top canonical pathways had the same pattern in naïve and pre-infected kids when comparing different dpi. The activation of almost all pathways was higher in pre-infected compared with naïve kids at 0 and 4 dpi. However, no difference in pathways activated at 6 dpi was observed between pre-infected and naïve kids. The fatty acid β-oxidation pathway was highly downregulated in pre-infected animals when comparing 6 and 0 dpi, while it was upregulated comparing pre-infected and naïve at 0 and 4 dpi.

The top 20 upstream regulators for different comparisons are presented in **Figure 3**. At 0 dpi pre-infected kids had already activated upstream regulators immune genes like TNF, IL1B, IL-2, IL-5, IFNG, CD3, TCR, IL-8, IL-6, and IL-4. On the other hand, these upstream regulators genes showed high activation in naïve animals when comparing 6 and 0 dpi, with activation score higher than those observed in pre-infected animals comparing the same days.

<table><caption>**TABLE 1** | The number of differentially expressed genes (DEGs) for the different comparisons including log₂-fold change and the proportion of genes with human orthologs.</caption><thead><tr><th rowspan="2">Comparison</th><th colspan="2">DEG</th><th rowspan="2">Log₂-fold change</th><th rowspan="2">Human orthologs</th></tr><tr><th>Up-regulated</th><th>Down-regulated</th></tr></thead><tbody><tr><td>Pre-inf vs. Naïve 0 dpi</td><td>1618</td><td>1175</td><td>-12.03, 25.04</td><td>78.2%</td></tr><tr><td>Pre-inf vs. Naïve 4 dpi</td><td>1013</td><td>218</td><td>-7.40, 18.63</td><td>75.6%</td></tr><tr><td>Pre-inf vs. Naïve 6 dpi</td><td>8</td><td>10</td><td>-8.04, 2.81</td><td>77.8%</td></tr><tr><td>Naïve 4 vs. 0 dpi</td><td>11</td><td>1</td><td>-1.37, 6.99</td><td>91.6%</td></tr><tr><td>Naïve 6 vs. 0 dpi</td><td>1136</td><td>1101</td><td>-9.60, 13.24</td><td>83.4%</td></tr><tr><td>Pre-inf 4 vs. 0 dpi</td><td>9</td><td>0</td><td>3.18, 9.24</td><td>88.9%</td></tr><tr><td>Pre-inf 6 vs. 0 dpi</td><td>1331</td><td>1875</td><td>-25.93, 16.61</td><td>79.1%</td></tr></tbody></table>

Pre-inf, Pre-infected kids.
dpi, days post-infection.
<table><caption>**TABLE 2** | Differentially expressed genes (DEGs) and log₂-fold change for pre-infected and naïve kids at 6 days post-infection (dpi), naïve and pre-infected comparing 0 and 4 dpi.</caption><thead><tr><th colspan="2">Pre-inf vs. Naïve at 6 dpi</th><th colspan="2">Naïve 4 vs. 0 dpi</th><th colspan="2">Pre-inf 4 vs. 0 dpi</th></tr><tr><th>Gene</th><th>Log₂ -old change</th><th>Gene</th><th>Log₂-fold change</th><th>Gene</th><th>Log₂-fold change</th></tr></thead><tbody><tr><td>RPL9</td><td>-8.04</td><td>PARVB</td><td>-1.37</td><td>FER1L6</td><td>3.18</td></tr><tr><td>ACTB-1</td><td>-3.44</td><td>ITGA2</td><td>2.06</td><td>B3GNT6</td><td>4.65</td></tr><tr><td>MUC5AC</td><td>-3.37</td><td>TNFRSF12A</td><td>2.92</td><td>LOC106502720</td><td>4.70</td></tr><tr><td>CCL20</td><td>-3.05</td><td>PLAUR</td><td>3.00</td><td>VSIG1</td><td>4.80</td></tr><tr><td>LOC102180547</td><td>-2.64</td><td>IL1B</td><td>3.14</td><td>SLC6A14</td><td>5.48</td></tr><tr><td>RAB3B</td><td>-1.95</td><td>IFIT1</td><td>3.30</td><td>CHST4</td><td>5.49</td></tr><tr><td>FAM3B</td><td>-1.40</td><td>CDA</td><td>3.33</td><td>SOX21</td><td>7.78</td></tr><tr><td>TMC5</td><td>-1.36</td><td>MMP3</td><td>3.46</td><td>GKN1</td><td>8.83</td></tr><tr><td>DGAT2</td><td>-1.13</td><td>CXCL8</td><td>3.50</td><td>GKN2</td><td>9.24</td></tr><tr><td>VILL</td><td>-1.09</td><td>DIO2</td><td>3.62</td><td></td><td></td></tr><tr><td>STMN2</td><td>1.15</td><td>MMP1</td><td>4.69</td><td></td><td></td></tr><tr><td>IGF1</td><td>1.21</td><td>OLFM4</td><td>6.99</td><td></td><td></td></tr><tr><td>MCTP1</td><td>1.29</td><td></td><td></td><td></td><td></td></tr><tr><td>CPXM2</td><td>1.30</td><td></td><td></td><td></td><td></td></tr><tr><td>TRPC4</td><td>1.36</td><td></td><td></td><td></td><td></td></tr><tr><td>SLC2A3</td><td>1.63</td><td></td><td></td><td></td><td></td></tr><tr><td>P2RX1</td><td>1.76</td><td></td><td></td><td></td><td></td></tr><tr><td>KLRD1</td><td>2.81</td><td></td><td></td><td></td><td></td></tr></tbody></table>

Pre-inf, Pre-infected kids.
dpi, days post-infection.

<figure><img src="image_4.png" /></figure>
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Top 20 upstream regulators of differently expressed genes for naïve and pre-infected (Pre-inf) kids comparing 0 vs. 4 and 6 days post-infection (dpi) and pre-infected vs. naïve at 0, 4, and 6 dpi.</figcaption></figure>

DISCUSSION

Numerous studies, mainly in sheep, suggested that the protective response against GIN infection would be a set of several parasite stage-specific mechanisms that differ both between breeds and divergently selected lines (resistant and susceptible) (25, 26). Thus, to better characterize the mechanism involved at the early stage of infection in goats, here were developed a model to monitor the same animals the response at 0, 3, and 6 dpi with *H. contortus*, and we compared naïve and pre-infected kids.

The lower establishment rate of infective larvae at 6 dpi observed in pre-infected kids suggested that a stage-specific mechanism developed after repeated infections targeted the early stage of infection. A similar study, comparing two lines of Merino sheep, divergently selected for resistance and susceptibility to GIN infections, was conducted recently in naïve and pre-infected lambs (27). Unfortunately, no measure of worm counts was performed, but the authors showed, that the number of DEG between the abomasal transcriptome of resistant and susceptible lambs was higher in the naïve compared with the pre-infected ones, on the third day after infection. In addition, the primary host response against *H. contortus* infection in the resistant line was characterized by the activation of several B cells-related canonical pathways. Indeed, many studies on sheep have shown that the major host-immune mechanism was acquired and not innate, with protective immunity developing over time in response to infection (15, 19, 28). In the present study, our results showed that the number of DEG between abomasal

transcriptome of pre-infected and naïve kid goats was higher at the beginning of infection (day 0 and 4 post-infection), while there were very fewer DEGs between pre-infected and naïve goats at 6 dpi. These results are in accordance with previous findings showing a very low number of DEG at 3 dpi (ranging from 2 to 134), between primary and tertiary challenged sheep from two flocks selected for divergent responses to *H. contortus* and *Trichostrongylus colubriformis* (29). The slow development of the adaptive immune response or the high similarity of the first steps of the adaptive and primary immune responses have been hypothesized by the authors. Here, the first hypothesis would be unlikely given the lower establishment rate found in pre-infected kids. The high similarity of the adaptive and the primary immune response at 6 dpi is in fact observed, however, it is important to specify that it is at the transcriptional level. Further functional analysis (e.g., immunohistochemistry, proteomics) would have been necessary for a more detailed characterization of the mechanisms involved.

In pre-infected goats compared to the naïve ones, different immune pathways and upstream regulators were activated at the beginning of infection (0 and 4 dpi). Thereafter, at 6 dpi the same immune pathways and upstream regulators were activated in naïve goats. Our results comparing 6 and 0 dpi, showed that naïve and pre-infected animals had a high number of DEG involved in the activation of different immune responses including Th1, Th2, and B cell receptor signaling. The role of the Th1/Th2 balance in the development of a protective immune response against GIN infection has been previously discussed in sheep and goats (16, 30, 31). Sheep studies suggested that the Th1 response would be rather associated with susceptibility and the Th2 response with resistance (30, 31). Meanwhile, in goats, a concomitant activation of the Th1 and Th2 responses was observed both in resistant and susceptible animals, with only differences in the timing of the immune response between the two lines (16). Here, in line with this previous study, both Th1 and Th2 responses were activated in naïve and pre-infected animals, with earlier activation in pre-infected animals compared with the naïve ones. A similar pattern was recorded for upstream regulator genes which were related to immunity like TNF, IL1β, IL2, IL5, TGFβ1, IFNγ, TCR, IL18, IL6, and IL4.

Moreover, we recently identified genomic variants between resistant and susceptible goats to GIN, in 78% of genes controlling the T-cell receptor-signaling pathway (32). The T-cell receptor-signaling pathway would play a role in the development of host immunity against GIN infection. In this study, we found a B-cell receptor signaling pathway in the top activated pathway in naïve and pre-infected goats comparing 6 and 0 dpi. Thus, our results strongly suggested that in our biological model of Creole goats, receptor signaling pathways would be one key mechanism of the immune response against *H. contortus* activated both after the primary, or repeated infections. One difference between the immune responses after primary or repeated infection would be the kinetics of the mechanisms involved, with an earlier activation in pre-infected animals. Hence, we conclude that repeated *H. contortus* infection in kid goats induced a concomitant early activation of a Th1 and Th2 immune response resulting in the regulation of worm establishment.
DATA AVAILABILITY STATEMENT

The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih. gov/biosample?LinkName=bioproject_biosample_all&from_ uid=808243.

ETHICS STATEMENT

The animal study was reviewed and approved by Comité d'Ethique en Matière d'Expérimentation Animale des Antilles et de la Guyane, C2EA-69.

AUTHOR CONTRIBUTIONS

HA and JCB conceived and designed the experiments, performed the bioinformatic and statistical analysis, and wrote the paper. JCB, YF, JH, and CB collected samples for the hematological

and parasitological analysis. YF and JH performed the laboratory analysis. All authors read and approved the final manuscript.

FUNDING

This study was funded by the Project MALIN (Maladies Infectieuses en Milieu Insulaire Tropical, La Région Guadeloupe and Fonds Européens FEDER). HA was supported by a doctoral fellowship from the Animal Genetics Division of INRAE and the European Project EGS-ABG (European Graduate School in Animal Breeding and Genetics).

ACKNOWLEDGMENTS

The authors want to thank the PTEA Petits-ruminants team of Duclos for the care and handling of the animals.

REFERENCES

1. Mavrot F, Hertzberg H, Torgerson P. Effect of gastro-intestinal nematode infection on sheep performance: a systematic review and meta-analysis. *Parasit Vectors*. (2015) 8:557. doi: 10.1186/s13071- 015-1164-z

2. Adler N, Bachmann J, Blanckenhorn WU, Floate KD, Jensen J, Rombke J. Effects of ivermectin application on the diversity and function of dung and soil fauna: regulatory and scientific background information. *Environ Toxicol Chem*. (2016) 35:1914-23. doi: 10.1002/etc.3308

3. Jabbar A, Iqbal Z, Kerboeuf D, Muhammad G, Khan MN, Afaq M. Anthelmintic resistance: the state of play revisited. *Life Sci*. (2006) 79:2413-31. doi: 10.1016/j.lfs.2006.08.010

4. Geurden T, Hoste H, Jacquiet P, Traversa D, Sotiraki S, Frangipane A, et al. Anthelmintic resistance and multidrug resistance in sheep gastro-intestinal nematodes in France, Greece and Italy. *Vet Parasitol*. (2014) 201:59-66. doi: 10.1016/j.vetpar.2014.01.016

5. Bishop SC. Possibilities to breed for resistance to nematode parasite infections in small ruminants in tropical production systems. *Animal*. (2012) 6:741-7. doi: 10.1017/S1751731111000681

6. McRae KM, McEwan JC, Dodds KG, Gemmell NJ. Signatures of selection in sheep bred for resistance or susceptibility to gastrointestinal nematodes. *BMC Genomics*. (2014) 15:637. doi: 10.1186/1471-2164-15-637

7. Aboshady HM, Johansson A, Jonas E, Stear MJ, Bambou J-C. Immunoglobulins as biomarkers for gastrointestinal nematodes resistance in small ruminants: a systematic review and meta-analysis. *Sci Rep*. (2020) 10:7765. doi: 10.1038/s41598-020-64775-x

8. Escribano C, Saravia A, Costa M, Castells D, Ciappesoni G, Riet-Correa F, et al. Resistance to Haemonchus contortus in Corriedale sheep is associated to high parasite-specific IgA titer and a systemic Th2 immune response. *Sci Rep*. (2019) 9:1-10. doi: 10.1038/s41598-019-55447-6

9. Pomroy WE, Lambert MG, Betteridge K. Comparison of faecal strongylate egg counts of goats and sheep on the same pasture. *N Z Vet J*. (1986) 34:36-7. doi: 10.1080/00480169.1986.35272

10. Hoste H, Torres-Acosta JF, Aguilar-Caballero AJ. Nutrition-parasite interactions in goats: is immunoregulation involved in the control of gastrointestinal nematodes? *Parasite Immunol*. (2008) 30:79-88. doi: 10.1111/j.1365-3024.2007.00987.x

11. Costa V, Aprile M, Esposito R, Ciccodicola A. RNA-Seq and human complex diseases: recent accomplishments and future perspectives. *Eur J Hum Genet*. (2013) 21:134-42. doi: 10.1038/ejhg.2012.129

12. Gossner A, Wilkie H, Joshi A, Hopkins J. Exploring the abomasal lymph node transcriptome for genes associated with resistance to

the sheep nematode Teladorsagia circumcincta. *Vet Res*. (2013) 44:1-13. doi: 10.1186/1297-9716-44-68

13. Ahmed AM, Good B, Hanrahan JP, McGettigan P, Browne J, Keane OM, et al. Variation in the Ovine abomasal lymph node transcriptome between breeds known to differ in resistance to the gastrointestinal nematode. *PLoS ONE*. (2015) 10:e0124823. doi: 10.1371/journal.pone.0124823

14. Guo Z, González JF, Hernandez JN, McNeilly TN, Corripio-Miyar Y, Frew D, et al. Possible mechanisms of host resistance to Haemonchus contortus infection in sheep breeds native to the Canary Islands. *Sci Rep*. (2016) 6:1-14. doi: 10.1038/srep26200

15. McRae KM, Good B, Hanrahan JP, McCabe MS, Cormican P, Sweeney T, et al. Transcriptional profiling of the ovine abomasal lymph node reveals a role for timing of the immune response in gastrointestinal nematode resistance. *Vet Parasitol*. (2016) 224:96-108. doi: 10.1016/j.vetpar.2016.05.014

16. Aboshady HM, Mandonnet N, Félicité Y, Hira J, Fourcot A, Barbier C, et al. Dynamic transcriptomic changes of goat abomasal mucosa in response to Haemonchus contortus infection. *Vet Res*. (2020) 51:44. doi: 10.1186/s13567-020-00768-y

17. Gomez-Munoz MT, Cuquerella M, Gomez-Iglesias LA, Mendez S, Fernandez-Perez FJ, De La Fuente C, et al. Serum antibody response of Castellana sheep to Haemonchus contortus infection and challenge: relationship to abomasal worm burdens. *Vet Parasitol*. (1999) 81:281-93. doi: 10.1016/S0304-4017(98)00260-X

18. Schallig HDFH, van Leeuwen MAW, Hendrikx WML. Isotype-specific serum antibody responses of sheep to Haemonchus contortus antigens. *Vet Parasitol*. (1995) 56:149-62. doi: 10.1016/0304-4017(94)00675-3

19. Lacroux C, Nguyen THC, Andreoletti O, Prevot F, Grisez C, Bergeaud J-P, et al. Haemonchus contortus (Nematoda: Trichostrongylidae) infection in lambs elicits an unequivocal Th2 immune response. *Vet Res*. (2006) 37:607-22. doi: 10.1051/vetres:2006022

20. El-Ashram S, Al Nasr I, El-Kemary M, Mehmood R, Hu M, Suo X. Early and late gene expression profiles of the ovine mucosa in response to Haemonchus contortus infection employing Illumina RNA-seq technology. *Parasitol Int*. (2017) 66:681-92. doi: 10.1016/j.parint.2017.05.007

21. Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon: fast and bias-aware quantification of transcript expression using dual-phase inference. *Nat Methods*. (2017) 14:417-9. doi: 10.1038/nmeth.4197

22. Soneson C, Love MI, Robinson MD. Differential analyses for RNA-seq : transcript-level estimates improve gene-level inferences. *F1000Research*. (2015) 4:1521. doi: 10.12688/f1000research.7563.1

23. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biol*. (2014) 15:550. doi: 10.1186/s13059-014-0550-8
24. Bambou JC, Larcher T, Cei W, Dumoulin PJ, Mandonnet N. Effect of Experimental Infection with Haemonchus contortus on parasitological and local cellular responses in resistant and susceptible young creole goats. *Biomed Res Int.* (2013) 2013:902759. doi: 10.1155/2013/902759

25. Gonzalez JF, Hernandez A, Molina JM, Fernandez A, Raadsma HW, Meeusen ENT, et al. Comparative experimental Haemonchus contortus infection of two sheep breeds native to the Canary Islands. *Vet Parasitol.* (2008) 153:374–8. doi: 10.1016/j.vetpar.2008.02.019

26. Balic A, Bowles VM, Meeusen ENT. Mechanisms of immunity to Haemonchus contortus infection in sheep. *Parasite Immunol.* (2002) 24:39–46. doi: 10.1046/j.0141-9838.2001.00432.x

27. Zhang R, Liu F, Hunt P, Li C, Zhang L, Ingham A, et al. Transcriptome analysis unraveled potential mechanisms of resistance to Haemonchus contortus infection in Merino sheep populations bred for parasite resistance. *Vet Res.* (2019) 50:1–13. doi: 10.1186/s13567-019-0622-6

28. Stear MJ, Park M, Bishop SC. The Key Components of Resistance to Ostertagia circumcincta in Lambs. *Parasitol Today*. (1996) 12:438–41. doi: 10.1016/0169-4758(96)10069-7

29. Liu J, Tan M, Xu X, Shen T, Zhou Z, Hunt PW, Zhang R. From innate to adaptive immunity: Abomasal transcriptomic responses of merino sheep to Haemonchus contortus infection. *Mol Biochem Parasitol.* (2021) 246:111424. doi: 10.1016/j.molbiopara.2021.111424

30. Zaros LG, Neves MRM, Benvenuti CL, Navarro AMC, Sider LH, Coutinho LL, et al. Response of resistant and susceptible Brazilian Somalis crossbreed sheep naturally infected by Haemonchus contortus. *Parasitol Res.* (2014) 113:1155–61. doi: 10.1007/s00436-014-3753-8

31. Venturina VM, Gossner AG, Hopkins J. The immunology and genetics of resistance of sheep to Teladorsagia circumcincta. *Vet Res Commun.* (2013) 37:171–81. doi: 10.1007/s11259-013-9559-9

32. Aboshady HM, Mandonnet N, Johansson AM, Jonas E, Bambou J-C. Genomic variants from RNA-seq for goats resistant or susceptible to gastrointestinal nematode infection. *PLoS ONE*. (2021) 16:e248405. doi: 10.1371/journal.pone.0248405

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2022 Aboshady, Félicité, Hira, Barbier and Bambou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.