Assessing the potential of genotyping-by-sequencing-derived single nucleotide polymorphisms to identify the geographic origins of intercepted gypsy moth (Lymantria dispar) specimens: A proof-of-concept study

Employing hydrogen stable isotope analysis for traceability of invasive pine wood nematode vector insects

BACKGROUND

Bursaphelenchus xylophilus, the pathogen responsible for pine wilt disease, significantly threatens pine forests globally, with rapid infection leading to tree mortality within 40 days. This disease spreads efficiently through vector beetles, with Monochamus alternatus and Monochamus saltuarius serving as the primary vectors and Monochamus galloprovincialis being considered a potential vector in China. Molecular traceability and stable isotope traceability are both important techniques for insect tracing, each with its own advantages and limitations. Molecular traceability, such as DNA-based methods, offers high specificity and accuracy in identifying insect species, but it may require more complex procedures and higher costs due to the need for advanced equipment and skilled personnel. Stable isotope tracing, however, provides conservative and high-precision signatures, and is label-free, making it suitable for tracing the origins and pathways of insects in complex environments. Using the method of degreasing, the fat of the sample epidermis was removed, only the bone part was left, and the abdominal cavity was emptied to prevent the food from affecting the determination of hydrogen isotope. The present study aimed to trace the geographic origins of these vector insects in China using hydrogen stable isotope ratios, providing a novel tool for disease management.

RESULTS

Data processing and modeling were conducted using R and Origin software, yielding the equation y = 0.9822x + 19.1765. Statistical tests confirmed the model's significance and fit. Results indicate a significant positive correlation of temperature and precipitation with atmospheric hydrogen isotopic composition.

CONCLUSION

Our findings demonstrate the model's applicability in tracking the spread of pine wilt vectors, offering a promising approach to forest pathogen management. If the sample isotope test value is brought into the model and marked in the partition diagram, it can be proved that the model can effectively identify the source information of the sample. © 2025 Society of Chemical Industry.

Tracking the North American Asian Longhorned Beetle Invasion With Genomics

Biological invasions pose significant threats to ecological and economic stability, with invasive pests like the Asian longhorned beetle (Anoplophora glabripennis Motschulsky, ALB) causing substantial damage to forest ecosystems. Effective pest management relies on comprehensive knowledge of the insect's biology and invasion history. This study uses genomics to address these knowledge gaps and inform existing biosurveillance frameworks. We used 2768 genome-wide single nucleotide polymorphisms to compare invasive A. glabripennis populations in North America, using genomic variation to trace their sources of invasion and spread patterns, thereby refining our understanding of this species' invasion history. We found that most North American A. glabripennis infestations were distinct, resulting from multiple independent introductions from the native range. Following their introduction, all invasive populations experienced a genetic bottleneck which was followed by a population expansion, with a few also showing secondary spread to satellite infestations. Our study provides a foundation for a genome-based biosurveillance tool that can be used to clarify the origin of intercepted individuals, allowing regulatory agencies to strengthen biosecurity measures against this invasive beetle.

Complex European invasion history of Anoplophora glabripennis (Motschulsky): new insights in its population genomic differentiation using genotype-by-sequencing

Anthropogenic activities like trade facilitate increasing rates of biological invasions. Asian long-horned beetle (ALB), which is naturally distributed in eastern Asia (China, Korean peninsula), was introduced via wood packing materials (WPM) used in trade to North America (1996) and Europe (2001). We used 7810 single nucleotide polymorphisms (SNPs) derived by a genotype-by-sequencing (GBS) approach to decipher the introduction patterns into Europe. This is applied for the first time on European ALB outbreaks from Germany, Switzerland, and Italy, both from still active and already eradicated infestations. The genome-wide SNPs detected signs of small and highly structured populations within Europe, showing clear founder effects. The very high population differentiation is presumably derived from multiple independent introductions to Europe, which are spatially restricted in mating. By admixture and phylogenetic analyses, some cases of secondary dispersal were observed. Furthermore, some populations suggest admixture, which might have been originated by either multiple introductions from different sources into the new sites or recurrent introductions from an admixed source population. Our results confirmed a complex invasion history of the ALB into Europe and the usability of GBS obtained SNPs in invasion science even without source populations.

Lymantria Dispar Iflavirus 1 RNA Comprises a Large Proportion of RNA in Adult L. dispar Moths

The spongy moth virus Lymantria dispar iflavirus 1 (LdIV1), originally identified from a Lymantria dispar cell line, was detected in 24 RNA samples from female moths of four populations from the USA and China. Genome-length contigs were assembled for each population and compared with the reference genomes of the first reported LdIV1 genome (Ames strain) and two LdIV1 sequences available in GenBank originating from Novosibirsk, the Russian Federation. A whole-genome phylogeny was generated for these sequences, indicating that LdIV1 viruses observed in North American (flightless) and Asian (flighted) spongy moth lineages indeed partition into clades as would be expected per their host's geographic origin and biotype. A comprehensive listing of synonymous and non-synonymous mutations, as well as indels, among the polyprotein coding sequences of these seven LdIV1 variants was compiled and a codon-level phylogram was computed using polyprotein sequences of these, and 50 additional iflaviruses placed LdIV1 in a large clade consisting mostly of iflaviruses from other species of Lepidoptera. Of special note, LdIV1 RNA was present at very high levels in all samples, with LdIV1 reads accounting for a mean average of 36.41% (ranging from 1.84% to 68.75%, with a standard deviation of 20.91) of the total sequenced volume.

Analyses of adult transcriptomes from four different populations of the spongy moth, Lymantria dispar L., from China and the USA

The spongy moth Lymantria dispar, formerly known as the gypsy moth, is a forest pest that occurs as two different biotypes: the European spongy moth (ESM), Lymantria dispar dispar, which is distributed in Europe and North America; and the Asian spongy moth (ASM), which consists of subspecies Lymantria dispar asiatica and Lymantria dispar japonica and is distributed in China, Russia, Korea, and Japan. The Asian biotype is classified as a quarantine pest by the U.S. Department of Agriculture because of the superior flight ability of adult females compared to females of the European biotype. To identify genes that potentially account for differences in female flight capability between the two biotypes, we assembled and compared transcriptional profiles of two North American populations of ESM and two Chinese populations of ASM, including samples of unmated female adults and females after mating and oviposition. Of 129,286 unigenes identified, 306 were up-regulated in ASM samples relative to ESM, including genes involved in egg production. In contrast, 2309 unigenes were down-regulated in ASM samples, including genes involved in energy production. Although a previous study found that ASM female flight was reduced after oviposition, a comparison of gene expression before and after mating and oviposition did not reveal any genes which were consistently up- or down-regulated in the two ASM populations.

Lymantria dispar (L.) (Lepidoptera: Erebidae): Current Status of Biology, Ecology, and Management in Europe with Notes from North America

The European Spongy moth, Lymantria dispar (L.) (Lepidoptera: Erebidae), is an abundant species found in oak woods in Central and Southern Europe, the Near East, and North Africa and is an important economic pest. It is a voracious eater and can completely defoliate entire trees; repeated severe defoliation can add to other stresses, such as weather extremes or human activities. Lymantria dispar is most destructive in its larval stage (caterpillars), stripping away foliage from a broad variety of trees (>500 species). Caterpillar infestation is an underestimated problem; medical literature reports that established populations of caterpillars may cause health problems to people and animals. Inflammatory reactions may occur in most individuals after exposure to setae, independent of previous exposure. Currently, chemical and mechanical methods, natural predators, and silvicultural practices are included for the control of this species. Various insecticides have been used for its control, often through aerial sprayings, which negatively affect biodiversity, frequently fail, and are inappropriate for urban/recreational areas. However, bioinsecticides based on various microorganisms (e.g., entomopathogenic viruses, bacteria, and fungi) as well as technologies such as mating disruption using sex pheromone traps have replaced insecticides for the management of L. dispar.

Tracking global invasion pathways of the spongy moth (Lepidoptera: Erebidae) to the United States using stable isotopes as endogenous biomarkers

The spread of invasive insect species causes enormous ecological damage and economic losses worldwide. A reliable method that tracks back an invaded insect's origin would be of great use to entomologists, phytopathologists, and pest managers. The spongy moth (Lymantria dispar, Linnaeus 1758) is a persistent invasive pest in the Northeastern United States and periodically causes major defoliations in temperate forests. We analyzed field-captured (Europe, Asia, United States) and laboratory-reared L. dispar specimens for their natal isotopic hydrogen and nitrogen signatures imprinted in their biological tissues (δ2H and δ15N) and compared these values to the long-term mean δ2H of regional precipitation (Global Network of Isotopes in Precipitation) and δ15N of regional plants at the capture site. We established the percentage of hydrogen-deuterium exchange for L. dispar tissue (Pex = 8.2%) using the comparative equilibration method and two-source mixing models, which allowed the extraction of the moth's natal δ2H value. We confirmed that the natal δ2H and δ15N values of our specimens are related to the environmental signatures at their geographic origins. With our regression models, we were able to isolate potentially invasive individuals and give estimations of their geographic origin. To enable the application of these methods on eggs, we established an egg-to-adult fraction factor for L. dispar (Δegg-adult = 16.3 ± 4.3‰). Our models suggested that around 25% of the field-captured spongy moths worldwide were not native in the investigated capture sites. East Asia was the most frequently identified location of probable origin. Furthermore, our data suggested that eggs found on cargo ships in the United States harbors in Alaska, California, and Louisiana most probably originated from Asian L. dispar in East Russia. These findings show that stable isotope biomarkers give a unique insight into invasive insect species pathways, and thus, can be an effective tool to monitor the spread of insect pest epidemics.

Genome-scale phylogeography resolves the native population structure of the Asian longhorned beetle, Anoplophora glabripennis (Motschulsky)

Human-assisted movement has allowed the Asian longhorned beetle (ALB, Anoplophora glabripennis (Motschulsky)) to spread beyond its native range and become a globally regulated invasive pest. Within its native range of China and the Korean peninsula, human-mediated dispersal has also caused cryptic translocation of insects, resulting in population structure complexity. Previous studies used genetic methods to detangle this complexity but were unable to clearly delimit native populations which is needed to develop downstream biosurveillance tools. We used genome-wide markers to define historical population structure in native ALB populations and contemporary movement between regions. We used genotyping-by-sequencing to generate 6102 single-nucleotide polymorphisms (SNPs) and amplicon sequencing to genotype 53 microsatellites. In total, we genotyped 712 individuals from ALB's native distribution. We observed six distinct population clusters among native ALB populations, with a clear delineation between northern and southern groups. Most of the individuals from South Korea were distinct from populations in China. Our results also indicate historical divergence among populations and suggest limited large-scale admixture, but we did identify a restricted number of cases of contemporary movement between regions. We identified SNPs under selection and describe a clinal allele frequency pattern in a missense variant associated with glycerol kinase, an important enzyme in the utilization of an insect cryoprotectant. We further demonstrate that small numbers of SNPs can assign individuals to geographic regions with high probability, paving the way for novel ALB biosurveillance tools.

Potential Differences and Methods of Determining Gypsy Moth Female Flight Capabilities: Implications for the Establishment and Spread in Novel Habitats

The introduction of the Asian gypsy moth into novel environments continues with frequent interceptions in North America. There is a concern that these subspecies will pose a greater threat to the forests and urban environments of North America than the established gypsy moths (Lymantria dispardispar L.), due to their greater capacity for female flight. Asian gypsy moth populations vary in many key traits, including female flight capabilities. The potential impacts of female flight, in combination with the other key traits, on the ecology and spread of this insect are first discussed in this communication. This also provides the first review of most of the current literature on the variations in flight capability and flight distance of gypsy moth populations, as well as variation in other traits of concern and the potential methods of identification, with special attention paid to the Asian subspecies Lymantria dispar japonica Motschulsky and Lymantria dispar asiatica Vinkovskij. There are currently good tools for identifying the general origin of introduced gypsy moth populations, but these do not provide enough information to effectively manage introductions. Gypsy moth key traits differ among populations, even within each subspecies of the gypsy moth, so introduction of gypsy moths from other world areas into locations where the gypsy moth is already present could result in unwanted changes in gypsy moth biology. It also appears that the introduction of flight-capable females could enhance a population’s dispersal capability and require modifications to management protocols used for flightless females. Therefore, rapid tools to assess key traits in introduced populations are needed to adequately plan for, or deal with, new introductions into novel habitats.

Genetic Rescue and the Plight of Ponui Hybrids

Long-term sustainable and resilient populations is a key goal of conservation. How to best achieve this is controversial. There are, for instance, polarized views concerning the fitness and conservation value of hybrid populations founded through multi-origin translocations. A classic example concerns Apteryx (kiwi) in New Zealand. The A. mantelli of Ponui Island constitute a hybrid population where the birds are highly successful in their island habitat. A key dilemma for managers is understanding the reason for this success. Are the hybrid birds of Ponui Island of “no future conservation value” as recently asserted, or do they represent an outstanding example of genetic rescue and an important resource for future translocations? There has been a paradigm shift in scientific thinking concerning hybrids, but the ecological significance of admixed genomes remains difficult to assess. This limits what we can currently predict in conservation science. New understanding from genome science challenges the sufficiency of population genetic models to inform decision making and suggests instead that the contrasting outcomes of hybridization, “outbreeding depression” and “heterosis,” require understanding additional factors that modulate gene and protein expression and how these factors are influenced by the environment. We discuss these findings and the investigations that might help us to better understand the birds of Ponui, inform conservation management of kiwi and provide insight relevant for the future survival of Apteryx.

Disjunction between canola distribution and the genetic structure of its recently described pest, the canola flower midge (Contarinia brassicola)

Population genomics is a useful tool to support integrated pest management as it can elucidate population dynamics, demography, and histories of invasion. Here, we use a restriction site-associated DNA sequencing approach combined with whole-genome amplification (WGA) to assess genomic population structure of a newly described pest of canola, the diminutive canola flower midge, Contarinia brassicola. Clustering analyses recovered little geographic structure across the main canola production region but differentiated several geographically disparate populations at edges of the agricultural zone. Given a lack of alternative hypotheses for this pattern, we suggest these data support alternative hosts for this species and thus our canola-centric view of this midge as a pest has limited our understanding of its biology. These results speak to the need for increased surveying efforts across multiple habitats and other potential hosts within Brassicaceae to improve both our ecological and evolutionary knowledge of this species and contribute to effective management strategies. We additionally found that use of WGA prior to library preparation was an effective method for increasing DNA quantity of these small insects prior to restriction site-associated DNA sequencing and had no discernible impact on genotyping consistency for population genetic analysis; WGA is therefore likely to be tractable for other similar studies that seek to randomly sample markers across the genome in small organisms.

Mycorrhiza: genotype assignment using phylogenetic networks

MOTIVATION

The genotype assignment problem consists of predicting, from the genotype of an individual, which of a known set of populations it originated from. The problem arises in a variety of contexts, including wildlife forensics, invasive species detection and biodiversity monitoring. Existing approaches perform well under ideal conditions but are sensitive to a variety of common violations of the assumptions they rely on.

RESULTS

In this article, we introduce Mycorrhiza, a machine learning approach for the genotype assignment problem. Our algorithm makes use of phylogenetic networks to engineer features that encode the evolutionary relationships among samples. Those features are then used as input to a Random Forests classifier. The classification accuracy was assessed on multiple published empirical SNP, microsatellite or consensus sequence datasets with wide ranges of size, geographical distribution and population structure and on simulated datasets. It compared favorably against widely used assessment tests or mixture analysis methods such as STRUCTURE and Admixture, and against another machine-learning based approach using principal component analysis for dimensionality reduction. Mycorrhiza yields particularly significant gains on datasets with a large average fixation index (FST) or deviation from the Hardy-Weinberg equilibrium. Moreover, the phylogenetic network approach estimates mixture proportions with good accuracy.

AVAILABILITY AND IMPLEMENTATION

Mycorrhiza is released as an easy to use open-source python package at github.com/jgeofil/mycorrhiza.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genomic biosurveillance of forest invasive alien enemies: A story written in code

The world's forests face unprecedented threats from invasive insects and pathogens that can cause large irreversible damage to the ecosystems. This threatens the world's capacity to provide long-term fiber supply and ecosystem services that range from carbon storage, nutrient cycling, and water and air purification, to soil preservation and maintenance of wildlife habitat. Reducing the threat of forest invasive alien species requires vigilant biosurveillance, the process of gathering, integrating, interpreting, and communicating essential information about pest and pathogen threats to achieve early detection and warning and to enable better decision-making. This process is challenging due to the diversity of invasive pests and pathogens that need to be identified, the diverse pathways of introduction, and the difficulty in assessing the risk of establishment. Genomics can provide powerful new solutions to biosurveillance. The process of invasion is a story written in four chapters: transport, introduction, establishment, and spread. The series of processes that lead to a successful invasion can leave behind a DNA signature that tells the story of an invasion. This signature can help us understand the dynamic, multistep process of invasion and inform management of current and future introductions. This review describes current and future application of genomic tools and pipelines that will provide accurate identification of pests and pathogens, assign outbreak or survey samples to putative sources to identify pathways of spread, and assess risk based on traits that impact the outbreak outcome.

Range-wide population genomics of the Mexican fruit fly: Toward development of pathway analysis tools

Recurrently invading pests provide unique challenges for pest management, but also present opportunities to utilize genomics to understand invasion dynamics and inform regulatory management through pathway analysis. In the southern United States, the Mexican fruit fly Anastrepha ludens is such a pest, and its incursions into Texas and California represent major threats to the agricultural systems of those regions. We developed a draft genome assembly for A. ludens, conducted range-wide population genomics using restriction site-associated DNA sequencing, and then developed and demonstrated a panel of highly differentiated diagnostic SNPs for source determination of intercepted flies in this system. Using 2,081 genomewide SNPs, we identified four populations across the range of A. ludens, corresponding to western Mexico, eastern Mexico/Texas, Guatemala/Belize/Honduras, and Costa Rica/Panama, with some intergradation present between clusters, particularly in Central America. From this population genomics framework, we developed a diagnostic panel of 28 highly differentiated SNPs that were able to recreate the genomewide population structure in this species. We demonstrated this panel on a set of test specimens, including specimens intercepted as part of regular trapping surveillance in Texas and California, and we were able to predict populations of origin for these specimens. This methodology presents a highly applied use of genomic techniques and can be implemented in any group of recurrently invading pests.

Genetic evidence of broad spreading of Lymantria dispar in the West Siberian Plain

Gypsy moth Lymantria dispar L. 1758 (Lepidoptera: Erebidae) is one of the most dangerous forest pests of the Holarctic region. Outbreaks of gypsy moth populations lead to significant defoliation of local forests. Within the vast territory of the West Siberian Plain, we noted the outbreak front movement in the north-east direction with a speed 100-200 km per year. The reason for the outbreak's movement is still unclear because L. dispar females are characterised by limited flight ability, which is not enough to support that movement per se. Herein, we analysed the mtDNA divergence pattern among L. dispar populations collected from the vast territory of the West Siberian Plain to determine the boundaries of populations and reveal the effect of the outbreak's front movement on mtDNA patterns of populations. The 590-bp region of the cytochrome oxidase subunit I gene of the mitochondrial genome was sequenced for 220 specimens that were collected from 18 localities along a transect line (~ 1400 km). Our results clearly show that the gypsy moth populations of the vast Siberian territory are not subdivided. This result can be explained by extensive genetic exchange among local populations. Taking into account that the flight ability of L. dispar females is rather limited, we suggest that spreading occurs through ballooning of early instar larvae. This hypothesis was confirmed by the coincidence of the outbreaks' movement direction with that of the dominant winds, complemented by the observation of ballooned larvae far from a forest edge.

Evolutionary genomics of gypsy moth populations sampled along a latitudinal gradient

The European gypsy moth (Lymantria dispar L.) was first introduced to Massachusetts in 1869 and within 150 years has spread throughout eastern North America. This large-scale invasion across a heterogeneous landscape allows examination of the genetic signatures of adaptation potentially associated with rapid geographical spread. We tested the hypothesis that spatially divergent natural selection has driven observed changes in three developmental traits that were measured in a common garden for 165 adult moths sampled from six populations across a latitudinal gradient covering the entirety of the range. We generated genotype data for 91,468 single nucleotide polymorphisms based on double digest restriction-site associated DNA sequencing and used these data to discover genome-wide associations for each trait, as well as to test for signatures of selection on the discovered architectures. Genetic structure across the introduced range of gypsy moth was low in magnitude (F_ST  = 0.069), with signatures of bottlenecks and spatial expansion apparent in the rare portion of the allele frequency spectrum. Results from applications of Bayesian sparse linear mixed models were consistent with the presumed polygenic architectures of each trait. Further analyses indicated spatially divergent natural selection acting on larval development time and pupal mass, with the linkage disequilibrium component of this test acting as the main driver of observed patterns. The populations most important for these signals were two range-edge populations established less than 30 generations ago. We discuss the importance of rapid polygenic adaptation to the ability of non-native species to invade novel environments.

Global phylogeography reveals the origin and the evolutionary history of the gypsy moth (Lepidoptera, Erebidae)

We examined the global phylogeography of gypsy moth (Lymantria dispar L.) using molecular data based on mitochondrial and nuclear genes. Populations from all biogeographic regions of the native and introduced range of L. dispar, were sampled to fully document intraspecific and subspecies variation, identify potential cryptic species, and to clarify the relationships among major phylogeographic lineages. We recovered three major mtDNA lineages of L. dispar: Transcaucasia; East Asia + Japan; and Europe + Central Asia. The circumscription of these lineages is only partially consistent with the current taxonomic concept (i.e., L. dispar dispar; L. dispar asiatica; L. dispar japonica), with the following important discrepancies: (1) north-central Asian populations, including topotypical populations of L. dispar asiatica, may be more closely related to European rather than Asian segregates, which would require the synonymization of the taxon asiatica and establishment of a new name; (2) the Japanese populations (L. d. japonica) are not distinct from east Asian populations; (3) the presence of a distinct, unnamed mitogenomic lineage endemic to the Trancaucasus region. We demonstrated that the population from Transcaucasia contains the highest mitochondrial haplotype diversity among L. dispar, potentially indicative of an ancestral area for the entire dispar-group. Our study corroborates the endemic Hokkaido, Japan taxon Lymantria umbrosa (Butler) as the sister group to all other L. dispar populations, but the applicability of the names umbrosa versus hokkaidoensis Goldschmidt needs to be re-evaluated. The ancestral area analysis suggest that Japan was likely colonized via Sakhalin ∼1 Mya, in contrast to previous studies which have suggested colonization of the Japanese archipelago via the Korean Peninsula. Lastly, mitogenomic variation within L. dispar is incongruent with phylogenies based on nuclear DNA, as nDNA gene phylogenies did not recover the three major mtDNA lineages, and also failed to recover L. dispar and L. umbrosa as reciprocally monophyletic.

Insights into the Structure of the Spruce Budworm (Choristoneura fumiferana) Genome, as Revealed by Molecular Cytogenetic Analyses and a High-Density Linkage Map

Genome structure characterization can contribute to a better understanding of processes such as adaptation, speciation, and karyotype evolution, and can provide useful information for refining genome assemblies. We studied the genome of an important North American boreal forest pest, the spruce budworm, Choristoneura fumiferana, through a combination of molecular cytogenetic analyses and construction of a high-density linkage map based on single nucleotide polymorphism (SNP) markers obtained through a genotyping-by-sequencing (GBS) approach. Cytogenetic analyses using fluorescence in situ hybridization methods confirmed the haploid chromosome number of n = 30 in both sexes of C. fumiferana and showed, for the first time, that this species has a WZ/ZZ sex chromosome system. Synteny analysis based on a comparison of the Bombyx mori genome and the C. fumiferana linkage map revealed the presence of a neo-Z chromosome in the latter species, as previously reported for other tortricid moths. In this neo-Z chromosome, we detected an ABC transporter C2 (ABCC2) gene that has been associated with insecticide resistance. Sex-linkage of the ABCC2 gene provides a genomic context favorable to selection and rapid spread of resistance against Bacillus thuringiensis serotype kurstaki (Btk), the main insecticide used in Canada to control spruce budworm populations. Ultimately, the linkage map we developed, which comprises 3586 SNP markers distributed over 30 linkage groups for a total length of 1720.41 cM, will be a valuable tool for refining our draft assembly of the spruce budworm genome.

Comparative analysis of mitochondrial genomes of geographic variants of the gypsy moth, Lymantria dispar, reveals a previously undescribed genotypic entity

The gypsy moth, Lymantria dispar L., is one of the most destructive forest pests in the world. While the subspecies established in North America is the European gypsy moth (L. dispar dispar), whose females are flightless, the two Asian subspecies, L. dispar asiatica and L. dispar japonica, have flight-capable females, enhancing their invasiveness and warranting precautionary measures to prevent their permanent establishment in North America. Various molecular tools have been developed to help distinguish European from Asian subspecies, several of which are based on the mitochondrial barcode region. In an effort to identify additional informative markers, we undertook the sequencing and analysis of the mitogenomes of 10 geographic variants of L. dispar, including two or more variants of each subspecies, plus the closely related L. umbrosa as outgroup. Several regions of the gypsy moth mitogenomes displayed nucleotide substitutions with potential usefulness for the identification of subspecies and/or geographic origins. Interestingly, the mitogenome of one geographic variant displayed significant divergence relative to the remaining variants, raising questions about its taxonomic status. Phylogenetic analyses placed this population from northern Iran as basal to the L. dispar clades. The present findings will help improve diagnostic tests aimed at limiting risks of AGM invasions.