Multigene Family Evolution: Perspectives from Insect Chemoreceptors

Genomic content of chemosensory genes correlates with host range in wood-boring beetles (Dendroctonus ponderosae, Agrilus planipennis, and Anoplophora glabripennis)

Background

Olfaction and gustation underlie behaviors that are crucial for insect fitness, such as host and mate selection. The detection of semiochemicals is mediated via proteins from large and rapidly evolving chemosensory gene families; however, the links between a species’ ecology and the diversification of these genes remain poorly understood. Hence, we annotated the chemosensory genes from genomes of select wood-boring coleopterans, and compared the gene repertoires from stenophagous species with those from polyphagous species.

Results

We annotated 86 odorant receptors (ORs), 60 gustatory receptors (GRs), 57 ionotropic receptors (IRs), 4 sensory neuron membrane proteins (SNMPs), 36 odorant binding proteins (OBPs), and 11 chemosensory proteins (CSPs) in the mountain pine beetle (Dendroctonus ponderosae), and 47 ORs, 30 GRs, 31 IRs, 4 SNMPs, 12 OBPs, and 14 CSPs in the emerald ash borer (Agrilus planipennis). Four SNMPs and 17 CSPs were annotated in the polyphagous wood-borer Anoplophora glabripennis. The gene repertoires in the stenophagous D. ponderosae and A. planipennis are reduced compared with those in the polyphagous A. glabripennis and T. castaneum, which is largely manifested through small gene lineage expansions and entire lineage losses. Alternative splicing of GR genes was limited in D. ponderosae and apparently absent in A. planipennis, which also seems to have lost one carbon dioxide receptor (GR1). A. planipennis has two SNMPs, which are related to SNMP3 in T. castaneum. D. ponderosae has two alternatively spliced OBP genes, a novel OBP “tetramer”, and as many as eleven IR75 members. Simple orthology was generally rare in beetles; however, we found one clade with orthologues of putative bitter-taste GRs (named the “GR215 clade”), and conservation of IR60a from Drosophila melanogaster.

Conclusions

Our genome annotations represent important quantitative and qualitative improvements of the original datasets derived from transcriptomes of D. ponderosae and A. planipennis, facilitating evolutionary analysis of chemosensory genes in the Coleoptera where only a few genomes were previously annotated. Our analysis suggests a correlation between chemosensory gene content and host specificity in beetles. Future studies should include additional species to consolidate this correlation, and functionally characterize identified proteins as an important step towards improved control of these pests.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6054-x) contains supplementary material, which is available to authorized users.

Downregulation of Orco and 5-HTT Alters Nestmate Discrimination in the Subterranean Termite Odontotermes formosanus (Shiraki)

Nestmate discrimination allows social insects to recognize nestmates from non-nestmates using colony-specific chemosensory cues, which typically evoke aggressive behavior toward non-nestmates. Functional analysis of genes associated with nestmate discrimination has been primarily focused on inter-colonial discrimination in Hymenopterans, and parallel studies in termites, however, are grossly lacking. To fill this gap, we investigated the role of two genes, Orco and 5-HTT, associated with chemosensation and neurotransmission respectively, in nestmate discrimination in a highly eusocial subterranean termite, Odontotermes formosanus (Shiraki). We hypothesized that knocking down of these genes will compromise the nestmate recognition and lead to the antagonistic behavior. To test this hypothesis, we carried out (1) an in vivo RNAi to suppress the expression of Orco and 5-HTT, respectively, (2) a validation study to examine the knockdown efficiency, and finally, (3) a behavioral assay to document the phenotypic impacts/behavioral consequences. As expected, the suppression of either of these two genes elevated stress level (e.g., vibrations and retreats), and led to aggressive behaviors (e.g., biting) in O. formosanus workers toward their nestmates, suggesting both Orco and 5-HTT can modulate nestmate discrimination in termites. This research links chemosensation and neurotransmission with nestmate discrimination at the genetic basis, and lays the foundation for functional analyses of nestmate discrimination in termites.

U, Kulkarni R, Löfstedt C, Sowdhamini R, Olsson SB. A Functional Agonist of Insect Olfactory Receptors: Behavior, Physiology and Structure

Chemical signaling is ubiquitous and employs a variety of receptor types to detect the cacophony of molecules relevant for each living organism. Insects, our most diverse taxon, have evolved unique olfactory receptors with as little as 10% sequence identity between receptor types. We have identified a promiscuous volatile, 2-methyltetrahydro-3-furanone (coffee furanone), that elicits chemosensory and behavioral activity across multiple insect orders and receptors. In vivo and in vitro physiology showed that coffee furanone was detected by roughly 80% of the recorded neurons expressing the insect-specific olfactory receptor complex in the antenna of Drosophila melanogaster, at concentrations similar to other known, and less promiscuous, ligands. Neurons expressing specialized receptors, other chemoreceptor types, or mutants lacking the complex entirely did not respond to this compound. This indicates that coffee furanone is a promiscuous ligand for the insect olfactory receptor complex itself and did not induce non-specific cellular responses. In addition, we present homology modeling and docking studies with selected olfactory receptors that suggest conserved interaction regions for both coffee furanone and known ligands. Apart from its physiological activity, this known food additive elicits a behavioral response for several insects, including mosquitoes, flies, and cockroaches. A broad-scale behaviorally active molecule non-toxic to humans thus has significant implications for health and agriculture. Coffee furanone serves as a unique tool to unlock molecular, physiological, and behavioral relationships across this diverse receptor family and animal taxa.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Understanding insect foraging in complex habitats by comparing trophic levels: insights from specialist host-parasitoid-hyperparasitoid systems

Insects typically forage in complex habitats in which their resources are surrounded by non-resources. For herbivores, pollinators, parasitoids, and higher level predators research has focused on how specific trophic levels filter and integrate information from cues in their habitat to locate resources. However, these insights frequently build specific theory per trophic level and seldom across trophic levels. Here, we synthesize advances in understanding of insect foraging behavior in complex habitats by comparing trophic levels in specialist host-parasitoid-hyperparasitoid systems. We argue that resources may become less apparent to foraging insects when they are member of higher trophic levels and hypothesize that higher trophic level organisms require a larger number of steps in their foraging decisions. We identify important knowledge gaps of information integration strategies by insects that belong to higher trophic levels.

Molecular Evolution of the Major Arthropod Chemoreceptor Gene Families

The evolutionary origins of the three major families of chemoreceptors in arthropods-the odorant receptor (OR), gustatory receptor (GR), and ionotropic receptor (IR) families-occurred at the base of the Insecta, Animalia, and Protostomia, respectively. Comparison of receptor family sizes across arthropods reveals a generally positive correlation with their widely disparate complexity of chemical ecology. Closely related species reveal the ongoing processes of gene family evolution, including gene duplication, divergence, pseudogenization, and loss, that mediate these larger patterns. Sets of paralogous receptors within species reveal positive selection on amino acids in regions likely to contribute to ligand binding and specificity. Ligands of many ORs and some GRs and IRs have been identified; however, ligand identification for many more chemoreceptors is needed, as are structures for the OR/GR superfamily, to improve our understanding of the molecular evolution of these ecologically important receptors in arthropods.

The genome of the water strider Gerris buenoi reveals expansions of gene repertoires associated with adaptations to life on the water

BACKGROUND

Having conquered water surfaces worldwide, the semi-aquatic bugs occupy ponds, streams, lakes, mangroves, and even open oceans. The diversity of this group has inspired a range of scientific studies from ecology and evolution to developmental genetics and hydrodynamics of fluid locomotion. However, the lack of a representative water strider genome hinders our ability to more thoroughly investigate the molecular mechanisms underlying the processes of adaptation and diversification within this group.

RESULTS

Here we report the sequencing and manual annotation of the Gerris buenoi (G. buenoi) genome; the first water strider genome to be sequenced thus far. The size of the G. buenoi genome is approximately 1,000 Mb, and this sequencing effort has recovered 20,949 predicted protein-coding genes. Manual annotation uncovered a number of local (tandem and proximal) gene duplications and expansions of gene families known for their importance in a variety of processes associated with morphological and physiological adaptations to a water surface lifestyle. These expansions may affect key processes associated with growth, vision, desiccation resistance, detoxification, olfaction and epigenetic regulation. Strikingly, the G. buenoi genome contains three insulin receptors, suggesting key changes in the rewiring and function of the insulin pathway. Other genomic changes affecting with opsin genes may be associated with wavelength sensitivity shifts in opsins, which is likely to be key in facilitating specific adaptations in vision for diverse water habitats.

CONCLUSIONS

Our findings suggest that local gene duplications might have played an important role during the evolution of water striders. Along with these findings, the sequencing of the G. buenoi genome now provides us the opportunity to pursue exciting research opportunities to further understand the genomic underpinnings of traits associated with the extreme body plan and life history of water striders.

Whole-head recording of chemosensory activity in the marine annelid Platynereis dumerilii

Chemical detection is key to various behaviours in both marine and terrestrial animals. Marine species, though highly diverse, have been underrepresented so far in studies on chemosensory systems, and our knowledge mostly concerns the detection of airborne cues. A broader comparative approach is therefore desirable. Marine annelid worms with their rich behavioural repertoire represent attractive models for chemosensation. Here, we study the marine worm Platynereis dumerilii to provide the first comprehensive investigation of head chemosensory organ physiology in an annelid. By combining microfluidics and calcium imaging, we record neuronal activity in the entire head of early juveniles upon chemical stimulation. We find that Platynereis uses four types of organs to detect stimuli such as alcohols, esters, amino acids and sugars. Antennae are the main chemosensory organs, compared to the more differentially responding nuchal organs or palps. We report chemically evoked activity in possible downstream brain regions including the mushroom bodies (MBs), which are anatomically and molecularly similar to insect MBs. We conclude that chemosensation is a major sensory modality for marine annelids and propose early Platynereis juveniles as a model to study annelid chemosensory systems.

Antennal Transcriptome Analysis of the Chemosensory Gene Families From Trichoptera and Basal Lepidoptera

The chemosensory gene families of insects encode proteins that are crucial for host location, mate finding, oviposition, and avoidance behaviors. The insect peripheral chemosensory system comprises odorant receptors (ORs), gustatory receptors (GRs), ionotropic receptors (IRs), odorant binding proteins (OBPs), chemosensory proteins (CSPs), and sensory neuron membrane proteins (SNMPs). These protein families have been identified from a large number of insect species, however, they still remain unidentified from several taxa that could provide important clues to their evolution. These taxa include older lepidopteran lineages and the sister order of Lepidoptera, Trichoptera (caddisflies). Studies of these insects should improve evolutionary analyses of insect chemoreception, and in particular shed light on the origin of certain lepidopteran protein subfamilies. These include the pheromone receptors (PRs) in the "PR clade", the pheromone binding proteins (PBPs), general odorant binding proteins (GOBPs), and certain presumably Lepidoptera-specific IR subfamilies. Hence, we analyzed antennal transcriptomes from Rhyacophila nubila (Trichoptera), Eriocrania semipurpurella, and Lampronia capitella (representing two old lepidopteran lineages). We report 37 ORs, 17 IRs, 9 GRs, 30 OBPs, 7 CSPs, and 2 SNMPs in R. nubila; 37 ORs, 17 IRs, 3 GRs, 23 OBPs, 14 CSPs, and 2 SNMPs in E. semipurpurella; and 53 ORs, 20 IRs, 5 GRs, 29 OBPs, 17 CSPs, and 3 SNMPs in L. capitella. We identified IR members of the "Lepidoptera-specific" subfamilies IR1 and IR87a also in R. nubila, demonstrating that these IRs also occur in Trichoptera. Members of the GOBP subfamily were only found in the two lepidopterans. ORs grouping within the PR clade, as well as PBPs, were only found in L. capitella, a species that in contrast to R. nubila and E. semipurpurella uses a so-called Type I pheromone similar to the pheromones of most species of the derived Lepidoptera (Ditrysia). Thus, in addition to providing increased coverage for evolutionary analyses of chemoreception in insects, our findings suggest that certain subfamilies of chemosensory genes have evolved in parallel with the transition of sex pheromone types in Lepidoptera. In addition, other chemoreceptor subfamilies show a broader taxonomic occurrence than hitherto acknowledged.

Chemoreceptor proteins in the Caribbean spiny lobster, Panulirus argus: Expression of Ionotropic Receptors, Gustatory Receptors, and TRP channels in two chemosensory organs and brain

The spiny lobster, Panulirus argus, has two classes of chemosensilla representing “olfaction” and “distributed chemoreception,” as is typical for decapod crustaceans. Olfactory sensilla are found exclusively on antennular lateral flagella and are innervated only by olfactory receptor neurons (ORNs) that project into olfactory lobes organized into glomeruli in the brain. Distributed chemoreceptor sensilla are found on all body surfaces including the antennular lateral flagella (LF) and walking leg dactyls (dactyls), and are innervated by both chemoreceptor neurons (CRNs) and mechanoreceptor neurons that project into somatotopically organized neuropils. Here, we examined expression of three classes of chemosensory genes in transcriptomes of the LF (with ORNs and CRNs), dactyls (with only CRNs), and brain of P. argus: Ionotropic Receptors (IRs), which are related to ionotropic glutamate receptors and found in all protostomes including crustaceans; Gustatory Receptors (GRs), which are ionotropic receptors that are abundantly expressed in insects but more restricted in crustaceans; and Transient Receptor Potential (TRP) channels, a diverse set of sensor-channels that include several chemosensors in diverse animals. We identified 108 IRs, one GR, and 18 homologues representing all seven subfamilies of TRP channels. The number of IRs expressed in the LF is far greater than in dactyls, possibly reflecting the contribution of receptor proteins associated with the ORNs beyond those associated with CRNs. We found co-receptor IRs (IR8a, IR25a, IR76b, IR93a) and conserved IRs (IR21a, IR40a) in addition to the numerous divergent IRs in the LF, dactyl, and brain. Immunocytochemistry showed that IR25a is expressed in ORNs, CRNs, and a specific type of cell located in the brain near the olfactory lobes. While the function of IRs, TRP channels, and the GR was not explored, our results suggest that P. argus has an abundance of diverse putative chemoreceptor proteins that it may use in chemoreception.

Positive selection at sites of chemosensory genes is associated with the recent divergence and local ecological adaptation in cactophilic Drosophila

Background

Adaptation to new hosts in phytophagous insects often involves mechanisms of host recognition by genes of sensory pathways. Most often the molecular evolution of sensory genes has been explained in the context of the birth-and-death model. The role of positive selection is less understood, especially associated with host adaptation and specialization. Here we aim to contribute evidence for this latter hypothesis by considering the case of Drosophila mojavensis, a species with an evolutionary history shaped by multiple host shifts in a relatively short time scale, and its generalist sister species, D. arizonae.

Results

We used a phylogenetic and population genetic analysis framework to test for positive selection in a subset of four chemoreceptor genes, one gustatory receptor (Gr) and three odorant receptors (Or), for which their expression has been previously associated with host shifts. We found strong evidence of positive selection at several amino acid sites in all genes investigated, most of which exhibited changes predicted to cause functional effects in these transmembrane proteins. A significant portion of the sites identified as evolving positively were largely found in the cytoplasmic region, although a few were also present in the extracellular domains.

Conclusions

The pattern of substitution observed suggests that some of these changes likely had an effect on signal transduction as well as odorant recognition and protein-protein interactions. These findings support the role of positive selection in shaping the pattern of variation at chemosensory receptors, both during the specialization onto one or a few related hosts, but as well as during the evolution and adaptation of generalist species into utilizing several hosts.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1250-x) contains supplementary material, which is available to authorized users.

Enormous expansion of the chemosensory gene repertoire in the omnivorous German cockroach Blattella germanica

The acquisition of genome sequences from a wide range of insects and other arthropods has revealed a broad positive correlation between the complexity of their chemical ecology and the size of their chemosensory gene repertoire. The German cockroach Blattella germanica is an extreme omnivore and has the largest chemosensory gene repertoire known for an arthropod, exceeding even the highly polyphagous spider mite Tetranychus urticae. While the Odorant Receptor family is not particularly large, with 123 genes potentially encoding 134 receptors (105 intact), the Gustatory Receptor family is greatly expanded to 431 genes potentially encoding 545 receptors (483 intact), the largest known for insects and second only to the spider mite. The Ionotropic Receptor family of olfactory and gustatory receptors is vastly expanded to at least 897 genes (604 intact), the largest size known in arthropods, far surpassing the 150 known from the dampwood termite Zootermopsis nevadensis. Commensurately, the Odorant Binding Protein family is expanded to the largest known for insects at 109 genes (all intact). Comparison with the far more specialized, but phylogenetically related termite, within the Dictyoptera, reveals considerable gene losses from the termite, and massive species-specific gene expansions in the cockroach. The cockroach has lost function of 11%-41% of these three chemoreceptor gene families to pseudogenization, and most of these are young events, implying rapid turnover of genes along with these major expansions, presumably in response to changes in its chemical ecology.

A foreleg transcriptome for Ixodes scapularis ticks: Candidates for chemoreceptors and binding proteins that might be expressed in the sensory Haller's organ

Little is known about the molecular basis for the olfactory capabilities of the sensory Haller's organ on the forelegs of ticks. We first expanded the known repertoire of Ionotropic Receptors (IRs), a variant lineage of the ionotropic glutamate receptors, encoded by the black-legged Ixodes scapularis genome from 15 to 125. We then undertook a transcriptome study of fore- and hind-legs of this tick in an effort to identify candidate chemoreceptors differentially expressed in forelegs as likely to be involved in Haller's organ functions. We primarily identified members of the IR family, specifically Ir25a and Ir93a, as highly and differentially expressed in forelegs. Several other IRs, as well as a few members of the gustatory receptor family, were expressed at low levels in forelegs and might contribute to the sensory function of Haller's organ. In addition, we identified eight small families of secreted proteins, with sets of conserved cysteines, which might function as binding proteins. The genes encoding these Microplusin-Like proteins and two previously described Odorant Binding Protein-Like proteins share a common exon-intron structure, suggesting that they all evolved from a common ancestor and represent an independent origin of binding proteins with potential roles comparable to the ChemoSensory Proteins and Odorant Binding Proteins of insects. We also found two Niemann-Pick Type C2 proteins with foreleg-biased expression, however we were unable to detect foreleg-biased expression of a G-Protein-Coupled pathway previously proposed to mediate olfaction in the tick Haller's organ.

Characterization of Odorant Receptors from a Non-ditrysian Moth, Eriocrania semipurpurella Sheds Light on the Origin of Sex Pheromone Receptors in Lepidoptera

Pheromone receptors (PRs) are essential in moths to detect sex pheromones for mate finding. However, it remains unknown from which ancestral proteins these specialized receptors arose. The oldest lineages of moths, so-called non-ditrysian moths, use short-chain pheromone components, secondary alcohols, or ketones, so called Type 0 pheromones that are similar to many common plant volatiles. It is, therefore, possible that receptors for these ancestral pheromones evolved from receptors detecting plant volatiles. Hence, we identified the odorant receptors (ORs) from a non-ditrysian moth, Eriocrania semipurpurella (Eriocraniidae, Lepidoptera), and performed functional characterization of ORs using HEK293 cells. We report the first receptors that respond to Type 0 pheromone compounds; EsemOR3 displayed highest sensitivity toward (2S, 6Z)-6-nonen-2-ol, whereas EsemOR5 was most sensitive to the behavioral antagonist (Z)-6-nonen-2-one. These receptors also respond to plant volatiles of similar chemical structures, but with lower sensitivity. Phylogenetically, EsemOR3 and EsemOR5 group with a plant volatile-responding receptor from the tortricid moth Epiphyas postvittana (EposOR3), which together reside outside the previously defined lepidopteran PR clade that contains the PRs from more derived lepidopteran families. In addition, one receptor (EsemOR1) that falls at the base of the lepidopteran PR clade, responded specifically to β-caryophyllene and not to any other additional plant or pheromone compounds. Our results suggest that PRs for Type 0 pheromones have evolved from ORs that detect structurally-related plant volatiles. They are unrelated to PRs detecting pheromones in more derived Lepidoptera, which, in turn, also independently may have evolved a novel function from ORs detecting plant volatiles.

The genomic and functional landscapes of developmental plasticity in the American cockroach

Many cockroach species have adapted to urban environments, and some have been serious pests of public health in the tropics and subtropics. Here, we present the 3.38-Gb genome and a consensus gene set of the American cockroach, Periplaneta americana. We report insights from both genomic and functional investigations into the underlying basis of its adaptation to urban environments and developmental plasticity. In comparison with other insects, expansions of gene families in P. americana exist for most core gene families likely associated with environmental adaptation, such as chemoreception and detoxification. Multiple pathways regulating metamorphic development are well conserved, and RNAi experiments inform on key roles of 20-hydroxyecdysone, juvenile hormone, insulin, and decapentaplegic signals in regulating plasticity. Our analyses reveal a high level of sequence identity in genes between the American cockroach and two termite species, advancing it as a valuable model to study the evolutionary relationships between cockroaches and termites.

Comparative genomic analysis of SET domain family reveals the origin, expansion, and putative function of the arthropod-specific SmydA genes as histone modifiers in insects

The SET domain is an evolutionarily conserved motif present in histone lysine methyltransferases, which are important in the regulation of chromatin and gene expression in animals. In this study, we searched for SET domain–containing genes (SET genes) in all of the 147 arthropod genomes sequenced at the time of carrying out this experiment to understand the evolutionary history by which SET domains have evolved in insects. Phylogenetic and ancestral state reconstruction analysis revealed an arthropod-specific SET gene family, named SmydA, that is ancestral to arthropod animals and specifically diversified during insect evolution. Considering that pseudogenization is the most probable fate of the new emerging gene copies, we provided experimental and evolutionary evidence to demonstrate their essential functions. Fluorescence in situ hybridization analysis and in vitro methyltransferase activity assays showed that the SmydA-2 gene was transcriptionally active and retained the original histone methylation activity. Expression knockdown by RNA interference significantly increased mortality, implying that the SmydA genes may be essential for insect survival. We further showed predominantly strong purifying selection on the SmydA gene family and a potential association between the regulation of gene expression and insect phenotypic plasticity by transcriptome analysis. Overall, these data suggest that the SmydA gene family retains essential functions that may possibly define novel regulatory pathways in insects. This work provides insights into the roles of lineage-specific domain duplication in insect evolution.

Mass Spectrometry-Based Screening Platform Reveals Orco Interactome in Drosophila melanogaster

Animals use their odorant receptors to receive chemical information from the environment. Insect odorant receptors differ from the G protein-coupled odorant receptors in vertebrates and nematodes, and very little is known about their protein-protein interactions. Here, we introduce a mass spectrometric platform designed for the large-scale analysis of insect odorant receptor protein-protein interactions. Using this platform, we obtained the first Orco interactome from Drosophila melanogaster. From a total of 1,186 identified proteins, we narrowed the interaction candidates to 226, of which only two-thirds have been named. These candidates include the known olfactory proteins Or92a and Obp51a. Around 90% of the proteins having published names likely function inside the cell, and nearly half of these intracellular proteins are associated with the endomembrane system. In a basic loss-of-function electrophysiological screen, we found that the disruption of eight (i.e., Rab5, CG32795, Mpcp, Tom70, Vir-1, CG30427, Eaat1, and CG2781) of 28 randomly selected candidates affects olfactory responses in vivo. Thus, because this Orco interactome includes physiologically meaningful candidates, we anticipate that our platform will help guide further research on the molecular mechanisms of the insect odorant receptor family.

Endogenous insensitivity to the Orco agonist VUAA1 reveals novel olfactory receptor complex properties in the specialist fly Mayetiola destructor

Insect olfactory receptors are routinely expressed in heterologous systems for functional characterisation. It was recently discovered that the essential olfactory receptor co-receptor (Orco) of the Hessian fly, Mayetiola destructor (Mdes), does not respond to the agonist VUAA1, which activates Orco in all other insects analysed to date. Here, using a mutagenesis-based approach we identified three residues in MdesOrco, located in different transmembrane helices as supported by 3D modelling, that confer sensitivity to VUAA1. Reciprocal mutations in Drosophila melanogaster (Dmel) and the noctuid moth Agrotis segetum (Aseg) Orcos diminish sensitivity of these proteins to VUAA1. Additionally, mutating these residues in DmelOrco and AsegOrco compromised odourant receptor (OR) dependent ligand-induced Orco activation. In contrast, both wild-type and VUAA1-sensitive MdesOrco were capable of forming functional receptor complexes when coupled to ORs from all three species, suggesting unique complex properties in M. destructor, and that not all olfactory receptor complexes are “created” equal.

FUSTr: a tool to find gene families under selection in transcriptomes

Background

The recent proliferation of large amounts of biodiversity transcriptomic data has resulted in an ever-expanding need for scalable and user-friendly tools capable of answering large scale molecular evolution questions. FUSTr identifies gene families involved in the process of adaptation. This is a tool that finds genes in transcriptomic datasets under strong positive selection that automatically detects isoform designation patterns in transcriptome assemblies to maximize phylogenetic independence in downstream analysis.

Results

When applied to previously studied spider transcriptomic data as well as simulated data, FUSTr successfully grouped coding sequences into proper gene families as well as correctly identified those under strong positive selection in relatively little time.

Conclusions

FUSTr provides a useful tool for novice bioinformaticians to characterize the molecular evolution of organisms throughout the tree of life using large transcriptomic biodiversity datasets and can utilize multi-processor high-performance computational facilities.

Identification of odorant-binding and chemosensory protein genes and the ligand affinity of two of the encoded proteins suggest a complex olfactory perception system in Periplaneta americana

The American cockroach (Periplaneta americana) is an urban pest with a precise chemosensory system that helps it achieve complex physiological behaviours, including locating food and mating. However, its chemosensory mechanisms have not been well studied. Here, we identified 71 putative odorant carrier protein genes in P. americana, including 57 new odorant-binding proteins (OBPs) and 11 chemosensory proteins (CSPs). To identify their physiological functions, we investigated their tissue expression patterns in antennae, mouthparts, legs, and the remainder of the body of both sexes, and determined that most of these genes were expressed in chemosensory organs. A phylogenetic tree showed that the putative pheromone-binding proteins of P. americana were in different clades from those of moths. Two genes, PameOBP24 and PameCSP7, were expressed equally in antennae of both sexes and highly expressed amongst the OBPs and CSPs. These genes were expressed in Escherichia coli and the resultant proteins were purified. The binding affinities of 74 common odorant compounds were tested with recombinant PameOBP24 and PameCSP7. Both proteins bound a variety of ligands. Our findings provide a foundation for future research into the chemosensory mechanisms of P. americana and help in identifying potential target genes for managing this pest.

Current source density mapping of antennal sensory selectivity reveals conserved olfactory systems between tephritids and Drosophila

Ecological specialization of insects involves the functional and morphological reshaping of olfactory systems. Little is known about the degree to which insect sensitivity to odorant compounds is conserved between genera, tribes, or families. Here we compared the olfactory systems of six tephritid fruit fly species spanning two tribes and the distantly related Drosophila melanogaster at molecular, functional, and morphological levels. Olfaction in these flies is mediated by a set of olfactory receptors (ORs) expressed in different functional classes of neurons located in distinct antennal regions. We performed a phylogenetic analysis that revealed both family-specific OR genes and putative orthologous OR genes between tephritids and Drosophila. With respect to function, we then used a current source density (CSD) analysis to map activity across antennae. Functional maps mirrored the intrinsic structure of antennae observed with scanning electron microscopy. Together, the results revealed partial conservation of the olfactory systems between tephritids and Drosophila. We also demonstrate that the mapping of olfactory responses is necessary to decipher antennal sensory selectivity to olfactory compounds. CSD analysis can be easily applied to map antennae of other species and therefore enables the rapid deriving of olfactory maps and the reconstructing of the target organisms' history of evolution.

Genomic Features of the Damselfly Calopteryx splendens Representing a Sister Clade to Most Insect Orders

Insects comprise the most diverse and successful animal group with over one million described species that are found in almost every terrestrial and limnic habitat, with many being used as important models in genetics, ecology, and evolutionary research. Genome sequencing projects have greatly expanded the sampling of species from many insect orders, but genomic resources for species of certain insect lineages have remained relatively limited to date. To address this paucity, we sequenced the genome of the banded demoiselle, Calopteryx splendens, a damselfly (Odonata: Zygoptera) belonging to Palaeoptera, the clade containing the first winged insects. The 1.6 Gbp C. splendens draft genome assembly is one of the largest insect genomes sequenced to date and encodes a predicted set of 22,523 protein-coding genes. Comparative genomic analyses with other sequenced insects identified a relatively small repertoire of C. splendens detoxification genes, which could explain its previously noted sensitivity to habitat pollution. Intriguingly, this repertoire includes a cytochrome P450 gene not previously described in any insect genome. The C. splendens immune gene repertoire appears relatively complete and features several genes encoding novel multi-domain peptidoglycan recognition proteins. Analysis of chemosensory genes revealed the presence of both gustatory and ionotropic receptors, as well as the insect odorant receptor coreceptor gene (OrCo) and at least four partner odorant receptors (ORs). This represents the oldest known instance of a complete OrCo/OR system in insects, and provides the molecular underpinning for odonate olfaction. The C. splendens genome improves the sampling of insect lineages that diverged before the radiation of Holometabola and offers new opportunities for molecular-level evolutionary, ecological, and behavioral studies.

Getting Nervous: An Evolutionary Overhaul for Communication

The evolution of a nervous system as a control system of the body's functions is a key innovation of animals. Its fundamental units are neurons, highly specialized cells dedicated to fast cell-cell communication. Neurons pass signals to other neurons, muscle cells, or gland cells at specialized junctions, the synapses, where transmitters are released from vesicles in a Ca²⁺-dependent fashion to activate receptors in the membrane of the target cell. Reconstructing the origins of neuronal communication out of a more simple process remains a central challenge in biology. Recent genomic comparisons have revealed that all animals, including the nerveless poriferans and placozoans, share a basic set of genes for neuronal communication. This suggests that the first animal, the Urmetazoan, was already endowed with neurosecretory cells that probably started to connect into neuronal networks soon afterward. Here, we discuss scenarios for this pivotal transition in animal evolution.

Evolutionary History of Chemosensory-Related Gene Families across the Arthropoda

Chemosensory-related gene (CRG) families have been studied extensively in insects, but their evolutionary history across the Arthropoda had remained relatively unexplored. Here, we address current hypotheses and prior conclusions on CRG family evolution using a more comprehensive data set. In particular, odorant receptors were hypothesized to have proliferated during terrestrial colonization by insects (hexapods), but their association with other pancrustacean clades and with independent terrestrial colonizations in other arthropod subphyla have been unclear. We also examine hypotheses on which arthropod CRG family is most ancient. Thus, we reconstructed phylogenies of CRGs, including those from new arthropod genomes and transcriptomes, and mapped CRG gains and losses across arthropod lineages. Our analysis was strengthened by including crustaceans, especially copepods, which reside outside the hexapod/branchiopod clade within the subphylum Pancrustacea. We generated the first high-resolution genome sequence of the copepod Eurytemora affinis and annotated its CRGs. We found odorant receptors and odorant binding proteins present only in hexapods (insects) and absent from all other arthropod lineages, indicating that they are not universal adaptations to land. Gustatory receptors likely represent the oldest chemosensory receptors among CRGs, dating back to the Placozoa. We also clarified and confirmed the evolutionary history of antennal ionotropic receptors across the Arthropoda. All antennal ionotropic receptors in E. affinis were expressed more highly in males than in females, suggestive of an association with male mate-recognition behavior. This study is the most comprehensive comparative analysis to date of CRG family evolution across the largest and most speciose metazoan phylum Arthropoda.

Contact chemosensation of phytochemicals by insect herbivores

Contact chemosensation, or tasting, is a complex process governed by nonvolatile phytochemicals that tell host-seeking insects whether they should accept or reject a plant. During this process, insect gustatory receptors (GRs) contribute to deciphering a host plant's metabolic code. GRs recognise many different classes of nonvolatile compounds; some GRs are likely to be narrowly tuned and others, broadly tuned. Although primary and/or secondary plant metabolites influence the insect's feeding choice, their decoding by GRs is challenging, because metabolites in planta occur in complex mixtures that have additive or inhibitory effects; in diverse forms composed of structurally unrelated molecules; and at different concentrations depending on the plant species, its tissue and developmental stage. Future studies of the mechanism of insect herbivore GRs will benefit from functional characterisation taking into account the spatio-temporal dynamics and diversity of the plant's metabolome. Metabolic information, in turn, will help to elucidate the impact of single ligands and complex natural mixtures on the insect's feeding choice.

Antennal transcriptomes of three tortricid moths reveal putative conserved chemosensory receptors for social and habitat olfactory cues

Insects use chemical signals to find mates, food and oviposition sites. The main chemoreceptor gene families comprise odorant receptors (ORs), ionotropic receptors (IRs) and gustatory receptors (GRs). Understanding the evolution of these receptors as well as their function will assist in advancing our knowledge of how chemical stimuli are perceived and may consequently lead to the development of new insect management strategies. Tortricid moths are important pests in horticulture, forestry and agriculture around the globe. Here, we characterize chemoreceptors from the three main gene families of three economically important tortricids, based on male antennal transcriptomes using an RNA-Seq approach. We identified 49 ORs, 11 GRs and 23 IRs in the green budworm moth, Hedya nubiferana; 49 ORs, 12 GRs and 19 IRs in the beech moth, Cydia fagiglandana; and 48 ORs, 11 GRs and 19 IRs in the pea moth, Cydia nigricana. Transcript abundance estimation, phylogenetic relationships and molecular evolution rate comparisons with deorphanized receptors of Cydia pomonella allow us to hypothesize conserved functions and therefore candidate receptors for pheromones and kairomones.

Finding food: how marine invertebrates use chemical cues to track and select food

This review covers recent research on how marine invertebrates use chemical cues to find and select food.

Complex Evolutionary Dynamics of Massively Expanded Chemosensory Receptor Families in an Extreme Generalist Chelicerate Herbivore

While mechanisms to detoxify plant produced, anti-herbivore compounds have been associated with plant host use by herbivores, less is known about the role of chemosensory perception in their life histories. This is especially true for generalists, including chelicerate herbivores that evolved herbivory independently from the more studied insect lineages. To shed light on chemosensory perception in a generalist herbivore, we characterized the chemosensory receptors (CRs) of the chelicerate two-spotted spider mite, Tetranychus urticae, an extreme generalist. Strikingly, T. urticae has more CRs than reported in any other arthropod to date. Including pseudogenes, 689 gustatory receptors were identified, as were 136 degenerin/Epithelial Na+ Channels (ENaCs) that have also been implicated as CRs in insects. The genomic distribution of T. urticae gustatory receptors indicates recurring bursts of lineage-specific proliferations, with the extent of receptor clusters reminiscent of those observed in the CR-rich genomes of vertebrates or C. elegans Although pseudogenization of many gustatory receptors within clusters suggests relaxed selection, a subset of receptors is expressed. Consistent with functions as CRs, the genomic distribution and expression of ENaCs in lineage-specific T. urticae expansions mirrors that observed for gustatory receptors. The expansion of ENaCs in T. urticae to > 3-fold that reported in other animals was unexpected, raising the possibility that ENaCs in T. urticae have been co-opted to fulfill a major role performed by unrelated CRs in other animals. More broadly, our findings suggest an elaborate role for chemosensory perception in generalist herbivores that are of key ecological and agricultural importance.

Positive Allosteric Modulation of Insect Olfactory Receptor Function by ORco Agonists

Insect olfactory receptors (ORs) are heteromeric ligand-gated cation channels composed of a common olfactory receptor subunit (ORco) and a variable subunit (ORx) of as yet unknown structures and undetermined stoichiometries. In this study, we examined the allosteric modulation exerted on Anopheles gambiae heteromeric ORx/ORco olfactory receptors in vitro by a specific class of ORco agonists (OAs) comprising ORcoRAM2 and VUAA1. High OA concentrations produced stronger functional responses in cells expressing heteromeric receptor channels relative to cells expressing ORco alone. These OA-induced responses of ORx/ORco channels were also notably much stronger than those obtained upon administration of ORx-specific ligands to the same receptors. Most importantly, small concentrations of OAs were found to act as strong potentiators of ORx/ORco function, increasing dramatically both the efficacy and potency of ORx-specific odorants. These results suggest that insect heteromeric ORs are highly dynamic complexes adopting different conformations that change in a concerted fashion as a result of the interplay between the subunits of the oligomeric assemblies, and that allosteric modulation may constitute an important element in the modulation and fining tuning of olfactory reception function.

Under-Expression of Chemosensory Genes in Domiciliary Bugs of the Chagas Disease Vector Triatoma brasiliensis

BACKGROUND

In Latin America, the bloodsucking bugs Triatominae are vectors of Trypanosoma cruzi, the parasite that causes Chagas disease. Chemical elimination programs have been launched to control Chagas disease vectors. However, the disease persists because native vectors from sylvatic habitats are able to (re)colonize houses-a process called domiciliation. Triatoma brasiliensis is one example. Because the chemosensory system allows insects to interact with their environment and plays a key role in insect adaption, we conducted a descriptive and comparative study of the chemosensory transcriptome of T. brasiliensis samples from different ecotopes.

METHODOLOGY/PRINCIPAL FINDING

In a reference transcriptome built using de novo assembly, we found transcripts encoding 27 odorant-binding proteins (OBPs), 17 chemosensory proteins (CSPs), 3 odorant receptors (ORs), 5 transient receptor potential channel (TRPs), 1 sensory neuron membrane protein (SNMPs), 25 takeout proteins, 72 cytochrome P450s, 5 gluthatione S-transferases, and 49 cuticular proteins. Using protein phylogenies, we showed that most of the OBPs and CSPs for T. brasiliensis had well supported orthologs in the kissing bug Rhodnius prolixus. We also showed a higher number of these genes within the bloodsucking bugs and more generally within all Hemipterans compared to the other species in the super-order Paraneoptera. Using both DESeq2 and EdgeR software, we performed differential expression analyses between samples of T. brasiliensis, taking into account their environment (sylvatic, peridomiciliary and domiciliary) and sex. We also searched clusters of co-expressed contigs using HTSCluster. Among differentially expressed (DE) contigs, most were under-expressed in the chemosensory organs of the domiciliary bugs compared to the other samples and in females compared to males. We clearly identified DE genes that play a role in the chemosensory system.

CONCLUSION/SIGNIFICANCE

Chemosensory genes could be good candidates for genes that contribute to adaptation or plastic rearrangement to an anthropogenic system. The domiciliary environment probably includes less diversity of xenobiotics and probably has more stable abiotic parameters than do sylvatic and peridomiciliary environments. This could explain why both detoxification and cuticle protein genes are less expressed in domiciliary bugs. Understanding the molecular basis for how vectors adapt to human dwellings may reveal new tools to control disease vectors; for example, by disrupting chemical communication.

A Sex Pheromone Receptor in the Hessian Fly Mayetiola destructor (Diptera, Cecidomyiidae)

The Hessian fly, Mayetiola destructor Say (Diptera, Cecidomyiidae), is a pest of wheat and belongs to a group of gall-inducing herbivores. This species has a unique life history and several ecological features that differentiate it from other Diptera such as Drosophila melanogaster and blood-feeding mosquitoes. These features include a short, non-feeding adult life stage (1-2 days) and the use of a long-range sex pheromone produced and released by adult females. Sex pheromones are detected by members of the odorant receptor (OR) family within the Lepidoptera, but no receptors for similar long-range sex pheromones have been characterized from the Diptera. Previously, 122 OR genes have been annotated from the Hessian fly genome, with many of them showing sex-biased expression in the antennae. Here we have expressed, in HEK293 cells, five MdesORs that display male-biased expression in antennae, and we have identified MdesOR115 as a Hessian fly sex pheromone receptor. MdesOR115 responds primarily to the sex pheromone component (2S,8E,10E)-8,10-tridecadien-2-yl acetate, and secondarily to the corresponding Z,E-isomer. Certain sensory neuron membrane proteins (i.e., SNMP1) are important for responses of pheromone receptors in flies and moths. The Hessian fly genome is unusual in that it encodes six SNMP1 paralogs, of which five are expressed in antennae. We co-expressed each of the five antennal SNMP1 paralogs together with each of the five candidate sex pheromone receptors from the Hessian fly and found that they do not influence the response of MdesOR115, nor do they confer responsiveness in any of the non-responsive ORs to any of the sex pheromone components identified to date in the Hessian fly. Using Western blots, we detected protein expression of MdesOrco, all MdesSNMPs, and all MdesORs except for MdesOR113, potentially explaining the lack of response from this OR. In conclusion, we report the first functional characterization of an OR from the Cecidomyiidae, extending the role of ORs as long-range sex pheromone detectors from the Lepidoptera into the Diptera.

Multifaceted biological insights from a draft genome sequence of the tobacco hornworm moth, Manduca sexta

Manduca sexta, known as the tobacco hornworm or Carolina sphinx moth, is a lepidopteran insect that is used extensively as a model system for research in insect biochemistry, physiology, neurobiology, development, and immunity. One important benefit of this species as an experimental model is its extremely large size, reaching more than 10 g in the larval stage. M. sexta larvae feed on solanaceous plants and thus must tolerate a substantial challenge from plant allelochemicals, including nicotine. We report the sequence and annotation of the M. sexta genome, and a survey of gene expression in various tissues and developmental stages. The Msex_1.0 genome assembly resulted in a total genome size of 419.4 Mbp. Repetitive sequences accounted for 25.8% of the assembled genome. The official gene set is comprised of 15,451 protein-coding genes, of which 2498 were manually curated. Extensive RNA-seq data from many tissues and developmental stages were used to improve gene models and for insights into gene expression patterns. Genome wide synteny analysis indicated a high level of macrosynteny in the Lepidoptera. Annotation and analyses were carried out for gene families involved in a wide spectrum of biological processes, including apoptosis, vacuole sorting, growth and development, structures of exoskeleton, egg shells, and muscle, vision, chemosensation, ion channels, signal transduction, neuropeptide signaling, neurotransmitter synthesis and transport, nicotine tolerance, lipid metabolism, and immunity. This genome sequence, annotation, and analysis provide an important new resource from a well-studied model insect species and will facilitate further biochemical and mechanistic experimental studies of many biological systems in insects.

Antennal transcriptome analysis and expression profiles of odorant binding proteins in Eogystia hippophaecolus (Lepidoptera: Cossidae)

Background

Eogystia hippophaecolus (Hua et al.) (Lepidoptera: Cossidae) is the major threat to seabuckthorn plantations in China. Specific and highly efficient artificial sex pheromone traps was developed and used to control it. However, the molecular basis for the pheromone recognition is not known. So we established the antennal transcriptome of E. hippophaecolus and characterized the expression profiles of odorant binding proteins. These results establish and improve the basis knowledge of the olfactory receptive system, furthermore provide a theoretical basis for the development of new pest control method.

Results

We identified 29 transcripts encoding putative odorant-binding proteins (OBPs), 18 putative chemosensory proteins (CSPs), 63 odorant receptors (ORs), 13 gustatory receptors (GRs), 12 ionotropic receptors (IRs), and two sensory neuron membrane proteins (SNMPs). Based on phylogenetic analysis, we found one Orco and three pheromone receptors of E. hippophaecolus and found that EhipGR13 detects sugar, EhipGR11 and EhipGR3 detect bitter. Nine OBPs expression profile indicated that most were the highest expression in antennae, consistent with functions of OBPs in binding and transporting odors during the antennal recognition process. OBP6 was external expressed in male genital-biased in, and this locus may be responsible for pheromone binding and recognition as well as mating. OBP1 was the highest and biased expressed in the foot and may function as identification of host plant volatiles.

Conclusions

One hundred thirty-seven chemosensory proteins were identified and the accurate functions and groups of part proteins were obtained by phylogenetic analysis. The most OBPs were antenna-biased expressed, which are involved in antennal recognition. However, few OBP was detected biased expression in the foot and external genitalia, and these loci may function in pheromone recognition, mating, and the recognition of plant volatiles.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3008-4) contains supplementary material, which is available to authorized users.

Identification and tissue expression profile of genes from three chemoreceptor families in an urban pest, Periplaneta americana

Periplaneta americana is a notorious urban pest prevalent in human habitats; very little is known about its chemosensory mechanism. Employing the advanced next-generation sequencing technique, in the present study, we conducted transcriptome sequencing and analysis of the antennae of the adult males and females as well as their mouthparts using an Illumina platform. This resulted in the discovery of a huge number of the members of all major known chemosensory receptor families in P. americana, including 96 odorant receptors (ORs), 53 ionotropic receptors (IRs), and 33 gustatory receptors (GRs). Tissue expression profiles showed most of them mainly expressed in antennae and phylogenetic analysis demonstrated the expansion in the clade distinguishing them from other functionally well-known Lepidoptera species. A high percentage of chemosensory receptor genes (ORs in particular) showing female antenna bias in mRNA expression was observed. Our results provide a basis for further investigations on how P. americana coordinates its chemosensory receptor genes in chemical communication with environments, and for development of novel pest management approaches.

Molecular Mechanisms of Reception and Perireception in Crustacean Chemoreception: A Comparative Review

This review summarizes our present knowledge of chemoreceptor proteins in crustaceans, using a comparative perspective to review these molecules in crustaceans relative to other metazoan models of chemoreception including mammals, insects, nematodes, and molluscs. Evolution has resulted in unique expansions of specific gene families and repurposing of them for chemosensation in various clades, including crustaceans. A major class of chemoreceptor proteins across crustaceans is the Ionotropic Receptors, which diversified from ionotropic glutamate receptors in ancient protostomes but which are not present in deuterostomes. Representatives of another major class of chemoreceptor proteins-the Grl/GR/OR family of ionotropic 7-transmembrane receptors-are diversified in insects but to date have been reported in only one crustacean species, Daphnia pulex So far, canonic 7-transmembrane G-protein coupled receptors, the principal chemoreceptors in vertebrates and reported in a few protostome clades, have not been identified in crustaceans. More types of chemoreceptors are known throughout the metazoans and might well be expected to be discovered in crustaceans. Our review also provides a comparative coverage of perireceptor events in crustacean chemoreception, including molecules involved in stimulus acquisition, stimulus delivery, and stimulus removal, though much less is known about these events in crustaceans, particularly at the molecular level.