Expansion of the fatty acyl reductase gene family shaped pheromone communication in Hymenoptera

Post-transcriptional regulation of Dufour's gland reproductive signals in bumble bees

Pheromone communication is a key mechanism by which the reproductive division of labor is maintained within insect communities. Understanding how pheromones evolved to regulate social behavior requires knowledge of the molecular regulation of their production. However, even in cases where pheromones were identified, our understanding of their biosynthesis and molecular regulation remains limited. Bumble bees provide a unique system to explore pheromone biosynthesis since workers produce ester sterility signals in their Dufour's gland that differ from gyne-specific esters and are not produced by queens. These esters are hypothesized to be produced in the exocrine gland where they are stored, and indeed queens, gynes and workers differ significantly in the expression of Dufour's gland genes coding to enzymes involved in the biosynthesis of esters. However, a previous transcriptome analysis revealed no gene expression differences in the Dufour's gland of workers despite differences in both ester production and ovarian activation, suggesting that ester production may be regulated lower down. Proteomics of the Dufour's gland of queens, gynes, and workers recovered over 2400 proteins and broadly matched the previous RNAseq data. However, more than 100 differentially expressed proteins were found between the worker groups, including key enzymes in fatty acid biosynthesis, indicating that the regulation of reproductive signal biosynthesis in workers is done post-transcription. Overall, our data provide evidence that pheromone biosynthesis in the Dufour's gland is caste specific, that gynes and workers are likely using different enzymes to make their respective wax esters, and that the regulation on pheromone production in queens, gynes and workers is likely done at different regulatory levels, with workers signals being subjected to regulation at the protein level.

Causes and consequences of a complex recombinational landscape in the ant Cardiocondyla obscurior

Eusocial Hymenoptera have the highest recombination rates among all multicellular animals studied so far, but it is unclear why this is and how this affects the biology of individual species. A high-resolution linkage map for the ant Cardiocondyla obscurior corroborates genome-wide high recombination rates reported for ants (8.1 cM/Mb). However, recombination is locally suppressed in regions that are enriched with TEs, that have strong haplotype divergence, or that show signatures of epistatic selection in C. obscurior The results do not support the hypotheses that high recombination rates are linked to phenotypic plasticity or to modulating selection efficiency. Instead, genetic diversity and the frequency of structural variants correlate positively with local recombination rates, potentially compensating for the low levels of genetic variation expected in haplodiploid social Hymenoptera with low effective population size. Ultimately, the data show that recombination contributes to within-population polymorphism and to the divergence of the lineages within C. obscurior.

The genomic and cellular basis of biosynthetic innovation in rove beetles

How evolution at the cellular level potentiates macroevolutionary change is central to understanding biological diversification. The >66,000 rove beetle species (Staphylinidae) form the largest metazoan family. Combining genomic and cell type transcriptomic insights spanning the largest clade, Aleocharinae, we retrace evolution of two cell types comprising a defensive gland-a putative catalyst behind staphylinid megadiversity. We identify molecular evolutionary steps leading to benzoquinone production by one cell type via a mechanism convergent with plant toxin release systems, and synthesis by the second cell type of a solvent that weaponizes the total secretion. This cooperative system has been conserved since the Early Cretaceous as Aleocharinae radiated into tens of thousands of lineages. Reprogramming each cell type yielded biochemical novelties enabling ecological specialization-most dramatically in symbionts that infiltrate social insect colonies via host-manipulating secretions. Our findings uncover cell type evolutionary processes underlying the origin and evolvability of a beetle chemical innovation.

'Junk' that matters: the role of transposable elements in bumblebee genome evolution

Transposable elements (TEs) are mobile DNA sequences that are widely distributed in eukaryotic genomes, where they are known to serve as a major force in genome evolution. The phenotypic impacts of TEs, while less well-studied, have also been discovered. Bumblebees are globally important pollinators in natural ecosystems and agriculture. Although TEs comprise a small fraction of bumblebee genomes, emerging evidence suggests that TEs are the major contributor of genome size variation across species and are involved in the formation of new coding and regulatory sequences. We review recent discoveries related to TEs in bumblebees, as well as outlining three key questions for the future of the field. In the future, we argue long-read sequencing technologies and genome editing techniques will help us identify TEs in bumblebees, unveil mechanisms that could account for their silencing and limited abundance, and uncover their contributions to phenotypic diversification, ecological adaptation, and speciation.

Evolution of a fatty acyl-CoA elongase underlies desert adaptation in Drosophila

Traits that allow species to survive in extreme environments such as hot-arid deserts have independently evolved in multiple taxa. However, the genetic and evolutionary mechanisms underlying these traits have thus far not been elucidated. Here, we show that Drosophila mojavensis, a desert-adapted fruit fly species, has evolved high desiccation resistance by producing long-chain methyl-branched cuticular hydrocarbons (mbCHCs) that contribute to a cuticular lipid layer reducing water loss. We show that the ability to synthesize these longer mbCHCs is due to evolutionary changes in a fatty acyl-CoA elongase (mElo). mElo knockout in D. mojavensis led to loss of longer mbCHCs and reduction of desiccation resistance at high temperatures but did not affect mortality at either high temperatures or desiccating conditions individually. Phylogenetic analysis showed that mElo is a Drosophila-specific gene, suggesting that while the physiological mechanisms underlying desert adaptation may be similar between species, the genes involved in these mechanisms may be species or lineage specific.

A novel domain-duplicated SlitFAR3 gene involved in sex pheromone biosynthesis in Spodoptera litura

Fatty acyl reductases (FARs) are key enzymes that participate in sex pheromone biosynthesis by reducing fatty acids to fatty alcohols. Lepidoptera typically harbor numerous FAR gene family members. Although FAR genes are involved in the biosynthesis of sex pheromones in moths, the key FAR gene of Spodoptera litura remains unclear. In this work, we predicted 30 FAR genes from the S. litura genome and identified a domain duplication within gene SlitFAR3, which exhibited high and preferential expression in the sexually mature female pheromone glands (PGs) and a rhythmic expression pattern during the scotophase of sex pheromone production. The molecular docking of SlitFAR3, as predicted using a 3D model, revealed a co-factor NADPH binding cavity and 2 substrate binding cavities. Functional expression in yeast cells combined with comprehensive gas chromatography indicated that the SlitFAR3 gene could produce fatty alcohol products. This study is the first to focus on the special phenomenon of FAR domain duplication, which will advance our understanding of biosynthesis-related genes from the perspective of evolutionary biology.

The genomic and cellular basis of biosynthetic innovation in rove beetles

How evolution at the cellular level potentiates change at the macroevolutionary level is a major question in evolutionary biology. With >66,000 described species, rove beetles (Staphylinidae) comprise the largest metazoan family. Their exceptional radiation has been coupled to pervasive biosynthetic innovation whereby numerous lineages bear defensive glands with diverse chemistries. Here, we combine comparative genomic and single-cell transcriptomic data from across the largest rove beetle clade, Aleocharinae. We retrace the functional evolution of two novel secretory cell types that together comprise the tergal gland-a putative catalyst behind Aleocharinae's megadiversity. We identify key genomic contingencies that were critical to the assembly of each cell type and their organ-level partnership in manufacturing the beetle's defensive secretion. This process hinged on evolving a mechanism for regulated production of noxious benzoquinones that appears convergent with plant toxin release systems, and synthesis of an effective benzoquinone solvent that weaponized the total secretion. We show that this cooperative biosynthetic system arose at the Jurassic-Cretaceous boundary, and that following its establishment, both cell types underwent ∼150 million years of stasis, their chemistry and core molecular architecture maintained almost clade-wide as Aleocharinae radiated globally into tens of thousands of lineages. Despite this deep conservation, we show that the two cell types have acted as substrates for the emergence of adaptive, biochemical novelties-most dramatically in symbiotic lineages that have infiltrated social insect colonies and produce host behavior-manipulating secretions. Our findings uncover genomic and cell type evolutionary processes underlying the origin, functional conservation and evolvability of a chemical innovation in beetles.

Intrasexual cuticular hydrocarbon dimorphism in a wasp sheds light on hydrocarbon biosynthesis genes in Hymenoptera

Cuticular hydrocarbons (CHCs) cover the cuticle of insects and serve as desiccation barrier and as semiochemicals. While the main enzymatic steps of CHC biosynthesis are well understood, few of the underlying genes have been identified. Here we show how exploitation of intrasexual CHC dimorphism in a mason wasp, Odynerus spinipes, in combination with whole-genome sequencing and comparative transcriptomics facilitated identification of such genes. RNAi-mediated knockdown of twelve candidate gene orthologs in the honey bee, Apis mellifera, confirmed nine genes impacting CHC profile composition. Most of them have predicted functions consistent with current knowledge of CHC metabolism. However, we found first-time evidence for a fatty acid amide hydrolase also influencing CHC profile composition. In situ hybridization experiments furthermore suggest trophocytes participating in CHC biosynthesis. Our results set the base for experimental CHC profile manipulation in Hymenoptera and imply that the evolutionary origin of CHC biosynthesis predates the arthropods' colonization of land.

Ancient origin and conserved gene function in terpene pheromone and defense evolution of stink bugs and hemipteran insects

Insects use diverse arrays of small molecules such as metabolites of the large class of terpenes for intra- and inter-specific communication and defense. These molecules are synthesized by specialized metabolic pathways; however, the origin of enzymes involved in terpene biosynthesis and their evolution in insect genomes is still poorly understood. We addressed this question by investigating the evolution of isoprenyl diphosphate synthase (IDS)-like genes with terpene synthase (TPS) function in the family of stink bugs (Pentatomidae) within the large order of piercing-sucking Hemipteran insects. Stink bugs include species of global pest status, many of which emit structurally related 15-carbon sesquiterpenes as sex or aggregation pheromones. We provide evidence for the emergence of IDS-type TPS enzymes at the onset of pentatomid evolution over 100 million years ago, coinciding with the evolution of flowering plants. Stink bugs of different geographical origin maintain small IDS-type families with genes of conserved TPS function, which stands in contrast to the diversification of TPS genes in plants. Expanded gene mining and phylogenetic analysis in other hemipteran insects further provides evidence for an ancient emergence of IDS-like genes under presumed selection for terpene-mediated chemical interactions, and this process occurred independently from a similar evolution of IDS-type TPS genes in beetles. Our findings further suggest differences in TPS diversification in insects and plants in conjunction with different modes of gene functionalization in chemical interactions.

Differential gene expression underlying the biosynthesis of Dufour's gland signals in Bombus impatiens

Pheromones regulating social behavior are one of the most explored phenomena in social insects. However, compound identity, biosynthesis and their genetic basis are known in only a handful of species. Here we examined the gene expression associated with pheromone biosynthesis of two main chemical classes: esters and terpenes, using the social bee Bombus impatiens. We conducted chemical and RNA-seq analyses of the Dufour's gland, an exocrine gland producing a plethora of pheromones regulating social behavior in hymenopteran species. The Dufour's gland contains mostly long-chained hydrocarbons, terpenes and esters that signal reproductive and social status in several bee species. In bumble bees, the Dufour's gland contains queen- and worker-specific esters, in addition to terpenes and terpene-esters only found in gynes and queens. These compounds are assumed to be synthesized de novo in the gland, however, their genetic basis is unknown. A whole transcriptome gene expression analysis of the gland in queens, gynes, queenless and queenright workers showed distinct transcriptomic profiles, with thousands of differentially expressed genes between the groups. Workers and queens express genes associated with key enzymes in the biosynthesis of wax esters, while queens and gynes preferentially express key genes in terpene biosynthesis. Overall, our data demonstrate gland-specific regulation of chemical signals associated with social behavior and identifies candidate genes and pathways regulating caste-specific chemical signals in social insects.

Genome-wide identification of FAR gene family and functional analysis of NlFAR10 during embryogenesis in the brown planthopper Nilaparvata lugens

Fatty acyl-CoA reductases (FARs) catalyze the synthesis of fatty alcohols from corresponding fatty acid precursors in organisms. However, the function of FARs in insect fecundity and embryogenesis remains largely unclear. Here, a total of 22 putative FAR proteins were identified in the brown planthopper Nilaparvata lugens, a hemipteran insect pest of rice, and most of them were highly expressed in embryonic stages. Among them, NlFAR10 was specifically and highly expressed in the later embryogenesis, but was promiscuously expressed in tissues of adults. The heterologously expressed NlFAR10 was able to produce the intermediate fatty acid alcohols from the corresponding acyl-CoA precursors. When NlFAR10 was silenced through RNAi in vivo, the embryogenesis was obviously inhibited, resulting in low hatching rates. Moreover, the metabolome analyses indicated that loss of NlFAR10 affected lipid metabolism and purine metabolism during embryogenesis. To the best of our knowledge, this is the first report of a FAR member affecting insect embryogenesis, thus providing a new target for future pest management.

Mosaic Evolution of Molecular Pathways for Sex Pheromone Communication in a Butterfly

Unraveling the origin of molecular pathways underlying the evolution of adaptive traits is essential for understanding how new lineages emerge, including the relative contribution of conserved ancestral traits and newly evolved derived traits. Here, we investigated the evolutionary divergence of sex pheromone communication from moths (mostly nocturnal) to butterflies (mostly diurnal) that occurred ~119 million years ago. In moths, it is the females that typically emit pheromones to attract male mates, but in butterflies males emit pheromones that are used by females for mate choice. The molecular bases of sex pheromone communication are well understood in moths, but they have remained relatively unexplored in butterflies. We used a combination of transcriptomics, real time qPCR, and phylogenetics to identify genes involved in the different steps (i.e., production, regulation, and reception) of sex pheromone communication of the butterfly Bicyclus anynana. Our results show that the biosynthesis and reception of sex pheromones relies both on moth-specific gene families (reductases) and on more ancestral insect gene families (desaturases, olfactory receptors, odorant binding proteins). Interestingly, B. anynana appears to use what was believed to be the moth-specific neuropeptide Pheromone Biosynthesis Activating Neuropeptide (PBAN) for regulating sex pheromone production. Altogether, our results suggest that a mosaic pattern best explains how sex pheromone communication evolved in butterflies, with some molecular components derived from moths, and others conserved from more ancient insect ancestors. This is the first large-scale investigation of the genetic pathways underlying sex pheromone communication in a butterfly.

Developments in Fatty Acid-Derived Insect Pheromone Production Using Engineered Yeasts

The use of traditional chemical insecticides for pest control often leads to environmental pollution and a decrease in biodiversity. Recently, insect sex pheromones were applied for sustainable biocontrol of pests in fields, due to their limited adverse impacts on biodiversity and food safety compared to that of other conventional insecticides. However, the structures of insect pheromones are complex, and their chemical synthesis is not commercially feasible. As yeasts have been widely used for fatty acid-derived pheromone production in the past few years, using engineered yeasts may be promising and sustainable for the low-cost production of fatty acid-derived pheromones. The primary fatty acids produced by Saccharomyces cerevisiae and other yeasts are C16 and C18, and it is also possible to rewire/reprogram the metabolic flux for other fatty acids or fatty acid derivatives. This review summarizes the fatty acid biosynthetic pathway in S. cerevisiae and recent progress in yeast engineering in terms of metabolic engineering and synthetic biology strategies to produce insect pheromones. In the future, insect pheromones produced by yeasts might provide an eco-friendly pest control method in agricultural fields.

Independent variations in genome-wide expression, alternative splicing, and DNA methylation in brain tissues among castes of the buff-tailed bumblebee, Bombus terrestris

Caste differentiation in social hymenopterans is an intriguing example of phenotypic plasticity. However, the co-ordination among gene regulatory factors to mediate caste differentiation remains inconclusive. In this study, we determined the role of gene regulation and related epigenetic processes in pre-imaginal caste differentiation in the primitively eusocial bumblebee Bombus terrestris. By combining RNA-Seq data from Illumina and PacBio and accurately quantifying methylation at whole-genomic base pair resolution, we found that queens, workers, and drones mainly differentiate in gene expression but not in alternative splicing and DNA methylation. Gynes are the most distinct with the lowest global level of whole-genomic methylation and with the largest number of caste-specific transcripts and alternative splicing events. By contrast, workers exhibit few uniquely expressed or alternatively spliced genes. Moreover, several genes involved in hormone and neurotransmitter metabolism are related to caste differentiation, whereas several neuropeptides are linked with sex differentiation. Despite little genome-wide association among differential gene expression, splicing, and differential DNA methylation, the overlapped gene ontology (GO) terms point to nutrition-related activity. Therefore, variations in gene regulation correlate with the behavioral differences among castes and highlight the specialization of toolkit genes in bumblebee gynes at the beginning of the adult stage.

Integrated miRNA and transcriptome profiling to explore the molecular determinism of convergent adaptation to corn in two lepidopteran pests of agriculture

Background

The degree to which adaptation to same environment is determined by similar molecular mechanisms, is a topic of broad interest in evolutionary biology, as an indicator of evolutionary predictability. We wished to address if adaptation to the same host plant in phytophagous insects involved related gene expression patterns. We compared sRNA-Seq and RNA-Seq data between two pairs of taxa of Ostrinia and Spodoptera frugiperda sharing maize as host-plant. For the latter, we had previously carried out a reciprocal transplant experiment by feeding of the larvae of the Corn strain (Sf-C) and the Rice strain (Sf-R) on corn versus rice and characterized the mRNA and miRNA responses.

Results

First, we predicted the genes encoding miRNA in Ostrinia nubilalis (On) and O. scapulalis (Os). Respectively 67 and 65 known miRNA genes, as well as 196 and 190 novel ones were predicted with Os genome using sncRNAs extracted from whole larvae feeding on corn or mugwort. In On, a read counts analysis showed that 37 (55.22%) known miRNAs and 19 (9.84%) novel miRNAs were differentially expressed (DE) on mugwort compared to corn (in Os, 25 known miRs (38.46%) and 8 novel ones (4.34%)). Between species on corn, 8 (12.5%) known miRNAs and 8 (6.83%) novel ones were DE while only one novel miRNA showed expression variation between species on mugwort. Gene target prediction led to the identification of 2953 unique target genes in On and 2719 in Os, among which 11.6% (344) were DE when comparing species on corn. 1.8% (54) of On miR targets showed expression variation upon a change of host-plant.

We found molecular changes matching convergent phenotype, i.e., a set of nine miRNAs that are regulated either according to the host-plant both in On and Sf-C or between them on the same plant, corn. Among DE miR target genes between taxa, 13.7% shared exactly the same annotation between the two pairs of taxa and had function related to insect host-plant interaction.

Conclusion

There is some similarity in underlying genetic mechanisms of convergent evolution of two distant Lepidopteran species having adopted corn in their host range, highlighting possible adaptation genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07905-7.

Candidate genes involved in cuticular hydrocarbon differentiation between cryptic, parabiotic ant species

Insect cuticular hydrocarbons (CHCs) are highly diverse and have multiple functions, including communication and waterproofing. CHC profiles form species-specific, complex blends of up to 150 compounds. Especially in ants, even closely related species can have largely different profiles, raising the question how CHC differences are mirrored in the regulation of biosynthetic pathways. The neotropical ants Crematogaster levior and Camponotus femoratus both consist of two cryptic species each that are morphologically similar, but express strongly different CHC profiles. This is ideal to study the molecular basis of CHC differences. We thus investigated gene expression differences in fat-body transcriptomes of these ants. Despite common garden conditions, we found several thousand differentially expressed transcripts within each cryptic species pair. Many of these were related to metabolic processes, probably accounting for physiological differences. Moreover, we identified candidate genes from five gene families involved in CHC biosynthesis. By assigning candidate transcripts to orthologs in Drosophila, we inferred which CHCs might be influenced by differential gene expression. Expression of these candidate genes was often mirrored in the CHC profiles. For example, Cr. levior A, which has longer CHCs than its cryptic sister species, had a higher expression of elongases and a lower expression of fatty acyl- CoA reductases. This study is one of the first to identify CHC candidate genes in ants and will provide a basis for further research on the genetic basis of CHC biosynthesis.

Cell type innovation at the tips of the animal tree

Understanding how organs originate is challenging due to the twin problems of explaining how new cell types evolve and how collective interactions between cell types arise and become selectively advantageous. Animals are assemblages of organs and cell types of different antiquities, and among the most rapidly and convergently evolving are exocrine glands and their constituent secretory cell types. Such structures have arisen independently thousands of times across the Metazoa, impacting how animals chemically interact with their environments. The recurrent evolution of exocrine systems provides a paradigm for examining how qualitative phenotypic novelties arise from variation at the cellular level. Here, we take a hierarchical perspective, focusing on the evolutionary assembly of novel biosynthetic pathways and secretory cell types, and how both selection and non-adaptive molecular processes may combine to build the complex, modular architectures of many animal glands.

Production of Long Chain Fatty Alcohols Found in Bumblebee Pheromones by Yarrowia lipolytica

Fatty alcohols (FA-OH) are aliphatic unbranched primary alcohols with a chain of four or more carbon atoms. Besides potential industrial applications, fatty alcohols have important biological functions as well. In nature, fatty alcohols are produced as a part of a mixture of pheromones in several insect species, such as moths, termites, bees, wasps, etc. In addition, FA-OHs have a potential for agricultural applications, for example, they may be used as a suitable substitute for commercial insecticides. The insecticides have several drawbacks associated with their preparation, and they exert a negative impact on the environment. Currently, pheromone components are prepared mainly through the catalytic hydrogenation of plant oils and petrochemicals, which is an unsustainable, ecologically unfriendly, and highly expensive process. The biotechnological production of the pheromone components using engineered microbial strains and through the expression of the enzymes participating in the biosynthesis of these components is a promising approach that ensures ecological sustenance as well. The present study was aimed at evaluating the production of FA-OHs in the oleaginous yeast, Yarrowia lipolytica, with different lengths of fatty-acyl chains by expressing the fatty acyl-CoA reductase (FAR) BlapFAR4 from B. lapidarius, producing C16:0-OH, C16:1Δ⁹-OH, and lower quantities of both C14:0-OH and C18:1Δ⁹-OH, and BlucFAR1 from B. lucorum, producing FA-OHs with a chain length of 18-26 carbon atoms, in this yeast. Among the different novel Y. lipolytica strains used in the present study, the best results were obtained with JMY7086, which carried several lipid metabolism modifications and expressed the BlucFAR1 gene under the control of a strong constitutive promoter 8UAS-pTEF. JMY7086 produced only saturated fatty alcohols with chain lengths from 18 to 24 carbon atoms. The highest titer and accumulation achieved were 166.6 mg/L and 15.6 mg/g DCW of fatty alcohols, respectively. Unlike JMY7086, the BlapFAR4-expressing strain JMY7090 produced only 16 carbon atom-long FA-OHs with a titer of 14.6 mg/L.

The evolution of species recognition labels in insects

The evolution of pre-zygotic reproductive isolation is a key step in the process of speciation. In many organisms, particularly insects, chemical labels are used as pheromones for species-specific mate recognition. Although an enormous body of knowledge exists regarding the patterns of pheromone chemical ecology, much less is known about the evolutionary processes that underlie the origin of new mating pheromones. Here, we examine case studies that have illuminated the origins of species-specific mating pheromones and suggest future directions for productive research. This article is part of the theme issue 'Signal detection theory in recognition systems: from evolving models to experimental tests'.

Phylogenomics Identifies an Ancestral Burst of Gene Duplications Predating the Diversification of Aphidomorpha

Aphids (Aphidoidea) are a diverse group of hemipteran insects that feed on plant phloem sap. A common finding in studies of aphid genomes is the presence of a large number of duplicated genes. However, when these duplications occurred remains unclear, partly due to the high relatedness of sequenced species. To better understand the origin of aphid duplications we sequenced and assembled the genome of Cinara cedri, an early branching lineage (Lachninae) of the Aphididae family. We performed a phylogenomic comparison of this genome with 20 other sequenced genomes, including the available genomes of five other aphids, along with the transcriptomes of two species belonging to Adelgidae (a closely related clade to the aphids) and Coccoidea. We found that gene duplication has been pervasive throughout the evolution of aphids, including many parallel waves of recent, species-specific duplications. Most notably, we identified a consistent set of very ancestral duplications, originating from a large-scale gene duplication predating the diversification of Aphidomorpha (comprising aphids, phylloxerids, and adelgids). Genes duplicated in this ancestral wave are enriched in functions related to traits shared by Aphidomorpha, such as association with endosymbionts, and adaptation to plant defenses and phloem-sap-based diet. The ancestral nature of this duplication wave (106-227 Ma) and the lack of sufficiently conserved synteny make it difficult to conclude whether it originated from a whole-genome duplication event or, alternatively, from a burst of large-scale segmental duplications. Genome sequencing of other aphid species belonging to different Aphidomorpha and related lineages may clarify these findings.

Molecular evolution of gland cell types and chemical interactions in animals

Across the Metazoa, the emergence of new ecological interactions has been enabled by the repeated evolution of exocrine glands. Specialized glands have arisen recurrently and with great frequency, even in single genera or species, transforming how animals interact with their environment through trophic resource exploitation, pheromonal communication, chemical defense and parental care. The widespread convergent evolution of animal glands implies that exocrine secretory cells are a hotspot of metazoan cell type innovation. Each evolutionary origin of a novel gland involves a process of 'gland cell type assembly': the stitching together of unique biosynthesis pathways; coordinated changes in secretory systems to enable efficient chemical release; and transcriptional deployment of these machineries into cells constituting the gland. This molecular evolutionary process influences what types of compound a given species is capable of secreting, and, consequently, the kinds of ecological interactions that species can display. Here, we discuss what is known about the evolutionary assembly of gland cell types and propose a framework for how it may happen. We posit the existence of 'terminal selector' transcription factors that program gland function via regulatory recruitment of biosynthetic enzymes and secretory proteins. We suggest ancestral enzymes are initially co-opted into the novel gland, fostering pleiotropic conflict that drives enzyme duplication. This process has yielded the observed pattern of modular, gland-specific biosynthesis pathways optimized for manufacturing specific secretions. We anticipate that single-cell technologies and gene editing methods applicable in diverse species will transform the study of animal chemical interactions, revealing how gland cell types are assembled and functionally configured at a molecular level.