Next Article in Journal
The hEag1 K+ Channel Inhibitor Astemizole Stimulates Ca2+ Deposition in SaOS-2 and MG-63 Osteosarcoma Cultures
Next Article in Special Issue
The Third Extracellular Loop of Mammalian Odorant Receptors Is Involved in Ligand Binding
Previous Article in Journal
Vitamin D as a Possible COVID-19 Prevention Strategy
Previous Article in Special Issue
Taste Receptors beyond Taste Buds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Candidate Chemosensory Gene Families by Head Transcriptomes Analysis in the Mexican Fruit Fly, Anastrepha ludens Loew (Diptera: Tephritidae)

by
Obdulia L. Segura-León
1,*,
Brenda Torres-Huerta
1,*,
Alan Rubén Estrada-Pérez
2,
Juan Cibrián-Tovar
1,
Fidel de la Cruz Hernandez-Hernandez
3,
José Luis Cruz-Jaramillo
4,
José Salvador Meza-Hernández
5 and
Fabian Sánchez-Galicia
6
1
Entomology and Acarology Program, Colegio de Postgraduados Campus Montecillo, Km. 36.5, Texcoco 56230, Mexico
2
Escuela Superior de Medicina (ESM), Intituto Politécnico Nacional, Casco de Santo Tomas, Miguel Hidalgo, Mexico City 11340, Mexico
3
Department of Infectomics and Molecular Pathogenesis, Center for Research and Advanced Studies (CINVESTAV-IPN), Mexico City 07360, Mexico
4
Bioinformatics and Technologies Department, Solaria Biodata, Col del Valle Centro, Benito Juarez, Mexico City 03100, Mexico
5
Programa Operativo de Moscas, SENASICA-SADER/IICA, Metapa de Dominguez 30860, Mexico
6
Servicio Nacional de Sanidad Inocuidad y Calidad Agroalimentaria, Secretaría de Agricultura y Desarrollo Rural, Cuauhtemoc, Mexico City 06100, Mexico
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(18), 10531; https://doi.org/10.3390/ijms231810531
Submission received: 6 August 2022 / Revised: 31 August 2022 / Accepted: 6 September 2022 / Published: 11 September 2022
(This article belongs to the Special Issue Molecular and Cellular Mechanisms Underlying Taste, Smell and Beyond)

Abstract

:
Insect chemosensory systems, such as smell and taste, are mediated by chemosensory receptor and non-receptor protein families. In the last decade, many studies have focused on discovering these families in Tephritidae species of agricultural importance. However, to date, there is no information on the Mexican fruit fly Anastrepha ludens Loew, a priority pest of quarantine importance in Mexico and other countries. This work represents the first effort to identify, classify and characterize the six chemosensory gene families by analyzing two head transcriptomes of sexually immature and mature adults of A. ludens from laboratory-reared and wild populations, respectively. We identified 120 chemosensory genes encoding 31 Odorant-Binding Proteins (OBPs), 5 Chemosensory Proteins (CSPs), 2 Sensory Neuron Membrane Proteins (SNMPs), 42 Odorant Receptors (ORs), 17 Ionotropic Receptors (IRs), and 23 Gustatory Receptors (GRs). The 120 described chemosensory proteins of the Mexican fruit fly significantly contribute to the genetic databases of insects, particularly dipterans. Except for some OBPs, this work reports for the first time the repertoire of olfactory proteins for one species of the genus Anastrepha, which provides a further basis for studying the olfactory system in the family Tephritidae, one of the most important for its economic and social impact worldwide.

1. Introduction

The chemical sensory systems such as smell and taste are part of an insect’s extraordinary capacity to find food, hosts, sexual partners, oviposition sites, avoid natural enemies, and other dangers [1]. The discrimination and interpretation of chemical information in the environment begins in the insect’s olfactory and gustatory organs, with a series of nerve impulses leading to a change in the insect behavior. Before reaching olfactory receptor neurons, there are different biochemical interactions called perireceptor events, that involve the transport of the odorant molecules across the sensillar lymph and the activation of specific transmembrane receptors for the transduction of the chemical signal [2].
This robust and sophisticated system includes several gene families encoding chemosensory proteins involved in complex biochemical reactions in various insect tissues [3,4]. Extracellular Odorant-Binding Proteins (OBPs) and Chemosensory Proteins (CSPs) are the first filters in olfactory processing; they transport hydrophobic compounds from the external environment to specific transmembrane receptors. Odorant Receptors (OR), Ionotropic Receptors (IR), and Gustatory Receptors (GR) initiate downstream signaling related to behavioral responses. Finally, the Sensory Neuron Membrane Proteins (SNMP) play a role in mediating ligand-OR interactions [5,6,7].
The demands of modern agriculture have generated several different challenges, such as creating specific and environmentally friendly alternatives for integrated pest management that reduce or displace the application of insecticides in the future [8]. In this sense, the study of chemosensory genes has increased in insects of agricultural importance [9,10]. Understanding the olfactory system of pests and its role in processes such as host colonization and mating are the basis for the potential development of biotechnological tools such as the creation of biosensors, the computational design of behaviorally active chemicals, and their manipulation by gene silencing techniques [11,12,13,14,15,16,17].
Phytophagous flies of the family Tephritidae, commonly called fruit flies, are among the most economically significant insect groups, causing primary fruit and horticultural losses worldwide in tropical and subtropical regions [18,19]. In the Americas, the largest genus of Tephritidae is Anastrepha Schiner, with more than 250 species distributed from the southern United States to northern Argentina [20,21]. In Mexico, A. ludens, or the “Mexican fruit fly”, is a pest of quarantine importance that represents a major phytosanitary problem in various fruits, particularly citrus and mango [22,23,24]. It caused direct damage due to crop yield losses, increased control costs, and international marketing restrictions [25].
In the last decade, an increasing number of studies have focused on discovering the chemosensory genes of pest tephritid species such as Bactrocera dorsalis, B. minax, Zeugodacus cucurbitae, Z. tau, Ceratitis capitata, A. fraterculus, and A. obliqua, to understand the olfactory pathways and their role in fruit flies’ behavior [26,27,28,29,30,31,32]. To date, there is no information on the multigene families of the olfactory system of A. ludens [33,34], but new sequencing technologies and advances in model insect communication systems offer the opportunity to understand similar mechanisms in the Mexican fruit fly.
Identifying putative chemosensory genes in A. ludens is the first step to exploring crucial gene functions in the communication process. Their knowledge can provide further insight into the evolutionary history of the genus Anastrepha, related to the search behavior, host selection, and mating in recently diverged species with a limited number of distinctive morphological and genetic characters [27,35]. In this study, we identified, analyzed, and characterized the six chemosensory gene families of A. ludens by analyzing two head transcriptomes of sexually immature and mature adults from laboratory-reared and wild populations, respectively.

2. Results

2.1. Transcriptome Assembly and Annotation

As a result of sequencing five Illumina HiSeq2000 libraries of sexually immature male and female heads of A. ludens, we obtained an average of 6.82 million clean reads, 92.37% with a Phred quality score Q20. The clean reads were assembled by de novo into 24,368 transcripts, and after removing redundancies and isoforms, 21,103 unigenes were obtained with a mean length of 359 bp and an N50 of 880 (Supplementary Table S1). On the other hand, with the six paired libraries generated with the BGISEQ-500 platform from sexually mature male heads of A. ludens, an average of 47.94 million clean reads were obtained with a length of 99 bp, and more than 97% of reads had a Phred quality score Q20 and low duplication rates. The de novo assembly resulted in 269,924 transcripts and after removing redundancies, we obtained 158,693 unigenes with a mean length of 358 bp and N50 of 1245 (Supplementary Table S2).
We used Gene Ontology (GO), UniProtKB, and InterPro databases (accessed on 1 October 2021) to annotate unigenes obtained from A. ludens head transcriptomes from Illumina HiSeq2000 and BGISEQ-500 platforms. A total of 13,723 (65.02%) unigenes obtained with Illumina (Figure 1a) and 39,696 (25.01%) unigenes generated with BGI were recorded in at least one database (Figure 1). In both cases, of the total number of annotated unigenes, more than 80% were recorded in the Insecta-UniprotKB, more than 50% in the InterPro, and less than 40% in the GO databases (accessed on 1 October 2021). In both A. ludens transcriptomes, more than 90% of the genes had hits with sequences from the family Tephritidae, mostly with sequences from B. dorsalis, B. tryoni, C. capitata, B. latifrons, and Z. cucurbitae, and only 1.05% matched sequences from seven species of the genus Anastrepha. In addition, more than 50% of the genes had E-value lower than 1 × 10−50 and similarities higher than 70% (Figure 1).
Transcriptomes generated with Illumina HiSeq2000 and BGISEQ-500 were annotated in the three main GO categories. The molecular function category was the most representative, and genes enriched with more than 50% and 30% were binding and catalytic activity. Cellular processes and metabolic terms were the most represented subcategories in biological processes. In contrast, we classified more than 50% of unigenes in a cellular component in the subcategories of cell part, cell, and membrane (Supplementary Figure S1).

2.2. Mining of Chemosensory Genes

We mined the protein-coding transcripts with sequence similarities to the six chemosensory families using the Insecta base of UniProtKB and identified a total of 224 transcripts homologous to communication genes; 95 belong to non-receptor families, and 132 to receptor families. We obtained the open reading frames (ORFs) of the 224 transcripts and performed Blastp analysis with protein domain searches as retention criteria to maximize sensitivity and obtain functionally significant sequences. After functional analysis and elimination of redundant sequences, the number of filtered genes was reduced to 120 transcripts encoding 31 AludOBPs, 5 AludCSPs, 2 AludSNMPs, 42 AludORs, 23 AludGRs, and 17 AludIRs.

2.3. Odorant Binding Proteins (OBPs)

We identified thirty-one candidate OBP transcripts encoding proteins ranging in length from 102 to 306 amino acids. Of the 31 sequences, twenty-one had complete ORFs, seven partial ORFs with truncation in the 5′ region, three internal partials, and 19 AludOBPs had signal peptides. All proteins were functionally annotated by searching for domains within the insect pheromone superfamily/odor-binding proteins (SSF47565) and PBP/GOBP family (PF01395). All 31 OBPs had homologs with odor proteins of the family Tephritidae, 19 AludOBPs had similarities to A. obliqua and A. fraterculus with values mostly higher than 90%, and sequences with similarities to OBPs from B. dorsalis, C. capitata, and Z. tau were also recorded (Supplementary Table S3).
The classification of the AludOBPs was performed based on the number and location of the conserved cysteines; we grouped 31 proteins into classical, Minus-C, Plus-C, and dimer subfamilies (Table 1). The subfamily with the highest number of AludOBPs (17) was classical, which presented the general pattern of the six conserved cysteines and some with additional cysteines upstream of C1 (Supplementary Figure S2). We also identified eight AludOBPs Minus-C, which only conserved four cysteines and had a different number of residues between them. AludOBP49a1, AludOBP49a2, AludOBP50a, and AludOBP50e were classified within the Plus-C subfamily, which presented additional cysteines at the N-terminal end, two or three additional cysteines downstream of C6 and a conserved proline after the seventh cysteine. AludOBP83ef and AludOBP83cd presented two six-cysteine motifs classified within the dimeric subfamily (Supplementary Figure S3).
On the other hand, we constructed a phylogenetic tree by Bayesian Inference of candidate OBPs of A. ludens with proteins from seven Tephritidae species and the model organism Drosophila melanogaster (Figure 2). Most of the AludOBPs clustered into different subclades with their closest homologs from A. fraterculus and A. obliqua with 100% Bayesian posterior probabilities (BPP), related to proteins from other fruit flies species and D. melanogaster, showing highly supported terminal clusters of proteins with specific lineage expansions.
Four groups representing the subfamilies identified in the AludOBPs were defined in the phylogenetic tree; the Minus-C and dimer members were more closely related and clustered within the same clade with high supports (>90% BPP). Most Minus-C in Tephritidae are classified as OBP99a; however, in our analysis, several of them are related to OBPs from D. melanogaster OBP44a, OBP83g, OBP99b, OBP8a, and OBP99c, suggesting a classification that has reduced the representation of members of this subfamily and their diversification into other fruit fly genera.
The analysis grouped the classical family into three major clades; the group containing the proteins AludOBP73a and AludOBP59a was related to members of the Plus-C subfamily, with a higher number of residues (≈171) between C1 and C2. On the other hand, the phylogeny clustered two genes encoding OBP19a (AludOBP19a1 and AludOBP19a2) with DmelOBP19a and the OBP56h were divided into two groups within the same clade, each related to different DmelOBPs (100% BPP). Likewise, within the classical AludOBPs, we identified proteins related to D. melanogaster proteins as AludOBP83a and B. dorsalis and Z. tau proteins as AludOBP73a and AludOBP83a (Figure 2).

2.4. Chemosensory Proteins (CSPs)

We identified five candidate CSP transcripts from A. ludens, four encode proteins with complete ORFs and lengths between 120 and 149 aa, while AludCSP3-1 had a partial ORF in the 3′ region and length of 51 aa. All AludCSPs had annotations by conserved domains within the chemosensory protein superfamily CSP2 (SSF100910) and insect odorant-binding domains such as A10/Ejaculate Bulb Specific Protein (PF03392), related to chemosensory proteins (PTHR11257). The five AludCSPs had homologs with odor proteins of fruit flies B. minax, B. dorsalis, and Z. tau with identities higher than 79% (Supplementary Table S4).
The predicted AludCSPs proteins have all four conserved cysteine residues (C1-X6-C2-X18-C3-X2-C4). The phylogenetic analysis grouped the AludCSPs in two major clades with 100% BPP. The five AludCSPs were dispersed into five well-defined subgroups related to their orthologues from tephritid species, and except for the clade of AludCSP3-2, all CSPs from tephritid were clustered with one or two DmelCSPs with BPP > 90% (Figure 3). Sixteen of the tephritid and three D. melanogaster sequences homologous to AludCSPs, are not characterized in the NCBI non-redundant database (accessed on 1 August 2022) but have identity percentages >80% with AludCSPs, and present the characteristic protein domains of the CSP family (A10/PEBIII/OSD).

2.5. Sensory Neuron Membrane Proteins (SNMPs)

Two SNMP transcripts were identified in A. ludens head transcriptomes encoding 509 and 630 residue proteins with complete ORF and two transmembrane regions at the C- and N-terminal ends (Supplementary Table S5). The AludSNMPs had homologs with SNMP1 and SNMP2 from C. capitata and B. dorsalis with an E-value of zero and identities higher than 90% (Supplementary Table S5). The annotations based on conserved domains correspond to the CD36 family (PF01130). In the phylogenetic analysis of AludSNMPs with proteins from seven dipteran species (Figure 4), SNMPs were separated into SNMP2 and SNMP1 with 100% BPP. The AludSNMPs were more closely related to C. capitata proteins within each group. Furthermore, SNMP1 split into two groups, and the one identified for A. ludens grouped within SNMP1b.

2.6. Odorant Receptors (ORs)

We identified 42 OR genes, 35-recorded complete ORFs encoding proteins with lengths above 357 residues, and seven genes had partial ORFs in the 5′ or 3′ regions; however, only three of them had lengths less than 300 amino acids. We annotated all AludORs within the insect odorant receptor family (PF0249; Pfam: PF02949) and odorant receptor subfamilies (PANTHER: PTHR11857). AludORs had five to seven transmembrane domains, except for AludOR2a, AludOR7a6, AludOR71a1, AludOR85d, and AludOR88a; however, only three of them have partial sequences (Supplementary Table S6).
We performed two homology analyses for receptor annotation with the NCBI and UniprotKB-Insecta non-redundant protein databases (accessed on 1 February 2022). In both cases, all AludORs had homologs with dipteran proteins. With the UniprotKB base, almost 90% of the AludORs matched with proteins from seven tephritid species with identities of 30–96% (Supplementary Table S6). In contrast, with the non-redundant protein base, more than 50% of the AludORs had homologs with putative odorant receptors from R. pomonella, B. oleae, C. capitata, and B. tryoni, with higher similarity percentages (>70%) (Supplementary Table S6). One of the AludOR proteins, AludORCO, shares 96.19% and 95.77% of identity with a co-receptor ORCO/OR83b from Procecidochares utilis and B. tryoni, respectively.
To gain insight into the evolutionary relationships of the 42 AludORs, we performed a phylogenetic analysis by Bayesian Inference using a dataset of 283 odorant receptor sequences from eight tephritid species and the model fly D. melanogaster (Figure 5). ORs clustered into three major clades with 90% BPP; AludORs dispersed into subgroups related to their orthologues from tephritid species, most closely related to Rhagoletis species. AludORCO is clustered in a highly conserved clade with other olfactory co-receptors (OR83b) closely related to OR83a.
AludORs were organized with homologous sequences into 36 orthogroups with lineage-specific expansions; 29 consisted of well-defined single clusters (no duplications), of which 26 were related to homologous D. melanogaster receptors with >90% BPP. Three gene clusters (OR67c1, OR42a, and OR43a) did not show clustering with the model fly but with genes from at least six tephritid species (Figure 5). The remaining clusters included specific duplications such as OR71a, OR69a, and OR67d and specific expansions such as OR7a and OR74a (Figure 5).

2.7. Gustatory Receptors (GRs)

Twenty-three transcripts similar to GR were identified, encoding proteins with an average length of 350 residues. Ten had complete ORFs with five to eight transmembrane domains, whereas three to seven transmembrane domains were recorded in partial transcripts (Supplementary Table S7). The AludGRs had homologs with UniprotKB-base dipteran sequences, 20 AludGRs with similarity to B. correcta, B. dorsalis, B. latifrons, C. capitata, and Procecidochares utilis with identities mostly higher than 70% and three with proteins from D. ananassae, D. melanogaster, and D. mojavensis with similarities less than 60%.
AludGRs had similarities to different taste receptors of B. oleae, B. tryoni, and R. pomonella, more than 60% with an E-value of zero and identity percentages higher than 55% (Supplementary Table S7). InterPro analysis classified A. ludens GRs into the invertebrate gustatory receptor (PANTHER: PTHR21143) and chemosensory receptor (Pfam: PF08395) families. The phylogenetic analysis clustered the 23 AludGRs into three clades clustered with homologous sequences from seven tephritid species and the model species D. melanogaster with 100% BPP, most closely related to R. pomonella GRs, and other groups to C. capitata and B. latifrons (Figure 6).

2.8. Ionotropic Receptors (IRs)

We identified seventeen iGluRs/IRs-related transcripts in A. ludens head transcriptomes, ten with complete ORFs encoding proteins of 626–952 residues and seven partial ORFs with lengths above 180 residues (Supplementary Table S8). All AludIRs had transmembrane domains and annotation within the ionotropic receptor family (PANTHER: PTHR42643) and ion channel ligand-gated L-glutamate and glycine binding site (Pfam: PF10613). AludIRs had homology with proteins from B. correcta, B. dorsalis, B. latifrons, C. capitata, Z. tau, and Z. cucurbitae with identities of higher than 60%, except AludIR7c and AludIR31a that had similarities of 40% and 49% with IRs from D. lebanonensis and Musca domestica (Supplementary Table S8).
On the other hand, homology analyses with the NCBI non-redundant protein database (accessed on 1 February 2022) yielded different results, AludiGluRs/IRs had homologs only with proteins from the Tephritidae family with an E-value of zero and similarities higher than 70%, except for AludIR31a, and in contrast to UniProtKB, most of those proteins were not annotated (Supplementary Table S8). Phylogenetic analysis using Bayesian Inference corroborates the annotation of AludiGluRs/IRs with previously described homologous receptors from six tephritid species and D. melanogaster (Figure 7).
Three AludIRs members of the divergent IR clade (IR7c, 56c, and 100a) were identified; the group of antennal IRs was the largest, with ten members. The analysis clustered AludIR40a, AludIR41a, AludIR76c, and AludIR92a with their orthologues within the divergent clade. In addition, we identified a specific group for the antennal IRs AludIR31a, AludIR75a, AludIR75a2, and AludIR75d, while the receptors AludIR21a and AludIR93a were related to the putative co-receptors IR25a/IR8a and the non-NMD iGluRs AludGluIID and AludKaiR1D.

3. Discussion

The study of molecular communication in Tephritidae has boomed in the last fifteen years, for the potential benefits that understanding the mechanisms of chemosensory perception offers in developing strategies for the management of these insects [32]. Currently, the UniprotKB protein database (accessed on 1 February 2022) contains 1558 chemosensory proteins belonging to thirteen species of seven tephritid genera. Still, only 3% correspond to sequences from two species of the genus Anastrepha (A. fraterculus and A. obliqua). This work is the first report on the assembly and analysis of two transcriptomes of adult females and males of the Mexican fruit fly A. ludens and identifying 120 putative genes encoding proteins of the six chemosensory receptors and non-receptor families.
We performed A. ludens head transcriptomes with two different sequencing platforms, the number of unigenes obtained was higher with BGISEQ-500, but in general, both transcriptomes had a low percentage of annotated transcripts with respect to the total number of assembled unigenes. The differences in the number of unigenes obtained and annotation rates may be due to the specimen collection conditions and the number of libraries assembled with each platform. Likewise, other work has suggested that the low percentage of transcripts annotated in insect tissues may be due to genes that are not homologous to those deposited in the databases, high sequence divergence, or partial transcripts with no match [36,37].
Even with these differences, in both transcriptomes, more than 80% of unigenes translations shared significant similarities with entries from the UniprotKB database (accessed on 1 November 2021) of the order Diptera, followed by InterPro (≈50%), and less than 40% of transcripts fell into all three categories of GO gene ontology terms. The most abundant GO functional groups in both transcriptomes showed similar frequency in binding and catalytic activity, cellular, and metabolic process. Comparable to that reported in antennal and head transcriptomes of other dipteran species and insect orders, the similarity in the annotations of the two data sets indicates some level of conservation in gene expression in sensory tissues [36,38,39,40].
Different olfactory-related genes have been explored in tephritid species as potential candidates for developing biotechnological applications for pest management [26,41,42,43,44,45,46,47]. OBPs are small water-soluble proteins that transport odors to receptors in olfactory neurons and play a fundamental role in natural and sexual selection [48,49]. There are more than 15,000 amino acid sequences of insect OBPs in the NCBI database (accessed on 1 February 2022); however, for the genus Anastrepha, there is only information on 47 OBPs [27,50,51].
We identified 31 A. ludens genes encoding OBPs, this is higher than reported in head and reproductive tissue transcriptomes of the sister species A. obliqua (24) and A. fraterculus (23) but similar to OBPs identified in B. dorsalis (49), B. papayae (35), B. correcta (34), Z. cucurbitae (33), Z. tau (33), and C. capitata (34) [31,52,53]. The differences in putative OBPs described for the three Anastrepha species may be due to the physiological stage at which the transcriptomes were performed, related to the abundance of the genes of interest and the coverage obtained with the sequencing platforms [27,31]. Although this study increases the number of OBPs described for this genus, further studies in different tissues may increase the number of OBPs of the three Anastrepha species.
Comparative genomic analysis of different hexapod species has shown a very dynamic evolution of OBPs, with a wide range of protein lengths with varying profiles of cysteine [54,55,56]. The classification of the AludOBPs was identical to that described for A. fraterculus and A. obliqua; the largest subfamily was the classical subfamily with 17 proteins, followed by Minus-C (8), Plus-C (4), and we also identified two members of the dimer subfamily [42]. Most AludOBPs clustered with the closest orthologues of Anastrepha, supported with D. melanogaster OBPs with BPP >90%. Based on the relationships of AludOBPs with D. melanogaster orthologues, we made a name adjustment to make the proteins comparable with other studies.
In contrast to Campanini and Brito (2016) [27], we did not observe a basal division into four clades corresponding to each subfamily; the phylogenetic analyses clustered the members of the classical subfamily in three clades, one of them related to the Plus-C subfamily. As in other studies, OBPs of the dimer and Minus-C subfamily seem to have had independent origins, as they clustered in an external clade with 100% BPP supports [49]. In contrast to A. fraterculus and A. obliqua, we identified six additional AludOBPs in the classical subfamily (AludOBP28a, AludOBP73a, AludOBP84a, AludOBP19a1, AludOBP83a, and AludOBP19d2), all related to orthologues of D. melanogaster and tephritid except AludOBP19b-like, which was related to the OBP19b group and shared a 37.04% similarity to AludOBP19b.
It is common in OBPs to find homologous copies of a gene in the same species derived by the duplication process [55,56]. Campanini et al. (2017) [51] mentioned that the OBP56h group in different tephritid presents a duplication event and reported positive selection between AfraOBP56h-1/AoblOBP56h-1 and AfraOBP56h-2/AoblOBP56h-2. Initially, we identified two homologous members to OBP56h in A. ludens OBPs; however, their identity was less than 20%, so we corroborated the annotation against the D. melanogaster reference protein and observed identities of both proteins with OBP56h and OBP56g. Likewise, phylogenetic analysis showed the separation of tephritid OBP56h into two clades, each related to DmelOBP56g and DmelOBP56h with 100% BPP, suggesting two closely related sister groups with specific diversification.
On the other hand, OBP73a, which was related to the OBP59a, was characterized by having the presence of a large number of residues (>100) between C1 and C2, which had not been described in other Anastrepha transcriptomes [49,57,58,59]. Interestingly both proteins clustered within the same clade with the Plus-C subfamily; however, their unusual lengths could generate branch attraction, which may not reveal a true evolutionary history [60]. As in A. obliqua, we identified two copies of the Plus-C OBP49a proteins; the gene encoding this protein was not identified in A. fraterculus, suggesting diversification after the separation of the species [27]. The high divergence at the amino acid level and its identification in the Mexican fruit fly could indicate that this gene is present but not detected in A. fraterculus, probably because of low expression levels during the collection stages.
The CSP gene family is highly conserved in insects; comparative analyses of the genomes of different insect orders reported that, in contrast to OBPs, the repertoire of CSPs is markedly reduced from three members in the genome of D. melanogaster to 20 in Tribolium castaneum [55,61,62,63,64,65,66]. For Tephritidae, there are few published papers for these non-receptor families; in general, three to five members of CSPs have been reported in B. dorsalis and B. minax [29,30,31,67,68].
We obtained equivalent results in the transcriptome of A. ludens, with five members recorded for CSPs. Similar to Cheng et al. (2020) [31], we identified two genes that encode CSP3. According to the phylogenetic tree, AludCSP3-1 and AludCSP3-2 belong to two subgroups with a specific diversification, related to orthologs of different tephritid species, of which the CSPs of the database for B. dorsalis, B. tryoni, Z. tau, and Z. cucurbitae, also have two different genes encoding CSP3. However, only the cluster AludCSP3-1 is related to an orthologue with D. melanogaster A10/OSD.
Although this family is more conserved than OBPs and is more reduced in Diptera, knowledge about it is evolution and function in Tephritidae is limited to a few species [29,69,70], and most of the sequences for fruit flies are not reported and characterized. Even for D. melanogaster, three sequences homologous to the AludCSPs were not classified within the CSPs, even though they present the characteristic protein domains of this family. Therefore, it is necessary to carry out comparative studies of the information deposited in the world databases to unify the classification of CSP members in Tephritidae and define the phylogenetic relationships with the functional characterization.
On the other hand, we identified two genes encoding SNMPs in the head transcriptomes, both of which presented the characteristics of this family, including the two transmembrane regions and six conserved Cys residues [5,70]. Phylogenetic analysis corroborated the division of AludSNMPs and their orthologs into SNMP1 and SNMP2 subgroups with 100% BPP support [71,72,73]. SNMP1 splits into SNMP1a and SNMP1b subtypes, in four species of Bactrocera and two species of Zeugodacus; however, in the transcriptome of A. ludens, we only identified SNMP1b. In addition, we identified a gene encoding SNMP2, also described in antennae and non-olfactory tissues, with high expression in legs in B. minax [53,73].
Three gene families encode chemosensory receptors; insect ORs detect volatile chemical information and are very different from mammalian GPCR with an inverse heptahelical topology; the GRs that are part of the gustatory system perceive a broad spectrum of ligands, and those related to ionotropic glutamate receptors are called ionotropic receptors (IRs) [4]. In recent years, knowledge of chemosensory receptors in agriculturally important tephritid has increased; however, there are no reports or information in genetic databases (accessed on 1 September 2022) of these receptors in the genus Anastrepha [74,75,76]. This work reports for the first time the identification and analysis of eighty-two genes encoding the three-receptor families in the Mexican fruit fly.
The number of putative ORs we identified in A. ludens (42) is similar to that reported in B. dorsalis (43) and higher than other Bactrocera and Zeugodacus species that vary in the range of 39 to 41 genes. This result suggests that the ORs not only show conservation in gene number within the Dacinae group but is comparable to other tephritid subfamilies [30,77]. As in other insects, A. ludens has one ORCO gene that shares a high similarity (>90%) to co-receptors of the eight tephritid species analyzed. ORCO is coexpressed with a specific OR in nearly all olfactory neurons [11] and plays an essential role in response to common odors, including esters, aldehydes, ketones, aromatics, and terpenes; in the male oriental fruit fly B. dorsalis, Orco participates in the taxis to methyl eugenol [78].
A. ludens has expansions of AludOR7a/AludOR74a and specific duplications of AludOR71a, AludOR69a, and AludOR67d. Different studies report OR7a and OR74a as positive signaling receptors to aggregation compounds such as 9-tricosene and 1-nonanol/6-oxo-1-nonanol. Diversification of these genes in tephritid may enhance the perception of specific odorants or the combination of similar odorants, which may be related to species differences in the detection of host volatiles and pheromones [32,53,79,80].
Although the number of A. ludens ORs is identical to that of another Tephritidae, we found differences in the previously reported homologous. The putative receptors AludOR30a-like1, AludOR30a-like2, AludOR59a-like, and AludOR63a-like, only clustered with receptors from R. pomonella, R. zephyria derived from an annotated genomic sequence using the gene prediction method, whereas AludOR59b-like and AludOR5 clustered with two orthologues from B. minax and C. capitata. Although the six receptors presented ambiguous classification and did not cluster with receptors from other tephritid species or D. melanogaster, all, except AludOR59a-like and AludOR30a-like1, present complete ORFs and five to seven predicted transmembrane regions. Furthermore, the phylogenetic analysis clustered AludOR30a-like1 and AludOR30a-like2 in the same clade related to OR49b and present an identity of 30.77%.
On the other hand, GRs repertoires vary depending on the insects’ order, from 10 genes in the honeybee Apis melifera to 68 and 76 genes in D. melanogaster and the mosquito Anopheles gambiae [81,82,83,84]. In this work, we identified 24 GRs in the transcriptomes of A. ludens, which is higher than that reported in Tephritidae, from six GRs in Z. tau, to 15 in B. minax. Several studies have shown associations in the number of GRs and host range, where specialist insects conserve a lower number of GRs than generalists [74,85,86]. However, similar to the Mexican fly, the reported Tephritidae species are characterized by a broad host range [23]. Given that GRs present a high sequence divergence, making their identification difficult, it is necessary to perform comparative genome analyses to corroborate the GRs’ repertoire [32,74].
Phylogenetic analysis grouped the twenty-three AludGRs and their orthologues into three clades representing the different lineages conserved in insects with specific putative functions; AludGR63a, AludGR21a, and AludGR22 belong to CO2 receptors. The most conserved in insects. GR22, homologous to GR21a, was found expanded in Tephritidae, which could enhance their olfactory sensitivity to CO2 [53,87,88]. Seven AludGRs (AludGR64a-f and AludGR5a) were classified within the sugar receptors, and we identified a duplication of the AludGR64b gene related to the Dmel64b protein. The remaining 13 AludGRs group into the bitter/aversive receptor clade, representing a high proportion within this family with specific expansions. Some of them, such as AludGR98a, AludGR98b, AludGR2a, and AludGR22e, are only related to R. pomonella and D. melanogaster GRs [89,90,91].
We identified 17 putative IR/iGluR-encoding genes in A. ludens; this is lower than reported in six fruit fly species with an average of 23 IRs and D. melanogaster with 66 IRs [30,53,92]. Based on phylogenetic analysis and Croset et al., 2010 [93], the 17 AludIRs were classified into the subfamilies divergent IRs, antennal IRs, non-NMDA (N-Methyl-D-aspartic acid), iGluRs, and two putative co-receptors. Although the antennal subfamily represents only a fraction of the IR repertoire in most insects, in A. ludens transcriptomes, it was the most representative group with ten members. All AludIRs formed orthologous groups with D. melanogaster receptors, except for IR75a, which presents an expansion in tephritid and could enhance their ability to perceive fermented fruits [30,53,94].
Croset et al. (2010) [93] mentioned that dipteran insects have broad expansions of divergent IRs; however, in A. ludens, we only identified three members of this subfamily (AludIR7c, AludIR56c, and AludIR100a), which were closely related to three clusters of antennal IRs with 90% BPP. Several evolutionary and gene expression studies report the IR25a cluster as a putative co-receptor, analogous to the heteromeric assembly of iGluR subunits into functional complexes [95]. In addition, IR25 appears to be an atypical member that, together with its homolog IR8a, shows high conservation with the primary sequence to iGluRs. The above is consistent with the phylogenetic relationship of AludIR25a and AludIR8a, located within the iGluR family of non-NMDA receptors, suggesting that these IRs are among the ancestral groups [92,93,96].
In conclusion, this work represents the first effort to identify, classify and functionally characterize the six chemosensory gene families from the analysis of head transcriptomes of male and female adults of A. ludens in different physiological stages. The 154 described genes encoding putative proteins play an essential role in the chemosensory perception of the Mexican fruit fly. They are a significant contribution to the genetic databases of insects, particularly dipterans, which will provide a further basis for the study of molecular communication in the Tephritidae family, one of the most important for its economic and social impact worldwide. Moreover, the identification of olfactory genes provides information for the study and development of alternative strategies for the specific management of pests of agricultural importance.

4. Materials and Methods

4.1. Collection of Biological Material

In 2017, we collected sexually immature adult females and males between four and 11 days of the emergence of A. ludens from the laboratory-reared of the Colegio de Postgraduados, Campus Montecillo; heads were dissected and stored in tubes with 600 µL of Thermo Fisher RNAlater™ buffer (Thermo Fisher Scientific, Inc., Lithuania). On the other hand, in 2020, a direct collection of mangoes infested with wild A. ludens larvae was carried out in Tapachula, Chiapas. The fruits were placed in cages until the emergence of adults, which were kept in rearing with a 12 L (07:00–19:00 h):12 O photoperiod and fed on sweetened water: hydrolyzed protein mixture until sexual maturity. Then, sexually mature males of eighteen-day-old were collected at resting (10:00–12:00 a.m.) and calling (6:30–7:00 p.m.) physiology stages; the specie was taxonomically corroborated by the technical experts of the MOSCAMED program, and heads were dissected and preserved in 200 mL of RNAlater™.

4.2. Preparation and Sequencing of cDNA Libraries

For the breeding material collected in 2017, we performed four total RNA extractions with the SV Total RNA Isolation System kit (Promega, Madison, WI, USA), two replicates of twenty-five male and female heads (1:1) for each emergence time (4 and 11 days), and we evaluated the extraction quality with a Qubit™ fluorometer (Thermo Fisher Scientific, Inc., Waltham, MA, USA). Four paired libraries were constructed with the TruSeq RNA Sample Preparation Kit V3 (Illumina Inc., San Diego, CA, USA) and sequenced with the MiSeq platform at the Centro Nacional de Referencia Fitosanitaria, SENASICA, Tecámac, State of Mexico. On the other hand, the heads of wild males collected in mango were sent to BGI Hong Kong, China. We performed RNA extraction in three replicates of twenty-five male heads for each physiological stage; the quality and concentration were verified with an Agilent 2100 Bioanalyzer (Agilent Technologies, Inc., Santa Clara, CA, USA) and constructed six paired cDNA libraries with the MGIEasy-RNA Library-Prep kit (MGI Tech Co., Ltd., Hong Kong, China). Sequencing was performed with the BGISEQ-500 platform (BGI Group, Hong Kong, China), with a paired read length of 100 bp and an output of 4 GB of clean data.

4.3. Cleaning, Assembly, and Annotation of the De Novo Transcriptome

Ten paired libraries were sequenced with Illumina and BGI platforms and evaluated the quality with FastQC v0.10.1 (Babraham Institute, Cambridge) [97]. The cleanup of adapters and low-quality reads was performed based on the characteristics of each platform with FastP (https://github.com/OpenGene/fastp, accessed on 1 September 2021) [98]. The libraries were assembled de novo based on the sequencing platform with Trinity v2.06 (https://github.com/trinityrnaseq/trinityrnaseq, accessed on 1 September 2021) and a default value of kmer = 25. The quality of the assemblies was evaluated with QUAST v.0.3.0. (https://github.com/ablab/quast, accessed on 1 September 2021) [99]. We recorded the total number of reads generated, length N50, mean length, and the GC% and filtered out redundancies to obtain the final number of unigenes for each assembly. We performed the annotation of the unigenes of the two assemblies by homology analysis with BLASTx [100] against the Insecta database of UniprotKB (3,914,796 sequences) with an E-value of 1 × 10−7. We added a conserved protein domain search with InteproScan workflow [87] and obtained the Gene Ontology (GO) terms with HMMER2GO (https://github.com/sestaton/HMMER2GO, accessed on 1 October 2021).

4.4. Mining and Functional Analysis of Chemosensory Genes

We filtered the sequences of A. ludens chemosensory receptor and non-receptor protein families from the results obtained in BLASTx homology analysis against the UniprotKB base and conserved InterProScan domains. We verified the open reading frames of the sequences with ORFinder (https://www.ncbi.nlm.nih.gov/orffinder/, accessed on 1 November 2021) and performed Blastp analysis against the UniprotKB database (accessed on 1 February 2022), and for the receptors, with the NCBI non-redundant protein base (accessed on 1 February 2022). To predict the protein domains of the receptor proteins (including SNMPs), we used TMHMM 2.0 (DTU Health Tech, Lyngby, Denmark) [101] and SignalP-5.0. (DTU Health Tech, Lyngby, Denmark) [102] to predict the presence of signal peptides in OBPs and CSPs. We evaluated the sequences by multialignment. For OBPs, we generated a classification based on cysteine profiles and obtained the average amino acid identities with the MAFFT percent identity matrix. The final sequences of the six chemosensory protein families were named Alud, followed by the corresponding protein name based on the best match in the Blastp analysis.

4.5. Sequence Comparison and Phylogenetic Analysis

We generated six different datasets for each candidate chemosensory families of A. ludens and homologous protein sequences from diverse tephritid species (B. dorsalis, B. latifrons, B. correcta, B. minax, B. tryoni, B. oleae, Z. cucurbitae, Z. tau, C. capitata, R. pomonella, R. zephyria, P. utilis, A. obliqua, and A. fraterculus) and the model species D. melanogaster. In addition, we included proteins from the mosquito Anopheles gambiae for SNMPs analysis (Supplements Tables S9 and S10). For multiple alignments of OBPs and CSPs, Clustal Omega (UCD Dublin, Dublin, Ireland) [103] was used with default parameters, while for SNMPs and receptors; we made alignments with the iterative method with MAFFT refinement [104]. We visualized the alignments in AliView (Uppsala University, Uppsala, Sweden) [105]) and eliminated the extremes in short or highly divergent length sequences. For each dataset, we obtained the best protein evolutionary model in ModelTest-NG (https://github.com/ddarriba/modeltest, accessed on 1 April 2022) [106] using the AIC and BIC criteria and constructed six phylogenetic trees by Bayesian inference with the Markov chain Monte Carlo method in BEAST v1.10.4 (© BEAST Developers) [107] with a minimum of 15,000,000 generations (Supplements Table S11). We evaluated the MCM output in Tracer v1.7.1 (©BEAST Developers) and incremented the chains for those analyses that required it. We discarded 30% of the initial trees for tree annotation and calculated posterior probability with the remaining trees. Finally, we visualized and edited the consensus trees in iTOL v6.5.8 (European Molecular Biology Laboratory, © EMBL) [108].

4.6. Sequence Information

A. ludens chemosensory gene sequences were submitted to the NCBI to obtain accession numbers for AludOBPs (ON419948-ON419978), AludCSPs (ON419979-ON419981), AludSNMPs (ON419982 and ON419983), AludORs (ON419984-ON420025), AludGRs (ON420026-ON420055), and AludIRs (ON420056-ON420076).

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms231810531/s1.

Author Contributions

Conceptualization, supervision, writing—review and editing, funding acquisition, O.L.S.-L.; investigation, formal analysis, data curation and writing—original draft preparation, B.T.-H.; methodology, A.R.E.-P.; writing—review and editing, J.C.-T.; methodology, F.d.l.C.H.-H.; methodology, J.L.C.-J.; methodology, resources, J.S.M.-H.; methodology, F.S.-G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by CONACyT and an agreement between Colegio de Postgraduados Montecillo and Servicio Nacional de Sanidad, Inocuidad y Calidad Agroalimentaria (SENASICA) (CP-SENASICA-2019).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Colegio de Postgraduados, Campus Montecillo. Consejo Nacional de Ciencia y Tecnología (CONACyT). Francisco Javier Trujillo Arriaga and José Abel López Buenfil, for the support provided to have access to the Centro Nacional de Referencia Fitosanitaria SENASICA, for the preparation of RNA-Seq libraries.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hansson, B.S.; Stensmyr, M.C. Evolution of insect olfaction. Neuron 2011, 72, 698–711. [Google Scholar] [CrossRef] [PubMed]
  2. Sachse, S.; Krieger, J. Olfaction in insects—The primary processes of odor recognition and coding. E-Neuroforum 2011, 2, 49–60. [Google Scholar] [CrossRef]
  3. Tunstall, N.E.; Warr, C.G. Chemical communication in insects: The peripheral odour coding system of Drosophila melanogaster. Adv. Exp. Med. Biol. 2012, 739, 59–77. [Google Scholar]
  4. Leal, W.S. Odorant Reception in Insects: Roles of Receptors, Binding Proteins, and Degrading Enzymes. Annu. Rev. Entomol. 2013, 58, 373–391. [Google Scholar] [CrossRef]
  5. Benton, R.; Vannice, K.S.; Vosshall, L.B. An essential role for a CD36-related receptor in pheromone detection in Drosophila. Nature 2007, 450, 289–293. [Google Scholar] [CrossRef] [PubMed]
  6. Gomez-Diaz, C.; Bargeton, B.; Abuin, L.; Bukar, N.; Reina, J.H.; Bartoi, T.; Graf, M.; Ong, H.; Ulbrich, M.H.; Masson, J.-F.; et al. A CD36 ectodomain mediates insect pheromone detection via a putative tunnelling mechanism. Nat. Commun. 2016, 7, 11866. [Google Scholar] [CrossRef] [PubMed]
  7. Pelosi, P.; Iovinella, I.; Felicioli, A.; Dani, F.R. Soluble proteins of chemical communication: An overview across arthropods. Front. Physiol. 2014, 5, 320. [Google Scholar] [CrossRef] [PubMed]
  8. Isman, M.B. Challenges of Pest Management in the Twenty First Century: New Tools and Strategies to Combat Old and New Foes Alike. Front. Agron. 2019, 1, 2. [Google Scholar] [CrossRef]
  9. Bosa, C.F.; Cruz-López, L.; Zepeda-Cisneros, C.S.; Valle-Mora, J.; Guillén-Navarro, K.; Liedo, P. Sexual behavior and male volatile compounds in wild and mass-reared strains of the Mexican fruit fly Anastrepha ludens (Diptera: Tephritidae) held under different colony management regimes. Insect Sci. 2016, 23, 105–116. [Google Scholar] [CrossRef]
  10. Pelosi, P.; Iovinella, I.; Zhu, J.; Wang, G.; Dani, F.R. Beyond chemoreception: Diverse tasks of soluble olfactory proteins in insects. Biol. Rev. 2017, 93, 184–200. [Google Scholar] [CrossRef]
  11. Dimitratos, S.D.; Hommel, A.S.; Konrad, K.D.; Simpson, L.M.; Wu-Woods, J.J.; Woods, D.F. Biosensors to monitor water quality utilizing insect odorant-binding proteins as detector elements. Biosensors 2019, 9, 62. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Ingham, V.A.; Anthousi, A.; Douris, V.; Harding, N.J.; Lycett, G.; Morris, M.; Vontas, J.; Ranson, H. A sensory appendage protein protects malaria vectors from pyrethroids. Nature 2019, 577, 376–380. [Google Scholar] [CrossRef] [PubMed]
  13. Portilla-Pulido, J.S.; Castillo-Morales, R.M.; Barón-Rodríguez, M.A.; Duque, J.E.; Mendez-Sanchez, S.C. Design of a repellent against Aedes aegypti (Diptera: Culicidae) using in silico simulations with AaegOBP1 protein. J. Med. Entomol. 2019, 57, 463–476. [Google Scholar] [CrossRef] [PubMed]
  14. Cheema, J.A.; Carraher, C.; Plank, N.O.V.; Travas-Sejdic, J.; Kralicek, A. Insect odorant receptor-based biosensors: Current status and prospects. Biotechnol. Adv. 2021, 53, 107840. [Google Scholar] [CrossRef]
  15. Chen, X.; Lei, Y.; Li, H.; Xu, L.; Yang, H.; Wang, J.; Jiang, H. CRISPR/Cas9 mutagenesis abolishes odorant-binding protein BdorOBP56f-2 and impairs the perception of methyl eugenol in Bactrocera dorsalis (Hendel). Insect Biochem. Mol. Biol. 2021, 139, 103656. [Google Scholar] [CrossRef]
  16. Koutroumpa, F.; Monsempès, C.; Anton, S.; François, M.-C.; Montagné, N.; Jacquin-Joly, E. Pheromone Receptor Knock-Out Affects Pheromone Detection and Brain Structure in a Moth. Biomolecules 2022, 12, 341. [Google Scholar] [CrossRef]
  17. Martinazzo, J.; Ballen, C.S.; Steffens, J.; Steffens, C. Sensing of pheromones from Euschistus heros (F.) stink bugs by nanosensors. Sens. Actuators Rep. 2022, 4, 100071. [Google Scholar] [CrossRef]
  18. Trombik, J.; Ward, S.F.; Norrbom, A.L.; Liebhold, A.M. Global drivers of historical true fruit fly (Diptera: Tephritidae) invasions. J. Pest. Sci. 2022, 1–13. [Google Scholar] [CrossRef]
  19. Malacrida, A.R.; Gomulski, L.M.; Bonizzoni, M.; Bertin, S.; Gasperi, G.; Guglielmino, C.R. Globalization and fruitfly invasion and expansion: The medfly paradigm. Genetica 2007, 131, 1–9. [Google Scholar] [CrossRef]
  20. Vera, M.T.; Cáceres, C.; Wornoayporn, V.; Islam, A.; Robinson, A.S.; De La Vega, M.H.; Hendrichs, J.; Cayol, J.-P. Mating Incompatibility Among Populations of the South American Fruit Fly Anastrepha fraterculus (Diptera: Tephritidae). Ann. Entomol. Soc. Am. 2006, 99, 387–397. [Google Scholar] [CrossRef]
  21. Caceres, C.; Segura, D.F.; Vera, M.T.; Wornoayporn, V.; Cladera, J.L.; Teal, P.; Sapountzis, P.; Bourtzis, K.; Zacharopoulou, A.; Robinson, S.A. Incipient speciation revealed in Anastrepha fraterculus (Diptera: Tephritidae) by studies on mating compatibility, sex pheromones, hybridization and cytology. Biol. J. Linn. Soc. 2009, 97, 152–165. [Google Scholar] [CrossRef]
  22. Hernandez-Ortíz, V. Distribución y biogeografía del género Anastrepha en México. In Moscas de la Fruta en Latinoamérica, Diversidad, Biología y Manejo; Hernández-Ortiz, V., Ed.; S y G Editores: Mexico City, Mexico, 2007; pp. 53–76. ISBN 9687623284. [Google Scholar]
  23. Aluja, M.; Mangan, R.L. Fruit Fly (Diptera: Tephritidae) Host Status Determination: Critical Conceptual, Methodological, and Regulatory Considerations. Annu. Rev. Entomol. 2008, 53, 473–502. [Google Scholar] [CrossRef] [PubMed]
  24. Servicio Nacional de Sanidad, Inocuidad y Calidad Agroalimentaria (SENASICA). Available online: https://www.gob.mx/senasica/documentos/ficha-tecnica (accessed on 14 June 2022).
  25. Gazit, Y.; Rossler, Y.; Wang, S.; Tang, J.; Lurie, S. Thermal Death Kinetics of Egg and Third Instar Mediterranean Fruit Fly (Diptera: Tephritidae). J. Econ. Entomol. 2004, 97, 1540–1546. [Google Scholar] [CrossRef] [PubMed]
  26. Siciliano, P.; He, X.L.; Woodcock, C.; Pickett, J.A.; Field, L.M.; Birkett, M.A.; Kalinova, B.; Gomulski, L.M.; Scolari, F.; Gasperi, G.; et al. Identification of pheromone components and their binding affinity to the odorant binding protein CcapOBP83a-2 of the Mediterranean fruit fly, Ceratitis capitata. Insect Biochem. Mol. Biol. 2014, 48, 51–62. [Google Scholar] [CrossRef] [PubMed]
  27. Campanini, E.B.; de Brito, R.A. Molecular evolution of Odorant-binding proteins gene family in two closely related Anastrepha fruit flies. BMC Evol. Biol. 2016, 16, 198. [Google Scholar] [CrossRef]
  28. Elfekih, S.; Chen, C.-Y.; Hsu, J.-C.; Belcaid, M.; Haymer, D. Identification and preliminary characterization of chemosensory perception-associated proteins in the melon fly Bactrocera cucurbitae using RNA-seq. Sci. Rep. 2016, 6, 19112. [Google Scholar] [CrossRef] [PubMed]
  29. Wang, Z.; Qi, Z.; Chen, J.; Wang, F.; Gui, L.; Zhang, G. Molecular characterization of chemosensory protein genes in Bactrocera minax (Diptera: Tephritidae). Entomol. Res. 2021, 51, 349–356. [Google Scholar] [CrossRef]
  30. Wu, Z.; Kang, C.; Qu, M.; Chen, J.; Chen, M.; Bin, S.; Lin, J. Candidates for chemosensory genes identified in the Chinese citrus fly, Bactrocera minax, through a transcriptomic analysis. BMC Genom. 2019, 20. [Google Scholar] [CrossRef]
  31. Cheng, J.-F.; Yu, T.; Chen, Z.-J.; Chen, S.; Chen, Y.-P.; Gao, L.; Zhang, W.-H.; Jiang, B.; Bai, X.; Walker, E.D.; et al. Comparative genomic and transcriptomic analyses of chemosensory genes in the citrus fruit fly Bactrocera (Tetradacus) minax. Sci. Rep. 2020, 10, 18068. [Google Scholar] [CrossRef]
  32. Ono, H.; Hee, A.K.W.; Jiang, H. Recent advancements in studies on chemosensory mechanisms underlying detection of semiochemicals in Dacini fruit flies of economic importance (Diptera: Tephritidae). Insects 2021, 12, 106. [Google Scholar] [CrossRef]
  33. Wachi, N.; Matsubayashi, K.W.; Maeto, K. Application of next-generation sequencing to the study of non-model insects. Entomol. Sci. 2017, 21, 3–11. [Google Scholar] [CrossRef] [Green Version]
  34. Ruiz-May, E.; Altúzar-Molina, A.; Elizalde-Contreras, J.M.; Arellano-de los Santos, J.; Monribot-Villanueva, J.; Guillén, L.; Vázquez-Rosas-Landa, M.; Ibarra-Laclette, E.; Ramírez-Vázquez, M.; Ortega, R.; et al. A first glimpse of the mexican fruit fly Anastrepha ludens (Diptera: Tephritidae) antenna morphology and proteome in response to a proteinaceous Attractant. Int. J. Mol. Sci. 2020, 21, 8086. [Google Scholar] [CrossRef] [PubMed]
  35. Perre, P.; Jorge, L.R.; Lewinsohn, T.M.; Zucchi, R.A. Morphometric differentiation of fruit fly pest species of the Anastrepha fraterculus group (Diptera: Tephritidae). Ann. Entomol. Soc. Am. 2014, 107, 490–495. [Google Scholar] [CrossRef]
  36. Li, X.-M.; Zhu, X.-Y.; He, P.; Xu, L.; Sun, L.; Chen, L.; Wang, Z.-Q.; Deng, D.-G.; Zhang, Y.-N. Molecular characterization and sex distribution of chemosensory receptor gene family based on transcriptome analysis of Scaeva pyrastri. PLoS ONE 2016, 11, e0155323. [Google Scholar] [CrossRef]
  37. Jia, X.; Zhang, X.; Liu, H.; Wang, R.; Zhang, T. Identification of chemosensory genes from the antennal transcriptome of Indian meal moth Plodia interpunctella. PLoS ONE 2018, 13, e0189889. [Google Scholar] [CrossRef] [PubMed]
  38. Gu, S.-H.; Yang, R.-N.; Guo, M.-B.; Wang, G.-R.; Wu, K.-M.; Guo, Y.-Y.; Zhou, J.-J.; Zhang, Y.-J. Molecular identification and differential expression of sensory neuron membrane proteins in the antennae of the black cutworm moth Agrotis ipsilon. J. Insect Physiol. 2013, 59, 430–443. [Google Scholar] [CrossRef]
  39. Wang, B.; Liu, Y.; Wang, G.-R. Chemosensory genes in the antennal transcriptome of two syrphid species, Episyrphus balteatus and Eupeodes corollae (Diptera: Syrphidae). BMC Genom. 2017, 18, 586. [Google Scholar] [CrossRef]
  40. Torres-Huerta, B.; Segura-León, O.L.; Aragón-Magadan, M.A.; González-Hernández, H. Identification and motif analyses of candidate nonreceptor olfactory genes of Dendroctonus adjunctus Blandford (Coleoptera: Curculionidae) from the head transcriptome. Sci. Rep. 2020, 10, 20695. [Google Scholar] [CrossRef]
  41. Archunan, G. Odorant Binding Proteins: A key player in the sense of smell. Bioinformation 2018, 14, 36–37. [Google Scholar] [CrossRef]
  42. Chen, X.-F.; Xu, L.; Zhang, Y.-X.; Wei, D.; Wang, J.-J.; Jiang, H.-B. Genome-wide identification and expression profiling of odorant-binding proteins in the oriental fruit fly, Bactrocera dorsalis. Comp. Biochem. Physiol. Part D Genom. Proteom. 2019, 31, 100605. [Google Scholar] [CrossRef]
  43. Liu, Z.; Liang, X.-F.; Xu, L.; Keesey, I.W.; Lei, Z.-R.; Smagghe, G.; Wang, J.-J. An antennae-specific odorant-binding protein is involved in Bactrocera dorsalis olfaction. Front. Ecol. Evol. 2020, 8, 63. [Google Scholar] [CrossRef] [Green Version]
  44. Liu, Y.; Cui, Z.; Si, P.; Liu, Y.; Zhou, Q.; Wang, G. Characterization of a specific odorant receptor for linalool in the Chinese citrus fly Bactrocera minax (Diptera: Tephritidae). Insect Biochem. Mol. Biol. 2020, 122, 103389. [Google Scholar] [CrossRef] [PubMed]
  45. Garczynski, S.F.; Martin, J.A.; Griste, M.; Willett, L.S.; Cooper, W.R.; Swisher, K.D.; Unruh, T.R. CRISPR/Cas9 Editing of the Codling Moth (Lepidoptera: Tortricidae) CpomOR1 Gene Affects Egg Production and Viability. J. Econ. Entomol. 2017, 110, 1847–1855. [Google Scholar] [CrossRef] [PubMed]
  46. Chen, J.; Yang, L.; Tian, X.-L.; Gui, L.-Y.; Wang, F.-L.; Zhang, G.-H. Functional characterization of two antenna-enriched odorant-binding proteins from Bactrocera minax (Diptera: Tephritidae). J. Econ. Entomol. 2021, 144, 2361–2369. [Google Scholar] [CrossRef]
  47. Yao, R.; Zhao, M.; Zhong, L.; Li, Y.; Li, D.; Deng, Z.; Ma, X. Characterization of the binding ability of the odorant binding protein BminOBP9 of Bactrocera minax to citrus volatiles. Pest Manag. Sci. 2021, 77, 1214–1225. [Google Scholar] [CrossRef]
  48. Hekmat-Scafe, D.S.; Scafe, C.R.; McKinney, A.J.; Tanouye, M.A. Genome-wide analysis of the odorant-binding protein gene family in Drosophila melanogaster. Genome Res. 2002, 12, 1357–1369. [Google Scholar] [CrossRef]
  49. Vieira, F.G.; Rozas, J. Comparative genomics of the odorant-binding and chemosensory protein gene families across the Arthropoda: Origin and evolutionary history of the chemosensory system. Genome Biol. Evol. 2011, 3, 476–490. [Google Scholar] [CrossRef]
  50. Nakamura, A.M.; Chahad-Ehlers, S.; Lima, A.L.A.; Taniguti, C.H.; Sobrinho, I.; Torres, F.R.; de Brito, R.A. Reference genes for accessing differential expression among developmental stages and analysis of differential expression of OBP genes in Anastrepha obliqua. Sci. Rep. 2016, 6, 17480. [Google Scholar] [CrossRef]
  51. Campanini, E.B.; Congrains, C.; Torres, F.R.; de Brito, R.A. Odorant-binding proteins expression patterns in recently diverged species of Anastrepha fruit flies. Sci. Rep. 2017, 7, 2194. [Google Scholar] [CrossRef]
  52. Cheng, J.-F.; Yu, T.; Chen, Z.-J.; Chen, Y.-P.; Gao, L.; Zhang, W.-H.; Jiang, B.; Bai, X.; Liu, J.; Lu, Y. Identification and expression analysis of chemosensory genes in the citrus fruit fly Bactrocera (Tetradacus) minax. PeerJ Prepr. 2018, 6, e27297v1. [Google Scholar] [CrossRef]
  53. Wu, Z.; Cui, Y.; Ma, J.; Qu, M.; Lin, J. Analyses of chemosensory genes provide insight into the evolution of behavioral differences to phytochemicals in Bactrocera species. Mol. Phylogenet. Evol. 2020, 151, 106858. [Google Scholar] [CrossRef] [PubMed]
  54. Pelosi, P.; Zhou, J.-J.; Ban, L.P.; Calvello, M. Soluble proteins in insect chemical communication. Cell Mol. Life Sci. 2006, 63, 1658–1676. [Google Scholar] [CrossRef] [PubMed]
  55. Vieira, F.G.; Sánchez-Gracia, A.; Rozas, J. Comparative genomic analysis of the odorant-binding protein family in 12 Drosophila genomes: Purifying selection and birth-and-death evolution. Genome Biol. 2007, 8, R235. [Google Scholar] [CrossRef]
  56. Rihani, K.; Ferveur, J.-F.; Briand, L. The 40-Year Mystery of Insect Odorant-Binding Proteins. Biomolecules 2021, 11, 509. [Google Scholar] [CrossRef]
  57. Missbach, C.; Vogel, H.; Hansson, B.S.; Groβe-Wilde, E. Identification of odorant binding proteins and chemosensory proteins in antennal transcriptomes of the jumping bristletail Lepismachilis y-signata and the firebrat Thermobia domestica: Evidence for an independent OBP–OR origin. Chem. Senses 2015, 40, 615–626. [Google Scholar] [CrossRef] [PubMed]
  58. He, X.; He, Z.-B.; Zhang, Y.-J.; Zhou, Y.; Xian, P.-J.; Qiao, L.; Chen, B. Genome-wide identification and characterization of odorant-binding protein (OBP) genes in the malaria vector Anopheles sinensis (Diptera: Culicidae). Insect Sci. 2016, 23, 366–376. [Google Scholar] [CrossRef]
  59. Sun, J.S.; Xiao, S.; Carlson, J.R. The diverse small proteins called odorant-binding proteins. Open Biol. 2018, 8, 180208. [Google Scholar] [CrossRef]
  60. Simpson, M.G. Phylogenetic Systematics. In Plant Systematics; Elsevier-Academic Press: San Diego, CA, USA, 2010; pp. 17–52. [Google Scholar] [CrossRef]
  61. Pelosi, P.; Calvello, M.; Ban, L. Diversity of odorant-binding proteins and chemosensory proteins in insects. Chem. Senses 2005, 30, i291–i292. [Google Scholar] [CrossRef]
  62. Forêt, S.; Maleszka, R. Function and evolution of a gene family encoding odorant binding-like proteins in a social insect, the honey bee (Apis mellifera). Genome Res. 2006, 16, 1404–1413. [Google Scholar] [CrossRef]
  63. Forêt, S.; Wanner, K.W.; Maleszka, R. Chemosensory proteins in the honey bee: Insights from the annotated genome, comparative analyses and expressional profiling. Insect Biochem. Mol. Biol. 2007, 37, 19–28. [Google Scholar] [CrossRef]
  64. Gong, D.-P.; Zhang, H.-J.; Zhao, P.; Lin, Y.; Xia, Q.-Y.; Xiang, Z.-H. Identification and expression pattern of the chemosensory protein gene family in the silkworm, Bombyx mori. Insect Biochem. Mol. Biol. 2007, 37, 266–277. [Google Scholar] [CrossRef] [PubMed]
  65. Gong, D.-P.; Zhang, H.-J.; Zhao, P.; Xia, Q.-Y.; Xiang, Z.-H. The odorant binding protein gene family from the genome of silkworm, Bombyx mori. BMC Genom. 2009, 10, 332. [Google Scholar] [CrossRef] [PubMed]
  66. Dippel, S.; Oberhofer, G.; Kahnt, J.; Gerischer, L.; Opitz, L.; Schachtner, J.; Stanke, M.; Schütz, S.; Wimmer, E.A.; Angeli, S. Tissue-specific transcriptomics, chromosomal localization, and phylogeny of chemosensory and odorant binding proteins from the red flour beetle Tribolium castaneum reveal subgroup specificities for olfaction or more general functions. BMC Genom. 2014, 15, 1141. [Google Scholar] [CrossRef]
  67. Wu, Z.; Zhang, H.; Wang, Z.; Bin, S.; He, H.; Lin, J. Discovery of chemosensory genes in the oriental fruit fly, Bactrocera dorsalis. PLoS ONE 2015, 10, e0129794. [Google Scholar] [CrossRef] [PubMed]
  68. Liu, Z.; Smagghe, G.; Lei, Z.; Wang, J.-J. Identification of male- and female-specific olfaction genes in antennae of the oriental fruit fly (Bactrocera dorsalis). PLoS ONE 2016, 11, e0147783. [Google Scholar] [CrossRef]
  69. Xu, P.; Wang, Y.; Akami, M.; Niu, C.-Y. Identification of olfactory genes and functional analysis of BminCSP and BminOBP21 in Bactrocera minax. PLoS ONE 2019, 14, e0222193. [Google Scholar] [CrossRef]
  70. Yi, X.; Zhao, H.; Dong, X.; Wang, P.; Hu, M.; Zhong, G. BdorCSP2 Is Important for Antifeed and Oviposition-Deterring Activities Induced by Rhodojaponin-III against Bactrocera dorsalis. PLoS ONE 2013, 8, e77295. [Google Scholar] [CrossRef]
  71. Vogt, R.G. Biochemical diversity of odor detection: OBPs, ODEs and SNMPs. In Insect Pheromone Biochemistry and Molecular Biology: The Biosynthesis and Detection of Pheromones and Plant Volatiles; Blomquist, G., Vogt, R.G., Eds.; Elsevier Academic Press: San Diego, CA, USA, 2003; pp. 391–445. ISBN 9780121071516. [Google Scholar]
  72. Nichols, Z.; Vogt, R.G. The SNMP/CD36 gene family in Diptera, Hymenoptera and Coleoptera: Drosophila melanogaster, D. pseudoobscura, Anopheles gambiae, Aedes aegypti, Apis mellifera, and Tribolium castaneum. Insect Biochem. Mol. Biol. 2008, 38, 398–415. [Google Scholar] [CrossRef]
  73. Liu, C.; Zhang, J.; Liu, Y.; Wang, G.; Dong, S. Expression of SNMP1 and SNMP2 genes in antennal sensilla of Spodoptera exigua (Hübner). Arch. Insect Biochem. Physiol. 2014, 85, 114–126. [Google Scholar] [CrossRef]
  74. Isono, K.; Morita, H. Molecular and cellular designs of insect taste receptor system. Front. Cell Neurosci. 2010, 4, 20. [Google Scholar] [CrossRef] [Green Version]
  75. Rytz, R.; Croset, V.; Benton, R. Ionotropic receptors (IRs): Chemosensory ionotropic glutamate receptors in Drosophila and beyond. Insect Biochem. Mol. Biol. 2013, 43, 888–897. [Google Scholar] [CrossRef] [PubMed]
  76. Wicher, D. Tuning insect odorant receptors. Front. Cell Neurosci. 2018, 12, 94. [Google Scholar] [CrossRef] [PubMed]
  77. Jin, S.; Zhou, X.; Gu, F.; Zhong, G.; Yi, X. Olfactory plasticity: Variation in the expression of chemosensory receptors in Bactrocera dorsalis in different physiological states. Front. Physiol. 2017, 8, 672. [Google Scholar] [CrossRef] [PubMed]
  78. Zheng, W.; Zhu, C.; Peng, T.; Zhang, H. Odorant receptor co-receptor Orco is upregulated by methyl eugenol in male Bactrocera dorsalis (Diptera: Tephritidae). J. Insect Physiol. 2012, 58, 1122–1127. [Google Scholar] [CrossRef]
  79. Lin, C.-C.; Prokop-Prigge, K.A.; Preti, G.; Potter, C.J. Food odors trigger Drosophila males to deposit a pheromone that guides aggregation and female oviposition decisions. Neuroscience 2015, 4, e08688. [Google Scholar] [CrossRef]
  80. Miyazaki, H.; Otake, J.; Mitsuno, H.; Ozaki, K.; Kanzaki, R.; Chieng, A.C.-T.; Hee, A.K.-W.; Nishida, R.; Ono, H. Functional characterization of olfactory receptors in the oriental fruit fly Bactrocera dorsalis that respond to plant volatiles. Insect Biochem. Mol. Biol. 2018, 101, 32–46. [Google Scholar] [CrossRef]
  81. Hill, C.A.; Fox, A.N.; Pitts, R.J.; Kent, L.B.; Tan, P.L.; Chrystal, M.A.; Cravchik, A.; Collins, F.H.; Robertson, H.M.; Zwiebel, L.J. G protein-coupled receptors in Anopheles gambiae. Science 2002, 298, 176–178. [Google Scholar] [CrossRef]
  82. Robertson, H.M.; Warr, C.G.; Carlson, J.R. Molecular evolution of the insect chemoreceptor gene superfamily in Drosophila melanogaster. Proc. Natl. Acad. Sci. USA 2003, 100, 14537–14542. [Google Scholar] [CrossRef]
  83. Robertson, H.M.; Wanner, K.W. The chemoreceptor superfamily in the honey bee, Apis mellifera: Expansion of the odorant, but not gustatory, receptor family. Genome Res. 2006, 16, 1395–1403. [Google Scholar] [CrossRef]
  84. Sánchez-Gracia, A.; Vieira, F.G.; Rozas, J. Molecular evolution of the major chemosensory gene families in insects. Heredity 2009, 103, 208–216. [Google Scholar] [CrossRef]
  85. Gouin, A.; Bretaudeau, A.; Nam, K.; Gimenez, S.; Aury, J.-M.; Duvic, B.; Hilliou, F.; Durand, N.; Montagné, N.; Darboux, I.; et al. Two genomes of highly polyphagous lepidopteran pests (Spodoptera frugiperda, Noctuidae) with different host-plant ranges. Sci. Rep. 2017, 7, 11816. [Google Scholar] [CrossRef] [PubMed]
  86. Pearce, S.L.; Clarke, D.F.; East, P.D.; Elfekih, S.; Gordon, K.H.J.; Jermiin, L.S.; McGaughran, A.; Oakeshott, J.G.; Papanikolaou, A.; Perera, O.P.; et al. Genomic innovations, transcriptional plasticity and gene loss underlying the evolution and divergence of two highly polyphagous and invasive Helicoverpa pest species. BMC Biol. 2017, 15, 63. [Google Scholar] [CrossRef] [PubMed]
  87. Jones, P.; Binns, D.; Chang, H.-Y.; Fraser, M.; Li, W.; McAnulla, C.; McWilliam, H.; Maslen, J.; Mitchell, A.; Nuka, G.; et al. InterProScan 5: Genome-scale protein function classification. Bioinformatics 2014, 30, 1236–1240. [Google Scholar] [CrossRef] [PubMed]
  88. Kwon, J.Y.; Dahanukar, A.; Weiss, L.A.; Carlson, J.R. The molecular basis of CO2 reception in Drosophila. Proc. Natl. Acad. Sci. USA 2007, 104, 3574–3578. [Google Scholar] [CrossRef] [PubMed]
  89. Shim, J.; Lee, Y.; Jeong, Y.T.; Kim, Y.; Lee, M.G.; Montell, C.; Moon, S.J. The full repertoire of Drosophila gustatory receptors for detecting an aversive compound. Nat. Commun. 2015, 6, 8867. [Google Scholar] [CrossRef] [PubMed]
  90. Sung, H.Y.; Jeong, Y.T.; Lim, J.Y.; Kim, H.; Oh, S.M.; Hwang, S.W.; Oh, S.M.; Hwang, S.W.; Kwon, J.Y.; Moon, S.J. Heterogeneity in the Drosophila gustatory receptor complexes that detect aversive compounds. Nat. Commun. 2017, 8, 1484. [Google Scholar] [CrossRef]
  91. Aryal, B.; Lee, Y. Histamine gustatory aversion in Drosophila melanogaster. Insect Biochem. Mol. Biol. 2021, 134, 103586. [Google Scholar] [CrossRef]
  92. Benton, R.; Vannice, K.S.; Gomez-Diaz, C.; Vosshall, L.B. Variant ionotropic glutamate receptors as chemosensory receptors in Drosophila. Cell 2009, 136, 149–162. [Google Scholar] [CrossRef]
  93. Croset, V.; Rytz, R.; Cummins, S.F.; Budd, A.; Brawand, D.; Kaessmann, H.; Gibson, T.J.; Benton, R. Ancient protostome origin of chemosensory ionotropic glutamate receptors and the evolution of insect taste and olfaction. PLoS Genet. 2010, 6, e1001064. [Google Scholar] [CrossRef]
  94. Ai, M.; Min, S.; Grosjean, Y.; Leblanc, C.; Bell, R.; Benton, R.; Suh, G.S.B. Acid sensing by the Drosophila olfactory system. Nature 2010, 468, 691–695. [Google Scholar] [CrossRef] [Green Version]
  95. Gereau, R.W.; Swanson, G.T. The Glutamate Receptors; Humana Press: Totawa, NJ, USA, 2008; ISBN 9781588297921. [Google Scholar]
  96. Abuin, L.; Bargeton, B.; Ulbrich, M.H.; Isacoff, E.Y.; Kellenberger, S.; Benton, R. Functional architecture of olfactory ionotropic glutamate receptors. Neuron 2011, 69, 44–60. [Google Scholar] [CrossRef] [PubMed]
  97. Andrews, S. FastQC: A Quality Control Tool for High Troughput Sequence Data. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 14 June 2022).
  98. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. Fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef] [PubMed]
  99. Bushmanova, E.; Antipov, D.; Lapidus, A.; Suvorov, V.; Prjibelski, A.D. rnaQUAST: A quality assessment tool forde novotranscriptome assemblies. Bioinformatics 2016, 32, 2210–2212. [Google Scholar] [CrossRef] [PubMed]
  100. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST+: Architecture and applications. BMC Bioinform. 2009, 10, 421. [Google Scholar] [CrossRef]
  101. Krogh, A.; Larsson, B.; von Heijne, G.; Sonnhammer, E.L. Predicting transmembrane protein topology with a hidden markov model: Application to complete genomes. J. Mol. Biol. 2001, 305, 567–580. [Google Scholar] [CrossRef]
  102. Almagro, A.J.J.; Tsirigos, K.D.; Sønderby, C.K.; Petersen, T.N.; Winther, O.; Brunak, S.; von Heijne, G.; Nielsen, H. SignalP 5.0 improves signal peptide predictions using deep neural networks. Nat. Biotechnol. 2019, 37, 420–423. [Google Scholar] [CrossRef]
  103. Sievers, F.; Higgins, D.G. Clustal Omega for making accurate alignments of many protein sequences. Protein Sci. 2018, 27, 135–145. [Google Scholar] [CrossRef]
  104. Kazutaka, K.; Misakwa, K.; Kei-ichi, K.; Miyata, T. MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30, 3059–3066. [Google Scholar] [CrossRef]
  105. Larsson, A. AliView: A fast and lightweight alignment viewer and editor for large datasets. Bioinformatics 2014, 30, 3276–3278. [Google Scholar] [CrossRef]
  106. Darriba, D.; Posada, D.; Kozlov, A.M.; Stamatakis, A.; Morel, B.; Flouri, T. ModelTest-NG: A new and scalable tool for the selection of DNA and protein evolutionary models. Mol. Biol. Evol. 2020, 37, 291–294. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Bouckaert, R.; Vaughan, T.G.; Barido-Sottani, J.; Duchêne, S.; Fourment, M.; Gavryushkina, A.; Heled, J.; Jones, G.; Kühnert, D.; De Maio, N.; et al. BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis. PLoS Comput. Biol. 2019, 15, e1006650. [Google Scholar] [CrossRef] [PubMed]
  108. Letunic, I.; Bork, P. Interactive tree of life (iTOL) v3: An online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016, 44, 242–245. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Functional annotation in GO, InterPro, and UniprotKB databases (accessed on October 2021) of unigenes obtained from A. ludens head transcriptomes constructed with Illumina HiSeq2000 (a) and BGISEQ-500 (b) platforms; also showing species distribution, similarity percentages, and E-value from BLASTx analysis against the Insecta-UniprotKB database.
Figure 1. Functional annotation in GO, InterPro, and UniprotKB databases (accessed on October 2021) of unigenes obtained from A. ludens head transcriptomes constructed with Illumina HiSeq2000 (a) and BGISEQ-500 (b) platforms; also showing species distribution, similarity percentages, and E-value from BLASTx analysis against the Insecta-UniprotKB database.
Ijms 23 10531 g001
Figure 2. Bayesian phylogenetic analysis of candidate OBPs from Anastrepha ludens with homologs from seven Tephritidae species and the model organism Drosophila melanogaster. The clusters of subfamily OBPs identified in AludOBPs have distinct color bar marks. Gray bars delimit the branch expansions of the AludOBPs and their orthologues. Circles showed posterior probability values higher than 70%.
Figure 2. Bayesian phylogenetic analysis of candidate OBPs from Anastrepha ludens with homologs from seven Tephritidae species and the model organism Drosophila melanogaster. The clusters of subfamily OBPs identified in AludOBPs have distinct color bar marks. Gray bars delimit the branch expansions of the AludOBPs and their orthologues. Circles showed posterior probability values higher than 70%.
Ijms 23 10531 g002
Figure 3. Bayesian phylogenetic analysis of candidate CSPs of Anastrepha ludens with homologs from twelve dipteran species and the model organism Drosophila melanogaster. Blue and gray bars delimit the branch expansions of the AludCSPs and their orthologues. Circles showed posterior probability values over 70%. Letters “U-” and “P-” refer to sequences in the databases that are uncharacterized or were automatically annotated by the genomic sequence prediction method.
Figure 3. Bayesian phylogenetic analysis of candidate CSPs of Anastrepha ludens with homologs from twelve dipteran species and the model organism Drosophila melanogaster. Blue and gray bars delimit the branch expansions of the AludCSPs and their orthologues. Circles showed posterior probability values over 70%. Letters “U-” and “P-” refer to sequences in the databases that are uncharacterized or were automatically annotated by the genomic sequence prediction method.
Ijms 23 10531 g003
Figure 4. Bayesian phylogenetic analysis of candidate SNMPs from Anastrepha ludens with homologs from seven dipteran species and the model organism Drosophila melanogaster. Bars delimit clusters of AludSNMPs and their orthologues. Circles showed posterior probability values higher than 70%.
Figure 4. Bayesian phylogenetic analysis of candidate SNMPs from Anastrepha ludens with homologs from seven dipteran species and the model organism Drosophila melanogaster. Bars delimit clusters of AludSNMPs and their orthologues. Circles showed posterior probability values higher than 70%.
Ijms 23 10531 g004
Figure 5. Bayesian phylogenetic analysis of candidate ORs of Anastrepha ludens with homologs from eight dipteran species and the model organism Drosophila melanogaster. Circles showed posterior probability values higher than 70%. Pink bars delimit clusters of AludORs and their orthologues.
Figure 5. Bayesian phylogenetic analysis of candidate ORs of Anastrepha ludens with homologs from eight dipteran species and the model organism Drosophila melanogaster. Circles showed posterior probability values higher than 70%. Pink bars delimit clusters of AludORs and their orthologues.
Ijms 23 10531 g005
Figure 6. Bayesian phylogenetic analysis of candidate GRs from Anastrepha ludens with homologs from seven Tephritidae species and Drosophila melanogaster. Distinct colors represent clusters of AludGRs members based on putative functions reported in previous studies. Gray bars delimit branch expansions of AludGRs and their orthologues. Circles showed posterior probability values higher than 70%.
Figure 6. Bayesian phylogenetic analysis of candidate GRs from Anastrepha ludens with homologs from seven Tephritidae species and Drosophila melanogaster. Distinct colors represent clusters of AludGRs members based on putative functions reported in previous studies. Gray bars delimit branch expansions of AludGRs and their orthologues. Circles showed posterior probability values higher than 70%.
Ijms 23 10531 g006
Figure 7. Bayesian phylogenetic analysis of candidate iGluRs/IRs from Anastrepha ludens with homologs from seven Tephritidae species and Drosophila melanogaster. Distinct colors represent clusters of AludIRs members based on putative functions reported in previous studies. Gray bars delimit the branch expansions of the AludIRs and their orthologues. Circles showed posterior probability values higher than 70%.
Figure 7. Bayesian phylogenetic analysis of candidate iGluRs/IRs from Anastrepha ludens with homologs from seven Tephritidae species and Drosophila melanogaster. Distinct colors represent clusters of AludIRs members based on putative functions reported in previous studies. Gray bars delimit the branch expansions of the AludIRs and their orthologues. Circles showed posterior probability values higher than 70%.
Ijms 23 10531 g007
Table 1. Attributes of A. ludens Odorant-binding proteins subfamilies.
Table 1. Attributes of A. ludens Odorant-binding proteins subfamilies.
SubfamilyNumberLength (aa)Protein Core RegionSimilarity (%)
Classic17110–306C1-X36–171-C2-X3-C3-X32–50-C4-X8–21-C5-X8-C626.41
Minus-C8102–173C1-X16–30-C2-X38–39-C3-X18–19-C4 and C1-X28-C2-X3-C3-X39–40-C4-X10-C5-X8-C639.26
Plus-C4126–232C1a-X10–13-C1b-C1c-X13–19-C1-X9–48-C2-X3-C3-X43–44-C4-X19–33-C5-X9-C6-X8-C7-X10–11-C8-X9-C928.87
Dimer2271–282C1-X28–34-C2-X3-C3-X31-C4-X10–14-C5-X8-C6-X17–28-C1′-X24–25-C2′-X3-C3′-X35-C4′-X17–21-C5′-X8-C6′35.65
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Segura-León, O.L.; Torres-Huerta, B.; Estrada-Pérez, A.R.; Cibrián-Tovar, J.; Hernandez-Hernandez, F.d.l.C.; Cruz-Jaramillo, J.L.; Meza-Hernández, J.S.; Sánchez-Galicia, F. Identification of Candidate Chemosensory Gene Families by Head Transcriptomes Analysis in the Mexican Fruit Fly, Anastrepha ludens Loew (Diptera: Tephritidae). Int. J. Mol. Sci. 2022, 23, 10531. https://doi.org/10.3390/ijms231810531

AMA Style

Segura-León OL, Torres-Huerta B, Estrada-Pérez AR, Cibrián-Tovar J, Hernandez-Hernandez FdlC, Cruz-Jaramillo JL, Meza-Hernández JS, Sánchez-Galicia F. Identification of Candidate Chemosensory Gene Families by Head Transcriptomes Analysis in the Mexican Fruit Fly, Anastrepha ludens Loew (Diptera: Tephritidae). International Journal of Molecular Sciences. 2022; 23(18):10531. https://doi.org/10.3390/ijms231810531

Chicago/Turabian Style

Segura-León, Obdulia L., Brenda Torres-Huerta, Alan Rubén Estrada-Pérez, Juan Cibrián-Tovar, Fidel de la Cruz Hernandez-Hernandez, José Luis Cruz-Jaramillo, José Salvador Meza-Hernández, and Fabian Sánchez-Galicia. 2022. "Identification of Candidate Chemosensory Gene Families by Head Transcriptomes Analysis in the Mexican Fruit Fly, Anastrepha ludens Loew (Diptera: Tephritidae)" International Journal of Molecular Sciences 23, no. 18: 10531. https://doi.org/10.3390/ijms231810531

APA Style

Segura-León, O. L., Torres-Huerta, B., Estrada-Pérez, A. R., Cibrián-Tovar, J., Hernandez-Hernandez, F. d. l. C., Cruz-Jaramillo, J. L., Meza-Hernández, J. S., & Sánchez-Galicia, F. (2022). Identification of Candidate Chemosensory Gene Families by Head Transcriptomes Analysis in the Mexican Fruit Fly, Anastrepha ludens Loew (Diptera: Tephritidae). International Journal of Molecular Sciences, 23(18), 10531. https://doi.org/10.3390/ijms231810531

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop