The freeze-avoiding mountain pine beetle (Dendroctonus ponderosae) survives prolonged exposure to stressful cold by mitigating ionoregulatory collapse.

Insect performance is linked to environmental temperature, and surviving through winter represents a key challenge for temperate, alpine and polar species. To overwinter, insects have adapted a range of strategies to become truly cold hardy. However, although the mechanisms underlying the ability to avoid or tolerate freezing have been well studied, little attention has been given to the challenge of maintaining ion homeostasis at frigid temperatures in these species, despite this limiting cold tolerance for insects susceptible to mild chilling. Here, we investigated how prolonged exposure to temperatures just above the supercooling point affects ion balance in freeze-avoidant mountain pine beetle (Dendroctonus ponderosae) larvae in autumn, mid-winter and spring, and related it to organismal recovery times and survival. Hemolymph ion balance was gradually disrupted during the first day of exposure, characterized by hyperkalemia and hyponatremia, after which a plateau was reached and maintained for the rest of the 7-day experiment. The degree of ionoregulatory collapse correlated strongly with recovery times, which followed a similar asymptotical progression. Mortality increased slightly during extensive cold exposures, where hemolymph K+ concentration was highest, and a sigmoidal relationship was found between survival and hyperkalemia. Thus, the cold tolerance of the freeze-avoiding larvae of D. ponderosae appears limited by the ability to prevent ionoregulatory collapse in a manner similar to that of chill-susceptible insects, albeit at much lower temperatures. Based on these results, we propose that a prerequisite for the evolution of insect freeze avoidance may be a convergent or ancestral ability to maintain ion homeostasis during extreme cold stress.

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.

Stage Transitions in Lucilia sericata and Phormia regina (Diptera: Calliphoridae) and Implications for Forensic Science.

Blow fly development rates have become a key factor in estimating the postmortem interval where blow flies are among the first decomposers to occur on a body. Because the use of blow fly development requires short time durations and high accuracy, stage transition distributions are essential for proper development modeling. However, detailed examinations of stage transitions are not available for any blow fly species. Consequently, we examined this issue in two blow fly species: Lucilia sericata and Phormia regina. Transitions for all life stages across all measured temperatures were normally distributed. Use of probit analysis allowed determination of 50% transition points and associated measures of variation (i.e., standard errors). The greatest variation was noted for the L2-L3, L3-L3m, and L3m-P stage transitions. These results invalidate the notion that largest maggots should be preferentially collected for determining current maggot population stage, and further call into question the relationship between intrinsic variation and potential geographic variation in development rates.

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.

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.

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.

No time to delay reperfusion: A cross-sectional study of primary percutaneous coronary intervention times.

AIMS AND OBJECTIVES: To establish and report cross-sectional data of reperfusion times for emergency primary percutaneous coronary interventions (PPCI) and to examine factors associated with times to reperfusion. BACKGROUND: Rapid coronary reperfusion can salvage myocardial tissue, preserve left ventricular function and reduce mortality. PPCI is the gold standard of management. Researchers have reported on international median reperfusion times, but this is the first Irish study to do so. METHODS: This observational, prospective, cross-sectional study included patients diagnosed with ST-segment elevation myocardial infarction (STEMI) and admitted for emergency PPCI. Descriptive and inferential statistics were used. The study was ethically approved. We adopted the STROBE guidelines. RESULTS: All patients (N = 133) who met the inclusion criteria were included initially. Of these, 105 (79%) were diagnosed with STEMI and received emergency PPCI. The majority of STEMIs were diagnosed by paramedics and most (67%) were reperfused within 120 min, with a median time of 96 min. The results suggested that younger patients achieved timelier PPCI and source of referral was also significant in that more of those transferred directly to the coronary catheterisation laboratory achieved reperfusion within 120 min, compared with those who presented to the emergency department. CONCLUSION: A timely reperfusion service is achieved for the majority. Attention is needed in respect of the ageing and those admitted directly to the emergency departments with STEMI. RELEVANCE TO CLINICAL PRACTICE: Further international research is recommended to compare current reperfusion times against guidelines and to identify areas for improvement. Clinicians should be mindful of the importance of rapid reperfusion and the implications of its delay for patients with STEMI. Those presenting to emergency departments with chest pain should be prioritised.

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.