A blow to the fly — Lucilia cuprina draft genome and transcriptome to support advances in biology and biotechnology

Age determination of Chrysomya megacephala (Diptera: Calliphoridae) using lifespan patterns, gene expression, and pteridine concentration under constant and variable temperatures

Chrysomya megacephala (Fabricius, 1794) (Diptera: Calliphoridae), is a blowfly species widely studied in medical, veterinary, and entomological research. Our study examined the impact of constant (15, 20, 25, 30, and 35 °C) and variable (ranging from 21.0 to 25.4 °C, with an average of 23.31 °C) temperatures on the development and larval body length of C. megacephala. Additionally, we analyzed the age of the adult C. megacephala through pteridine content and related metabolic genes analysis. Our findings revealed three distinct growth patterns: isomorphen diagram, isomegalen diagram, and thermal accumulated models. At constant temperatures of 15, 20, 25, 30, and 35 °C, egg-hatching times were 44.5 ± 8.9, 26.7 ± 4.6, 12.6 ± 1.1, 11.0 ± 1.0, and 9.9 ± 1.9 h, respectively, while it was 15.3 ± 5.9 h at variable temperatures. The total development times from oviposition to adult eclosion in C. megacephala required 858.1 ± 69.2, 362.3 ± 5.9, 289.6 ± 17.8, 207.3 ± 9.3, and 184.7 ± 12.1 h at constant temperatures of 15, 20, 25, 30, and 35 °C, respectively. This duration was extended to 282.0 ± 64.1 h under variable temperatures. However, no significant differences were found in hatching times and the total developmental durations between 25 °C and variable temperatures. A developmental threshold temperature (D0) of 9.90 ± 0.77 °C and a thermal summation constant (K) of 4244.0 ± 347.0° hours were ascertained. Pteridine content patterns varied significantly across constant temperatures, but not between 25 °C and variable temperatures. Sex and temperature were identified as the primary factors influencing pteridine levels in the head of C. megacephala. Gene expression associated with pteridine metabolism decreased following adult eclosion, matching with increased pteridine concentration. Further investigations are needed to explore the use of pteridine cofactors for age-grading adult necrophagous flies. These findings provide valuable insights into the lifespan of C. megacephala, thereby offering valuable groundwork for forthcoming investigations and PMImin determination.

Harnessing model organism genomics to underpin the machine learning-based prediction of essential genes in eukaryotes - Biotechnological implications

The availability of high-quality genomes and advances in functional genomics have enabled large-scale studies of essential genes in model eukaryotes, including the 'elegant worm' (Caenorhabditis elegans; Nematoda) and the 'vinegar fly' (Drosophila melanogaster; Arthropoda). However, this is not the case for other, much less-studied organisms, such as socioeconomically important parasites, for which functional genomic platforms usually do not exist. Thus, there is a need to develop innovative techniques or approaches for the prediction, identification and investigation of essential genes. A key approach that could enable the prediction of such genes is machine learning (ML). Here, we undertake an historical review of experimental and computational approaches employed for the characterisation of essential genes in eukaryotes, with a particular focus on model ecdysozoans (C. elegans and D. melanogaster), and discuss the possible applicability of ML-approaches to organisms such as socioeconomically important parasites. We highlight some recent results showing that high-performance ML, combined with feature engineering, allows a reliable prediction of essential genes from extensive, publicly available 'omic data sets, with major potential to prioritise such genes (with statistical confidence) for subsequent functional genomic validation. These findings could 'open the door' to fundamental and applied research areas. Evidence of some commonality in the essential gene-complement between these two organisms indicates that an ML-engineering approach could find broader applicability to ecdysozoans such as parasitic nematodes or arthropods, provided that suitably large and informative data sets become/are available for proper feature engineering, and for the robust training and validation of algorithms. This area warrants detailed exploration to, for example, facilitate the identification and characterisation of essential molecules as novel targets for drugs and vaccines against parasitic diseases. This focus is particularly important, given the substantial impact that such diseases have worldwide, and the current challenges associated with their prevention and control and with drug resistance in parasite populations.

Phylogenomics of Opsin Genes in Diptera Reveals Lineage-Specific Events and Contrasting Evolutionary Dynamics in Anopheles and Drosophila

Diptera is one of the biggest insect orders and displays a large diversity of visual adaptations. Similarly to other animals, the dipteran visual process is mediated by opsin genes. Although the diversity and function of these genes are well studied in key model species, a comprehensive comparative genomic study across the dipteran phylogeny is missing. Here we mined the genomes of 61 dipteran species, reconstructed the evolutionary affinities of 528 opsin genes, and determined the selective pressure acting in different species. We found that opsins underwent several lineage-specific events, including an independent expansion of Long Wave Sensitive opsins in flies and mosquitoes, and numerous family-specific duplications and losses. Both the Drosophila and the Anopheles complement are derived in comparison with the ancestral dipteran state. Molecular evolutionary studies suggest that gene turnover rate, overall mutation rate, and site-specific selective pressure are higher in Anopheles than in Drosophila. Overall, our findings indicate an extremely variable pattern of opsin evolution in dipterans, showcasing how two similarly aged radiations, Anopheles and Drosophila, are characterized by contrasting dynamics in the evolution of this gene family. These results provide a foundation for future studies on the dipteran visual system.

Time Flies-Age Grading of Adult Flies for the Estimation of the Post-Mortem Interval

The estimation of the minimum time since death is one of the main applications of forensic entomology. This can be done by calculating the age of the immature stage of necrophagous flies developing on the corpse, which is confined to approximately 2–4 weeks, depending on temperature and species of the first colonizing wave of flies. Adding the age of the adult flies developed on the dead body could extend this time frame up to several weeks when the body is in a building or closed premise. However, the techniques for accurately estimating the age of adult flies are still in their beginning stages or not sufficiently validated. Here we review the current state of the art of analysing the aging of flies by evaluating the ovarian development, the amount of pteridine in the eyes, the degree of wing damage, the modification of their cuticular hydrocarbon patterns, and the increasing number of growth layers in the cuticula. New approaches, including the use of age specific molecular profiles based on the levels of gene and protein expression and the application of near infrared spectroscopy, are introduced, and the forensic relevance of these methods is discussed.

Evaluation of Reference Genes and Age Estimation of Forensically Useful Aldrichina grahami (Diptera: Calliphoridae) During Intrapuparial Period

The minimum postmortem interval (PMImin) could be evaluated from the developmental stage of forensically important insects colonize a corpse, such as blow flies (Diptera: Calliphoridae). Unlike larvae, the developmental stage of which is well established according to their morphology, estimating the age of pupae is proven to be challenging. Recently, several studies reported the regulation of special genes during the development of blow fly pupae. However, gene regulation in Aldrichina grahami during the intrapuparial period remains to be studied. Therefore, we set out to investigate the mRNA levels of heat shock protein 23 (Hsp23), heat shock protein 24 (Hsp24), and 1_16 during the metamorphosis of A. grahami pupae. First, we examined seven candidate reference genes (ribosomal protein 49 (RP49), 18S ribosomal RNA (18S rRNA), 28S ribosomal RNA (28S rRNA), beta-tubulin at 56D (β-tubulin), Ribosomal protein L23 (RPL23), glutathione S-transferase (GST1), and Actin. Three widely used algorithms (NormFinder, BestKeeper, and geNorm) were applied to evaluate the mRNA levels of reference gene candidates in puparium at three stable temperatures (15, 22, and 27°C). Next, mRNA expression of Hsp23, Hsp24, and 1_16 during A. grahami metamorphosis was examined. We demonstrated that mRNA expression levels of Hsp23, Hsp24, and 1_16 showed time-specific regulation. In summary, our study identified three gene markers for the intrapuparial period of A. grahami and might provide a potential application in PMImin estimation.

Molecular tools-advances, opportunities and prospects for the control of parasites of veterinary importance

The recent advancement in genome sequencing facilities, proteomics, transcriptomics, and metabolomics of eukaryotes have opened door for employment of molecular diagnostic techniques for early detection of parasites and determining target molecules for formulating control strategies. It further leads to the introduction of several purified vaccines in the field of veterinary parasitology. Earlier, the conventional diagnostic methods was entirely based upon morphological taxonomy for diagnosis of parasites but nowadays improved molecular techniques help in phylogenetic study and open an another area of molecular taxonomy of parasites with high precision. Control measures based upon targeting endosymbionts in parasites like Dirofilaria immitis is also under exploration in veterinary parasitology. Metagenomics have added an inside story of parasites bionomics which have created havoc in human and animals population since centuries. Omics era is playing a key role in opening the new approaches on parasite biology. Various newer generations of safer vaccines like edible vaccines and subunit vaccines and diagnostic techniques based upon purified immunologically active epitopes have become commercially available against the parasites (helminths, protozoa and arthropod borne diseases). Nowadays, a transgenic and gene knock out studies using RNA interference and CRISPR are also helping in understanding the functions of genes and screening of target genes, which are not available before the advent of molecular tools. Molecular techniques had paramount impact on increasing the sensitivity of diagnostic tools, epidemiological studies and more importantly in controlling these diseases. This review is about the advancements in veterinary parasitology and their impact on the control of these pathogens.

Major SCP/TAPS protein expansion in Lucilia cuprina is associated with novel tandem array organisation and domain architecture

Background

Larvae of the Australian sheep blowfly, Lucilia cuprina, parasitise sheep by feeding on skin excretions, dermal tissue and blood, causing severe damage known as flystrike or myiasis. Recent advances in -omic technologies and bioinformatic data analyses have led to a greater understanding of blowfly biology and should allow the identification of protein families involved in host-parasite interactions and disease. Current literature suggests that proteins of the SCP (Sperm-Coating Protein)/TAPS (Tpx-1/Ag5/PR-1/Sc7) (SCP/TAPS) superfamily play key roles in immune modulation, cross-talk between parasite and host as well as developmental and reproductive processes in parasites.

Methods

Here, we employed a bioinformatics workflow to curate the SCP/TAPS protein gene family in L. cuprina. Protein sequence, the presence and number of conserved CAP-domains and phylogeny were used to group identified SCP/TAPS proteins; these were compared to those found in Drosophila melanogaster to make functional predictions. In addition, transcription levels of SCP/TAPS protein-encoding genes were explored in different developmental stages.

Results

A total of 27 genes were identified as belonging to the SCP/TAPS gene family: encoding 26 single-domain proteins each with a single CAP domain and a solitary double-domain protein containing two conserved cysteine-rich secretory protein/antigen 5/pathogenesis related-1 (CAP) domains. Surprisingly, 16 SCP/TAPS predicted proteins formed an extended tandem array spanning a 53 kb region of one genomic region, which was confirmed by MinION long-read sequencing. RNA-seq data indicated that these 16 genes are highly transcribed in all developmental stages (excluding the embryo).

Conclusions

Future work should assess the potential of selected SCP/TAPS proteins as novel targets for the control of L. cuprina and related parasitic flies of major socioeconomic importance.

Combined use of feature engineering and machine-learning to predict essential genes in Drosophila melanogaster

Characterizing genes that are critical for the survival of an organism (i.e. essential) is important to gain a deep understanding of the fundamental cellular and molecular mechanisms that sustain life. Functional genomic investigations of the vinegar fly, Drosophila melanogaster, have unravelled the functions of numerous genes of this model species, but results from phenomic experiments can sometimes be ambiguous. Moreover, the features underlying gene essentiality are poorly understood, posing challenges for computational prediction. Here, we harnessed comprehensive genomic-phenomic datasets publicly available for D. melanogaster and a machine-learning-based workflow to predict essential genes of this fly. We discovered strong predictors of such genes, paving the way for computational predictions of essentiality in less-studied arthropod pests and vectors of infectious diseases.

Towards reconstructing the dipteran demise of an ancient essential gene: E3 ubiquitin ligase Murine double minute

Genome studies have uncovered many examples of essential gene loss, raising the question of how ancient genes transition from essentiality to dispensability. We explored this process for the deeply conserved E3 ubiquitin ligase Murine double minute (Mdm), which is lacking in Drosophila despite the conservation of its main regulatory target, the cellular stress response gene p53. Conducting gene expression and knockdown experiments in the red flour beetle Tribolium castaneum, we found evidence that Mdm has remained essential in insects where it is present. Using bioinformatics approaches, we confirm the absence of the Mdm gene family in Drosophila, mapping its loss to the stem lineage of schizophoran Diptera and Pipunculidae (big-headed flies), about 95-85 million years ago. Intriguingly, this gene loss event was preceded by the de novo origin of the gene Companion of reaper (Corp), a novel p53 regulatory factor that is characterized by functional similarities to vertebrate Mdm2 despite lacking E3 ubiquitin ligase protein domains. Speaking against a 1:1 compensatory gene gain/loss scenario, however, we found that hoverflies (Syrphidae) and pointed-wing flies (Lonchopteridae) possess both Mdm and Corp. This implies that the two p53 regulators have been coexisting for ~ 150 million years in select dipteran clades and for at least 50 million years in the lineage to Schizophora and Pipunculidae. Given these extensive time spans of Mdm/Corp coexistence, we speculate that the loss of Mdm in the lineage to Drosophila involved further acquisitions of compensatory gene activities besides the emergence of Corp. Combined with the previously noted reduction of an ancestral P53 contact domain in the Mdm homologs of crustaceans and insects, we conclude that the loss of the ancient Mdm gene family in flies was the outcome of incremental functional regression over long macroevolutionary time scales.

Gene content evolution in the arthropods

BACKGROUND

Arthropods comprise the largest and most diverse phylum on Earth and play vital roles in nearly every ecosystem. Their diversity stems in part from variations on a conserved body plan, resulting from and recorded in adaptive changes in the genome. Dissection of the genomic record of sequence change enables broad questions regarding genome evolution to be addressed, even across hyper-diverse taxa within arthropods.

RESULTS

Using 76 whole genome sequences representing 21 orders spanning more than 500 million years of arthropod evolution, we document changes in gene and protein domain content and provide temporal and phylogenetic context for interpreting these innovations. We identify many novel gene families that arose early in the evolution of arthropods and during the diversification of insects into modern orders. We reveal unexpected variation in patterns of DNA methylation across arthropods and examples of gene family and protein domain evolution coincident with the appearance of notable phenotypic and physiological adaptations such as flight, metamorphosis, sociality, and chemoperception.

CONCLUSIONS

These analyses demonstrate how large-scale comparative genomics can provide broad new insights into the genotype to phenotype map and generate testable hypotheses about the evolution of animal diversity.

Specific Gene Disruption in the Major Livestock Pests Cochliomyia hominivorax and Lucilia cuprina Using CRISPR/Cas9

Cochliomyia hominivorax and Lucilia cuprina are major pests of livestock. Their larvae infest warm-blooded vertebrates and feed on host's tissues, resulting in severe industry losses. As they are serious pests, considerable effort has been made to develop genomic resources and functional tools aiming to improve their management and control. Here, we report a significant addition to the pool of genome manipulation tools through the establishment of efficient CRISPR/Cas9 protocols for the generation of directed and inheritable modifications in the genome of these flies. Site-directed mutations were introduced in the C hominivorax and L cuprina yellow genes (ChY and LcY) producing lightly pigmented adults. High rates of somatic mosaicism were induced when embryos were injected with Cas9 ribonucleoprotein complexes (RNPs) pre-assembled with guide RNAs (sgRNAs) at high concentrations. Adult flies carrying disrupted yellow alleles lacked normal pigmentation (brown body phenotype) and efficiently transmitted the mutated alleles to the subsequent generation, allowing the rapid creation of homozygous strains for reverse genetics of candidate loci. We next used our established CRISPR protocol to disrupt the C hominivorax transformer gene (Chtra). Surviving females carrying mutations in the Chtra locus developed mosaic phenotypes of transformed ovipositors with characteristics of male genitalia while exhibiting abnormal reproductive tissues. The CRISPR protocol described here is a significant improvement on the existing toolkit of molecular methods in calliphorids. Our results also suggest that Cas9-based systems targeting Chtra and Lctra could be an effective means for controlling natural populations of these important pests.

Genome and Ontogenetic-Based Transcriptomic Analyses of the Flesh Fly, Sarcophaga bullata

The flesh fly, Sarcophaga bullata, is a widely-used model for examining the physiology of insect diapause, development, stress tolerance, neurobiology, and host-parasitoid interactions. Flies in this taxon are implicated in myiasis (larval infection of vertebrates) and feed on carrion, aspects that are important in forensic studies. Here we present the genome of S. bullata, along with developmental- and reproduction-based RNA-Seq analyses. We predict 15,768 protein coding genes, identify orthology in relation to closely related flies, and establish sex and developmental-specific gene sets based on our RNA-Seq analyses. Genomic sequences, predicted genes, and sequencing data sets have been deposited at the National Center for Biotechnology Information. Our results provide groundwork for genomic studies that will expand the flesh fly’s utility as a model system.

Insecticide activity of Curcuma longa (leaves) essential oil and its major compound α-phellandrene against Lucilia cuprina larvae (Diptera: Calliphoridae): Histological and ultrastructural biomarkers assessment

Lucilia cuprina, known as the Australian blowfly, is of high medico-sanitary and veterinary importance due to its ability to induce myiasis. Synthetic products are the most frequent form of fly control, but their indiscriminate use has selected for resistant populations and accounted for high levels of residues in animal products. This study aimed to assess the effect of essential oil from leaves of Curcuma longa (CLLEO), and its major compound α-phellandrene against L. cuprina L3. An additional goal was to determine the morphological alterations in target organs/tissues through ultrastructural assessment (SEM) and light microscopy, as well as macroscopic damage to cuticle induced by CLLEO. Groups of 20 L3 were placed on filter paper impregnated with increasing concentrations of CLLEO (0.15 to 2.86 μL/cm²) and α-phellandrene (0.29 to 1.47 μL/cm²). Efficacy was determined by quantifying L3 mortality 6, 24 and 48 h after contact with CLLEO and by measuring the structural damage to L3. CLLEO and α-phellandrene inhibited adult emergence by 96.22 and 100%, respectively. Macroscopic cuticle damage, appeared as diffuse pigment and darkening of larval body, was caused by both extracts. The SEM revealed dryness on the cuticle surface, distortion of the sensorial structures and general degeneration in treated L3. Furthermore, alterations in target organs (digestive tract, fat body and brain) were noticed and shall be used as biomarkers in future attempts to elucidate the mechanism of action of these compounds. The vacuolar degeneration and pyknotic profiles observed in the brain tissue of treated larvae with both extracts and the decreased motility within <6 h after treatment leads us to suggest a neurotoxic activity of the products. This work demonstrates the potential use of CLLEO and α-phellandrene as bioinsecticides to be used against L. cuprina, representing an ecofriendly alternative for myiasis control in humans and animals.

The Battle Against Flystrike - Past Research and New Prospects Through Genomics

Flystrike, or cutaneous myiasis, is caused by blow fly larvae of the genus Lucilia. This disease is a major problem in countries with large sheep populations. In Australia, Lucilia cuprina (Wiedemann, 1830) is the principal fly involved in flystrike. While much research has been conducted on L. cuprina, including physical, chemical, immunological, genetic and biological investigations, the molecular biology of this fly is still poorly understood. The recent sequencing, assembly and annotation of the draft genome and analyses of selected transcriptomes of L. cuprina have given a first global glimpse of its molecular biology and insights into host-fly interactions, insecticide resistance genes and intervention targets. The present article introduces L. cuprina, flystrike and associated issues, details past control efforts and research foci, reviews salient aspects of the L. cuprina genome project and discusses how the new genomic and transcriptomic resources for this fly might accelerate fundamental molecular research of L. cuprina towards developing new methods for the treatment and control of flystrike.