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

Osiris gene family defines the cuticle nanopatterns of Drosophila

Nanostructures of pores and protrusions in the insect cuticle modify molecular permeability and surface wetting and help insects sense various environmental cues. However, the cellular mechanisms that modify cuticle nanostructures are poorly understood. Here, we elucidate how insect-specific Osiris family genes are expressed in various cuticle-secreting cells in the Drosophila head during the early stages of cuticle secretion and cover nearly the entire surface of the head epidermis. Furthermore, we demonstrate how each sense organ cell with various cuticular nanostructures expressed a unique combination of Osiris genes. Osiris gene mutations cause various cuticle defects in the corneal nipples and pores of the chemosensory sensilla. Thus, our study emphasizes on the importance of Osiris genes for elucidating cuticle nanopatterning in insects.

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.

Osiris17 is indispensable for morphogenesis of intestinal tract in Locusta migratoria

The Osiris gene family is believed to play important roles in insect biology. Previous studies mainly focused on the roles of Osiris in Drorophila, how Osiris operates during the development of other species remains largely unknown. Here, we investigated the role of LmOsi17 in development of the hemimetabolous insect Locusta migratoria. LmOsi17 was highly expressed in the intestinal tract of nymphs. Knockdown of LmOsi17 by RNA interference (RNAi) in nymphs resulted in growth defects. The dsLmOsi17-injected nymphs did not increase in body weight or size and eventually died. Immunohistochemical analysis showed that LmOsi17 was localized to the epithelial cells of the foregut and the gastric caecum. Histological observation and hematoxylin-eosin staining indicate that the foregut and gastric caecum are deformed in dsLmOsi17 treated nymphs, suggesting that LmOsi17 is involved in morphogenesis of foregut and gastric caecum. In addition, we observed a significant reduction in the thickness of the new cuticle in dsLmOsi17-injected nymphs compared to control nymphs. Taken together, these results suggest that LmOsi17 contributes to morphogenesis of intestinal tract that affects growth and development of nymphs in locusts.

Osiris17 is essential for stable integrin localization and function during insect wing epithelia remodeling

The dynamic adhesion between cells and their extracellular matrix is essential for the development and function of organs. During insect wing development, two epithelial sheets contact each other at their basal sites through the interaction of βPS integrins with the extracellular matrix. We report that Osiris17 contributes to the maintenance of βPS integrins localization and function in developing wing of Drosophila and locust. In flies with reduced Osiris17 expression the epithelia sheets fail to maintain the integrity of basal cytoplasmic junctional bridges and basal adhesion. In contrast to the continuous basal integrin localization in control wings, this localization is disrupted during late stages of wing development in Osiris17 depleted flies. In addition, the subcellular localization revealed that Osiris17 co-localizes with the endosomal markers Rab5 and Rab11. This observation suggests an involvement of Osiris17 in endosomal recycling of integrins. Indeed, Osiris17 depletion reduced the numbers of Rab5 and Rab11 positive endosomes. Moreover, overexpression of Osiris17 increased co-localization of Rab5 and βPS integrins and partially rescued the detachment phenotype in flies with reduced βPS integrins. Taken together, our data suggest that Osiris17 is an endosome related protein that contributes to epithelial remodeling and morphogenesis by assisting basal integrins localization in insects.

Evidence of Nonrandom Patterns of Functional Chromosome Organization in Danaus plexippus

Our understanding on the interplay between gene functionality and gene arrangement at different chromosome scales relies on a few Diptera and the honeybee, species with quality reference genome assemblies, accurate gene annotations, and abundant transcriptome data. Using recently generated 'omic resources in the monarch butterfly Danaus plexippus, a species with many more and smaller chromosomes relative to Drosophila species and the honeybee, we examined the organization of genes preferentially expressed at broadly defined developmental stages (larva, pupa, adult males, and adult females) at both fine and whole-chromosome scales. We found that developmental stage-regulated genes do not form more clusters, but do form larger clusters, than expected by chance, a pattern consistent across the gene categories examined. Notably, out of the 30 chromosomes in the monarch genome, 12 of them, plus the fraction of the chromosome Z that corresponds to the ancestral Z in other Lepidoptera, were found enriched for developmental stage-regulated genes. These two levels of nonrandom gene organization are not independent as enriched chromosomes for developmental stage-regulated genes tend to harbor disproportionately large clusters of these genes. Further, although paralogous genes were overrepresented in gene clusters, their presence is not enough to explain two-thirds of the documented cases of whole-chromosome enrichment. The composition of the largest clusters often included paralogs from more than one multigene family as well as unrelated single-copy genes. Our results reveal intriguing patterns at the whole-chromosome scale in D. plexippus while shedding light on the interplay between gene expression and chromosome organization beyond Diptera and Hymenoptera.

Developmental transcriptomes predict adult social behaviours in the socially flexible sweat bee, Lasioglossum baleicum

Natural variation can provide important insights into the genetic and environmental factors that shape social behaviour and its evolution. The sweat bee, Lasioglossum baleicum, is a socially flexible bee capable of producing both solitary and eusocial nests. We demonstrate that within a single nesting aggregation, soil temperatures are a strong predictor of the social structure of nests. Sites with warmer temperatures in the spring have a higher frequency of social nests than cooler sites, perhaps because warmer temperatures provide a longer reproductive window for those nests. To identify the molecular correlates of this behavioural variation, we generated a de novo genome assembly for L. baleicum, and we used transcriptomic profiling to compare adults and developing offspring from eusocial and solitary nests. We find that adult, reproductive females have similar expression profiles regardless of social structure in the nest, but that there are strong differences between reproductive females and workers from social nests. We also find substantial differences in the transcriptomic profiles of stage-matched pupae from warmer, social-biased sites compared to cooler, solitary-biased sites. These transcriptional differences are strongly predictive of adult reproductive state, suggesting that the developmental environment may set the stage for adult behaviours in L. baleicum. Together, our results help to characterize the molecular mechanisms shaping variation in social behaviour and highlight a potential role of environmental tuning during development as a factor shaping adult behaviour and physiology in this socially flexible bee.

Single-cell transcriptomics identifies new blood cell populations in Drosophila released at the onset of metamorphosis

Drosophila blood cells called hemocytes form an efficient barrier against infections and tissue damage. During metamorphosis, hemocytes undergo tremendous changes in their shape and behavior, preparing them for tissue clearance. Yet, the diversity and functional plasticity of pupal blood cells have not been explored. Here, we combine single-cell transcriptomics and high-resolution microscopy to dissect the heterogeneity and plasticity of pupal hemocytes. We identified undifferentiated and specified hemocytes with different molecular signatures associated with distinct functions such as antimicrobial, antifungal immune defense, cell adhesion or secretion. Strikingly, we identified a highly migratory and immune-responsive pupal cell population expressing typical markers of the posterior signaling center (PSC), which is known to be an important niche in the larval lymph gland. PSC-like cells become restricted to the abdominal segments and are morphologically very distinct from typical Hemolectin (Hml)-positive plasmatocytes. G-TRACE lineage experiments further suggest that PSC-like cells can transdifferentiate to lamellocytes triggered by parasitoid wasp infestation. In summary, we present the first molecular description of pupal Drosophila blood cells, providing insights into blood cell functional diversification and plasticity during pupal metamorphosis.

Developmental transcriptomes predict adult social behaviors in the socially flexible sweat bee, Lasioglossum baleicum

Natural variation can provide important insights into the genetic and environmental factors that shape social behavior and its evolution. The sweat bee, Lasioglossum baleicum , is a socially flexible bee capable of producing both solitary and eusocial nests. We demonstrate that within a single nesting aggregation, soil temperatures are a strong predictor of the social structure of nests. Sites with warmer temperatures in the spring have a higher frequency of social nests than cooler sites, perhaps because warmer temperatures provide a longer reproductive window for those nests. To identify the molecular correlates of this behavioral variation, we generated a de novo genome assembly for L. baleicum , and we used transcriptomic profiling to compare adults and developing offspring from eusocial and solitary nests. We find that adult, reproductive females have similar expression profiles regardless of social structure in the nest, but that there are strong differences between reproductive females and workers from social nests. We also find substantial differences in the transcriptomic profiles of stage-matched pupae from warmer, social-biased sites compared to cooler, solitary-biased sites. These transcriptional differences are strongly predictive of adult reproductive state, suggesting that the developmental environment may set the stage for adult behaviors in L. baleicum . Together, our results help to characterize the molecular mechanisms shaping variation in social behavior and highlight a potential role of environmental tuning during development as a factor shaping adult behavior and physiology in this socially flexible bee.

Sexual dimorphism as a facilitator of worker caste evolution in ants

Ant societies are primarily composed of females, whereby labor is divided into reproductive, or queen, and non-reproductive, or worker, castes. Workers and reproductive queens can differ greatly in behavior, longevity, physiology, and morphology, but queen-worker differences are usually modest relative to the differences in males. Males are short-lived, typically do not provide the colony with labor, often look like a different species, and only occur seasonally. It is these differences that have historically led to their neglect in social insect research, but also why they may facilitate novel phenotypic variation - by increasing the phenotypic variability that is available for selection. In this study, worker variation in multivariate size-shape space paralleled male-queen variation. As worker variation increased within species, so did sexual variation. Across species in two independent genera, using head width as a proxy for body size, sexual size dimorphism correlated with worker polymorphism regardless of whether the ancestral condition was large or small worker/sexual dimorphism. Mounting molecular data support the hypothesis that queen-worker caste determination has co-opted many genes/pathways from sex determination. The molecular evidence, coupled with the observations from this study, leads to the hypothesis that sexual selection and selection on colony-level traits are non-independent, and that sexual dimorphism may even have facilitated the evolution of the distinct worker caste.

In the battle of the disease: a transcriptomic analysis of European foulbrood-diseased larvae of the Western honey bee (Apis mellifera)

BACKGROUND

European foulbrood is a significant bacterial brood disease of Apis sp. and can cause severe and devastating damages in beekeeping operations. Nevertheless, the epidemiology of its causative agent Melissococcus plutonius has been begun to uncover but the underlying mechanisms of infection and cause of disease still is not well understood. Here, we sought to provide insight into the infection mechanism of EFB employing RNAseq in in vitro reared Apis mellifera larvae of two developmental stages to trace transcriptional changes in the course of the disease, including Paenibacillus alvei secondary infected individuals.

RESULTS

In consideration of the progressing development of the larva, we show that infected individuals incur a shift in metabolic and structural protein-encoding genes, which are involved in metabolism of crucial compounds including all branches of macronutrient metabolism, transport protein genes and most strikingly chitin and cuticle associated genes. These changes underpin the frequently observed developmental retardation in EFB disease. Further, sets of expressed genes markedly differ in different stages of infection with almost no overlap. In an earlier stage of infection, a group of regulators of the melanization response cascade and complement component-like genes, predominantly C-type lectin genes, are up-regulated while a differential expression of immune effector genes is completely missing. In contrast, late-stage infected larvae up-regulated the expression of antimicrobial peptides, lysozymes and prominent bacteria-binding haemocyte receptor genes compared to controls. While we clearly show a significant effect of infection on expressed genes, these changes may partly result from a shift in expression timing due to developmental alterations of infection. A secondary infection with P. alvei elicits a specific response with most of the M. plutonius associated differential immune effector gene expression missing and several immune pathway genes even down-regulated.

CONCLUSION

We conclude that with progressing infection diseased individuals undergo a systemic response with a change of metabolism and their activated immune defence repertoire. Moreover, larvae are capable of adjusting their response to a secondary invasion in late stage infections.

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.

BmSuc1 Affects Silk Properties by Acting on Sericin1 in Bombyx mori

BmSuc1, a novel animal-type β-fructofuranosidase (β-FFase, EC 3.2.1.26) encoding gene, was cloned and identified for the first time in the silkworm, Bombyx mori. BmSuc1 was specifically and highly expressed in the midgut and silk gland of Bombyx mori. Until now, the function of BmSuc1 in the silk gland was unclear. In this study, it was found that the expression changes of BmSuc1 in the fifth instar silk gland were consistent with the growth rate of the silk gland. Next, with the aid of the CRISPR/Cas9 system, the BmSuc1 locus was genetically mutated, and homozygous mutant silkworm strains with truncated β-FFase (BmSUC1) proteins were established. BmSuc1 mutant larvae exhibited stunted growth and decreased body weight. Interestingly, the molecular weight of part of Sericin1 (Ser1) in the silk gland of the mutant silkworms was reduced. The knockout of BmSuc1 reduced the sericin content in the silkworm cocoon shell, and the mechanical properties of the mutant line silk fibers were also negatively affected. These results reveal that BmSUC1 is involved in the synthesis of Ser1 protein in silk glands and helps to maintain the homeostasis of silk protein content in silk fibers and the mechanical properties of silk fibers, laying a foundation for the study of BmSUC1 regulation of silk protein synthesis in silk glands.

Abnormal Antennal Olfactory Sensilla Phenotypes Involved in Olfactory Deficit in Bactrocera correcta (Diptera: Tephritidae)

Simple Summary

Tephritidae fruit flies sense odorants mainly through antennal olfactory sensilla with nanopores. Therefore, theoretically, the development of nanopore-targeted pest control technologies is an important direction in the future. Here, we report naturally occurring abnormal antennal trichoid and basiconic olfactory sensilla phenotypes consisting of abnormal bulges and reduced nanopore numbers in a long-term laboratory rearing colony of the guava fruit fly Bactrocera correcta, and further find that the reduction of nanopore numbers in these sensilla led to an olfactory deficit. Our findings provide a basis for developing nanopore-targeted pest control technologies in the future.

Abstract

The guava fruit fly, Bactrocera correcta, is one of the most destructive pests in the genus Bactrocera and detects environmental odorants mainly through antennal olfactory sensilla phenotypes with nanopores. However, it is unclear whether there are naturally occurring abnormal antennal olfactory sensilla phenotypes that affect olfaction. Here, we found that there were abnormal bulges besides nanopores on the surface of trichoid and basiconic olfactory sensilla in the antennal flagellum of long-term laboratory rearing colony (LTC), and that nanopore number in these olfactory sensilla was also remarkably reduced. Notably, the electroantennogram (EAG) responses of LTC insects to methyl eugenol or β-caryophyllene were inhibited, and their behavioral responses elicited by the same odorants were also impaired. These results revealed naturally occurring abnormal antennal olfactory sensilla phenotypes which were involved in olfactory deficit in B. correcta, providing a platform to further study nanopore-targeted pest control technologies in the future.

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.

Multiscale analysis of the randomization limits of the chromosomal gene organization between Lepidoptera and Diptera

How chromosome gene organization and gene content evolve among distantly related and structurally malleable genomes remains unresolved. This is particularly the case when considering different insect orders. We have compared the highly contiguous genome assemblies of the lepidopteran Danaus plexippus and the dipteran Drosophila melanogaster, which shared a common ancestor around 290 Ma. The gene content of 23 out of 30 D. plexippus chromosomes was significantly associated with one or two of the six chromosomal elements of the Drosophila genome, denoting common ancestry. Despite the phylogenetic distance, 9.6% of the 1-to-1 orthologues still reside within the same ancestral genome neighbourhood. Furthermore, the comparison D. plexippus-Bombyx mori indicated that the rates of chromosome repatterning are lower in Lepidoptera than in Diptera, although still within the same order of magnitude. Concordantly, 14 developmental gene clusters showed a higher tendency to retain full or partial clustering in D. plexippus, further supporting that the physical association between the SuperHox and NK clusters existed in the ancestral bilaterian. Our results illuminate the scope and limits of the evolution of the gene organization and content of the ancestral chromosomes to the Lepidoptera and Diptera while helping reconstruct portions of the genome in their most recent common ancestor.

Current understanding of immune priming phenomena in insects

It may seem that the most important issues related to insect immunity have already been described. However, novel phenomena observed in recent years shed new light on the understanding of the immune response in insects.The adaptive abilities of insects helped them to populate all ecological land niches.One important adaptive ability of insects that facilitates their success is the plasticity of their immune system. Although they only have innate immune mechanisms, insects can increase their resistance after the first encounter with the pathogen. In recent years, this phenomenon,namedimmunepriming, has become a "hot topic" in immunobiology.Priming can occur within or across generations. In the first case, the resistance of a given individual can increase after surviving a previous infection. Transstadial immune priming occurs when infection takes place at one of the initial developmental stages and increased resistance is observed at the pupal or imago stages. Priming across generations (transgenerationalimmune priming, TGIP) relies on the increased resistance of the offspring when one or both parents are infected during their lifetime.Despite the attention that immune priming has received, basic questions remain to be answered, such as regulation of immune priming at the molecular level. Research indicates that pathogen recognition receptors (PRRs) can be involved in the priming phenomenon. Recent studies have highlighted the special role of microRNAs and epigenetics, which can influence expression of genes that can be transmitted through generations although they are not encoded in the nucleotide sequence. Considerable amounts of research are required to fully understand the mechanisms that regulate priming phenomena. The aim of our work is to analyse thoroughly the most important information on immune priming in insects and help raise pertinent questions such that a greater understanding of this phenomenon can be obtained in the future.

Molecular underpinnings of the early brain developmental response to differential feeding in the honey bee Apis mellifera

Brain differential morphogenesis in females is one of the major phenotypic manifestations of caste development in honey bees. Brain diphenism appears at the fourth larval phase as a result of the differential feeding regime developing females are submitted during early phases of larval development. Here, we used a forward genetics approach to test the early brain molecular response to differential feeding leading to the brain diphenism observed at later developmental phases. Using RNA sequencing analysis, we identified 53 differentially expressed genes (DEGs) between the brains of queens and workers at the third larval phase. Since miRNAs have been suggested to play a role in caste differentiation after horizontal and vertical transmission, we tested their potential participation in regulating the DEGs. The miRNA-mRNA interaction network, including the DEGs and the royal- and worker-jelly enriched miRNA populations, revealed a subset of miRNAs potentially involved in regulating the expression of DEGs. The interaction of miR-34, miR-210, and miR-317 with Takeout, Neurotrophin-1, Forked, and Masquerade genes was experimentally confirmed using a luciferase reporter system. Taken together, our results reconstruct the regulatory network that governs the development of the early brain diphenism in honey bees.

In silico-driven analysis of the Glossina morsitans morsitans antennae transcriptome in response to repellent or attractant compounds

BACKGROUND

High-throughput sequencing generates large volumes of biological data that must be interpreted to make meaningful inference on the biological function. Problems arise due to the large number of characteristics p (dimensions) that describe each record [n] in the database. Feature selection using a subset of variables extracted from the large datasets is one of the approaches towards solving this problem.

METHODOLOGY

In this study we analyzed the transcriptome of Glossina morsitans morsitans (Tsetsefly) antennae after exposure to either a repellant (δ-nonalactone) or an attractant (ε-nonalactone). We identified 308 genes that were upregulated or downregulated due to exposure to a repellant (δ-nonalactone) or an attractant (ε-nonalactone) respectively. Weighted gene coexpression network analysis was used to cluster the genes into 12 modules and filter unconnected genes. Discretized and association rule mining was used to find association between genes thereby predicting the putative function of unannotated genes.

RESULTS AND DISCUSSION

Among the significantly expressed chemosensory genes (FDR < 0.05) in response to Ɛ-nonalactone were gustatory receptors (GrIA and Gr28b), ionotrophic receptors (Ir41a and Ir75a), odorant binding proteins (Obp99b, Obp99d, Obp59a and Obp28a) and the odorant receptor (Or67d). Several non-chemosensory genes with no assigned function in the NCBI database were co-expressed with the chemosensory genes. Exposure to a repellent (δ-nonalactone) did not show any significant change between the treatment and control samples. We generated a coexpression network with 276 edges and 130 nodes. Genes CAH3, Ahcy, Ir64a, Or67c, Ir8a and Or67a had node degree values above 11 and therefore could be regarded as the top hub genes in the network. Association rule mining showed a relation between various genes based on their appearance in the same itemsets as consequent and antecedent.

Transcriptomic, proteomic and ultrastructural studies on salinity-tolerant Aedes aegypti in the context of rising sea levels and arboviral disease epidemiology

BACKGROUND

Aedes aegypti mosquito, the principal global vector of arboviral diseases, lays eggs and undergoes larval and pupal development to become adult mosquitoes in fresh water (FW). It has recently been observed to develop in coastal brackish water (BW) habitats of up to 50% sea water, and such salinity tolerance shown to be an inheritable trait. Genomics of salinity tolerance in Ae. aegypti has not been previously studied, but it is of fundamental biological interest and important for controlling arboviral diseases in the context of rising sea levels increasing coastal ground water salinity.

RESULTS

BW- and FW-Ae. aegypti were compared by RNA-seq analysis on the gut, anal papillae and rest of the carcass in fourth instar larvae (L4), proteomics of cuticles shed when L4 metamorphose into pupae, and transmission electron microscopy of cuticles in L4 and adults. Genes for specific cuticle proteins, signalling proteins, moulting hormone-related proteins, membrane transporters, enzymes involved in cuticle metabolism, and cytochrome P450 showed different mRNA levels in BW and FW L4 tissues. The salinity-tolerant Ae. aegypti were also characterized by altered L4 cuticle proteomics and changes in cuticle ultrastructure of L4 and adults.

CONCLUSIONS

The findings provide new information on molecular and ultrastructural changes associated with salinity adaptation in FW mosquitoes. Changes in cuticles of larvae and adults of salinity-tolerant Ae. aegypti are expected to reduce the efficacy of insecticides used for controlling arboviral diseases. Expansion of coastal BW habitats and their neglect for control measures facilitates the spread of salinity-tolerant Ae. aegypti and genes for salinity tolerance. The transmission of arboviral diseases can therefore be amplified in multiple ways by salinity-tolerant Ae. aegypti and requires appropriate mitigating measures. The findings in Ae. aegypti have attendant implications for the development of salinity tolerance in other fresh water mosquito vectors and the diseases they transmit.

A Tale of Two Transcriptomic Responses in Agricultural Pests via Host Defenses and Viral Replication

Autographa californica Multiple Nucleopolyhedrovirus (AcMNPV) is a baculovirus that causes systemic infections in many arthropod pests. The specific molecular processes underlying the biocidal activity of AcMNPV on its insect hosts are largely unknown. We describe the transcriptional responses in two major pests, Spodoptera frugiperda (fall armyworm) and Trichoplusia ni (cabbage looper), to determine the host-pathogen responses during systemic infection, concurrently with the viral response to the host. We assembled species-specific transcriptomes of the hemolymph to identify host transcriptional responses during systemic infection and assessed the viral transcript abundance in infected hemolymph from both species. We found transcriptional suppression of chitin metabolism and tracheal development in infected hosts. Synergistic transcriptional support was observed to suggest suppression of immune responses and induction of oxidative stress indicating disease progression in the host. The entire AcMNPV core genome was expressed in the infected host hemolymph with a proportional high abundance detected for viral transcripts associated with replication, structure, and movement. Interestingly, several of the host genes that were targeted by AcMNPV as revealed by our study are also targets of chemical insecticides currently used commercially to control arthropod pests. Our results reveal an extensive overlap between biological processes represented by transcriptional responses in both hosts, as well as convergence on highly abundant viral genes expressed in the two hosts, providing an overview of the host-pathogen transcriptomic landscape during systemic infection.

Transcriptomic response of Anopheles gambiae sensu stricto mosquito larvae to Curry tree (Murraya koenigii) phytochemicals

BACKGROUND

Insect growth regulators (IGRs) can control insect vector populations by disrupting growth and development in juvenile stages of the vectors. We previously identified and described the curry tree (Murraya koenigii (L.) Spreng) phytochemical leaf extract composition (neplanocin A, 3-(1-naphthyl)-L-alanine, lumiflavine, terezine C, agelaspongin and murrayazolinol), which disrupted growth and development in Anopheles gambiae sensu stricto mosquito larvae by inducing morphogenetic abnormalities, reducing locomotion and delaying pupation in the mosquito. Here, we attempted to establish the transcriptional process in the larvae that underpins these phenotypes in the mosquito.

METHODS

We first exposed third-fourth instar larvae of the mosquito to the leaf extract and consequently the inherent phytochemicals (and corresponding non-exposed controls) in two independent biological replicates. We collected the larvae for our experiments sampled 24 h before peak pupation, which was 7 and 18 days post-exposure for controls and exposed larvae, respectively. The differences in duration to peak pupation were due to extract-induced growth delay in the larvae. The two study groups (exposed vs control) were consequently not age-matched. We then sequentially (i) isolated RNA (whole larvae) from each replicate treatment, (ii) sequenced the RNA on Illumina HiSeq platform, (iii) performed differential bioinformatics analyses between libraries (exposed vs control) and (iv) independently validated the transcriptome expression profiles through RT-qPCR.

RESULTS

Our analyses revealed significant induction of transcripts predominantly associated with hard cuticular proteins, juvenile hormone esterases, immunity and detoxification in the larvae samples exposed to the extract relative to the non-exposed control samples. Our analysis also revealed alteration of pathways functionally associated with putrescine metabolism and structural constituents of the cuticle in the extract-exposed larvae relative to the non-exposed control, putatively linked to the exoskeleton and immune response in the larvae. The extract-exposed larvae also appeared to have suppressed pathways functionally associated with molting, cell division and growth in the larvae. However, given the age mismatch between the extract-exposed and non-exposed larvae, we can attribute the modulation of innate immune, detoxification, cuticular and associated transcripts and pathways we observed to effects of age differences among the larvae samples (exposed vs control) and to exposures of the larvae to the extract.

CONCLUSIONS

The exposure treatment appears to disrupt cuticular development, immune response and oxidative stress pathways in Anopheles gambiae s.s larvae. These pathways can potentially be targeted in development of more efficacious curry tree phytochemical-based IGRs against An. gambiae s.s mosquito larvae.

CoFly: A gene coexpression database for the fruit fly Drosophila melanogaster

The fruit fly Drosophila melanogaster can be used as a model organism for studying various problems in biomedicine and pest management. A large number of fruit fly transcriptomes have been profiled in various cell types, tissues, development stages, toxicological exposures, and other conditions by microarray. Until now, there are still no database developed for exploring those precious data. Microarray data for 4,367 samples from National Center for Biotechnology Information Gene Expression Omnibus was collected, and analyzed by weighted gene coexpression network analysis algorithm. Fifty one gene coexpression modules that are related to cell types, tissues, development stages, and other experimental conditions were identified. The high dimensional gene expression was reduced to tens of modules that were associated with experiments/traits, representing signatures for phenotypes. Six modules were enriched with genomic regions of clustered genes. Hub genes could also be screened by intramodule connectivity. By analyzing higher order module networks, we found that cell signaling modules are more connected than other modules. Module-based gene function identification may help to discover novel gene function. An easy-to-use database was developed, which provides a new source for gene function study in the fruit fly (http://bioinformatics.fafu.edu.cn/fly/).

The nuclear and mitochondrial genomes of Frieseomelitta varia - a highly eusocial stingless bee (Meliponini) with a permanently sterile worker caste

BACKGROUND

Most of our understanding on the social behavior and genomics of bees and other social insects is centered on the Western honey bee, Apis mellifera. The genus Apis, however, is a highly derived branch comprising less than a dozen species, four of which genomically characterized. In contrast, for the equally highly eusocial, yet taxonomically and biologically more diverse Meliponini, a full genome sequence was so far available for a single Melipona species only. We present here the genome sequence of Frieseomelitta varia, a stingless bee that has, as a peculiarity, a completely sterile worker caste.

RESULTS

The assembly of 243,974,526 high quality Illumina reads resulted in a predicted assembled genome size of 275 Mb composed of 2173 scaffolds. A BUSCO analysis for the 10,526 predicted genes showed that these represent 96.6% of the expected hymenopteran orthologs. We also predicted 169,371 repetitive genomic components, 2083 putative transposable elements, and 1946 genes for non-coding RNAs, largely long non-coding RNAs. The mitochondrial genome comprises 15,144 bp, encoding 13 proteins, 22 tRNAs and 2 rRNAs. We observed considerable rearrangement in the mitochondrial gene order compared to other bees. For an in-depth analysis of genes related to social biology, we manually checked the annotations for 533 automatically predicted gene models, including 127 genes related to reproductive processes, 104 to development, and 174 immunity-related genes. We also performed specific searches for genes containing transcription factor domains and genes related to neurogenesis and chemosensory communication.

CONCLUSIONS

The total genome size for F. varia is similar to the sequenced genomes of other bees. Using specific prediction methods, we identified a large number of repetitive genome components and long non-coding RNAs, which could provide the molecular basis for gene regulatory plasticity, including worker reproduction. The remarkable reshuffling in gene order in the mitochondrial genome suggests that stingless bees may be a hotspot for mtDNA evolution. Hence, while being just the second stingless bee genome sequenced, we expect that subsequent targeting of a selected set of species from this diverse clade of highly eusocial bees will reveal relevant evolutionary signals and trends related to eusociality in these important pollinators.

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.

Development of the transcriptome for a sediment ecotoxicological model species, Chironomus dilutus

Chironomus dilutus is a prominent model species in conventional sediment toxicity testing and sediment contamination diagnosis. However, lack of genomic data significantly limited its application in identifying toxicological mode of action (MOA) and molecular biomarkers of toxicants. Here the transcriptome of C. dilutus in full life span and both sexes (1st, 2nd, 3rd and 4th instar larvae, pupae, and adults) were developed and temporal gene expression across adjacent life stages were investigated to understand the regulation of development. Furthermore, transcriptional response of Midges (the 4th instar larvae) exposed to chemicals of different MOAs (CdCl₂, nonylphenol and triclosan) were profiled based on the reference transcriptome. Consequently, a complete transcriptome of 31132 unigenes with N50 of 3117bp, covering 98.8% of the arthropod single-copy orthologs were assembled. While 364 genes were differentially expressed among adjacent larval stages, 7142 and 2127 of transcripts were significantly changed for the transition of larvae-pupae and pupae-adults, respectively. Finally, chemical-specific gene expression profile were identified in the midges, showed its potential in classifying distinct contaminants. Overall, the comprehensive transcriptome of C. dilutus developed here could not only facilitate the mechanistic understanding of environmental toxicants during critical life stage of aquatic insects, but also provide molecular diagnostic tools in sediment ecotoxicology.

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.

Improved Genome Assembly and Annotation of the Soybean Aphid (Aphis glycines Matsumura)

Aphids are an economically important insect group due to their role as plant disease vectors. Despite this economic impact, genomic resources have only been generated for a small number of aphid species. The soybean aphid (Aphis glycines Matsumura) was the third aphid species to have its genome sequenced and the first to use long-read sequence data. However, version 1 of the soybean aphid genome assembly has low contiguity (contig N50 = 57 Kb, scaffold N50 = 174 Kb), poor representation of conserved genes and the presence of genomic scaffolds likely derived from parasitoid wasp contamination. Here, I use recently developed methods to reassemble the soybean aphid genome. The version 2 genome assembly is highly contiguous, containing half of the genome in only 40 scaffolds (contig N50 = 2.00 Mb, scaffold N50 = 2.51 Mb) and contains 11% more conserved single-copy arthropod genes than version 1. To demonstrate the utility of this improved assembly, I identify a region of conserved synteny between aphids and Drosophila containing members of the Osiris gene family that was split over multiple scaffolds in the original assembly. The improved genome assembly and annotation of A. glycines demonstrates the benefit of applying new methods to old data sets and will provide a useful resource for future comparative genome analysis of aphids.

Molecular Insights into the Insensitivity of Lepidopteran Pests to Cycloxaprid

Cycloxaprid (CYC) is effective in the control of hemipteran pests, but its bioactivity against lepidopteran pests is still unclear. Here, the bioactivity of CYC against lepidopteran pests was found to be much worse than that against hemipteran insects. To reveal the mechanism, the transcriptomes of CYC-treated and untreated Ostrinia furnacalis larvae were compared. Among the top 20 differentially expressed genes, 11 encode proteins involved in cuticle formation, while only one encodes a detoxifying enzyme. Thus, the cuticle appears to be important for the insensitivity of O. furnacalis to CYC. A pretreatment of O. furnacalis larvae with methoprene enhanced the bioactivity of CYC by 1.12-fold. Moreover, mixtures of CYC with graphene oxide increased the bioactivity of CYC by 1.88-fold. Because lepidopteran and hemipteran insects often harm crops at the same time, the work can help make full use of CYC and reduce the environmental impacts of using multiple pesticides.

How interactions with plant chemicals shape insect genomes

The transition to herbivory by insects is associated with distinct genomic signatures. Sequenced genomes of extant herbivore species reveal the result of these transitions, but in lieu of comparisons between herbivorous and non-herbivorous lineages that diverged recently, such datasets have shed less light on the evolutionary genomic processes involved in diet shifts to or from herbivory. Here, we propose that the comparative genomics of diet shifts between closely related insect herbivores and non-herbivores, and within densely-sampled clades of herbivores, will help reveal the extent to which herbivory evolves through the co-option and subtle remodeling of widely-conserved gene families with functions ancestrally distinct from phytophagy.

Effect of heritable symbionts on maternally-derived embryo transcripts

Maternally-transmitted endosymbiotic bacteria are ubiquitous in insects. Among other influential phenotypes, many heritable symbionts of arthropods are notorious for manipulating host reproduction through one of four reproductive syndromes, which are generally exerted during early developmental stages of the host: male feminization; parthenogenesis induction; male killing; and cytoplasmic incompatibility (CI). Major advances have been achieved in understanding mechanisms and identifying symbiont factors involved in reproductive manipulation, particularly male killing and cytoplasmic incompatibility. Nonetheless, whether cytoplasmically-transmitted bacteria influence the maternally-loaded components of the egg or early embryo has not been examined. In the present study, we investigated whether heritable endosymbionts that cause different reproductive phenotypes in Drosophila melanogaster influence the mRNA transcriptome of early embryos. We used mRNA-seq to evaluate differential expression in Drosophila embryos lacking endosymbionts (control) to those harbouring the male-killing Spiroplasma poulsonii strain MSRO-Br, the CI-inducing Wolbachia strain wMel, or Spiroplasma poulsonii strain Hyd1; a strain that lacks a reproductive phenotype and is naturally associated with Drosophila hydei. We found no consistent evidence of influence of symbiont on mRNA composition of early embryos, suggesting that the reproductive manipulation mechanism does not involve alteration of maternally-loaded transcripts. In addition, we capitalized on several available mRNA-seq datasets derived from Spiroplasma-infected Drosophila melanogaster embryos, to search for signals of depurination of rRNA, consistent with the activity of Ribosome Inactivating Proteins (RIPs) encoded by Spiroplasma poulsonii. We found small but statistically significant signals of depurination of Drosophila rRNA in the Spiroplasma treatments (both strains), but not in the symbiont-free control or Wolbachia treatment, consistent with the action of RIPs. The depurination signal was slightly stronger in the treatment with the male-killing strain. This result supports a recent report that RIP-induced damage contributes to male embryo death.

Nanopore Formation in the Cuticle of an Insect Olfactory Sensillum

Nanometer-level patterned surface structures form the basis of biological functions, including superhydrophobicity, structural coloration, and light absorption [1-3]. In insects, the cuticle overlying the olfactory sensilla has multiple small (50- to 200-nm diameter) pores [4-8], which are supposed to function as a filter that admits odorant molecules, while preventing the entry of larger airborne particles and limiting water loss. However, the cellular processes underlying the patterning of extracellular matrices into functional nano-structures remain unknown. Here, we show that cuticular nanopores in Drosophila olfactory sensilla originate from a curved ultrathin film that is formed in the outermost envelope layer of the cuticle and secreted from specialized protrusions in the plasma membrane of the hair forming (trichogen) cell. The envelope curvature coincides with plasma membrane undulations associated with endocytic structures. The gore-tex/Osiris23 gene encodes an endosomal protein that is essential for envelope curvature, nanopore formation, and odor receptivity and is expressed specifically in developing olfactory trichogen cells. The 24-member Osiris gene family is expressed in cuticle-secreting cells and is found only in insect genomes. These results reveal an essential requirement for nanopores for odor reception and identify Osiris genes as a platform for investigating the evolution of surface nano-fabrication in insects.

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.