Chromosome-level de novo genome assembly of Sarcophaga peregrina provides insights into the evolutionary adaptation of flesh flies

The gut bacterial composition across life stages of Sarcophaga peregrina (Diptera: Sarcophagidae) and the effects of amikacin on their development

Insects and microorganisms, ubiquitous organisms in the natural world, have developed intricate relationships throughout their evolutionary histories. However, most studies have concentrated on specific time points or life stages, but some limited studies have investigated the dynamics of microbial diversity within insects across life stages. Here, 16S rDNA sequencing technology was used to investigate the gut bacterial community across the life stages of Sarcophaga peregrina (Robineau-Desvoidy) (Diptera: Sarcophagidae). The results revealed that the gut bacterial diversity of S. peregrina varied with life stage and showed similarity in the nearby life stages. Proteobacteria, Actinobacteria, Firmicutes, and Bacteroidetes were the dominant phyla in S. peregrina. Genera such as Providencia, Ignatzschineria, and Myroides are implicated in potentially pivotal roles during the developmental processes of this flesh fly. Furthermore, the effects of amikacin on the growth and development of S. peregrina were not statistically significant. However, we did observe significant changes at the protein level, which suggests a close association between protein-level alterations and growth and development. Additionally, we speculate that S. peregrina regulates its nutritional status during nonfeeding stages to meet the demands of eclosion. This study represents the first comprehensive examination of the intestinal bacterial composition across various life stages of S. peregrina. Our findings deepen our understanding of the gut microbiota in this flesh fly and lay the groundwork for further exploration into the intricate interactions between microorganisms and insects.

HOFE: an interactive forensic entomological database

The significance of entomological evidence in inferring the time, location and cause of death has been demonstrated both theoretically and practically. With the advancement of sequencing technologies, reports have emerged on necrophagous insects' nuclear genomes, transcriptomes, proteomes and mitochondrial genomes. However, within the field of forensic entomology, there is currently no available database that can integrate, store and share the resources of necrophagous insects. The absence of a database poses an inconvenience to the application of entomological evidence in judicial practice and hampers the development of the forensic entomology discipline. Given this, we have developed the Home Of Forensic Entomology database, encompassing 10 core functional modules: Home, Browse, Mitochondria, Proteome, JBrowse, Search, BLAST, Tools, Case base and Maps. Notably, the 'Tools' module enables multiple sequence alignment analysis (Muscle), homologous protein prediction (Genewise), primer design (Primer), large-scale genomic analysis (Lastz), Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis, as well as expression profiling (PCA Analysis, Hcluster and Correlation Heatmap). In addition, the present database also works as an interactive platform for researchers by sharing forensic entomological case reports and uploading data and material. This database provides potential visitors with a comprehensive function for multi-omics data analysis, offers substantial references to researchers and criminal scene investigators and facilitates the utilization of entomological evidence in court. Database URL: http://ihofe.com/.

New insights into immune genes and other expanded gene families of the house fly, Musca domestica, from an improved whole genome sequence

The house fly, Musca domestica, is a pest of livestock, transmits pathogens of human diseases, and is a model organism in multiple biological research areas. The first house fly genome assembly was published in 2014 and has been of tremendous use to the community of house fly biologists, but that genome is discontiguous and incomplete by contemporary standards. To improve the house fly reference genome, we sequenced, assembled, and annotated the house fly genome using improved techniques and technologies that were not available at the time of the original genome sequencing project. The new genome assembly is substantially more contiguous and complete than the previous genome. The new genome assembly has a scaffold N50 of 12.46 Mb, which is a 50-fold improvement over the previous assembly. In addition, the new genome assembly is within 1% of the estimated genome size based on flow cytometry, whereas the previous assembly was missing nearly one-third of the predicted genome sequence. The improved genome assembly has much more contiguous scaffolds containing large gene families. To provide an example of the benefit of the new genome, we used it to investigate tandemly arrayed immune gene families. The new contiguous assembly of these loci provides a clearer picture of the regulation of the expression of immune genes, and it leads to new insights into the selection pressures that shape their evolution.

Effects of heat tolerance on the gut microbiota of Sarcophaga peregrina (Diptera: Sarcophagidae) and impacts on the life history traits

BACKGROUND

Heat tolerance is a distinct abiotic factor affecting the distribution and abundance of insects. Gut microbiota can contribute to host fitness, thereby increasing resistance to abiotic stress conditions. In this study, Sarcophaga peregrina is closely associated with human life in ecological habits and shows remarkable adaptability to daily and seasonal temperature fluctuations. To date, the role of gut microbiota in S. peregrina response to heat stress and its influence on the host phenotypic variability remain poorly studied.

METHODS

We exposed S. peregrina to heat stress at 40 °C for 3 h every day throughout the developmental stages from newly hatched larva to adult, after which gut DNA was extracted from third-instar larvae, early pupal stage, late pupal stage, and newly emerged adults, respectively. Then, 16S rRNA microbial community analyses were performed.

RESULTS

Firstly, we analyzed whether heat stress could have an impact on the life history traits of S. peregrina and showed that the growth rate of larvae was higher and the developmental time was significantly shorter after heat stress. We then proposed the role of the gut microbiota in the heat tolerance of S. peregrina, which indicated that the bacterial abundance and community structure changed significantly after heat tolerance. In particular, the relative abundance of Wohlfahrtiimonas and Ignatzschineria was higher in the third-instar larval larvae; the former increased and the latter decreased significantly after heat stress. To further explore the effect of disturbing the microbial community on thermotolerant phenotype, newly hatched larvae were fed with amikacin under heat stress, which indicated that the larval length and the whole developmental cycle was significantly shorter.

CONCLUSION

This study indicated that Wohlfahrtiimonas and Ignatzschineria should play an important role in the post-feeding stage under heat stress, but further study is still needed. In general, heat tolerance can affect the gut microbial community structure, which in turn affects the fitness of the host.

Application of omics techniques in forensic entomology research

With the advent of the post-genome era, omics technologies have developed rapidly and are widely used, including in genomics, transcriptomics, proteomics, metabolomics, and microbiome research. These omics techniques are often based on comprehensive and systematic analysis of biological samples using high-throughput analysis methods and bioinformatics, to provide new insights into biological phenomena. Currently, omics techniques are gradually being applied to forensic entomology research and are useful in species identification, phylogenetics, screening for developmentally relevant differentially expressed genes, and the interpretation of behavioral characteristics of forensic-related species at the genetic level. These all provide valuable information for estimating the postmortem interval (PMI). This review mainly discusses the available omics techniques, summarizes the application of omics techniques in forensic entomology, and their future in the field.

Genome-wide analysis of long noncoding RNAs and their association in regulating the metamorphosis of the Sarcophaga peregrina (Diptera: Sarcophagidae)

BACKGROUND

The flesh fly, Sarcophaga peregrina (Diptera: Sarcophagidae), is an important hygiene pest, that causes myiasis in humans and other mammals, typically livestock, and as a vector for various parasitic agents, including bacteria, viruses, and parasites. The role of long non-coding RNAs (lncRNAs) in regulating gene expression during metamorphosis of the flesh fly has not been well established.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we performed genome-wide identification and characterization of lncRNAs from the early pupal stage (1-days pupae), mid-term pupal stage (5-days pupae), and late pupal stage (9-days pupae) of S. peregrina by RNA-seq, and a total of 6921 lncRNAs transcripts were identified. RT-qPCR and enrichment analyses revealed the differentially expressed lncRNAs (DE lncRNAs) that might be associated with insect metamorphosis development. Furthermore, functional analysis revealed that the DE lncRNA (SP_lnc5000) could potentially be involved in regulating the metamorphosis of S. peregrina. RNA interference of SP_lnc5000 caused reduced expression of metamorphosis-related genes in 20-hydroxyecdysone (20E) signaling (Br-c, Ftz-F1), cuticle tanning pathway (TH, DOPA), and chitin related pathway (Cht5). Injection of dsSP_lnc5000 in 3rd instar larvae of S. peregrina resulted in deformed pupae, stagnation of pupal-adult metamorphosis, and a decrease in development time of pupal, pupariation rates and eclosion rates. Hematoxylin-eosin staining (H&E), scanning electron microscope (SEM) observation and cuticle hydrocarbons (CHCs) analysis indicated that SP_lnc5000 had crucial roles in the metamorphosis developmental by modulating pupal cuticular development.

CONCLUSIONS/SIGNIFICANCE

We established that the lncRNA SP_lnc5000 potentially regulates the metamorphosis of S. peregrina by putatively affecting the structure and composition of the pupal cuticle. This study enhances our understanding of lncRNAs as regulators of metamorphosis in S. peregrina, and provide valuable insights into the identification of potential targets for vector control and the development of effective strategies for controlling the spread of myiasis and parasitic diseases.

Development of Forensically Important Sarcophaga peregrina (Diptera: Sarcophagidae) and Intra-Puparial Age Estimation Utilizing Multiple Methods at Constant and Fluctuating Temperatures

Sarcophaga peregrina (Robineau-Desvoidy, 1830) has the potential to estimate the minimum postmortem interval (PMI_min). Development data and intra-puparial age estimation are significant for PMI_min estimation. Previous research has focused on constant temperatures, although fluctuating temperatures are a more real scenario at a crime scene. The current study examined the growth patterns of S. peregrina under constant (25.75 °C) and fluctuating temperatures (18-36 °C; 22-30 °C). Furthermore, differentially expressed genes, attenuated total reflectance Fourier-transform infrared spectroscopy, and cuticular hydrocarbons of S. peregrina during the intra-puparial period were used to estimate age. The results indicated that S. peregrina at fluctuating temperatures took longer to develop and had a lower pupariation rate, eclosion rate, and pupal weight than the group at constant temperatures did. Moreover, we found that six DEG expression profiles and ATR-FTIR technology, CHCs detection methods, and chemometrics can potentially estimate the intra-puparial age of S. peregrina at both constant and fluctuating temperatures. The findings of the study support the use of S. peregrina for PMI_min estimation and encourage the use of entomological evidence in forensic practice.

Chromosome-level genome assembly of Dastarcus helophoroides provides insights into CYP450 genes expression upon insecticide exposure

BACKGROUND

Dastarcus helophoroides is an important natural enemy of cerambycids, and is wildly used in biological control of pests. Nevertheless, the absence of complete genomic information limits the investigation of the underlying molecular mechanisms. Here, a chromosome-level of Dastarcus helophoroides genome is assembled using a combination strategy of Illumina, PacBio, 10x™ Genomics, and Hi-C.

RESULTS

The final assembly is 609.09 Mb with contig N50, scaffold N50 and GC content of 5.46 Mb, 42.56 Mb and 31.50%, respectively, and 95.25% of the contigs anchor into 13 chromosomes. In total 14 890 protein-coding genes and 65.37% repeat sequences are predicted in the assembly genome. The phylogenetic analysis of single-copy gene families shared among 20 insect species indicates that Dastarcus helophoroides is placed as the sister species to clade (Nitidulidae+Curculionoidea+Chrysomeloidea) + Tenebrionoidea, and diverges from the related species ~242.9 Mya. In total 36 expanded gene families are identified in Dastarcus helophoroides genome, and are functionally related to drug metabolism and metabolism of xenobiotics by cytochrome P450. Some members of CYP4 Clade and CYP6 Clade are up-regulated in Dastarcus helophoroides adults upon insecticide exposure, of which expressions of DhCYP4Q, DhCYP6A14X1 and DhCYP4C1 are significantly up-regulated. The silencing of the three genes leads to adults more sensitive to insecticide and increased knocked-down rate, which may indicate their critical roles in stress resistance and detoxication.

CONCLUSION

Our study systematically integrated the chromosome-level genome, transcriptome and gene expression of Dastarcus helophoroides, which will provide valuable resources for understanding mechanisms of pesticide metabolism, growth and development, and utilization of the natural enemy in integrated control. © 2022 Society of Chemical Industry.

Multimethod combination for age estimation of Sarcophaga peregrina (Diptera: Sarcophagidae) with implications for estimation of the postmortem interval

Sarcophaga peregrina (Robineau-Desvoidy, 1830) (Diptera: Sarcophagidae) is a forensically important flesh fly with potential value for estimating the minimum postmortem interval (PMI_min). Here, the developmental patterns of S. peregrina were investigated at 5 constant temperatures (15-35 °C). Morphological changes at different developmental stages of the pupa were observed at 4 constant temperatures (15-30 °C) by removing the puparium and staining the pupa with hematoxylin and eosin. Furthermore, differentially expressed genes (DEGs) were analyzed at 25 °C in the intrapuparial period to estimate the age of S. peregrina during the intrapuparial stage. S. peregrina completed development from larviposition to adult eclosion at 15 °C, 20 °C, 25 °C, and 30 °C; the developmental durations were 1090.3 ± 30.6 h, 566.6 ± 21.9 h, 404.6 ± 13.01 h, and 280.3 ± 4.5 h, respectively, while the development could not be completed at 35 °C. The intrapuparial period of S. peregrina was divided into 12 sub-stages on the basis of the overall external morphological changes; 6 sub-stages on the basis of individual morphological structures such as the compound eyes, antennae, thorax, legs, wings, and abdomen; and 10 sub-stages on the basis of internal morphological changes detected using histological analysis. The period of each sub-stage or structure that appeared was determined. Moreover, we found that 6 genes (NDUFS2, CPAMD8, NDUFV2, Hsp27, Hsp23, and TPP) with differential expression can be used for the precise age estimation of S. peregrina during the intrapuparial period. This study provided basic developmental data for the use of S. peregrina in PMI_min estimation, and we successfully estimated PMI_min in a real forensic case by using a multimethod combination.

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.

Chromosome-level genome assembly reveals potential epigenetic mechanisms of the thermal tolerance in the oriental fruit fly, Bactrocera dorsalis

The oriental fruit fly, Bactrocera dorsalis (Hendel), has very strong ecological adaptability and phenotypic plasticity. Here, the genome of B. dorsalis was assembled into 549.45 Mb sequences with a contig N50 length of 12.81 Mb. Among, 95.67 % assembled genome sequences were anchored on six chromosomes with an N50 length of 94.63 Mb. According to the basic characteristics of the sex chromosomes of Tephritidae, the X chromosome of B. dorsalis was identified. Significant gene expansions were detected in several important gene families related to adaptability. In particular, we annotated 50 histone modification enzymes (HMEs) in this genome. A comparative transcriptome analysis indicated that 12 HME genes were differentially expressed in two thermo-tolerant strains (heat and cold). Interestingly, four and seven of the 12 HME genes responded to heat shock or cold hardening, respectively. These evidences suggested that the histone modification as an epigenetic modification may be involved in the thermal tolerance of B. dorsalis, but with different regulation mechanisms in thermal acclimation and hardening. The high quality genome of B. dorsalis provides an invaluable resource for further functional genomic study. Moreover, comparative genomic analysis will shed insights on revealing the mechanisms of adaptive evolution in this fly.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

Predicting the Weathering Time by the Empty Puparium of Sarcophaga peregrina (Diptera: Sarcophagidae) with the ANN Models

Empty puparium are frequently collected at crime scenes and may provide valuable evidence in cases with a long postmortem interval (PMI). Here, we collected the puparium of Sarcophaga peregrina (Diptera: Sarcophagidae) (Robineau-Desvoidy, 1830) for 120 days at three temperatures (10 °C, 25 °C, and 40 °C) with the aim to estimate the weathering time of empty puparium. The CHC profiles were analyzed by gas chromatography-mass spectrometry (GC-MS). The partial least squares (PLS), support vector regression (SVR), and artificial neural network (ANN) models were used to estimate the weathering time. This identified 49 CHCs with a carbon chain length between 10 and 33 in empty puparium. The three models demonstrate that the variation tendency of hydrocarbon could be used to estimate the weathering time, while the ANN models show the best predictive ability among these three models. This work indicated that puparial hydrocarbon weathering has certain regularity with weathering time and can gain insight into estimating PMI in forensic investigations.

Temporal Expression Profiles Reveal Potential Targets during Postembryonic Development of Forensically Important Sarcophaga peregrina (Diptera: Sarcophagidae)

Sarcophaga peregrina (Robineau-Desvoidy, 1830) is a species of medical and forensic importance. In order to investigate the molecular mechanism during postembryonic development and identify specific genes that may serve as potential targets, transcriptome analysis was used to investigate its gene expression dynamics from the larval to pupal stages, based on our previous de novo-assembled genome of S. peregrina. Totals of 2457, 3656, 3764, and 2554 differentially expressed genes were identified. The specific genes encoding the structural constituent of cuticle were significantly differentially expressed, suggesting that degradation and synthesis of cuticle-related proteins might actively occur during metamorphosis. Molting (20-hydroxyecdysone, 20E) and juvenile (JH) hormone pathways were significantly enriched, and gene expression levels changed in a dynamic pattern during the developmental stages. In addition, the genes in the oxidative phosphorylation pathway were significantly expressed at a high level during the larval stage, and down-regulated from the wandering to pupal stages. Weighted gene co-expression correlation network analysis (WGCNA) further demonstrated the potential regulation mechanism of tyrosine metabolism in the process of puparium tanning. Moreover, 10 consistently up-regulated genes were further validated by qRT-PCR. The utility of the models was then examined in a blind study, indicating the ability to predict larval development. The developmental, stage-specific gene profiles suggest novel molecular markers for age prediction of forensically important flies.

A Study of Cuticular Hydrocarbons of All Life Stages in Sarcophaga peregrina (Diptera: Sarcophagidae)

Sarcophaga peregrina (Robineau-Desvoidy, 1830), a synanthropic flesh fly species found in different parts of the world, is of medical and forensic importance. Traditional methods of inferring developmental age rely on the life stage of insects and morphological changes. However, once the larvae reach the pupal and adult stage, morphological changes would become barely visible, so that the classic method would be invalid. Here, we studied the cuticular hydrocarbon profile of S. peregrina of the whole life cycle from larval stage to adult stage by GC-MS. Sixty-three compounds with carbon chain length ranging from 8 to 36 were detected, which could be categorized into four classes: n-alkanes, branched alkanes, alkenes, and unknowns. As developmental increased, branched alkanes dominant, and the content of high-molecular-weight hydrocarbons is variable, especially for 2-methyl C19, DiMethyl C21, docosane (C22), and tricosane (C23). This study shows that the composition of CHC could be used to determine the developmental age of S. peregrina and aid in postmortem interval estimations in forensic science.