Diapause survival requires a temperature-sensitive preparatory period.

Diapause is a form of internally-controlled dormancy that allows insects to avoid stressful conditions and periods of low food availability. Eastern spruce budworm (Choristoneura fumiferana Clemens), like many cold-adapted insects, enter diapause well in advance of winter conditions, thus exposing them to elevated temperatures during fall that can deplete energy stores and impact post-diapause survival. We explored the impact of fall conditions on C. fumiferana by manipulating the length of the fall period and exposure temperatures during the diapause initiation phase of second instar larvae in a factorial design. We exposed second instar larvae to four fall temperatures (10, 15, 20, and 25°C) and five exposure times (1, 2, 4, 6, and 10 weeks) prior to standardized diapause conditions. We measured metabolites (glycogen, glycerol, and protein) prior to and during diapause for a subset of individuals. We also measured post-diapause survival by quantifying emergence following diapause conditions for a subset of individuals. We found that long, warm fall conditions depleted glycogen content and lowered post-diapause survival. We also found that short, cool conditions impacted post-diapause survival, although glycogen content remained high. Our results showed that fall conditions have substantial fitness consequences to overwintering insects. Optimal fall conditions struck a balance between exposure time and temperature. Our findings point to a potentially adaptive reason for early diapause onset: that an undescribed, but temperature-sensitive process is occurring in C. fumiferana larvae during the diapause initiation period that is essential for overwintering survival and successful post-diapause emergence.

Transcriptome profiling of Fraxinus excelsior genotypes infested by emerald ash borer.

European ash, Fraxinus excelsior is facing the double threat of ongoing devastation by the invasive fungal pathogen, Hymenoscyphus fraxineus and the imminent arrival of the non-native emerald ash borer (EAB), Agrilus planipennis. The spread of EAB which is currently moving westwards from European Russia and Ukraine into central Europe, poses an additional substantial threat to European ash, F. excelsior. While the molecular basis for resistance or variation in resistance among European ash genotypes is heavily investigated, comparatively little is known about the molecular ash traits involved in resistance against EAB. In this study we have gathered transcriptomic data from EAB inoculated genotypes of F. excelsior that have previously shown different levels of susceptibility to EAB. Resultant datasets show differential gene expression in susceptible and resistant genotypes in response to EAB infestation. This data will provide important information on the molecular basis of resistance to the EAB and allow the development of management plans to combat a pending threat of a culturally and ecologically important European tree species.

A glimmer of hope - ash genotypes with increased resistance to ash dieback pathogen show cross-resistance to emerald ash borer.

Plants rely on cross-resistance traits to defend against multiple, phylogenetically distinct enemies. These traits are often the result of long co-evolutionary histories. Biological invasions can force naïve plants to cope with novel, coincident pests, and pathogens. For example, European ash (Fraxinus excelsior) is substantially threatened by the emerald ash borer (EAB), Agrilus planipennis, a wood-boring beetle, and the ash dieback (ADB) pathogen, Hymenoscyphus fraxineus. Yet, plant cross-resistance traits against novel enemies are poorly explored and it is unknown whether naïve ash trees can defend against novel enemy complexes via cross-resistance mechanisms. To gain mechanistic insights, we quantified EAB performance on grafted replicates of ash genotypes varying in ADB resistance and characterized ash phloem chemistry with targeted and untargeted metabolomics. Emerald ash borer performed better on ADB-susceptible than on ADB-resistant genotypes. Moreover, changes in EAB performance aligned with differences in phloem chemical profiles between ADB-susceptible and ADB-resistant genotypes. We show that intraspecific variation in phloem chemistry in European ash can confer increased cross-resistance to invasive antagonists from different taxonomic kingdoms. Our study suggests that promotion of ADB-resistant ash genotypes may simultaneously help to control the ADB disease and reduce EAB-caused ash losses, which may be critical for the long-term stability of this keystone tree species.

Molecular signatures of diapause in the Asian longhorned beetle: Gene expression.

Most previous studies on gene expression during insect diapause do not address among-tissue variation in physiological processes. We measured transcriptomic changes during larval diapause in the Asian longhorned beetle, Anoplophora glabripennis (Coleoptera: Cerambycidae). We conducted RNA-seq on fat body, the supraesophageal ganglion, midgut, hindgut, and Malpighian tubules during pre-diapause, diapause maintenance, post-diapause quiescence, and post-diapause development. We observed a small, but consistent, proportion of genes within each gene expression profile that were shared among tissues, lending support for a core set of diapause-associated genes whose expression is tissue-independent. We evaluated the overarching hypotheses that diapause would be associated with cell cycle arrest, developmental arrest, and increased stress tolerance and found evidence of repressed TOR and insulin signaling, reduced cell cycle activity and increased capacity of stress response via heat shock protein expression and remodeling of the cytoskeleton. However, these processes varied among tissues, with the brain and fat body appearing to maintain higher levels of cellular activity during diapause than the midgut or Malpighian tubules. We also observed temperature-dependent changes in gene expression during diapause maintenance, particularly in genes related to the heat shock response and MAPK, insulin, and TOR signaling pathways. Additionally, we provide evidence for epigenetic reorganization during the diapause/post-diapause quiescence transition and expression of genes involved in post-translational modification, highlighting the need for investigations of the protein activity of these candidate genes and processes. We conclude that diapause development is coordinated via diverse tissue-specific gene expression profiles and that canonical diapause phenotypes vary among tissues.

Host association, environment, and geography underlie genomic differentiation in a major forest pest.

Diverse geographic, environmental, and ecological factors affect gene flow and adaptive genomic variation within species. With recent advances in landscape ecological modelling and high-throughput DNA sequencing, it is now possible to effectively quantify and partition their relative contributions. Here, we use landscape genomics to identify determinants of genomic differentiation in the forest tent caterpillar, Malacosoma disstria, a widespread and irruptive pest of numerous deciduous tree species in North America. We collected larvae from multiple populations across Eastern Canada, where the species experiences a diversity of environmental gradients and feeds on a number of different host tree species, including trembling aspen (Populus tremuloides), sugar maple (Acer saccharum), red oak (Quercus rubra), and white birch (Betula papyrifera). Using a combination of reciprocal causal modelling (RCM) and distance-based redundancy analyses (dbRDA), we show that differentiation of thousands of genome-wide single nucleotide polymorphisms (SNPs) among individuals is best explained by a combination of isolation by distance, isolation by environment (spatial variation in summer temperatures and length of the growing season), and differences in host association. Configuration of suitable habitat inferred from ecological niche models was not significantly related to genomic differentiation, suggesting that M. disstria dispersal is agnostic with respect to habitat quality. Although population structure was not discretely related to host association, our modelling framework provides the first molecular evidence of host-associated differentiation in M. disstria, congruent with previous documentation of reduced growth and survival of larvae moved between natal host species. We conclude that ecologically mediated selection is contributing to variation within M. disstria, and that divergent adaptation related to both environmental conditions and host association should be considered in ongoing research and management of this important forest pest.

Plasticity drives extreme cold tolerance of emerald ash borer (Agrilus planipennis) during a polar vortex.

Invasive species must often survive combinations of environmental conditions that differ considerably from their native range; however, for a given species it is unclear whether improved tolerance is the result of phenotypic plasticity or genetic adaptation (or both). Agrilus planipennis (Coleoptera: Buprestidae; the emerald ash borer) is an invasive pest of Fraxinus trees in North America and Europe. Previous studies in SW Ontario, Canada, showed that A. planipennis is freeze avoidant, preventing internal ice formation by accumulating Molar concentrations of glycerol in its hemolymph and depressing its supercooling point (SCP, the temperature at which it freezes). The cold tolerance of these SW Ontario animals was used to predict potential distribution, revealing that some Canadian cities should be too cold to allow populations to persist. However, a small population of A. planipennis has persisted in Winnipeg, Manitoba, Canada, through several severe 'polar vortex' events. In 2018/19, we collected A. planipennis larvae and prepupae from Winnipeg, MB and Southern Ontario, and found that individuals from Winnipeg were extremely cold tolerant - with SCPs as low as -52°C in prepupae (compared to -32°C in SW Ontario), and observed survival of unfrozen individuals exposed to -50°C for one hour. This cold tolerance was accompanied by higher hemolymph osmolality and glycerol concentration than in the SW Ontario individuals. To distinguish between phenotypic plasticity and local adaptation, in 2020/21 we overwintered Winnipeg-sourced individuals either outdoors in SW Ontario or in a simulated Winnipeg winter. Simulated Winnipeg winter individuals had cold tolerance similar to those overwintered in Winnipeg, while SW Ontario overwintered individuals had cold tolerance similar to those collected previously in the region. The simulated winter individuals had higher hemolymph glycerol concentrations than SW Ontario overwintered animals, at least in part due to greater dehydration. Thus, A. planipennis are cold-tolerant enough to survive some of the harshest winters where their host trees can grow, and most likely attain this cold tolerance via phenotypic plasticity. These findings raise the importance of delineating sensitivity of conclusions to unexpected phenotypic plasticity when predicting potential distributions of new invasives or responses to climate change.

The Spruce Budworm Genome: Reconstructing the Evolutionary History of Antifreeze Proteins.

Insects have developed various adaptations to survive harsh winter conditions. Among freeze-intolerant species, some produce "antifreeze proteins" (AFPs) that bind to nascent ice crystals and inhibit further ice growth. Such is the case of the spruce budworm, Choristoneura fumiferana (Lepidoptera: Tortricidae), a destructive North American conifer pest that can withstand temperatures below -30°C. Despite the potential importance of AFPs in the adaptive diversification of Choristoneura, genomic tools to explore their origins have until now been limited. Here we present a chromosome-scale genome assembly for C. fumiferana, which we used to conduct comparative genomic analyses aimed at reconstructing the evolutionary history of tortricid AFPs. The budworm genome features 16 genes homologous to previously reported C. fumiferana AFPs (CfAFPs), 15 of which map to a single region on chromosome 18. Fourteen of these were also detected in five congeneric species, indicating Choristoneura AFP diversification occurred before the speciation event that led to C. fumiferana. Although budworm AFPs were previously considered unique to the genus Choristoneura, a search for homologs targeting recently sequenced tortricid genomes identified seven CfAFP-like genes in the distantly related Notocelia uddmanniana. High structural similarity between Notocelia and Choristoneura AFPs suggests a common origin, despite the absence of homologs in three related tortricids. Interestingly, one Notocelia AFP formed the C-terminus of a "zonadhesin-like" protein, possibly representing the ancestral condition from which tortricid AFPs evolved. Future work should clarify the evolutionary path of AFPs between Notocelia and Choristoneura and assess the role of the "zonadhesin-like" protein as precursor of tortricid AFPs.

Plasticity of cold hardiness in the eastern spruce budworm, Choristoneura fumiferana.

High latitude insect populations must cope with extreme conditions, particularly low temperatures. Insects use a variety of cold hardiness mechanisms to withstand this temperature stress, and these can drive geographic distributions through overwintering mortality. The degree of cold hardiness can be altered by two evolved responses: phenotypic plasticity and local adaptation. Phenotypic plasticity can occur within or between generations (transgenerational plasticity; TGP), and local adaptation can evolve through directional selection in response to regional climatic differences. We used the eastern spruce budworm, Choristoneura fumiferana (Lepidoptera: Tortricidae) as a model to explore the role that variable winter temperatures play in inducing two aspects of plasticity in cold hardiness: TGP and local adaptation in phenotypic plasticity. This species is one of the most destructive boreal forest pests in North America, therefore accurately predicting overwintering survival is essential for effective management. While we found no evidence of TGP in cold hardiness, there was a long term fitness cost to larvae that experienced repeated cold exposures. We also found evidence of local adaptation in both seasonal and short-term plasticity of cold hardiness, as our more northerly populations that would experience lower overwintering temperatures had more plastic responses to cold exposure. These findings provide evidence for the importance of phenotypic plasticity and local adaptation when modelling species distributions.

Cold tolerance of laboratory-reared Asian longhorned beetles.

Low winter temperatures in temperate climates can limit the success of non-native species. The Asian longhorned beetle, Anoplophora glabripennis, is an invasive wood-boring pest of hardwood trees in North America and Europe. Native A. glabripennis populations are spread across several climate zones in China and the Korean Peninsula and are likely to encounter low temperatures in at least some of this range. Understanding the lethal limits of the overwintering life stages of A. glabripennis is essential for accurately modeling the risk that invasive populations pose to non-native environments. In this study, we provide the first systematic characterization of the cold tolerance strategy and lower lethal limits of A. glabripennis eggs, larvae, and pupae. In diapausing larvae, the most common overwintering stage in this species, we measure hemolymph glycerol and osmolality and identify the effects of prolonged low temperature exposure. In developing pupae, we identify sublethal effects caused by low temperature exposure before freezing. Eggs and larvae were the most cold-tolerant life stages; eggs were freeze-avoidant with an average supercooling point of -25.8 °C and larvae were freeze tolerant with an LT90 of -25 °C. Hemolymph osmolality of freeze-tolerant larvae, on average, increased to 811 mOsm during chilling. This increase was primarily driven by a concurrent, average increase of 232 mM hemolymph glycerol. Pupae died upon exposure to freezing temperatures, but accumulate strong sublethal effects prior to freezing, indicating that they are chill susceptible. Taken together, these data will be useful to inform species distribution modeling in A. glabripennis.

Surviving in a Frozen Forest: the Physiology of Eastern Spruce Budworm Overwintering.

The eastern spruce budworm, Choristoneura fumiferana, is one of North America's most destructive forest insects. It survives the harsh winters by deploying both a sophisticated diapause program and a complex suite of cryoprotective molecules. The spruce budworm's cryoprotective biochemistry could revolutionize organ storage and transplants. Here we review the latest in C. fumiferana overwintering physiology and identify emerging theoretical and practical questions that are open for exploration.

Dormancy in laboratory-reared Asian longhorned beetles, Anoplophora glabripennis.

An insect's capacity to survive winter is critical for range expansion in temperate regions. The Asian longhorned beetle (Anoplophora glabripennis) is a polyphagous wood-boring insect native to China and the Korean peninsula and poses a high risk of invasion in North America and Europe. It is unclear whether A. glabripennis enters diapause, which means that diapause cannot be included in assessments of the risk of this species invading forests in temperate regions. Using a laboratory colony, we examine larval developmental arrest, metabolic rates, gas exchange patterns, thermal sensitivity, and body composition to characterize larval dormancy. Chilled larvae entered a temperature-independent developmental arrest which usually required more than four weeks of chilling to break, decreased their metabolic rate by as much as 63%, and maintained energy stores throughout the chilling period - results consistent with an obligate diapause. We also observed a switch to discontinuous gas exchange at low temperatures. Thermal sensitivity of metabolic rate did not differ between chilled and non-chilled larvae. Taken together, we conclude that A. glabripennis enters a larval diapause during chilling and terminates diapause after a requisite chilling period. These results will enhance our ability to predict phenology and potential distribution of current and future invasions of A. glabripennis.

Assessing Niche Shifts and Conservatism by Comparing the Native and Post-Invasion Niches of Major Forest Invasive Species.

Invasive species experience biotic and abiotic conditions that may (or may not) resemble their native environment. We explored the methodology of determining climatic niches and compared the native and post-invasion niches of four invasive forest pests to determine if these species experienced shifts or changes in their new climatic niches. We used environmental principle components analysis (PCA-env) method to quantify climatic niche shifts, expansions, and temporal changes. Furthermore, we assessed the effect of variable selection in the delineation and comparison of niche space. We found that variable selection influenced the delineation and overlap of each niche, whereas the subset of climatic variables selected from the first two PCA-env axes explained more variance in environmental conditions than the complete set of climatic variables for all four species. Most focal species showed climatic niche shifts in their invasive range and had not yet fully occupied the available niche within the invaded range. Our species varied the proportion of niche overlap between the native and invasive ranges. By comparing native and invasive niches, we can help predict a species' potential range expansion and invasion potential. Our results can guide monitoring and help inform management of these and other invasive species.

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.

Origins and History of Laboratory Insect Stocks in a Multispecies Insect Production Facility, With the Proposal of Standardized Nomenclature and Designation of Formal Standard Names.

Laboratory insect colonies are an essential part of experimental insect science. Formalized naming of laboratory stocks is standard practice in model organisms such as mice and fruit flies, but crucial details such as colony origin and standard names are often lacking in nonmodel systems. For institutions involved in rearing multiple nonmodel species, effective monitoring requires standardized naming and nomenclature, from establishment to production, distribution, and publication. Insect rearing has been the cornerstone of the Insect Production and Quarantine Laboratories (IPQL) at the Great Lakes Forestry Centre for over 70 yr, but the histories of the insect colonies in this facility have not been adequately documented and formal, standardized names do not exist. We propose a standardized naming framework that we applied to the eight species reared at the IPQL to rectify these deficiencies. We also present the origin and history of each colony, essential information that is challenging to obtain post hoc. We suggest that other research institutions consider developing similar standards, so they can accurately document, communicate, and track laboratory insect their within the facilities and through the scientific literature.

Genetic and genomic evidence of niche partitioning and adaptive radiation in mountain pine beetle fungal symbionts.

Bark beetles form multipartite symbiotic associations with blue stain fungi (Ophiostomatales, Ascomycota). These fungal symbionts play an important role during the beetle's life cycle by providing nutritional supplementation, overcoming tree defences and modifying host tissues to favour brood development. The maintenance of stable multipartite symbioses with seemingly less competitive symbionts in similar habitats is of fundamental interest to ecology and evolution. We tested the hypothesis that the coexistence of three fungal species associated with the mountain pine beetle is the result of niche partitioning and adaptive radiation using SNP genotyping coupled with genotype-environment association analysis and phenotypic characterization of growth rate under different temperatures. We found that genetic variation and population structure within each species is best explained by distinct spatial and environmental variables. We observed both common (temperature seasonality and the host species) and distinct (drought, cold stress, precipitation) environmental and spatial factors that shaped the genomes of these fungi resulting in contrasting outcomes. Phenotypic intraspecific variations in Grosmannia clavigera and Leptographium longiclavatum, together with high heritability, suggest potential for adaptive selection in these species. By contrast, Ophiostoma montium displayed narrower intraspecific variation but greater tolerance to extreme high temperatures. Our study highlights unique phenotypic and genotypic characteristics in these symbionts that are consistent with our hypothesis. By maintaining this multipartite relationship, the bark beetles have a greater likelihood of obtaining the benefits afforded by the fungi and reduce the risk of being left aposymbiotic. Complementarity among species could facilitate colonization of new habitats and survival under adverse conditions.

Modeling Landscape-Level Spatial Variation in Sex Ratio Skew in the Mountain Pine Beetle (Coleoptera: Curculionidae).

Through their influence on effective population sizes, sex ratio skew affects population dynamics. We examined spatial variation in female-biased sex ratios in the mountain pine beetle (Dendroctonus ponderosae Hopkins) outbreak in western Canada to better understand how environmental context affects sex ratio skew. Our specific objectives were to: 1) characterize spatial variation in mountain pine beetle sex ratio; 2) test previously asserted hypotheses that beetle sex ratio varies with tree diameter and year in outbreak; and 3) develop predictive models of sex ratio skew for larval and adult populations. Using logistic regression, we modeled the probability that an individual beetle (n = 2,369) was female as a function of multiple environmental variables across 34 stands in British Columbia and Alberta, Canada. We identified a consistent female-biased sex ratio with significantly greater skew in adults (2:1, n = 713) than in larvae (1.76:1, n = 1,643). We found that the proportion of larval females increased with decreasing tree size and with outbreak age. However, adults did not respond to tree size and larvae did not respond to outbreak age. Predictive models differed between larvae and adults. All identified models perform well and included predictors related to weather, tree diameter, and year in outbreak. Female-biased sex ratios appear to originate from differential male mortality during development rather than from sex-biased oviposition, suggesting sex ratio skew is not the cause of outbreaks, but rather a consequence.

Population structure of mountain pine beetle symbiont Leptographium longiclavatum and the implication on the multipartite beetle-fungi relationships.

Over 18 million ha of forests have been destroyed in the past decade in Canada by the mountain pine beetle (MPB) and its fungal symbionts. Understanding their population dynamics is critical to improving modeling of beetle epidemics and providing potential clues to predict population expansion. Leptographium longiclavatum and Grosmannia clavigera are fungal symbionts of MPB that aid the beetle to colonize and kill their pine hosts. We investigated the genetic structure and demographic expansion of L. longiclavatum in populations established within the historic distribution range and in the newly colonized regions. We identified three genetic clusters/populations that coincide with independent geographic locations. The genetic profiles of the recently established populations in northern British Columbia (BC) and Alberta suggest that they originated from central and southern BC. Approximate Bayesian Computation supports the scenario that this recent expansion represents an admixture of individuals originating from BC and the Rocky Mountains. Highly significant correlations were found among genetic distance matrices of L. longiclavatum, G. clavigera, and MPB. This highlights the concordance of demographic processes in these interacting organisms sharing a highly specialized niche and supports the hypothesis of long-term multipartite beetle-fungus co-evolutionary history and mutualistic relationships.

Fitness dynamics within a poplar hybrid zone: I. Prezygotic and postzygotic barriers impacting a native poplar hybrid stand.

Hybridization and introgression are pervasive evolutionary phenomena that provide insight into the selective forces that maintain species boundaries, permit gene flow, and control the direction of evolutionary change. Poplar trees (Populus L.) are well known for their ability to form viable hybrids and maintain their distinct species boundaries despite this interspecific gene flow. We sought to quantify the hybridization dynamics and postzygotic fitness within a hybrid stand of balsam poplar (Populus balsamifera L.), eastern cottonwood (P. deltoides Marsh.), and their natural hybrids to gain insight into the barriers maintaining this stable hybrid zone. We observed asymmetrical hybrid formation with P. deltoides acting as the seed parent, but with subsequent introgression biased toward P. balsamifera. Native hybrids expressed fitness traits intermediate to the parental species and were not universally unfit. That said, native hybrid seedlings were absent from the seedling population, which may indicate additional selective pressures controlling their recruitment. It is imperative that we understand the selective forces maintaining this native hybrid zone in order to quantify the impact of exotic poplar hybrids on this native system.

Fitness dynamics within a poplar hybrid zone: II. Impact of exotic sex on native poplars in an urban jungle.

Trees bearing novel or exotic gene components are poised to contribute to the bioeconomy for a variety of purposes such as bioenergy production, phytoremediation, and carbon sequestration within the forestry sector, but sustainable release of trees with novel traits in large-scale plantations requires the quantification of risks posed to native tree populations. Over the last century, exotic hybrid poplars produced through artificial crosses were planted throughout eastern Canada as ornamentals or windbreaks and these exotics provide a proxy by which to examine the fitness of exotic poplar traits within the natural environment to assess risk of exotic gene escape, establishment, and spread into native gene pools. We assessed postzygotic fitness traits of native and exotic poplars within a naturally regenerated stand in eastern Canada (Quebec City, QC). Pure natives (P. balsamifera and P. deltoides spp. deltoides), native hybrids (P. deltoides × P. balsamifera), and exotic hybrids (trees bearing Populus nigra and P. maximowiczii genetic components) were screened for reproductive biomass, yield, seed germination, and fungal disease susceptibility. Exotic hybrids expressed fitness traits intermediate to pure species and were not significantly different from native hybrids. They formed fully viable seed and backcrossed predominantly with P. balsamifera. These data show that exotic hybrids were not unfit and were capable of establishing and competing within the native stand. Future research will seek to examine the impact of exotic gene regions on associated biotic communities to fully quantify the risk exotic poplars pose to native poplar forests.

Multi-locus species delimitation in closely related animals and fungi: one marker is not enough.

Despite taxonomy's 250-year history, the past 20 years have borne witness to remarkable advances in technology and techniques, as well as debate. DNA barcoding has generated a substantial proportion of this debate, with its proposition that a single mitochondrial sequence will consistently identify and delimit species, replacing more evidence-rich and time-intensive methods. Although mitochondrial DNA (mtDNA) has since been the focus of voluminous discussion and case studies, little effort has been made to comprehensively evaluate its success in delimiting closely related species. We have conducted the first broadly comparative literature review addressing the efficacy of molecular markers for delimiting such species over a broad taxonomic range. By considering only closely related species, we sought to avoid confusion of success rates with those due to deeply divergent taxa. We also address whether increased population-level or geographic sampling affects delimitation success. Based on the results from 101 studies, we found that all marker groups had approximately equal success rates (~70%) in delimiting closely related species and that the use of additional loci increased average delimitation success. We also found no relationship between increased sampling of intraspecific variability and delimitation success. Ultimately, our results support a multi-locus integrative approach to species delimitation and taxonomy.