Thrips Species Composition in Ontario Greenhouse Floriculture: Innovative Identification Tools and Implications for Integrated Pest Management.

Proper species identification is the keystone of successful integrated pest management (IPM). However, efforts to identify thrips species in Canadian greenhouses have not been formally made since the 1980s. In response to recent increases in crop damage, we sampled thrips communities from eight commercial floriculture greenhouses in the Niagara region (Ontario, Canada) from May until August 2016. Selected sites were revisited in 2017, 2018, and 2019 to determine changes in species composition over time. Western flower thrips (Frankliniella occidentalis (Pergande)), along with onion thrips (Thrips tabaci Lindeman), constituted the majority of species found. Other pest species (less than 8% of specimens across all sampling years) included poinsettia thrips (Echinothrips americanus Morgan), chrysanthemum thrips (Thrips nigropilosus Uzel), and Frankliniella fusca (Hinds). Further investigations of thrips outbreaks in Ontario from 2016 to 2023 revealed other important species, including Thrips parvispinus (Karny), Hercinothrips femoralis (Reuter), and Scirtothrips dorsalis Hood. The current biocontrol strategies used in Ontario floriculture crops for western flower thrips do not adequately control onion thrips or other thrips pests in ornamental crops, making identification a fundamental step in determining whether biocontrol or chemical control strategies should be implemented. However, traditional taxonomic keys are inaccessible to non-specialists due to their technical difficulty. Using the data gathered in these surveys, we developed a simplified, illustrated identification key for use by growers and IPM consultants.

Whole-genome sequence of the Cooley spruce gall adelgid, Adelges cooleyi (Hemiptera: Sternorrhyncha: Adelgidae).

The adelgids (Adelgidae) are a small family of sap-feeding insects, which, together with true aphids (Aphididae) and phylloxerans (Phylloxeridae), make up the infraorder Aphidomorpha. Some adelgid species are highly destructive to forest ecosystems such as Adelges tsugae, Adelges piceae, Adelges laricis, Pineus pini, and Pineus boerneri. Despite this, there are no high-quality genomic resources for adelgids, hindering advanced genomic analyses within Adelgidae and among Aphidomorpha. Here, we used PacBio continuous long-read and Illumina RNA-sequencing to construct a high-quality draft genome assembly for the Cooley spruce gall adelgid, Adelges cooleyi (Gillette), a gall-forming species endemic to North America. The assembled genome is 270.2 Mb in total size and has scaffold and contig N50 statistics of 14.87 and 7.18 Mb, respectively. There are 24,967 predicted coding sequences, and the assembly completeness is estimated at 98.1 and 99.6% with core BUSCO gene sets of Arthropoda and Hemiptera, respectively. Phylogenomic analysis using the A. cooleyi genome, 3 publicly available adelgid transcriptomes, 4 phylloxera transcriptomes, the Daktulosphaira vitifoliae (grape phylloxera) genome, 4 aphid genomes, and 2 outgroup coccoid genomes fully resolves adelgids and phylloxerans as sister taxa. The mitochondrial genome is 24 kb, among the largest in insects sampled to date, with 39.4% composed of noncoding regions. This genome assembly is currently the only genome-scale, annotated assembly for adelgids and will be a valuable resource for understanding the ecology and evolution of Aphidomorpha.

Discovery of the primary aphid (Hemiptera: Aphidomorpha) and scale insect (Hemiptera: Coccomorpha) type specimens from the collection of Theodor Hartig (18051880).

Theodor Hartigs aphid and scale insect type specimens have been presumed lost or destroyed for the last 140 years. Here we document their discovery at the Bavarian State Collection of Zoology (Zoologische Staatssammlung Mnchen, ZSM), in Munich, Germany. These specimens include primary types for 24 aphid, three adelgid, and two armored scale insect species named by Hartig between 1834 and 1851, as well as other specimens of unknown importance.

Geographic isolation drives speciation in Nearctic aphids.

Across herbivorous insect clades, species richness and host-use diversity tend to positively covary. This could be because host-use divergence drives speciation, or because it raises the ecological limits on species richness. To evaluate these hypotheses, we performed phylogenetic path model analyses of the species diversity of Nearctic aphids. Here, we show that variation in the species richness of aphid clades is caused mainly by host-use divergence, whereas variation in speciation rates is caused more by divergence in non-host-related niche variables. Aphid speciation is affected by both the evolution of host and non-host-related niche components, but the former is largely caused by the latter. Thus, our analyses suggest that host-use divergence can both raise the ecological limits on species richness and drive speciation, although in the latter case, host-use divergence tends to be a step along the causal path leading from non-host-related niche evolution to speciation.

Continent-wide population genomic structure and phylogeography of North America's most destructive conifer defoliator, the spruce budworm (Choristoneura fumiferana).

The spruce budworm, Choristoneura fumiferana, is presumed to be panmictic across vast regions of North America. We examined the extent of panmixia by genotyping 3,650 single nucleotide polymorphism (SNP) loci in 1975 individuals from 128 collections across the continent. We found three spatially structured subpopulations: Western (Alaska, Yukon), Central (southeastern Yukon to the Manitoba-Ontario border), and Eastern (Manitoba-Ontario border to the Atlantic). Additionally, the most diagnostic genetic differentiation between the Central and Eastern subpopulations was chromosomally restricted to a single block of SNPs that may constitute an island of differentiation within the species. Geographic differentiation in the spruce budworm parallels that of its principal larval host, white spruce (Picea glauca), providing evidence that spruce budworm and spruce trees survived in the Beringian refugium through the Last Glacial Maximum and that at least two isolated spruce budworm populations diverged with spruce/fir south of the ice sheets. Gene flow in the spruce budworm may also be affected by mountains in western North America, habitat isolation in West Virginia, regional adaptations, factors related to dispersal, and proximity of other species in the spruce budworm species complex. The central and eastern geographic regions contain individuals that assign to Eastern and Central subpopulations, respectively, indicating that these barriers are not complete. Our discovery of previously undetected geographic and genomic structure in the spruce budworm suggests that further population modelling of this ecologically important insect should consider regional differentiation, potentially co-adapted blocks of genes, and gene flow between subpopulations.

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.

Convergent herbivory on conifers by Choristoneura moths after boreal forest formation.

Mitogenomes are useful markers for phylogenetic studies across a range of taxonomic levels. Here, we focus on mitogenome variation across the tortricid moth genus Choristoneura and particularly the spruce budworm (Choristoneura fumiferana) species complex, a notorious pest group of North American conifer forests. Phylogenetic relationships of Tortricidae, representing two subfamilies, four tribes and nine genera, were analyzed using 21 mitogenomes. These included six newly-sequenced mitogenomes for species in the spruce budworm complex plus three additional Choristoneura species and 12 previously published mitogenomes from other tortricids and one from the Cossidae. We evaluated the phylogenetic informativeness of the mitogenomes and reconstructed a time-calibrated tree with fossil and secondary calibrations. We found that tortricid mitogenomes had conserved protein and ribosomal regions, and analysis of all protein-coding plus ribosomal genes together provided an efficient marker at any taxonomic rank. The time-calibrated phylogeny showed evolutionary convergence of conifer feeding within Choristoneura, with two independent lineages, the Nearctic spruce budworm complex and the Palearctic species Choristoneura murinana, both shifting onto conifers about 11 million years ago from angiosperms. These two host-plant shifts both occurred after the formation of boreal forest in the late Miocene. Haplotype diversification within the spruce budworm complex occurred in the last 4 million years, and is probably linked to the initial cooling cycles of the Northern Hemisphere in the Pliocene.

Distinct sources of gene flow produce contrasting population genetic dynamics at different range boundaries of a Choristoneura budworm.

Populations are often exposed to multiple sources of gene flow, but accounts are lacking of the population genetic dynamics that result from these interactions or their effects on local evolution. Using a genomic clines framework applied to 1,195 single nucleotide polymorphisms, we documented genomewide, locus-specific patterns of introgression between Choristoneura occidentalis biennis spruce budworms and two ecologically divergent relatives, C. o. occidentalis and Choristoneura fumiferana, that it interacts with at alternate boundaries of its range. We observe contrasting hybrid indexes between the two hybrid zones, no overlap in "gene-flow outliers" (clines showing relatively extreme extents or rates of locus-specific introgression) and variable linkage disequilibrium among those outliers. At the same time, correlated genomewide rates of introgression between zones suggest the presence of processes common to both boundaries. These findings highlight the contrasting population genetic dynamics that can occur at separate frontiers of a single population, while also suggesting that shared patterns may frequently accompany cases of divergence-with-gene-flow that involve a lineage in common. Our results point to potentially complex evolutionary outcomes for populations experiencing multiple sources of gene flow.

Life-stage differences in spatial genetic structure in an irruptive forest insect: implications for dispersal and spatial synchrony.

Dispersal determines the flux of individuals, energy and information and is therefore a key determinant of ecological and evolutionary dynamics. Yet, it remains difficult to quantify its importance relative to other factors. This is particularly true in cyclic populations in which demography, drift and dispersal contribute to spatio-temporal variability in genetic structure. Improved understanding of how dispersal influences spatial genetic structure is needed to disentangle the multiple processes that give rise to spatial synchrony in irruptive species. In this study, we examined spatial genetic structure in an economically important irruptive forest insect, the spruce budworm (Choristoneura fumiferana) to better characterize how dispersal, demography and ecological context interact to influence spatial synchrony in a localized outbreak. We characterized spatial variation in microsatellite allele frequencies using 231 individuals and seven geographic locations. We show that (i) gene flow among populations is likely very high (Fst  ≈ 0); (ii) despite an overall low level of genetic structure, important differences exist between adult (moth) and juvenile (larvae) life stages; and (iii) the localized outbreak is the likely source of moths captured elsewhere in our study area. This study demonstrates the potential of using molecular methods to distinguish residents from migrants and for understanding how dispersal contributes to spatial synchronization. In irruptive populations, the strength of genetic structure depends on the timing of data collection (e.g. trough vs. peak), location and dispersal. Taking into account this ecological context allows us to make more general characterizations of how dispersal can affect spatial synchrony in irruptive populations.

Transverse and vortex domain wall structure in magnetic nanowires with uniaxial in-plane anisotropy.

Micromagnetic and analytical models are used to investigate how in-plane uniaxial anisotropy affects transverse and vortex domain walls in nanowires where shape anisotropy dominates. The effect of the uniaxial anisotropy can be interpreted as a modification of the effective wire dimensions. When the anisotropy axis is aligned with the wire axis (θ(a) = 0), the wall width is narrower than when no anisotropy is present. Conversely, the wall width increases when the anisotropy axis is perpendicular to the wire axis (θ(a) = π/2). The anisotropy also affects the nanowire dimensions at which transverse walls become unstable. This phase boundary shifts to larger widths or thicknesses when θ(a) = 0, but smaller widths or thicknesses when θ(a) = π/2.