Birth-and-Death Evolution of the Fatty Acyl-CoA Reductase (FAR) Gene Family and Diversification of Cuticular Hydrocarbon Synthesis in Drosophila

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.

Fatty Acid Origin of Insect Pheromones

Pheromones are utilized to a great extent in insects. Many of these pheromones are biosynthesized through a pathway involving fatty acids. This chapter will provide examples where the biosynthetic pathways of fatty acid-derived pheromones have been studied in detail. These include pheromones from Lepidoptera, Coleoptera, and Hymenoptera. Many species of Lepidoptera utilize fatty acids as precursors to pheromones with a functional group that include aldehydes, alcohols, and acetate esters. In addition, the biosynthesis of hydrocarbons will be briefly examined because many insects utilize hydrocarbons or modified hydrocarbons as pheromones.

Strong resistance to β-cyfluthrin in a strain of the beetle Alphitobius diaperinus: a de novo transcriptome analysis

The lesser mealworm, Alphitobius diaperinus, is an invasive tenebrionid beetle and a vector of pathogens. Due to the emergence of insecticide resistance and consequent outbreaks that generate significant phytosanitary and energy costs for poultry farmers, it has become a major insect pest worldwide. To better understand the molecular mechanisms behind this resistance, we studied a strain of A. diaperinus from a poultry house in Brittany that was found to be highly resistant to the β-cyfluthrin. The strain survived β-cyfluthrin exposures corresponding to more than 100 times the recommended dose. We used a comparative de novo RNA-Seq approach to explore genes expression in resistant versus sensitive strains. Our de novo transcriptomic analyses showed that responses to β-cyfluthrin likely involved a whole set of resistance mechanisms. Genes related to detoxification, metabolic resistance, cuticular hydrocarbon biosynthesis and proteolysis were found to be constitutively overexpressed in the resistant compared to the sensitive strain. Follow-up enzymatic assays confirmed that the resistant strain exhibited high basal activities for detoxification enzymes such as cytochrome P450 monooxygenase and glutathione-S-transferase. The in-depth analysis of differentially expressed genes suggests the involvement of complex regulation of signaling pathways. Detailed knowledge of these resistance mechanisms is essential for the establishment of effective pest control.

Adaptation of an arthropod predator to a challenging environment is associated with a loss of a genome-wide plastic transcriptional response

BACKGROUND

Structural and chemical plant defence traits may reduce the efficacy of biological control agents in integrated pest management. Breeding programmes have shown arthropod predators' potential to acclimate to challenging host plants. However, whether and how these predators adapt to novel plant environments remain unclear. Using the predatory mite Phytoseiulus persimilis - herbivorous mite Tetranychus urticae system in an experimental evolution setup, we studied the adaptation mechanisms to tomato and cucumber, plants that possess a distinct repertoire of defensive traits.

RESULTS

Experimental evolution experiments on whole plants revealed that allowing P. persimilis to adapt to tomatoes led to an ~100% larger population size. Independent feeding assays showed that tomato- and cucumber-adapted prey reduced predator fecundity. The deleterious effect of ingesting low-quality prey persisted after adaptation of the predator to both cucumber and tomato. We demonstrated that jasmonic acid (JA)-dependent defences reduce prey quality by evaluating predator performance on prey fed on JA defence-deficient tomato plants. Transcriptomic profiling of the replicated P. persimilis lines showed that long-term propagation on tomato and cucumber plants produces distinctive gene-expression levels. Predator adaptation to tomatoes results in the loss of a large transcriptional response, in which predicted cuticle-building rather than detoxification pathways are affected.

CONCLUSION

We showed that the adaptation of predatory arthropods to a novel, challenging plant does not necessarily occur via the prey, but rather through the physical environment of the plant. We provided first insights into the underlying molecular mechanisms. © 2023 Society of Chemical Industry.

Phylogenomics of the Ecdysteroid Kinase-like (EcKL) Gene Family in Insects Highlights Roles in Both Steroid Hormone Metabolism and Detoxification

The evolutionary dynamics of large gene families can offer important insights into the functions of their individual members. While the ecdysteroid kinase-like (EcKL) gene family has previously been linked to the metabolism of both steroid molting hormones and xenobiotic toxins, the functions of nearly all EcKL genes are unknown, and there is little information on their evolution across all insects. Here, we perform comprehensive phylogenetic analyses on a manually annotated set of EcKL genes from 140 insect genomes, revealing the gene family is comprised of at least 13 subfamilies that differ in retention and stability. Our results show the only two genes known to encode ecdysteroid kinases belong to different subfamilies and therefore ecdysteroid metabolism functions must be spread throughout the EcKL family. We provide comparative phylogenomic evidence that EcKLs are involved in detoxification across insects, with positive associations between family size and dietary chemical complexity, and we also find similar evidence for the cytochrome P450 and glutathione S-transferase gene families. Unexpectedly, we find that the size of the clade containing a known ecdysteroid kinase is positively associated with host plant taxonomic diversity in Lepidoptera, possibly suggesting multiple functional shifts between hormone and xenobiotic metabolism. Our evolutionary analyses provide hypotheses of function and a robust framework for future experimental studies of the EcKL gene family. They also open promising new avenues for exploring the genomic basis of dietary adaptation in insects, including the classically studied coevolution of butterflies with their host plants.

Fatty Acid Elongase Gene PcELO7 is Essential for Lipid Accumulation and Fecundity of Panonychus citri (Acari: Tetranychidae)

RNA interference (RNAi) has been proposed as a promising strategy for sustainable and ecofriendly pest control. The insect cuticle lipids were deposited on the body surface and functioned as a defense against chemical xenobiotics. They consisted of aliphatic compounds, including free fatty acids (FFAs). However, elongase of very long chain fatty acids (ELOs) is essential for FFA biosynthesis; the function of ELO is still unknown in many arthropods, including Panonychus citri (P. citri). In this study, three ELOs were cloned. Developmental-specific mRNA expression results revealed that three PcELOs were highly expressed in egg and adult females. Whereas PcELO7 was dominantly expressed in adult females. Under spirobudiclofen stress, ELOs mRNA expression had different changes, and PcELO7 was down-regulated. The silencing of PcELO7 resulted in a dramatic reduction of oviposition and hatchability. Significant reduction of FFA contents was also examined within PcELO7-repressed P. citri. In addition, we found that PcELO7 mRNA levels were related to fecundity and could affect triacylglycerol (TG) contents. The findings demonstrated that the introduction of dsPcELO7 via oral feeding induced the RNA interference-mediated silencing of a special target gene and could result in mortality and reproduction. In conclusion, PcELO7 is a special RNAi target for P. citri control, and its lethal mechanism might be disturbing lipids biosynthesis.

Knockdown of the ABCG23 Gene Disrupts the Development and Lipid Accumulation of Panonychus citri (Acari/Tetranychidae)

Panonychus citri is a worldwide citrus pest that is currently controlled through the use of insecticides. However, alternative strategies are required to manage P. citri. Recent studies suggest that the ATP-binding cassette (ABC) transporter G subfamily plays a crucial role in transporting cuticular lipids, which are essential for the insect's barrier function against microbial penetration. Therefore, investigating the potential of the ABC transporter G subfamily as a control measure for P. citri could be a promising approach. Based on the genome database, the gene was cloned, and the transcriptional response of ABCG23 for the different developmental stages of P. citri and under spirobudiclofen stress was investigated. Our results showed that the expression level of ABCG23 was significantly lower in adult females exposed to treatment compared to the control and was higher in females than males. The knockdown of ABCG23 using RNAi led to a decrease in the survival rate, fecundity, and TG contents of P. citri. Additionally, a lethal phenotype was characterized by body wrinkling and darkening. These results indicate that ABCG23 may be involved in cuticular lipid transportation and have adverse effects on the development and reproduction of P. citri, providing insight into the discovery of new targets for pest management based on the insect cuticle's penetration barrier function.

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.

Symbiosis preservation: Putative regulation of fatty acyl-CoA reductase by miR-31a within the symbiont harboring bacteriome through tsetse evolution

Tsetse flies are the sole vectors of African trypanosomes. In addition to trypanosomes, tsetse harbor obligate Wigglesworthia glossinidia bacteria that are essential to tsetse biology. The absence of Wigglesworthia results in fly sterility, thus offering promise for population control strategies. Here, microRNA (miRNAs) and mRNA expression are characterized and compared between the exclusive Wigglesworthia-containing bacteriome and adjacent aposymbiotic tissue in females of two evolutionarily distant tsetse species (Glossina brevipalpis and G. morsitans). A total of 193 miRNAs were expressed in either species, with 188 of these expressed in both species, 166 of these were novel to Glossinidae, and 41 miRNAs exhibited comparable expression levels between species. Within bacteriomes, 83 homologous mRNAs demonstrated differential expression between G. morsitans aposymbiotic and bacteriome tissues, with 21 of these having conserved interspecific expression. A large proportion of these differentially expressed genes are involved in amino acid metabolism and transport, symbolizing the essential nutritional role of the symbiosis. Further bioinformatic analyses identified a sole conserved miRNA::mRNA interaction (miR-31a::fatty acyl-CoA reductase) within bacteriomes likely catalyzing the reduction of fatty acids to alcohols which comprise components of esters and lipids involved in structural maintenance. The Glossina fatty acyl-CoA reductase gene family is characterized here through phylogenetic analyses to further understand its evolutionary diversification and the functional roles of members. Further research to characterize the nature of the miR-31a::fatty acyl-CoA reductase interaction may find novel contributions to the symbiosis to be exploited for vector control.

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.

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.

Natural variation at a single gene generates sexual antagonism across fitness components in Drosophila

Mutations with conflicting fitness effects in males and females accumulate in sexual populations, reducing their adaptive capacity.^1,2 Although quantitative genetic studies indicate that sexually antagonistic polymorphisms are common,^3-5 their molecular basis and population genetic properties remain poorly understood.^6,7 Here, we show in fruit flies how natural variation at a single gene generates sexual antagonism through phenotypic effects on cuticular hydrocarbon (CHC) traits that function as both mate signals and protectors against abiotic stress⁸ across a latitudinal gradient. Tropical populations of Drosophila serrata have polymorphic CHCs producing sexual antagonism through opposing but sex-limited effects on these two fitness-related functions. We dissected this polymorphism to a single fatty-acyl CoA reductase gene, DsFAR2-B, that is expressed in oenocyte cells where CHCs are synthesized. RNAi-mediated disruption of the DsFAR2-B ortholog in D. melanogaster oenocytes affected CHCs in a similar way to that seen in D. serrata. Population genomic analysis revealed that balancing selection likely operates at the DsFAR2-B locus in the wild. Our study provides insights into the genetic basis of sexual antagonism in nature and connects sexually varying antagonistic selection on phenotypes with balancing selection on genotypes that maintains molecular variation.

Genetic and genomic architecture of species-specific cuticular hydrocarbon variation in parasitoid wasps

Cuticular hydrocarbons (CHCs) serve two fundamental functions in insects: protection against desiccation and chemical signalling. How the interaction of genes shapes CHC profiles, which are essential for insect survival, adaptation and reproductive success, is still poorly understood. Here we investigate the genetic and genomic basis of CHC biosynthesis and variation in parasitoid wasps of the genus Nasonia. We mapped 91 quantitative trait loci (QTL) explaining the variation of a total of 43 CHCs in F₂ hybrid males from interspecific crosses between three Nasonia species. To identify candidate genes, we localized orthologues of CHC biosynthesis-related genes in the Nasonia genomes. We discovered multiple genomic regions where the location of QTL coincides with the location of CHC biosynthesis-related candidate genes. Most conspicuously, on a region close to the centromere of chromosome 1, multiple CHC biosynthesis-related candidate genes co-localize with several QTL explaining variation in methyl-branched alkanes. The genetic underpinnings behind this compound class are not well understood so far, despite their high potential for encoding chemical information as well as their prevalence in hymenopteran CHC profiles. Our study considerably extends our knowledge on the genetic architecture governing this important compound class, establishing a model for methyl-branched alkane genetics in the Hymenoptera in general.

Differential gene expression in aphids following virus acquisition from plants or from an artificial medium

BACKGROUND

Poleroviruses, such as turnip yellows virus (TuYV), are plant viruses strictly transmitted by aphids in a persistent and circulative manner. Acquisition of either virus particles or plant material altered by virus infection is expected to induce gene expression deregulation in aphids which may ultimately alter their behavior.

RESULTS

By conducting an RNA-Seq analysis on viruliferous aphids fed either on TuYV-infected plants or on an artificial medium containing purified virus particles, we identified several hundreds of genes deregulated in Myzus persicae, despite non-replication of the virus in the vector. Only a few genes linked to receptor activities and/or vesicular transport were common between the two modes of acquisition with, however, a low level of deregulation. Behavioral studies on aphids after virus acquisition showed that M. persicae locomotion behavior was affected by feeding on TuYV-infected plants, but not by feeding on the artificial medium containing the purified virus particles. Consistent with this, genes potentially involved in aphid behavior were deregulated in aphids fed on infected plants, but not on the artificial medium.

CONCLUSIONS

These data show that TuYV particles acquisition alone is associated with a moderate deregulation of a few genes, while higher gene deregulation is associated with aphid ingestion of phloem from TuYV-infected plants. Our data are also in favor of a major role of infected plant components on aphid behavior.

Phenotypic disruption of cuticular hydrocarbon production in hybrids between sympatric species of Hawaiian picture-wing Drosophila

Interspecies hybrids can express phenotypic traits far outside the range of parental species. The atypical traits of hybrids provide insight into differences in the factors that regulate the expression of these traits in the parental species. In some cases, the unusual phenotypic traits of hybrids can lead to phenotypic dysfunction with hybrids experiencing reduced survival or reproduction. Cuticular hydrocarbons (CHCs) in insects are important phenotypic traits that serve several functions, including desiccation resistance and pheromones for mating. We used gas chromatography mass spectrometry to investigate the differences in CHC production between two closely related sympatric Hawaiian picture-wing Drosophila species, Drosophila heteroneura and D. silvestris, and their F1 and backcross hybrid offspring. CHC profiles differed between males of the two species, with substantial sexual dimorphism in D. silvestris but limited sexual dimorphism in D. heteroneura. Surprisingly, F1 hybrids did not produce three CHCs, and the abundances of several other CHCs occurred outside the ranges present in the two parental species. Backcross hybrids produced all CHCs with greater variation than observed in F1 or parental species. Overall, these results suggest that the production of CHCs was disrupted in F1 and backcross hybrids, which may have important consequences for their survival or reproduction.

CYP4G8 is responsible for the synthesis of methyl-branched hydrocarbons in the polyphagous caterpillar of Helicoverpa armigera

Insect cuticular hydrocarbons (CHCs) have dual functions as physical barrier and chemical signals. The last step of CHC biosynthesis is known to be catalyzed by cytochrome P450 CYP4G in a number of insects. Until recently, studies on CYP4Gs in the context of functional evolution are rare. In this study, we analyzed sequence similarity and temporal-spatial expression patterns of the five CYP4G genes in the cotton bollworm Helicoverpa armigera, an important agricultural pest and also typical representative of lepidopteran insects. Moreover, the CRISPR/Cas9-induced knockout was used to clarify the roles of the five CYP4Gs in CHC biosynthesis. Temporal-spatial expression patterns revealed that CYP4G8 was highly expressed at all developmental stages and in most tissues examined. Larvae with CYP4G8 knocked out could not produce methyl-branched CHCs and failed to pupate, while larvae with the other four CYP4G genes knocked out (4G1-type-KO) showed no significant changes in their CHC profiles, weight gain and survival. Comparative transcriptomics revealed that knocking out CYP4G8 affected the global gene expression in larvae, especially down-regulated the expression of genes in the fatty acid biosynthetic pathway, while no significant change in 4G1-type-KO transcriptome was observed. These findings indicate that the five members of the CYP4G subfamily have undergone functional divergence: CYP4G8 maintains the essential function in CHC biosynthesis, while the function of the other four CYP4G genes remains unclear. Intriguingly, CYP4G8 has evolved to be a P450 enzyme responsible for the synthesis of larval methyl-branched hydrocarbons. The observation that CYP4G8 knockout is lethal strongly suggest that CYP4G8 may serve as a candidate target for the development of insecticidal agents for the control of cotton bollworms.

Evolutionary assembly of cooperating cell types in an animal chemical defense system

How the functions of multicellular organs emerge from the underlying evolution of cell types is poorly understood. We deconstructed evolution of an organ novelty: a rove beetle gland that secretes a defensive cocktail. We show how gland function arose via assembly of two cell types that manufacture distinct compounds. One cell type, comprising a chemical reservoir within the abdomen, produces alkane and ester compounds. We demonstrate that this cell type is a hybrid of cuticle cells and ancient pheromone and adipocyte-like cells, executing its function via a mosaic of enzymes from each parental cell type. The second cell type synthesizes benzoquinones using a chimera of conserved cellular energy and cuticle formation pathways. We show that evolution of each cell type was shaped by coevolution between the two cell types, yielding a potent secretion that confers adaptive value. Our findings illustrate how cooperation between cell types arises, generating new, organ-level behaviors.

The Birth-and-Death Evolution of Cytochrome P450 Genes in Bees

The birth-and-death model of multigene family evolution describes how gene families evolve and diversify through duplication and deletion. The cytochrome P450s are one of the most diverse and well-studied multigene families, involved in both physiological and xenobiotic functions. Extensive studies of insect P450 genes have demonstrated their role in insecticide resistance. Bees are thought to experience toxin exposure through their diet of nectar and pollen, as well as the resin-collecting behavior exhibited by some species. Here, we describe the repertoire of P450 genes in the orchid bee Euglossa dilemma. Male orchid bees form perfume bouquets used in courtship displays by collecting volatile compounds, resulting in exposure to compounds known to be toxic. In addition, we conducted phylogenetic and selection analyses across ten bee species encompassing three bee families. We find that social behavior and resin collection are not correlated with the repertoire of P450 present in a bee species. However, our analyses revealed that P450 clades can be classified as stable and unstable, and that genes involved in xenobiotic metabolism are more likely to belong to unstable clades. Furthermore, we find that unstable clades are under more dynamic evolutionary pressures and exhibit signals of adaptive evolution. This work highlights the complexity of multigene family evolution, revealing that multiple factors contribute to the diversification, stability, and dynamics of this gene family. Furthermore, we provide a resource for future detailed studies investigating the function of different P450s in economically important bee species.

DrosoPhyla: Resources for Drosophilid Phylogeny and Systematics

The vinegar fly Drosophila melanogaster is a pivotal model for invertebrate development, genetics, physiology, neuroscience, and disease. The whole family Drosophilidae, which contains over 4,400 species, offers a plethora of cases for comparative and evolutionary studies. Despite a long history of phylogenetic inference, many relationships remain unresolved among the genera, subgenera, and species groups in the Drosophilidae. To clarify these relationships, we first developed a set of new genomic markers and assembled a multilocus data set of 17 genes from 704 species of Drosophilidae. We then inferred a species tree with highly supported groups for this family. Additionally, we were able to determine the phylogenetic position of some previously unplaced species. These results establish a new framework for investigating the evolution of traits in fruit flies, as well as valuable resources for systematics.

Large-scale characterization of sex pheromone communication systems in Drosophila

Insects use sex pheromones as a reproductive isolating mechanism to attract conspecifics and repel heterospecifics. Despite the profound knowledge of sex pheromones, little is known about the coevolutionary mechanisms and constraints on their production and detection. Using whole-genome sequences to infer the kinship among 99 drosophilids, we investigate how phylogenetic and chemical traits have interacted at a wide evolutionary timescale. Through a series of chemical syntheses and electrophysiological recordings, we identify 52 sex-specific compounds, many of which are detected via olfaction. Behavioral analyses reveal that many of the 43 male-specific compounds are transferred to the female during copulation and mediate female receptivity and/or male courtship inhibition. Measurement of phylogenetic signals demonstrates that sex pheromones and their cognate olfactory channels evolve rapidly and independently over evolutionary time to guarantee efficient intra- and inter-specific communication systems. Our results show how sexual isolation barriers between species can be reinforced by species-specific olfactory signals.

Despite the profound knowledge of sex pheromones, little is known about the coevolutionary mechanisms and constraints on their production and detection. Whole-genome sequences from 99 drosophilids, with chemical and behavioural data, show that sex pheromones and their cognate olfactory channels evolve rapidly and independently.

Conversion between 100-million-year-old duplicated genes contributes to rice subspecies divergence

BACKGROUND

Duplicated gene pairs produced by ancient polyploidy maintain high sequence similarity over a long period of time and may result from illegitimate recombination between homeologous chromosomes. The genomes of Asian cultivated rice Oryza sativa ssp. indica (XI) and Oryza sativa ssp. japonica (GJ) have recently been updated, providing new opportunities for investigating ongoing gene conversion events and their impact on genome evolution.

RESULTS

Using comparative genomics and phylogenetic analyses, we evaluated gene conversion rates between duplicated genes produced by polyploidization 100 million years ago (mya) in GJ and XI. At least 5.19-5.77% of genes duplicated across the three rice genomes were affected by whole-gene conversion after the divergence of GJ and XI at ~ 0.4 mya, with more (7.77-9.53%) showing conversion of only portions of genes. Independently converted duplicates surviving in the genomes of different subspecies often use the same donor genes. The ongoing gene conversion frequency was higher near chromosome termini, with a single pair of homoeologous chromosomes, 11 and 12, in each rice genome being most affected. Notably, ongoing gene conversion has maintained similarity between very ancient duplicates, provided opportunities for further gene conversion, and accelerated rice divergence. Chromosome rearrangements after polyploidization are associated with ongoing gene conversion events, and they directly restrict recombination and inhibit duplicated gene conversion between homeologous regions. Furthermore, we found that the converted genes tended to have more similar expression patterns than nonconverted duplicates. Gene conversion affects biological functions associated with multiple genes, such as catalytic activity, implying opportunities for interaction among members of large gene families, such as NBS-LRR disease-resistance genes, contributing to the occurrence of the gene conversion.

CONCLUSION

Duplicated genes in rice subspecies generated by grass polyploidization ~ 100 mya remain affected by gene conversion at high frequency, with important implications for the divergence of rice subspecies.

Eco-genetics of desiccation resistance in Drosophila

Climate change globally perturbs water circulation thereby influencing ecosystems including cultivated land. Both harmful and beneficial species of insects are likely to be vulnerable to such changes in climate. As small animals with a disadvantageous surface area to body mass ratio, they face a risk of desiccation. A number of behavioural, physiological and genetic strategies are deployed to solve these problems during adaptation in various Drosophila species. Over 100 desiccation-related genes have been identified in laboratory and wild populations of the cosmopolitan fruit fly Drosophila melanogaster and its sister species in large-scale and single-gene approaches. These genes are involved in water sensing and homeostasis, and barrier formation and function via the production and composition of surface lipids and via pigmentation. Interestingly, the genetic strategy implemented in a given population appears to be unpredictable. In part, this may be due to different experimental approaches in different studies. The observed variability may also reflect a rich standing genetic variation in Drosophila allowing a quasi-random choice of response strategies through soft-sweep events, although further studies are needed to unravel any underlying principles. These findings underline that D. melanogaster is a robust species well adapted to resist climate change-related desiccation. The rich data obtained in Drosophila research provide a framework to address and understand desiccation resistance in other insects. Through the application of powerful genetic tools in the model organism D. melanogaster, the functions of desiccation-related genes revealed by correlative studies can be tested and the underlying molecular mechanisms of desiccation tolerance understood. The combination of the wealth of available data and its genetic accessibility makes Drosophila an ideal bioindicator. Accumulation of data on desiccation resistance in Drosophila may allow us to create a world map of genetic evolution in response to climate change in an insect genome. Ultimately these efforts may provide guidelines for dealing with the effects of climate-related perturbations on insect population dynamics in the future.

Transcriptomic analyses of the termite, Cryptotermes secundus, reveal a gene network underlying a long lifespan and high fecundity

Organisms are typically characterized by a trade-off between fecundity and longevity. Notable exceptions are social insects. In insect colonies, the reproducing caste (queens) outlive their non-reproducing nestmate workers by orders of magnitude and realize fecundities and lifespans unparalleled among insects. How this is achieved is not understood. Here, we identified a single module of co-expressed genes that characterized queens in the termite species Cryptotermes secundus. It encompassed genes from all essential pathways known to be involved in life-history regulation in solitary model organisms. By manipulating its endocrine component, we tested the recent hypothesis that re-wiring along the nutrient-sensing/endocrine/fecundity axis can account for the reversal of the fecundity/longevity trade-off in social insect queens. Our data from termites do not support this hypothesis. However, they revealed striking links to social communication that offer new avenues to understand the re-modelling of the fecundity/longevity trade-off in social insects.

Advances in deciphering the genetic basis of insect cuticular hydrocarbon biosynthesis and variation

Cuticular hydrocarbons (CHCs) have two fundamental functions in insects. They protect terrestrial insects against desiccation and serve as signaling molecules in a wide variety of chemical communication systems. It has been hypothesized that these pivotal dual traits for adaptation to both desiccation and signaling have contributed to the considerable evolutionary success of insects. CHCs have been extensively studied concerning their variation, behavioral impact, physiological properties, and chemical compositions. However, our understanding of the genetic underpinnings of CHC biosynthesis has remained limited and mostly biased towards one particular model organism (Drosophila). This rather narrow focus has hampered the establishment of a comprehensive view of CHC genetics across wider phylogenetic boundaries. This review attempts to integrate new insights and recent knowledge gained in the genetics of CHC biosynthesis, which is just beginning to incorporate work on more insect taxa beyond Drosophila. It is intended to provide a stepping stone towards a wider and more general understanding of the genetic mechanisms that gave rise to the astonishing diversity of CHC compounds across different insect taxa. Further research in this field is encouraged to aim at better discriminating conserved versus taxon-specific genetic elements underlying CHC variation. This will be instrumental in greatly expanding our knowledge of the origins and variation of genes governing the biosynthesis of these crucial phenotypic traits that have greatly impacted insect behavior, physiology, and evolution.

The evolution of insect metallothioneins

Metallothioneins (MTs) are a family of cysteine-rich metal-binding proteins that are important in the chelating and detoxification of toxic heavy metals. Until now, the short length and the low sequence complexity of MTs have hindered the inference of robust phylogenies, hampering the study of their evolution. To address this longstanding question, we applied an iterative BLAST search pipeline that allowed us to build a unique dataset of more than 300 MT sequences in insects. By combining phylogenetics and synteny analysis, we reconstructed the evolutionary history of MTs in insects. We show that the MT content in insects has been shaped by lineage-specific tandem duplications from a single ancestral MT. Strikingly, we also uncovered a sixth MT, MtnF, in the model organism Drosophila melanogaster. MtnF evolves faster than other MTs and is characterized by a non-canonical length and higher cysteine content. Our methodological framework not only paves the way for future studies on heavy metal detoxification but can also allow us to identify other previously unidentified genes and other low complexity genomic features.

Knockdown of LmCYP303A1 alters cuticular hydrocarbon profiles and increases the susceptibility to desiccation and insecticides in Locusta migratoria

Cytochrome P450 monooxygenases (CYPs) serve many functions in insects, from the regulation of development to xenobiotic detoxification. Several conserved CYPs have been shown to play a role in insect growth and development. CYP303A1 is a highly conserved CYP with a single ortholog in most insects, but its underlying molecular characteristics and specific physiological functions remain poorly understood. In Drosophila melanogaster and Locusta migratoria, CYP303A1 is indispensable for eclosion to adult. Here, we report additional functions of the locust gene LmCYP303A1 in nymphal molts, cuticular lipid deposition and insecticide penetration. RT-qPCR revealed that LmCYP303A1 had a high expression level before ecdysis and was highly expressed in integument, wing pads, foregut and hindgut. Suppression of LmCYP303A1 expression by RNA interference (RNAi) caused a lethal phenotype with molting defect from nymph to nymph. In addition, LmCYP303A1 RNAi resulted in locusts being more susceptible to desiccation and to insecticide toxicity. Furthermore, knockdown of LmCYP303A1 efficiently suppressed the transcript level of key genes (ELO7, FAR15 and CYP4G102) responsible for cuticular hydrocarbon (CHC) synthesis, which led to a decrease in some CHC levels. Taken together, our results suggest that one of the functions of LmCYP303A1 is to regulate the biosynthesis of CHC, which plays critical roles in protecting locusts from water loss and insecticide penetration.

The fatty acid elongase gene LmELO7 is required for hydrocarbon biosynthesis and cuticle permeability in the migratory locust, Locusta migratoria

Insect cuticular lipids are a complex cocktail of highly diverse cuticular hydrocarbons (CHCs), which form a hydrophobic surface coat to maintain water balance and to prevent desiccation and penetration of exogenous substances. Fatty acid elongases (ELOs) are key enzymes that participate in a common CHC synthesis pathway in insects. However, the importance of ELOs for CHC synthesis and function remains understudied. Using transcriptomic data, we have identified seven ELO genes (LmELO1-7) in the migratory locust Locusta migratoria. We determined their tissue-specific and temporal expression profiles in fifth instar nymphs. As we are interested in cuticle barrier formation, we performed RNA interference against LmELO7, which is mainly expressed in the integument. Suppression of LmELO7 significantly decreased its expression and caused lethality during or shortly after molting. CHC quantification by GC-MS analysis indicated that suppression of LmELO7 resulted in a decrease in total CHC amounts. By consequence, CHC deficiency reduced desiccation resistance and enhanced cuticle permeability in LmELO7-suppressed L. migratoria. Interestingly, LmELO7 expression is induced at low air humidity. Our results indicate that LmELO7 plays a vital role in the production of CHCs and, hence, cuticle permeability. Induction of LmELO7 expression in drought conditions suggests a key role of this gene in regulating desiccation resistance. This work is expected to help developing new strategies for insect pest management based on CHC function.

Climate change and the genetics of insecticide resistance

Changes in global temperature and humidity as a result of climate change are producing rapid evolutionary changes in many animal species, including agricultural pests and disease vectors, leading to changes in allele frequencies of genes involved in thermotolerance and desiccation resistance. As some of these genes have pleiotropic effects on insecticide resistance, climate change is likely to affect insecticide resistance in the field. In this review, we discuss how the interactions between adaptation to climate change and resistance to insecticides can affect insecticide resistance in the field using examples in phytophagous and hematophagous pest insects, focusing on the effects of increased temperature and increased aridity. We then use detailed genetic and mechanistic studies in the model insect, Drosophila melanogaster, to explain the mechanisms underlying this phenomenon. We suggest that tradeoffs or facilitation between adaptation to climate change and resistance to insecticides can alter insecticide resistance allele frequencies in the field. The dynamics of these interactions will need to be considered when managing agricultural pests and disease vectors in a changing climate. © 2019 Society of Chemical Industry.

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.

Mating-Induced Differential Expression in Genes Related to Reproduction and Immunity in Spodoptera litura (Lepidoptera: Noctuidae) Female Moths

Mating promotes reproductive activity, which may impact immune performance. Paradoxically, mating frequently challenges females' immunity (e.g., infections). Therefore, studies of postmating resource allocation between reproduction and survival are likely to shed new light on life-history trade-off and sexual selection. Here, we used RNAseq to test whether and how mating affected mRNA expression in genes related to reproduction and immunity in Spodoptera litura female moths. Results show a divergent change in the differentially expressed genes (DEGs) between reproduction and immunity: the immune response was largely downregulated shortly after mating (~6 h postmating), which has some recovery at 24 h postmating; reproductive response is trivial shortly after mating (~6 h postmating), but it largely upregulated at 24 h postmating (e.g., egg maturation related genes were highly upregulated). Considering the fact that most of the total DEGs downregulated from 0 to 6 h postmating (from 51/68 to 214/260) but most of the total DEGs upregulated at 24 h postmating (816/928), it is possible that trade-offs between reproduction and immunity occurred in mated females. For example, they may shut down immunity to favor sperm storage and save limited resources to support the increased energy required in reproduction (e.g., egg maturation and oviposition). Mating-induced infections should be trivial due to low polyandry in S. litura. A reduced immune defense may have no threat to S. litura survival but may benefit reproduction significantly. Furthermore, obvious expression changes were detected in genes related to hormone production, suggesting that endocrine changes could play important roles in postmating responses.