Global characterization of gene expression in the brain of starved immature Rhodnius prolixus

Abstract

Background

Rhodnius prolixus is a vector of Chagas disease and has become a model organism to study physiology, behavior, and pathogen interaction. The publication of its genome allowed initiating a process of comparative characterization of the gene expression profiles of diverse organs exposed to varying conditions. Brain processes control the expression of behavior and, as such, mediate immediate adjustment to a changing environment, allowing organisms to maximize their chances to survive and reproduce. The expression of fundamental behavioral processes like feeding requires fine control in triatomines because they obtain their blood meals from potential predators. Therefore, the characterization of gene expression profiles of key components modulating behavior in brain processes, like those of neuropeptide precursors and their receptors, seems fundamental. Here we study global gene expression profiles in the brain of starved R. prolixus fifth instar nymphs by means of RNA sequencing (RNA-Seq).

Results

The expression of neuromodulatory genes such as those of precursors of neuropeptides, neurohormones, and their receptors; as well as the enzymes involved in the biosynthesis and processing of neuropeptides and biogenic amines were fully characterized. Other important gene targets such as neurotransmitter receptors, nuclear receptors, clock genes, sensory receptors, and takeouts genes were identified and their gene expression analyzed.

Conclusion

We propose that the set of neuromodulatory-related genes highly expressed in the brain of starved R. prolixus nymphs deserves functional characterization to allow the subsequent development of tools targeting them for bug control. As the brain is a complex structure that presents functionally specialized areas, future studies should focus on characterizing gene expression profiles in target areas, e.g. mushroom bodies, to complement our current knowledge.

Introduction

Triatomines are hematophagous insects that can transmit Trypanosoma cruzi, the etiological agent of Chagas disease. It is estimated that this neglected disease affects 7 million people, located mostly in Central and South America. The study of their biology is relevant because T. cruzi transmission is mostly controlled by eliminating domiciliated bugs [1].

Triatomines are nocturnal insects that assume an akinetic state while hidden in shelters during daylight hours. At nightfall, they eventually start a non-oriented locomotor activity, outside shelters to search for hosts. For host recognition, starved bugs detect cues released by vertebrates, such as radiant heat, water vapor, carbon dioxide, and other odorants [2]. The decision to leave a shelter and engage in foraging is risky, as triatomine hosts are often predators as well. For this reason, starved bugs mostly leave the protection of the shelters when a robust set of host clues is present [3–5]. Bugs of all nymphal instars and adults of both sexes feed on blood and can tolerate long starvation. Whereas nymphs have to feed to be able to molt, and adult females require nutrients to produce eggs [6]. Starved insects orientate towards host-emitted stimuli, while fed insects can remain indifferent or avoid these cues depending on the time elapsed after feeding [7].

The central nervous system (CNS) is the main regulator of physiology and behavior. Besides processing sensory information, the brain is the major accumulation of neuropiles integrating neural activity of sensory, memory, and proprioceptive nature [8]. As such, it has a main role in the coordination of motor responses, adjusting their proper timing through a set of clock neurons [9–11]. Signal transfer and modulation of neural processes in the CNS depend on neuroactive compounds, including neurotransmitters of diverse chemical nature like biogenic amines and neuropeptides, and their receptors [12, 13]. Neuropeptides and biogenic amines can also act as endocrine factors mediating signaling processes in multicellular organisms, and in the case of insects, they are fundamental in coordinating growth and development, as well as physiological processes such as metabolism, diuresis, digestion, reproduction, and behavior [14, 15]. Rhodnius prolixus Stål, 1859 (Hemiptera, Reduviidae, Triatominae) is considered an important vector of Chagas disease in Colombia and Venezuela due to its adaptability to both domestic and peridomestic environments, its rapid developmental cycle and the great population density it reaches in human dwellings [16–19]. Furthermore, R. prolixus has been widely used as a model for insect physiology studies, including research on reproduction, development, immunology, and vector-parasite interactions [20]. After the publication of its genome sequence [21], and the introduction of next-generation sequencing (NGS) methods, several studies have described genetic and molecular components underlying the physiology of R. prolixus [14, 15, 22–30]. Transcriptomic studies allowed the discovery of new genes and transcripts, the identification of differentially expressed genes, and determining targets for broader functional analyses. Some transcriptomes have analyzed gene expression in different R. prolixus tissues such as salivary glands [31], ovaries [32, 33], gut [29], testicles [34], and antennae [26]. Furthermore, this technique allowed defining the molecular bases of female reproductive physiology under differing nutritional states [22], as well as characterizing the innate immune system of these bugs at the molecular level [35]. Several of these bioinformatic analyses have recently shown that the genome of R. prolixus has many missing or miss-annotated genes [22, 33, 36], highlighting the importance of transcriptomes for improving the quality of the annotated genome of R. prolixus. Therefore, the present study aims to describe the genetic components that serve as the molecular neural bases controlling behavior in the brain of unfed R. prolixus nymphs.

Materials and methods

Insects

Rhodnius prolixus were obtained from a colony derived from insects collected in Honduras around 1990 and maintained by the Vector Behavior and Pathogen Interaction Group at the René Rachou Institute, Belo Horizonte, Brazil. Insects were monthly fed on citrated rabbit blood obtained from CECAL (Centro de Criação de Animais de Laboratório, FIOCRUZ, Rio de Janeiro, Brazil) offered through an artificial feeder at 37°C, alternating with feeding on anesthetized chicken and mice. Chickens were anesthetized with intraperitoneal injections of a mixture of ketamine (20 mg/kg; Cristália, Brazil) and detomidine (0.3 mg/kg; Syntec, Brazil), and mice with ketamine (150 mg/kg; Cristália, Brazil) and xylazine (10 mg/kg; Bayer, Brazil). Insects were reared in the insectary under 27 ± 2°C, 51 ± 7% of relative humidity and natural illumination. For this study, seven-day-old 4^th instar nymphs were fed on citrated rabbit blood using an artificial feeder. Insects were kept unfed for 30 days after their ecdysis to the 5^th instar and subsequently dissected to collect their brains. This study was carried out in strict accordance with the recommendations in CONCEA/MCT (http://www.cobea.org.br/), which is associated with the American Association for Animal Science (AAAS), the Federation of European Laboratory Animal Science Associations (FELASA), the International Council for Animal Science (ICLAS) and the Association for Assessment and Accreditation of Laboratory Animal Care International (AAALAC). The protocol was approved by the Committee for Ethics in the Use of Animals, CEUA, of the FIOCRUZ (Protocol Number: LW-61/2012).

RNA extraction and illumina sequencing

Insect brains were dissected on a freeze cold dissecting dish (BioQuip, Gardena, CA, US), collected with forceps, and immediately transferred to a microtube immersed in dry ice and added 1 mL of TRIzol™ reagent (Invitrogen, Thermo Fisher Scientific, MA, USA). For sample completion, dissections occurred along three days, exclusively between 2 and 4 pm. Two separate experiments were performed, with three independent replicates each. Samples for replicates were, each composed of a pool of 20 brains. RNA extraction was performed with TRIzol™ according to the manufacturer’s instructions. Total RNA concentrations were determined using a Qubit 2.0 Fluorometer (Life Technologies, Carlsbad, CA, US). Six libraries were constructed using the TruSeq Stranded mRNA Sample Preparation Kit (Illumina, San Diego, CA) and sequenced on an Illumina HiSeq 2500 platform at the Max Planck Genome Center in Cologne (Germany). Approximately 15 million reads were obtained for each library, using 150 base-pair (bp) paired-end reads. The raw sequence dataset is available with the NCBI-SRA Bioproject number PRJNA853796 at NCBI.

Bioinformatic analysis

Raw reads were filtered and trimmed for low-quality bases using Trimmomatic (v0.36) [37], according to standard quality score parameters (Phred-33 (>15); and 50 base-pair minimum length). Then, STAR v2.6.0 [38] was used with default parameters to map reads to the R. prolixus reference genome (version RproC3.3) accessed through the VectorBase website [39]. Mapped reads were assigned to each gene through the coverage tool in BEDTools (v2.29.2) based on an updated gene annotation file [26]. A Principal Component Analysis (PCA) was performed through the plotPCA function in DESeq2 on RStudio to analyze the variation between the six libraries. Gene length and total counts of mapped reads were used for calculating Transcripts per Kilobase per Million mapped reads (TPM) values for target genes in each library. Subsequently, target gene expression (as Log₁₀ TPM+1) was depicted in heatmaps built using the pheatmap R package (v1.0.12). Finally, all genes were ranked according to the highest expression using TPM values, and the top 50 that were present in at least four of the six libraries, and which presented annotation in VectorBase [39] and/or a positive hit in BLASTp searches were selected. BLASTp searches were performed against Insecta class sequences in GenBank to identify putative functions of these highly expressed genes.

Results and discussion

Overall analysis

RNA-Seq data from starved 5^th instar R. prolixus brain transcriptomes were summarized in Table 1. After filtering and trimming raw reads, all libraries showed coverage of around 13 million reads. The number of uniquely mapped reads against the R. prolixus genome ranged from 8,8 M to 10,3 M reads. According to the PCA graph, three libraries clustered together (Rep2, Rep3 from Experiment 1 and Rep4 from Experiment 2), and three segregated apart (S1 Fig).

[Figure omitted. See PDF.]

Based on our brain library outputs, six lists of the expressed genes were obtained after ranking their TPM values. To characterize the set of genes most highly expressed in the brain, we compared these ranks and built a consensus list depicting genes that ranked top 50 in at least four out of six libraries (Table 2). Several initiation and elongation factors together with ribosomal proteins were identified among top expressed genes, probably reflecting protein biosynthesis induced by starvation. Several heat shock proteins presented high expression in the brain, also probably due to starvation-generated stress. Different types of soluble carrier proteins, like odorant binding proteins, takeouts and lipocalins were very abundantly expressed in the CNS, suggesting roles other than odor transportation or detection. Finally, other highly expressed genes were the neuroendocrine secretory protein 7B2 which functions as a specific chaperone for the prohormone convertase 2 (PC2) [40], an enzyme required for the maturation of neuropeptide and peptide hormone precursors; and the Glutamine synthetase that catalyzes the synthesis of glutamine, which has a central role in nitrogen metabolism and the regulation of neurotransmitter production [41].

[Figure omitted. See PDF.]

Neuropeptide precursor genes

Neuropeptide precursor gene (NPG) expression patterns in insects tend to be stereotyped, and each neuropeptide may be involved in neurotransmission, synaptic neuromodulation, or in conveying neuroendocrine signals at peripheral targets [12]. Currently, structural and functional studies on insect neuropeptides are being performed to develop new insect control approaches. This is especially true for neuropeptides involved in developmental, nutritional, and survival processes [13]. For this reason, we focused on the neuropeptides showing the highest expression in our study.

The four most highly expressed NPGs (Neuropeptide-like precursor 1—NPLP1, Tachykinin—TK, ITG-like, and NVP-like) presented expression values > 3 (Log₁₀ TPM+1). Most genes coding for neuropeptides (26 genes) had expression values between two and three, while eight had values between 1 and 2. A few of them, like ecdysis triggering hormone (ETH), elevenin 1 (Ele1), eclosion hormone (EH), sulphakinin (SK), and sifamide (SIFa) had expression values < 1. Adipokinetic hormone/corazonin-related peptide (ACP) was the only NPG out of 44 annotated for R. prolixus [26] that was not expressed in any of the six libraries here analyzed (Fig 1A).

[Figure omitted. See PDF.]

(A) Heatmap depicting the expression level of neuropeptide precursor genes, and (B) neuropeptide receptor genes in the brain of R. prolixus nymphs. Expression (displayed as Log₁₀ TPM+1) is represented by means of a color scale in which blue/red represent the lowest/highest expression. Each column represents the expression of one library.

Neuropeptide-like precursor 1 (NPLP1) was the NPG showing the highest expression in our transcriptomic database. Furthermore, it ranked amongst the 10 most highly expressed genes in the brain, presenting a mean of 3.90±0.03 Log₁₀ TPM+1 among our libraries. There were two different isoforms with distinct expression patterns previously detected in R. prolixus, isoform A expressed in CNS, ovaries, testes, and antennae, while expression of isoform B was only detected in antennae. The presence of mature peptides from NPLP1 precursors was also previously detected in the salivary glands of R. prolixus using proteomics [27]. The high expression observed for this neuropeptide precursor in our work suggests a relevant role in the CNS. Indeed, Sterkel et al. [30] showed that the levels of two mature peptides encoded by the NPLP1 precursor gene significantly decreased in the CNS 24 hours after feeding, suggesting a role connected to starvation. A decrease in NPLP1 expression was also observed in the antennal lobe of mosquitoes after blood or sugar ingestion [42].

According to our brain dataset, TK was the NPG showing the second-highest expression level (3.22 ± 0.07 Log₁₀ TPM+1). Sterkel et al. [30] have also detected TK in the CNS of R. prolixus by means of peptidomics, not observing changes in its abundance after blood ingestion. This NPG presents expression in a wide set of R. prolixus tissues, as its transcripts have been detected in salivary glands, fat body, dorsal vessel, the intestinal tract of 5^th instar nymphs (by RT-qPCR) [43], and in the antennae of 5^th instar nymphs and adults (using RNA-Seq) [26]. Studies on other insects have shown that TKs can modulate early olfactory processing at the olfactory lobes, circuits controlling locomotion and food search, aggression, metabolic stress, and nociception [44]. Similar roles could be expected in triatomines based on the high expression of the TK gene observed in our study.

ITG-like and NVP-like genes also presented high expression values (3.21±0.08 and 3.10±0.04 Log₁₀ (TPM+1), respectively) in our database; nonetheless, to our knowledge, a functional characterization is not available for these peptides so far. Latorre-Estivalis et al. [26] detected high expression of ITG-like in the antennae of immature and adults of R. prolixus and proposed a modulatory role for this neuropeptide at the peripheral level. In its turn, Leyria et al. [22] indicated that blood ingestion induced ovarian downregulation of the NVP-like gene, suggesting a role in R. prolixus reproduction [22]. The quantitative peptidomics analysis of the R. prolixus CNS by Sterkel et al. [30] showed that the abundance of ITG-like and NVP-like neuropeptides significantly decreases a few hours post-blood meal, indicating their implication in a neuroendocrine response to feeding. Considering this and the very high expression observed for ITG-like and NVP-like in our dataset from starved bugs, we suggest that they might act by signaling starvation status. Functional genetic studies based on gene silencing should be implemented to verify potential behavioral phenotypes that could offer evidence about the roles of these neuropeptides in R. prolixus.

Insulin-like peptide (ILP) presented a high expression in our database (2.90±0.02 Log₁₀ TPM+1). This result was consistent with the characteristics of our samples (brains from starved bugs) and the putative function of this neuropeptide as a modulator of lipid and carbohydrate metabolism [45, 46]. Indeed, in vitro immunofluorescence studies in R. prolixus brains demonstrated strong and abundant fluorescence of ILP neurons in unfed nymphs, followed by an acute decrease 4 hours after ingestion of a food meal, indicating transport and release of these signaling molecules into the hemolymph soon after feeding [47]. The high expression of this neuropeptide in the CNS of unfed insects was also observed by Leyria et al. [22]. Interestingly, these authors observed that blood ingestion did not affect ILP gene expression in the CNS of R. prolixus [22].

Neuropeptide processing enzymes

Mature neuropeptides are synthesized by a series of enzymatic steps that sequentially cleave and modify larger precursor molecules. One step of neuropeptide biosynthesis involves peptide amidation, a process that occurs on half of the known bioactive neuropeptides [48, 49]. The expression pattern observed for these enzymes might be a good proxy to estimate their abundance and activity [50], deserving their subsequent functional characterization.

A set of eleven neuropeptide processing enzymes was previously annotated for R. prolixus [26], using sequences of D. melanogaster as references [51]. All of these enzymes showed relative expression values higher than 1 (Log₁₀ TPM+1) (S2 Fig). Peptidyl alpha hydroxyglycine alpha amidatinglyase 2 (PAL2) and peptidylglycine α-amidating monooxygenase (PHM), the enzymes involved in amidation reactions, were the most highly expressed genes of this set (PAL2–2.92±0.07, PHM– 2.77±0.11).

Neuropeptide receptors

Whether a tissue is targeted by a certain neuropeptide is defined by the presence of the corresponding neuropeptide receptor on the surface of its cells. Neuropeptide receptor genes showed much lower expression values than neuropeptide gene precursors. Only two out of the 48 neuropeptide receptor genes analyzed showed a relative expression higher than 2 (Fig 1B); allatostatin C (AstC) receptor with 2.38±0.05 Log₁₀ TPM+1 and calcitonin-like diuretic hormone receptor 3 -CT/DH-R3 with 2.12±0.03 Log₁₀ TPM+1. Twenty-six genes (54%) presented a relative expression between 2 and 1, and twenty genes (41%) had values lower than 1.

The expression of the AstC receptor in the R. prolixus CNS had been previously reported by Ons et al. [15] using RNA-Seq. These authors did not see expression changes for this receptor at the CNS after blood ingestion. Furthermore, Villalobos-Sambucaro et al. [52] observed the presence of this receptor in the hindgut, midgut and dorsal vessel, and showed that the receptor and its ligand play a key myoregulatory and cardioregulatory role in R. prolixus. Our results suggest that the AstC receptor may also have a fundamental role at the central level. Two receptors for CT/DH (named R1 and R2) were previously described in R. prolixus, their transcript expression being detected in the CNS and reproductive tissues [53]. The existence of a third CT/DH receptor in R. prolixus was suggested by Ons et al. [15] and confirmed by an antennal transcriptome [26]. The latter showed that this receptor had increased expression in male antennae when compared to those of nymphs, suggesting a stage-enriched role for this gene [26]. Our results suggest that this receptor may also act at the central level.

Neurotransmitter receptors

As expected, neurotransmitter receptors tended to present higher expression values than the other receptor gene families studied here. N-methyl-D-aspartate receptor type 2A (NMDAr2a) was the second gene with the highest expression in the whole brain transcriptome (Fig 2; 4.05±0.03 Log₁₀ TPM+1). Muscarinic acetylcholine receptor type A (AChR-A) and dopamine 1-like receptor 1 (DOP1) were the other two genes that presented FPKM values higher than 2 in this gene set (2.16±0.12 Log₁₀ TPM+1 and 2.15±0.01 Log₁₀ TPM+1) for AChR-A and DOP1, respectively). Only two genes, 5-hydroxytryptamine (serotonin) receptors 1A (5HT-1A-R) and 7A (5HT-7A-R) showed values lower than 1 (0.62±0.08 and 0.82+0.08 Log₁₀ TPM+1) for 5HT-1A-R and 5HT-7A-R, respectively.

[Figure omitted. See PDF.]

Heatmap depicting the expression level of neurotransmitter receptor genes in the brain of R. prolixus nymphs. Expression level (displayed as Log₁₀ TPM+1) is represented by means of a color scale in which blue/red represent the lowest/highest expression. Each column represents the expression of one library.

N-methyl-D-aspartate (NMDA) receptors are one of the subtypes of ionotropic receptors that bind to L-glutamate, mediating an excitatory activity in the CNS of insects. The NMDARs are usually constituted of two subunits NR1 and NR2 [54]. Even though these receptors were characterized in the brain of several invertebrate species, their functions in insects are poorly understood. However, their involvement in behavioral plasticity is already known [55, 56]. Similarly, a study of the evaluation of NMDAR expression in different tissues from female Dactyola punctata, showed that DpunNR1A, DpunNR1B and DpunNR2 were highly expressed in the brain [57]. In Drosophila melanogaster, both NMDA receptors called DmelNR1 and DmelNR2 were weakly expressed throughout the entire brain, with higher expression observed in some scattered cell bodies [58]. As far as we know, this is the first report on the expression of this receptor in the brain of R. prolixus; its high expression suggests a very relevant role in the neural physiology of these insects. Immunostaining experiments to characterize brain neuropiles depicting NMDAr2a expression will be required to initiate functional studies to uncover its putative function.

Nuclear receptors

Most nuclear receptor genes presented low expression in the brain of unfed nymphs (Fig 3). Out of this gene set, the ecdysone-induced protein 75B (Eip75B) receptor was the only gene that presented TPM values higher than 2 (2.24±0.07 Log₁₀ TPM+1). Ten out of 23 genes (43%) presented values between 1 and 2 Log₁₀ TPM+1). The Eip75B and HR51 transcripts (the latter also known as unfulfilled) have been identified in central clock cells of D. melanogaster and control the expression of clock genes, playing an important role in the maintenance of locomotor rhythms [59–61]. Similar roles could be proposed for R. prolixus, however, functional information is not available for Eip75B and HR51 genes in this species; only Latorre-Estivalis et al. [26] reported similar expression values for the Eip75B gene in kissing bug antennae.

[Figure omitted. See PDF.]

Heatmap depicting the expression level of nuclear receptor genes in the brain of R. prolixus nymphs. Expression level (displayed as Log₁₀ TPM+1) is represented by means of a color scale in which blue/red represent the lowest/highest expression. Each column represents the expression of one library.

Clock and behavior-related genes

Clock genes are responsible for controlling circadian rhythms, and they can cycle in a synchronized way according to daily oscillations of environmental cues such as light and temperature [62]. Even though most available information on clock gene function has been generated using D. melanogaster as a model, their fundamental roles and their high level of sequence homology suggest that their functions should be conserved across insect orders. A total of 31 clock genes have been previously described and annotated in the R. prolixus genome [21]; however, as far as we know, this is the first time that the expression of the whole clock gene set is studied in this insect. As all our brain samples were generated at the same interval of the daily cycle (2–4 PM), the expression profiles obtained here define the levels of expression of clock genes at this time. Except for vrille (vri), cycle (cyc), single-minded (sim) and timeless (tim) genes, the rest of the clock genes had expression higher than 1 Log₁₀ TPM+1 values (Fig 4). The most highly expressed genes (> 2 Log₁₀ TPM+1) were: casein kinase 2 (ck2), protein phosphatases (Pp) 1a and 2a, no circadian temperature entrainment (nocte), poly(A) binding proteins (pAbp) 1 and 2, and no receptor potential A (norpA). These genes showing high expression probably have key roles in the brain of unfed kissing bugs. Based on D. melanogaster studies, pAbps genes form a complex with twenty-four (tyf) and Ataxin-2 (Atx2) that maintain circadian rhythms in locomotor behavior [63]. The Pp2A gene, which is also highly expressed, controls the cyclic expression of the PER protein [64], which was also detected in our dataset. The nocte gene encodes a protein involved in temperature compensation of the circadian clock in Drosophila. It would be relevant to characterize 24h expression profiles of clock genes in the R. prolixus CNS to identify genes with cycling expression.

[Figure omitted. See PDF.]

Heatmap depicting the expression level of clock-related genes in the brain of R. prolixus nymphs. Expression level (displayed as Log₁₀ TPM+1) is represented by means of a color scale in which blue/red represent the lowest/highest expression. Each column represents the expression of one library.

Other genes related to the control of insect behavior include foraging (for) whose expression was detected in our brain transcriptome (S1 Table). This gene encodes a cGMP-dependent protein kinase and plays an essential role in modulating food search in different species of insects, such as D. melanogaster [65–67], locusts [68], ants [69], honeybees [70], and social wasps [71]. In R. prolixus, for has been shown to participate in the modulation of locomotory activity [4, 72], and its expression in the brain and fat body changes depending on the nutritional status of the insect, increasing with starvation [72]. Our transcriptomic data seem to reinforce its relevance in the brain of R. prolixus.

Sensory-related genes

Chemosensory proteins (CSPs) and odorant binding proteins (OBPs) can bind, solubilize and transport hydrophobic molecules [73]. The role of these transporters has been mainly studied in insect antennae and other sensory tissues where they bind odor molecules [73, 74].

Nevertheless, the expression of certain CSPs and OBPs has been reported in non-sensory tissues and involved in different functions, like releasing semichemicals in pheromone glands or associated to insecticide resistance, among others (revised by Pelosi et al. 2018) [75]. For this reason, we decided to data mine our database to characterize whether representatives of these protein families showed expression in bug brains (S2 Fig). Interestingly, high expression levels of several OBP and CSP transcripts, such as those of RproCsp3, RproCsp5, RproCsp7, RproObp1, RproObp3, RproObp11, RproObp20, and RproObp26, were detected in the brain of R. prolixus. The presence of these carrier proteins was previously reported in the brain of other insects [76–78]; however, its functional role in this tissue is still unknown. As proposed by Walker et al. (2019), these carriers could be monitoring internal chemical signals in the R. prolixus brain [76].

Drosophila melanogaster takeout 1 (to1) gene (DmelTo1) has been related to the regulation of feeding behavior and locomotor activity, and its expression has been detected in various fly structures and tissues, including the head, fat body, crop, and antennae [79, 80]. DmelTo1 also affects male courtship behavior [81]. The role of these proteins has been poorly studied in insects other than Dipterans, even though to date the scarce evidence also points to behavioral roles in a locust and a moth [82, 83]. Regarding triatomines, the expression of TO genes has been reported in the digestive tract [29] and antennae of R. prolixus [26] and T. brasiliensis [84]. In our brain transcriptome, RproTo1, RproTo2, RproTo4, and RproTo6 were highly expressed (S4 Fig), suggesting that they may have relevant behavioral functions at the central level, as observed for D. melanogaster.

Final remarks

Expression datasets obtained using transcriptomes represent powerful tools to uncover molecular targets for functional studies. Furthermore, these data allow improving automatically predicted gene models of interest through the manual curation of sequences [21, 22]. Therefore, they also increase the chances of performing successful functional experiments based on more trustable gene models. A drawback associated with whole tissue transcriptomes is that complex structures like the brain can present an intricate organization with specialized areas having very differentiated functional roles, and consequently, specific gene expression profiles. Therefore, brain transcriptome studies should acknowledge that expression profiles represent averages of neuropiles having differentiated properties. This is especially true for clock genes or NPGs which can be expressed in very restricted sets of neurons. Therefore, any lack of differential expression observed in studies comparing levels of expression in different developmental or physiological conditions should be later validated with tissue-specific or single-cell sequencing methods, when available. Still, the high levels of expression on which we decided to focus here seem to denote relevant functions that deserve attention, as they might guide research toward specific targets allowing the development of more rational control methods.

Supporting information

S1 Fig. Principal component analysis graph of the six libraries.

Rep 1, Rep 2, Rep 3 (experiment 1) and Rep 4, Rep 5, Rep 6 (experiment 2).

https://doi.org/10.1371/journal.pone.0282490.s001

(TIF)

S2 Fig. Expression profiles of neuropeptide processing enzymes genes.

Heatmap depicting the expression level of neuropeptide processing enzymes genes in the brain of R. prolixus nymphs. Expression level (displayed as Log₁₀ TPM+1) is represented by means of a color scale in which blue/red represent the lowest/highest expression. Each column represents the expression of one library.

https://doi.org/10.1371/journal.pone.0282490.s002

(TIF)

S3 Fig. Expression profiles of chemosensory proteins and odorant binding proteins.

(A) Heatmap depicting the expression level of chemosensory proteins (CSPs), and (B) odorant binding protein (OBPs) genes in the brain of R. prolixus nymphs. Expression level (displayed as Log₁₀ TPM+1) is represented by means of a color scale in which blue/red represent the lowest/highest expression. Each column represents the expression of one library.

https://doi.org/10.1371/journal.pone.0282490.s003

(TIF)

S4 Fig. Expression profiles of takeout genes.

Heatmap depicting the expression level of takeout genes in the brain of R. prolixus nymphs. Expression level (displayed as Log₁₀ TPM+1) is represented by means of a color scale in which blue/red represent the lowest/highest expression. Each column represents the expression of one library.

https://doi.org/10.1371/journal.pone.0282490.s004

(TIF)

S1 Table. Details of the mRNA expression of Figs 1–4, S2–S4 Figs.

Columns are: the abbreviation of the gene assigned; gene name according to the annotation; VectorBase code–the official gene number in the RproC3 genome assembly; values of TPM in each library (replicate); values of (Log₁₀ TPM+1) in each library. ND: not determined.

https://doi.org/10.1371/journal.pone.0282490.s005

(XLSX)

Citation: Nevoa JC, Latorre-Estivalis JM, Pais FS-M, Marliére NP, Fernandes GdR, Lorenzo MG, et al. (2023) Global characterization of gene expression in the brain of starved immature Rhodnius prolixus. PLoS ONE 18(3): e0282490. https://doi.org/10.1371/journal.pone.0282490

About the Authors:

Jessica Coraiola Nevoa

Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Vector Behaviour and Pathogen Interaction Group, Instituto René Rachou – FIOCRUZ, Belo Horizonte, Minas Gerais, Brazil

Jose Manuel Latorre-Estivalis

Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Writing – original draft, Writing – review & editing

Affiliation: Laboratorio de Insectos Sociales, Instituto de Fisiología, Biología Molecular y Neurociencias, Universidad de Buenos Aires - CONICET, Buenos Aires, Argentina

Fabiano Sviatopolk-Mirsky Pais

Roles: Data curation, Formal analysis, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Bioinformatics Plataform, Instituto René Rachou - FIOCRUZ, Belo Horizonte, Minas Gerais, Brazil

Newmar Pinto Marliére

Roles: Conceptualization, Methodology

Affiliation: Vector Behaviour and Pathogen Interaction Group, Instituto René Rachou – FIOCRUZ, Belo Horizonte, Minas Gerais, Brazil

Gabriel da Rocha Fernandes

Roles: Data curation, Formal analysis, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Bioinformatics Plataform, Instituto René Rachou - FIOCRUZ, Belo Horizonte, Minas Gerais, Brazil

Marcelo Gustavo Lorenzo

Roles: Conceptualization, Funding acquisition, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Vector Behaviour and Pathogen Interaction Group, Instituto René Rachou – FIOCRUZ, Belo Horizonte, Minas Gerais, Brazil

Alessandra Aparecida Guarneri

Roles: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing

E-mail: [email protected]

Affiliation: Vector Behaviour and Pathogen Interaction Group, Instituto René Rachou – FIOCRUZ, Belo Horizonte, Minas Gerais, Brazil

ORICD: https://orcid.org/0000-0002-9162-2731

References

1. Peterson JK, Bartsch SM, Lee BY, Dobson AP. Broad patterns in domestic vector-borne Trypanosoma cruzi transmission dynamics: synanthropic animals and vector control Quantitative analysis of strategies to achieve the 2020 goals for neglected tropical diseases: where are we now? Parasites and Vectors. 2015;8(1):1–1 Available from: https://parasitesandvectors.biomedcentral.com/articles/10.1186/s13071-015-1146-1

2. Lazzari CR. The Behaviour of Kissing Bugs. In: Guarneri AA, Lorenzo MG, editors. Triatominae—The Biology of Chagas Disease Vectors. Springer Nature. 2021. p. 215 https://link.springer.com/10.1007/978-3-030-64548-9_10

3. Ferreira R, Guarneri AA, Lorenzo MG. Activity and shelter-related behavior in Rhodnius prolixus: The role of host odours. Acta Trop. 2019;196:150–4. Available from: https://pubmed.ncbi.nlm.nih.gov/31116968/

4. Marliére NP, Latorre-Estivalis JM, Lorenzo MG, Carrasco D, Alves-Silva J, Rodrigues J de O, et al. Trypanosomes modify the behavior of their insect hosts: effects on locomotion and on the expression of a related gene. PLoS Negl Trop Dis. 2015; 9(8). pmid:26291723

5. Marliére NP, Lorenzo MG, Guarneri AA. Trypanosoma cruzi-infected Rhodnius prolixus endure increased predation facilitating parasite transmission to mammal hosts. PLoS Negl Trop Dis. 2021; 15(7). Available from: /pmc/articles/PMC8279422/

6. Friend WG, Smith JJ. Factors affecting feeding by bloodsucking insects. Annu Rev Entomol. 1977; 22:309–31. pmid:319741

7. Bodin A, Vinauger C, Lazzari CR. State-dependency of host-seeking in Rhodnius prolixus: The post-ecdysis time. J Insect Physiol. 2009;55(6).

8. Chapman RF. The insects: structure and function. 6th ed. Cambridge: Harvard University Press; 1982.

9. Giebultowicz JM. Molecular mechanism and cellular distribution of insect circadian clocks. Annu Rev Entomol. 2000;45: 769–93. pmid:10761596

10. Numata H, Miyazaki Y, Ikeno T. Common features in diverse insect clocks. Zool Lett. 2015;1(10):1–17. pmid:26605055

11. Saunders DS, Steel CGH, Vafopoulou X, Lewis RD. Insect Clocks. Elsevier. 2002. https://linkinghub.elsevier.com/retrieve/pii/B9780444504074X50009

12. Nässel DR. Neuropeptides in the nervous system of Drosophila and other insects: multiple roles as neuromodulators and neurohormones. Prog Neurobiol. 2002;68: 1–84.

13. Nässel DR. Functional roles of neuropeptides in the insect central nervous system. Naturwissenschaften. 2000;87: 439–49. pmid:11129943

14. Ons S. Neuropeptides in the regulation of Rhodnius prolixus physiology. J Insect Physiol. 2017;97: 77–92. Available from: http://dx.doi.org/10.1016/j.jinsphys.2016.05.003

15. Ons S, Lavore A, Sterkel M, Wulff JP, Sierra I, Martínez-Barnetche J, et al. Identification of G protein coupled receptors for opsines and neurohormones in Rhodnius prolixus. Genomic and transcriptomic analysis. Insect Biochem Mol Biol. 2016;69: 34–50.

16. Angulo VM, Esteban L. Nueva trampa para la captura de triatominos en hábitats silvestres y peridomésticos. Biomédica. 2011 Mar 7;31(2).

17. Angulo VM, Esteban L, Luna KP. Attalea butyracea próximas a las viviendas como posible fuente de infestación domiciliaria por Rhodnius prolixus (Hemiptera: Reduviidae) en los Llanos Orientales de Colombia. Biomédica. 2012 Jun 30;32(2).

18. Esteban L, Montes JM, Angulo VM. Diversity of Triatominae (Hemiptera: Reduviidae) in Santander, Colombia: Epidemiological implications. Biomédica. 2017 Jan 24;37(1). pmid:28527247

19. Urbano P, Hincapié E, Angulo VM, Esteban L, Núñez-Avellaneda LA. Variación poblacional de Rhodnius prolixus (Reduviidae: Triatominae) en Attalea butyracea (Arecaceae) en la Orinoquia colombiana. Rev Colomb Entomol. 2018 Dec 31;44(2).

20. Nunes-da-Fonseca R, Berni M, Tobias-Santos V, Pane A, Araujo HM. Rhodnius prolixus: From classical physiology to modern developmental biology. Genesis. 2017 May 1;55(5).

21. Mesquita RD, Vionette-Amaral RJ, Lowenberger C, Rivera-Pomar R, Monteiro FA, Minx P, et al. Genome of Rhodnius prolixus, an insect vector of Chagas disease, reveals unique adaptations to hematophagy and parasite infection. Proc Natl Acad Sci U S A. 2015;112(48).

22. Leyria J, Orchard I, Lange AB. Transcriptomic analysis of regulatory pathways involved in female reproductive physiology of Rhodnius prolixus under different nutritional states. Sci Rep. 2020;10(1).

23. Capriotti N, Ianowski JP, Gioino P, Ons S. The neuropeptide CCHamide2 regulates diuresis in the Chagas disease vector Rhodnius prolixus. J Exp Biol. 2019;222(10).

24. Latorre-Estivalis JM, de Oliveira ES, Beiral Esteves B, Santos Guimarães L, Neves Ramos M, Lorenzo MG. Patterns of expression of odorant receptor genes in a Chagas disease vector. Insect Biochem Mol Biol. 2016;69: 71–81. pmid:26003917

25. Latorre-Estivalis JM, Robertson HM, Walden KKO, Ruiz J, Gonçalves LO, Guarneri AA, et al. The molecular sensory machinery of a Chagas disease vector: Expression changes through imaginal moult and sexually dimorphic features. Sci Rep. 2017;7. pmid:28059141

26. Latorre-Estivalis JM, Sterkel M, Ons S, Lorenzo MG. Transcriptomics supports local sensory regulation in the antenna of the kissing-bug Rhodnius prolixus. BMC Genomics. 2020;21(1).

27. Ons S, Richter F, Urlaub H, Pomar RR. The neuropeptidome of Rhodnius prolixus brain. Proteomics. 2009;9(3): 788–92.

28. Ons S, Sterkel M, Diambra L, Urlaub H, Rivera-Pomar R. Neuropeptide precursor gene discovery in the Chagas disease vector Rhodnius prolixus. Insect Mol Biol. 2011;20(1):29–44.

29. Ribeiro JMC, Genta FA, Sorgine MHF, Logullo R, Mesquita RD, Paiva-Silva GO, et al. An insight into the transcriptome of the digestive tract of the bloodsucking bug, Rhodnius prolixus. PLoS Negl Trop Dis. 2014 Jan 9;8(1).

30. Sterkel M, Urlaub H, Rivera-Pomar R, Ons S. Functional proteomics of neuropeptidome dynamics during the feeding process of Rhodnius prolixus. J Proteome Res. 2011 Aug 5;10(8):3363–71.

31. Ribeiro JMC, Andersen J, Silva-Neto MAC, Pham VM, Garfield MK, Valenzuela JG. Exploring the sialome of the blood-sucking bug Rhodnius prolixus. Insect Biochem Mol Biol. 2004 Jan 1;34(1):61–79.

32. Medeiros MN, Logullo R, Ramos IB, Sorgine MHF, Paiva-Silva GO, Mesquita RD, et al. Transcriptome and gene expression profile of ovarian follicle tissue of the triatomine bug Rhodnius prolixus. Insect Biochem Mol Biol [Internet]. 2011 Oct;41(10):823–31. Available from: https://linkinghub.elsevier.com/retrieve/pii/S0965174811001342

33. Coelho VL, De Brito TF, Brito IADA, Cardoso MA, Berni MA, Maria H, et al. Analysis of ovarian transcriptomes reveals thousands of novel genes in the insect vector Rhodnius prolixus. Sci Rep [Internet]. 2021;11:1918. Available from: https://doi.org/10.1038/s41598-021-81387-1

34. Cardoso JC, Ribeiro JMC, Santos DV dos, Pereira MH, Araújo RN, Gontijo NF, et al. Analysis of the testicle’s transcriptome of the Chagas disease vector Rhodnius prolixus. bioRxiv [Internet]. 2020 Apr 13 [cited 2022 May 4];616193. https://www.biorxiv.org/content/10.1101/616193v2

35. Salcedo-Porras N, Guarneri A, Oliveira PL, Lowenberger C. Rhodnius prolixus: Identification of missing components of the IMD immune signaling pathway and functional characterization of its role in eliminating bacteria. PLoS One 2019 Apr;14(4). Available from: https://dx.plos.org/10.1371/journal.pone.0214794

36. Latorre-Estivalis JM, Große-Wilde E, da Rocha Fernandes G, Hansson BS, Lorenzo MG. Changes in antennal gene expression underlying sensory system maturation in Rhodnius prolixus. Insect Biochem Mol Biol. 2022 Jan;140:103704.

37. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 2014 Aug 1;30(15):2114–20. pmid:24695404

38. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013 Jan 1 [cited 2021 Nov 1];29(1):15–21. Available from: https://academic.oup.com/bioinformatics/article/29/1/15/272537 pmid:23104886

39. Giraldo-Calderón G, Emrich S, MacCallum R, 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 Jan 28 [cited 2021 Nov 1];43(Database issue):D707–13. Available from: https://pubmed.ncbi.nlm.nih.gov/25510499/ pmid:25510499

40. Hwang JR, Siekhaus DE, Fuller RS, Taghert PH, Lindberg I. Interaction of Drosophila melanogaster Prohormone Convertase 2 and 7B2. J Biol Chem. 2000 Jun;275(23):17886–93. Available from: https://linkinghub.elsevier.com/retrieve/pii/S0021925819833215

41. Caggese C, Barsanti P, Viggiano L, Bozzetti MP, Caizzi R. Genetic, molecular and developmental analysis of the glutamine synthetase isozymes of Drosophila melanogaster. Genetica. 1994 Jun;94(2–3):275–81. Available from: http://link.springer.com/10.1007/BF01443441

42. Christ P, Reifenrath A, Kahnt J, Hauser F, Hill SR, Schachtner J, et al. Feeding-induced changes in allatostatin-A and short neuropeptide F in the antennal lobes affect odor-mediated host seeking in the yellow fever mosquito, Aedes aegypti. PLoS One. 2017 Nov 1;12(11): e0188243. Available from: https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0188243

43. Haddad ANS, Defferrari MS, Hana S, Szeto SG, Lange AB. Expression and functional characterization of tachykinin-related peptides in the blood-feeding bug, Rhodnius prolixus. Peptides. 2018 Jan 1;99:247–54.

44. Nässel DR, Zandawala M. Recent advances in neuropeptide signaling in Drosophila, from genes to physiology and behavior. Prog Neurobiol. 2019 Aug;179:101607. Available from: https://linkinghub.elsevier.com/retrieve/pii/S0301008218302144

45. Claeys I, Simonet G, Poels J, Van Loy T, Vercammen L, De Loof A, et al. Insulin-related peptides and their conserved signal transduction pathway. Peptides. 2002 Apr;23(4). pmid:11897402

46. Nässel DR, Kubrak OI, Liu Y, Luo J, Lushchak O V. Factors that regulate insulin producing cells and their output in Drosophila. Front Physiol. 2013;4.

47. Vafopoulou X, Steel CGH. Insulin-like and testis ecdysiotropin neuropeptides are regulated by the circadian timing system in the brain during larval-adult development in the insect Rhodnius prolixus (Hemiptera). Gen Comp Endocrinol. 2012 Nov 1;179(2):277–88.

48. Lindberg I, Zhou Y. Overexpression of neuropeptide precursors and processing enzymes. Methods Neurosci. 1995;23:94–108. Available from: https://linkinghub.elsevier.com/retrieve/pii/S1043947106801163

49. Sossin W, Fisher J, Scheller R. Cellular and molecular biology of neuropeptide processing and packaging. Neuron. 1989;2:1407–17. pmid:2697236

50. Nässel DR, Homberg U. Neuropeptides in interneurons of the insect brain. Cell Tissue Res. 2006 Oct 8;326(1):1–24. Available from: http://link.springer.com/10.1007/s00441-006-0210-8 pmid:16761145

51. Pauls D, Chen J, Reiher W, Vanselow JT, Schlosser A, Kahnt J, et al. Peptidomics and processing of regulatory peptides in the fruit fly Drosophila. EuPA Open Proteomics. 2014 Jun;3:114–27. Available from: https://linkinghub.elsevier.com/retrieve/pii/S2212968514000117

52. Villalobos-Sambucaro MJ, Diambra LA, Noriega FG, Ronderos JR. Allatostatin-C antagonizes the synergistic myostimulatory effect of allatotropin and serotonin in Rhodnius prolixus (Stal). Gen Comp Endocrinol. 2016 Jul 1;233:1–7. Available from: https://pubmed.ncbi.nlm.nih.gov/27174746/

53. Zandawala M, Li S, Hauser F, Grimmelikhuijzen CJP, Orchard I. Isolation and functional characterization of calcitonin-like diuretic hormone receptors in Rhodnius prolixus. PLoS One. 2013 Nov 29;8(11).

54. Mcbain CJ, Mayer ML. N-methyl-D-aspartic acid receptor structure and function. Physiol Rev. 1994;74(3). pmid:8036251

55. Ultsch A, Schuster CM, Laube B, Betz H, Schmitt B. Glutamate receptors of Drosophila melanogaster Primary structure of a putative NMDA receptor protein expressed in the head of the adult fly. Fed Eur Biochem Soc. 1993;324(2).

56. Xia S, Chiang AS. NMDA Receptors in Drosophila. Eur PMC. 2009.

57. Huang J, Hult EF, Marchal E, Tobe SS. Identification and characterization of the NMDA receptor and its role in regulating reproduction in the cockroach Diploptera punctata. J Exp Biol. 2015 Apr 1;218(7):983–90.

58. Xia S, Miyashita T, Fu TF, Lin WY, Wu CL, Pyzocha L, et al. NMDA receptors mediate olfactory learning and memory in Drosophila. Curr Biol. 2005 Apr 12;15(7):603–15.

59. Jaumouillé E, Machado Almeida P, Stähli P, Koch R, Nagoshi E. Transcriptional regulation via nuclear receptor crosstalk required for the drosophila circadian clock. Curr Biol. 2015 Jun 1;25(11):1502–8. Available from: http://dx.doi.org/10.1016/j.cub.2015.04.017 pmid:26004759

60. Beuchle D, Jaumouillé E, Nagoshi E. The nuclear receptor unfulfilled is required for free-running clocks in Drosophila pacemaker neurons. Curr Biol. 2012 Jul 10;22(13):1221–7.

61. Kumar S, Chen D, Jang C, Nall A, Zheng X, Sehgal A. An ecdysone-responsive nuclear receptor regulates circadian rhythms in Drosophila. Nat Commun. 2014; Available from: www.nature.com/naturecommunications

62. Crane BR, Young MW. Interactive features of proteins composing eukaryotic circadian clocks. Annu Rev Biochem. 2014;83:191–219. pmid:24905781

63. Lee J, Yoo E, Lee H, Park K, Hur J-H, Lim C. LSM12 and ME31B/DDX6 define distinct modes of posttranscriptional regulation by ATAXIN-2 protein complex in Drosophila circadian pacemaker neurons. Mol Cell. 2017 Apr;66(1):129–140.e7. Available from: https://linkinghub.elsevier.com/retrieve/pii/S1097276517301715

64. Sathyanarayanan S, Zheng X, Xiao R, Sehgal A. Posttranslational regulation of Drosophila PERIOD protein by protein phosphatase 2A. Cell. 2004 Feb;116(4):603–15. Available from: https://linkinghub.elsevier.com/retrieve/pii/S009286740400128X

65. Allen AM, Anreiter I, Neville MC, Sokolowski MB. Feeding-related traits are affected by dosage of the foraging gene in Drosophila melanogaster. Genetics. 2017 Feb 1;205(2):761–73.

66. Anreiter I, Sokolowski MB. The foraging gene and its behavioral effects: Pleiotropy and plasticity. Annu Rev Genet. 2019;53:10–20. Available from: https://doi.org/10.1146/annurev-genet-112618- pmid:31487469

67. Osborne KA. Natural behavior polymorphism due to a cGMP-Dependent protein kinase of Drosophila. Science (80-). 1997 Aug 8;277(5327).

68. Lucas C, Kornfein R, Chakaborty-Chatterjee M, Schonfeld J, Geva N, Sokolowski MB, et al. The locust foraging gene. Arch Insect Biochem Physiol. 2010 May;74(1). pmid:20422718

69. Ingram KK, Oefner P, Gordon DM. Task-specific expression of the foraging gene in harvester ants. Mol Ecol. 2005 Feb 7;14(3). pmid:15723672

70. Ben-Shahar Y. Influence of Gene Action Across Different Time Scales on Behavior. Science (80-). 2002 Apr 26;296(5568). pmid:11976457

71. Tobback J, Mommaerts V, Vandersmissen HP, Smagghe G, Huybrechts R. Age- and task-dependent foraging gene expression in the bumblebee Bombus terrestris. Arch Insect Biochem Physiol. 2011 Jan;76(1).

72. Marliére NP, Lorenzo MG, Martínez Villegas LE, Guarneri AA. Co-existing locomotory activity and gene expression profiles in a kissing-bug vector of Chagas disease. J Insect Physiol. 2020 Apr 1;122.

73. Rihani K, Ferveur J-F, Briand L. The 40-year mystery of insect odorant-binding proteins. Biomolecules. 2021 Mar 30;11(4):509. Available from: https://www.mdpi.com/2218-273X/11/4/509 pmid:33808208

74. Wang Q, Zhou J-J, Liu J-T, Huang G-Z, Xu W-Y, Zhang Q, et al. Integrative transcriptomic and genomic analysis of odorant binding proteins and chemosensory proteins in aphids. Insect Mol Biol. 2019 Feb;28(1):1–22. Available from: https://onlinelibrary.wiley.com/doi/10.1111/imb.12513 pmid:29888835

75. Pelosi P, Iovinella I, Zhu J, Wang G, Dani FR. Beyond chemoreception: diverse tasks of soluble olfactory proteins in insects. Biol Rev [Internet]. 2018 [cited 2023 Feb 3];93:184–200. Available from: https://onlinelibrary.wiley.com/doi/10.1111/brv.12339 pmid:28480618

76. Walker WB, Roy A, Anderson P, Schlyter F, Hansson BS, Larsson MC. Transcriptome Analysis of Gene Families Involved in Chemosensory Function in Spodoptera littoralis (Lepidoptera: Noctuidae). BMC Genomics [Internet]. 2019 Dec 28;20(1):428. Available from: https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-019-5815-x pmid:31138111

77. Chen W-B, Du L-X, Gao X-Y, Sun L-L, Chen L-L, Xie G-Y, et al. Identification of Odorant-Binding and Chemosensory Protein Genes in Mythimna separata Adult Brains Using Transcriptome Analyses. Front Physiol [Internet]. 2022 Feb 28;13. Available from: https://www.frontiersin.org/articles/10.3389/fphys.2022.839559/full pmid:35295575

78. Bian L, Li Z-Q, Ma L, Cai X-M, Luo Z-X, Chen Z-M. Identification of the genes in tea leafhopper, Empoasca onukii (Hemiptera: Cicadellidae), that encode odorant-binding proteins and chemosensory proteins using transcriptome analyses of insect heads. Appl Entomol Zool [Internet]. 2018 Feb 22;53(1):93–105. Available from: http://link.springer.com/10.1007/s13355-017-0533-9

79. Meunier N, Belgacem YH, Martin JR. Regulation of feeding behaviour and locomotor activity by takeout in Drosophila. J Exp Biol. 2007 Apr;210(8):1424–34. pmid:17401125

80. Sarov-Blat L, So WV, Liu L, Rosbash M. The Drosophila takeout gene is a novel molecular link between circadian rhythms and feeding behavior. Cell. 2000;101(6).

81. Dauwalder B, Tsujimoto S, Moss J, Mattox W. The Drosophila takeout gene is regulated by the somatic sex-determination pathway and affects male courtship behavior. Genes Dev. 2002 Nov 15;16(22):2879–92.

82. Guo W, Wang X, Ma Z, Xue L, Han J, Yu D, et al. CSP and takeout genes modulate the switch between attraction and repulsion during behavioral phase change in the migratory locust. Copenhaver GP, editor. PLoS Genet. 2011 Feb 3;7(2):e1001291. Available from: https://dx.plos.org/10.1371/journal.pgen.1001291 pmid:21304893

83. Saito K, Su ZH, Emi A, Mita K, Takeda M, Fujiwara Y. Cloning and expression analysis of takeout/JHBP family genes of silkworm, Bombyx mori. Insect Mol Biol. 2006 Jun;15(3):245–51. pmid:16756543

84. Marchant A, Mougel F, Jacquin-Joly E, Costa J, Almeida CE, Harry M. Under-expression of chemosensory genes in domiciliary bugs of the Chagas disease vector Triatoma brasiliensis. PLoS Negl Trop Dis [Internet]. 2016 Oct 28 [cited 2022 Jun 1];10(10):e0005067. Available from: https://journals.plos.org/plosntds/article?id=10.1371/journal.pntd.0005067