The genome of the stable fly, Stomoxys calcitrans, reveals potential mechanisms underlying reproduction, host interactions, and novel targets for pest control

Morphological analysis and lipid composition of the stable fly Stomoxys calcitrans salivary glands reinforce important features of feeding behaviour

Stomoxys calcitrans L. (Diptera: Muscidae), the stable fly, is a hematophagous insect of great veterinary importance, because it is a mechanical vector of diverse pathogens in livestock. The saliva of blood-feeding insects presents important pharmacologically active molecules that impair blood clotting, promote vasodilation and modulate the host immune system response, crucial processes for successful feeding. These properties also enable pathogens' transmission. In the present work, we describe an efficient protocol to dissect S. calcitrans salivary glands, their morphological characteristics and lipid profile. The mean length of the tubular gland is 3.23 mm with a bulbous posterior end and a narrow anterior end. Histological analysis revealed a monolayer of large polygonal epithelial cells with voluminous nuclei and high lipid content in their cytoplasm. Ultrastructural analysis showed that the epithelium is rich in mitochondria, free ribosomes, Golgi complex cisternae, presenting a great extension of rough endoplasmic reticulum that contains an electron-dense material. Lipid analysis by thin-layer chromatography showed that neutral fatty acids and phosphatidylcholine are predominant in the fly salivary glands. Lysophosphatidylcholine, an important signalling biomolecule involved in different metabolic processes, including host's immunomodulation and pathogens proliferation and differentiation, is also present.

Defensins of the Stable Fly (Stomoxys calcitrans) have Developmental-Specific Regulation and Evolve at Different Rates

Organisms produce antimicrobial peptides (AMPs) either in response to infection (induced) or continuously (constitutively) to combat microbes encountered during normal trophic activities and/or through pathogenic infections. The expression of AMPs is tightly regulated, often with specificity to particular tissues or developmental stages. As a result, AMPs face varying selective pressures based on the microbes the organism's tissue or developmental stage encounters. Here, we analyzed the evolution and developmental-specific expression of Defensins, which are ancient AMPs in insects, in the stable fly (Stomoxys calcitrans). Stable fly larvae inhibit microbe-rich environments, whereas adult flies, as blood-feeders, experience comparatively fewer encounters with diverse microbial communities. Using existing RNA-seq datasets, we identified six Defensins that were only expressed in larvae (larval Defensins) and five that were not expressed in larvae (non-larval Defensins). Each of the non-larval Defensins was expressed in at least one adult tissue sample. Half of the larval Defensins were induced by mating or feeding in adults, and all three of the induced Defensins were located downstream of canonical binding sites for an Imd transcription factor involved in the highly conserved NF-κB signaling that regulates induction of AMPs. The larval and non-larval Defensins were located in distinct genomic regions, and the amino acid sequences of the larval Defensins formed a monophyletic clade. There were more amino acid substitutions across non-larval Defensins, with multiple genes losing a highly conserved furin cleavage site thought to be required for the removal of the amino terminus from the mature Defensin domain. However, larval Defensins had a higher proportion of radical amino acid substitutions, altering amino acid size and polarity. Our results reveal insights into the developmental stage-specific regulation of AMPs, and they suggest different regulatory regimes impose unique selection pressures on AMPs, possibly as a result of variation in exposure to microbial communities across development.

Spiroplasma endosymbiont reduction of host lipid synthesis and Stomoxyn-like peptide contribute to trypanosome resistance in the tsetse fly Glossina fuscipes

Tsetse flies (Glossina spp.) vector African trypanosomes that cause devastating diseases in humans and domestic animals. Within the Glossina genus, species in the Palpalis subgroup exhibit greater resistance to trypanosome infections compared to those in the Morsitans subgroup. Varying microbiota composition and species-specific genetic traits can significantly influence the efficiency of parasite transmission. Notably, infections with the endosymbiotic bacterium Spiroplasma have been documented in several Palpalis subgroup species, including Glossina fuscipes fuscipes (Gff). While Spiroplasma infections in Gff are known to hinder trypanosome transmission, the underlying mechanisms remain unknown. To investigate Spiroplasma-mediated factors affecting Gff vector competence, we conducted high-throughput RNA sequencing of the midgut tissue along with functional assays. Our findings reveal elevated oxidative stress in the midgut environment in the presence of Spiroplasma, evidenced by increased expression of nitric oxide synthase, which catalyzes the production of trypanocidal nitric oxide. Additionally, we observed impaired lipid biosynthesis leading to a reduction of this important class of nutrients essential for parasite and host physiologies. In contrast, trypanosome infections in Gff's midgut significantly upregulated various immunity-related genes, including a small peptide, Stomoxyn-like, homologous to Stomoxyns first discovered in the stable fly Stomoxys calcitrans. We observed that the Stomoxyn-like locus is exclusive to the genomes of Palpalis subgroup tsetse species. GffStomoxyn is constitutively expressed in the cardia (proventriculus) and synthetic GffStomoxyn exhibits potent activity against Escherichia coli and bloodstream form of Trypanosoma brucei parasites, while showing no effect against insect stage procyclic forms or tsetse's commensal endosymbiont Sodalis in vitro. Reducing GffStomoxyn levels significantly increased trypanosome infection prevalence, indicating its potential trypanocidal role in vivo. Collectively, our results suggest that the enhanced resistance to trypanosomes observed in Spiroplasma-infected Gff may be due to the reduced lipid availability necessary for parasite metabolic maintenance. Furthermore, GffStomoxyn could play a crucial role in the initial immune response(s) against mammalian parasites early in the infection process in the midgut and prevent gut colonization. We discuss the molecular characteristics of GffStomoxyn, its spatial and temporal expression regulation and its microbicidal activity against Trypanosome parasites. Our findings reinforce the nutritional influences of microbiota on host physiology and host-pathogen dynamics.

Trans regulation of an odorant binding protein by a proto-Y chromosome affects male courtship in house fly

The male-limited inheritance of Y chromosomes favors alleles that increase male fitness, often at the expense of female fitness. Determining the mechanisms underlying these sexually antagonistic effects is challenging because it can require studying Y-linked alleles while they still segregate as polymorphisms. We used a Y chromosome polymorphism in the house fly, Musca domestica, to address this challenge. Two male determining Y chromosomes (YM and IIIM) segregate as stable polymorphisms in natural populations, and they differentially affect multiple traits, including male courtship performance. We identified differentially expressed genes encoding odorant binding proteins (in the Obp56h family) as candidate agents for the courtship differences. Through network analysis and allele-specific expression measurements, we identified multiple genes on the house fly IIIM chromosome that could serve as trans regulators of Obp56h gene expression. One of those genes is homologous to Drosophila melanogaster CG2120, which encodes a transcription factor that binds near Obp56h. Upregulation of CG2120 in D. melanogaster nervous tissues reduces copulation latency, consistent with this transcription factor acting as a negative regulator of Obp56h expression. The transcription factor gene, which we name speed date, demonstrates a molecular mechanism by which a Y-linked gene can evolve male-beneficial effects.

Identifying the sex chromosome and sex determination genes in the cattle tick, Rhipicephalus (Boophilus) microplus

Rhipicephalus (Boophilus) microplus is globally one of the most economically important ectoparasites of cattle costing the agriculture industry billions of dollars annually. Resistance to chemical control measures has prompted the development of novel methods of control. Recent advancements in genetic control measures for human and other animal vectors have utilized sex determination research to manipulate sex ratios, which have shown promising results in mosquitoes namely Aedes aegypti and Anopheles stephensi. Here, we use R. (B.) microplus as a model to provide foundational research to allow similar avenues of investigation in ticks using R. (B.) microplus as a model. Karyotypes for R. (B.) microplus show an XX:XO sex determining system with the largest chromosome being the sex chromosome. Using flow cytometric methods, the size of the sex chromosome was estimated at 526.91 Mb. All measures to identify the sex chromosome within the cattle tick genome assembly associated sex chromosomal characteristics to two chromosomes. This discrepancy between the assembly and karyotypes of the tick led to generating a new genome assembly with a single adult male specimen. The two chromosomes in question aligned with a single scaffold within the new genome that had a length of 513.29 Mb and was the first time the sex chromosome was identified in an Ixodid genome assembly.

Diaci v3.0: chromosome-level assembly, de novo transcriptome, and manual annotation of Diaphorina citri, insect vector of Huanglongbing

BACKGROUND

Diaphorina citri is an insect vector of "Candidatus Liberibacter asiaticus" (CLas), the gram-negative bacterial pathogen associated with citrus greening disease. Control measures rely on pesticides with negative impacts on the environment, natural ecosystems, and human and animal health. In contrast, gene-targeting methods have the potential to specifically target the vector species and/or reduce pathogen transmission.

RESULTS

To improve the genomic resources needed for targeted pest control, we assembled a D. citri genome based on PacBio long reads followed by proximity ligation-based scaffolding. The 474-Mb genome has 13 chromosomal-length scaffolds. In total, 1,036 genes were manually curated as part of a community annotation project, composed primarily of undergraduate students. We also computationally identified a total of 1,015 putative transcription factors (TFs) and were able to infer motifs for 337 TFs (33%). In addition, we produced a genome-independent transcriptome and genomes for D. citri endosymbionts.

CONCLUSIONS

Manual annotation provided more accurate gene models for use by researchers and provided an excellent training opportunity for students from multiple institutions. All resources are available on CitrusGreening.org and NCBI. The chromosomal-length D. citri genome assembly serves as a blueprint for the development of collaborative genomics projects for other medically and agriculturally significant insect vectors.

Close to complete conservation of the brachyceran opsin repertoire in the stalk-eyed fly Teleopsis dalmanni

Due to the unique morphology of their adult visual system, stalk-eyed flies represent an important model of exaggerated trait evolution through sexual selection. Early physiological measurements indicated wavelength sensitivity peaks in the ultraviolet (360 nm), blue (450), blue-green (490 nm), and red (>550 nm) ranges in the compound eye retina of the stalk-eyed fly Teleopsis dalmanni, consistent with the trichromatic color and broad range motion detection vision system of brachyceran Diptera. A previous study of dipteran opsin gene diversification, however, detected only homologs of members of the long wavelength range sensitive opsin subfamilies Rh2 and Rh6 in T. dalmanni. Here, I report findings from analyzing the most recent T. dalmanni genome assembly, which revealed the conservation of most brachyceran opsin homologs except for the UV wavelength range-sensitive homolog Rh4. These results and other examples highlight the caution that needs to be applied to gene loss conclusions.

The genome of sheep ked (Melophagus ovinus) reveals potential mechanisms underlying reproduction and narrower ecological niches

BACKGROUND

Melophagus ovinus is considered to be of great veterinary health significance. However, little is known about the information on genetic mechanisms of the specific biological characteristics and novel methods for controlling M. ovinus.

RESULTS

In total, the de novo genome assembly of M. ovinus was 188.421 Mb in size (330 scaffolds, N50 Length: 10.666 Mb), with a mean GC content of 27.74%. A total of 13,372 protein-coding genes were functionally annotated. Phylogenetic analysis indicated that the diversification of M. ovinus and Glossina fuscipes took place 72.76 Mya within the Late Cretaceous. Gene family expansion and contraction analysis revealed that M. ovinus has 65 rapidly-evolving families (26 expansion and 39 contractions) mainly involved DNA metabolic activity, transposases activity, odorant receptor 59a/67d-like, IMD domain-containing protein, and cuticle protein, etc. The universal and tightly conserved list of milk protein orthologues has been assembled from the genome of M. ovinus. Contractions and losses of sensory receptors and vision-associated Rhodopsin genes were significant in M. ovinus, which indicate that the M. ovinus has narrower ecological niches.

CONCLUSIONS

We sequenced, assembled, and annotated the whole genome sequence of M. ovinus, and launches into the preliminary genetic mechanisms analysis of the adaptive evolution characteristics of M. ovinus. These resources will provide insights to understand the biological underpinnings of this parasite and the disease control strategies.

Constitutively-expressed and induced immune effectors in the house fly (Musca domestica) and the transcription factors that may regulate them

Insects possess both infection-induced and constitutively expressed innate immune defences. Some effectors, such as lysozymes and antimicrobial peptides (AMPs), are constitutively expressed in flies, but expression patterns vary across tissues and species. The house fly (Musca domestica L.) has an impressive immune repertoire, with more effector genes than any other flies. We used RNA-seq to explore both constitutive and induced expression of immune effectors in flies. House flies were fed either Pseudomonas aeruginosa or Escherichia coli, or sterile control broth, and gene expression in the gut and carcass was analysed 4 h post-feeding. Flies fed either bacterium did not induce AMP expression, but some lysozyme and AMP genes were constitutively expressed. Prior transcriptome data from flies injected with bacteria also were analysed, and these constitutively expressed genes differed from those induced by bacterial injection. Binding sites for the transcription factor Myc were enriched upstream of constitutively expressed AMP genes, while upstream regions of induced AMPs were enriched for NF-κB binding sites resembling those of the Imd-responsive transcription factor Relish. Therefore, we identified at least two expression repertoires for AMPs in the house fly: constitutively expressed genes that may be regulated by Myc, and induced AMPs likely regulated by Relish.

Characterized Gene Repertoires and Functional Gene Reference for Forensic Entomology: Genomic and Developmental Transcriptomic Analysis of Aldrichina grahami (Diptera: Calliphoridae)

Many flies of Diptera are common entomological evidence employed in forensic investigation. Exploring the existence of inter- and intra-species genomic differences of forensically relevant insects is of great importance. Aldrichina grahami is a common blow fly species of forensic importance. The present study characterized the gene repertoires of A. grahami, and provides insights into issues related to forensic entomology, such as necrophagous behavior, gene family features, and developmental patterns. Gene families were clustered and classified according to their function in different aspects of the necrophagous lifestyle, generating several gene repertoires. The genes under positive selection pressure and evolutionary changes were screen and identified. Moreover, genes that exhibited potential prediction value in the post mortem interval (PMI) estimation and development of immature stages were subjected to analysis based on the developmental transcriptome. Related insect species were compared at the genomic level to reveal the genes associated with necrophagous behaviors. The expression of selected genes in separated repositories was verified using qPCR. This work was conducted using a high-quality chromosome-level genome assembly of A. grahami and its developmental transcriptome. The findings will facilitate future research on A. grahami and the other forensically important species.

Odorant-binding protein from the stable fly (Stomoxys calcitrans) has a high-histidine N-terminal extension that binds transition metals

The role of odorant- and pheromone-binding proteins (OBPs) in olfactory function is not fully understood. We found an OBP sequence from the stable fly, Stomoxys calcitrans, ScalOBP60, that has a 25 amino acid N-terminal extension with a high content of histidine and acidic amino acids, suggesting a possible metal binding activity. A search of public databases revealed a large number of other fly OBPs with histidine-rich N-terminal extensions, as well as beetle, wasp and ant OBPs with histidine-rich C-terminal extensions. We recombinantly expressed ScalOBP60, as well as a truncated sequence which lacks the histidine-rich N-terminal region, tScalOBP60. Using fluorescence quenching and electrospray quadrupole time-of-flight mass spectrometry (ESI-QTOF), we detected two different types of metal-binding sites. Divalent copper, nickel and zinc bind to the N-terminal histidine-rich region, and divalent copper binds to an internal sequence position. Comparison of the ESI-QTOF spectra of ScalOBP60 and tScalOBP60 showed that the histidine-rich sequence is structurally disordered, but it becomes more ordered in the presence of divalent metal. When copper is bound to the internal site, binding of a hydrophobic ligand to ScalOBP60 is inhibited. The internal and N-terminal metal sites interact allosterically, possibly through a conformational equilibrium, suggesting a mechanism for metal regulation of ligand binding to ScalOBP60. Based on our studies of ScalOBP60, we propose several possible olfactory and non-olfactory functions for this OBP.

Ensembl Genomes 2022: an expanding genome resource for non-vertebrates

Ensembl Genomes (https://www.ensemblgenomes.org) provides access to non-vertebrate genomes and analysis complementing vertebrate resources developed by the Ensembl project (https://www.ensembl.org). The two resources collectively present genome annotation through a consistent set of interfaces spanning the tree of life presenting genome sequence, annotation, variation, transcriptomic data and comparative analysis. Here, we present our largest increase in plant, metazoan and fungal genomes since the project's inception creating one of the world's most comprehensive genomic resources and describe our efforts to reduce genome redundancy in our Bacteria portal. We detail our new efforts in gene annotation, our emerging support for pangenome analysis, our efforts to accelerate data dissemination through the Ensembl Rapid Release resource and our new AlphaFold visualization. Finally, we present details of our future plans including updates on our integration with Ensembl, and how we plan to improve our support for the microbial research community. Software and data are made available without restriction via our website, online tools platform and programmatic interfaces (available under an Apache 2.0 license). Data updates are synchronised with Ensembl's release cycle.

Genome and transcriptome sequencing of the green bottle fly, Lucilia sericata, reveals underlying factors of sheep flystrike and maggot debridement therapy

The common green bottle blow fly Lucilia sericata (family, Calliphoridae) is widely used for maggot debridement therapy, which involves the application of sterile maggots to wounds. The larval excretions and secretions are important for consuming necrotic tissue and inhibiting bacterial growth in wounds of patients. Lucilia sericata is also of importance as a pest of sheep and in forensic studies to estimate a postmortem interval. Here we report the assembly of a 565.3 Mb genome from long read PacBio DNA sequencing of genomic DNA. The genome contains 14,704 predicted protein coding genes and 1709 non-coding genes. Targeted annotation and transcriptional analyses identified genes that are highly expressed in the larval salivary glands (secretions) and Malpighian tubules (excretions) under normal growth conditions and following heat stress. The genomic resources will underpin future genetic studies and in development of engineered strains for genetic control of L. sericata and for biotechnology-enhanced maggot therapy.

Genome of the pincer wasp Gonatopus flavifemur reveals unique venom evolution and a dual adaptation to parasitism and predation

Background

Hymenoptera comprise extremely diverse insect species with extensive variation in their life histories. The Dryinidae, a family of solitary wasps of Hymenoptera, have evolved innovations that allow them to hunt using venom and a pair of chelae developed from the fore legs that can grasp prey. Dryinidae larvae are also parasitoids of Auchenorrhyncha, a group including common pests such as planthoppers and leafhoppers. Both of these traits make them effective and valuable for pest control, but little is yet known about the genetic basis of its dual adaptation to parasitism and predation.

Results

We sequenced and assembled a high-quality genome of the dryinid wasp Gonatopus flavifemur, which at 636.5 Mb is larger than most hymenopterans. The expansion of transposable elements, especially DNA transposons, is a major contributor to the genome size enlargement. Our genome-wide screens reveal a number of positively selected genes and rapidly evolving proteins involved in energy production and motor activity, which may contribute to the predatory adaptation of dryinid wasp. We further show that three female-biased, reproductive-associated yellow genes, in response to the prey feeding behavior, are significantly elevated in adult females, which may facilitate the egg production. Venom is a powerful weapon for dryinid wasp during parasitism and predation. We therefore analyze the transcriptomes of venom glands and describe specific expansions in venom Idgf-like genes and neprilysin-like genes. Furthermore, we find the LWS2-opsin gene is exclusively expressed in male G. flavifemur, which may contribute to partner searching and mating.

Conclusions

Our results provide new insights into the genome evolution, predatory adaptation, venom evolution, and sex-biased genes in G. flavifemur, and present genomic resources for future in-depth comparative analyses of hymenopterans that may benefit pest control.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01081-6.