# Transcriptional profiling of midguts prepared from Trypanosoma/T. congolense-positive Glossina palpalis palpalis collected from two distinct Cameroonian foci: coordinated signatures of the midguts' remodeling as T. congolense-supportive niches

Jean Tsagmo Ngoue, Flobert Njiokou, Béatrice Loriod, Ginette

Kame-Ngasse, Nicolas Fernandez-Nunez, Claire Rioualen, Jacques van Helden,

Anne Geiger

► To cite this version:

Jean Tsagmo Ngoue, Flobert Njiokou, Béatrice Loriod, Ginette Kame-Ngasse, Nicolas Fernandez-Nunez, et al.. Transcriptional profiling of midguts prepared from Trypanosoma/T. congolense-positive Glossina palpalis palpalis collected from two distinct Cameroonian foci: coordinated signatures of the midguts' remodeling as T. congolense-supportive niches. Frontiers in Immunology, 2017, 8, 10.3389/fimmu.2017.00876. [hal-01874935](https://hal HAL Id: hal-01874935)

Submitted on 18 Oct 2018

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

Check for
updates

# Transcriptional Profiling of Midguts Prepared from *Trypanosoma* *T. congolense*-Positive *Glossina palpalis* Collected from Two Distinct Cameroonian Foci: Coordinated Signatures of the Midguts' Remodeling As *T. congolense*-Supportive Niches

OPEN ACCESS

Edited by:

Alexandre Morrot,
Federal University of
Rio de Janeiro, Brazil

Reviewed by:

Geneviève Milon,
Institut Pasteur, France
Carl De Trez,
Vrije Universiteit Brussel, Belgium

*Correspondence:

Anne Geiger
anne.geiger@ird.fr

Specialty section:

This article was submitted
to Microbial Immunology,
a section of the journal
Frontiers in Immunology

Received: 01 June 2017

Accepted: 10 July 2017

Published: 28 July 2017

Citation:

Tsagmo Ngoue JM, Njiokou F,
Loriod B, Kame-Ngasse G,
Fernandez-Nunez N, Rioualen C,
van Helden J and Geiger A (2017)
*Transcriptional Profiling of Midguts
Prepared from Trypanosoma/
T. congolense-Positive Glossina
palpalis palpalis Collected from
Two Distinct Cameroonian Foci:
Coordinated Signatures of the
Midguts' Remodeling As
T. congolense-Supportive Niches.*
Front. Immunol. 8:876.
doi: 10.3389/fimmu.2017.00876

Jean M. Tsagmo Ngoue<fnref n="1" />,², Flobert Njiokou<fnref n="1" />, Béatrice Loriod³, Ginette Kame-Ngasse<fnref n="1" />, Nicolas Fernandez-Nunez³, Claire Rioualen³, Jacques van Helden³ and Anne Geiger²*

Our previous transcriptomic analysis of *Glossina palpalis gambiensis* experimentally infected or not with *Trypanosoma brucei gambiense* aimed to detect differentially expressed genes (DEGs) associated with infection. Specifically, we selected candidate genes governing tsetse fly vector competence that could be used in the context of an anti-vector strategy, to control human and/or animal trypanosomiasis. The present study aimed to verify whether gene expression in field tsetse flies (*G. p. palpalis*) is modified in response to natural infection by trypanosomes (*T. congolense*), as reported when insectary-raised flies (*G. p. gambiensis*) are experimentally infected with *T. b. gambiense*. This was achieved using the RNA-seq approach, which identified 524 DEGs in infected vs. non-infected tsetse flies, including 285 downregulated genes and 239 upregulated genes (identified using DESeq2). Several of these genes were highly differentially expressed, with log2 fold change values in the vicinity of either +40 or -40. Downregulated genes were primarily involved in transcription/translation processes, whereas encoded upregulated genes governed amino acid and nucleotide biosynthesis pathways. The BioCyc metabolic pathways associated with infection also revealed that downregulated genes were mainly involved in fly immunity processes. Importantly, our study demonstrates that data on the molecular cross-talk between the host and the parasite (as well as the always present fly microbiome) recorded from an experimental biological model has a counterpart in field flies, which in turn validates the use of experimental host/parasite couples.

Keywords: field tsetse fly, Cameroonian foci, RNAseq, vector control, trypanosomiasis

<footnote n="1">¹Faculty of Science, University of Yaoundé I, Yaoundé, Cameroon, ²UMR 177, IRD-CIRAD, CIRAD TA A-17/G, Campus International de Baillarguet, Montpellier, France, ³Aix-Marseille University, INSERM, TAGC, Technological Advances for Genomics and Clinics, UMR S 1090, Marseille, France</footnote>
# INTRODUCTION

Human African trypanosomiasis [HAT or sleeping sickness; (1)] and animal African trypanosomiasis [AAT or nagana; (2)] are two vector-borne diseases that inflict heavy social and economic burdens on sub-Saharan African populations. Although the number of newly diagnosed HAT cases is decreasing (<10,000 per year) (3), more than 60 million people living in endemic areas are at risk of infection (4). In addition, AAT causes a large amount of livestock loss, which has been estimated as high as US$ 4.5 billion per year (5). HAT is due to either *Trypanosoma brucei gambiense* (Tbg; the chronic form of the disease in West and Central Africa) or *T. b. rhodesiense* (the acute form of the disease in East Africa), which are, respectively, transmitted by *Glossina palpalis* and *G. morsitans*. In contrast, AAT is caused by *T. b. brucei*, *T. congolense* (Tc; the forest or savannah type), or *T. vivax*, and is transmitted by *G. palpalis* or *G. morsitans*.

Despite differences between Tc and Tbg [reviewed in Ref. (6)], the parasites share several important characteristics. In particular, they are digenetic, meaning that they need to successively infect two different hosts to achieve their life cycle. One of these hosts, a *Glossina* fly, is the vector, whereas the other host is a vertebrate, typically a mammal. These parasites must accomplish a crucial part of their life cycle within their vector, namely their multiplication and maturation into the infectious form that can be transmitted to the vertebrate host while the tsetse fly ingests its blood meal. More specifically, Tc and Tbg undergo sequential differentiations after their ingestion by the fly, from the ingested blood stream form to the vertebrate-infective metacyclic form. The latter differentiation occurs either in the proboscis (for Tc) or in the salivary glands (for Tbg) (7), which is the basis for their respective classification into two different subgenera, *Nannomonas* and *Trypanozoon* (8). They also share the ability to excrete/secrete a number of proteins, some of which are considered to be involved in their establishment in the tsetse midgut and/or in the pathogenic process developed within the vertebrate host (9–12). Finally, the establishment of both Tc and Tbg in the *G. palpalis* vector is reported to be favored by *Sodalis glossinidius*, the secondary symbiont inhabitant of the tsetse gut (13). This finding demonstrates the occurrence in naturally infected field tsetse flies of a tripartite interaction (fly/trypanosome/gut bacteria) already reported to occur in experimentally infected insectary flies (14–18).

Another similarity between the parasites is that their mantle, which consists of a variant surface glycoprotein, allows them to evade the host's immune system by means of antigenic variations (19–21), thus rendering ineffective any vaccine approaches to fight HAT or AAT. Nevertheless, progress has been made in rapid diagnosis (22) and therapy that uses a nifurtimox-elfornitine combination in the treatment of the second phase of HAT (23). Furthermore, besides the use of trypanocidal drugs, the incidence of AAT can be lowered by introducing trypanotolerant cattle into AAT-infected area or through the antibody-mediated inhibition of trypanosome-secreted proteins involved in the parasite pathogenic process (24, 25).

Another approach to fight HAT or AAT is by vector control. Diverse strategies are available, including the application of

pesticides, the use of sterile males, and the development of para-transgenic approaches (26–32).

The normal status of tsetse flies is considered to be refractory to trypanosome infection, given that artificial or natural infection rates are always low (28, 33–36). Recently, a global transcriptomic analysis was performed (15–17) in the context of an anti-vector strategy, aimed at deciphering the molecular cross-talk occurring between the different participants involved in tsetse infection: the fly, the trypanosome, and the fly gut bacteria, especially the primary (*Wigglesworthia glossinidia*) and secondary (*S. glossinidius*) symbionts. The authors also focused on identifying differentially expressed genes (DEGs) associated with fly susceptibility or refractoriness as a result of fly infection by the trypanosome. These investigations were performed on insectary-raised *G. p. gambiensis* (Gpg) flies that were artificially infected (or not) by Tbg. This study raised the question of whether the results recorded under these experimental conditions could be transposed to what actually occurs under natural conditions in HAT and AAT foci.

To address this question we have conducted similar transcriptomic analyses on *G. p. palpalis* (Gpp) flies infected or not with Tc, collected in two HAT foci in southern Cameroon. Our experimental design involved a different host vector/parasite couple (Gpp/Tc) from what was used in the previous insectary-raised approach. However, as shown above and in support of this approach, several notable characteristics are shared between the Gpp/Tc couple and the previously used Gpg/Tbg couple, including life cycle, sequential differentiation within the vector, transmission modalities, host immune response escape, and pathological effects on susceptible vertebrate hosts, among others. Thus, the objectives of this study were to determine whether or not field-collected tsetse flies react to trypanosome infection under natural conditions similar to insectary flies under experimental conditions, and whether or not Tc induces molecular disruptions in Gpp similar to those provoked by Tbg in Gpg. Importantly, our approach provides novel evidence that validates the use of experimental host/parasite couples in the context of investigating anti-vector strategies.

# MATERIALS AND METHODS

## Sampling Areas

Tsetse flies were sampled in May and June 2015 in two active HAT foci (Campo and Bipindi), located in the Ocean Division of the southern region of Cameroon. The Campo focus (2°20′N, 9°52′E) is located on the Atlantic coast and extends along the Ntem river. The HAT National Control Program that visits Campo once per year diagnosed 61 novel HAT cases between 2001 and 2011. The passive identification of two cases in 2012 (37) indicates that HAT is still present. The Bipindi focus (3°2′N, 10°22′E) has a typical forest bioecological environment, including equatorial forest and farmland along roads and around villages. This focus has been recognized since 1920 (38) and includes several villages. Sleeping sickness is still present, since approximately 83 HAT cases were identified by the National Control Program in this focus between 1998 and 2011 (Ebo'o Eyenga, personal communication). In addition to HAT cases that involve *G. palpalis* and Tbg, regular global
surveys have identified the presence of several other Glossina (including Gpp) and *Trypanosoma* species (including Tc) in both foci. Surveys have also identified a variety of domestic and wild animals that serve as reservoirs for diverse *Trypanosoma* species (39–42). As described below, flies were trapped in these areas in order to select non-infected and Tc-infected individuals.

## Fly Sampling, Dissection, and Subsequent RNA Preservation

The May 2015 tsetse fly trapping campaign was conducted in three Campo villages (Ipono, Mabiogo, and Campo-Beach), and the June 2015 campaign was conducted in three Bipindi villages (Lambi, Bidjouka, and Ebiminbang). The geographical positions of the sampling sites were determined by GPS. Tsetse flies were captured using pyramidal traps (43) placed in suitable tsetse fly biotopes. Each trap was installed for four consecutive days, and the flies were collected twice per day.

Prior to handling samples, work stations and dissecting instruments were cleaned with RNase away (Ambion) in order to eliminate any RNases that could degrade sample RNA. Furthermore, tsetse flies were dissected alive to prevent RNA degradation by normal post mortem degradation processes. The first step in sample processing consisted in identifying the collected tsetse flies to the species level on the basis of morphological criteria and adapted taxonomic keys (44). Next, the samples were separated into two groups of teneral and non-teneral flies. The non-teneral Gpp flies were dissected in a drop of sterile 0.9% saline solution, according to the midgut dissection protocol developed by Penchenier and Itard (45). The organs were immediately transferred to tubes containing RNAlater (Ambion) for DNA and RNA extraction. These samples were then used for parasite identification by specific PCR amplification, and ultimately for transcriptomic analysis. All tools were carefully cleaned after the dissection of each fly to prevent cross-contamination. During field manipulations, the tubes containing the organs were stored at -20°C for 5 days; subsequently, they were stored in the laboratory at -80°C until use.

## DNA and RNA Extraction

To prepare for extraction, samples stored at -80°C were thawed and RNAlater was removed. The midguts were treated with the NucleoSpin TriPrep extraction kit (Macherey-Nagel) according to the manufacturer's instructions, which allow the separate extraction of DNA and RNA. RNase free water (40 µl) was added to elute the RNA, and 100 µl of DNA elute solution was added to recover the DNA. RNA quality and the absence of any DNA

contamination were checked on an Agilent RNA 6000 Bioanalyzer and quantified using the Agilent RNA 6000 Nano kit (Agilent Technologies, France).

## PCR Amplification

To identify which trypanosome species had infected the sampled tsetse flies, the isolated DNA samples stored at -80°C were thawed and used as a template for PCR amplification with specific primers (Table 1). PCR amplification of parasites was performed as described by Herder et al. (46) and consisted of a denaturing step at 94°C (5 min) followed by 44 amplification cycles, each comprising a denaturing step at 94°C (30 s), annealing at 55°C (30 s), and an extension step at 72°C (1 min). A final extension was performed at 72°C for 10 min. The amplified products were separated on a 2% agarose gel containing ethidium bromide and visualized under UV illumination. Positive (2 ng of reference DNA) and negative controls were included in each PCR amplification experiment. PCR amplifications that gave a positive result were repeated once for confirmation.

## RNA-Seq Processing

### Preparation of cDNA Libraries

Total RNA from 10 Gpp flies (5 non-infected flies and 5 flies infected by Tc s.l.) was assayed using the TruSeq mRNA-seq Stranded v2 Kit (Illumina), according to the manufacturer's instructions. Briefly, 4 µg of total RNA were used for poly(A)-selection to generate 120-210 bp cDNA fragments (mean size: 155 bp) after an 8-min elution-fragmentation incubation. Each library was barcoded using TruSeq Single Index (Illumina), according to the manufacturer's instructions. After library preparation, Agencourt AMPure XP (Beckman Coulter, Inc.) was used to select 200- to 400-bp size libraries. Each library size distribution was examined using the Bioanalyzer with a High Sensitivity DNA chip (Agilent) to ensure that the samples had the proper size and that they were devoid of any adaptor contamination. The sample concentration was quantified on Qubit with the Qubit® dsDNA HS Assay Kit (Life Technologies). Each library was then diluted to 4 nM and pooled at an equimolar ratio.

## NextSeq-500 Sequencing

For sequencing, 5 µL of pooled libraries (4 nM) were denatured with 5 µl NaOH (0.2 N) according to the manufacturer's instructions. Following a 5-min incubation, 5 µl of Tris-HCl (200 mM; pH 7) were added, and 20 pM of the pooled libraries were diluted with HT1 to a 1.6-pM final concentration. As a sequencing

<table><caption>**TABLE 1** | Primers used for PCR amplification of trypanosomes.</caption><thead><tr><td>Species</td><td>Primer sequence</td><td>Amplified product (bp)</td><td>Reference</td></tr></thead><tbody><tr><td>T. brucei s.l.</td><td>5'-CGAATGAATATTAAACAATGCGCAG-3'<br/>5'-AGAACCATTATTAGCTTTGTTGC-3'</td><td>164</td><td>Masiga et al. (47)</td></tr><tr><td>T. congolense ("forest" type)</td><td>5'-CGAATGAATATTAAACAATGCGCAG-3'<br/>5'-AGAACCATTATTAGCTTTGTTGC-3'</td><td>350</td><td>Masiga et al. (47)</td></tr><tr><td>T. congolense ("savannah" type)</td><td>5'-CGAATGAATATTAAACAATGCGCAG-3'<br/>5'-AGAACCATTATTAGCTTTGTTGC-3'</td><td>341</td><td>Moser et al. (48)</td></tr></tbody></table>

T. *Trypanosoma*; s.l., sensu lato.
control, a PhiX library was denatured and diluted according to the manufacturer's instructions, and 1.2 µl were added to the sample of denatured and diluted pooled libraries before loading. Finally, the libraries were sequenced on a high-output flow cell (400M clusters) using the NextSeq® 500/550 High Output Kit v2 (150 cycles; Illumina) in paired-end 75/75 nt mode, according to the manufacturer's instructions.

Datasets for the reads are available from the NCBI, GEO submission, accession number GSE98989.

# Bioinformatics Analysis

## Workflow

The successive tasks of the bioinformatics analysis were managed using a Snakemake workflow (49). This workflow enables the reproduction of all analyses from the raw read files and is available from the supporting Web site.<fnref n="1" />

## Reference Genomes

The *G. palpalis* genomic sequences (Glossina-palpalis-IAEA_SCAFFOLDS_GpapI1.fa) and annotations (Glossina-palpalis-IAEA_BASEFEATURES_GpapI1.1.gff3) were downloaded from VectorBase (50). For the annotation of Tc genes, the reference genome (TriTrypDB-9.0_TcongolenseIL3000.gff) was downloaded from TriTrypDB (51), whereas the *Drosophila melanogaster* reference genome (*Drosophila_melanogaster*.BDGP6.30.gff3) was downloaded from Flybase (52).

## Lane Merging

Since each sample was sequenced on four lanes, the original fastq-formatted read files were merged to produce two files per sample (one for each paired-end extremity).

## Read Quality Control

FastQC<fnref n="2" /> was run on the raw reads in order to check their quality.

## Read Mapping

Raw reads were mapped onto the genome with the local alignment algorithm Subread-align (53) in paired-ends mode with at most 10 mismatches. Read mapping statistics were computed using samtools flagstats (54) and are summarized in Table S1 in Supplementary Material.

## Read Counts per Gene

The number of read pairs (fragments) per gene was counted using the featureCounts tool from the Subread package (55), including the option “feature type” to only count reads overlapping transcripts.

## Visualization

Genome maps were generated using the Integrative Genomic Viewer (56).

## Detection of DEGs

Differential gene expression analysis was performed using the SARTools R package (57), which separately runs DESeq2 (58) and egdeR (59) as well as generates readable reports.

## Identification of Orthologs between Gpp and D. melanogaster

Because the *G. palpalis* genome is inadequately assembled and annotated, our functional interpretation of the DEGs relied on a comparative genomics approach. This was based on the identification of bidirectional best hits (BBH) between all sequences of *G. palpalis* and *D. melanogaster* (assembly BDGP6). BBH were identified using blastp (60) and the BLOSUM45 substitution matrix, and by setting a threshold of 10⁻⁵ on the expected score.

## Functional Enrichment of DEGs

Identification of functions associated with the DEGs was based on *Drosophila* orthologs of the DEGs (ortho-DEG). Functional enrichment was separately performed using the DAVID (61) and g:Profiler (62) tools. The Bonferroni correction was used to obtain the enrichments of these functions, with a threshold set at 10⁻³.

## Pathway Mapping of DEGs

*Drosophila* orthologs of the DEGs were loaded into the metabolic cellular overview of BioCyc (63) in order to highlight the pathways affected by the infection.

## Statistical Treatment of Entomological Data

Entomological data, as well as all other calculations, were evaluated using the statistical package SPSS version 2.0. Spreadsheets were made using Microsoft Office Excel 2007.

# RESULTS

## Entomological Data

A total of 1,991 tsetse flies were collected during the entomological survey (775 flies from Campo and 1,216 flies from Bipindi). The Campo fly population was composed of Gpp (95.61%), *Glossina caliginea* (2.06%), *Glossina palicera* (1.87%), and *Glossina nigrofusca* (0.52%). Two tsetse fly species were identified at the Bipindi focus: Gpp (99.34%) and *G. palicera* (0.66%). The mean apparent density was 4.24 flies per trap per day; however, this parameter was highly variable between the different villages and was higher in Bipindi (8.1) than in Campo (3.52) (Table 2). The frequency of teneral flies was typically low in both Bipindi (0.08%) and Campo (1.16%). These data are roughly in line with data recorded in 2007/2008 (13), although the rate of teneral flies was much lower in the present study. Only 1,245 of the trapped 1,991 tsetse flies were dissected, since 10 flies were teneral and 736 flies were desiccated.

## PCR Identification of Trypanosome Species in the Tsetse Midgut

The number of flies carrying single or mixed trypanosome infections is presented in Table 3. Of the 337 Campo flies analyzed, 25 (7.41%) were infected by the Tc “forest” type, 16 (4.74%) by the

<footnote n="1">¹https://github.com/rioualen/gene-regulation.</footnote>

<footnote n="2">²http://www.bioinformatics.babraham.ac.uk/projects/fastqc/.</footnote>
<table><caption>**TABLE 2** | Entomological field data from the Bipindi and Campo foci.</caption><thead><tr><th>Focus</th><th>Village</th><th>Number of traps</th><th>Number of tsetse flies captured</th><th>ADT</th><th>Number of teneral tsetse flies (%)</th><th>Number of tsetse flies dissected</th></tr></thead><tbody><tr><td rowspan="4">Campo</td><td>Ipono</td><td>15</td><td>161</td><td>2.68</td><td>4 (2.48)</td><td>110</td></tr><tr><td>Beach</td><td>15</td><td>341</td><td>4.55</td><td>3 (0.87)</td><td>264</td></tr><tr><td>Mabiogo</td><td>18</td><td>273</td><td>3.21</td><td>2 (0.73)</td><td>228</td></tr><tr><td>Total Campo</td><td>48</td><td>775</td><td>3.52</td><td>9 (1.16)</td><td>602</td></tr><tr><td rowspan="4">Bipindi</td><td>Bidjouka</td><td>23</td><td>608</td><td>5.28</td><td>1 (0.16)</td><td>278</td></tr><tr><td>Lambi</td><td>12</td><td>486</td><td>8.1</td><td>0 (0)</td><td>303</td></tr><tr><td>Ebimimbang</td><td>15</td><td>122</td><td>1.74</td><td>0 (0)</td><td>72</td></tr><tr><td>Total Bipindi</td><td>50</td><td>1,216</td><td>4.96</td><td>1 (0.50)</td><td>653</td></tr><tr><td>Total</td><td></td><td>98</td><td>1,991</td><td>4.24</td><td>10 (0.50)</td><td>1,255</td></tr></tbody></table>

ADT, apparent density per trap per day.

<table><caption>**TABLE 3** | Number of *Trypanosoma congolense* s.l. simple and mixed infections by village.</caption><thead><tr><th>Focus</th><th>Village</th><th>Number of tsetse flies analyzed</th><th>Number of flies infected with TcF</th><th>Number of flies infected with TcS</th><th>Number of flies carrying a mixed infection</th></tr></thead><tbody><tr><td rowspan="4">Campo</td><td>Ipono</td><td>63</td><td>2</td><td>1</td><td>1</td></tr><tr><td>Beach</td><td>170</td><td>15</td><td>8</td><td>7</td></tr><tr><td>Mabiogo</td><td>104</td><td>8</td><td>7</td><td>6</td></tr><tr><td>Total Campo</td><td>337</td><td>25</td><td>16</td><td>14</td></tr><tr><td rowspan="4">Bipindi</td><td>Bidjouka</td><td>40</td><td>1</td><td>0</td><td>0</td></tr><tr><td>Lambi</td><td>33</td><td>5</td><td>0</td><td>0</td></tr><tr><td>Ebimimbang</td><td>11</td><td>1</td><td>0</td><td>0</td></tr><tr><td>Total Bipindi</td><td>84</td><td>7</td><td>0</td><td>0</td></tr><tr><td>Total</td><td></td><td>421</td><td>32</td><td>16</td><td>14</td></tr></tbody></table>

TcF, T. congolense "forest" type; TcS, T. congolense "savannah" type.

Tc "savanah" type, and 14 (4.15%) by both parasites. In contrast, Bipindi flies only carried the Tc "forest" type (8.33%). Table S1 in Supplementary Material details the characteristics of the different samples including those used for transcriptomic analyses.

## Raw Data

The sequencing of libraries produced a total of 400 million reads (theoretically 40 million reads per sample), which represent a satisfactory sequencing depth for subsequent differential gene expression analysis (Figure 1). Out of the 328 raw clusters generated, 77.6% were successfully filtered, with each sample producing 50–72 million clusters (mean: 61 million clusters) (Figure 1). Sequencing also revealed a total of 31,320 contigs distributed in 3,926 scaffolds with a mean size of 96,817 bp (varying in size from 545 bp to 3.6 Mb).

## Mapping on Gpp

RNA-seq sequencing produced an average of 124 million reads per sample. From this, 111.45 million reads (83.6%) were mapped onto the genome of Gpp, 103.68 M (73.7%) of which were properly paired (Table 4).

## Mapping on Tc

Reads were also mapped onto the trypanosome genome in order to validate the infection status of the samples, as well as to

investigate the role played by the trypanosome in this molecular dialog. Since trypanosome cells represented a small fraction of the analyzed material, only a smaller fraction of the reads could be mapped. Specifically, 300,302.5 (0.24%) reads were mapped from an average of 124 million reads per sample. This resulted in 137,305.4 (0.11%) properly paired reads and 157,712.5 (0.12%) singletons (Table 4).

## DEGs between Infected and Non-Infected Flies

We used DESeq2 (Table 5) to detect genes that were differentially expressed between the five infected and five non-infected samples. When a Benjamini–Hochberg corrected *p*-value lower than 0.05 was applied, 524 genes were observed to be significantly differentially expressed in infected vs. non-infected flies, among which 285 genes were downregulated and 239 were upregulated. A similar DEG analysis was performed using edgeR (Table 5), which identified only 20 downregulated genes and 53 upregulated genes. Figure 2 presents the volcano plots produced by DESeq2 (Figure 2A) and edgeR (Figure 2B); genes that were significantly (*p*-value <0.05) differentially expressed and with a fold change of log2 (fold change) >2 (i.e., upregulated genes) or log2 (fold change) <-2 (i.e., downregulated genes) were considered relevant.

## Functional Annotation

To understand the roles of DEGs associated with tsetse fly infection by Tc, identifiers of the *Drosophila* orthologs of the *Glossina* DEGs were examined using the DAVID-functional enrichment tool. Separately, we analyzed the 290 downregulated and 213 upregulated genes reported by DESeq2 and identified 207 *Drosophila* best hits (121 downregulated and 86 upregulated genes). The same analysis by edgeR only resulted in 25 DEG orthologs (6 downregulated and 19 upregulated genes). This reduced number of orthologs could possibly be due to the incomplete assembly and annotation of the *Glossina* genome and/or the high stringency of the BBH criterion (which discards the case where several tsetse fly proteins have the same closest hit in *Drosophila*). These DEGs were examined using DAVID, which compares the list of input genes with a variety of functional annotations. This analysis was focused on the three primary categories of the Gene
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | The score quality of different sequences of clusters, presented base by base. Green indicates a good quality. (A) Total read count per sample (millions). (B) Number of clusters per sample.</figcaption></figure>

<table><caption>**TABLE 4** | Read mapping statistics.</caption><thead><tr><th>Mapping on</th><th>Mapped reads (%)</th><th>Properly paired reads (%)</th><th>Singleton reads (%)</th><th>QC-passed reads</th></tr></thead><tbody><tr><td>Glossina</td><td>111.45 M (89.87)</td><td>103.68 M (83.61)</td><td>3.62 M (2.92)</td><td>123.8 M</td></tr><tr><td>Trypanosoma</td><td>0.3 M (0.24)</td><td>0.13 M (0.11)</td><td>0.15 M (0.12)</td><td>123.8 M</td></tr></tbody></table>

M, millions of reads.

<table><caption>**TABLE 5** | Differentially expressed genes.</caption><thead><tr><th></th><th>Downregulated</th><th>Upregulated</th><th>Total</th></tr></thead><tbody><tr><td>DESeq2</td><td>285</td><td>239</td><td>524</td></tr><tr><td>edgeR</td><td>20</td><td>53</td><td>73</td></tr><tr><td>Common</td><td>20</td><td>46</td><td>66</td></tr><tr><td>Total</td><td>285</td><td>246</td><td>531</td></tr></tbody></table>

Ontology annotation: biological process (BP), molecular function (MF), and cellular component (CC). The list of different features is provided in **Table 6**.

The 285 downregulated genes identified by DESeq2 in flies infected with Tc mainly belonged to the BP category, in which the major functional classes were RNA processing (58 genes; 23.48%), ribosome biogenesis (24 genes; 9.71%), and translation (13 genes; 5.26%); the other genes corresponded to several poorly represented classes. The MF category included the functional RNA binding classes (20 genes; 8.09%) and catalytic activity (25 genes; 10.12%). Finally, the CC category included the intracellular lumen (48 genes; 19.43%), non-membrane-bound organelle (34 genes; 13.75%), and ribonucleoprotein complex (25 genes; 10.12%) functional classes (Figure 3).

In addition, 239 DEGs were overexpressed in Tc-infected tsetse flies. These DEGs encoded proteins corresponding to the same three primary ontology categories. The BP category includes

the neuron morphogenesis functional class (83 genes; 38.96%), amino acid biosynthesis (24 genes; 11.26%), and carboxylic acid biosynthesis (10 genes; 4.7%). The MF category included the iron binding (64 genes, 30%) and catalytic activity (24 genes; 11.26%) functional classes. Finally, the CC category only included the mitochondrion functional class (8 genes; 3.75%).

As already shown for DEGs, edgeR provides a much lower number of functional annotations than DESeq2. Here, using edgeR, fly genes that displayed an increased expression in response to Tc infection were found to belong to the MF category, with only three functional classes: phospholipase activity (2 genes; 10.5%), lipase activity (2 genes; 10.5%), and carboxylesterase activity (2 genes; 10.5%). Finally, expression was decreased for only two genes belonging to the one-carbon metabolic process term (BP category).

## Functional Enrichment of DEGs

To make our analysis more focused and efficient, a functional enrichment was performed to refine the list of tsetse fly DEGs in which expression was influenced by Tc infection. We therefore combined fold enrichment and the *p*-value at a 5% threshold, which allowed applying a Bonferroni correction to eliminate false positives. The Bonferroni correction threshold was fixed at $\alpha = 10^{-2}$, and all functionality with a Bonferroni value below this threshold was considered to be due to trypanosome infection.

Following this correction, 16 functional classes were found to be selectively altered by trypanosome infection. These classes are mainly involved in the transcription process, including (a) RNA related processes (rRNA processing, rRNA metabolic process, ncRNA processing, ncRNA metabolic processes, and RNA processing), involving 63 DEGs; (b) monitoring processes related to the synthesis of the ribonucleoprotein complex (ribonucleoprotein complex biogenesis and ribonucleoprotein
<figure><img src="image_4.png" /></figure>

**FIGURE 2** | Comparisons of significantly differentially expressed genes. DESeq2 (A) and edgeR (B) results are illustrated by volcano plots, in which the differentially expressed features are shown in red. Upregulated genes are thus observed as positive values, and downregulated genes as negative values. Triangles correspond to features where the log of the adjusted *p*-value is too low or too high to be displayed on the plot.

complex), involving 29 DEGs; (c) RNA binding, involving 24 DEGs; (d) nucleolus biogenesis (nucleolus and nuclear lumen), involving 23 DEGs; (e) ribosome synthesis, involving 17 DEGs; and (f) eukaryotic translation factor 3 complex, involving 4 DEGs. In contrast to the transcription process, which involved 160 DEGs, the BPs of organic acid synthesis (amine biosynthetic process, cellular amino acid biosynthetic process, carboxylic acid biosynthetic process, and organic acid biosynthetic process) that were found to be activated by trypanosome infection only involved 21 DEGs (Table 7).

## Associated Metabolic Pathways

The BioCyc metabolic map (Figure 4) and Table 8 illustrate the different pathways that the DEGs are involved in. Among these, the amino acid biosynthesis pathway (which includes the biosynthesis of L-glutamine, L-glutamate, L-serine, L-asparagine, L-aspartate, etc.) is controlled by genes that were shown to be overexpressed following trypanosome infection in the flies. Similarly, DEGs associated with the nucleotide biosynthesis pathway were overexpressed following trypanosome infection, especially uridine monophosphate, an RNA monomer. In contrast, genes involved in the pentose phosphate pathway, namely those implicated in the synthesis of D-ribose 5-phosphate, were downregulated. Phosphorylated pentose is converted by ribose phosphate diphosphokinase into phosphoribosylpyrophosphate, a precursor of nucleotide synthesis. Finally, regarding the carbohydrate biosynthesis pathway, we observed an overexpression of genes encoding malate dehydrogenase, which converts malate into pyruvate.

Interestingly, a large number of up- and downregulated DEGs were related to a given biosynthetic process meaning that

upregulated DEGs encoding amino acids (such as tyrosine, serine, glutamine, and several others) were found in the amino acid synthesis pathway. Other overexpressed genes that were identified are involved in the biosynthesis of galactosyltransferase, N-acetylglucosamine, and beta-1,4-manosylglycolipide, which are all molecules that interact with the immune system of the tsetse fly (64–66). In contrast, genes involved in the biosynthesis of cytokines were downregulated in trypanosome-infected flies, as compared to non-infected flies. Similarly, genes involved in folate metabolism (e.g., the biosynthesis of formyltetrahydrofolate dehydrogenase), the main source of energy in flies, were downregulated.

Finally, DEGs involved in the transport of several molecules from the extracellular space toward the cytosol compartment were upregulated. This transport includes molecules with a role in cell nutrition, and nutrients such as lipids, but also ATP, succinate, and L-carnitine (which participates in the degradation of fats).

## DISCUSSION

Understanding the mechanisms involved in tsetse fly susceptibility or refractoriness to trypanosome infection is crucial for developing a novel anti-vector based strategy to control the spread of sleeping sickness and nagana. One recent study was performed within this context to identify genes in Gpg associated with its susceptibility or refractoriness to Tbg infection, using an RNA-seq approach (17). The underlying hypothesis was that some of the genes involved in controlling fly susceptibility/resistance to trypanosome infection could be targeted in order to increase the refractoriness of the fly, thereby decreasing its vector competence while enabling the development of an anti-vector strategy against
<table><caption>TABLE 6 | Functional annotation of differentially expressed genes (DEGs).</caption><thead><tr><th>Category</th><th>Term</th><th>Number of DEG</th><th>p-value</th><th>Log2 (fold enrichment)</th></tr></thead><tbody><tr><td colspan="5"><strong>DEGs identified with DESeq2</strong></td></tr><tr><td>GOTERM_BP_FAT</td><td>Amine biosynthetic process</td><td>6</td><td>2.80E-05</td><td>15.70</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>Amino acid biosynthesis</td><td>3</td><td>2.00E-03</td><td>42.80</td></tr><tr><td>GOTERM_MF_FAT</td><td>ATPase activity, uncoupled</td><td>5</td><td>9.10E-02</td><td>-2.90</td></tr><tr><td>GOTERM_MF_FAT</td><td>ATP-dependent helicase activity</td><td>4</td><td>5.80E-02</td><td>-4.50</td></tr><tr><td>GOTERM_MF_FAT</td><td>ATP-dependent RNA helicase activity</td><td>4</td><td>1.10E-02</td><td>-8.40</td></tr><tr><td>GOTERM_BP_FAT</td><td>Axon guidance</td><td>6</td><td>2.90E-03</td><td>5.90</td></tr><tr><td>GOTERM_BP_FAT</td><td>Axonal defasciculation</td><td>2</td><td>8.70E-02</td><td>21.50</td></tr><tr><td>GOTERM_BP_FAT</td><td>Axonogenesis</td><td>6</td><td>1.20E-02</td><td>4.20</td></tr><tr><td>GOTERM_BP_FAT</td><td>Carboxylic acid biosynthetic process</td><td>5</td><td>9.90E-04</td><td>10.80</td></tr><tr><td>GOTERM_MF_FAT</td><td>Cation binding</td><td>18</td><td>6.10E-02</td><td>1.50</td></tr><tr><td>GOTERM_BP_FAT</td><td>Cell morphogenesis</td><td>7</td><td>8.10E-02</td><td>2.20</td></tr><tr><td>GOTERM_BP_FAT</td><td>Cell morphogenesis involved in differentiation</td><td>6</td><td>5.70E-02</td><td>2.80</td></tr><tr><td>GOTERM_BP_FAT</td><td>Cell morphogenesis involved in neuron differentiation</td><td>6</td><td>5.00E-02</td><td>2.90</td></tr><tr><td>GOTERM_BP_FAT</td><td>Cell motion</td><td>6</td><td>5.20E-02</td><td>2.90</td></tr><tr><td>GOTERM_BP_FAT</td><td>Cell part morphogenesis</td><td>6</td><td>7.60E-02</td><td>2.60</td></tr><tr><td>GOTERM_BP_FAT</td><td>Cell projection morphogenesis</td><td>6</td><td>7.20E-02</td><td>2.60</td></tr><tr><td>GOTERM_BP_FAT</td><td>Cell recognition</td><td>3</td><td>4.40E-02</td><td>8.70</td></tr><tr><td>GOTERM_BP_FAT</td><td>Cellular amino acid biosynthetic process</td><td>5</td><td>5.40E-05</td><td>22.40</td></tr><tr><td>GOTERM_BP_FAT</td><td>Chemical homeostasis</td><td>3</td><td>7.70E-02</td><td>6.30</td></tr><tr><td>GOTERM_MF_FAT</td><td>Coenzyme binding</td><td>4</td><td>6.40E-02</td><td>4.30</td></tr><tr><td>GOTERM_CC_FAT</td><td>Cytosol</td><td>6</td><td>9.30E-02</td><td>-2.40</td></tr><tr><td>INTERPRO</td><td>DEAD-like helicase, N-terminal</td><td>4</td><td>6.20E-02</td><td>-4.30</td></tr><tr><td>GOTERM_BP_FAT</td><td>Defasciculation of motor neuron axon</td><td>2</td><td>6.20E-02</td><td>30.70</td></tr><tr><td>SMART</td><td>DEXDc</td><td>4</td><td>4.30E-02</td><td>-4.90</td></tr><tr><td>GOTERM_BP_FAT</td><td>Di-, tri-valent inorganic cation transport</td><td>3</td><td>3.40E-02</td><td>10.10</td></tr><tr><td>INTERPRO</td><td>DNA/RNA helicase, C-terminal</td><td>4</td><td>6.20E-02</td><td>-4.30</td></tr><tr><td>INTERPRO</td><td>DNA/RNA helicase, DEAD/DEAH box type, N-terminal</td><td>4</td><td>1.60E-02</td><td>-7.40</td></tr><tr><td>GOTERM_MF_FAT</td><td>Electron carrier activity</td><td>5</td><td>2.50E-02</td><td>4.40</td></tr><tr><td>GOTERM_MF_FAT</td><td>Enzyme inhibitor activity</td><td>3</td><td>9.90E-02</td><td>5.50</td></tr><tr><td>GOTERM_CC_FAT</td><td>Eukaryotic translation initiation factor 3 complex</td><td>4</td><td>4.40E-04</td><td>-24.50</td></tr><tr><td>GOTERM_MF_FAT</td><td>Glutamate synthase activity</td><td>2</td><td>2.10E-02</td><td>93.60</td></tr><tr><td>GOTERM_BP_FAT</td><td>Glutamine family amino acid biosynthetic process</td><td>2</td><td>9.50E-02</td><td>19.60</td></tr><tr><td>GOTERM_BP_FAT</td><td>Glutamine family amino acid metabolic process</td><td>3</td><td>1.50E-02</td><td>15.40</td></tr><tr><td>GOTERM_BP_FAT</td><td>Glutamine metabolic process</td><td>2</td><td>6.20E-02</td><td>30.70</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>Helicase</td><td>4</td><td>6.80E-02</td><td>-4.20</td></tr><tr><td>INTERPRO</td><td>Helicase, superfamily 1 and 2, ATP-binding</td><td>4</td><td>6.00E-02</td><td>-4.40</td></tr><tr><td>SMART</td><td>HELICc</td><td>4</td><td>4.30E-02</td><td>-4.90</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>Heme</td><td>3</td><td>8.40E-02</td><td>6.10</td></tr><tr><td>GOTERM_MF_FAT</td><td>Heme binding</td><td>4</td><td>2.90E-02</td><td>5.90</td></tr><tr><td>GOTERM_BP_FAT</td><td>Homeostatic process</td><td>4</td><td>8.90E-02</td><td>3.70</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>Hydrolase</td><td>14</td><td>6.20E-02</td><td>1.70</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>Initiation factor</td><td>3</td><td>3.40E-02</td><td>-10.10</td></tr><tr><td>GOTERM_CC_FAT</td><td>Intracellular non-membrane-bounded organelle</td><td>17</td><td>3.60E-03</td><td>-2.00</td></tr><tr><td>GOTERM_CC_FAT</td><td>Intracellular organelle lumen</td><td>12</td><td>1.00E-02</td><td>-2.30</td></tr><tr><td>GOTERM_MF_FAT</td><td>Ion binding</td><td>18</td><td>6.40E-02</td><td>1.50</td></tr><tr><td>GOTERM_MF_FAT</td><td>Iron ion binding</td><td>6</td><td>1.10E-02</td><td>4.30</td></tr><tr><td>KEGG_PATHWAY</td><td>Limonene and pinene degradation</td><td>3</td><td>6.60E-02</td><td>6.60</td></tr><tr><td>GOTERM_MF_FAT</td><td>Lipase activity</td><td>3</td><td>9.30E-02</td><td>5.70</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>Lipid-binding</td><td>2</td><td>4.80E-02</td><td>39.90</td></tr><tr><td>GOTERM_CC_FAT</td><td>Membrane-enclosed lumen</td><td>12</td><td>1.20E-02</td><td>-2.20</td></tr><tr><td>GOTERM_MF_FAT</td><td>Metal ion binding</td><td>18</td><td>5.10E-02</td><td>1.50</td></tr><tr><td>INTERPRO</td><td>Mitochondrial substrate carrier</td><td>3</td><td>3.10E-02</td><td>10.60</td></tr><tr><td>INTERPRO</td><td>Mitochondrial substrate/solute carrier</td><td>3</td><td>3.30E-02</td><td>10.20</td></tr><tr><td>GOTERM_CC_FAT</td><td>Mitochondrion</td><td>8</td><td>6.00E-02</td><td>2.10</td></tr><tr><td>GOTERM_BP_FAT</td><td>Mitotic spindle elongation</td><td>4</td><td>6.50E-02</td><td>-4.20</td></tr><tr><td>GOTERM_MF_FAT</td><td>mRNA binding</td><td>6</td><td>3.70E-02</td><td>-3.20</td></tr><tr><td>GOTERM_BP_FAT</td><td>ncRNA metabolic process</td><td>12</td><td>3.10E-07</td><td>-7.50</td></tr><tr><td>GOTERM_BP_FAT</td><td>ncRNA processing</td><td>12</td><td>6.00E-09</td><td>-10.80</td></tr><tr><td>GOTERM_BP_FAT</td><td>Neuron development</td><td>6</td><td>9.30E-02</td><td>2.40</td></tr><tr><td>GOTERM_BP_FAT</td><td>Neuron projection development</td><td>6</td><td>5.00E-02</td><td>2.90</td></tr><tr><td>GOTERM_BP_FAT</td><td>Neuron projection morphogenesis</td><td>6</td><td>4.90E-02</td><td>2.90</td></tr><tr><td>GOTERM_BP_FAT</td><td>Neuron recognition</td><td>3</td><td>4.40E-02</td><td>8.70</td></tr><tr><td>GOTERM_BP_FAT</td><td>Nitrogen compound biosynthetic process</td><td>6</td><td>1.10E-02</td><td>4.20</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 6 | Continued</caption><thead><tr><th>Category</th><th>Term</th><th>Number of DEG</th><th>p-value</th><th>Log2 (fold enrichment)</th></tr></thead><tbody><tr><td>GOTERM_CC_FAT</td><td>Non-membrane-bounded organelle</td><td>17</td><td>3.60E-03</td><td>-2.00</td></tr><tr><td>GOTERM_CC_FAT</td><td>Nuclear lumen</td><td>12</td><td>4.10E-04</td><td>-3.30</td></tr><tr><td>GOTERM_CC_FAT</td><td>Nucleolus</td><td>11</td><td>2.20E-09</td><td>-13.50</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>Nucleus</td><td>18</td><td>1.90E-02</td><td>-1.80</td></tr><tr><td>GOTERM_CC_FAT</td><td>Organelle lumen</td><td>12</td><td>1.00E-02</td><td>-2.30</td></tr><tr><td>GOTERM_BP_FAT</td><td>Organic acid biosynthetic process</td><td>5</td><td>9.90E-04</td><td>10.80</td></tr><tr><td>GOTERM_BP_FAT</td><td>Oxidation reduction</td><td>8</td><td>3.10E-02</td><td>2.50</td></tr><tr><td>GOTERM_MF_FAT</td><td>Phospholipase activity</td><td>3</td><td>3.60E-02</td><td>9.70</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>Phosphoprotein</td><td>16</td><td>4.90E-02</td><td>-1.70</td></tr><tr><td>GOTERM_BP_FAT</td><td>Positive regulation of protein kinase cascade</td><td>2</td><td>8.50E-02</td><td>-22.30</td></tr><tr><td>GOTERM_CC_FAT</td><td>Preribosome</td><td>3</td><td>3.40E-03</td><td>-31.50</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>Protein biosynthesis</td><td>5</td><td>1.70E-02</td><td>-5.00</td></tr><tr><td>GOTERM_BP_FAT</td><td>Pseudouridine synthesis</td><td>2</td><td>9.90E-02</td><td>-19.10</td></tr><tr><td>GOTERM_MF_FAT</td><td>Purine NTP-dependent helicase activity</td><td>4</td><td>5.80E-02</td><td>-4.50</td></tr><tr><td>GOTERM_BP_FAT</td><td>Regulation of translational initiation</td><td>3</td><td>1.10E-02</td><td>-18.20</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>Ribonucleoprotein</td><td>6</td><td>1.30E-02</td><td>-4.20</td></tr><tr><td>GOTERM_CC_FAT</td><td>Ribonucleoprotein complex</td><td>13</td><td>6.50E-05</td><td>-3.70</td></tr><tr><td>GOTERM_BP_FAT</td><td>Ribonucleoprotein complex biogenesis</td><td>12</td><td>2.80E-09</td><td>-11.60</td></tr><tr><td>GOTERM_BP_FAT</td><td>Ribosome biogenesis</td><td>12</td><td>7.10E-11</td><td>-16.00</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>Ribosome biogenesis</td><td>5</td><td>1.20E-04</td><td>-18.60</td></tr><tr><td>GOTERM_MF_FAT</td><td>RNA binding</td><td>14</td><td>9.00E-05</td><td>-3.60</td></tr><tr><td>GOTERM_MF_FAT</td><td>RNA helicase activity</td><td>4</td><td>1.50E-02</td><td>-7.50</td></tr><tr><td>INTERPRO</td><td>RNA helicase, ATP-dependent, DEAD-box, conserved site</td><td>3</td><td>4.20E-02</td><td>-9.00</td></tr><tr><td>INTERPRO</td><td>RNA helicase, DEAD-box type, Q motif</td><td>4</td><td>6.50E-03</td><td>-10.20</td></tr><tr><td>GOTERM_BP_FAT</td><td>RNA modification</td><td>3</td><td>7.10E-02</td><td>-6.70</td></tr><tr><td>GOTERM_BP_FAT</td><td>RNA processing</td><td>14</td><td>1.40E-05</td><td>-4.30</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>RNA-binding</td><td>10</td><td>2.20E-05</td><td>-6.40</td></tr><tr><td>GOTERM_MF_FAT</td><td>RNA-dependent ATPase activity</td><td>4</td><td>1.10E-02</td><td>-8.40</td></tr><tr><td>GOTERM_BP_FAT</td><td>rRNA metabolic process</td><td>10</td><td>2.10E-10</td><td>-22.30</td></tr><tr><td>GOTERM_BP_FAT</td><td>rRNA modification</td><td>2</td><td>9.90E-02</td><td>-19.10</td></tr><tr><td>GOTERM_BP_FAT</td><td>rRNA processing</td><td>10</td><td>2.10E-10</td><td>-22.30</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>rRNA processing</td><td>5</td><td>1.50E-04</td><td>-17.70</td></tr><tr><td>COG_ONTOLOGY</td><td>Secondary metabolites biosynthesis, transport, and catabolism</td><td>3</td><td>3.40E-02</td><td>8.60</td></tr><tr><td>GOTERM_CC_FAT</td><td>Small nuclear ribonucleoprotein complex</td><td>4</td><td>1.80E-02</td><td>-6.80</td></tr><tr><td>GOTERM_CC_FAT</td><td>Small nucleolar ribonucleoprotein complex</td><td>4</td><td>4.20E-05</td><td>-49.00</td></tr><tr><td>GOTERM_CC_FAT</td><td>Small-subunit processome</td><td>3</td><td>1.70E-03</td><td>-44.10</td></tr><tr><td>GOTERM_BP_FAT</td><td>Spindle elongation</td><td>4</td><td>6.70E-02</td><td>-4.20</td></tr><tr><td>GOTERM_MF_FAT</td><td>Tetrapyrrole binding</td><td>4</td><td>2.90E-02</td><td>5.90</td></tr><tr><td>GOTERM_BP_FAT</td><td>Translation</td><td>9</td><td>2.10E-02</td><td>-2.60</td></tr><tr><td>GOTERM_BP_FAT</td><td>Translational initiation</td><td>4</td><td>1.80E-02</td><td>-7.00</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>Transport</td><td>7</td><td>5.60E-02</td><td>2.50</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>WD repeat</td><td>7</td><td>1.00E-02</td><td>-3.80</td></tr><tr><td>SMART</td><td>WD40</td><td>7</td><td>8.30E-03</td><td>-3.80</td></tr><tr><td>INTERPRO</td><td>WD40 repeat</td><td>7</td><td>1.80E-02</td><td>-3.30</td></tr><tr><td>INTERPRO</td><td>WD40 repeat, conserved site</td><td>5</td><td>5.30E-02</td><td>-3.50</td></tr><tr><td>INTERPRO</td><td>WD40 repeat, region</td><td>6</td><td>2.60E-02</td><td>-3.60</td></tr><tr><td>INTERPRO</td><td>WD40 repeat, subgroup</td><td>7</td><td>7.80E-03</td><td>-4.00</td></tr><tr><td>INTERPRO</td><td>WD40/YVTN repeat-like</td><td>8</td><td>5.90E-03</td><td>-3.60</td></tr><tr><td colspan="5"><strong>DEGs identified with DESeq2</strong></td></tr><tr><td>GOTERM_BP_FAT</td><td>One-carbon metabolic process</td><td>2</td><td>6.00E-02</td><td>-27.2</td></tr><tr><td>GOTERM_MF_FAT</td><td>Carboxylesterase activity</td><td>2</td><td>8.40E-02</td><td>20.4</td></tr><tr><td>GOTERM_MF_FAT</td><td>Lipase activity</td><td>2</td><td>7.90E-02</td><td>21.7</td></tr><tr><td>GOTERM_MF_FAT</td><td>Phospholipase activity</td><td>2</td><td>4.80E-02</td><td>36.6</td></tr></tbody></table>

Categories: GOTERM_BP_FAT, biological process; GOTERM_CC_FAT, cellular component; GOTERM_MF_FAT, molecular function.

SMART & INTERPRO, protein domains; SPIR_KEYWORD, protein information resource provided by SWISSPROT and UniProt.

CDG_ONTOLOGY, cluster orthology group.

Black fonts: downregulated DEGs; red fonts: upregulated DEGs.

the disease. As this analysis was performed with insectary flies artificially infected with trypanosomes, it was necessary to verify that similar molecular events occur in field flies naturally infected by trypanosome vs. non-infected flies. To accomplish this,

we have chosen the Gpp/Tc couple, whose prevalence (even in HAT foci) is often higher than observed with the Gpp/Tbg couple.

As in the previous study (17), we employed an RNA-seq approach. This provided satisfactory results regarding the
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Functional annotation of differentially expressed genes using DAVID. (A) Downregulated genes identified by DESeq2. (B) Upregulated genes identified by DESeq2. The x-axis indicates the number of genes enriched for the term, and the y-axis indicates the functional classes that were differentially expressed.</figcaption></figure>

<table><caption>**TABLE 7** | Bonferroni correction for differentially expressed genes (DEGs) enriched functionalities.</caption><thead><tr><th>Category</th><th>Term</th><th>Number of DEGs</th><th>Bonferroni</th></tr></thead><tbody><tr><td colspan="4" style="text-align: center; font-weight: bold;">Downregulated DESeq2</td></tr><tr><td>GOTERM_BP_FAT</td><td>Ribosome biogenesis</td><td>12</td><td>3.90E-08</td></tr><tr><td>GOTERM_BP_FAT</td><td>rRNA processing</td><td>10</td><td>5.80E-08</td></tr><tr><td>GOTERM_BP_FAT</td><td>rRNA metabolic process</td><td>10</td><td>5.80E-08</td></tr><tr><td>GOTERM_BP_FAT</td><td>Ribonucleoprotein complex biogenesis</td><td>12</td><td>5.10E-07</td></tr><tr><td>GOTERM_BP_FAT</td><td>ncRNA processing</td><td>12</td><td>8.30E-07</td></tr><tr><td>GOTERM_BP_FAT</td><td>ncRNA metabolic process</td><td>12</td><td>3.40E-05</td></tr><tr><td>GOTERM_BP_FAT</td><td>RNA processing</td><td>14</td><td>1.20E-03</td></tr><tr><td>GOTERM_CC_FAT</td><td>Nucleolus</td><td>11</td><td>2.00E-07</td></tr><tr><td>GOTERM_CC_FAT</td><td>Small nucleolar ribonucleoprotein complex</td><td>4</td><td>1.90E-03</td></tr><tr><td>GOTERM_CC_FAT</td><td>Ribonucleoprotein complex</td><td>13</td><td>2.00E-03</td></tr><tr><td>GOTERM_CC_FAT</td><td>Nuclear lumen</td><td>12</td><td>9.40E-03</td></tr><tr><td>GOTERM_CC_FAT</td><td>Eukaryotic transl. initiation factor 3 complex</td><td>4</td><td>8.00E-03</td></tr><tr><td>GOTERM_MF_FAT</td><td>RNA binding</td><td>14</td><td>1.40E-02</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>RNA-binding</td><td>10</td><td>1.80E-03</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>Ribosome biogenesis</td><td>5</td><td>4.70E-03</td></tr><tr><td>SP_PIR_KEYWORDS</td><td>rRNA processing</td><td>5</td><td>4.00E-03</td></tr><tr><td colspan="4" style="text-align: center; font-weight: bold;">Upregulated DESeq2</td></tr><tr><td>GOTERM_BP_FAT</td><td>Amine biosynthetic process</td><td>6</td><td>8.80E-03</td></tr><tr><td>GOTERM_BP_FAT</td><td>Cellular amino acid biosynthetic process</td><td>5</td><td>8.50E-03</td></tr><tr><td>GOTERM_BP_FAT</td><td>Carboxylic acid biosynthetic process</td><td>5</td><td>9.90E-02</td></tr><tr><td>GOTERM_BP_FAT</td><td>Organic acid biosynthetic process</td><td>5</td><td>9.90E-02</td></tr></tbody></table>

Categories: GOTERM_BP_FAT, biological process; GOTERM_CC_FAT, cellular component; GOTERM_MF_FAT, molecular function. SMART & INTERPRO, protein domains; SPIR_KEYWORD, protein information resource provided by SWISSPROT and UniProt.

mapping of reads on the fly genome, since nearly 75% of the 124 million reads (mean number per sample) were properly paired. However, this was not the case for the Tc genome, which displayed an average of less than 1% of properly mapped reads.

This result is not surprising, given that our pre-sequencing manipulations did not target the trypanosome genome. Other contributing factors include the preparation of libraries, which was based on poly(a) selection using Oligo(dT) beads (67) and
<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | BioCyc metabolic map illustrating the different differentially expressed genes involved in tsetse fly metabolic pathways. The genes activated by the infection are displayed in red, whereas the repressed genes are displayed in green.</figcaption></figure>

the *Trypanosoma* genome, which is organized in polycistron units (68).

The infection duration in artificial infection experiments was monitored in the previous report, revealing that the levels of over- or under-expression in DEGs at 3, 10, or 20 days after infection can vary largely (17). In contrast, the present study was performed on tsetse flies sampled in the field, thus neither their age- nor the time-elapsed post-fly infection could be measured. Consequently, the recorded results represent an average level of DEG expression in Gpp flies that may have been infected by trypanosomes (Tc) recently or in the past several days. Similarly, it is possible that non-infected samples could group together flies that were truly never infected with flies that have eliminated their ingested trypanosomes (i.e., “self-cured” or “refractory” flies).

Despite this uncertainty, the results clearly demonstrate a very strong interaction between the parasite and its host/vector, resulting in major transcriptomic changes in the fly. For instance, the level of the “rRNA processing” function in infected vs. non-infected flies was as low as $log2 = -22.3$. In other words, when the infected flies were captured and dissected, the “rRNA processing” function was $2^{22.3} = 5.16 \times 10^6$ fold lower than the value recorded in non-infected flies sampled at the same time and in the same areas. This indicates that the “rRNA processing” function was not effective at that time, and that at least 1 of the 10 DEGs shown to be involved in this function was essentially no longer expressed; however, this does not mean that it could not be reactivated at a later point in a fly’s life.

In this study, we reported that 290 fly genes were down-regulated and 213 were upregulated. This type of imbalance is expected to be induced either by a parasite or a symbiont, and to result in disturbing the host metabolism in such a way as to facilitate microorganism establishment (69, 70). In agreement with this, we observed the repression or non-activation of transcription genes that may allow the trypanosome to alter its host’s transcription steps. Furthermore, certain metabolic pathways were downregulated that can prevent the host from synthesizing factors (proteins or metabolites) needed to fight infection (71). In this context and concerning the “Biological Process,” “Cellular Component,” and “Molecular Function” categories, most of the functional classes were associated with the host transcription/translation machinery (translation, RNA binding, ribonucleoprotein complex, ribosome processing helicase, etc.). Only 10% of the DEGs were related to “Catalytic activities.” In contrast, overexpressed DEGs were involved with catalytic activities, cellular activities (morphogenesis, motion, and cell recognition) and, surprisingly, neuron activities (neuron development and neuron recognition). This is coherently illustrated in **Table 6**, where those “terms” that were over- or under-represented in DEGs (equal to or higher than a fourfold change) and that were identified through functional annotation on the *D. melanogaster* database have been alphabetically classified.

Our identification of the metabolic pathways associated with infection (**Table 8**) highlights the importance of the amino acid
<table><caption>**TABLE 8** | Metabolic pathways associated with fly infection.</caption><thead><tr><th>Pathways</th><th>Functions</th></tr></thead><tbody><tr><td>Amino acid biosynthesis pathway</td><td>Synthesis of L-glutamine<br>Synthesis of L-glutamate<br>Synthesis of L-serine<br>Synthesis of L-asparagine<br>Synthesis of L-aspartate<br>Synthesis of L-valine<br>Synthesis of L-proline<br>Synthesis of L-isoleucine</td></tr><tr><td>Nucleotide biosynthesis pathway</td><td>Biosynthesis of uridine monophosphate</td></tr><tr><td>Pentose phosphate pathway</td><td>Repression of D-ribose-5-phosphate<br>Synthesis of orotidine<br>Synthesis of pyruvate</td></tr><tr><td>Carbohydrate biosynthesis pathway<br>Isolated reactions</td><td>Biosynthesis of pyruvate<br>Synthesis of tyrosine<br>Synthesis of serine<br>Synthesis of L-glutamine<br>Synthesis of L-glutamate<br>Synthesis of L-serine<br>Synthesis of L-asparagine<br>Synthesis of L-methionine<br>Synthesis of L-aspartate<br>Synthesis of L-valine<br>Synthesis of L-isoleucine<br>Synthesis of N-acetylglucosamine<br>Synthesis of galactosyltransferase<br>Synthesis of beta-1,4-manosylglycolipid<br>Synthesis of S-adenosyl-L-homocysteine<br>Repression of immune cytokines<br>Repression of formyltetrahydrofolate DH</td></tr><tr><td>Transport</td><td>Transport of lipids<br>Transport of ATP<br>Transport of succinate<br>Transport of L-carnitine<br>Transport of GTP<br>Transport of acid dicarboxylic<br>Transport of acid monocarboxylic<br>Transport of L-tyrosine<br>Transport of L-serine<br>Transport of Ca<sup>2+</sup><br>Transport of nucleotide<br>Transport of Cyclic GMP<br>Transport of proteinogenic amino acid<br>Transport of NAD<sup>+</sup><br>Transport of L-fructose<br>Transport of GDP<br>Transport of fatty acid<br>Transporter activity<br>Calcium ion binding</td></tr></tbody></table>

biosynthesis pathway. This provides the parasite with a broad range of amino acids that serve as a valuable source of energy, as previously reported for *T. cruzi*, the parasite causing Chagas disease (72), and microsporidia, a parasite of fishes (73). One such amino acid that we identified is proline, whose synthesis was overexpressed in Gpg infected with Tbg in comparison to non-infected flies (17). We also observed an increase in the

biosynthesis of N-acetyl-glucosamine, a molecule that can affix itself to lectins that possess a sugar recognition area (74). This process inactivates tsetse fly lectins that are otherwise lethal to procyclic forms of trypanosomes (65), which consequently favor trypanosome installation in the fly vector. Interestingly, this mechanism has also been reported in Gpg infected with Tbg.

As reported in experimental Gpg insectary flies infected with Tbg, we have shown that field-collected Gpp naturally infected with Tc exhibit a strong cytokine repression in comparison to uninfected tsetse flies. This result indicates that strong alteration of the immune system occurred in infected flies, favoring parasite installation. In addition, *Trypanosoma* infection repressed 34 DEGs encoding non-membrane-bound organelles and 48 DEGs encoding expression of the intracellular lumen (an organelle consisting of chromatin). This type of scenario has also been described for the herpes simplex virus type 1, which can modify the structure and dynamics of chromatin through posttranscriptional modification of histone or other chromatin-forming proteins, contributing to their establishment within the host (75).

This is the first study to evaluate the transcriptomic events associated with infection by the Tc trypanosome in field Gpp tsetse flies. Our results establish that field flies naturally infected by trypanosomes display disruptions in their gene expression that result in either overexpression or under-expression of certain fly genes, as similarly observed in experimentally infected insectary flies. Furthermore, molecular disruptions occur in Gpp when infected with Tc, just as in Gpg that have been artificially infected with Tbg. Importantly, these findings indicate that different *Glossina* species infected with different trypanosome species under different conditions display comparable molecular reactions, which validate the use of experimental host/parasite couples for future research programs.

## AUTHOR CONTRIBUTIONS

AG conceived and designed the experiments. JMTN, FN, BL, GK-N, and NF-N performed the experiments. JMTN, CR, JH, and AG analyzed the data. BL and AG contributed reagents/ materials/analysis tools. JMTN, FN, JH, and AG wrote the paper.

## FUNDING

The authors thank the Région Languedoc-Roussillon-Appel d'Offre «Chercheur d' Avenir 2011» and the Institut de Recherche pour le Développement for their financial support. Highthroughput sequencing was performed at the TGML Platform, supported by grants from Inserm, GIS IBiSA, Aix-Marseille Université, and ANR-10-INBS-0009-10. JMTN is a PhD student supported by an ARTS fellowship from the Institut de Recherche pour le Développement (IRD).

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at http://journal.frontiersin.org/article/10.3389/fimmu.2017.00876/full#supplementary-material.
REFERENCES

1. Buguet A, Bert J, Tapie P, Bogui P, Doua F, Mouanga G, et al. The distribution of sleep and wakefulness in human African trypanosomiasis. *Bull Soc Pathol Exot* (1994) 87:36–7.

2. Mehlitz D, Brinkmann U, Haller L. Epidemiological studies on the animal reservoir of gambiense sleeping sickness. Part I. Review of literature and description of the study areas. *Tropenmed Parasitol* (1981) 32:129–33.

3. Simarro PP, Diarra A, Postigo JAR, Franco JR, Jannin JG. The human African trypanosomiasis control and surveillance programme of the World Health Organization 2000–2009: the way forward. *PLoS Negl Trop Dis* (2011) 5:e1007. doi:10.1371/journal.pntd.0001007

4. Louis FJ, Simarro PP, Lucas P. Sleeping sickness: one hundred years of control strategy evolution. *Bull Soc Pathol Exot* (2002) 95:331–6.

5. Reinhardt E. *Travailler ensemble: la mouche tsé-tsé et la pauvreté rurale*. Chronique ONU (2002). Available from: www.un.org/french/pubs/chronique/2002/numero2/0202p17_la_mouche_tsetse.html

6. Gibson W, Kay C, Peacock L. *Trypanosoma congolense*: molecular toolkit and resources for studying a major livestock pathogen and model trypanosome. *Advances in Parasitology*. (Vol. 98) (2017). doi:10.1016/bs.apar.2017.03.002

7. Van den Abbeele J, Claes Y, Bockstaale D, Ray D, Coosemans M. *Trypanosoma brucei* spp. development in the tsetse fly: characterization of the post-mesocyclic stages in the foregut and proboscis. *Parasitology* (1999) 118:469–78. doi:10.1017/S0031182099004217

8. Peacock L, Cook S, Ferris V, Bailey M, Gibson W. The life cycle of *Trypanosoma (Nannomonas) congolense* in the tsetse fly. *Parasit Vectors* (2012) 5:109. doi:10.1186/1756-3305-5-109

9. Holzmuller P, Grébaut P, Peltier JB, Brizard JP, Perrone T, Gonzatti M, et al. Secretome of animal trypanosomes. *Ann N Y Acad Sci* (2008) 1149:337–42. doi:10.1196/annals.1428.097

10. Grébaut P, Chuchana P, Brizard JP, Demettre E, Seveno M, Bossard G, et al. Identification of total and differentially expressed excreted-secreted proteins from *Trypanosoma congolense* strains exhibiting different virulence and pathogenicity. *Int J Parasitol* (2009) 39:1137–50. doi:10.1016/j.ijpara.2009.02.018

11. Atyame Nten CM, Sommerer N, Rofidal V, Hirtz C, Rossignol M, Cuny G, et al. Excreted/secreted proteins from trypanosome procyclic strains. *J Biomed Biotechnol* (2010) 2010:212817. doi:10.1155/2010/212817

12. Geiger A, Hirtz C, Bécue T, Bellard E, Centeno D, Gargani D, et al. Exocytosis and protein secretion in *Trypanosoma*. *BMC Microbiol* (2010) 10:20. doi:10.1186/1471-2180-10-20

13. Farikou O, Njiokou F, Mbida Mbida JA, Njitchouang GR, Djeunga HN, Asonganyi T, et al. Tripartite interactions between tsetse flies, *Sodalis glossinidius* and trypanosomes—An epidemiological approach in two historical human African trypanosomiasis foci in Cameroon. *Infect Genet Evol* (2010) 10:115–21. doi:10.1016/j.meegid.2009.10.008

14. Geiger A, Ravel S, Mateille T, Janelle J, Patrel D, Cuny G, et al. Vector competence of *Glossina palpalis gambiensis* for *Trypanosoma brucei* s.l. and genetic diversity of the symbiont *Sodalis glossinidius*. *Mol Biol Evol* (2007) 24:102–9. doi:10.1093/molbev/msl135

15. Hamidou Soumana I, Tchicaya B, Simo G, Geiger A. Comparative gene expression of *Wigglesworthia* inhabiting non-infected and *Trypanosoma brucei* gambiense-infected *Glossina palpalis gambiensis* flies. *Front Microbiol* (2014) 5:620. doi:10.3389/fmicb.2014.00620

16. Hamidou Soumana I, Loriod B, Ravel S, Tchicaya B, Simo G, Rihet P, et al. The transcriptional signatures of *Sodalis glossinidius* in the *Glossina palpalis gambiensis* flies negative for *Trypanosoma brucei gambiense* contrast with those of this symbiont in tsetse flies positive for the parasite: possible involvement of a *Sodalis*-hosted prophage in fly *Trypanosoma* refractoriness? *Infect Genet Evol* (2014) 24:41–56.

17. Hamidou Soumana I, Klopp C, Ravel S, Nabi houdine I, Tchicaya B, Parrinello H, et al. RNA-seq de novo assembly reveals differential gene expression in *Glossina palpalis gambiense* infected with *Trypanosoma brucei* gambiense vs. non-infected and self-cured flies. *Front Microbiol* (2015) 6:1259. doi:10.3389/fmicb.2015.01259

18. Hamidou Soumana I, Tchicaya B, Rialle S, Parinello H, Geiger A. Comparative genomics of *Glossina palpalis gambiensis* and *G. morsitans morsitans* to reveal gene orthologs involved in infection by *Trypanosoma brucei* gambiense. *Front Microbiol* (2017) 8:540. doi:10.3389/fmicb.2017.00540

19. Matthews KR. 25 years of African trypanosome research: from description to molecular dissection and new drug discovery. *Mol Biochem Parasitol* (2015) 200:30–40. doi:10.1016/j.molbiopara.2015.01.006

20. Mugnier MR, Stebbins CE, Papavasiliou FN. Masters of disguise: antigenic variation and the VSG coat in *Trypanosoma brucei*. *PLoS Pathog* (2016) 12:e1005784. doi:10.1371/journal.ppat.1005784

21. Tiengwe C, Muratore KA, Bangs JD. Surface proteins, ERAD and antigenic variation in *Trypanosoma brucei*. *Cell Microbiol* (2016) 18:1673–88. doi:10.1111/cmi.12605

22. Odongo S, Sterckx YG, Stijlemans B, Pillay D, Baltz T, Muyldermans S, et al. An anti-proteome nanobody library approach yields a specific immunoassay for *Trypanosoma congolense* diagnosis targeting glycosomal aldolase. *PLoS Negl Trop Dis* (2016) 10:e0004420. doi:10.1371/journal.pntd.0004420

23. WHO. *The World Health Report – Working Together for Health*. (2006). Available from: http://www.who.int/whr/2006/en/

24. Authié E, Boulangé A, Muteti D, Lalmanach G, Gauthier F, Musoke AJ. Immunisation of cattle with cysteine proteinases of *Trypanosoma congolense*: targeting the disease rather than the parasite. *Int J Parasitol* (2001) 31:1429–33. doi:10.1016/S0020-7519(01)00266-1

25. Huson LE, Authié E, Boulangé AF, Goldring JP, Coetzer TH. Modulation of the immunogenicity of the *Trypanosoma congolense* cysteine protease, con-gopain, through complexation with alpha(2)-macroglobulin. *Vet Res* (2009) 40:52. doi:10.1051/vetres/2009036

26. Aksoy S, Weiss B, Attardo G. Paratransgenesis applied for control of tsetse transmitted sleeping sickness. In: Aksoy S, editor. *Transgenesis and the Management of Vector-Borne Disease, Advances in Experimental Medicine and Biology*. New York: Springer (2008). p. 35–48.

27. Matetovici I, Caljon G, Van Den Abbeele J. Tsetse fly tolerance to *T. brucei* infection: transcriptome analysis of trypanosome-associated changes in the tsetse fly salivary gland. *BMC Genomics* (2016) 17:971. doi:10.1186/s12864-016-3283-0

28. Rio RVM, Hu Y, Aksoy S. Strategies of the home-team: symbioses exploited for vector-borne disease control. *Trends Microbiol* (2004) 12:325–36. doi:10.1016/j.tim.2004.05.001

29. Abd-Alla AMM, Bergoin M, Parker AG, Maniania NK, Vlak JM, Bourtzis K, et al. Improving Sterile Insect Technique (SIT) for tsetse flies through research on their symbionts and pathogens. *J Invertebr Pathol* (2013) 112:S2–10. doi:10.1016/j.jip.2012.07.009

30. Vreysen MJB, Saleh KM, Ali MY, Abdulla AM, Zhu Z-R, Juma KG, et al. *Glossina austeni* (Diptera: Glossinidae) eradicated on the island of Unguja, Zanzibar, using the sterile insect technique. *J Econ Entomol* (2000) 93:123–35. doi:10.1603/0022-0493-93.1-123

31. Geiger A, Ponton F, Simo G. Adult blood-feeding tsetse flies, trypanosomes, microbiota and the fluctuating environment in sub-Saharan Africa. *ISME J* (2015) 9:1496–507. doi:10.1038/ismej.2014.236

32. Geiger A, Bossard G, Sereno D, Pissarra J, Lemesre JL, Vincendeau P, et al. Escaping deleterious immune response in their hosts: lessons from trypanosomatids. *Front Immunol* (2016) 7:212. doi:10.3389/fimmu.2016.00212

33. Frezil J-L, Cuisance D. Trypanosomiases, maladies d'avenir: leurs perspectives et leurs inconnues. *Trypanosomiasis, diseases with future: prospects and uncertainty*. *Bull Soc Pathol Exot* (1994) 87:391–3.

34. Aksoy S, Gibson WC, Lehane MJ. Interactions between tsetse and trypanosomes with implications for the control of trypanosomiasis. *Adv Parasitol* (2003) 53:1–83. doi:10.1016/S0065-308X(03)53002-0

35. Lehane MJ, Aksoy S, Gibson W, Kerhornou A, Berriman M, Hamilton J, et al. Adult midgut expressed sequence tags from the tsetse fly *Glossina morsitans morsitans* and expression analysis of putative immune response genes. *Genome Biol* (2003) 4:R63. doi:10.1186/gb-2003-4-10-r63

36. Lehane MJ, Gibson W, Lehane SM. Differential expression of fat body genes in *Glossina morsitans morsitans* following infection with *Trypanosoma brucei brucei*. *Int J Parasitol* (2008) 38:93–101. doi:10.1016/j.ijpara.2007.06.004

37. Grébaut P, Melachio T, Nyangmang S, Eyenga VE, Njitchouang GR, Ofon E, et al. Xenomonitoring of sleeping sickness transmission in Campo (Cameroon). *Parasit Vectors* (2016) 9:201. doi:10.1186/s13071-016-1479-4

38. Grebaut P, Bodo JM, Assona A, Foumane Ngane V, Njiokou F, Ollivier G, et al. Risk factors for human African trypanosomiasis in the Bipindi region of Cameroon. *Med Trop Rev Corps Sante Colon* (2001) 61:377–83.
39. Farikou O, Njiokou F, Simo G, Asonganyi T, Cuny G, Geiger A. Tsetse fly blood meal modification and trypanosome identification in two sleeping sickness foci in the forest of southern Cameroon. *Acta Trop* (2010) 116:81–8. doi:10.1016/j.actatropica.2010.06.002

40. Njiokou F, Simo G, Nkinin SW, Laveissière C, Herder S. Infection rate of *Trypanosoma brucei* s.l., *T. vivax*, *T. congolense* "forest type", and *T. simiae* in small wild vertebrates in south Cameroon. *Acta Trop* (2004) 92:139–46. doi:10.1016/j.actatropica.2004.04.011

41. Njiokou F, Laveissière C, Simo G, Nkinin S, Grébaut P, Cuny G, et al. Wild fauna as a probable animal reservoir for *Trypanosoma brucei gambiense* in Cameroon. *Infect Genet Evol* (2006) 6:147–53. doi:10.1016/j.meegid.2005.04.003

42. Simo G, Njiokou F, Mbida Mbida JA, Njitchouang GR, Herder S, Asonganyi T, et al. Tsetse fly host preference from sleeping sickness foci in Cameroon: epidemiological implications. *Infect Genet Evol* (2008) 8:34–9. doi:10.1016/j. meegid.2007.09.005

43. Lancien J. Description du piège monoconique utilisé pour l'élimination des glossines en République Populaire du Congo. *Ent Méd Parasitol* (1981) 19:235–8.

44. Grébaut P, Mbida JAM, Kondjio CA, Njiokou F, Penchenier L, Laveissière C. Spatial and temporal patterns of human African trypanosomosis (HAT) transmission risk in the bipindi focus, in the forest zone of southern Cameroon. *Vector Borne Zoonotic Dis* (2004) 4:230–8. doi:10.1089/vbz.2004.4.230

45. Penchenier L, Itard J. Une nouvelle technique de dissection rapide des glandes salivaires et de l'intestin de glossines. *Ent Méd Parasitol* (1981) 19:55–7.

46. Herder S, Simo G, Nkinin S, Njiokou F. Identification of trypanosomes in wild animals from southern Cameroon using the polymerase chain reaction (PCR). *Parasite* (2002) 9:345–9. doi:10.1051/parasite/2002094345

47. Masiga DK, Smyth AJ, Ayes P, Bromidge TJ, Gibson WC. Sensitive detection of trypanosomes in tsetse flies by DNA amplification. *Int J Parasitol* (1992) 22:909–18. doi:10.1016/0020-7519(92)90047-0

48. Moser DR, Cook GA, Ochs DE, Bailey CP, McKane MR, Donelson JE. Detection of *Trypanosoma congolense* and *Trypanosoma brucei* subspecies by DNA amplification using the polymerase chain reaction. *Parasitology* (1989) 99:57–66. doi:10.1017/S0031182000061023

49. Köster J, Rahmann S. Snakemake – a scalable bioinformatics workflow engine. *Bioinformatics* (2012) 28:2520–2. doi:10.1093/bioinformatics/bts480

50. Giraldo-Calderón GI, Emrich SJ, MacCallum RM, Maslen G, Dialynas E, Topalis P, et al. VectorBase: an updated bioinformatics resource for invertebrate vectors and other organisms related with human diseases. *Nucleic Acids Res* (2015) 43:D707–13. doi:10.1093/nar/gku1117

51. Aslett M, Aurrecoechea C, Berriman M, Brestelli J, Brunk BP, Carrington M, et al. TriTrypDB: a functional genomic resource for the *Trypanosomatidae*. *Nucleic Acids Res* (2010) 38:D457–62. doi:10.1093/nar/gkp851

52. Attrill H, Falls K, Goodman JL, Millburn GH, Antonazzo G, Rey AJ, et al. FlyBase: establishing a gene group resource for *Drosophila melanogaster*. *Nucleic Acids Res* (2016) 44:D786–92. doi:10.1093/nar/gkv1046

53. Liao Y, Smyth GK, Shi W. The subread aligner: fast, accurate and scalable read mapping by seed-and-vote. *Nucleic Acids Res* (2013) 41:e108. doi:10.1093/nar/gkt214

54. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. *Bioinformatics* (2009) 25:2078–9. doi:10.1093/bioinformatics/btt352

55. Liao Y, Smyth GK, Shi W. FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features. *Bioinformatics* (2014) 30:923–30. doi:10.1093/bioinformatics/btt656

56. Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. *Nat Biotechnol* (2011) 29:24–6. doi:10.1038/nbt.1754

57. Varet H, Brillet-Guéguen L, Coppée J-Y, Dillies M-A. SARTools: a DESeq2- and EdgeR-based R pipeline for comprehensive differential analysis of RNA-Seq data. *PLoS One* (2016) 11:e0157022. doi:10.1371/journal.pone.0157022

58. 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

59. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. *Bioinformatics* (2010) 26:139–40. doi:10.1093/bioinformatics/btp616

60. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. *Nucleic Acids Res* (1997) 25:3389–402. doi:10.1093/nar/25.17.3389

61. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. *Nat Protoc* (2009) 4:44–57. doi:10.1038/nprot.2008.211

62. Reimand J, Kull M, Peterson H, Hansen J, Vilo J. g:Profiler—a web-based toolset for functional profiling of gene lists from large-scale experiments. *Nucleic Acids Res* (2007) 35:W193–200. doi:10.1093/nar/gkm226

63. Caspi R, Billington R, Ferrer L, Foerster H, Fulcher CA, Keseler IM, et al. The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of pathway/genome databases. *Nucleic Acids Res* (2016) 44:D471–80. doi:10.1093/nar/gkv1164

64. Maudlin I, Ellis DS. Association between intracellular *Rickettsia*-like infections of midgut cells and susceptibility to trypanosome infection in *Glossina* spp. *Z Parasitenkd* (1985) 71:683–7. doi:10.1007/BF00925601

65. Welburn SC, Maudlin I. Tsetse-trypanosome interactions: rites of passage. *Parasitol Today* (1999) 15:399–403. doi:10.1016/S0169-4758(99)01512-4

66. Dale C, Welburn SC. The endosymbionts of tsetse flies: manipulating host-parasite interactions. *Int J Parasitol* (2001) 31:628–31. doi:10.1016/S0020-7519(01)00151-5

67. Koch H, Raabe M, Urlaub H, Bindereif A, Preußer C. The polyadenylation complex of *Trypanosoma brucei*: characterization of the functional poly(A) polymerase. *RNA Biol* (2016) 13:221–31. doi:10.1080/15476286.2015.1130208

68. Kolev NG, Franklin JB, Carmi S, Shi H, Michaeli S, Tschudi C. The transcriptome of the human pathogen *Trypanosoma brucei* at single-nucleotide resolution. *PLoS Pathog* (2010) 6:e1001090. doi:10.1371/journal.ppat.1001090

69. Martinson EO, Wheeler D, Wright J, Mrinalini Siebert AL, Werren JH. *Nasonia vitripennis* venom causes targeted gene expression changes in its fly host. *Mol Ecol* (2014) 23:5918–30. doi:10.1111/mec.12967

70. Snyder AK, Rio RVM. "Wigglesworthia morsitans" folate (vitamin b9) biosynthesis contributes to tsetse host fitness. *Appl Environ Microbiol* (2015) 81:5375-86. doi:10.1128/AEM.00553-15

71. Zhang Q, Lai M-M, Lou Y-Y, Guo B-H, Wang H-Y, Zheng X-Q. Transcriptome altered by latent human cytomegalovirus infection on THP-1 cells using RNA-seq. *Gene* (2016) 594:144–50. doi:10.1016/j.gene.2016.09.014

72. Barisón MJ, Damasceno FS, Mantilla BS, Silber AM. The active transport of histidine and its role in ATP production in *Trypanosoma cruzi*. *J Bioenerg Biomembr* (2016) 48:437–49. doi:10.1007/s10863-016-9665-9

73. Dean P, Hirt RP, Embley TM. *Microsporidia*: why make nucleotides if you can steal them? *PLoS Pathog* (2016) 12:e1005870. doi:10.1371/journal.ppat.1005870

74. Kanost MR, Jiang H, Yu X-Q. Innate immune responses of a lepidopteran insect, *Manduca sexta*. *Immunol Rev* (2004) 198:97–105. doi:10.1111/j.0105-2896.2004.0121.x

75. Kent JR, Zeng P-Y, Atanasiu D, Gardner J, Fraser NW, Berger SL. During lytic infection herpes simplex virus type 1 is associated with histones bearing modifications that correlate with active transcription. *J Virol* (2004) 78:10178–86. doi:10.1128/JVI.78.18.10178-10186.2004

**Conflict of Interest Statement:** 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.

Copyright © 2017 Tsagmo Ngoune, Njiokou, Loriod, Kame-Ngasse, Fernandez-Nunez, Rioualen, van Helden and Geiger. 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) or licensor 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.