Sirex noctilio is a major international quarantine pest that recently emerged in northeast China to specifically invade conifers. During female oviposition, venom is injected into the host together with its symbiotic fungus to alter the normal Pinus physiology and weaken or even kill the tree. In China, the Mongolian pine (Pinus sylvestris var. mongolica), an important wind-proof and sand-fixing species, is the unique host of S. noctilio. To explore the interplay between S. noctilio venom and Mongolian pine, we performed a transcriptome comparative analysis of a 10-year-old Mongolian pine after wounding and inoculation with S. noctilio venom. The analysis was performed at 12 h, 24 h and 72 h. PacBio ISO-seq was used and integrated with RNA-seq to construct an accurate full-length transcriptomic database. We obtained 52,963 high-precision unigenes, consisting of 48,654 (91.86%) unigenes that were BLASTed to known sequences in the public database and 4309 unigenes without any annotation information, which were presumed to be new genes. The number of differentially expressed genes (DEGs) increased with the treatment time, and the DEGs were most abundant at 72 h. A total of 706 inoculation-specific DEGs (475 upregulated and 231 downregulated) and 387 wounding-specific DEGs (183 upregulated and 204 downregulated) were identified compared with the control. Under venom stress, we identified 6 DEGs associated with reactive oxygen species (ROS) and 20 resistance genes in Mongolian pine. Overall, 52 transcription factors (TFs) were found under venom stress, 45 of which belonged to the AP2/ERF TF family and were upregulated. A total of 13 genes related to the photosystem, 3 genes related photo-regulation, and 9 TFs were identified under wounding stress. In conclusion, several novel putative genes were found in Mongolian pine by PacBio ISO seq. Meanwhile, we also identified various genes that were resistant to S. noctilio venom, such as GAPDH, GPX, CAT, FL2, CERK1, and HSP83A, etc.

Obolodiplosis robiniae was discovered in Eurasia at the beginning of the 21st century. In this study, we explore the present and future (in the years 2050 and 2070) trends in the potential distribution of O. robiniae in Eurasia under diverse climate change scenarios based on a maximum entropy model. Our findings indicated that the current potential distribution area of O. robiniae is within the range of 21°34′ and 65°39′ N in the Eurasian continent. The primary factor controlling the distribution of O. robiniae is temperature. The highly and moderately suitable areas are mainly distributed in the semi-humid and semi-arid regions, which also happen to be the locations where the host black locust (Robinia pseudoacacia L.) grows at its fastest rate. The forecast of the potential distribution area of O. robiniae revealed that the species would benefit from global warming. The region suitable for the habitat of O. robiniae is characterized by a large-scale northward expansion trend and an increase in temperature. This information would help the forestry quarantine departments of Asian and European countries provide early warnings on the probable distribution areas of O. robiniae and provide a scientific basis for the prevention and control of O. robiniae spread and outbreaks.

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.

Dutch elm disease is one of the most devastating plant diseases, primarily spread through bark beetles. Scolytus scolytus is a key vector of this disease. In this study, distribution data of S. scolytus were collected and filtered. Combined with environmental and climatic variables, an ensemble model was developed using the Biomod2 platform to predict its potential geographical distribution in China. The selection of climate variables was critical for accurate prediction. Eight bioclimatic factors with high importance were selected from 19 candidate variables. Among these, the three most important factors are the minimum temperature of the coldest month (bio6), precipitation seasonality (bio15), and precipitation in the driest quarter (bio17). Under current climate conditions, suitable habitats for S. scolytus are mainly located in the temperate regions between 30° and 60° N latitude. These include parts of Europe, East Asia, eastern and northwestern North America, and southern and northeastern South America. In China, the low-suitability area was estimated at 37,883.39 km², and the medium-suitability area at 251.14 km². No high-suitability regions were identified. However, low-suitability zones were widespread across multiple provinces. Under future climate scenarios, low-suitability areas are still projected across China. Medium-suitability areas are expected to increase under SSP370 and SSP585, particularly along the eastern coastal regions, peaking between 2041 and 2060. High-suitability zones may also emerge under these two scenarios, again concentrated in coastal areas. These findings provide a theoretical basis for entry quarantine measures and early warning systems aimed at controlling the spread of S. scolytus in China.

In the Western Palaearctic, the congeneric bark beetles Hylurgus ligniperda and Hylurgus micklitzi are developing on pine trees. Although H. ligniperda has a larger native range in Europe, both species often occur in sympatry in the Mediterranean Basin. They share rather similar life histories, thus providing an excellent model for comparative studies on species differentiation and ecological traits. Moreover, along with increasing human activities and international trade, H. ligniperda has invaded all continents, whereas H. micklitzi remains confined to the Mediterranean Basin. Our research utilised COI markers to assess the genetic diversity and structure, alongside the demographic history, of both species and revealed marked differences. H. micklitzi exhibited low genetic diversity and shallow population structures with restricted expansion. In contrast, H. ligniperda displayed a longer and more complex post-glacial evolutionary route with four distinct clades originating from separate glacial refugia, coupled with considerable genetic variability across Europe. Comparative analyses about the two species suggest that species traits, such as larger body size, increased voltinism, broader ecological niches and frequent introduction events may have helped the fast global invasion and spread of H. ligniperda over the past century.

Hylurgus ligniperda (F.) is a highly successful invader among bark beetles (Scolytinae) and forest insects in general. Native to the western Palearctic region, it has become established in every continent where its host plants ( Pinus spp.) occur. Especially in southern hemisphere regions with large pine plantations, it is often highly abundant. As a repeat invader with a wealth of information on various aspects relevant for biological invasions, it is highly suitable as a model organism for studying the role of international trade, the planting of non-native trees, and the potential occurrence of bridgehead invasions (where abundant non-native populations precipitate further invasions). In the present study, our specific objectives were to reconstruct the worldwide invasions of H. ligniperda and the pathways involved by using a multi-pronged approach including population genetics, analysis of historic interception data generated from inspections of imports, and records of establishments in the literature. Our review of the native and non-native ranges of H. ligniperda and the chronology of establishments revealed at least 13 separate invasions of non-native regions, beginning with Madeira (Portugal) before 1850, and, most recently, eastern China in 2019. We compared the population genetics of 464 specimens from eight countries in the native range and eight countries in the non-native range. Sequencing of the mitochondrial COI gene revealed the presence of 29 haplotypes in six well-supported clades, based on a Bayesian analysis. Non-native populations had significantly lower haplotype diversity (mean h = 0.219) than populations in the native range (mean h = 0.691). Countries in the non-native range had an average of about two haplotypes compared with about four haplotypes in native countries. In the non-native range, only one or two haplotypes were dominant, and these differed among invaded regions except for haplotype HL-H3 which occurred in Australia, New Zealand and China as well as in four countries in southern Europe, and HL-H4 which was dominant in New Zealand, California, New York State, and eastern China as well as two countries in southern Europe. Analyses of interceptions of H. ligniperda with imports arriving in five countries revealed that between 74% and 99% of interceptions originated from other non-native regions while in the USA, most interceptions were linked to imports from the native range beginning in the 1970s. Based on the combined evidence of the chronology of invasions, interception data, and analysis of haplotype distribution, we conclude that the early invasions (before 1950) probably all originated from the native range, while several of the more recent invasions probably originated from parts of the non-native range (suggestive of a bridgehead effect). However, it cannot be determined with certainty what the original sources of each of the invading populations were.