<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Physiol.</journal-id>
<journal-title>Frontiers in Physiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Physiol.</abbrev-journal-title>
<issn pub-type="epub">1664-042X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fphys.2020.585883</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Physiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Red Imported Fire Ant (<italic>Solenopsis invicta</italic>) Chemosensory Proteins Are Expressed in Tissue, Developmental, and Caste-Specific Patterns</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Wanchoo</surname> <given-names>Arun</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn002"><sup>&#x2020;</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhang</surname> <given-names>Wei</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="author-notes" rid="fn002"><sup>&#x2020;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/729748/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Ortiz-Urquiza</surname> <given-names>Almudena</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<xref ref-type="author-notes" rid="fn002"><sup>&#x2020;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/79947/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Boswell</surname> <given-names>John</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Xia</surname> <given-names>Yuxian</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/538858/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Keyhani</surname> <given-names>Nemat O.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c002"><sup>&#x002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/57878/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Microbiology and Cell Science, Institute of Food and Agricultural Sciences, University of Florida</institution>, <addr-line>Gainesville, FL</addr-line>, <country>United States</country></aff>
<aff id="aff2"><sup>2</sup><institution>Genetic Engineering Research Center, School of Life Sciences, Chongqing University</institution>, <addr-line>Chongqing</addr-line>, <country>China</country></aff>
<aff id="aff3"><sup>3</sup><institution>Department of Biosciences, College of Science, Swansea University</institution>, <addr-line>Swansea</addr-line>, <country>United Kingdom</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Ya-Nan Zhang, Huaibei Normal University, China</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Herbert Venthur, University of La Frontera, Chile; Yang Liu, Chinese Academy of Agricultural Sciences, China; Sufang Zhang, Chinese Academy of Forestry, China</p></fn>
<corresp id="c001">&#x002A;Correspondence: Almudena Ortiz-Urquiza, <email>a.ortiz-urquiza@swansea.ac.uk</email></corresp>
<corresp id="c002">Nemat O. Keyhani, <email>keyhani@ufl.edu</email></corresp>
<fn fn-type="other" id="fn002"><p><sup>&#x2020;</sup>These authors have contributed equally to this work</p></fn>
<fn fn-type="other" id="fn004"><p>This article was submitted to Invertebrate Physiology, a section of the journal Frontiers in Physiology</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>23</day>
<month>10</month>
<year>2020</year>
</pub-date>
<pub-date pub-type="collection">
<year>2020</year>
</pub-date>
<volume>11</volume>
<elocation-id>585883</elocation-id>
<history>
<date date-type="received">
<day>21</day>
<month>07</month>
<year>2020</year>
</date>
<date date-type="accepted">
<day>22</day>
<month>09</month>
<year>2020</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2020 Wanchoo, Zhang, Ortiz-Urquiza, Boswell, Xia and Keyhani.</copyright-statement>
<copyright-year>2020</copyright-year>
<copyright-holder>Wanchoo, Zhang, Ortiz-Urquiza, Boswell, Xia and Keyhani</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p>The red imported fire ant, <italic>Solenopsis invicta</italic>, is a eusocial invasive insect that has spread worldwide. Chemosensory proteins (CSPs) are ligand-binding proteins that participate in a diverse range of physiological processes that include olfaction and chemical transport. Here, we performed a systematic survey of the expression of the 21 gene <italic>S. invicta CSP</italic> family that includes at least two groups of apparent <italic>S. invicta</italic>-specific gene expansions. These data revealed caste, tissue, and developmental stage-specific differential expression of the <italic>SiCSPs</italic>. In general, moderate to high <italic>SiCSP</italic> expression was seen in worker antennae and abdomen tissues with lower expression in head/thorax regions. Male and female alates showed high antennal expression of fewer <italic>SiCSPs</italic>, with the female alate thorax showing comparatively high <italic>SiCSP</italic> expression. <italic>SiCSP</italic> expression was lower in male alates tissues compared to workers and female alates, albeit with some highly expressed <italic>SiCSPs</italic>. <italic>SiCSP</italic> expression was low during development including in eggs, larvae (early and late instars), and pupae. Global analyses revealed examples of conserved, divergent, and convergent <italic>SiCSP</italic> expression patterns linked to phylogenetic relationships. The developmental and caste-specific variation seen in <italic>SiCSP</italic> expression patterns suggests specific functional diversification of CSPs that may translate into differential chemical recognition and communication among individuals and/or reflect other cellular roles of CSPs. Our results support a model for CSPs acting as general ligand carriers involved in a wide range of physiological processes beyond olfaction. As compared to the expression patterns of the <italic>S. invicta</italic> odorant binding proteins (OBPs), an inverse correlation between <italic>SiOBP</italic> and <italic>SiCSP</italic> expression was seen, suggesting potential complementary and/or compensatory functions between these two classes of ligand carriers.</p>
</abstract>
<kwd-group>
<kwd>social insects</kwd>
<kwd>chemical communication</kwd>
<kwd>ligand binding proteins</kwd>
<kwd>chemosensory proteins</kwd>
<kwd>caste</kwd>
<kwd>tissue</kwd>
<kwd>developmental expression</kwd>
</kwd-group>
<contract-sponsor id="cn001">National Science Foundation<named-content content-type="fundref-id">10.13039/100000001</named-content></contract-sponsor>
<contract-sponsor id="cn002">Institute of Food and Agricultural Sciences<named-content content-type="fundref-id">10.13039/100007861</named-content></contract-sponsor>
<counts>
<fig-count count="6"/>
<table-count count="0"/>
<equation-count count="0"/>
<ref-count count="60"/>
<page-count count="19"/>
<word-count count="0"/>
</counts>
</article-meta>
</front>
<body>
<sec id="S1">
<title>Introduction</title>
<p>Ants are the most numerous of the social insects and one of the most successful animals on earth (<xref ref-type="bibr" rid="B19">Holldobler and Wilson, 1990</xref>). The red imported fire ant, <italic>Solenopsis invicta</italic> (Buren) was originally endemic to South America (Northern Argentina, Southern Brazil, and parts of Paraguay). Sometime between 1933 and 1945, this ant likely found its way aboard ships carrying exotic fruits, lumber, and other goods from South America apparently establishing somewhere around Mobile, Alabama in the United States, from where it has spread to many other parts of the world (<xref ref-type="bibr" rid="B50">Tschinkel, 2006</xref>; <xref ref-type="bibr" rid="B46">Shoemaker et al., 2011</xref>). In this respect, <italic>S. invicta</italic> has demonstrated a remarkable plasticity in environmental adaptations making it one of the most successful invasive animal species on our planet. It is an important pest throughout the southern United States and has recently invaded and established itself in California, South East Asia, and even parts of Europe. <italic>S. invicta</italic> form complex societies that are territorial, and whose members display strong nest-mate recognition, and elaborate task specialization (<xref ref-type="bibr" rid="B50">Tschinkel, 2006</xref>). Their colonies are considered to form a multi-tiered &#x201C;caste&#x201D; system consisting of non-reproductive workers, that display task-differentiation, and reproductive winged males and females. Male and female reproductive forms are produced seasonally as the colony matures, a process that can take up to 5 years. Depending upon their genetic background <italic>S. invicta</italic> colonies can have one (monogyny) or multiple queens (polygyny) (<xref ref-type="bibr" rid="B45">Ross and Fletcher, 1985</xref>; <xref ref-type="bibr" rid="B13">Gotzek and Ross, 2008</xref>). These two differing colony organizations are genetically distinct and have consequent important differences in the biology and responses of <italic>S. invicta</italic>. Originally linked to a single gene, termed <italic>Gp-9</italic> (renamed OBP3), mono- versus polygyne colony organization appears to be linked to a larger 13 Mb non-recombining &#x201C;social&#x201D; chromosomal fragment (&#x223C;55% of the chromosome) containing at least 616 open reading frames (<xref ref-type="bibr" rid="B54">Wurm et al., 2011</xref>; <xref ref-type="bibr" rid="B52">Wang et al., 2013</xref>).</p>
<p>Chemosensory proteins (CSPs) are a protein family defined by amino acid homology and sequence motifs that include four cysteine residues (two disulfide bridges) with conserved spacing and a set of &#x03B1;-helices that form a hydrophobic binding cavity as part of a compact, small molecular weight (10&#x2013;15 kDa) structure (<xref ref-type="bibr" rid="B28">Lartigue et al., 2002</xref>; <xref ref-type="bibr" rid="B59">Zhou et al., 2006</xref>; <xref ref-type="bibr" rid="B42">Pelosi et al., 2014</xref>). CSPs were originally identified as a protein accumulating during leg regeneration in nymphal stages of the cockroach <italic>Periplaneta americana</italic> (<xref ref-type="bibr" rid="B40">Nomura et al., 1992</xref>), other CSPs were subsequently found enriched in sensory organs including the antennae, and hence presumed to function within the context of chemoreception or olfaction by acting as odorant (ligand) carrier proteins (<xref ref-type="bibr" rid="B1">Angeli et al., 1999</xref>; <xref ref-type="bibr" rid="B43">Picimbon et al., 2000</xref>; <xref ref-type="bibr" rid="B5">Ban et al., 2003</xref>). CSPs in the alfalfa plant bug (<italic>Adelphocoris lineolatus</italic>) have been implicated in mediating host recognition (<xref ref-type="bibr" rid="B14">Gu et al., 2012</xref>), and CSP/Takeout genes have been shown to be involved in olfactory-based behaviors including repulsion and attraction in the migratory locust (<italic>Locusta migratoria</italic>) (<xref ref-type="bibr" rid="B16">Guo et al., 2011</xref>). However, it has also been recognized that many CSPs are expressed in a variety of other tissues, e.g., in the pheromone glands of the cabbage moth, <italic>Mamestra brassicae</italic> and the silk moth, <italic>Bombyx mori</italic>, and in male and female reproductive organs, e.g., ejaculatory bulb of <italic>Drosophila melanogaster</italic> and female organs of <italic>L. migratoria manilensis</italic> (<xref ref-type="bibr" rid="B60">Zhou et al., 2013</xref>), where CSPs potentially function in pheromone detection and release (<xref ref-type="bibr" rid="B23">Jacquin-Joly et al., 2001</xref>; <xref ref-type="bibr" rid="B7">Dani et al., 2011</xref>). More recently, various CSPs have been implicated in a wide range of physiological processes beyond olfaction and/or chemical communication. RNAi knockdown of the honeybee (<italic>Apis mellifera</italic>) CSP5 indicated a role for this protein in embryonic integument development (<xref ref-type="bibr" rid="B36">Maleszka et al., 2007</xref>) and CSP3 (designed ASP3c), also from the honeybee, has been shown to act as a brood pheromone carrier protein (<xref ref-type="bibr" rid="B6">Briand et al., 2002</xref>). CSPs, therefore, should be considered as general ligand carrier proteins, some of which may function within the scope of chemical perception mediating binding of volatile and/or hydrophobic odorants, whereas others may ferry hormones and/or other chemical compounds participating in organismal homeostasis and/or developmental process. Participation in pheromone/semiochemical sequestration and/or release would still link CSPs to chemical communication, albeit in the broad sense, namely outside of sensillar mediated olfaction.</p>
<p>Chemosensory proteins are considered to have at least some parallels in function to another class of small molecular weight (&#x223C;130&#x2013;150 amino acids), soluble proteins known as odorant binding proteins (OBPs). OBPs, like CSPs, are thought to mediate signal transduction in insects by shuttling hydrophobic compounds (e.g., odorants) to odorant receptors expressed on the dendritic membrane of the olfactory neurons. Insects OBPs do not share homology to the vertebrate OBPs and are sometimes annotated as pheromone binding proteins (PBPs) or general odorant binding proteins (GOBPs) (<xref ref-type="bibr" rid="B58">Zhang et al., 2016</xref>). Expression profiling of the <italic>S. invicta</italic> 17-gene member OBP family revealed both antennal and non-antennal expression specific <italic>SiOBPs</italic>, as well as caste and tissue specificity in the expression of these proteins (<xref ref-type="bibr" rid="B58">Zhang et al., 2016</xref>). Intriguingly, the <italic>S. invicta</italic> &#x201C;social&#x201D; chromosomal fragment contains multiple members of both the CSP and OBP gene families amongst other genes.</p>
<p>Antenna-specific CSPs, enriched in specialized chemosensory sensilla, have been identified in the Argentine ant (<italic>Linepithema humile</italic>) and the Japanese carpenter ant (<italic>Camponotus japonicus</italic>), with CjCSP1 mediating cuticular hydrocarbon recognition and transport, a critical process for nestmate/non-nestmate discrimination in <italic>C. japonicus</italic> (<xref ref-type="bibr" rid="B22">Ishida et al., 2002</xref>; <xref ref-type="bibr" rid="B41">Ozaki et al., 2005</xref>). Comparative transcriptomics in a variety of ant species also indicated expression of CSPs in non-antennal tissues, further supporting their functioning beyond olfaction (<xref ref-type="bibr" rid="B37">McKenzie et al., 2014</xref>). Initially, fourteen CSP sequences were identified in <italic>S. invicta</italic> expressed sequence tag (EST) libraries, with one, designated SiCSP1(SiCSP19 in our study), as well as an apolipophorin-like protein, found to be highly expressed in the worker antennae (<xref ref-type="bibr" rid="B15">Guntur et al., 2004</xref>; <xref ref-type="bibr" rid="B12">Gonzalez et al., 2009</xref>; <xref ref-type="bibr" rid="B55">Xu et al., 2009</xref>). Subsequent genome mining and phylogenetic analyses revealed intriguing CSP gene expansions in the available ant genomes with the apparent evolution of ant-specific CSP clades (<xref ref-type="bibr" rid="B27">Kulmuni et al., 2013</xref>). Within these analyses 21 CSPs in total were identified in the <italic>S. invicta</italic> genome, representing one of the largest CSP sets annotated in eusocial Hymenoptera to date. Most <italic>Drosophila</italic> spp., by contrast, have only four CSP genes (<xref ref-type="bibr" rid="B51">Vieira and Rozas, 2011</xref>).</p>
<p>Here, the expression pattern of the 21 CSP genes identified in <italic>S. invicta</italic> was examined by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR). <italic>SiCSP</italic> gene expression patterns were examined as a function of caste, tissue distribution, and developmental stage. Developmental specificity of <italic>SiCSP</italic> was examined in eggs, early instar larvae, late instar larvae, and pupae. For adult ant (workers, male and female alates) samples were further subdivided into specific anatomical tissues that included the antenna, head, thorax, and abdomen. High to robust expression of most <italic>SiCSPs</italic> was found in the antennae of workers and male alates, and to a lesser extent, in female alate antennae. <italic>SiCSP</italic> expression was significantly lower in the developmental stages as compared to the adults, although robust expression of specific <italic>SiCSPs</italic> were found in early-instar larvae and pupae. These data revealed dynamic and differential <italic>SiCSP</italic> expression patterns among developmental stages and adult tissues with little correlation to their phylogenetic relatedness, although sub-groups of phylogenetically related <italic>SiCSPs</italic> showed similar patterns of expression. Both divergence between more closely related SiCSPs and convergence between more remotely related SiCSPs was observed. This study provides a basis for a further systematic characterization of the functions of the different CSPs in <italic>S. invicta</italic>.</p>
</sec>
<sec id="S2" sec-type="materials|methods">
<title>Materials and Methods</title>
<sec id="S2.SS1">
<title>Insects and Experimental Samples</title>
<p><italic>Solenopsis invicta</italic> laboratory colonies were collected from the field (Gainesville, FL, United States) and maintained as described in <xref ref-type="bibr" rid="B8">Fan et al. (2012)</xref>. The field colony was assessed to be polygyne with multiple queens evident in the founding colony. In addition, sequencing of the full-length cDNA of <italic>SiOBP3/Gp-9</italic> gene revealed the presence of both <italic>Gp-9B</italic> and <italic>Gp-9b</italic> alleles. Laboratory colonies were maintained at room temperature with &#x223C;70% relative humidity and a 16:8 dark:light photoperiod. Ants were fed with 300 mM sucrose solutions randomly dispersed throughout the trays and supplemented with freeze dried <italic>Galleria mellonella</italic> larvae. Dissections into four sections; antenna, head (without antennae), thorax and abdomen were performed using separately processed (&#x223C;200&#x2013;500 ants) workers, male and female alates that were immersed in RNALater (Invitrogen, Thermo Fisher Scientific, Waltham, MA, United States) and dissected under a stereomicroscope. Adult stages were not sampled as same-age cohorts. Eggs, collected from mated queens within 24 h, as well as small larvae, large larvae and pupa (approximately 200&#x2013;500 each) were immediately suspended in RNALater and stored at &#x2212;80&#x00B0;C until RNA extraction. Larvae were distinguished by their size and reflected a mixture of minor and major workers that could include male and female alates. Early instars (1st and 2nd) were selected based on size (&#x2264;1 mm) and color, as these often displayed some melanization, whereas late instars (3rd and 4th) were larger (1.1&#x2013;4.5 mm) and clear to whitish.</p>
</sec>
<sec id="S2.SS2">
<title>RNA Preparation and cDNA Library Construction</title>
<p>Samples (&#x223C;100 mg) were ground in liquid nitrogen using a mortar and pestle, after which Trizol reagent (1 ml) was added, and total RNA extracted following the manufacturers&#x2019; protocols (Invitrogen, Carlsbad, CA, United States). Genomic DNA in samples were digested using TURBODNase (Invitrogen, Carlsbad, CA, United States). Total RNA quality and quantity were analyzed by agarose gel electrophoresis and via NanoDrop 2000 spectrophotometric analyses. Quantification of RNA concentrations in samples was performed using a Qubit H 2.0 fluorometer (Invitrogen, Carlsbad, CA, United States) and cDNA libraries were constructed using 2 mg total RNA using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA, United States). At least three biological replicates were prepared for each sample.</p>
</sec>
<sec id="S2.SS3">
<title>Absolute Qualification of cDNAs</title>
<p>The amino acid sequences of the <italic>S. invicta</italic> chemosensory proteins (SiCSPs) were downloaded from NCBI (<xref ref-type="supplementary-material" rid="TS1">Supplementary Table S1</xref>). Absolute qualification of the <italic>SiCSPs</italic> were performed as previous described (<xref ref-type="bibr" rid="B53">Whelan et al., 2003</xref>). Briefly, target-specific qRT-PCR primers were designed using the Beacon designer 8.13 software program (Palo Alto, CA, United States) and synthesized commercially (Invitrogen, Carlsbad, CA, United States). The list of primers used is given in <xref ref-type="supplementary-material" rid="TS2">Supplementary Table S2</xref>. Primers were used in PCR reactions to amplify <italic>SiCSP</italic> target sequences using an <italic>S. invicta</italic> cDNA library constructed as described above as the template. After purification, PCR fragments were cloned into the pGEMT vector (Promega Corp., Madison, WI, United States). Positive clones were isolated, and the integrity of the inserts verified by sequencing (Eton Biosciences, San Diego, CA, United States). Concentrations of the pGEMT plasmids for each <italic>SiCSP</italic> were quantified using Qubit dsDNA BR Assay Kit (Invitrogen, Carlsbad, CA, United States) for use in constructing absolute quantification standard curves. RT-PCR primer sets for each <italic>SiCSP</italic> were validated for production of the correct amplicon size, optimal ratios of primers, <italic>T</italic><sub>m</sub>, and efficiency (<xref ref-type="supplementary-material" rid="TS3">Supplementary Table S3</xref>). The amplification efficiency (E) was calculated using the slope of a linear regression determined by the Ct values (<italic>Y</italic>-axis) and the log<sub>10</sub> concentration of the cDNA (<italic>X</italic>-axis). Slopes were used to calculated Efficiencies (<italic>E</italic>) using the formula: <italic>E</italic> = 10<sup>1/slope</sup>&#x2212;1. Plasmid constructs and optimized PCR conditions were used to acquire data for construction of absolute expression standard curves, using serial dilution of the plasmid templates (10<sup>&#x2013;5</sup> to 10<sup>&#x2013;9</sup>). The number of transcript copies was calculated using the molecular weight of plasmids and their empirically determined concentrations using the formula: copy number = (6.02 &#x00D7; 10<sup>23</sup>) &#x00D7; (amount ng &#x00D7; 10<sup>&#x2013;9</sup>)/(DNA length &#x00D7; 660).</p>
<p>Individual standard curves and empirically determined C<sub>t</sub> values (as derived from the qRT-PCR experiments) were used to calculate <italic>SiCSP</italic> absolute transcript expression values. qRT-PCR reactions were performed using 2XSYBR Green qPCR Master Mix (Biotools, Houston, TX, United States). In most cases, <italic>S. invicta</italic> cDNA libraries were diluted 40-fold in sterile RNAase-free H<sub>2</sub>O. Typically, reactions included: 2XMaster Mix (7.5 ml), 5 &#x03BC;L template (12.5 ng of cDNA) and 200 nM of each gene specific primer in a total volume of 15 ml. qRT-PCR reactions were performed for at least three RNA preparations (i.e., three independent biological samples) from each tissue sample. PCR reactions were performed using the Eco Real-Time qPCR System (Illumina, San Diego, CA, United States) with a thermo-profile of one cycle of 95&#x00B0;C 5 min, 95&#x00B0;C 2 min, then 45 cycles of 95&#x00B0;C 15 s, and 60/59&#x00B0;C 45 s, followed by a melting curve analysis from 55 to 95&#x00B0;C. Absolute quantification of the <italic>S. invicta</italic> elongation factor-&#x03B1; (<italic>EF1</italic>&#x03B1;) and glyceraldehyde 6-phosphate dehydrogenase (<italic>GAPDH</italic>) genes of <italic>S. invicta</italic> were used in normalization analyses, and the geometric mean of the absolute number of transcripts of <italic>EF1</italic>&#x03B1; and <italic>GAPDH</italic> in every examined tissue was used to obtain the values of normalized relative mRNA transcript abundance.</p>
</sec>
<sec id="S2.SS4">
<title>Data Analysis</title>
<p>Analysis of variance and multivariate analysis of the variance (ANOVA and MANCOVA) were performed to compare <italic>SiCSP</italic> expression levels among adult tissues and across developmental stages. The data, normalized relative mRNA transcript abundance, were transformed logarithmically (Log<sub>10</sub>) in order to correct normality and unequal variances. Means were compared with Tukey&#x2019;s HSD (honestly significant differences) test. Statistical analyses were performed using the IBM SPSS Statistics, version 26 (Armonk, NY, United States). Expression profiles for eggs, larvae and pupae were assessed using one-way ANOVA with <italic>post hoc</italic> comparisons using Tukey&#x2019;s HSD test within and across developmental stages. CSP expression in adults were analyzed via MANCOVA (two variables: caste and tissue) and Tukey&#x2019;s HSD test was performed for each examined tissue among castes. Co-expression patterns of CSPs genes among developmental stages and among adult tissues were inferred by unsupervised hierarchical clustering of the log-transformed normalized expression data. Clustering analyses were implemented by rows (i.e., <italic>CSP</italic>) and columns (i.e., tissue/developmental stage) using the web-based tool Morpheus<sup><xref ref-type="fn" rid="footnote1">1</xref></sup>. The similarity between the objects in the matrix was assessed with the one minus Pearson&#x2019;s correlation metric and the complete linkage approach.</p>
</sec>
<sec id="S2.SS5">
<title>Phylogenetic and Motif Analyses</title>
<p>The amino acid sequences of a total of 92 CSPs from <italic>S. invicta</italic>, <italic>Acromyrmex echinatior</italic>, <italic>C. floridanus</italic>, <italic>C. japonicus</italic>, <italic>Harpegnathos saltator</italic>, <italic>D. melanogaster</italic>, <italic>D. grimshawi</italic>, <italic>Apis mellifera</italic>, and <italic>Polistes canadensis</italic> were used for phylogenetic tree construction and motif analyses (<xref ref-type="supplementary-material" rid="TS4">Supplementary Table S4</xref>). Putative N-terminal signal peptides were identified using Signal P<sup><xref ref-type="fn" rid="footnote2">2</xref></sup> (<xref ref-type="supplementary-material" rid="TS1">Supplementary Table S1</xref>). Exon-intron splice positions were identified using the online Gene Structure Display Server (GSDS<sup><xref ref-type="fn" rid="footnote3">3</xref></sup>). Amino acid motifs were identified using MEME (version 4.11.2<sup><xref ref-type="fn" rid="footnote4">4</xref></sup>). The amino acid multiple sequence alignment (MSA) was generated with PRANK (<xref ref-type="bibr" rid="B33">L&#x00F6;ytynoja and Goldman, 2005</xref>, <xref ref-type="bibr" rid="B34">2010</xref>). The best fitting model of amino acid substitution was estimated with MEGA 6.0 (<xref ref-type="bibr" rid="B49">Tamura et al., 2013</xref>). MEGA optimizes the tree topology search starting with a Neighbor Joining tree and uses the likelihood function and three model criteria BIC (Bayesian Information Criterion), AIC (Akaike information criterion) and LnL (log likelihood) to find the best fitting model of amino acid substitution, which was: LG + G (G = Gamma shape parameters). The phylogenetic tree was built using RaxML at the CIPRES Science Gateway (<xref ref-type="bibr" rid="B38">Miller et al., 2010</xref>; <xref ref-type="bibr" rid="B48">Stamatakis, 2014</xref>) and the LG amino acid substitution model. G was estimated, branch lengths optimized, and branch support calculated by bootstrapping. RaxML was allowed to execute 1000 rapid bootstrap inferences and halt bootstrapping automatically after a thorough maximum likelihood search (893 bootstraps). The software MEGA 6.0 was used to draw the tree.</p>
</sec>
</sec>
<sec id="S3">
<title>Results</title>
<sec id="S3.SS1">
<title>The Chemosensory Protein Repertoire of <italic>Solenopsis invicta</italic></title>
<p>A set of 21 <italic>CSP</italic> genes has been previously identified in the genome of <italic>S. invicta</italic> (<xref ref-type="supplementary-material" rid="TS1">Supplementary Table S1</xref>). A limited protein sequence based phylogenetic tree illustrates the division of these CSPs into two discrete major branches, one containing those showing orthology to sequences found in other insects, i.e., CSPs 8, 7, 2, 6, 3, and 4 (boxed in red), which we henceforth term &#x201C;general CSPs,&#x201D; and those found essentially only in ants albeit with some exceptions: &#x201C;ant CSP expansion&#x201D; (<xref ref-type="fig" rid="F1">Figure 1</xref>). Within the ant expansion grouping, two <italic>S. invicta</italic>-specific gene expansions were apparent; namely one consisting of <italic>SiCSPs 19</italic> [identified as a major antennal SiCSP (<xref ref-type="bibr" rid="B12">Gonzalez et al., 2009</xref>)], <italic>9</italic>, <italic>10</italic>, <italic>12</italic>, <italic>13</italic>, <italic>20</italic> and <italic>22</italic> (boxed in green in <xref ref-type="fig" rid="F1">Figure 1</xref>) and another including <italic>SiCSPs 11</italic>, <italic>15</italic>, <italic>16</italic>, <italic>17</italic> and <italic>21</italic> (boxed in blue). Three additional <italic>SiCSPs</italic>: <italic>1</italic>, <italic>14</italic> and <italic>18</italic> (boxed in brown) were also found distributed within the ant expansion group forming separate sub-clades. Of note, SiCSP1 showed high similarity (83.7%) to the cuticular hydrocarbon recognition CSP identified in the Japanese carpenter ant (<italic>C. japonicus</italic>), labeled as CjCSP1 in <xref ref-type="fig" rid="F1">Figure 1</xref> (<xref ref-type="bibr" rid="B41">Ozaki et al., 2005</xref>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption><p>Phylogenetic analysis of <italic>S. invicta</italic> chemosensory proteins (SiCSPs). Limited maximum likelihood phylogeny of <italic>S. invicta</italic> CSPs compared to the CSP repertoires found in the ant species, <italic>A. echinatior</italic> (AeCSP), <italic>C. floridanus</italic> (CfCSP), <italic>C. japonicus</italic> (CjCSP), and <italic>H. saltator</italic> (HsCSP), the European honeybee <italic>Apis mellifera</italic>, the primitive eusocial wasp <italic>Polistes canadensis</italic> and in the fruit flies <italic>D. melanogaster</italic> (DmCSP) and <italic>D. grimshawi</italic> (DgCSP) (accession numbers given in <xref ref-type="supplementary-material" rid="TS4">Supplementary Table S4</xref>). Numbers at nodes indicate bootstrap values. The tree is midpoint-rooted in the absence of a suitable out-group.</p></caption>
<graphic xlink:href="fphys-11-585883-g001.tif"/>
</fig>
<p>Mapping of the <italic>SiCSP</italic> sequences to the <italic>S. invicta</italic> genome revealed their clustering into seven of the sixteen linkage groups reported by <xref ref-type="bibr" rid="B52">Wang et al. (2013)</xref> (<xref ref-type="fig" rid="F2">Figure 2A</xref>, CSPs are color coded as in <xref ref-type="fig" rid="F1">Figure 1</xref>). Of the general <italic>SiCSPs</italic>, <italic>SiCSPs 3</italic> and <italic>4</italic> (on same linkage group), <italic>SiCSP2</italic>, <italic>SiCSP8</italic>, and <italic>SiCSP7</italic> were located on linkage groups (Lg-) 1, 3, 8, and 16 respectively. <italic>SiCSP7</italic> was found on the same linkage group (Lg16) as the member of the ant expansion <italic>SiCSP14</italic> and four representative of one of the fire ant specific CSP gene expansion (boxed in green), namely <italic>SiCSPs 13</italic>, <italic>12</italic>, <italic>22</italic> and <italic>9</italic>. The remaining members of the latter fire ant specific expansion, i.e., <italic>SiCSPs19</italic>, <italic>10</italic> and <italic>20</italic>, were localized to linkage group Lg4. Four (out of five) members of the second fire ant specific CSP gene expansion (boxed in blue), namely, <italic>SiCSPs 17</italic>, <italic>16</italic>, <italic>15</italic>, and <italic>11</italic> were clustered on linkage group Lg9, along with the two remaining members of the ant CSP expansion, <italic>CSPs 1</italic> and <italic>18</italic>. The remaining member of the (blue boxed) fire ant specific CSP gene expansion, <italic>CSP21</italic>, was localized to linkage group Lg11. <italic>SiCSP6</italic> could not be assigned to any of the 16 linkage groups due to incomplete assembly.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption><p><bold>(A)</bold> Mapping of <italic>SiCSP</italic> nucleotide sequences to 16 linkage groups corresponding to the 16 <italic>S. invicta</italic> chromosomes (<xref ref-type="bibr" rid="B52">Wang et al., 2013</xref>). <italic>SiCSP</italic> positions are not drawn to scale but indicate approximate locations on the linkage groups. <italic>SiCSP</italic> exon boxes are color coded to match their phylogenetic placement as in <xref ref-type="fig" rid="F1">Figure 1</xref>. <bold>(B)</bold> Intron&#x2013;exon structures of <italic>SiCSPs</italic>.</p></caption>
<graphic xlink:href="fphys-11-585883-g002.tif"/>
</fig>
<p>The intron/exon structure of the <italic>SiCSPs</italic>, was determined using the online Gene Structure Display Server (GSDS, see text footnote 3)<underline>,</underline> and correlated with <italic>SiCSP</italic> phylogenetic relatedness (<xref ref-type="fig" rid="F2">Figure 2B</xref>). Members of the general <italic>SiCSPs</italic> (<italic>2</italic>, <italic>3</italic>, <italic>4</italic>, <italic>6</italic>, <italic>7</italic>, and <italic>8</italic>) showed significant variation in intron/exon structure (the intron/exon structure of <italic>SiCSP8</italic> could not be definitively assigned due to ambiguities in the published sequence). Aside from <italic>SiCSP7</italic>, which contained two introns, all other <italic>SiCSPs</italic> contained a single intron. The overall genomic sequences of the general <italic>SiCSPs</italic> ranged from &#x223C;0.75 to &#x003E;7 kb (the open reading frame sequence for each <italic>SiCSP</italic> is &#x223C;300&#x2013;330 bp coding for proteins &#x223C;100&#x2013;110 amino acids in length). <italic>SiCSPs 1</italic>, <italic>14</italic>, and <italic>18</italic> (distributed within the ant CSP expansion group) also varied in intron/exon structure. However, the overall genomic context ranged from 1 to 2.5 kb. The fire ant specific clade consisting of <italic>SiCSPs 9, 10, 12</italic>, <italic>13</italic>, <italic>19</italic>, <italic>20</italic>, and <italic>22</italic> were very similar in intron/exon structure, ranging in size from only &#x223C;0.7 to 1.2 kb. Within the subgroup of <italic>SiCSPs 10</italic>, <italic>19</italic>, and <italic>20</italic>, which clustered together in the phylogenetic tree and were also found tandemly arrayed on linkage group Lg4, little variation in intron/exon structure was seen. Similarly, the variation of intron/exon structure seen for the second fire ant CSP gene expansion (<italic>SiCSPs 11</italic>, <italic>15</italic>, <italic>16</italic>, <italic>17</italic>, and <italic>21</italic>) was also low, ranging from &#x223C;1.1 to 2.5 kb, with <italic>SiCSP11</italic> and <italic>17</italic> nearly identical.</p>
</sec>
<sec id="S3.SS2">
<title>Chemosensory Protein Expression in <italic>S. invicta</italic> Adults and Developmental Stages</title>
<p>Oligonucleotide primers were validated for quantitative RT-PCR as described in the &#x201C;Materials and Methods&#x201D; section. Primer efficiency and amplification of single bands corresponding to the predicted size of each amplicon were verified (<xref ref-type="supplementary-material" rid="TS2">Supplementary Tables S2</xref>, <xref ref-type="supplementary-material" rid="TS3">S3</xref>). All amplicons were cloned and used to construct standard curves for absolute quantification of each respective transcript number in the cDNA samples. In addition to primers designed to the 21 <italic>SiCSPs</italic>, primers were also designed to two different &#x201C;housekeeping&#x201D; genes for use as references. These included primers targeting transcripts for glyceraldehyde-6-phosphate dehydrogenase (<italic>GAPDH</italic>) and elongation factors 1a (<italic>EF1</italic>&#x03B1;). Multivariate analysis of the variance of the normalized gene expression data showed that the examined tissues and castes, as well as the interaction between both independent variables, had a significant effect on the expression of <italic>S. invicta CSPs</italic> (<italic>P</italic> &#x003C; 0.001).</p>
<p>In order to obtain a clearer picture of the expression patterns of the <italic>SiCSPs</italic>, two separate analyses are presented for the adult ants: (A) the absolute number of transcripts as the fractions in each tissue relative to the total concentration of each <italic>SiCSP</italic> across all of the examined tissues (<xref ref-type="fig" rid="F3">Figure 3</xref>) and (B) the normalized relative transcript abundance (with respect to the geometric mean of <italic>EF1&#x03B1;</italic> and <italic>GAPDH</italic> expression levels) (<xref ref-type="fig" rid="F4">Figure 4</xref>). The former analysis gives the relative tissue distribution of each <italic>SiCSP</italic> and is a reflection of <italic>SiCSP</italic> absolute abundance in each tissue relative to the total concentration of the <italic>SiCSP</italic> in workers, female, or male alates. For the developmental stages only (B) the normalized relative transcript abundance in the four examined developmental stages (eggs, 1st&#x2013;2st instar larvae, 3rd&#x2013;4th instar larvae, and pupae) are shown (<xref ref-type="fig" rid="F5">Figure 5</xref>).</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption><p><italic>SiCSP</italic> relative expression distribution across tissues in <bold>(A)</bold> workers, <bold>(B)</bold> male, and <bold>(C)</bold> female alates. The percent of individual <italic>SiCSP</italic> expression in each tissue was calculated as follows: [<italic>SiCSP</italic> expression in tissue]/[total <italic>SiCSP</italic> expression across all four tissues] &#x00D7; 100. Data are representative of at least two independent preparations. Clusters indicate antenna-, head-, thorax-, and abdomen- biased <italic>SiCSPs</italic>. <italic>SiCSPs</italic> are color coded to match their phylogenetic relation.</p></caption>
<graphic xlink:href="fphys-11-585883-g003.tif"/>
</fig>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption><p>Expression of <italic>SiCSPs</italic> in fire ant workers, males and female alates. <bold>(A)</bold> Normalized <italic>SiCSP</italic> expression levels in the worker, male and female antennae. <bold>(B)</bold> Normalized <italic>SiCSP</italic> expression levels in the worker, male and female head. <bold>(C)</bold> Normalized <italic>SiCSP</italic> expression levels in the worker, male and female thorax. <bold>(D)</bold> Normalized <italic>SiCSP</italic> expression levels in the worker, male and female abdomen. Data are representative of at least three independent preparations. Clusters indicate co-expression networks of <italic>SiCSPs</italic> across the examined tissues. Within each row (i.e., <italic>SiCSP</italic>), different letters indicate significant differences at <italic>P</italic> &#x003C; 0.05. The following terms are broadly used to describe the results in terms of number of normalized relative mRNA transcripts (NRTs): very low expression &#x2264;10<sup>&#x2013;4</sup> to 10<sup>&#x2013;7</sup> NRTs, slightly low to low &#x2264;10<sup>&#x2013;2</sup> to &#x003E;10<sup>&#x2013;4</sup> NRTs, moderately robust &#x003C;1 to &#x003E;10<sup>&#x2013;2</sup> NRTs, robust to very robust &#x2265;1 to &#x003C;10 NRTs, slightly high to high expression &#x2265;10 to &#x003C;10<sup>3</sup> NRTs, and very high expression &#x2265;10<sup>3</sup> to 10<sup>5</sup> NRTs.</p></caption>
<graphic xlink:href="fphys-11-585883-g004.tif"/>
</fig>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption><p>Expression of <italic>SiCSPs</italic> in fire ant eggs, larvae and pupae. Normalized <italic>SiCSP</italic> expression levels in 1st&#x2013;2nd instar larvae, 3rd&#x2013;4th instar larvae, eggs and pupae. Data are representative of at least two independent preparations. Clusters indicate co-expression networks of <italic>SiCSPs</italic> across the examined tissues. Within each row (i.e., <italic>SiCSP</italic>), different letters indicate significant differences at <italic>P</italic> &#x003C; 0.05. The following terms are broadly used to describe the results in terms of number of normalized relative mRNA transcripts (NRTs): very low expression &#x2264;10<sup>&#x2013;4</sup> to 10<sup>&#x2013;7</sup> NRTs, slightly low to low &#x2264;10<sup>&#x2013;2</sup> to &#x003E;10<sup>&#x2013;4</sup> NRTs, moderately robust &#x003C;1 to &#x003E;10<sup>&#x2013;2</sup> NRTs, robust to very robust &#x2265;1 to &#x003C;10 NRTs, slightly high to high expression &#x2265;10 to &#x003C;10<sup>3</sup> NRTs, and very high expression &#x2265;10<sup>3</sup> to 10<sup>5</sup> NRTs.</p></caption>
<graphic xlink:href="fphys-11-585883-g005.tif"/>
</fig>
<sec id="S3.SS2.SSS1">
<title>Relative Tissue Distribution of <italic>S. invicta</italic> Chemosensory Proteins in Workers, and Male and Female Alates</title>
<p>Analyses of the tissue distributions of the <italic>SiCSPs</italic>, i.e., <inline-formula><mml:math id="INEQ8"><mml:mfrac><mml:msub><mml:mrow><mml:mo stretchy="false">[</mml:mo><mml:mrow><mml:mi>S</mml:mi><mml:mi>I</mml:mi><mml:mi>C</mml:mi><mml:mi>S</mml:mi><mml:mi>P</mml:mi><mml:mi>X</mml:mi></mml:mrow><mml:mo stretchy="false">]</mml:mo></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mi>u</mml:mi><mml:mi>m</mml:mi><mml:mi>b</mml:mi><mml:mi>e</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>r</mml:mi></mml:mpadded><mml:mi>o</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>f</mml:mi></mml:mpadded><mml:mi>t</mml:mi><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>n</mml:mi><mml:mi>s</mml:mi><mml:mi>c</mml:mi><mml:mi>r</mml:mi><mml:mi>i</mml:mi><mml:mi>p</mml:mi><mml:mi>t</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>s</mml:mi></mml:mpadded><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>r</mml:mi></mml:mpadded><mml:mi>u</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>g</mml:mi></mml:mpadded><mml:mi>o</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>f</mml:mi></mml:mpadded><mml:mi>R</mml:mi><mml:mi>N</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>A</mml:mi></mml:mpadded><mml:mi>i</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>n</mml:mi></mml:mpadded><mml:mi>s</mml:mi><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>c</mml:mi><mml:mi>i</mml:mi><mml:mi>f</mml:mi><mml:mi>i</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>c</mml:mi></mml:mpadded><mml:mi>t</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi><mml:mi>u</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mo stretchy="false">[</mml:mo><mml:mrow><mml:mi>S</mml:mi><mml:mi>I</mml:mi><mml:mi>C</mml:mi><mml:mi>S</mml:mi><mml:mi>P</mml:mi><mml:mi>X</mml:mi></mml:mrow><mml:mo stretchy="false">]</mml:mo></mml:mrow><mml:mrow><mml:mi>t</mml:mi><mml:mi>o</mml:mi><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>l</mml:mi></mml:mpadded><mml:mi>n</mml:mi><mml:mi>u</mml:mi><mml:mi>m</mml:mi><mml:mi>b</mml:mi><mml:mi>e</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>r</mml:mi></mml:mpadded><mml:mi>o</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>f</mml:mi></mml:mpadded><mml:mi>t</mml:mi><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>n</mml:mi><mml:mi>s</mml:mi><mml:mi>c</mml:mi><mml:mi>r</mml:mi><mml:mi>p</mml:mi><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>s</mml:mi></mml:mpadded><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>r</mml:mi></mml:mpadded><mml:mi>u</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>g</mml:mi></mml:mpadded><mml:mi>o</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>f</mml:mi></mml:mpadded><mml:mi>R</mml:mi><mml:mi>N</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>A</mml:mi></mml:mpadded><mml:mi>a</mml:mi><mml:mi>c</mml:mi><mml:mi>r</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>s</mml:mi></mml:mpadded><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mpadded width="+2.8pt"><mml:mi>l</mml:mi></mml:mpadded><mml:mi>t</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi><mml:mi>u</mml:mi><mml:mi>e</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mfrac></mml:math></inline-formula>, coupled with hierarchical clustering analyses allowed for the identification of groups of <italic>SiCSPs</italic> with similar patterns of transcript abundance across tissues in workers, males and female alates (<xref ref-type="fig" rid="F3">Figure 3</xref>). These analyses also revealed that some CSPs appeared to be preferentially expressed in some tissues, with these tissue-biased expression patterns differing across castes and being somewhat phylogenetically related (<xref ref-type="fig" rid="F3">Figure 3</xref>). <italic>SiCSPs</italic> preferentially expressed in the antennae of workers, males and females tended to belong to the group of &#x2018;general CSPs&#x2019; (boxed in red in <xref ref-type="fig" rid="F1">Figure 1</xref>), except for <italic>SiCSP19</italic> in workers, males and females, and <italic>SiCSP11</italic> in males, which were members of the <italic>S. invicta</italic> expansions (green and blue boxes in <xref ref-type="fig" rid="F1">Figure 1</xref>). On the other hand, <italic>SiCSPs</italic> exclusively found in <italic>S. invicta</italic> (green and blue boxes in <xref ref-type="fig" rid="F1">Figure 1</xref>) were seen most abundant in the worker abdomen and the male and female thorax.</p>
<p>In workers, <italic>SiCSP</italic> expression tended to be higher in the antennae and abdomen than in the head and the thorax (<xref ref-type="fig" rid="F3">Figure 3A</xref>). <italic>SiCSPs 6</italic>, <italic>2</italic>, <italic>8</italic> and <italic>19</italic> were found almost exclusively in the worker antennae (&#x003E;91.7% abundance), with <italic>SiCSP1</italic> also mainly distributed in the antennae but in lesser amounts (53.8%). In worker ants, abdomen-biased <italic>SiCSPs</italic> included <italic>4</italic>, <italic>3</italic>, <italic>13</italic>, <italic>14</italic>, <italic>7</italic>, <italic>22</italic>, <italic>10</italic>, <italic>9</italic>, <italic>15</italic>, <italic>17</italic> and <italic>16</italic> (56.0&#x2013;95.8% abundance). <italic>SiCSP18</italic> and <italic>11</italic> were seen most abundant in the worker head, while no thorax-biased <italic>SiCSPs</italic> were found in workers. Notably, <italic>SiCSP 12</italic>, <italic>20</italic> and <italic>21</italic> showed comparable transcript abundance across all the examined tissues in workers.</p>
<p>To a lesser extent, the tissue distribution of <italic>SiCSPs</italic> in males mainly favored thorax tissues although subsets of antenna-biased, head-biased and abdomen-biased <italic>SiCSPs</italic> were also identified in males (<xref ref-type="fig" rid="F3">Figure 3B</xref>). Transcript abundance in the male thorax of <italic>SiCSPs 12, 3, 22, 16, 21, 9, 17, 8</italic> and <italic>13</italic> varied between 44.7% and 66.8%. <italic>SiCSPs 18, 7, 14</italic> and <italic>15</italic> exhibited preferences for the male abdomen (46.2%&#x2013;92.2% abundance), while <italic>SiCSPs 2, 6, 19, 11</italic> and <italic>4</italic> appeared most found in the male antennae, with transcript abundance between 34.9% and 76.2%. Conversely, <italic>SiCSPs 20, 1</italic> and <italic>10</italic> were most expressed in the head, with <italic>SiCSP1</italic> found almost solely in the male head (99.9% abundance).</p>
<p>Similar to males, the fractional tissue distribution of the <italic>SiCSPs</italic> in females showed a clear preference for the thorax, with 15 out of 21 <italic>SiCSPs</italic> being thorax-biased (<xref ref-type="fig" rid="F3">Figure 3C</xref>). Notable exceptions to this were <italic>SiCSPs 19, 8, 1</italic> and <italic>4</italic>, that were preferentially found in the female antennae, showing transcript abundance between 70.8 and 98.0%. In addition, <italic>SiCSP8</italic> was found mainly in the head (66.9% abundance), and <italic>SiCSP2</italic> was found mostly in the female abdomen (72.6% abundance).</p>
</sec>
<sec id="S3.SS2.SSS2">
<title>Chemosensory Protein Expression in Antennae of <italic>S. invicta</italic> Workers, and Male and Female Alates</title>
<p>In order to better present the results, the following terms are broadly used to describe the results in terms of number of normalized relative mRNA transcripts (NRTs): <italic>very low</italic> expression &#x2264;10<sup>&#x2013;4</sup> to 10<sup>&#x2013;7</sup> NRTs, <italic>slightly low to low</italic> &#x2264;10<sup>&#x2013;2</sup> to &#x003E;10<sup>&#x2013;4</sup> NRTs, <italic>moderately robust</italic> &#x003C;1 to &#x003E;10<sup>&#x2013;2</sup> NRTs, <italic>robust to very robust</italic> &#x2265;1 to &#x003C;10 NRTs, <italic>slightly high to high expression</italic> &#x2265;10 to &#x003C;10<sup>3</sup> NRTs, and <italic>very high expression</italic> &#x2265;10<sup>3</sup> to 10<sup>5</sup> NRTs. Overall <italic>SiCSP</italic> expression in the antennae across workers, males and females ranged between robust and very high, with the exception of some <italic>SiCSPs</italic> in female antennae whose expression varied between slightly low and low (<xref ref-type="fig" rid="F4">Figure 4A</xref>). Hierarchical clustering by caste (i.e., columns) revealed that the overall <italic>SiCSP</italic> expression in workers was more similar to female alates than to males, while hierarchical clustering by <italic>SiCSP</italic> (rows) yielded four antennal co-expression clusters, namely groupings of CSPs with similar expression patterns. These clusters of co-expressed genes consisted of subgroups of <italic>SiCSPs</italic> with no apparent phylogenetic relationship: I) <italic>SiCSPs 12, 8</italic> and <italic>17</italic>, II) <italic>SiCSPs 9, 19</italic> and <italic>11</italic>, III) <italic>SiCSPs 18, 10, 15, 7</italic> and <italic>1</italic>, and IV) <italic>SiCSPs 22, 2, 6, 16, 3, 20, 21, 14, 4</italic>, and <italic>13</italic> (<xref ref-type="fig" rid="F4">Figure 4A</xref>). In cluster I, <italic>SiCSP12</italic> showed high expression in males and workers and very high expression in female alates, although these differences were not statistically significant. <italic>SiCSP8</italic>, which was differentially expressed across castes, displayed robust expression in males, high expression in workers and very high expression in female alates. On the other hand, <italic>SiCSP17</italic> exhibited similar expression levels among castes although expression of this gene was lower than that seen for <italic>SiCSPs 12</italic> and <italic>8</italic>.</p>
<p>Within the antennal co-expression cluster II, <italic>SiCSP19</italic> and <italic>11</italic> were equally expressed across caste, although <italic>SiCSP11</italic> showed overall greater expression (i.e., high) than <italic>SiCSP19</italic>, which was only robustly expressed in adult antennae. Expression of <italic>SiCSP19</italic> in the antennae varied significantly between males and workers and males and females alates, ranging between slightly high (males) and high (workers and female alates). Expression levels in the antennal co-expression cluster III fluctuated between slightly low to very high. However, expression values within the same <italic>SiCSP</italic> remained not statistically distinct among male, worker and female antennae for the majority of the <italic>SiCSPs</italic> in grouping III (i.e., <italic>SiCSPs 18, 10, 7</italic> and <italic>1</italic>). Only <italic>SiCSP15</italic>, which was high to very highly expressed in adult antennae, showed significantly different expression among castes. Antennal co-expression cluster IV included ten <italic>SiCSPs</italic> which included members of the general, ant expansion and <italic>S. invicta</italic>-specific <italic>CSPs</italic>. <italic>SiCSP</italic> expression in male and worker antennae within grouping IV was overall robust, while the expression of <italic>SiCSPs 22, 2, 6, 16, 20, 14, 4</italic> and <italic>13</italic> in female antennae varied between low (<italic>SiCSP2</italic>) and slightly robust (<italic>SiCSPs 16</italic> and <italic>4</italic>). Within this group of co-expressed genes, only <italic>SiCSPs 3</italic> and <italic>21</italic> showed no significant differences among castes, while the rest of the <italic>SiCSPs</italic> within grouping IV was significantly lower expressed in female antennae than in male and worker antennae.</p>
</sec>
<sec id="S3.SS2.SSS3">
<title>Chemosensory Protein Expression in the Head of <italic>S. invicta</italic> Workers, and Male and Female Alates</title>
<p>In general, expression of the <italic>S. invicta CSPs</italic> in the adult head was lower than in adult antennae and was similar between workers and female alates. Furthermore, the overall expression of <italic>SiCSPs</italic> in the male head was somewhat lower than in the other two examined castes (<xref ref-type="fig" rid="F4">Figure 4B</xref>). Hierarchical clustering generated five clusters of co-expressed <italic>SiCSPs</italic> which included subgroups of phylogenetically similar and more distant <italic>SiCSP</italic>s. Co-expression cluster I incorporated <italic>SiCSPs 19</italic> and <italic>8</italic>, which differed in expression pattern across castes. <italic>SiCSP19</italic> expression was slightly low in the head of all the studied castes, whereas <italic>SiCSP8</italic> expression was low in the male and worker head, with no statistical differences seen between the two tissues. However, <italic>SiCSP8</italic> was significantly much more highly expressed in the head of female alates as compared to the other two samples.</p>
<p>Co-expression cluster II encompassed two sub-clusters, II<sub>a</sub> (<italic>SiCSPs 10, 18</italic> and <italic>9</italic>) and II<sub>b</sub> (<italic>SiCSPs 22, 12, 3, 7, 17</italic> and <italic>15</italic>). Within sub-cluster II<sub>a</sub>, <italic>SiCSP 10</italic> and <italic>18</italic> expression levels were low to robust, respectively, and did not exhibit significant differences among castes. Conversely, SiCSP9 exhibited robust expression in the female head, which was statistically higher (<italic>P</italic> &#x003C; 0.05) than in the male head. Expression levels among the <italic>SiCSPs</italic> from sub-cluster II<sub>b</sub> varied between low and slightly high. This sub-cluster comprised general <italic>SiCSPs</italic> (i.e., <italic>SiCSPs 3</italic> and <italic>7</italic> boxed in red in <xref ref-type="fig" rid="F1">Figure 1</xref>) and <italic>S. invicta</italic>-specific <italic>SiCSPs</italic> (<italic>SiCSPs 22</italic> and <italic>12</italic> - boxed in green in <xref ref-type="fig" rid="F1">Figure 1</xref> &#x2013; and <italic>SiCSPs 17</italic> and <italic>15</italic> &#x2013; boxed in blue in <xref ref-type="fig" rid="F1">Figure 1</xref>) which, in turn, displayed similar co-expression patterns between them. <italic>SiCSP22</italic> expression was low in the head of all ant castes, and <italic>SiCSP12</italic> expression was significantly lower (<italic>P</italic> &#x003C; 0.05) in the male head than in worker and female heads. Expression levels of <italic>SiCSPs 3</italic> and <italic>7</italic> was low across adult heads, with no significant differences among castes. Similar to <italic>SiCSP12</italic>, the expression of <italic>SiCSP17</italic> and <italic>15</italic> was statistically lower (<italic>P</italic> &#x003C; 0.05) in male heads than in worker and female heads. Significant differences were observed between female and male heads in the expression of <italic>SiCSP17</italic>, but not in <italic>SiCSP15</italic>.</p>
<p>Similar to co-expression cluster II, co-expression cluster III could be divided into two sub-clusters that were phylogenetically related. Sub-cluster III<sub>a</sub> included <italic>SiCSP2</italic> and <italic>4</italic> (general CSPs &#x2013; boxed in red in <xref ref-type="fig" rid="F1">Figure 1</xref>) whose expression levels ranged from very low (<italic>SiCSP2</italic> &#x2013; worker, male and female heads) to slightly low (<italic>SiCSP4</italic> &#x2013; worker head). Significant differences in expression were seen in <italic>SiCSP4</italic> between worker (higher) and male/female heads (lower). Sub-cluster III<sub>b</sub> contained only <italic>SiCSPs</italic> found exclusively in <italic>S. invicta</italic> (i.e., <italic>SiCSPs 16</italic>, <italic>21</italic>, and <italic>11</italic> &#x2013; boxed in blue in <xref ref-type="fig" rid="F1">Figure 1</xref> and <italic>SiCSP 20</italic> &#x2013; boxed in green in <xref ref-type="fig" rid="F1">Figure 1</xref>). Within this sub-cluster, the expression of the majority of the <italic>SiCSPs</italic> (i.e., <italic>SiCSP16</italic>, <italic>21</italic>, and <italic>20</italic>) was in general not statistically different and low across worker, male and female heads, with the exception of <italic>SiCSP11</italic> which was somewhat robust expressed in worker and female heads and significantly higher (<italic>P</italic> &#x003C; 0.05) in these two tissues than in the male head.</p>
<p>One outlier, namely <italic>SiCSP1</italic>, was found within co-expression cluster IV. <italic>SiCSP1</italic> was differentially expressed among castes and was very highly expressed in the male head. It showed slightly low expression in the female head, and low expression in worker heads. Conversely, the other members of co-expression cluster IV (<italic>SiCSP6</italic>, <italic>13</italic>, and <italic>14</italic>) were expressed at somewhat low levels and did not show statistical differences across castes.</p>
</sec>
<sec id="S3.SS2.SSS4">
<title>Chemosensory Protein Expression in the Thorax of <italic>S. invicta</italic> Workers, and Male and Female Alates</title>
<p>For the most part, <italic>SiCSP</italic> expression in worker, male and female thoraces ranged between very low to robust, with the exception of <italic>SiCSPs 22</italic> and <italic>9</italic> which were highly expressed in the female thorax (<xref ref-type="fig" rid="F4">Figure 4C</xref>). Unlike for the rest of the tissues examined (i.e., antennae, head, and abdomen), the <italic>SiCSP</italic> expression profile of adult thoraces was more similar between males and females than between workers and both reproductive castes, as revealed by the hierarchical clustering analysis. Clustering by <italic>SiCSPs</italic> generated five co-expression clusters that included sub-clusters of phylogenetically related and distant CSPs sharing similar expression profiles across castes: (I) <italic>SiCSPs 2</italic>, <italic>3</italic>, and <italic>16</italic>, (II) <italic>SiCSPs 22</italic>, <italic>12</italic>, <italic>10</italic>, <italic>9</italic> and <italic>7</italic>, (III) <italic>SiCSPs 8</italic>, <italic>18</italic> and <italic>13</italic>, (IV) <italic>SiCSPs 13</italic>, <italic>19, 11</italic>, <italic>17</italic>, <italic>15</italic>, and (V) <italic>SiCSPs 6</italic>, <italic>4</italic>, <italic>21</italic>, <italic>14</italic>, <italic>1</italic>, <italic>20</italic> (<xref ref-type="fig" rid="F4">Figure 4C</xref>).</p>
<p>Expression of <italic>SiCSP2</italic> in the thorax varied between very low in workers and females to low in males. The expression of <italic>SiCSP3</italic> in the thoraces of the three castes was slightly robust and, although <italic>SiCSP3</italic> followed the trend of <italic>SiCSP2</italic>, being lower in workers and females than in male, these differences were not found to be statistically significant. <italic>SiCSP16</italic> was the outlier in co-expression cluster II, showing similar slightly robust expression across adult thoraces with no statistical differences. Co-expression cluster II consisted of two sub-clusters: II<sub>a</sub> (<italic>SiCSPs 22, 12</italic> and <italic>10</italic> &#x2013; formed by <italic>SiCSPs</italic> found exclusively in <italic>S. invicta</italic>, which are boxed in green in <xref ref-type="fig" rid="F1">Figure 1</xref>) and II<sub>b</sub> (<italic>SiCSps 7</italic> and <italic>9</italic>). Within sub-cluster II<sub>a</sub>, <italic>SiCSPs</italic> expression fluctuated between low and somewhat high. Among worker, male and female thoraces, no significance differences in expression were seen in <italic>SiCSP22</italic> and <italic>12</italic>, which exhibited slightly robust and somewhat high expression respectively in adult thoraces. <italic>SiCSP10</italic>, however, was significantly lower expressed in males than females but not in workers, showing low expression in the male and worker thoraces and slightly robust expression in the female thorax. <italic>SiCSP7</italic> and <italic>9</italic> appeared to be more expressed in the female thorax (slightly robust and high, respectively) than in the worker and male thoraces (slightly low), although only <italic>SiCSP9</italic> was found to be significantly (<italic>P</italic> &#x003C; 0.05) higher expressed in the female head.</p>
<p>The <italic>SiCSPs</italic> included in co-expression cluster III were not phylogenetically related, and their expression varied between low to robust. <italic>SiCSP8</italic>, whose expression in the thorax was low in the males, slightly robust in workers and robust in females, was statically differentially expressed across castes. <italic>SiCSP18</italic> in the thorax was low expressed in workers and somewhat robustly expressed in males and females and only showed statistical differences between workers and males and workers and females. The expression of <italic>SiCSP13</italic> in adult thoraces was slightly robust, and no significant differences were seen among castes. Co-expression cluster IV consisted of members of the two <italic>S. invicta</italic> expansions (i.e., <italic>SiCSP19</italic> &#x2013; colored in green in <xref ref-type="fig" rid="F1">Figure 1</xref> &#x2013; and <italic>SiCSPs 11</italic>, <italic>17</italic>, and <italic>15</italic> &#x2013; colored in blue in <xref ref-type="fig" rid="F1">Figure 1</xref>). <italic>SiCSPs 19</italic> and <italic>11</italic> were expressed at very low levels in the male thorax and were significantly more highly expressed in the thoraces of workers and females. <italic>SiCSP17</italic> was somewhat robustly expressed throughout adult thoraces and showed no differential expression. In contrast, <italic>SiCSP15</italic> was differentially expressed in the male thorax and exhibited low expression, as opposed to workers and females which showed moderate low expression. <italic>SiCSPs 6</italic>, <italic>4</italic>, <italic>21</italic>, <italic>14</italic>, and <italic>20</italic> formed co-expression cluster V and mainly appeared to be more highly expressed in the worker thorax than in male and female thoraces. The expression of <italic>SiCSPs 4</italic>, <italic>14</italic>, <italic>1</italic>, and <italic>20</italic> was statistically higher (<italic>P</italic> &#x003C; 0.05) in the worker thorax than in the thoraces of males and females, exhibiting moderate low (e.g., <italic>SiCSP4</italic>) and slightly robust (e.g., <italic>SiCSPs 14</italic>, <italic>1</italic> and <italic>20</italic>) expression. <italic>SiCSPs 6</italic> and <italic>21</italic> showed low to slightly robust expression, respectively, in adult thoraces with no statistical differences found among castes.</p>
</sec>
<sec id="S3.SS2.SSS5">
<title>Chemosensory Protein Expression in the Abdomen of <italic>S. invicta</italic> Workers, and Male and Female Alates</title>
<p>Overall SiCSP expression in the abdomens of <italic>S. invicta</italic> was slightly low in males and females and robust to high in workers (<xref ref-type="fig" rid="F4">Figure 4D</xref>). Despite some apparent differences in expression, the hierarchical clustering analysis exposed similar co-expression profiles between worker and female abdomens, with a more distant expression profile seen for the male abdomen. Three groupings or clusters of co-expressed <italic>CSPs</italic> were inferred by hierarchical clustering, showing phylogenetic relatedness: (I) <italic>SiCSPs 18</italic>, <italic>2</italic>, and <italic>8</italic>, (II) <italic>SiCSPs 19</italic>, <italic>6</italic>, <italic>7</italic>, <italic>14</italic>, <italic>1</italic>, <italic>22</italic>, <italic>17</italic>, <italic>16</italic>, <italic>3</italic>, <italic>15</italic>, <italic>4</italic>, <italic>13</italic>, and (III) <italic>SiCSPs 10</italic>, <italic>12</italic>, <italic>21</italic>, <italic>11</italic>, <italic>9</italic>, and <italic>20</italic> (<xref ref-type="fig" rid="F4">Figure 4D</xref>). <italic>SiCSPs 18</italic>, <italic>2</italic> and <italic>8</italic>, in co-expression cluster I, showed significantly lower expression in the female abdomen than in worker and male abdomens. <italic>SiCSP8</italic> and <italic>2</italic> expression was low in worker and male abdomens, and very low in the female abdomen, while <italic>SiCSP18</italic> was robustly expressed in workers and males and showed low expression in the female abdomen.</p>
<p>Co-expression cluster II consisted of two sub-cluster, II<sub>a</sub> (<italic>SiCSPs 19</italic>, <italic>6</italic>, <italic>7</italic>, <italic>14</italic> and <italic>1</italic>) and II<sub>b</sub> (<italic>SiCSPs 22</italic>, <italic>17</italic>, <italic>16</italic>, <italic>3</italic>, <italic>15</italic>, <italic>4</italic> and <italic>13</italic>). Within sub-cluster II<sub>a</sub>, <italic>SiCSPs 19</italic>, <italic>6</italic>, <italic>14</italic> and <italic>1</italic> exhibited significant (<italic>P</italic> &#x003C; 0.05) differential expression across all examined castes, with <italic>SiCSP7</italic> showing significant differential expression between the worker abdomen and the female abdomen. <italic>SiCSPs 19</italic>, <italic>14</italic> and <italic>1</italic> showed moderately low expression in the abdomen of males, with high (e.g., <italic>SiCSP19</italic>) to very high expression (e.g., <italic>SiCSPs 14</italic> and <italic>1</italic>) in the abdomen of workers but lower expression in the abdomen of female alates. SiCSP7 displayed moderately low expression in the female abdomen with higher expression in male and worker abdomen. Among the examined abdomen, <italic>SiCSP6</italic> expression was somewhat robust in workers, slightly low in males and low in females. In sub-cluster II<sub>b</sub>, the expression of <italic>SiCSP15</italic>, <italic>4</italic> and <italic>13</italic> in the abdomen showed significant differences across caste, while the expression of <italic>SiCSPs 22</italic>, <italic>17</italic>, <italic>16</italic> and <italic>3</italic> only presented differences between the worker and male abdomen and between the worker and female abdomen. <italic>SiCSP15</italic> and <italic>13</italic> exhibited moderately low expression in males, low expression in females but high expression in workers. Similarly, <italic>SiCSP4</italic> also showed high expression in the worker abdomen, moderately low expression in males, and low expression in females, although its expression in male and female abdomen was much lower than that seen for <italic>SiCSPs 15</italic> and <italic>13</italic>. Comparatively, <italic>SiCSPs 22</italic>, <italic>17</italic>, and <italic>16</italic> displayed quasi-identical co-expression pattern across castes, with moderately low expression in the abdomen of males and females and high expression in worker abdomen. <italic>SiCSP3</italic> also appeared more highly expressed in worker abdomen than in male and female abdomen.</p>
<p>Co-expression cluster III included only members of the two <italic>S. invicta</italic> expansions (i.e., <italic>SiCSPs 10, 12, 9</italic> and <italic>20</italic> &#x2013; boxed in green in <xref ref-type="fig" rid="F1">Figure 1</xref> &#x2013; and <italic>SiCSP21</italic> and <italic>11</italic> &#x2013; boxed in blue in <xref ref-type="fig" rid="F1">Figure 1</xref>). All the <italic>SiCSPs</italic> within this co-expression cluster showed significantly (<italic>P</italic> &#x003C; 0.05) higher levels of expression in the worker abdomen than in the male and female abdomen. <italic>SiCSP10</italic> displayed robust expression in the worker abdomen and low expression in the male and female abdomen. Expression of <italic>SiCSP12</italic> in the abdomen was much higher across the examined castes than that seen for <italic>SiCSP10</italic>, with robust levels of expression in males and females and high expression in workers. <italic>SiCSPs 21, 11</italic>, and <italic>9</italic> showed very similar expression patterns across adult abdomen samples, namely high in worker abdomen and somewhat robust in male and female abdomen. <italic>SiCSP20</italic> was also highly expressed in the worker abdomen but, unlike <italic>SiCSPs 21, 11</italic> and <italic>9, SiCSP20</italic> appeared low expressed in the male and female abdomens.</p>
</sec>
<sec id="S3.SS2.SSS6">
<title>Chemosensory Protein Expression in <italic>S. invicta</italic> Developmental Stages</title>
<p><italic>SiCSP</italic> expression was examined in <italic>S. invicta</italic> eggs, larvae, and pupae (<xref ref-type="fig" rid="F5">Figure 5</xref>). Larvae were separated into 1st/2nd (early) and 3rd/4th (late) instars based on size and morphology. One-way ANOVA yielded significant differential expression of the <italic>SiCSPs</italic> within each developmental stage examined (<italic>P</italic> &#x003C; 0.001). Overall, expression of <italic>SiCSPs</italic> was low within all the developmental stages examined. Hierarchical clustering analysis by column (i.e., developmental phase) revealed similar co-expression patterns between early- and late-instar larvae, and a much distant and reduced <italic>SiCSP</italic> expression profiles in pupae and eggs (<xref ref-type="fig" rid="F5">Figure 5</xref>).</p>
<p>In early-instar larvae, <italic>SiCSPs 14</italic>, <italic>11</italic> and <italic>20</italic> showed somewhat high expression, <italic>SiCSPs 21</italic>, <italic>17</italic>, <italic>7</italic>, <italic>1</italic>, <italic>15</italic>, <italic>13</italic>, and <italic>4</italic> robust to slightly robust expression, while the rest of the <italic>SiCSPs</italic> (<italic>SiCSPs</italic> 6, 3, 18, 10, 22, 16, 2, 8, 19, 12, 9) exhibited low to very low expression (<xref ref-type="supplementary-material" rid="FS1">Supplementary Figure S1</xref>). Although early- and late-instar larvae showed similar expression profiles, overall <italic>SiCSP</italic> expression in late-instars was much lower than in early instars, with all <italic>SiCSPs</italic> exhibiting slightly low to very low abundance of normalized transcripts. Late-instar larvae displayed low expression of <italic>SiCSPs 11</italic>, <italic>14</italic>, <italic>7</italic>, <italic>20</italic>, <italic>21</italic>, <italic>1</italic>, <italic>4</italic>, <italic>17</italic>,<italic>15</italic>, <italic>13</italic>, <italic>8</italic>, and <italic>6</italic>, and low to very low expression of the remaining <italic>SiCSPs</italic> (<italic>SiCSP19</italic>, <italic>2</italic>, <italic>12</italic>, <italic>16</italic>, <italic>18</italic>, <italic>3</italic>, <italic>22</italic>, <italic>10</italic>, <italic>9</italic>). In eggs, SiCSP expression was low to very low, while <italic>SiCSP</italic> expression levels appeared to increase for pupae, which showed significant differential expression that ranged between robust (<italic>SiCS7</italic>) to very low (<italic>SiCSP10</italic>) (<xref ref-type="supplementary-material" rid="FS1">Supplementary Figure S1</xref>).</p>
<p>One-way ANOVA analysis of the normalized expression data showed that <italic>S. invicta CSPs</italic> are significantly (<italic>P</italic> &#x003C; 0.001) differentially expressed across developmental stages (<xref ref-type="fig" rid="F5">Figure 5</xref>). Hierarchical clustering by <italic>SiCSP</italic> allowed for the identification of five co-expression CSP clusters: (I) <italic>SiCSPs 19</italic>, <italic>12</italic> and <italic>9</italic>, (II) <italic>SiCSPs 7</italic> and <italic>16</italic>, (III) <italic>SiCSPs 18</italic>, <italic>3</italic> and <italic>6</italic>, (IV) <italic>SiCSPs 22</italic>, <italic>10</italic>, <italic>20</italic>, <italic>13</italic> and <italic>17</italic>, and (V) <italic>SiCSPs 4</italic>, <italic>1</italic>, <italic>21</italic>, <italic>11</italic>, <italic>15</italic> and <italic>14</italic>. These groupings of co-expressed <italic>CSPs</italic> revealed shared expression within more closely related SiCSPs, but also cases among more distantly related SiCSPs (<xref ref-type="fig" rid="F5">Figure 5</xref>). The expression patterns of <italic>SiCSP8</italic> and <italic>2</italic> across developmental stages were not shared with the other <italic>SiCSPs</italic>. <italic>SiCSP8</italic> expression was low in early-instar larvae and pupae and moderately low in late-instar larvae and egg, with significance differences found between both larval stages, early-instar larvae and eggs, and pupae and eggs. <italic>SiCSP2</italic>, on the other hand, was more highly expressed in pupae than in the other developmental phases, which showed no significant differences among them.</p>
<p>Within co-expression cluster I, expression of <italic>SiCSPs 19, 12</italic> and <italic>9</italic> (members of one of the <italic>S. invicta</italic> expansions boxed in green in <xref ref-type="fig" rid="F1">Figure 1</xref>) ranged from low to very low among developmental stages, with expression levels of <italic>SiCSP19</italic> being statistically different across all developmental stages, and expression of <italic>SiCSP12</italic> being not statistically different between larvae (early and late instars) but distinct when comparing both larval stages (early and late instars) to pupae, and larvae (early and late instars) to eggs, and pupae to eggs. The expression of <italic>SiCSP9</italic> was very low and similar across all the developmental stages.</p>
<p>Expression of <italic>SiCSPs 7</italic> and <italic>16</italic> (co-expression cluster II) in developmental stages varied between very low (<italic>SiCSP16</italic> in late-instar larvae and eggs) to robust (<italic>SiCSP7</italic> in pupae). Significant different levels of expression were observed in <italic>SiCSP7</italic> across all the examined developmental stages, showing more robust expression in early-instar larvae and pupae than in late-instar larvae and eggs. <italic>SiCSP16</italic> displayed low expression in early-instar larvae and pupae and very low expression in late-instar larvae and eggs. While no significant differences in expression levels of <italic>SiCSP16</italic> were noticed among early and late larval instars and pupae and late larval instars and eggs, significant differences were detected between both larval instar and between early-instar larvae and eggs.</p>
<p>Within co-expression cluster III (i.e., <italic>SiCSPs 18, 3</italic> and <italic>6</italic>) expression levels ranged from moderately low to very low. <italic>SiCSPs 18</italic> and <italic>3</italic> displayed similar expression profiles, with very low expression in late-instar larvae and eggs, and low expression in early-instar larvae and pupae. Significantly different levels of expression were observed for <italic>SiCSPs 18</italic> and <italic>3</italic> between both larval stages and late-instar larvae and pupae. In addition, a significant increase of expression was observed between early-instar larvae and pupae for <italic>SiCSP18</italic> but not for <italic>SiCSP3</italic>. <italic>SiCSP6</italic> showed a comparable trend of expression to <italic>SiCSPs 18</italic> and <italic>3</italic>, although with slightly higher expression values. Also, <italic>SiCSP6</italic> displayed low levels of expression throughout all developmental stages, although expression was significantly lower in late-instar larvae and eggs than in the other two phases.</p>
<p>Co-expression cluster IV, which includes <italic>SiCSPs 22, 10, 20, 13</italic> (members of one of the <italic>S. invicta</italic> expansion boxed in green in <xref ref-type="fig" rid="F1">Figure 1</xref>) and <italic>17</italic> (member of the <italic>S. invicta</italic> expansion boxed in blue in <xref ref-type="fig" rid="F1">Figure 1</xref>), displayed somewhat high to very low expression values. <italic>SiCSP22</italic> and <italic>10</italic> showed almost identical expression profiles across developmental stages, with lower expression in eggs, pupae, and late-instar larvae than in early-instar larvae. In addition, significant differences were seen comparing eggs to pupae, and eggs to both larval stages. <italic>SiCSP20</italic> was slightly highly expressed in early-instar larvae and showed low expression in late-instar larvae, pupae and eggs. Expression of <italic>SiCSP20</italic> was significantly lower in eggs, late-instar larvae and pupae than in early-instar larvae. <italic>SiCSP13</italic> and <italic>17</italic> showed a similar trend as <italic>SiCSP20</italic>, with higher levels of expression in early-instar larvae and much (significantly) as compared to lower expression levels seen in eggs, late-instar larvae, and pupae.</p>
<p>Co-expression cluster V, comprising <italic>SiCSPs 11, 21, 15</italic> (members of the <italic>S. invicta</italic> expansion boxed in blue in <xref ref-type="fig" rid="F1">Figure 1</xref>), <italic>1, 4</italic> and <italic>14</italic>, exhibited similarities with cluster IV since the overall expression of this group of <italic>SiCSPs</italic> was robust in early-instar larvae but low/very low in the other examined stages. <italic>SiCSP4, 1</italic>, and <italic>21</italic> were statistically differentially expressed (<italic>P</italic> &#x003C; 0.05) across all developmental phases, whereas <italic>SiCSP11, 15</italic> and <italic>14</italic> expression levels were statistically different (<italic>P</italic> &#x003C; 0.05) comparing both larval stages (early to late instars), larvae to eggs, and larvae to pupae, and not significantly different when comparing eggs to pupae.</p>
</sec>
</sec>
<sec id="S3.SS3">
<title>Protein Motifs Identified in the <italic>S. invicta</italic> CSPs</title>
<p>The MEME server (<xref ref-type="bibr" rid="B4">Bailey and Elkan, 1994</xref>) was used to identify conserved amino acid motifs present in the 21 <italic>S. invicta</italic> CSPs, along with the other homologs used to build the phylogenetic tree. Eight distinct motifs were identified with e-values &#x003C; 1.2 e<sup>&#x2013;37</sup>, and their distribution within <italic>S. invicta</italic> (<xref ref-type="fig" rid="F6">Figure 6</xref>) and the 92 CSP set (<xref ref-type="supplementary-material" rid="FS2">Supplementary Figure S2</xref>) determined. Motifs 1 and 2 (red and turquoise, respectively) were found in all of the CSPs examined, including the full 21 SiCSPs, and contain the four conserved cysteines. Motif 3 (neon green), at the C-terminus, was also found in the majority of SiCSPs, with the exception of SiCSPs 10, 12, 13, 19, 20, 22, and 9 (<italic>S. invicta</italic> CSP expansion boxed in green in <xref ref-type="fig" rid="F1">Figure 1</xref>) which contained a separate motif (8, in pink) instead, that appeared to be unique to <italic>S. invicta</italic> and mutually exclusive with motif 3 (<xref ref-type="supplementary-material" rid="FS2">Supplementary Figure S2</xref>). Motif 4 (purple) was also widely distributed within the majority of the CSPs examined. In <italic>S. invicta</italic>, motif 4 was missing in SiCSPs 12, 22, and 2, although in the latter, this motif appeared to be substituted by a shuffled motif 6 (green) (<xref ref-type="fig" rid="F6">Figure 6</xref>). Motif 6 (green) was found as an extended C-terminus in SiCSPs 3, 4, and 7 but not in SiCSP2. Motif 5 and 7 (orange and dark blue, respectively) were found at the N-terminus of the CSPs and are both part of the signal peptide sequence. Motif 5 was found in most of the CSPs, with only few exceptions (e.g., SiCSP8), while Motif 7 only appeared in the CSPs that formed part of the CSP ant expansion (<xref ref-type="supplementary-material" rid="FS2">Supplementary Figure S2</xref>, CSP ant expansion illustrated within the gray shaded area in <xref ref-type="fig" rid="F1">Figure 1</xref>).</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption><p>Motif analyses of <italic>S. invicta</italic> CSPs. Position of the 7 most common amino acid motifs identified in <italic>S. invicta</italic> and 71 more CSPs using the MEME server (version 4.11.2, <ext-link ext-link-type="uri" xlink:href="http://meme-suite.org/tools/meme">http://meme-suite.org/tools/meme</ext-link>).</p></caption>
<graphic xlink:href="fphys-11-585883-g006.tif"/>
</fig>
</sec>
</sec>
<sec id="S4">
<title>Discussion</title>
<p>Chemosensory proteins represent one branch of at least three phylogenetically un-related insect protein families, that are however, considered partially functionally related. The other two members are the arthropod OBPs and the lipid transporters classified as Niemann-Pick type C<sub>2</sub> proteins (NPC<sub>2</sub>) (<xref ref-type="bibr" rid="B42">Pelosi et al., 2014</xref>). All three classes of proteins are small (&#x223C;12&#x2013;20 kDa), soluble, ligand-binding proteins, that are often considered capable of solubilizing hydrophobic compounds including odorants, volatiles, and lipids within the aqueous milieu of the organism and/or cell. These molecules can encompass pheromones, semio-chemicals, hydrocarbons, fatty acids, amongst other hydrophobic compounds, but can also potentially include soluble substrates such as carbohydrates and other molecules. In arthropods, CSPs exist as multigene families, and have been associated with chemoreception of both olfactory and gustatory cues. Aside from numerous studies showing expression of CSPs in antennae (for recent reports see <xref ref-type="bibr" rid="B2">Antony et al., 2016</xref>; <xref ref-type="bibr" rid="B21">Hu et al., 2016</xref>; <xref ref-type="bibr" rid="B56">Xue et al., 2016</xref>; <xref ref-type="bibr" rid="B25">Jing et al., 2020</xref>), these proteins have also been found in many other tissues some of which contain specialized sensilla, e.g., in tarsal chemosensilla that respond to non-volatile surface chemicals and are involved in gustatory detection (<xref ref-type="bibr" rid="B35">Ma et al., 2016</xref>). In the latter case, transcriptomic analyses in the tea geometrid, <italic>Ectropis obliqua</italic>, revealed sex and leg-specific expression of CSPs, suggesting their involvement in taste detection. From these and other studies, it is now well established that many CSPs are expressed outside olfactory organs, functioning in chemosensory (e.g., taste) and non-chemosensory roles (e.g., development). Transcriptomic analyses in several ant species revealed that many OBPs and CSPs are expressed primarily in non-olfactory tissues (<xref ref-type="bibr" rid="B37">McKenzie et al., 2014</xref>), strongly suggesting non-olfactory roles for many of these proteins in ants. In <italic>A. mellifera</italic>, <italic>AmCSP5</italic> is not expressed in the antennae, but instead has been implicated in embryonic development, with high expression noted in queen ovaries and eggs (<xref ref-type="bibr" rid="B36">Maleszka et al., 2007</xref>). In addition, proteomic analyses have indicated the presence of CSPs in insect hemolymph (<xref ref-type="bibr" rid="B20">Hou et al., 2016</xref>).</p>
<p>Our data indicate caste-specific expression patterns of the <italic>S. invicta</italic> CSPs. Transcripts for most of the <italic>SiCSPs</italic> could be detected at robust and high levels in worker and male antennae. Levels of <italic>SiCSP</italic> expression in female alate antennae were notably lower, with the expression of at least four <italic>SiCSPs</italic>, i.e., <italic>22</italic>, <italic>2</italic>, <italic>6</italic>, and <italic>14</italic>, being low in female alate antennae but robustly expressed in workers and male alate antennae. This suggests that in <italic>S. invicta</italic>, similar to <italic>H. saltator</italic>, olfactory responses depend on reproductive status (<xref ref-type="bibr" rid="B11">Ghaninia et al., 2017</xref>). Pheromonal inhibition by the reproductive queen(s) and/or ecological adaptations - since female alates remain most of the time in the nest and are not involved in colony maintenance &#x2013; might account for <italic>S. invicta</italic> female alates (i.e., virgin females) being less responsive to olfactory cues than male alates and workers. Generalized robust expression of <italic>CSPs</italic> in male alates, on the other hand, might reflect the need for male alates to respond to olfactory signals in the context of mating (<xref ref-type="bibr" rid="B10">Ghaninia et al., 2018</xref>).</p>
<p>In workers, <italic>SiCSP19</italic> was the most highly expressed <italic>CSP</italic>, consistent with its original discovery in <italic>S. invicta</italic> antennae (<xref ref-type="bibr" rid="B29">Leal and Ishida, 2008</xref>; <xref ref-type="bibr" rid="B12">Gonzalez et al., 2009</xref>). The protein product of <italic>SiCSP19</italic> has been localized to the tip (A9-A10 segments) of the antennal club (<xref ref-type="bibr" rid="B29">Leal and Ishida, 2008</xref>). As the porous sensilla contain self-enclosed neuronal compartments located in the <italic>S. invicta</italic> antennal club that are not contiguous with the cuticle and hemolymph of the rest of the antennae (<xref ref-type="bibr" rid="B44">Renthal et al., 2003</xref>), it was concluded that SiCSP19 is most likely found in the olfactory sensilla (<xref ref-type="bibr" rid="B12">Gonzalez et al., 2009</xref>). Intriguingly, SiCSP19 is not an ortholog of the <italic>C. japonicus</italic> CSP1 (CjCSP1) that has been shown to be localized to chemosensory sensilla and implicated as critical for binding of hydrocarbons implicated in nestmate/non-nestmate discrimination (<xref ref-type="bibr" rid="B41">Ozaki et al., 2005</xref>). The ortholog of <italic>CjCSP1</italic> appears to be <italic>SiCSP1</italic>, which, while equally highly expressed across female alate, male and worker antennae, is more preferentially expressed in the male head. Proteomic analyses of <italic>S. invicta</italic> worker antennae also revealed the presence of SiCSP12 (<xref ref-type="bibr" rid="B12">Gonzalez et al., 2009</xref>), that was not only highly expressed in worker antennae in our dataset, but also in male and female alate antennae. In addition, robust to high expression of <italic>SiCSPs 8</italic> and <italic>11</italic> were noted in worker, male and female alate antennae. Whether these proteins (and/or other SiCSPs) are made in the antennae will require further deeper probing of the proteomic content of these structures.</p>
<p>A limited number of downstream behavioral effects have been attributed to CSPs involved in odorant detection. The <italic>B. dorsalis BdorCSP2</italic> has been implicated in insect responses to rhodojaponin-III, a non-volatile diterpene which shows antifeedant and anti-ovipositing activities against insects (<xref ref-type="bibr" rid="B57">Yi et al., 2013</xref>). BdorCSP2 was shown to be capable of binding rhodojaponin-III <italic>in vitro</italic>, and RNAi mediated knockdown of <italic>BdorCSP2</italic> decreased the effects of the chemical compound on the insect. <italic>SiCSP4</italic> is the ortholog of <italic>BdorCSP2</italic> and is robustly expressed in the antennae of fire ant males and workers and in much higher levels in the worker abdomen. Several CSPs, highly expressed in <italic>A. cerana</italic> antennae, have been shown to bind various floral chemicals (<xref ref-type="bibr" rid="B30">Li H.L. et al., 2016</xref>) and two CSPs in the mite <italic>Tyrophagus putrescentiae</italic> have been implicated in host recognition (<xref ref-type="bibr" rid="B31">Li X.M. et al., 2016</xref>), although these studies lack experimental designs that would directly test the role(s) of these CSPs in mediating specific behaviors. Differential male and female-specific expression of olfaction genes, has been reported in a number of insects including the oriental fruit fly (<italic>Bactrocera dorsalis</italic>) and the hover fly <italic>Scaeva pyrastri</italic> (<xref ref-type="bibr" rid="B31">Li X.M. et al., 2016</xref>; <xref ref-type="bibr" rid="B32">Liu et al., 2016</xref>). Caste and sex-specific expression of chemosensory genes has also been reported in termites and honeybees (<xref ref-type="bibr" rid="B39">Mitaka et al., 2016</xref>; <xref ref-type="bibr" rid="B24">Jain and Brockmann, 2020</xref>). When the absolute number of <italic>SiCSP</italic> transcripts is examined by tissue distribution, clear caste- and sex-specific differences are apparent. In workers, <italic>SiCSP</italic> transcripts are distributed mainly in the antennae and the abdomen. This suggests these SiCSPs may function as ligand carriers in antennal chemosensation, and that in the abdomen these proteins may act to bind/sequester endogenous hydrophobic compounds that can include pheromones (for release), hydrocarbons and lipids (for delivery to the cuticle), and/or endocrine molecules (endogenous hormones and other compounds). In contrast, in female alates, <italic>SiCSP</italic> transcripts appear to be widely distributed in the thorax, and to a lesser extent, the male thorax appeared also enriched in <italic>SiCSP</italic> expression. These tissues include the legs and the wings that contain tarsal and wing chemosensory sensilla, some of which have gustatory roles (<xref ref-type="bibr" rid="B9">Fleischer and Krieger, 2018</xref>; <xref ref-type="bibr" rid="B17">He et al., 2019</xref>). As reveal by the clustering analyses, similar patterns of CSP expression in <italic>S. invicta</italic> male and female thoraces could indicate roles for these CSPs in gustatory and/or taste functions in this ant species. Gustatory pathways in <italic>Drosophila</italic> have been reported to be activated by pheromonal cuticular hydrocarbons (CHCs) (<xref ref-type="bibr" rid="B9">Fleischer and Krieger, 2018</xref>). In ants, select CHCs also mediate social interaction between colony members and can act as mating pheromones (<xref ref-type="bibr" rid="B26">Kather and Martin, 2015</xref>; <xref ref-type="bibr" rid="B47">Slone et al., 2017</xref>). In <italic>Drosophila</italic>, some pheromonal CHCs are primarily detected via taste organs on legs and wings (<xref ref-type="bibr" rid="B9">Fleischer and Krieger, 2018</xref>). It is intriguing to speculate that in male and female <italic>S. invicta</italic> alates, some <italic>SiCSPs</italic>, (particularly, <italic>S. invicta</italic>-specific CSPs) expressed in the thorax may function to bind mating CHCs, and that sensilla located in the male and female thorax participate in mate recognition.</p>
<p><italic>SiCSP</italic> expression levels and tissue distribution appear to also reflect the task distribution of the various castes. Workers are &#x201C;pluripotent&#x201D; in task and are hence exposed to a significantly larger range of environmental compounds than female or male alates. Female alates need to mate, establish the colony, and subsequently act as the reproductive unit of the colony, and thus antennal chemosensation is likely to be far less important than in workers. In addition, queen pheromones act to regulate offspring behavior, but the queen herself may be largely non-responsive to these chemicals. Finally, male alates have only one task, and that is to mate. Therefore comparatively, in males, the requirements for and functions of the CSPs may be more limited, consistent with their lower expression and scattered tissue distribution.</p>
<p>A lack of congruent patterns expression of orthologous CSPs was generally seen when our data are compared to CSP expression analyses in <italic>C. japonicus</italic>. For example, the most highly expressed CSPs in <italic>C. japonicus</italic> were <italic>CjCSPs 1</italic>, <italic>3</italic>, <italic>7</italic>, <italic>12</italic>, and <italic>13</italic>, with only <italic>CjCSPs1</italic>, <italic>3</italic>, and <italic>7</italic> having true orthologs in <italic>S. invicta</italic> (<xref ref-type="bibr" rid="B18">Hojo et al., 2015</xref>). Of the corresponding <italic>SiCSPs</italic> only <italic>SiCSP1</italic> was highly expressed in worker antennae. These data suggest little cross-species conservation of the expression patterns of orthologous CSPs.</p>
<p>In general, the <italic>SiCSPs</italic> were poorly expressed in the developmental stages examined and were especially low in eggs and pupae. Some expression was seen in the various larval stages, and when compared across the developmental stages, <italic>SiCSP</italic> expression appeared to peak in 1st/2nd instar larvae. These data suggest the potential for a role for some <italic>SiCSPs</italic> as ligand carriers in larvae, potentially for distribution and delivery of cuticular hydrocarbons and lipids and/or of developmental hormones and other molecules. Apparent co-expression of some <italic>CSPs</italic> has been reported in the antennae of <italic>C. japonicus</italic> (<xref ref-type="bibr" rid="B18">Hojo et al., 2015</xref>). Hierarchical clustering analyses examining potential co-regulation of phylogenetically related <italic>SiCSPs</italic> revealed patterns of co-expression of <italic>SiCSPs</italic> expression within caste tissues and between developmental stages. However, notable divergences were also seen. No robust co-expression consistent with phylogeny was observed among <italic>SiCSPs</italic> across worker, male and female alate antennae. In contrast, in the head, thorax, abdomen and developmental stages co-expression of phylogenetic sub-groups could be discerned. Our data revealed potential examples of co-regulation of both divergent and convergent lines of <italic>SiCSPs</italic>.</p>
<p>These data, together with the characterization of the expression profiles of the <italic>S. invicta</italic> OBPs (<xref ref-type="bibr" rid="B58">Zhang et al., 2016</xref>), indicate that both OBPs and CSPs show discrete expression profiles within <italic>S. invicta</italic>. These profiles extended between workers, male and female alates, and across various developmental stages and supports an emerging model of these proteins as general ligand carriers involved in a wide range of physiological processes beyond olfaction and/or chemosensation. Comparison between <italic>SiOBP</italic> expression distribution and <italic>SiCSP</italic> expression revealed a number of interesting insights. Similar to the <italic>SiCSPs</italic>, only a subset of the <italic>SiOBPs</italic> were highly expressed in the worker, male and female alate antennae. However, <italic>SiCSP</italic> expression was more enriched in the antennae than <italic>SiOBP</italic> expression (<xref ref-type="bibr" rid="B58">Zhang et al., 2016</xref>). In the abdomen, <italic>SiCSP</italic> expression was enhanced in the worker but not in male and female alates, while <italic>SiOBP</italic> expression showed the opposite trend (more highly expressed in the male and female abdomen than in the worker abdomen). The fire ant head and thorax tissues showed lower expression of the <italic>SiCSPs</italic> than of the <italic>SiOBPs</italic>. In contrast, <italic>SiCSP</italic> expression was much higher than <italic>SiOBP</italic> expression in the various developmental stages (<xref ref-type="bibr" rid="B58">Zhang et al., 2016</xref>). These comparisons indicate a general inverse correlation between <italic>SiCSP</italic> and <italic>SiOBP</italic> expression. This may suggest compensatory functions and/or that these proteins might perform unique tasks within different tissues and, in some instance, may complement each other. Thus, the differential regulation of <italic>SiCSPs</italic> and <italic>SiOBPs</italic> may be linked to the substrates they may recognize.</p>
<p>Motif analyses of the SiCSPs revealed a core set of three motifs (i.e., motifs 4, 1, and 2) found in almost all of the SiCSPs as well as two slightly different conserved signal peptide sequences consisting on motif 5 alone and a combination of motifs 7 and 5. Additional motifs (i.e., 3, 6, and 8) were found at the C-terminal end in subsets of the SiCSPs. Essentially, three different C-terminal ends were observed in SiCSPs, which yielded three distinct groups of SiCSPs. Group I, which included one of the fire ant-specific CSP gene expansion (SiCSPs 9, 22, 10, 19, 20,13, and 12) and the shortest C-terminus, containing motif 8. Group II containing motif 3 and including SiCSPs 1, 14, 18 and the second fire ant-specific CSP gene expansion (i.e., SiCSPs 21, 15, 11, 17, and 16). And Group III (SiCSPs 8, 6, 2, 3, 4, and 7) with the longest C-terminal end, which comprised motifs 3 and 6. Within Group III, SiCSP8 and 6 exhibited a degenerate motif 6 and SiCSP2 a much shorter C-terminus as motif 6 appeared shuffled toward the N-terminus, where it replaced motif 4. The similarities and differences observed amongst SiCSPs within the core motifs, the N-terminal and C-terminal ends suggest that SiCSP functional diversification was mainly driven by gradual changes occurring at the periphery of the protein core via point mutations, indels and domain shuffling (<xref ref-type="bibr" rid="B3">Bagowski et al., 2010</xref>). Nevertheless, the consequences of divergences and convergences of SiCSP expression, as well as the motif analyses that likely reflects structural and functional properties, will only be understood as more information concerning their ligand binding affinities and downstream interacting partners are uncovered.</p>
<p>Our data provide a foundation for future exploration of the functional roles of insect CSPs and suggest that like OBPs, CSPs participate in diverse biological processes as ligand carriers. The systematic expression profiles of both families of SiCSPs and SiOBPs are now available and can be used to identify targets that may participate in specific physiological processes. In addition, these data can be used as a basis for exploring factors that may contribute to ant colony organization and development, self/non-self-recognition and communication, and/or identify potential targets for ant control, i.e., attractants and/or repellents. Combining these expression data with the ligand binding profiles of these proteins is a much-needed future step in their functional characterization (e.g., behavioral studies combined with RNAi). Such ligands may be derived from the environment during chemosensation or endogenously within the organism. In addition, the differential expression of the <italic>CSPs</italic> and <italic>OBPs</italic> in ant castes may be one defining feature or signature that can account for the different behavioral responsiveness of each caste in terms of reacting to chemical cues that lead to task allocation, mating, and other features of the social nature of this ant species.</p>
</sec>
<sec id="S5">
<title>Data Availability Statement</title>
<p>All the data supporting the findings is contained within the article and/or <xref ref-type="supplementary-material" rid="FS1">Supplementary Materials</xref>.</p>
</sec>
<sec id="S6">
<title>Author Contributions</title>
<p>NK, AW, WZ, AO-U, and JB initiated the project, conceived and designed the study. WZ, AW, AO-U, and JB performed the samples collection, library constructions, q-RT-PCR, data processing, bioinformatic analyses, data processing, and contributed to revisions of the manuscript. JB and YX provided resources and/or materials for the study. NK, AO-U, and YX wrote the manuscript, contributed to bioinformatic analyses, data processing, and data interpretation. All authors contributed to the article and approved the submitted version.</p>
</sec>
<sec id="conf1">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
</body>
<back>
<fn-group>
<fn fn-type="financial-disclosure">
<p><bold>Funding.</bold> This work was supported by funds from the University of Florida Institute of Food and Agricultural Sciences, USDA-NIFA award 2019-05150 and NSF grant IOS-1557704 to NK, and from funds from the National Natural Science Foundation of China (award 32001961) and the Natural Science Foundation of Chongqin (award cstc2019jcyj-bsh0070) to WZ.</p>
</fn>
</fn-group>
<ack>
<p>Photographs of worker ants, male, and female alates were taken by AO-U and NK with kind assistance from Tyler Jones, UF-IFAS Communications Photography Department. Ant anatomy images are open access and were originally made by Robert Svensson (<ext-link ext-link-type="uri" xlink:href="http://www.msitua.net">msitua.net</ext-link>).</p>
</ack>
<sec id="S9" sec-type="supplementary material"><title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fphys.2020.585883/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fphys.2020.585883/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Image_1.TIFF" id="FS1" mimetype="image/tiff" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Supplementary Figure 1</label>
<caption><p>Homogenous subsets yielded by the comparison of means via HSD Tukey&#x2019;s test within developmental stage. Within each developmental stage, different letters indicate significant differences at <italic>P</italic> &#x003C; 0.05. <italic>SiCSPs</italic> for each developmental stage are arranged from low to high levels of expression.</p></caption>
</supplementary-material>
<supplementary-material xlink:href="Image_2.TIFF" id="FS2" mimetype="image/tiff" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Supplementary Figure 2</label>
<caption><p>Motif analyses of CSPs. <bold>(A)</bold> Positions of indicated motifs in the amino acid sequences of 92 CSPs. <bold>(B)</bold> Most common amino acid motifs identified in 92 CSPs using the MEME server (version 4.11.2, <ext-link ext-link-type="uri" xlink:href="http://meme-suite.org/tools/meme">http://meme-suite.org/tools/meme</ext-link>). The size (height) of the amino acid represents the degree of conservation of that amino acid in the consensus sequence.</p></caption>
</supplementary-material>
<supplementary-material xlink:href="Table_1.DOCX" id="TS1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Supplementary Table 1</label>
<caption><p>List of SiCSP genes.</p></caption>
</supplementary-material>
<supplementary-material xlink:href="Table_1.DOCX" id="TS2" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Supplementary Table 2</label>
<caption><p>List of primers.</p></caption>
</supplementary-material>
<supplementary-material xlink:href="Table_1.DOCX" id="TS3" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Supplementary Table 3</label>
<caption><p>Validation of primer pairs.</p></caption>
</supplementary-material>
<supplementary-material xlink:href="Table_1.DOCX" id="TS4" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Supplementary Table 4</label>
<caption><p>Proteins used in phylogenetic analyses.</p></caption>
</supplementary-material>
<supplementary-material xlink:href="Table_1.DOCX" id="SD1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Supplementary Datasheet 1</label>
<caption><p>Ct values yielded by the RT-qPCR analyses, and absolute and normalised number of transcripts for the developmental stages, workers, males and female alates.</p></caption>
</supplementary-material>
</sec>
<ref-list>
<title>References</title>
<ref id="B1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Angeli</surname> <given-names>S.</given-names></name> <name><surname>Ceron</surname> <given-names>F.</given-names></name> <name><surname>Scaloni</surname> <given-names>A.</given-names></name> <name><surname>Monti</surname> <given-names>M.</given-names></name> <name><surname>Monteforti</surname> <given-names>G.</given-names></name> <name><surname>Minnocci</surname> <given-names>A.</given-names></name><etal/></person-group> (<year>1999</year>). <article-title>Purification, structural characterization, cloning and immunocytochemical localization of chemoreception proteins from <italic>Schistocerca gregaria</italic>.</article-title> <source><italic>Eur. J. Biochem.</italic></source> <volume>262</volume> <fpage>745</fpage>&#x2013;<lpage>754</lpage>. <pub-id pub-id-type="doi">10.1046/j.1432-1327.1999.00438.x</pub-id> <pub-id pub-id-type="pmid">10411636</pub-id></citation></ref>
<ref id="B2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Antony</surname> <given-names>B.</given-names></name> <name><surname>Soffan</surname> <given-names>A.</given-names></name> <name><surname>Jakse</surname> <given-names>J.</given-names></name> <name><surname>Abdelazim</surname> <given-names>M. M.</given-names></name> <name><surname>Aldosari</surname> <given-names>S. A.</given-names></name> <name><surname>Aldawood</surname> <given-names>A. S.</given-names></name><etal/></person-group> (<year>2016</year>). <article-title>Identification of the genes involved in odorant reception and detection in the palm weevil <italic>Rhynchophorus ferrugineus</italic>, an important quarantine pest, by antennal transcriptome analysis.</article-title> <source><italic>BMC Genomics</italic></source> <volume>17</volume>:<issue>69</issue>. <pub-id pub-id-type="doi">10.1186/s12864-016-2362-6</pub-id> <pub-id pub-id-type="pmid">26800671</pub-id></citation></ref>
<ref id="B3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bagowski</surname> <given-names>C. P.</given-names></name> <name><surname>Bruins</surname> <given-names>W.</given-names></name> <name><surname>Te Velthuis</surname> <given-names>A. J. W.</given-names></name></person-group> (<year>2010</year>). <article-title>The nature of protein domain evolution: shaping the interaction network.</article-title> <source><italic>Curr. Genom.</italic></source> <volume>11</volume> <fpage>368</fpage>&#x2013;<lpage>376</lpage>. <pub-id pub-id-type="doi">10.2174/138920210791616725</pub-id> <pub-id pub-id-type="pmid">21286315</pub-id></citation></ref>
<ref id="B4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bailey</surname> <given-names>T. L.</given-names></name> <name><surname>Elkan</surname> <given-names>C.</given-names></name></person-group> (<year>1994</year>). <article-title>Fitting a mixture model by expectation maximization to discover motifs in biopolymers.</article-title> <source><italic>Proc. Int. Conf. Intell. Syst. Mol. Biol.</italic></source> <volume>2</volume> <fpage>28</fpage>&#x2013;<lpage>36</lpage>.</citation></ref>
<ref id="B5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ban</surname> <given-names>L.</given-names></name> <name><surname>Scaloni</surname> <given-names>A.</given-names></name> <name><surname>Brandazza</surname> <given-names>A.</given-names></name> <name><surname>Angeli</surname> <given-names>S.</given-names></name> <name><surname>Zhang</surname> <given-names>L.</given-names></name> <name><surname>Yan</surname> <given-names>Y.</given-names></name><etal/></person-group> (<year>2003</year>). <article-title>Chemosensory proteins of <italic>Locusta migratoria</italic>.</article-title> <source><italic>Insect. Mol. Biol.</italic></source> <volume>12</volume> <fpage>125</fpage>&#x2013;<lpage>134</lpage>. <pub-id pub-id-type="doi">10.1046/j.1365-2583.2003.00394.x</pub-id> <pub-id pub-id-type="pmid">12653934</pub-id></citation></ref>
<ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Briand</surname> <given-names>L.</given-names></name> <name><surname>Swasdipan</surname> <given-names>N.</given-names></name> <name><surname>Nespoulous</surname> <given-names>C.</given-names></name> <name><surname>Bezirard</surname> <given-names>V.</given-names></name> <name><surname>Blon</surname> <given-names>F.</given-names></name> <name><surname>Huet</surname> <given-names>J. C.</given-names></name><etal/></person-group> (<year>2002</year>). <article-title>Characterization of a chemosensory protein (ASP3c) from honeybee (<italic>Apis mellifera</italic> L.) as a brood pheromone carrier.</article-title> <source><italic>Eur. J. Biochem.</italic></source> <volume>269</volume> <fpage>4586</fpage>&#x2013;<lpage>4596</lpage>.</citation></ref>
<ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dani</surname> <given-names>F. R.</given-names></name> <name><surname>Michelucci</surname> <given-names>E.</given-names></name> <name><surname>Francese</surname> <given-names>S.</given-names></name> <name><surname>Mastrobuoni</surname> <given-names>G.</given-names></name> <name><surname>Cappellozza</surname> <given-names>S.</given-names></name> <name><surname>La Marca</surname> <given-names>G.</given-names></name><etal/></person-group> (<year>2011</year>). <article-title>Odorant-binding proteins and chemosensory proteins in pheromone detection and release in the silkmoth <italic>Bombyx mori</italic>.</article-title> <source><italic>Chem. Sens.</italic></source> <volume>36</volume> <fpage>335</fpage>&#x2013;<lpage>344</lpage>. <pub-id pub-id-type="doi">10.1093/chemse/bjq137</pub-id> <pub-id pub-id-type="pmid">21220518</pub-id></citation></ref>
<ref id="B8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fan</surname> <given-names>Y.</given-names></name> <name><surname>Pereira</surname> <given-names>R. M.</given-names></name> <name><surname>Kilic</surname> <given-names>E.</given-names></name> <name><surname>Casella</surname> <given-names>G.</given-names></name> <name><surname>Keyhani</surname> <given-names>N. O.</given-names></name></person-group> (<year>2012</year>). <article-title>Pyrokinin beta-neuropeptide affects necrophoretic behavior in fire ants (S. invicta), and expression of beta-NP in a mycoinsecticide increases its virulence.</article-title> <source><italic>PLoS One</italic></source> <volume>7</volume>:<issue>e26924</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0026924</pub-id> <pub-id pub-id-type="pmid">22238569</pub-id></citation></ref>
<ref id="B9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fleischer</surname> <given-names>J.</given-names></name> <name><surname>Krieger</surname> <given-names>J.</given-names></name></person-group> (<year>2018</year>). <article-title>Insect pheromone receptors - key elements in sensing intraspecific chemical signals.</article-title> <source><italic>Front. Cell. Neurosci.</italic></source> <volume>12</volume>:<issue>425</issue>. <pub-id pub-id-type="doi">10.3389/fphys.2014.00425</pub-id> <pub-id pub-id-type="pmid">25400590</pub-id></citation></ref>
<ref id="B10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ghaninia</surname> <given-names>M.</given-names></name> <name><surname>Berger</surname> <given-names>S. L.</given-names></name> <name><surname>Reinberg</surname> <given-names>D.</given-names></name> <name><surname>Zwiebel</surname> <given-names>L. J.</given-names></name> <name><surname>Ray</surname> <given-names>A.</given-names></name> <name><surname>Liebig</surname> <given-names>J.</given-names></name></person-group> (<year>2018</year>). <article-title>Antennal olfactory physiology and behavior of males of the ponerine ant harpegnathos saltator.</article-title> <source><italic>J. Chem. Ecol.</italic></source> <volume>44</volume> <fpage>999</fpage>&#x2013;<lpage>1007</lpage>. <pub-id pub-id-type="doi">10.1007/s10886-018-1013-6</pub-id> <pub-id pub-id-type="pmid">30191433</pub-id></citation></ref>
<ref id="B11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ghaninia</surname> <given-names>M.</given-names></name> <name><surname>Haight</surname> <given-names>K.</given-names></name> <name><surname>Berger</surname> <given-names>S. L.</given-names></name> <name><surname>Reinberg</surname> <given-names>D.</given-names></name> <name><surname>Zwiebel</surname> <given-names>L. J.</given-names></name> <name><surname>Ray</surname> <given-names>A.</given-names></name><etal/></person-group> (<year>2017</year>). <article-title>Chemosensory sensitivity reflects reproductive status in the ant <italic>Harpegnathos saltator</italic>.</article-title> <source><italic>Sci. Rep.</italic></source> <volume>7</volume>:<issue>3732</issue>.</citation></ref>
<ref id="B12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gonzalez</surname> <given-names>D.</given-names></name> <name><surname>Zhao</surname> <given-names>Q.</given-names></name> <name><surname>McMahan</surname> <given-names>C.</given-names></name> <name><surname>Velasquez</surname> <given-names>D.</given-names></name> <name><surname>Haskins</surname> <given-names>W. E.</given-names></name> <name><surname>Sponsel</surname> <given-names>V.</given-names></name><etal/></person-group> (<year>2009</year>). <article-title>The major antennal chemosensory protein of red imported fire ant workers.</article-title> <source><italic>Insect. Mol. Biol.</italic></source> <volume>18</volume> <fpage>395</fpage>&#x2013;<lpage>404</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-2583.2009.00883.x</pub-id> <pub-id pub-id-type="pmid">19523071</pub-id></citation></ref>
<ref id="B13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gotzek</surname> <given-names>D.</given-names></name> <name><surname>Ross</surname> <given-names>K. G.</given-names></name></person-group> (<year>2008</year>). <article-title>Experimental conversion of colony social organization in fire ants (<italic>Solenopsis invicta</italic>): worker genotype manipulation in the absence of queen effects.</article-title> <source><italic>J. Insect. Behav.</italic></source> <volume>21</volume> <fpage>337</fpage>&#x2013;<lpage>350</lpage>. <pub-id pub-id-type="doi">10.1007/s10905-008-9130-7</pub-id></citation></ref>
<ref id="B14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gu</surname> <given-names>S. H.</given-names></name> <name><surname>Wang</surname> <given-names>S. Y.</given-names></name> <name><surname>Zhang</surname> <given-names>X. Y.</given-names></name> <name><surname>Ji</surname> <given-names>P.</given-names></name> <name><surname>Liu</surname> <given-names>J. T.</given-names></name> <name><surname>Wang</surname> <given-names>G. R.</given-names></name><etal/></person-group> (<year>2012</year>). <article-title>Functional characterizations of chemosensory proteins of the alfalfa plant bug <italic>Adelphocoris lineolatus</italic> indicate their involvement in host recognition.</article-title> <source><italic>PLoS One</italic></source> <volume>7</volume>:<issue>e42871</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0042871</pub-id> <pub-id pub-id-type="pmid">22900060</pub-id></citation></ref>
<ref id="B15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Guntur</surname> <given-names>K. V. P.</given-names></name> <name><surname>Velasquez</surname> <given-names>D.</given-names></name> <name><surname>Chadwell</surname> <given-names>L.</given-names></name> <name><surname>Carroll</surname> <given-names>C.</given-names></name> <name><surname>Weintraub</surname> <given-names>S.</given-names></name> <name><surname>Cassill</surname> <given-names>J. A.</given-names></name><etal/></person-group> (<year>2004</year>). <article-title>Apolipophorin-III-like protein expressed in the antenna of the red imported fire ant, <italic>Solenopsis invicta</italic> Buren (Hymenoptera : Formicidae).</article-title> <source><italic>Arch. Insect. Biochem. Physiol.</italic></source> <volume>57</volume> <fpage>101</fpage>&#x2013;<lpage>110</lpage>. <pub-id pub-id-type="doi">10.1002/arch.20019</pub-id> <pub-id pub-id-type="pmid">15484261</pub-id></citation></ref>
<ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Guo</surname> <given-names>W.</given-names></name> <name><surname>Wang</surname> <given-names>X.</given-names></name> <name><surname>Ma</surname> <given-names>Z.</given-names></name> <name><surname>Xue</surname> <given-names>L.</given-names></name> <name><surname>Han</surname> <given-names>J.</given-names></name> <name><surname>Yu</surname> <given-names>D.</given-names></name><etal/></person-group> (<year>2011</year>). <article-title><italic>CSP</italic> and <italic>takeout</italic> genes modulate the switch between attraction and repulsion during behavioral phase change in the migratory locust.</article-title> <source><italic>PLoS Genet.</italic></source> <volume>7</volume>:<issue>e1001291</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.1001291</pub-id></citation></ref>
<ref id="B17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>He</surname> <given-names>Z.</given-names></name> <name><surname>Luo</surname> <given-names>Y.</given-names></name> <name><surname>Shang</surname> <given-names>X.</given-names></name> <name><surname>Sun</surname> <given-names>J. S.</given-names></name> <name><surname>Carlson</surname> <given-names>J. R.</given-names></name></person-group> (<year>2019</year>). <article-title>Chemosensory sensilla of the <italic>Drosophila</italic> wing express a candidate ionotropic pheromone receptor.</article-title> <source><italic>PLoS Biol.</italic></source> <volume>17</volume>:<issue>e2006619</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.02006619</pub-id></citation></ref>
<ref id="B18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hojo</surname> <given-names>M. K.</given-names></name> <name><surname>Ishii</surname> <given-names>K.</given-names></name> <name><surname>Sakura</surname> <given-names>M.</given-names></name> <name><surname>Yamaguchi</surname> <given-names>K.</given-names></name> <name><surname>Shigenobu</surname> <given-names>S.</given-names></name> <name><surname>Ozaki</surname> <given-names>M.</given-names></name></person-group> (<year>2015</year>). <article-title>Antennal RNA-sequencing analysis reveals evolutionary aspects of chemosensory proteins in the carpenter ant, <italic>Camponotus japonicus</italic>.</article-title> <source><italic>Sci. Rep.</italic></source> <volume>5</volume>:<issue>13541</issue>.</citation></ref>
<ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Holldobler</surname> <given-names>B.</given-names></name> <name><surname>Wilson</surname> <given-names>E. O.</given-names></name></person-group> (<year>1990</year>). <source><italic>The Ants.</italic></source> <publisher-loc>Cambridge, MA</publisher-loc>: <publisher-name>Belknap Press</publisher-name>.</citation></ref>
<ref id="B20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hou</surname> <given-names>Y.</given-names></name> <name><surname>Zhang</surname> <given-names>Y.</given-names></name> <name><surname>Gong</surname> <given-names>J.</given-names></name> <name><surname>Tian</surname> <given-names>S.</given-names></name> <name><surname>Li</surname> <given-names>J. W.</given-names></name> <name><surname>Dong</surname> <given-names>Z. M.</given-names></name><etal/></person-group> (<year>2016</year>). <article-title>Comparative proteomics analysis of silkworm hemolymph during the stages of metamorphosis via liquid chromatography and mass spectrometry.</article-title> <source><italic>Proteomics</italic></source> <volume>16</volume> <fpage>1421</fpage>&#x2013;<lpage>1431</lpage>. <pub-id pub-id-type="doi">10.1002/pmic.201500427</pub-id> <pub-id pub-id-type="pmid">26960012</pub-id></citation></ref>
<ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hu</surname> <given-names>P.</given-names></name> <name><surname>Tao</surname> <given-names>J.</given-names></name> <name><surname>Cui</surname> <given-names>M. M.</given-names></name> <name><surname>Gao</surname> <given-names>C. L.</given-names></name> <name><surname>Lu</surname> <given-names>P. F.</given-names></name> <name><surname>Luo</surname> <given-names>Y. Q.</given-names></name></person-group> (<year>2016</year>). <article-title>Antennal transcriptome analysis and expression profiles of odorant binding proteins in <italic>Eogystia hippophaecolus</italic> (Lepidoptera: Cossidae).</article-title> <source><italic>BMC Genom.</italic></source> <volume>17</volume>:<issue>651</issue>. <pub-id pub-id-type="doi">10.1186/s12864-016-3008-4</pub-id> <pub-id pub-id-type="pmid">27538507</pub-id></citation></ref>
<ref id="B22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ishida</surname> <given-names>Y.</given-names></name> <name><surname>Chiang</surname> <given-names>V.</given-names></name> <name><surname>Leal</surname> <given-names>W. S.</given-names></name></person-group> (<year>2002</year>). <article-title>Protein that makes sense in the Argentine ant.</article-title> <source><italic>Naturwissenschaften</italic></source> <volume>89</volume> <fpage>505</fpage>&#x2013;<lpage>507</lpage>. <pub-id pub-id-type="doi">10.1007/s00114-002-0368-1</pub-id> <pub-id pub-id-type="pmid">12451452</pub-id></citation></ref>
<ref id="B23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jacquin-Joly</surname> <given-names>E.</given-names></name> <name><surname>Vogt</surname> <given-names>R. G.</given-names></name> <name><surname>Francois</surname> <given-names>M. C.</given-names></name> <name><surname>Nagnan-Le Meillour</surname> <given-names>P.</given-names></name></person-group> (<year>2001</year>). <article-title>Functional and expression pattern analysis of chemosensory proteins expressed in antennae and pheromonal gland of <italic>Mamestra brassicae</italic>.</article-title> <source><italic>Chem. Sens.</italic></source> <volume>26</volume> <fpage>833</fpage>&#x2013;<lpage>844</lpage>. <pub-id pub-id-type="doi">10.1093/chemse/26.7.833</pub-id> <pub-id pub-id-type="pmid">11555479</pub-id></citation></ref>
<ref id="B24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jain</surname> <given-names>R.</given-names></name> <name><surname>Brockmann</surname> <given-names>A.</given-names></name></person-group> (<year>2020</year>). <article-title>Sex-specific molecular specialization and activity rhythm dependent gene expression in honey bee antennae.</article-title> <source><italic>J. Exper. Biol.</italic></source> <volume>2020</volume>:<issue>jeb.217406</issue>.</citation></ref>
<ref id="B25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jing</surname> <given-names>D.</given-names></name> <name><surname>Zhang</surname> <given-names>T.</given-names></name> <name><surname>Bai</surname> <given-names>S.</given-names></name> <name><surname>He</surname> <given-names>K.</given-names></name> <name><surname>Prabu</surname> <given-names>S.</given-names></name> <name><surname>Luan</surname> <given-names>J.</given-names></name><etal/></person-group> (<year>2020</year>). <article-title>Sexual-biased gene expression of olfactory-related genes in the antennae of <italic>Conogethes pinicolalis</italic> (Lepidoptera: Crambidae).</article-title> <source><italic>BMC Genomics</italic></source> <volume>21</volume>:<issue>244</issue>. <pub-id pub-id-type="doi">10.1186/1471-2164-15-244</pub-id> <pub-id pub-id-type="pmid">24678773</pub-id></citation></ref>
<ref id="B26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kather</surname> <given-names>R.</given-names></name> <name><surname>Martin</surname> <given-names>S. J.</given-names></name></person-group> (<year>2015</year>). <article-title>Evolution of cuticular hydrocarbons in the hymenoptera: a meta-analysis.</article-title> <source><italic>J. Chem. Ecol.</italic></source> <volume>41</volume> <fpage>871</fpage>&#x2013;<lpage>883</lpage>. <pub-id pub-id-type="doi">10.1007/s10886-015-0631-5</pub-id> <pub-id pub-id-type="pmid">26410609</pub-id></citation></ref>
<ref id="B27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kulmuni</surname> <given-names>J.</given-names></name> <name><surname>Wurm</surname> <given-names>Y.</given-names></name> <name><surname>Pamilo</surname> <given-names>P.</given-names></name></person-group> (<year>2013</year>). <article-title>Comparative genomics of chemosensory protein genes reveals rapid evolution and positive selection in ant-specific duplicates.</article-title> <source><italic>Heredity</italic></source> <volume>110</volume> <fpage>538</fpage>&#x2013;<lpage>547</lpage>. <pub-id pub-id-type="doi">10.1038/hdy.2012.122</pub-id> <pub-id pub-id-type="pmid">23403962</pub-id></citation></ref>
<ref id="B28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lartigue</surname> <given-names>A.</given-names></name> <name><surname>Campanacci</surname> <given-names>V.</given-names></name> <name><surname>Roussel</surname> <given-names>A.</given-names></name> <name><surname>Larsson</surname> <given-names>A. M.</given-names></name> <name><surname>Jones</surname> <given-names>T. A.</given-names></name> <name><surname>Tegoni</surname> <given-names>M.</given-names></name><etal/></person-group> (<year>2002</year>). <article-title>X-ray structure and ligand binding study of a moth chemosensory protein.</article-title> <source><italic>J. Biol. Chem.</italic></source> <volume>277</volume> <fpage>32094</fpage>&#x2013;<lpage>32098</lpage>. <pub-id pub-id-type="doi">10.1074/jbc.m204371200</pub-id> <pub-id pub-id-type="pmid">12068017</pub-id></citation></ref>
<ref id="B29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Leal</surname> <given-names>W. S.</given-names></name> <name><surname>Ishida</surname> <given-names>Y.</given-names></name></person-group> (<year>2008</year>). <article-title>GP-9s are ubiquitous proteins unlikely involved in olfactory mediation of social organization in the red imported fire ant, <italic>Solenopsis invicta</italic>.</article-title> <source><italic>PLoS One</italic></source> <volume>3</volume>:<issue>e003762</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0003762</pub-id> <pub-id pub-id-type="pmid">19018280</pub-id></citation></ref>
<ref id="B30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>H. L.</given-names></name> <name><surname>Ni</surname> <given-names>C. X.</given-names></name> <name><surname>Tan</surname> <given-names>J.</given-names></name> <name><surname>Zhang</surname> <given-names>L. Y.</given-names></name> <name><surname>Hu</surname> <given-names>F. L.</given-names></name></person-group> (<year>2016</year>). <article-title>Chemosensory proteins of the eastern honeybee, <italic>Apis cerana</italic>: identification, tissue distribution and olfactory related functional characterization.</article-title> <source><italic>Compar. Biochem. Physiol. B Biochem. Mol. Biol.</italic></source> <volume>194</volume> <fpage>11</fpage>&#x2013;<lpage>19</lpage>. <pub-id pub-id-type="doi">10.1016/j.cbpb.2015.11.014</pub-id> <pub-id pub-id-type="pmid">26773657</pub-id></citation></ref>
<ref id="B31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>X. M.</given-names></name> <name><surname>Zhu</surname> <given-names>X. Y.</given-names></name> <name><surname>He</surname> <given-names>P.</given-names></name> <name><surname>Xu</surname> <given-names>L.</given-names></name> <name><surname>Sun</surname> <given-names>L.</given-names></name> <name><surname>Chen</surname> <given-names>L.</given-names></name><etal/></person-group> (<year>2016</year>). <article-title>Molecular characterization and sex distribution of chemosensory receptor gene family based on transcriptome analysis of <italic>Scaeva pyrastri</italic>.</article-title> <source><italic>PLoS One</italic></source> <volume>11</volume>:<issue>e0155323</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0155323</pub-id> <pub-id pub-id-type="pmid">27171401</pub-id></citation></ref>
<ref id="B32"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>Z.</given-names></name> <name><surname>Smagghe</surname> <given-names>G.</given-names></name> <name><surname>Lei</surname> <given-names>Z. R.</given-names></name> <name><surname>Wang</surname> <given-names>J. J.</given-names></name></person-group> (<year>2016</year>). <article-title>Identification of male- and female-specific olfaction genes in antennae of the oriental fruit fly (<italic>Bactrocera dorsalis</italic>).</article-title> <source><italic>PLoS One</italic></source> <volume>11</volume>:<issue>e0147783</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0147783</pub-id> <pub-id pub-id-type="pmid">26845547</pub-id></citation></ref>
<ref id="B33"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>L&#x00F6;ytynoja</surname> <given-names>A.</given-names></name> <name><surname>Goldman</surname> <given-names>N.</given-names></name></person-group> (<year>2005</year>). <article-title>An algorithm for progressive multiple alignment of sequences with insertions.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>102</volume> <fpage>10557</fpage>&#x2013;<lpage>10562</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.0409137102</pub-id> <pub-id pub-id-type="pmid">16000407</pub-id></citation></ref>
<ref id="B34"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>L&#x00F6;ytynoja</surname> <given-names>A.</given-names></name> <name><surname>Goldman</surname> <given-names>N.</given-names></name></person-group> (<year>2010</year>). <article-title>webPRANK: a phylogeny-aware multiple sequence aligner with interactive alignment browser.</article-title> <source><italic>BMC Bioinform.</italic></source> <volume>11</volume>:<issue>579</issue>. <pub-id pub-id-type="doi">10.1186/1471-2164-15-579</pub-id> <pub-id pub-id-type="pmid">25011954</pub-id></citation></ref>
<ref id="B35"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ma</surname> <given-names>L.</given-names></name> <name><surname>Li</surname> <given-names>Z. Q.</given-names></name> <name><surname>Bian</surname> <given-names>L.</given-names></name> <name><surname>Cai</surname> <given-names>X. M.</given-names></name> <name><surname>Luo</surname> <given-names>Z. X.</given-names></name> <name><surname>Zhang</surname> <given-names>Y. J.</given-names></name><etal/></person-group> (<year>2016</year>). <article-title>Identification and comparative study of chemosensory genes related to host selection by legs transcriptome analysis in the tea geometrid <italic>Ectropis obliqua</italic>.</article-title> <source><italic>PLoS One</italic></source> <volume>11</volume>:<issue>e0149591</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0149591</pub-id> <pub-id pub-id-type="pmid">26930056</pub-id></citation></ref>
<ref id="B36"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Maleszka</surname> <given-names>J.</given-names></name> <name><surname>Foret</surname> <given-names>S.</given-names></name> <name><surname>Saint</surname> <given-names>R.</given-names></name> <name><surname>Maleszka</surname> <given-names>R.</given-names></name></person-group> (<year>2007</year>). <article-title>RNAi-induced phenotypes suggest a novel role for a chemosensory protein CSP5 in the development of embryonic integument in the honeybee (<italic>Apis mellifera</italic>).</article-title> <source><italic>Dev. Genes Evol.</italic></source> <volume>217</volume> <fpage>189</fpage>&#x2013;<lpage>196</lpage>. <pub-id pub-id-type="doi">10.1007/s00427-006-0127-y</pub-id> <pub-id pub-id-type="pmid">17216269</pub-id></citation></ref>
<ref id="B37"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>McKenzie</surname> <given-names>S. K.</given-names></name> <name><surname>Oxley</surname> <given-names>P. R.</given-names></name> <name><surname>Kronauer</surname> <given-names>D. J. C.</given-names></name></person-group> (<year>2014</year>). <article-title>Comparative genomics and transcriptomics in ants provide new insights into the evolution and function of odorant binding and chemosensory proteins.</article-title> <source><italic>BMC Genomics</italic></source> <volume>15</volume>:<issue>718</issue>. <pub-id pub-id-type="doi">10.1186/1471-2164-15-718</pub-id> <pub-id pub-id-type="pmid">25159315</pub-id></citation></ref>
<ref id="B38"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Miller</surname> <given-names>M. A.</given-names></name> <name><surname>Pfeiffer</surname> <given-names>W.</given-names></name> <name><surname>Schwartz</surname> <given-names>T.</given-names></name></person-group> (<year>2010</year>). &#x201C;<article-title>Creating the CIPRES Science Gateway for inference of large phylogenetic trees</article-title>,&#x201D; in <source><italic>Proceedings of the Gateway Computing Environments Workshop (GCE)</italic></source>, <publisher-loc>New Orleans, LA</publisher-loc>.</citation></ref>
<ref id="B39"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mitaka</surname> <given-names>Y.</given-names></name> <name><surname>Kobayashi</surname> <given-names>K.</given-names></name> <name><surname>Mikheyev</surname> <given-names>A.</given-names></name> <name><surname>Tin</surname> <given-names>M. M. Y.</given-names></name> <name><surname>Watanabe</surname> <given-names>Y.</given-names></name> <name><surname>Matsuura</surname> <given-names>K.</given-names></name></person-group> (<year>2016</year>). <article-title>Caste-specific and sex-specific expression of chemoreceptor genes in a termite.</article-title> <source><italic>PLoS One</italic></source> <volume>11</volume>:<issue>e0146125</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0146125</pub-id> <pub-id pub-id-type="pmid">26760975</pub-id></citation></ref>
<ref id="B40"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nomura</surname> <given-names>A.</given-names></name> <name><surname>Kawasaki</surname> <given-names>K.</given-names></name> <name><surname>Kubo</surname> <given-names>T.</given-names></name> <name><surname>Natori</surname> <given-names>S.</given-names></name></person-group> (<year>1992</year>). <article-title>Purification and localization of P10, a novel protein that increases in nymphal regenerating legs of <italic>Periplaneta americana</italic> (American Cockroach).</article-title> <source><italic>Int. J. Dev. Biol.</italic></source> <volume>36</volume> <fpage>391</fpage>&#x2013;<lpage>398</lpage>.</citation></ref>
<ref id="B41"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ozaki</surname> <given-names>M.</given-names></name> <name><surname>Wada-Katsumata</surname> <given-names>A.</given-names></name> <name><surname>Fujikawa</surname> <given-names>K.</given-names></name> <name><surname>Iwasaki</surname> <given-names>M.</given-names></name> <name><surname>Yokohari</surname> <given-names>F.</given-names></name> <name><surname>Satoji</surname> <given-names>Y.</given-names></name><etal/></person-group> (<year>2005</year>). <article-title>Ant nestmate and non-nestmate discrimination by a chemosensory sensillum.</article-title> <source><italic>Science</italic></source> <volume>309</volume> <fpage>311</fpage>&#x2013;<lpage>314</lpage>. <pub-id pub-id-type="doi">10.1126/science.1105244</pub-id> <pub-id pub-id-type="pmid">15947139</pub-id></citation></ref>
<ref id="B42"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pelosi</surname> <given-names>P.</given-names></name> <name><surname>Iovinella</surname> <given-names>I.</given-names></name> <name><surname>Felicioli</surname> <given-names>A.</given-names></name> <name><surname>Dani</surname> <given-names>F. R.</given-names></name></person-group> (<year>2014</year>). <article-title>Soluble proteins of chemical communication: an overview across arthropods.</article-title> <source><italic>Front. Physiol.</italic></source> <volume>5</volume>:<issue>320</issue>. <pub-id pub-id-type="doi">10.3389/fphys.2014.00320</pub-id> <pub-id pub-id-type="pmid">25221516</pub-id></citation></ref>
<ref id="B43"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Picimbon</surname> <given-names>J. F.</given-names></name> <name><surname>Dietrich</surname> <given-names>K.</given-names></name> <name><surname>Angeli</surname> <given-names>S.</given-names></name> <name><surname>Scaloni</surname> <given-names>A.</given-names></name> <name><surname>Krieger</surname> <given-names>J.</given-names></name> <name><surname>Breer</surname> <given-names>H.</given-names></name><etal/></person-group> (<year>2000</year>). <article-title>Purification and molecular cloning of chemosensory proteins from <italic>Bombyx mori</italic>.</article-title> <source><italic>Arch. Insect. Biochem. Physiol.</italic></source> <volume>44</volume> <fpage>120</fpage>&#x2013;<lpage>129</lpage>. <pub-id pub-id-type="doi">10.1002/1520-6327(200007)44:3&#x003C;120::aid-arch3&#x003E;3.0.co;2-h</pub-id></citation></ref>
<ref id="B44"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Renthal</surname> <given-names>R.</given-names></name> <name><surname>Velasquez</surname> <given-names>D.</given-names></name> <name><surname>Olmos</surname> <given-names>D.</given-names></name> <name><surname>Hampton</surname> <given-names>J.</given-names></name> <name><surname>Wergin</surname> <given-names>W. P.</given-names></name></person-group> (<year>2003</year>). <article-title>Structure and distribution of antennal sensilla of the red imported fire ant.</article-title> <source><italic>Micron</italic></source> <volume>34</volume> <fpage>405</fpage>&#x2013;<lpage>413</lpage>. <pub-id pub-id-type="doi">10.1016/s0968-4328(03)00050-7</pub-id></citation></ref>
<ref id="B45"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ross</surname> <given-names>K. G.</given-names></name> <name><surname>Fletcher</surname> <given-names>D. J. C.</given-names></name></person-group> (<year>1985</year>). <article-title>Comparative study of genetic and social structure in two forms of the Fire Ant <italic>Solenopsis invicta</italic> (Hymenoptera, Formicidae).</article-title> <source><italic>Behav. Ecol. Sociobiol.</italic></source> <volume>17</volume> <fpage>349</fpage>&#x2013;<lpage>356</lpage>. <pub-id pub-id-type="doi">10.1007/bf00293212</pub-id></citation></ref>
<ref id="B46"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shoemaker</surname> <given-names>D.</given-names></name> <name><surname>Ascunce</surname> <given-names>M. S.</given-names></name> <name><surname>Yang</surname> <given-names>C. C.</given-names></name> <name><surname>Oakey</surname> <given-names>J.</given-names></name> <name><surname>Calcaterra</surname> <given-names>L.</given-names></name> <name><surname>Wu</surname> <given-names>W. J.</given-names></name><etal/></person-group> (<year>2011</year>). <article-title>Global invasion history of the fire ant <italic>Solenopsis invicta</italic>.</article-title> <source><italic>Science</italic></source> <volume>331</volume> <fpage>1066</fpage>&#x2013;<lpage>1068</lpage>.</citation></ref>
<ref id="B47"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Slone</surname> <given-names>J. D.</given-names></name> <name><surname>Pask</surname> <given-names>G. M.</given-names></name> <name><surname>Ferguson</surname> <given-names>S. T.</given-names></name> <name><surname>Millar</surname> <given-names>J. G.</given-names></name> <name><surname>Berger</surname> <given-names>S. L.</given-names></name> <name><surname>Reinberg</surname> <given-names>D.</given-names></name><etal/></person-group> (<year>2017</year>). <article-title>Functional characterization of odorant receptors in the ponerine ant, <italic>Harpegnathos saltator</italic>.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>114</volume> <fpage>8586</fpage>&#x2013;<lpage>8591</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1704647114</pub-id> <pub-id pub-id-type="pmid">28696298</pub-id></citation></ref>
<ref id="B48"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Stamatakis</surname> <given-names>A.</given-names></name></person-group> (<year>2014</year>). <article-title>RAxML Version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies.</article-title> <source><italic>Bioinformatics</italic></source> <volume>30</volume> <fpage>1312</fpage>&#x2013;<lpage>1313</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btu033</pub-id> <pub-id pub-id-type="pmid">24451623</pub-id></citation></ref>
<ref id="B49"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tamura</surname> <given-names>K.</given-names></name> <name><surname>Stecher</surname> <given-names>G.</given-names></name> <name><surname>Peterson</surname> <given-names>D.</given-names></name> <name><surname>Filipski</surname> <given-names>A.</given-names></name> <name><surname>Kumar</surname> <given-names>S.</given-names></name></person-group> (<year>2013</year>). <article-title>MEGA6: molecular evolutionary genetics analysis version 6.0.</article-title> <source><italic>Mol. Biol. Evol.</italic></source> <volume>30</volume> <fpage>2725</fpage>&#x2013;<lpage>2729</lpage>. <pub-id pub-id-type="doi">10.1093/molbev/mst197</pub-id> <pub-id pub-id-type="pmid">24132122</pub-id></citation></ref>
<ref id="B50"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tschinkel</surname> <given-names>W. R.</given-names></name></person-group> (<year>2006</year>). <source><italic>The Fire Ants.</italic></source> <publisher-loc>Cambridge, MA</publisher-loc>: <publisher-name>Belknap Press of Harvard University Press</publisher-name>.</citation></ref>
<ref id="B51"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Vieira</surname> <given-names>F. G.</given-names></name> <name><surname>Rozas</surname> <given-names>J.</given-names></name></person-group> (<year>2011</year>). <article-title>Comparative genomics of the odorant-binding and chemosensory protein gene families across the arthropoda: origin and evolutionary history of the chemosensory system.</article-title> <source><italic>Genome Biol. Evol.</italic></source> <volume>3</volume> <fpage>476</fpage>&#x2013;<lpage>490</lpage>. <pub-id pub-id-type="doi">10.1093/gbe/evr033</pub-id> <pub-id pub-id-type="pmid">21527792</pub-id></citation></ref>
<ref id="B52"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>J.</given-names></name> <name><surname>Wurm</surname> <given-names>Y.</given-names></name> <name><surname>Nipitwattanaphon</surname> <given-names>M.</given-names></name> <name><surname>Riba-Grognuz</surname> <given-names>O.</given-names></name> <name><surname>Huang</surname> <given-names>Y. C.</given-names></name> <name><surname>Shoemaker</surname> <given-names>D.</given-names></name><etal/></person-group> (<year>2013</year>). <article-title>A Y-like social chromosome causes alternative colony organization in fire ants.</article-title> <source><italic>Nature</italic></source> <volume>493</volume> <fpage>664</fpage>&#x2013;<lpage>668</lpage>. <pub-id pub-id-type="doi">10.1038/nature11832</pub-id> <pub-id pub-id-type="pmid">23334415</pub-id></citation></ref>
<ref id="B53"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Whelan</surname> <given-names>J. A.</given-names></name> <name><surname>Russell</surname> <given-names>N. B.</given-names></name> <name><surname>Whelan</surname> <given-names>M. A.</given-names></name></person-group> (<year>2003</year>). <article-title>A method for the absolute quantification of cDNA using real-time PCR.</article-title> <source><italic>J. Immunol. Methods</italic></source> <volume>278</volume> <fpage>261</fpage>&#x2013;<lpage>269</lpage>. <pub-id pub-id-type="doi">10.1016/s0022-1759(03)00223-0</pub-id></citation></ref>
<ref id="B54"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wurm</surname> <given-names>Y.</given-names></name> <name><surname>Wang</surname> <given-names>J.</given-names></name> <name><surname>Riba-Grognuz</surname> <given-names>O.</given-names></name> <name><surname>Corona</surname> <given-names>M.</given-names></name> <name><surname>Nygaard</surname> <given-names>S.</given-names></name> <name><surname>Hunt</surname> <given-names>B. G.</given-names></name><etal/></person-group> (<year>2011</year>). <article-title>The genome of the fire ant <italic>Solenopsis invicta</italic>.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>108</volume> <fpage>5679</fpage>&#x2013;<lpage>5684</lpage>.</citation></ref>
<ref id="B55"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xu</surname> <given-names>Y. L.</given-names></name> <name><surname>He</surname> <given-names>P.</given-names></name> <name><surname>Zhang</surname> <given-names>L.</given-names></name> <name><surname>Fang</surname> <given-names>S. Q.</given-names></name> <name><surname>Dong</surname> <given-names>S. L.</given-names></name> <name><surname>Zhang</surname> <given-names>Y. J.</given-names></name><etal/></person-group> (<year>2009</year>). <article-title>Large-scale identification of odorant-binding proteins and chemosensory proteins from expressed sequence tags in insects.</article-title> <source><italic>BMC Genomics</italic></source> <volume>10</volume>:<issue>632</issue>. <pub-id pub-id-type="doi">10.1186/1471-2164-15-632</pub-id> <pub-id pub-id-type="pmid">25070353</pub-id></citation></ref>
<ref id="B56"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xue</surname> <given-names>W. X.</given-names></name> <name><surname>Fan</surname> <given-names>J.</given-names></name> <name><surname>Zhang</surname> <given-names>Y.</given-names></name> <name><surname>Xu</surname> <given-names>Q. X.</given-names></name> <name><surname>Han</surname> <given-names>Z. L.</given-names></name> <name><surname>Sun</surname> <given-names>J. R.</given-names></name><etal/></person-group> (<year>2016</year>). <article-title>Identification and expression analysis of candidate odorant-binding protein and chemosensory protein genes by antennal transcriptome of <italic>Sitobion avenae</italic>.</article-title> <source><italic>PLoS One</italic></source> <volume>11</volume>:<issue>e0161839</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0161839</pub-id> <pub-id pub-id-type="pmid">27561107</pub-id></citation></ref>
<ref id="B57"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yi</surname> <given-names>X.</given-names></name> <name><surname>Zhao</surname> <given-names>H. M.</given-names></name> <name><surname>Dong</surname> <given-names>X. L.</given-names></name> <name><surname>Wang</surname> <given-names>P. D.</given-names></name> <name><surname>Hu</surname> <given-names>M. Y.</given-names></name> <name><surname>Zhong</surname> <given-names>G. H.</given-names></name></person-group> (<year>2013</year>). <article-title>BdorCSP2 Is Important for antifeed and oviposition-deterring activities induced by rhodojaponin-III against <italic>Bactrocera dorsalis</italic>.</article-title> <source><italic>PLoS One</italic></source> <volume>8</volume>:<issue>e077295</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0077295</pub-id> <pub-id pub-id-type="pmid">24155937</pub-id></citation></ref>
<ref id="B58"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>W.</given-names></name> <name><surname>Wanchoo</surname> <given-names>A.</given-names></name> <name><surname>Ortiz-Urquiza</surname> <given-names>A.</given-names></name> <name><surname>Xia</surname> <given-names>Y.</given-names></name> <name><surname>Keyhani</surname> <given-names>N. O.</given-names></name></person-group> (<year>2016</year>). <article-title>Tissue, developmental, and caste-specific expression of odorant binding proteins in a eusocial insect, the red imported fire ant, <italic>Solenopsis invicta</italic>.</article-title> <source><italic>Sci. Rep.</italic></source> <volume>6</volume>:<issue>35452</issue>.</citation></ref>
<ref id="B59"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname> <given-names>J. J.</given-names></name> <name><surname>Kan</surname> <given-names>Y. C.</given-names></name> <name><surname>Antoniw</surname> <given-names>J.</given-names></name> <name><surname>Pickett</surname> <given-names>J. A.</given-names></name> <name><surname>Field</surname> <given-names>L. M.</given-names></name></person-group> (<year>2006</year>). <article-title>Genome and EST analyses and expression of a gene family with putative functions in insect chemoreception.</article-title> <source><italic>Chem. Sens.</italic></source> <volume>31</volume> <fpage>453</fpage>&#x2013;<lpage>465</lpage>. <pub-id pub-id-type="doi">10.1093/chemse/bjj050</pub-id> <pub-id pub-id-type="pmid">16581978</pub-id></citation></ref>
<ref id="B60"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname> <given-names>X. H.</given-names></name> <name><surname>Ban</surname> <given-names>L. P.</given-names></name> <name><surname>Iovinella</surname> <given-names>I.</given-names></name> <name><surname>Zhao</surname> <given-names>L. J.</given-names></name> <name><surname>Gao</surname> <given-names>Q.</given-names></name> <name><surname>Felicioli</surname> <given-names>A.</given-names></name><etal/></person-group> (<year>2013</year>). <article-title>Diversity, abundance, and sex-specific expression of chemosensory proteins in the reproductive organs of the locust <italic>Locusta migratoria manilensis</italic>.</article-title> <source><italic>Biol. Chem.</italic></source> <volume>394</volume> <fpage>43</fpage>&#x2013;<lpage>54</lpage>. <pub-id pub-id-type="doi">10.1515/hsz-2012-0114</pub-id> <pub-id pub-id-type="pmid">23096575</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn id="footnote1">
<label>1</label>
<p><ext-link ext-link-type="uri" xlink:href="https://software.broadinstitute.org/Morpheus">https://software.broadinstitute.org/Morpheus</ext-link></p></fn>
<fn id="footnote2">
<label>2</label>
<p><ext-link ext-link-type="uri" xlink:href="http://www.cbs.dtu.dk/services/SignalP/">http://www.cbs.dtu.dk/services/SignalP/</ext-link></p></fn>
<fn id="footnote3">
<label>3</label>
<p><ext-link ext-link-type="uri" xlink:href="http://gsds.cbi.pku.edu.cn">http://gsds.cbi.pku.edu.cn</ext-link></p></fn>
<fn id="footnote4">
<label>4</label>
<p><ext-link ext-link-type="uri" xlink:href="http://meme-suite.org/tools/meme">http://meme-suite.org/tools/meme</ext-link></p></fn>
</fn-group>
</back>
</article>
