Blowflies of the genus Lucilia (Insecta: Diptera: Calliphoridae) are important in forensic science, medical research, and agriculture, particularly due to their role as facultative ectoparasites of domesticated sheep. The Australian sheep blowfly, L. cuprina dorsalis, is a primary causative agent of cutaneous myiasis, leading to substantial economic losses to the wool industry. In contrast, the closely related subspecies L. cuprina cuprina does not induce myiasis in sheep, but functions as a necrophagous subspecies in urban regions of Australia. Differences in the physiology, behavior, and host interactions of these subspecies might be attributed to their bacterial community structures as demonstrated in other blowfly species. However, our understanding of the bacterial communities associated with these flies is limited. This study aimed to characterize the microbiomes of two geographically distinct subspecies of L. cuprina—L. c. dorsalis (rural) and L. c. cuprina (urban) using 16S rRNA gene sequencing. The findings indicate that, while both subspecies harbor broadly similar microbial compositions, their respective urban and rural environments influence the prevalence of specific bacterial genera. The bacterial communities of these flies were predominantly composed of members of the phyla Proteobacteria, Firmicutes, Bacteroidetes, and Actinobacteria. The genus Ignatzschineria was highly abundant in L. c. dorsalis, whereas Pseudomonas was more prevalent in L. c. cuprina, highlighting distinct microbial compositions between the two subspecies. These findings provide fundamental insights into the microbial composition of Lucilia species, contributing to a broader understanding of their ecological and biological roles, with potential applications in forensic entomology, medical research, and agricultural pest management. • The microbial communities of Lucilia cuprina dorsalis and Lucilia cuprina cuprina have a broadly similar composition, but the abundance of particular bacterial genera is notably different. • Proteobacteria is the most abundant phylum in both subspecies. • The genus Ignatzschineria is notably abundant in Lucilia cuprina dorsalis, while Pseudomonas is predominant in Lucilia cuprina cuprina, emphasizing differences between the microbial community members between these subspecies.

Bactrocera zonata is a highly invasive agricultural pest that causes extensive damage to fruit crops. The Sterile Insect Technique (SIT), a species-specific and environmentally friendly pest control method, significantly benefits from the availability of Genetic Sexing Strains (GSSs) that enable efficient mass production of males for sterile release. However, no GSS currently exists for B. zonata limiting SIT applications targeting this important invasive pest. Here, we report two key advancements toward GSS development in this species. First, we present a high-quality, chromosome-level genome assembly from male B. zonata, identifying two scaffolds derived from the Y chromosome, which represent potential targets for future male-specific genetic engineering. Second, we demonstrate the feasibility of CRISPR/Cas9 genome editing in B. zonata by generating stable, homozygous white-eye mutants through targeted disruption of the conserved white-eye gene. This visible, recessive phenotype serves as a proof-of-concept for developing selectable markers in this species. Together, these results provide foundational genomic and genetic tools to support the development of GSSs in B. zonata, advancing the potential for sustainable, genetics-based pest control strategies.

Flystrike (cutaneous myiasis) is caused by blowfly larvae of the genus Lucilia. This disease is a major obstacle to sustainable global sheep and wool production. Flystrike control relies primarily on breech modification surgery (mulesing) and insecticidal treatment; however, control is constantly compromised by the emergence and spread of insecticide-resistance. Preventing the spread of resistance is severely hindered by a limited understanding of genetic variation, structure and gene flow within and among Lucilia populations. Australia is one of the world’s largest producers of sheep and wool products, where Lucilia cuprina dorsalis is the major cause of flystrike. Here, we collected 2,034 Lucilia cuprina dorsalis individuals among 86 populations from sheep-grazing regions across the continent. Each fly was genetically characterised at 20,000 loci using DArTseq, a reduced complexity genome sequencing strategy. Three genetically distinct population clusters (i.e., Western Australia, Eastern Australia and Tasmania) were revealed through population structure analyses. This investigation into population structure and gene flow yields significant insights into the genetic composition of diverse L. c. dorsalis populations throughout Australia. These findings will be essential for the sustainable management of flystrike on a global scale and for addressing the ongoing challenge of insecticide resistance.

Bactrocera tryoni, the Queensland fruit fly, is among the most damaging insect pests to the Australian horticultural industry as larvae infest ripening fruits or vegetables prior to harvest. Genetic biocontrol using Sterile Insect Technique (SIT) programs have been used to successfully suppress populations, via mass release of factory-reared sterile males that mate with wild females. Bi-sex flies are currently used for releases, although the efficiency of these control programs could be improved through using genetic sexing strains that eliminate females early during development, as they are not required for SIT. Here we used CRISPR/Cas9 mutagenesis to modify two nucleotides in the B. tryoni gene shibire, which created a proline to serine amino acid substitution and produced a temperature sensitive phenotype. Shibire is an essential GTPase required in endocytosis and synaptic vesicle recycling, and classical mutagenic screens in the vinegar fly Drosophila melanogaster previously identified temperature sensitive alleles including shits4 that results in adult paralysis. In B. tryoni, the shits4 mutant strain exhibited similar adult paralytic phenotypes when exposed to high temperatures, as well as temperature dependent lethality at egg, larval and pupal stages when subjected to heat treatment above standard rearing temperatures. These temperature sensitive phenotypes could be adapted to develop a SIT genetic sexing strain for conditional elimination of females prior to sterile releases, to improve efficiency and reduce costs.

The remarkable diversity of insect pigmentation offers a captivating avenue for studying evolution and genetics. In tephritids, understanding the molecular basis of mutant traits is also crucial for applied entomology, enabling the creation of genetic sexing strains through genome editing, thus facilitating sex-sorting before sterile insect releases. Here, we present evidence from classical and modern genetics showing that the black pupae (bp) phenotype in the GUA10 strain of Anastrepha ludens is caused by a large deletion at the ebony locus, removing the gene’s entire coding region. Targeted knockout of ebony induced analogous bp phenotypes across six major tephritid agricultural pests, demonstrating that disruption of Ebony alone is sufficient to produce the mutant trait in distantly related species. This functional characterization further allowed a deeper exploration of Ebony’s role in pigmentation and development across life stages in diverse species. Our findings offer key insights for molecular engineering of sexing strains based on the bp marker and for future evolutionary developmental biology studies in tephritids. This study integrates classical and modern genetics to show that the ebony gene is linked to the black pupae phenotype in tephritids, offering new insights into Ebony’s role in the evolutionary developmental biology of this diverse group of flies

The sterile insect technique can suppress and eliminate population outbreaks of the Australian horticultural pest, Bactrocera tryoni, the Queensland fruit fly. Sterile males mate with wild females that produce inviable embryos, causing population suppression or elimination. Current sterile insect releases are mixed sex, as the efficient removal of unrequired factory-reared females is not yet possible. In this paper, we assessed the known Drosophila melanogaster temperature-sensitive embryonic lethal alleles shibire (G268D, shits1) and RNA polymerase II 215 (R977C, RpII215ts) for potential use in developing B. tryoni genetic sexing strains (GSS) for the conditional removal of females. Complementation tests in D. melanogaster wild-type or temperature-sensitive genetic backgrounds were performed using the GAL4–UAS transgene expression system. A B. tryoni wild-type shibire isoform partially rescued Drosophila temperature lethality at 29°C by improving survivorship to pupation, while expressing B. tryoni shits1 failed to rescue the lethality, supporting a temperature-sensitive phenotype. Expression of the B. tryoni RpII215 wild-type protein rescued the lethality of D. melanogaster RpII215ts flies at 29°C. Overexpressing the B. tryoni RpII215ts allele in the D. melanogaster wild-type background unexpectedly produced a dominant lethal phenotype at 29°C. The B. tryoni shibire and RpII215 wild-type alleles were able to compensate, to varying degrees, for the function of the D. melanogaster temperature-sensitive proteins, supporting functional conservation across species. Shibire and RpII215 hold potential for developing insect strains that can selectively kill using elevated temperatures; however, alleles with milder effects than shits1 will need to be considered.

A synthetic gene drive that targets haplolethal genes on the X chromosome can skew the sex ratio toward males. Like an "X-shredder," it does not involve "homing," and that has advantages including the reduction of gene drive resistance allele formation. We examine this "X-poisoning" strategy by targeting 4 of the 11 known X-linked haplolethal/haplosterile genes of Drosophila melanogaster with CRISPR/Cas9. We find that targeting the wupA gene during spermatogenesis skews the sex ratio so fewer than 14% of progeny are daughters. That is unless we cross the mutagenic males to X^XY female flies that bear attached-X chromosomes, which reverses the inheritance of the poisoned X chromosome so that sons inherit it from their father, in which case only 2% of the progeny are sons. These sex ratio biases suggest that most of the CRISPR/Cas9 mutants we induced in the wupA gene are haplolethal but some are recessive lethal. The males generating wupA mutants do not suffer from reduced fertility; rather, the haplolethal mutants arrest development in the late stages of embryogenesis well after fertilized eggs have been laid. This provides a distinct advantage over genetic manipulation strategies involving sterility which can be countered by the remating of females. We also find that wupA mutants that destroy the nuclear localization signal of shorter isoforms are not haplolethal as long as the open reading frame remains intact. Like D. melanogaster, wupA orthologs of Drosophila suzukii and Anopheles mosquitos are found on X chromosomes making wupA a viable X-poisoning target in multiple species.

The remarkable diversity of insect pigmentation offers a captivating avenue for exploring evolution and genetics. In tephritid fruit flies, decoding the molecular pathways underlying pigmentation traits also plays a central role in applied entomology. Mutant phenotypes like the black pupae (bp) have long been used as a component of genetic sexing strains, allowing male-only release in tephritid sterile insect technique applications. However, the genetic basis of bp remains largely unknown. Here, we present independent evidence from classical and modern genetics showing that the bp phenotype in the GUA10 strain of the Mexican fruit fly, Anastrepha ludens , is caused by a large deletion at the ebony locus resulting in the removal of the entire protein-coding region of the gene. Targeted knockout of ebony induced analogous bp phenotypes across six tephritid species spanning over 50 million years of divergent evolution. This functionally validated our findings and allowed for a deeper investigation into the role of Ebony in pigmentation and development in these species. Our study offers fundamental knowledge for developing new sexing strains based on the bp marker and for future evolutionary developmental biology studies in tephritid fruit flies. ### Competing Interest Statement The authors have declared no competing interest.

BACKGROUND One of the proposed applications of gene drives has been to revert pesticide resistant mutations back to the ancestral susceptible state. Insecticides that have become ineffective because of the rise of resistance could have reinvigorated utility and be used to suppress pest populations again, perhaps at lower application doses. RESULTS We have created a laboratory model for susceptibility gene drives that replaces field-selected resistant variants of the acetylcholine esterase ( Ace ) locus of Drosophila melanogaster with ancestral susceptible variants. We constructed a CRISPR/Cas9 homing drive and found that homing occurred in many genetic backgrounds with varying efficiencies. While the drive itself could not be homozygosed, it converted resistant alleles into susceptible ones and produced recessive lethal alleles that could suppress populations. Our studies provided evidence for two distinct classes of Gene Drive Resistance (GDR): rather than being mediated by the conventional Non-Homologous End-joining (NHEJ) pathway, one seemed to involve short homologous repair and the other was defined by genetic background. Additionally, we used simulations to explore a distinct application of susceptibility drives; the use of chemicals to prevent the spread of synthetic gene drives into protected areas. CONCLUSIONS Insecticide susceptibility gene drives could be useful tools to control pest insects however problems associated with particularities of the target loci and GDR will need to be overcome for them to be effective. Furthermore, realistic patterns of pest dispersal and high insecticide exposure rates would be required if susceptibility were to be useful as a ‘safety-switch’ to prevent the unwanted spread of gene drives.

Bacillus thuringiensis (Bt) three-domain Cry toxins are highly successful biological pesticides; however, the mechanism through which they cause death to targeted larval midgut cells is not fully understood. Herein, we challenged transgenic Bt-susceptible Drosophila melanogaster larvae with moderate doses of activated Cry1Ac toxin and assessed the midgut tissues after one, three, and five hours using transmission electron microscopy and transcriptome sequencing. Larvae treated with Cry1Ac showed dramatic changes to their midgut morphology, including shortened microvilli, enlarged vacuoles, thickened peritrophic membranes, and swelling of the basal labyrinth, suggesting water influx. Transcriptome analysis showed that innate immune responses were repressed, genes involved with cell death pathways were largely unchanged, and mitochondria-related genes were strongly upregulated following toxin exposure. Defective mitochondria produced after toxin exposure were likely to contribute to significant levels of oxidative stress, which represent a common physiological response to a range of toxic chemicals. Significant reductions in both mitochondrial aconitase activity and ATP levels in the midgut tissue supported a rapid increase in reactive oxygen species (ROS) following exposure to Cry1Ac. Overall, these findings support the role of water influx, midgut cell swelling, and ROS activity in response to moderate concentrations of Cry1Ac.

Mass releases of sterilized male insects, in the frame of sterile insect technique programs, have helped suppress insect pest populations since the 1950s. In the major horticultural pests Bactrocera dorsalis, Ceratitis capitata , and Zeugodacus cucurbitae , a key phenotype white pupae (wp) has been used for decades to selectively remove females before releases, yet the gene responsible remained unknown. Here, we use classical and modern genetic approaches to identify and functionally characterize causal wp − mutations in these distantly related fruit fly species. We find that the wp phenotype is produced by parallel mutations in a single, conserved gene. CRISPR/Cas9-mediated knockout of the wp gene leads to the rapid generation of white pupae strains in C. capitata and B. tryoni . The conserved phenotype and independent nature of wp − mutations suggest this technique can provide a generic approach to produce sexing strains in other major medical and agricultural insect pests.

BACKGROUND Although decoding the molecular mechanisms underlying insecticide resistance has often proven difficult, recent progress has revealed that specific mutations in the ryanodine receptor (RyR) of the diamondback moth, Plutella xylostella, can confer resistance to diamide insecticides. The extent to which specific RyR mutations contribute to the diamide resistance phenotype, the associated genetic traits and fitness costs remain limited. RESULTS Three field-evolved PxRyR mutations (G4946E, I4790 M, and I4790 K) were respectively introgressed into a common susceptible background strain (IPP-S) of P. xylostella with marker-assisted backcrossing. The mutations alone can result in moderate to high levels of resistance to five commercial diamides (flubendiamide, chlorantraniliprole, cyantraniliprole, tetraniliprole, and cyclaniliprole), and the resistance intensity mediated by the three mutations was hierarchical in order of I4790 K (1199- to >2778-fold) > G4946E (39- to 739-fold) > I4790 M (16- to 57-fold). Flubendiamide resistance was autosomal and incompletely recessive, and was significantly linked with the introgressed mutations in the three constructed strains. In addition, the resistance levels to flubendiamide of hybrid progeny from any two resistant strains fell in between the status of their parents. Furthermore, by comparing the net replacement rate, the fitness of 4946E, 4790 M and 4790 K strains were 0.77, 0.93 and 0.92 relative to the IPP-S strain, respectively. CONCLUSION Three independent PxRyR mutations confer varying degrees of resistance to diamides in P. xylostella. Among the three mutations, I4790 K confers highest levels of resistance (> 1000-fold) to all five commercial diamides. The findings can guide resistance management practices for diamides in P. xylostella and other arthropods.

The diamondback moth, Plutella xylostella (L.), is a highly mobile brassica crop pest with worldwide distribution and can rapidly evolve resistance to insecticides, including group 28 diamides. Reference genomes assembled using Illumina sequencing technology have provided valuable resources to advance our knowledge regarding the biology, origin and movement of diamondback moth, and more recently with its sister species, Plutella australiana. Here we apply a trio binning approach to sequence and annotate a chromosome level reference genome of P. xylostella using PacBio Sequel and Dovetail Hi-C sequencing technology and identify a point mutation that causes resistance to commercial diamides. A P. xylostella population collected from brassica crops in the Lockyer Valley, Australia (LV-R), was reselected for chlorantraniliprole resistance then a single male was crossed to a P. australiana female and a hybrid pupa sequenced. A chromosome level 328 Mb P. xylostella genome was assembled with 98.1% assigned to 30 autosomes and the Z chromosome. The genome was highly complete with 98.4% of BUSCO Insecta genes identified and RNAseq informed protein prediction annotated 19,002 coding genes. The LV-R strain survived recommended field application doses of chlorantraniliprole, flubendiamide and cyclaniliprole. Some hybrids also survived these doses, indicating significant departure from recessivity, which has not been previously documented for diamides. Diamide chemicals modulate insect Ryanodine Receptors (RyR), disrupting calcium homeostasis, and we identified an amino acid substitution (I4790K) recently reported to cause diamide resistance in a strain from Japan. This chromosome level assembly provides a new resource for insect comparative genomics and highlights the emergence of diamide resistance in Australia. Resistance management plans need to account for the fact that resistance is not completely recessive.

A major obstacle of sterile insect technique (SIT) programs is the availability of robust sex-separation systems for conditional removal of females. Sterilized male-only releases improve SIT efficiency and cost-effectiveness for agricultural pests, whereas it is critical to remove female disease-vector pests prior to release as they maintain the capacity to transmit disease. Some of the most successful Genetic Sexing Strains (GSS) reared and released for SIT control were developed for Mediterranean fruit fly (Medfly), Ceratitis capitata, and carry a temperature sensitive lethal (tsl) mutation that eliminates female but not male embryos when heat treated. The Medfly tsl mutation was generated by random mutagenesis and the genetic mechanism causing this valuable heat sensitive phenotype remains unknown. Conditional temperature sensitive lethal mutations have also been developed using random mutagenesis in the insect model, Drosophila melanogaster, and were used for some of the founding genetic research published in the fields of neuro- and developmental biology. Here we review mutations in select D. melanogaster genes shibire, Notch, RNA polymerase II 215kDa, pale, transformer-2, Dsor1 and CK2α that cause temperature sensitive phenotypes. Precise introduction of orthologous point mutations in pest insect species with CRISPR/Cas9 genome editing technology holds potential to establish GSSs with embryonic lethality to improve and advance SIT pest control.

The analysis and interpretation of datasets generated through sequencing large numbers of individual genomes is becoming commonplace in population and evolutionary genetic studies. Here we introduce geaR, a modular R package for evolutionary analysis of genome-wide genotype data. The package leverages the Genomic Data Structure (GDS) format, which enables memory and time efficient querying of genotype datasets compared to standard VCF genotype files. geaR utilizes GRange object classes to partition an analysis based on features from GFF annotation files, select codons based on position or degeneracy, and construct both positional and coordinate genomic windows. Tests of genetic diversity (eg. , ) and admixture along with tree building and sequence output, can be carried out on partitions using a single function regardless of sample ploidy or number of observed alleles. The package and associated documentation are available on GitHub at .

The diamondback moth, Plutella xylostella, is a cosmopolitan pest and the first species to develop field resistance to toxins from the gram-positive bacterium Bacillus thuringiensis (Bt). Although previous work has suggested that mutations of ATP-binding cassette transporter subfamily C2 (ABCC2) or C3 (ABCC3) genes can confer Cry1Ac resistance, here we reveal that P. xylostella requires combined mutations in both PxABCC2 and PxABCC3 to achieve high-level Cry1Ac resistance, rather than simply a mutation of either gene. We identified natural mutations of PxABCC2 and PxABCC3 that concurrently occurred in a Cry1Ac-resistant strain (Cry1S1000) of P. xylostella, with a mutation (RA2) causing the mis-splicing of PxABCC2 and another mutation (RA3) leading to the premature termination of PxABCC3. Genetic linkage analysis showed that RA2 and RA3 were tightly linked to Cry1Ac resistance. Introgression of RA2 and RA3 enabled a susceptible strain (G88) of P. xylostella to obtain high resistance to Cry1Ac, confirming that these genes confer resistance. To further support the role of PxABCC2 and PxABCC3 in Cry1Ac resistance, frameshift mutations were introduced into PxABCC2 and PxABCC3 singly and in combination in the G88 strain with CRISPR/Cas9 mediated mutagenesis. Bioassays of CRISPR-based mutant strains, plus genetic complementation tests, demonstrated that the deletion of PxABCC2 or PxABCC3 alone provided < 4-fold tolerance to Cry1Ac, while disruption of both genes together conferred >8,000-fold resistance to Cry1Ac, suggesting the redundant/complementary roles of PxABCC2 and PxABCC3. This work advances our understanding of Bt resistance in P. xylostella by demonstrating mutations within both PxABCC2 and PxABCC3 genes are required for high-level Cry1Ac resistance.

Meadow brown butterflies (Maniola jurtina) on the Isles of Scilly represent an ideal model in which to dissect the links between genotype, phenotype and long-term patterns of selection in the wild - a largely unfulfilled but fundamental aim of modern biology. To meet this aim, a clear description of genotype is required. Here we present the draft genome sequence of M. jurtina to serve as a founding genetic resource for this species. Seven libraries were constructed using pooled DNA from five wild caught spotted females and sequenced using Illumina, PacBio RSII and MinION technology. A novel hybrid assembly approach was employed to generate a final assembly with an N50 of 214 kb (longest scaffold 2.9 Mb). The sequence assembly described here predicts a gene count of 36,294 and includes variants and gene duplicates from five genotypes. Core BUSCO (Benchmarking Universal Single-Copy Orthologs) gene sets of Arthropoda and Insecta recovered 90.5% and 88.7% complete and single-copy genes respectively. Comparisons with 17 other Lepidopteran species placed 86.5% of the assembled genes in orthogroups. Our results provide the first high-quality draft genome and annotation of the butterfly M. jurtina.

The diamondback moth, Plutella xylostella is a cosmopolitan pest that has evolved resistance to all classes of insecticide, and costs the world economy an estimated US $4-5 billion annually. We analyse patterns of variation among 532 P. xylostella genomes, representing a worldwide sample of 114 populations. We find evidence that suggests South America is the geographical area of origin of this species, challenging earlier hypotheses of an Old-World origin. Our analysis indicates that Plutella xylostella has experienced three major expansions across the world, mainly facilitated by European colonization and global trade. We identify genomic signatures of selection in genes related to metabolic and signaling pathways that could be evidence of environmental adaptation. This evolutionary history of P. xylostella provides insights into transoceanic movements that have enabled it to become a worldwide pest.

Background Pest eradication using the Sterile Insect Technique (SIT) involves high-density releases of sterilized males that mate with wild females and ultimately suppress the population. Sterilized females are not required for SIT and their removal or separation from males prior to release remains challenging. In order to develop genetic sexing strains (GSS), conditional traits such as temperature sensitive lethality are required. Results Here we introduce a known Drosophila melanogaster temperature sensitive embryonic lethal mutation into Bactrocera tryoni , a serious horticultural pest in Australia. A non-synonymous point mutation in the D. melanogaster gene shibire causes embryonic lethality at 29 °C and we successfully used CRISPR/Cas9 technology to recreate the orthologous shibire temperature sensitive-1 ( shi ts1 ) mutation in B. tryoni . Genotypic analyses over three generations revealed that a high fitness cost was associated with the shi ts1 mutant allele and shi ts1 homozygotes were not viable at 21 °C, which is a more severe phenotype than that documented in D. melanogaster . Conclusions We have demonstrated the first successful use of CRISPR/Cas9 to introduce precise single base substitutions in an endogenous gene via homology-directed repair in an agricultural pest insect and this technology can be used to trial other conditional mutations for the ultimate aim of generating genetic sexing strains for SIT.

Irradiated Queensland fruit flies ( Bactrocera tryoni ) used in Sterile Insect Technique (SIT) programmes are marked with fluorescent dyes to distinguish them from wild flies when recaptured in monitoring traps. However, coating sterile pupae with powdered dyes can reduce emergence rates and fly quality and can sometimes produce insufficiently certain discrimination through inadequate coating or because the dye is transferred to wild flies through contact. Here we created a phenotypically distinct B. tryoni strain that lacks typical melanisation patterns through CRISPR/Cas9-mediated mutagenesis of tandemly duplicated yellow-y genes and then assessed effects of this visible trait on fly performance. Recessive mutations are only required in one of these copies to restrict melanisation and generate a phenotype clearly distinguished from wild type. The yellow strain showed significant declines in eclosion rates and in the percentage of fliers directly after emergence. Locomotor activity was greater in the yellow strain, and these mutations did not generally affect mating probability, copula latency, or copula duration. The longevity of yellow flies was approximately 10 days shorter than wild-type flies in both sexes. Overall, replacing dyes with yellow body marker for SIT can simplify production, eliminate a step that is known to reduce fly quality, remove potentially hazardous dyes from production, enable accurate discrimination from wild flies, and improve cost-effectiveness; however, direct comparisons of the decrements in performance associated with dyes on mass-reared wild-type flies and disruption of yellow-y genes are now required to determine the relative suitability of these marking methods for B. tryoni SIT.