Genetic basis of octanoic acid resistance in Drosophila sechellia: functional analysis of a fine-mapped region

Genetics of flight in spongy moths (Lymantria dispar ssp.): functionally integrated profiling of a complex invasive trait

BACKGROUND

Flight can drastically enhance dispersal capacity and is a key trait defining the potential of exotic insect species to spread and invade new habitats. The phytophagous European spongy moths (ESM, Lymantria dispar dispar) and Asian spongy moths (ASM; a multi-species group represented here by L. d. asiatica and L. d. japonica), are globally invasive species that vary in adult female flight capability-female ASM are typically flight capable, whereas female ESM are typically flightless. Genetic markers of flight capability would supply a powerful tool for flight profiling of these species at any intercepted life stage. To assess the functional complexity of spongy moth flight and to identify potential markers of flight capability, we used multiple genetic approaches aimed at capturing complementary signals of putative flight-relevant genetic divergence between ESM and ASM: reduced representation genome-wide association studies, whole genome sequence comparisons, and developmental transcriptomics. We then judged the candidacy of flight-associated genes through functional analyses aimed at addressing the proximate demands of flight and salient features of the ecological context of spongy moth flight evolution.

RESULTS

Candidate gene sets were typically non-overlapping across different genetic approaches, with only nine gene annotations shared between any pair of approaches. We detected an array of flight-relevant functional themes across gene sets that collectively suggest divergence in flight capability between European and Asian spongy moth lineages has coincided with evolutionary differentiation in multiple aspects of flight development, execution, and surrounding life history. Overall, our results indicate that spongy moth flight evolution has shaped or been influenced by a large and functionally broad network of traits.

CONCLUSIONS

Our study identified a suite of flight-associated genes in spongy moths suited to exploration of the genetic architecture and evolution of flight, or validation for flight profiling purposes. This work illustrates how complementary genetic approaches combined with phenotypically targeted functional analyses can help to characterize genetically complex traits.

Evolution of fatty acid taste in drosophilids

Comparative studies of related but ecologically distinct species can reveal how the nervous system evolves to drive behaviors that are particularly suited to certain environments. Drosophila melanogaster is a generalist that feeds and oviposits on most overripe fruits. A sibling species, D. sechellia, is an obligate specialist of Morinda citrifolia (noni) fruit, which is rich in fatty acids (FAs). To understand evolution of noni taste preference, we characterized behavioral and cellular responses to noni-associated FAs in three related drosophilids. We find that mixtures of sugar and noni FAs evoke strong aversion in the generalist species but not in D. sechellia. Surveys of taste sensory responses reveal noni FA- and species-specific differences in at least two mechanisms-bitter neuron activation and sweet neuron inhibition-that correlate with shifts in noni preference. Chemoreceptor mutant analysis in D. melanogaster predicts that multiple genetic changes account for evolution of gustatory preference in D. sechellia.

Genetic variation in disease resistance in Drosophila spp. is mitigated in Drosophila sechellia by specialization to a toxic host

We found that Drosophila species vary in their susceptibility to the broad-spectrum entomopathogen, Metarhizium anisopliae (strain Ma549). Generalist species were generally more resistant than dietary specialists, with the cactophilic Drosophila buzzatii and Drosophila sechellia, a specialist of the Morinda citrifolia (Morinda) fruit, being most susceptible. Morinda fruit is reported to be toxic to most herbivores because it contains Octanoic Acid (OA). We confirmed that OA is toxic to Drosophila spp., other than D. sechellia, and we also found that OA is highly toxic to entomopathogenic fungi including Ma549 and Beauveria bassiana. Drosophila sechellia fed a diet containing OA, even at levels much less than found in Morinda fruit, had greatly reduced susceptibility to Ma549. This suggests that specializing to Morinda may have provided an enemy-free space, reducing adaptive prioritization on a strong immune response. Our results demonstrate that M. anisopliae and Drosophila species with divergent lifestyles provide a versatile model system for understanding the mechanisms of host-pathogen interactions at different scales and in environmental context.

Pheromone sensing in Drosophila requires support cell-expressed Osiris 8

Background

The nose of most animals comprises multiple sensory subsystems, which are defined by the expression of different olfactory receptor families. Drosophila melanogaster antennae contain two morphologically and functionally distinct subsystems that express odorant receptors (Ors) or ionotropic receptors (Irs). Although these receptors have been thoroughly characterized in this species, the subsystem-specific expression and roles of other genes are much less well-understood.

Results

Here we generate subsystem-specific transcriptomic datasets to identify hundreds of genes, encoding diverse protein classes, that are selectively enriched in either Or or Ir subsystems. Using single-cell antennal transcriptomic data and RNA in situ hybridization, we find that most neuronal genes—other than sensory receptor genes—are broadly expressed within the subsystems. By contrast, we identify many non-neuronal genes that exhibit highly selective expression, revealing substantial molecular heterogeneity in the non-neuronal cellular components of the olfactory subsystems. We characterize one Or subsystem-specific non-neuronal molecule, Osiris 8 (Osi8), a conserved member of a large, insect-specific family of transmembrane proteins. Osi8 is expressed in the membranes of tormogen support cells of pheromone-sensing trichoid sensilla. Loss of Osi8 does not have obvious impact on trichoid sensillar development or basal neuronal activity, but abolishes high sensitivity responses to pheromone ligands.

Conclusions

This work identifies a new protein required for insect pheromone detection, emphasizes the importance of support cells in neuronal sensory functions, and provides a resource for future characterization of other olfactory subsystem-specific genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01425-w.

Potential Action Mechanism and Inhibition Efficacy of Morinda citrifolia Essential Oil and Octanoic Acid against Stagonosporopsis cucurbitacearum Infestations

The use of plant-based products has been shown to efficiently inhibit fungi-mediated diseases in agricultural crops. Here, we extracted and evaluated the composition of noni, Morinda citrifolia L., essential oil and assessed its activities against Stagonosporopsis cucurbitacearum in Cucumis melo L. Using in silico molecular approaches, potential interactions between the essential oil major components and S. cucurbitacearum tyrosine-tRNA ligase were predicted. Finally, we also measured the potential interference of plant physiology (the stomatal conductance and net photosynthesis) mediated by the application of the M. citrifolia essential oil. Chromatographic analysis revealed that octanoic acid (75.8%), hexanoic acid (12.8%), and isobutyl pent-4-enyl carbonate (3.1%) were the major essential oil compounds. Octanoic acid and noni essential oil, when used as preventive measures, reduce fungal mycelial growth at a concentration of 5 mg/mL without causing significant damage to the treated leaves, which reinforces their efficacies as preventive tools against S. cucurbitacearum. Molecular docking analyses predicted very stable interactions between the major essential oil constituents and S. cucurbitacearum tyrosine-tRNA ligase, suggesting the interference of these plant-based molecules upon enzyme activation. Octanoic acid and M. citrifolia essential oil at concentrations of 20 mg/mL decreased the stomatal conductance and net photosynthesis rate of melon plants, resulting in robust phytotoxicity. Collectively, our findings indicated that despite the phytotoxicity risks at higher concentrations, M. citrifolia essential oil and octanoic acid, have potential as alternative tools for the integrative management of S. cucurbitacearum.

Ionotropic receptors in the turnip moth Agrotis segetum respond to repellent medium-chain fatty acids

Background

In insects, airborne chemical signals are mainly detected by two receptor families, odorant receptors (ORs) and ionotropic receptors (IRs). Functions of ORs have been intensively investigated in Diptera and Lepidoptera, while the functions and evolution of the more ancient IR family remain largely unexplored beyond Diptera.

Results

Here, we identified a repertoire of 26 IRs from transcriptomes of female and male antennae, and ovipositors in the moth Agrotis segetum. We observed that a large clade formed by IR75p and IR75q expansions is closely related to the acid-sensing IRs identified in Diptera. We functionally assayed each of the five AsegIRs from this clade using Xenopus oocytes and found that two receptors responded to the tested ligands. AsegIR75p.1 responded to several compounds but hexanoic acid was revealed to be the primary ligand, and AsegIR75q.1 responded primarily to octanoic acid, and less so to nonanoic acid. It has been reported that the C₆-C₁₀ medium-chain fatty acids repel various insects including many drosophilids and mosquitos. We show that the C₆-C₁₀ medium-chain fatty acids elicited antennal responses of both sexes of A. segetum, while only octanoic acid had repellent effect to the moths in a behavioral assay. In addition, using fluorescence in situ hybridization, we demonstrated that the five IRs and their co-receptor AsegIR8a are not located in coeloconic sensilla as found in Drosophila, but in basiconic or trichoid sensilla.

Conclusions

Our results significantly expand the current knowledge of the insect IR family. Based on the functional data in combination with phylogenetic analysis, we propose that subfunctionalization after gene duplication plays an important role in the evolution of ligand specificities of the acid-sensing IRs in Lepidoptera.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01235-0.

Gene expression in Lucilia sericata (Diptera: Calliphoridae) larvae exposed to Pseudomonas aeruginosa and Acinetobacter baumannii identifies shared and microbe-specific induction of immune genes

Antibiotic resistance is a continuing challenge in medicine. There are various strategies for expanding antibiotic therapeutic repertoires, including the use of blow flies. Their larvae exhibit strong antibiotic and antibiofilm properties that alter microbiome communities. One species, Lucilia sericata, is used to treat problematic wounds due to its debridement capabilities and its excretions and secretions that kill some pathogenic bacteria. There is much to be learned about how L. sericata interacts with microbiomes at the molecular level. To address this deficiency, gene expression was assessed after feeding exposure (1 h or 4 h) to two clinically problematic pathogens: Pseudomonas aeruginosa and Acinetobacter baumannii. The results identified immunity-related genes that were differentially expressed when exposed to these pathogens, as well as non-immune genes possibly involved in gut responses to bacterial infection. There was a greater response to P. aeruginosa that increased over time, while few genes responded to A. baumannii exposure, and expression was not time-dependent. The response to feeding on pathogens indicates a few common responses and features distinct to each pathogen, which is useful in improving the wound debridement therapy and helps to develop biomimetic alternatives.

Drosophila sechellia: A Genetic Model for Behavioral Evolution and Neuroecology

Defining the mechanisms by which animals adapt to their ecological niche is an important problem bridging evolution, genetics, and neurobiology. We review the establishment of a powerful genetic model for comparative behavioral analysis and neuroecology, Drosophila sechellia. This island-endemic fly species is closely related to several cosmopolitan generalists, including Drosophila melanogaster, but has evolved extreme specialism, feeding and reproducing exclusively on the noni fruit of the tropical shrub Morinda citrifolia. We first describe the development and use of genetic approaches to facilitate genotype/phenotype associations in these drosophilids. Next, we survey the behavioral, physiological, and morphological adaptations of D. sechellia throughout its life cycle and outline our current understanding of the genetic and cellular basis of these traits. Finally, we discuss the principles this knowledge begins to establish in the context of host specialization, speciation, and the neurobiology of behavioral evolution and consider open questions and challenges in the field.

Genomics analysis of hexanoic acid exposure in Drosophila species

Drosophila sechellia is a dietary specialist endemic to the Seychelles islands that has evolved to consume the fruit of Morinda citrifolia. When ripe, the fruit of M. citrifolia contains octanoic acid and hexanoic acid, two medium-chain fatty acid volatiles that deter and are toxic to generalist insects. Drosophila sechellia has evolved resistance to these volatiles allowing it to feed almost exclusively on this host plant. The genetic basis of octanoic acid resistance has been the focus of multiple recent studies, but the mechanisms that govern hexanoic acid resistance in D. sechellia remain unknown. To understand how D. sechellia has evolved to specialize on M. citrifolia fruit and avoid the toxic effects of hexanoic acid, we exposed adult D. sechellia, D. melanogaster and D. simulans to hexanoic acid and performed RNA sequencing comparing their transcriptional responses to identify D. sechellia specific responses. Our analysis identified many more genes responding transcriptionally to hexanoic acid in the susceptible generalist species than in the specialist D. sechellia. Interrogation of the sets of differentially expressed genes showed that generalists regulated the expression of many genes involved in metabolism and detoxification whereas the specialist primarily downregulated genes involved in the innate immunity. Using these data, we have identified interesting candidate genes that may be critically important in aspects of adaptation to their food source that contains high concentrations of HA. Understanding how gene expression evolves during dietary specialization is crucial for our understanding of how ecological communities are built and how evolution shapes trophic interactions.

A knockout screen of genes expressed specifically in Ae. aegypti pupae reveals a critical role for stretchin in mosquito flight

Aedes aegypti is a critical vector for transmitting Zika, dengue, chikungunya, and yellow fever viruses to humans. Genetic strategies to limit mosquito survival based upon sex distortion or disruption of development may be valuable new tools to control Ae. aegypti populations. We identified six genes with expression limited to pupal development; osi8 and osi11 (Osiris protein family), CPRs and CPF (cuticle protein family), and stretchin (a muscle protein). Heritable CRISPR/Cas9-mediated gene knockout of these genes did not reveal any defects in pupal development. However, stretchin-null mutations (strn^Δ35/Δ41) resulted in flightless mosquitoes with an abnormal open wing posture. The inability of adult strn^Δ35/Δ41 mosquitoes to fly restricted their escape from aquatic rearing media following eclosion, and substantially reduced adult survival rates. Transgenic strains which contain the EGFP marker gene under the control of strn regulatory regions (0.8 kb, 1.4 kb, and 2.2 kb upstream, respectively), revealed the gene expression pattern of strn in muscle-like tissues in the thorax during late morphogenesis from L₄ larvae to young adults. We demonstrated that Ae. aegypti pupae-specific strn is critical for adult mosquito flight capability and a key late-acting lethal target for mosquito-borne disease control.

Seven Questions on the Chemical Ecology and Neurogenetics of Resource-Mediated Speciation

Adaptation to different environments can result in reproductive isolation between populations and the formation of new species. Food resources are among the most important environmental factors shaping local adaptation. The chemosensory system, the most ubiquitous sensory channel in the animal kingdom, not only detects food resources and their chemical composition, but also mediates sexual communication and reproductive isolation in many taxa. Chemosensory divergence may thus play a crucial role in resource-mediated adaptation and speciation. Understanding how the chemosensory system can facilitate resource-mediated ecological speciation requires integrating mechanistic studies of the chemosensory system with ecological studies, to link the genetics and physiology of chemosensory properties to divergent adaptation. In this review, we use examples of insect research to present seven key questions that can be used to understand how the chemosensory system can facilitate resource-mediated ecological speciation in consumer populations.

Does Drosophila sechellia escape parasitoid attack by feeding on a toxic resource?

Host shifts can drastically change the selective pressures that animals experience from their environment. Drosophila sechellia is a species restricted to the Seychelles islands, where it specializes on the fruit Morinda citrifolia (noni). This fruit is known to be toxic to closely related Drosophila species, including D. melanogaster and D. simulans, releasing D. sechellia from interspecific competition when breeding on this substrate. Previously, we showed that larvae of D. sechellia are unable to mount an effective immunological response against wasp attack, while larvae of closely-related species can defend themselves from parasitoid attack by melanotic encapsulation. We hypothesized that this inability constitutes a trait loss due to a reduced risk of parasitoid attack in noni. Here we present a lab experiment and field survey aimed to test the hypothesis that specialization on noni has released D. sechellia from the antagonistic interaction with its larval parasitoids. Our results from the lab experiment suggest that noni may be harmful to parasitoid wasps. Our results from the field survey indicate that D. sechellia was found in ripe noni, whereas another Drosophila species, D. malerkotliana, was present in unripe and overripe stages. Parasitic wasps of the species Leptopilina boulardi emerged from overripe noni, where D. malerkotliana was the most abundant host, but not from ripe noni. These results indicate that the specialization of D. sechellia on noni has indeed drastically altered its ecological interactions, leading to a relaxation in the selection pressure to maintain parasitoid resistance.

Functional insights from the GC-poor genomes of two aphid parasitoids, Aphidius ervi and Lysiphlebus fabarum

BACKGROUND

Parasitoid wasps have fascinating life cycles and play an important role in trophic networks, yet little is known about their genome content and function. Parasitoids that infect aphids are an important group with the potential for biological control. Their success depends on adapting to develop inside aphids and overcoming both host aphid defenses and their protective endosymbionts.

RESULTS

We present the de novo genome assemblies, detailed annotation, and comparative analysis of two closely related parasitoid wasps that target pest aphids: Aphidius ervi and Lysiphlebus fabarum (Hymenoptera: Braconidae: Aphidiinae). The genomes are small (139 and 141 Mbp) and the most AT-rich reported thus far for any arthropod (GC content: 25.8 and 23.8%). This nucleotide bias is accompanied by skewed codon usage and is stronger in genes with adult-biased expression. AT-richness may be the consequence of reduced genome size, a near absence of DNA methylation, and energy efficiency. We identify missing desaturase genes, whose absence may underlie mimicry in the cuticular hydrocarbon profile of L. fabarum. We highlight key gene groups including those underlying venom composition, chemosensory perception, and sex determination, as well as potential losses in immune pathway genes.

CONCLUSIONS

These findings are of fundamental interest for insect evolution and biological control applications. They provide a strong foundation for further functional studies into coevolution between parasitoids and their hosts. Both genomes are available at https://bipaa.genouest.org.

Effects of Osiris9a on Silk Properties in Bombyx mori Determined by Transgenic Overexpression

Osiris is an insect-specific gene family with multiple biological roles in development, phenotypic polymorphism, and protection. In the silkworm, we have previously identified twenty-five Osiris genes with high evolutionary conservation and remarkable synteny among several insects. Bombxy moriOsiris9a (BmOsi9a) is expressed only in the silk gland, particularly in the middle silk gland (MSG). However, the biological function of BmOsi9a is still unknown. In this study, we overexpressed BmOsi9a in the silk gland by germline transgene expression. BmOsi9a was overexpressed not only in the MSG but also in the posterior silk gland (PSG). Interestingly, BmOsi9a could be secreted into the lumen in the MSG but not in the PSG. In the silk fiber, overexpressed BmOsi9a interacted with Sericin1 in the MSG, as confirmed by a co-immunoprecipitation assay. The overexpression of BmOsi9a altered the secondary structure and crystallinity of the silk fiber, thereby changing the mechanical properties. These results provide insight into the mechanisms underlying silk proteins secretion and silk fiber formation.

Evolutionary genomics of host plant adaptation: insights from Drosophila

Variation in gene expression in response to the use of alternate host plants can reveal genetic and physiological mechanisms explaining why insect-host relationships vary from host specialism to generalism. Interpreting transcriptome variation relies on well-annotated genomes, making drosophilids valuable model systems, particularly those species with tractable ecological associations. Patterns of whole genome expression and alternate gene splicing in response to growth on different hosts have revealed expression of gene networks of known detoxification genes as well as novel functionally enriched genes of diverse metabolic and structural functions. Integrating trancriptomic responses with fitness differences and levels of phenotypic plasticity in response to alternate hosts will help to reveal the general nature of genotype-phenotype relationships.

Derived esterase activity in Drosophila sechellia contributes to evolved octanoic acid resistance

The dietary specialist fruit fly Drosophila sechellia has evolved resistance to the secondary defence compounds produced by the fruit of its host plant, Morinda citrifolia. The primary chemicals that contribute to lethality of M. citrifolia are the medium-chain fatty acids octanoic acid (OA) and hexanoic acid. At least five genomic regions contribute to this adaptation in D. sechellia and whereas the fine-mapped major effect locus for OA resistance on chromosome 3R has been thoroughly analysed, the remaining four genomic regions that contribute to toxin resistance remain uncharacterized. To begin to identify the genetic basis of toxin resistance in this species, we removed the function of well-known detoxification gene families to determine whether they contribute to toxin resistance. Previous work found that evolution of cytochrome P450 enzymatic activity or expression is not responsible for the OA resistance in D. sechellia. Here, we tested the role of the two other major detoxification gene families in resistance to Morinda fruit toxins - glutathione-S-transferases and esterases - through the use of the pesticide synergists diethyl maleate and tribufos that inhibit the function of these gene families. This work suggests that one or more esterase(s) contribute to evolved OA resistance in D. sechellia.

Genomics Analysis of L-DOPA Exposure in Drosophila sechellia

Drosophila sechellia is a dietary specialist fruit fly that evolved from a generalist ancestor to specialize on the toxic fruit of Morinda citrifolia This species pair has been the subject of numerous studies where the goal has largely been to determine the genetic basis of adaptations associated with host specialization. Because one of the most striking features of M. citrifolia fruit is the production of toxic volatile compounds that kill insects, most genomic studies in D. sechellia to date have focused on gene expression responses to the toxic compounds in its food. In this study, we aim to identify new genes important for host specialization by profiling gene expression response to 3,4-dihydroxyphenylalanine (L-DOPA). Recent work found it to be highly abundant in M. citrifolia, critical for reproductive success of D. sechellia, and supplementation of diet with the downstream pathway product dopamine can influence toxin resistance phenotypes in related species. Here we used a combination of functional genetics and genomics techniques to identify new genes that are important for D. sechellia ecological adaptation to this new niche. We show that L-DOPA exposure can affect toxin resistance phenotypes, identify genes with plastic responses to L-DOPA exposure, and functionally test an identified candidate gene. We found that knock-down of Esterase 6 (Est6) in a heterologous species alters toxin resistance suggesting Est6 may play an important role in D. sechellia host specialization.

Host use and host shifts in Drosophila

Over a thousand Drosophila species have radiated onto a wide range of feeding and breeding sites. These radiations involve adaptations for locating, accepting, and growing in hosts with highly differing characteristics. In a number of species, owing to the availability of sequenced genomes, particular steps in host specialization and genes that control them, are being identified. Many cases of specialization involve the ability to detoxify some component of the host. Examples include Drosophila sechellia and the octanoic acid in Morinda citrifolia, alpha-amanitin in mycophagous drosophilids, and the alkaloids in cactophilic species. Owing to the known ecologies of many species for which genomes exist, the Drosophila model system provides an unprecedented opportunity to simultaneously examine the genes underlying HOST LOCATION, HOST ACCEPTANCE and HOST USE, the types of selection acting upon them and any coevolutionary interactions among the genes underlying these steps.

Investigating the role of Osiris genes in Drosophila sechellia larval resistance to a host plant toxin

The underlying genetic basis of adaptive phenotypic changes is generally poorly understood, yet a growing number of case studies are beginning to shed light on important questions about the molecular nature and pleiotropy of such changes. We use Drosophila sechellia, a dietary specialist fruit fly that evolved to specialize on a single toxic host plant, Morinda citrifolia, as a model for adaptive phenotypic change and seek to determine the genetic basis of traits associated with host specialization in this species. The fruit of M. citrifolia is toxic to other drosophilids, primarily due to high levels of the defense chemical octanoic acid (OA), yet D. sechellia has evolved resistance to OA. Our prior work identified three Osiris family genes that reside in a fine-mapped QTL for OA resistance: Osiris 6 (Osi6), Osi7, and Osi8, which can alter OA resistance in adult D. melanogaster when knocked down with RNA interference suggesting they may contribute to OA resistance in D. sechellia. Genetic mapping identified overlapping genomic regions involved in larval and adult OA resistance in D. sechellia, yet it remains unknown whether Osiris genes contribute to resistance in both life stages. Furthermore, because multiple genomic regions contribute to OA resistance, we aim to identify other gene(s) involved in this adaptation. Here, we identify candidate larval OA resistance genes using RNA sequencing to measure genome-wide differential gene expression in D. sechellia larvae after exposure to OA and functionally test identified genes for a role in OA resistance. We then test the Osiris genes previously shown to alter adult OA resistance for effects on OA resistance in larvae. We found that Osi8 knockdown decreased OA resistance in D. melanogaster larvae. These data suggest that evolved changes in Osi8 could impact OA resistance in multiple life stages while Osi6 and Osi7 may only impact adult resistance to OA.

Characterization and expression profiling of microRNAs in response to plant feeding in two host-plant strains of the lepidopteran pest Spodoptera frugiperda

Background

A change in the environment may impair development or survival of living organisms leading them to adapt to the change. The resulting adaptation trait may reverse, or become fixed in the population leading to evolution of species. Deciphering the molecular basis of adaptive traits can thus give evolutionary clues. In phytophagous insects, a change in host-plant range can lead to emergence of new species. Among them, Spodoptera frugiperda is a major agricultural lepidopteran pest consisting of two host-plant strains having diverged 3 MA, based on mitochondrial markers. In this paper, we address the role of microRNAs, important gene expression regulators, in response to host-plant change and in adaptive evolution.

Results

Using small RNA sequencing, we characterized miRNA repertoires of the corn (C) and rice (R) strains of S. frugiperda, expressed during larval development on two different host-plants, corn and rice, in the frame of reciprocal transplant experiments. We provide evidence for 76 and 68 known miRNAs in C and R strains and 139 and 171 novel miRNAs. Based on read counts analysis, 34 of the microRNAs were differentially expressed in the C strain larvae fed on rice as compared to the C strain larvae fed on corn. Twenty one were differentially expressed on rice compared to corn in R strain. Nine were differentially expressed in the R strain compared to C strain when reared on corn. A similar ratio of microRNAs was differentially expressed between strains on rice. We could validate experimentally by QPCR, variation in expression of the most differentially expressed candidates. We used bioinformatics methods to determine the target mRNAs of known microRNAs. Comparison with the mRNA expression profile during similar reciprocal transplant experiment revealed potential mRNA targets of these host-plant regulated miRNAs.

Conclusions

In the current study, we performed the first systematic analysis of miRNAs in Lepidopteran pests feeding on host-plants. We identified a set of the differentially expressed miRNAs that respond to the plant diet, or differ constitutively between the two host plant strains. Among the latter, the ones that are also deregulated in response to host-plant are molecular candidates underlying a complex adaptive trait.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5119-6) contains supplementary material, which is available to authorized users.

The molecular genetic basis of herbivory between butterflies and their host plants

Interactions between herbivorous insects and their host-plants are a central component of terrestrial food webs and a critical topic in agriculture, where a substantial fraction of potential crop yield is lost annually to pests. Important insights into plant-insect interactions have come from research on specific plant defenses and insect detoxification mechanisms. Yet, much remains unknown about the molecular mechanisms that mediate plant-insect interactions. Here we use multiple genome-wide approaches to map the molecular basis of herbivory from both plant and insect perspectives, focusing on butterflies and their larval host-plants. Parallel genome-wide association studies in the Cabbage White butterfly, Pieris rapae, and its host-plant, Arabidopsis thaliana, pinpointed a small number of butterfly and plant genes that influenced herbivory. These genes, along with much of the genome, were regulated in a dynamic way over the time course of the feeding interaction. Comparative analyses, including diverse butterfly/plant systems, showed a variety of genome-wide responses to herbivory, yet a core set of highly conserved genes in butterflies as well as their host-plants. These results greatly expand our understanding of the genomic causes and evolutionary consequences of ecological interactions across two of nature’s most diverse taxa, butterflies and flowering plants.

Supervised machine learning reveals introgressed loci in the genomes of Drosophila simulans and D. sechellia

Hybridization and gene flow between species appears to be common. Even though it is clear that hybridization is widespread across all surveyed taxonomic groups, the magnitude and consequences of introgression are still largely unknown. Thus it is crucial to develop the statistical machinery required to uncover which genomic regions have recently acquired haplotypes via introgression from a sister population. We developed a novel machine learning framework, called FILET (Finding Introgressed Loci via Extra-Trees) capable of revealing genomic introgression with far greater power than competing methods. FILET works by combining information from a number of population genetic summary statistics, including several new statistics that we introduce, that capture patterns of variation across two populations. We show that FILET is able to identify loci that have experienced gene flow between related species with high accuracy, and in most situations can correctly infer which population was the donor and which was the recipient. Here we describe a data set of outbred diploid Drosophila sechellia genomes, and combine them with data from D. simulans to examine recent introgression between these species using FILET. Although we find that these populations may have split more recently than previously appreciated, FILET confirms that there has indeed been appreciable recent introgression (some of which might have been adaptive) between these species, and reveals that this gene flow is primarily in the direction of D. simulans to D. sechellia.

Understanding the extent to which species or diverged populations hybridize in nature is crucially important if we are to understand the speciation process. Accordingly numerous research groups have developed methodology for finding the genetic evidence of such introgression. In this report we develop a supervised machine learning approach for uncovering loci which have introgressed across species boundaries. We show that our method, FILET, has greater accuracy and power than competing methods in discovering introgression, and in addition can detect the directionality associated with the gene flow between species. Using whole genome sequences from Drosophila simulans and Drosophila sechellia we show that FILET discovers quite extensive introgression between these species that has occurred mostly from D. simulans to D. sechellia. Our work highlights the complex process of speciation even within a well-studied system and points to the growing importance of supervised machine learning in population genetics.

Tracheal expression of Osiris gene family in Drosophila

The Drosophila trachea is a premier genetic system to investigate the fundamental mechanisms of tubular organ formation. Development of the trachea consists of the well understood early branch specification and migration processes, and the less clear later branch maturation process including the apical membrane expansion, cytoskeleton rearrangement, luminal matrix clearance, and air-filling. We identified seven members of the Osiris (Osi) gene family with obvious tracheal expression in Drosophila. In addition, HA-tagged Osi proteins are highly concentrated in vesicle-like structures at and near the apical membrane. Osi proteins are predicted to contain endocytic signals and transmembrane domains. The localization of Osi proteins is consistent with these predictions. Interestingly, the Drosophila tracheal tube maturation process also occurs at the apical membrane. Taken together, the localization of Osi proteins suggest that these proteins are likely involved in tube maturation through vesicular trafficking or interacting with other apical membrane proteins.

Conserved roles of Osiris genes in insect development, polymorphism and protection

Much of the variation among insects is derived from the different ways that chitin has been moulded to form rigid structures, both internal and external. In this study, we identify a highly conserved expression pattern in an insect-only gene family, the Osiris genes, that is essential for development, but also plays a significant role in phenotypic plasticity and in immunity/toxicity responses. The majority of Osiris genes exist in a highly syntenic cluster, and the cluster itself appears to have arisen very early in the evolution of insects. We used developmental gene expression in the fruit fly, Drosophila melanogaster, the bumble bee, Bombus terrestris, the harvester ant, Pogonomyrmex barbatus, and the wood ant, Formica exsecta, to compare patterns of Osiris gene expression both during development and between alternate caste phenotypes in the polymorphic social insects. Developmental gene expression of Osiris genes is highly conserved across species and correlated with gene location and evolutionary history. The social insect castes are highly divergent in pupal Osiris gene expression. Sets of co-expressed genes that include Osiris genes are enriched in gene ontology terms related to chitin/cuticle and peptidase activity. Osiris genes are essential for cuticle formation in both embryos and pupae, and genes co-expressed with Osiris genes affect wing development. Additionally, Osiris genes and those co-expressed seem to play a conserved role in insect toxicology defences and digestion. Given their role in development, plasticity, and protection, we propose that the Osiris genes play a central role in insect adaptive evolution.

Transcriptomic Analysis of Octanoic Acid Response in Drosophila sechellia Using RNA-Sequencing

The dietary specialist fruit fly Drosophila sechellia has evolved to specialize on the toxic fruit of its host plant Morinda citrifolia. Toxicity of Morinda fruit is primarily due to high levels of octanoic acid (OA). Using RNA interference (RNAi), prior work found that knockdown of Osiris family genes Osiris 6 (Osi6), Osi7, and Osi8 led to increased susceptibility to OA in adult D. melanogaster flies, likely representing genes underlying a Quantitative Trait Locus (QTL) for OA resistance in D. sechellia. While genes in this major effect locus are beginning to be revealed, prior work has shown at least five regions of the genome contribute to OA resistance. Here, we identify new candidate OA resistance genes by performing differential gene expression analysis using RNA-sequencing (RNA-seq) on control and OA-exposed D. sechellia flies. We found 104 significantly differentially expressed genes with annotated orthologs in D. melanogaster, including six Osiris gene family members, consistent with previous functional studies and gene expression analyses. Gene ontology (GO) term enrichment showed significant enrichment for cuticle development in upregulated genes and significant enrichment of immune and defense responses in downregulated genes, suggesting important aspects of the physiology of D. sechellia that may play a role in OA resistance. In addition, we identified five candidate OA resistance genes that potentially underlie QTL peaks outside of the major effect region, representing promising new candidate genes for future functional studies.

Analysis of cytochrome P450 contribution to evolved plant toxin resistance in Drosophila sechellia

Drosophila sechellia is a dietary specialist species of fruit fly that has evolved resistance to the toxic secondary defence compounds produced by the fruit of its preferred host plant Morinda citrifolia. The genetic basis of adult toxin resistance is the result of evolution at five loci across the genome. Genetic mapping between D. sechellia and Drosophila simulans and subsequent functional studies in Drosophila melanogaster have identified candidate genes potentially underlying one locus involved in toxin resistance but the remainder of the genes involved are unknown. Genes in the mixed function oxidase or cytochrome P450 gene family are frequently utilized in evolved toxin resistance in insects, yet whether they play a role in D. sechellia's resistance to the toxins found in its host plant is unknown. Here we test the role of cytochrome P450 enzymatic activity in evolved resistance to the two primary toxins found in M. citrifolia fruit: octanoic acid and hexanoic acid. We found that although cytochrome P450 enzymatic activity is involved in basal resistance it is not involved in derived toxin resistance in D. sechellia.

Osiris9a is a major component of silk fiber in lepidopteran insects

In a previous high-throughput proteomics study, it was found that the silkworm cocoon contains hundreds of complex proteins, many of which have unknown functions, in addition to fibroins, sericins, and some protease inhibitors. Osiris was one of the proteins with no known function. In this study, we identified the Osiris gene family members and constructed a phylogenetic tree based on the sequences from different species. Our results indicate that the Osiris9 gene subfamily contains six members; it is specifically expressed in lepidopteran insects and has evolved by gene duplication. An Osiris gene family member from Bombyx mori was designated as BmOsiris9a (BmOsi9a) on the basis of its homology to Drosophila melanogaster Osiris9. The expression pattern of BmOsi9a showed that it was highly expressed only in the middle silk gland of silkworm larvae, similar to Sericin1 (Ser1). BmOsi9a was visualized as two bands in western blot analysis, implying that it probably undergoes post-translational modifications. Immunohistochemistry analysis revealed that BmOsi9a was synthesized and secreted into the lumen of the middle silk gland, and was localized in the sericin layer in the silk fiber. BmOsi9a was found in the silk fibers of not only three Bombycidae species, viz. B. mori, B. mandarina, and B. huttoni, but also in the fibers collected from Saturniidae species, including Antheraea assama, Antheraea mylitta, and Samia cynthia. Although the exact biological function of Osi9a in the silk fibers is unknown, our results are important because they demonstrate that Osi9a is a common structural component of silk fiber and is expressed widely among the silk-producing Bombycidae and Saturniidae insects. Our results should help in understanding the role of Osi9a in silk fibers.

Defensive repertoire of Drosophila larvae in response to toxic fungi

Chemical warfare including insecticidal secondary metabolites is a well-known strategy for environmental microbes to monopolize a food source. Insects in turn have evolved behavioural and physiological defences to eradicate or neutralize the harmful microorganisms. We studied the defensive repertoire of insects in this interference competition by combining behavioural and developmental assays with whole-transcriptome time-series analysis. Confrontation with the toxic filamentous fungus Aspergillus nidulans severely reduced the survival of Drosophila melanogaster larvae. Nonetheless, the larvae did not behaviourally avoid the fungus, but aggregated at it. Confrontation with fungi strongly affected larval gene expression, including many genes involved in detoxification (e.g., CYP, GST and UGT genes) and the formation of the insect cuticle (e.g., Tweedle genes). The most strongly upregulated genes were several members of the insect-specific gene family Osiris, and CHK-kinase-like domains were over-represented. Immune responses were not activated, reflecting the competitive rather than pathogenic nature of the antagonistic interaction. While internal microbes are widely acknowledged as important, our study emphasizes the underappreciated role of environmental microbes as fierce competitors.