A male-killing Wolbachia endosymbiont is concealed by another endosymbiont and a nuclear suppressor.

Bacteria that live inside the cells of insect hosts (endosymbionts) can alter the reproduction of their hosts, including the killing of male offspring (male killing, MK). MK has only been described in a few insects, but this may reflect challenges in detecting MK rather than its rarity. Here, we identify MK Wolbachia at a low frequency (around 4%) in natural populations of Drosophila pseudotakahashii. MK Wolbachia had a stable density and maternal transmission during laboratory culture, but the MK phenotype which manifested mainly at the larval stage was lost rapidly. MK Wolbachia occurred alongside a second Wolbachia strain expressing a different reproductive manipulation, cytoplasmic incompatibility (CI). A genomic analysis highlighted Wolbachia regions diverged between the 2 strains involving 17 genes, and homologs of the wmk and cif genes implicated in MK and CI were identified in the Wolbachia assembly. Doubly infected males induced CI with uninfected females but not females singly infected with CI-causing Wolbachia. A rapidly spreading dominant nuclear suppressor genetic element affecting MK was identified through backcrossing and subsequent analysis with ddRAD SNPs of the D. pseudotakahashii genome. These findings highlight the complexity of nuclear and microbial components affecting MK endosymbiont detection and dynamics in populations and the challenges of making connections between endosymbionts and the host phenotypes affected by them.

A decade of stability for wMel Wolbachia in natural Aedes aegypti populations.

Mosquitoes carrying Wolbachia endosymbionts are being released in many countries for arbovirus control. The wMel strain of Wolbachia blocks Aedes-borne virus transmission and can spread throughout mosquito populations by inducing cytoplasmic incompatibility. Aedes aegypti mosquitoes carrying wMel were first released into the field in Cairns, Australia, over a decade ago, and with wider releases have resulted in the near elimination of local dengue transmission. The long-term stability of Wolbachia effects is critical for ongoing disease suppression, requiring tracking of phenotypic and genomic changes in Wolbachia infections following releases. We used a combination of field surveys, phenotypic assessments, and Wolbachia genome sequencing to show that wMel has remained stable in its effects for up to a decade in Australian Ae. aegypti populations. Phenotypic comparisons of wMel-infected and uninfected mosquitoes from near-field and long-term laboratory populations suggest limited changes in the effects of wMel on mosquito fitness. Treating mosquitoes with antibiotics used to cure the wMel infection had limited effects on fitness in the next generation, supporting the use of tetracycline for generating uninfected mosquitoes without off-target effects. wMel has a temporally stable within-host density and continues to induce complete cytoplasmic incompatibility. A comparison of wMel genomes from pre-release (2010) and nine years post-release (2020) populations show few genomic differences and little divergence between release locations, consistent with the lack of phenotypic changes. These results indicate that releases of Wolbachia-infected mosquitoes for population replacement are likely to be effective for many years, but ongoing monitoring remains important to track potential evolutionary changes.

Close Kin Dyads Indicate Intergenerational Dispersal and Barriers.

AbstractThe movement of individuals through continuous space is typically constrained by dispersal ability and dispersal barriers. A range of approaches have been developed to investigate these. Kindisperse is a new approach that infers recent intergenerational dispersal (σ) from close kin dyads and appears particularly useful for investigating taxa that are difficult to observe individually. This study, focusing on the mosquito Aedes aegypti, shows how the same close kin data can also be used for barrier detection. We empirically demonstrate this new extension of the method using genome-wide sequence data from 266 Ae. aegypti. First, we use the spatial distribution of full-sib dyads collected within one generation to infer past movements of ovipositing female mosquitoes. These dyads indicated the relative barrier strengths of two roads and performed favorably against alternative genetic methods for detecting barriers. We then use Kindisperse to quantify recent intergenerational dispersal (σ=81.5-197.1 m generation-1/2) from the difference in variance between the sib and the first cousin spatial distributions and, from this, estimate effective population density (ρ=833-4,864 km-2). Dispersal estimates showed general agreement with those from mark-release-recapture studies. Barriers, σ, ρ, and neighborhood size (331-526) can inform forthcoming releases of dengue-suppressing Wolbachia bacteria into this mosquito population.

Heterogeneous patterns of heterozygosity loss in isolated populations of the threatened eastern barred bandicoot (Perameles gunnii).

Identifying and analysing isolated populations is critical for conservation. Isolation can make populations vulnerable to local extinction due to increased genetic drift and inbreeding, both of which should leave imprints of decreased genome-wide heterozygosity. While decreases in heterozygosity among populations are frequently investigated, fewer studies have analysed how heterozygosity varies among individuals, including whether heterozygosity varies geographically along lines of discrete population structure or with continuous patterns analogous to isolation by distance. Here we explore geographical patterns of differentiation and individual heterozygosity in the threatened eastern barred bandicoot (Perameles gunnii) in Tasmania, Australia, using genomic data from 85 samples collected between 2008 and 2011. Our analyses identified two isolated demes undergoing significant genetic drift, and several areas of fine-scale differentiation across Tasmania. We observed discrete genetic structures across geographical barriers and continuous patterns of isolation by distance, with little evidence of recent or historical migration. Using a recently developed analytical pipeline for estimating autosomal heterozygosity, we found individual heterozygosities varied within demes by up to a factor of two, and demes with low-heterozygosity individuals also still contained those with high heterozygosity. Spatial interpolation of heterozygosity scores clarified these patterns and identified the isolated Tasman Peninsula as a location where low-heterozygosity individuals were more common than elsewhere. Our results provide novel insights into the relationship between isolation-driven genetic structure and local heterozygosity patterns. These may help improve translocation efforts, by identifying populations in need of assistance, and by providing an individualised metric for identifying source animals for translocation.

Urban population structure and dispersal of an Australian mosquito (Aedes notoscriptus) involved in disease transmission.

Dispersal is a critical parameter for successful pest control measures as it determines the rate of movement across target control areas and influences the risk of human exposure. We used a fine-scale spatial population genomic approach to investigate the dispersal ecology and population structure of Aedes notoscriptus, an important disease transmitting mosquito at the Mornington Peninsula, Australia. We sampled and reared Ae. notoscriptus eggs at two time points from 170 traps up to 5 km apart and generated genomic data from 240 individuals. We also produced a draft genome assembly from a laboratory colony established from mosquitoes sampled near the study area. We found low genetic structure (Fst) and high coancestry throughout the study region. Using genetic data to identify close kin dyads, we found that mosquitoes had moved distances of >1 km within a generation, which is further than previously recorded. A spatial autocorrelation analysis of genetic distances indicated genetic similarity at >1 km separation, a tenfold higher distance than for a comparable population of Ae. aegypti, from Cairns, Australia. These findings point to high mobility of Ae. notoscriptus, highlighting challenges of localised intervention strategies. Further sampling within the same area 6 and 12 months after initial sampling showed that egg-counts were relatively consistent across time, and that spatial variation in egg-counts covaried with spatial variation in Wright's neighbourhood size (NS). As NS increases linearly with population density, egg-counts may be useful for estimating relative density in Ae. notoscriptus. The results highlight the importance of acquiring species-specific data when planning control measures.

Genetic stability of Aedes aegypti populations following invasion by wMel Wolbachia.

BACKGROUND: Wolbachia wMel is the most commonly used strain in rear and release strategies for Aedes aegypti mosquitoes that aim to inhibit the transmission of arboviruses such as dengue, Zika, Chikungunya and yellow fever. However, the long-term establishment of wMel in natural Ae. aegypti populations raises concerns that interactions between Wolbachia wMel and Ae. aegypti may lead to changes in the host genome, which could affect useful attributes of Wolbachia that allow it to invade and suppress disease transmission. RESULTS: We applied an evolve-and-resequence approach to study genome-wide genetic changes in Ae. aegypti from the Cairns region, Australia, where Wolbachia wMel was first introduced more than 10 years ago. Mosquito samples were collected at three different time points in Gordonvale, Australia, covering the phase before (2010) and after (2013 and 2018) Wolbachia releases. An additional three locations where Wolbachia replacement happened at different times across the last decade were also sampled in 2018. We found that the genomes of mosquito populations mostly remained stable after Wolbachia release, with population differences tending to reflect the geographic location of the populations rather than Wolbachia infection status. However, outlier analysis suggests that Wolbachia may have had an influence on some genes related to immune response, development, recognition and behavior. CONCLUSIONS: Ae. aegypti populations remained geographically distinct after Wolbachia wMel releases in North Australia despite their Wolbachia infection status. At some specific genomic loci, we found signs of selection associated with Wolbachia, suggesting potential evolutionary impacts can happen in the future and further monitoring is warranted.

Anthropogenic and natural barriers affect genetic connectivity in an Alpine butterfly.

Dispersal is a key biological process serving several functions including connectivity among populations. Habitat fragmentation caused by natural or anthropogenic structures may hamper dispersal, thereby disrupting genetic connectivity. Investigating factors affecting dispersal and gene flow is important in the current era of anthropogenic global change, as dispersal comprises a vital part of a species' resilience to environmental change. Using finescale landscape genomics, we investigated gene flow and genetic structure of the Sooty Copper butterfly (Lycaena tityrus) in the Alpine Ötz valley system in Austria. We found surprisingly high levels of gene flow in L. tityrus across the region. Nevertheless, ravines, forests, and roads had effects on genetic structure, while rivers did not. The latter is surprising as roads and rivers have a similar width and run largely in parallel in our study area, pointing towards a higher impact of anthropogenic compared with natural linear structures. Additionally, we detected eleven loci potentially under thermal selection, including ones related to membranes, metabolism, and immune function. This study demonstrates the usefulness of molecular approaches in obtaining estimates of dispersal and population processes in the wild. Our results suggest that, despite high gene flow in the Alpine valley system investigated, L. tityrus nevertheless seems to be vulnerable to anthropogenically-driven habitat fragmentation. With anthropogenic rather than natural linear structures affecting gene flow, this may have important consequences for the persistence of species such as the butterfly studied here in altered landscapes.

Spatial population genomics of a recent mosquito invasion.

Population genomic approaches can characterize dispersal across a single generation through to many generations in the past, bridging the gap between individual movement and intergenerational gene flow. These approaches are particularly useful when investigating dispersal in recently altered systems, where they provide a way of inferring long-distance dispersal between newly established populations and their interactions with existing populations. Human-mediated biological invasions represent such altered systems which can be investigated with appropriate study designs and analyses. Here we apply temporally restricted sampling and a range of population genomic approaches to investigate dispersal in a 2004 invasion of Aedes albopictus (the Asian tiger mosquito) in the Torres Strait Islands (TSI) of Australia. We sampled mosquitoes from 13 TSI villages simultaneously and genotyped 373 mosquitoes at genome-wide single nucleotide polymorphisms (SNPs): 331 from the TSI, 36 from Papua New Guinea (PNG) and four incursive mosquitoes detected in uninvaded regions. Within villages, spatial genetic structure varied substantially but overall displayed isolation by distance and a neighbourhood size of 232-577. Close kin dyads revealed recent movement between islands 31-203 km apart, and deep learning inferences showed incursive Ae. albopictus had travelled to uninvaded regions from both adjacent and nonadjacent islands. Private alleles and a co-ancestry matrix indicated direct gene flow from PNG into nearby islands. Outlier analyses also detected four linked alleles introgressed from PNG, with the alleles surrounding 12 resistance-associated cytochrome P450 genes. By treating dispersal as both an intergenerational process and a set of discrete events, we describe a highly interconnected invasive system.

Heterogeneous genetic invasions of three insecticide resistance mutations in Indo-Pacific populations of Aedes aegypti (L.).

Nations throughout the Indo-Pacific region use pyrethroid insecticides to control Aedes aegypti, the mosquito vector of dengue, often without knowledge of pyrethroid resistance status of the pest or origin of resistance. Two mutations (V1016G + F1534C) in the sodium channel gene (Vssc) of Ae. aegypti modify ion channel function and cause target-site resistance to pyrethroid insecticides, with a third mutation (S989P) having a potential additive effect. Of 27 possible genotypes involving these mutations, some allelic combinations are never seen whereas others predominate. Here, five allelic combinations common in Ae. aegypti from the Indo-Pacific region are described and their geographical distributions investigated using genome-wide SNP markers. We tested the hypothesis that resistance allele combinations evolved de novo in populations versus the alternative that dispersal of Ae. aegypti between populations facilitated genetic invasions of allele combinations. We used latent factor mixed-models to detect SNPs throughout the genome that showed structuring in line with resistance allele combinations and compared variation at SNPs within the Vssc gene with genome-wide variation. Mixed-models detected an array of SNPs linked to resistance allele combinations, all located within or in close proximity to the Vssc gene. Variation at SNPs within the Vssc gene was structured by resistance profile, whereas genome-wide SNPs were structured by population. These results demonstrate that alleles near to resistance mutations have been transferred between populations via linked selection. This indicates that genetic invasions have contributed to the widespread occurrence of Vssc allele combinations in Ae. aegypti in the Indo-Pacific region, pointing to undocumented mosquito invasions between countries.

Mechanical Resonances of Mobile Impurities in a One-Dimensional Quantum Fluid.

We study a one-dimensional interacting quantum liquid hosting a pair of mobile impurities causing backscattering. We determine the effective retarded interaction between the two impurities mediated by the liquid. We show that for strong backscattering this interaction gives rise to resonances and antiresonances in the finite-frequency mobility of the impurity pair. At the antiresonances, the two impurities remain at rest even when driven by a (small) external force. At the resonances, their synchronous motion follows the external drive in phase and reaches maximum amplitude. Using a perturbative renormalization group analysis in quantum tunneling across the impurities, we study the range of validity of our model. We predict that these mechanical antiresonances are observable in experiments on ultracold atom gases confined to one dimension.

Fine-scale landscape genomics helps explain the slow spatial spread of Wolbachia through the Aedes aegypti population in Cairns, Australia.

The endosymbiotic bacterium Wolbachia suppresses the capacity for arbovirus transmission in the mosquito Aedes aegypti, and can spread spatially through wild mosquito populations following local introductions. Recent introductions in Cairns, Australia have demonstrated slower than expected spatial spread. Potential reasons for this include: (i) barriers to Ae. aegypti dispersal; (ii) higher incidence of long-range dispersal; and (iii) intergenerational loss of Wolbachia. We investigated these three potential factors using genome-wide single-nucleotide polymorphisms (SNPs) and an assay for the Wolbachia infection wMel in 161 Ae. aegypti collected from Cairns in 2015. We detected a small but significant barrier effect of Cairns highways on Ae. aegypti dispersal using distance-based redundancy analysis and patch-based simulation analysis. We detected a pair of putative full-siblings in ovitraps 1312 m apart, indicating long-distance female movement likely mediated by human transport. We also found a pair of full-siblings of different infection status, indicating intergenerational loss of Wolbachia in the field. These three factors are all expected to contribute to the slow spread of Wolbachia through Ae. aegypti populations, though from our results it is unclear whether Wolbachia loss and long-distance movement are sufficiently common to reduce the speed of spatial spread appreciably. Our findings inform the strategic deployment of Wolbachia-infected mosquitoes during releases, and show how parameter estimates from laboratory studies may differ from those estimated using field data. Our landscape genomics approach can be extended to other host/symbiont systems that are being considered for biocontrol.

How Can Genomics Help or Hinder Wildlife Conservation?

Genomic data are becoming increasingly affordable and easy to collect, and new tools for their analysis are appearing rapidly. Conservation biologists are interested in using this information to assist in management and planning but are typically limited financially and by the lack of genomic resources available for non-model taxa. It is therefore important to be aware of the pitfalls as well as the benefits of applying genomic approaches. Here, we highlight recent methods aimed at standardizing population assessments of genetic variation, inbreeding, and forms of genetic load and methods that help identify past and ongoing patterns of genetic interchange between populations, including those subjected to recent disturbance. We emphasize challenges in applying some of these methods and the need for adequate bioinformatic support. We also consider the promises and challenges of applying genomic approaches to understand adaptive changes in natural populations to predict their future adaptive capacity.

Diversity and distribution of sodium channel mutations in Aedes albopictus (Diptera: Culicidae).

There is growing interest in insecticide resistance in the mosquito, Aedes albopictus (Skuse), as its potential for spreading diseases is increasing as urbanization and control efforts intensify. Here we review the presence and diversity of mutations in the voltage-sensitive sodium channel (Vssc) gene associated with pyrethroid resistance and report on additional surveys of these mutations in new populations with an analysis of their spread. The known diversity of these mutations has increased in recent years including the identification of 26 non-synonymous mutations, although phenotypic data associating mutations with resistance remain limited. We provide data on mutations in several new locations including those in Timor Leste, Indonesia, and Vanuatu. We use population genomic data from ddRAD analyses of target populations with the 1534C mutation to identify single nucleotide polymorphisms (SNPs) associated with the mutant to test for clustering of SNPs based on the presence of the 1534C mutation rather than population origin. Our findings suggest spread of resistance alleles via genetic invasion, which is further supported by patterns from a genome-wide principal components analysis. These data point to movement of resistance alleles across wide areas with likely impacts on local control options.

Global, asynchronous partial sweeps at multiple insecticide resistance genes in Aedes mosquitoes.

Aedes aegypti (yellow fever mosquito) and Ae. albopictus (Asian tiger mosquito) are globally invasive pests that confer the world's dengue burden. Insecticide-based management has led to the evolution of insecticide resistance in both species, though the genetic architecture and geographical spread of resistance remains incompletely understood. This study investigates partial selective sweeps at resistance genes on two chromosomes and characterises their spread across populations. Sweeps at the voltage-sensitive sodium channel (VSSC) gene on chromosome 3 correspond to one resistance-associated nucleotide substitution in Ae. albopictus and three in Ae. aegypti, including two substitutions at the same nucleotide position (F1534C) that have evolved and spread independently. In Ae. aegypti, we also identify partial sweeps at a second locus on chromosome 2. This locus contains 15 glutathione S-transferase (GST) epsilon class genes with significant copy number variation among populations and where three distinct genetic backgrounds have spread across the Indo-Pacific region, the Americas, and Australia. Local geographical patterns and linkage networks indicate VSSC and GST backgrounds probably spread at different times and interact locally with different genes to produce resistance phenotypes. These findings highlight the rapid global spread of resistance and are evidence for the critical importance of GST genes in resistance evolution.

Majorana qubit rotations in microwave cavities.

Majorana bound states have been proposed as building blocks for qubits on which certain operations can be performed in a topologically protected way using braiding. However, the set of these protected operations is not sufficient to realize universal quantum computing. We show that the electric field in a microwave cavity can induce Rabi oscillations between adjacent Majorana bound states. These oscillations can be used to implement an additional single-qubit gate. Supplemented with one braiding operation, this gate allows us to perform arbitrary single-qubit operations.

Evaluation of In2Care mosquito stations for suppression of the Australian backyard mosquito, Aedes notoscriptus (Diptera: Culicidae).

Aedes notoscriptus (Skuse) is a container-inhabiting mosquito endemic to Australia that vectors arboviruses and is suspected to transmit Mycobacterium ulcerans, the cause of Buruli ulcer. We evaluated the effectiveness of the In2Care station, which suppresses mosquito populations via the entomopathogenic fungus, Beauveria bassiana, and the insect growth regulator pyriproxyfen, the latter of which is autodisseminated among larval habitats by contaminated mosquitoes. A field trial was conducted using 110 In2Care stations in a 50,000 m2 area and results were compared to 4 control areas that did not receive the treatment. Efficacy was evaluated by comparing egg counts and measuring larvicidal impact in surrounding breeding sites. Laboratory experiments validated the effect of B. bassiana on adult survival. Results of this field trial indicate that, 6 wk after the In2Care stations were deployed, treatment site ovitraps contained 43% fewer eggs than control site ovitraps, and 33% fewer eggs after 10 wk, suggesting that the In2Care station was able to reduce the egg density of Ae. notoscriptus. Population reduction remained evident for up to 3 wk after In2Care stations were removed. Treatment site ovitraps had significantly fewer Ae. notoscriptus eclosing than control site ovitraps, confirming the pyriproxyfen autodissemination feature of the stations. An average reduction of 50% in adult eclosion was achieved. Exposure to B. bassiana resulted in four-times higher mortality among adult mosquitoes. Additionally, using fresh In2Care nettings led to an 88% decrease in average survival compared to 4-wk-old nettings. The use of In2Care stations has potential for suppressing Ae. notoscriptus egg density.

Inelastic electron backscattering in a generic helical edge channel.

We evaluate the low-temperature conductance of a weakly interacting one-dimensional helical liquid without axial spin symmetry. The lack of that symmetry allows for inelastic backscattering of a single electron, accompanied by forward scattering of another. This joint effect of weak interactions and potential scattering off impurities results in a temperature-dependent deviation from the quantized conductance, δG ∝ T4. In addition, δG is sensitive to the position of the Fermi level. We determine numerically the parameters entering our generic model for the Bernevig-Hughes-Zhang Hamiltonian of a HgTe/CdTe quantum well in the presence of Rashba spin-orbit coupling.

Estimating dispersal using close kin dyads: The kindisperse R package.

Investigating dispersal in animal populations can be difficult, particularly for taxa that are hard to directly observe such as those that are small or rare. A promising solution may come from new approaches that use genome-wide sequence data to detect close kin dyads and estimate dispersal parameters from the distribution of these dyads. These methods have so far only been applied to mosquito populations. However, they should have broad applicability to a range of taxa, although no assessment has yet been made on their performance under different dispersal conditions and study designs. Here we develop an R package and shiny app, kindisperse, that can be used to estimate dispersal parameters from the spatial distribution of close kin. kindisperse can handle study designs that target different life stages and allows for a range of dispersal kernel shapes and organismal life histories; we provide implementation examples for a vertebrate (Antechinus) and an invertebrate (Aedes). We use simulations run in kindisperse to compare the performance of two published close kin methodologies, showing that one method produces unbiased estimates whereas the other produces downward-biased estimates. We also use kindisperse simulations to investigate how study design affects dispersal estimates, and we provide guidelines for the size and shape of sample sites as well as the number of close kin needed for accurate estimates. kindisperse is easily adaptable for application to a variety of research contexts ranging from invasive pests to threatened species where noninvasive DNA sampling can be used to detect close kin.

Sex-specific distribution and classification of Wolbachia infections and mitochondrial DNA haplogroups in Aedes albopictus from the Indo-Pacific.

The arbovirus vector Aedes albopictus (Asian tiger mosquito) is common throughout the Indo-Pacific region, where most global dengue transmission occurs. We analysed population genomic data and tested for cryptic species in 160 Ae. albopictus sampled from 16 locations across this region. We found no evidence of cryptic Ae. albopictus but found multiple intraspecific COI haplotypes partitioned into groups representing three Asian lineages: East Asia, Southeast Asia and Indonesia. Papua New Guinea (PNG), Vanuatu and Christmas Island shared recent coancestry, and Indonesia and Timor-Leste were likely invaded from East Asia. We used a machine learning trained on morphologically sexed samples to classify sexes using multiple genetic features and then characterized the wAlbA and wAlbB Wolbachia infections in 664 other samples. The wAlbA and wAlbB infections as detected by qPCR showed markedly different patterns in the sexes. For females, most populations had a very high double infection incidence, with 67% being the lowest value (from Timor-Leste). For males, the incidence of double infections ranged from 100% (PNG) to 0% (Vanuatu). Only 6 females were infected solely by the wAlbA infection, while rare uninfected mosquitoes were found in both sexes. The wAlbA and wAlbB densities varied significantly among populations. For mosquitoes from Torres Strait and Vietnam, the wAlbB density was similar in single-infected and superinfected (wAlbA and wAlbB) mosquitoes. There was a positive association between wAlbA and wAlbB infection densities in superinfected Ae. albopictus. Our findings provide no evidence of cryptic species of Ae. albopictus in the region and suggest site-specific factors influencing the incidence of Wolbachia infections and their densities. We also demonstrate the usefulness of ddRAD tag depths as sex-specific mosquito markers. The results provide baseline data for the exploitation of Wolbachia-induced cytoplasmic incompatibility (CI) in dengue control.

Genome-wide SNPs of vegetable leafminer, Liriomyza sativae: Insights into the recent Australian invasion.

Liriomyza sativae, the vegetable leafminer, is an important agricultural pest originally from the Americas, which has now colonized all continents except Antarctica. In 2015, L. sativae arrived on the Australian mainland and established on the Cape York Peninsula in the northeast of the country near the Torres Strait, which provides a possible pathway for pests to enter Australia and evade biosecurity efforts. Here, we assessed genetic variation in L. sativae based on genome-wide single nucleotide polymorphisms (SNPs) generated by double digest restriction-site-associated DNA sequencing (ddRAD-seq), aiming to uncover the potential origin(s) of this pest in Australia and contribute to reconstructing its global invasion history. Our fineRADstructure results and principal component analysis suggest Australian mainland populations were genetically close to populations from the Torres Strait, whereas populations from Asia, Africa, and Papua New Guinea (PNG) were more distantly related. Hawaiian populations were genetically distinct from all other populations of L. sativae included in our study. Admixture analyses further revealed that L. sativae from the Torres Strait may have genetic variation originating from multiple sources including Indonesia and PNG, and which has now spread to the Australian mainland. The L. sativae lineages from Asia and Africa appear closely related. Isolation-by-distance (IBD) was found at a broad global scale, but not within small regions, suggesting that human-mediated factors likely contribute to the local spread of this pest. Overall, our findings suggest that an exotic Liriomyza pest invaded Australia through the Indo-Papuan conduit, highlighting the importance of biosecurity programs aimed at restricting the movement of pests and diseases through this corridor.