Repurposing population genetics data to discern genomic architecture: A case study of linkage cohort detection in mountain pine beetle (Dendroctonus ponderosae)

Genetic management on the brink of extinction: sequencing microsatellites does not improve estimates of inbreeding in wild and captive Vancouver Island marmots (Marmota vancouverensis)

Captive breeding is often a last resort management option in the conservation of endangered species which can in turn lead to increased risk of inbreeding depression and loss of genetic diversity. Thus, recording breeding events via studbook for the purpose of estimating relatedness, and facilitating mating pair selection to minimize inbreeding, is common practice. However, as founder relatedness is often unknown, loss of genetic variation and inbreeding cannot be entirely avoided. Molecular genotyping is slowly being adopted in captive breeding programs, however achieving sufficient resolution can be challenging in small, low diversity, populations. Here, we evaluate the success of the Vancouver Island marmot (Marmota vancouverensis; VIM; among the worlds most endangered mammals) captive breeding program in preventing inbreeding and maintaining genetic diversity. We explored the use of high-throughput amplicon sequencing of microsatellite regions to assay greater genetic variation in both captive and wild populations than traditional length-based fragment analysis. Contrary to other studies, this method did not considerably increase diversity estimates, suggesting: (1) that the technique does not universally improve resolution, and (2) VIM have exceedingly low diversity. Studbook estimates of pairwise relatedness and inbreeding in the current population were weakly, but positively, correlated to molecular estimates. Thus, current studbooks are moderately effective at predicting genetic similarity when founder relatedness is known. Finally, we found that captive and wild populations did not differ in allelic frequencies, and conservation efforts to maintain diversity have been successful with no significant decrease in diversity over the last three generations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10592-022-01429-7.

Chromosome-level genome assembly reveals genomic architecture of northern range expansion in the mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae)

Genome sequencing methods and assembly tools have improved dramatically since the 2013 publication of draft genome assemblies for the mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae). We conducted proximity ligation library sequencing and scaffolding to improve contiguity, and then used linkage mapping and recent bioinformatic tools for correction and further improvement. The new assemblies have dramatically improved contiguity and gaps compared to the originals: N50 values increased 26- to 36-fold, and the number of gaps were reduced by half. Ninety per cent of the content of the assemblies is now contained in 12 and 11 scaffolds for the female and male assemblies, respectively. Based on linkage mapping information, the 12 largest scaffolds in both assemblies represent all 11 autosomal chromosomes and the neo-X chromosome. These assemblies now have nearly chromosome-sized scaffolds and will be instrumental for studying genomic architecture, chromosome evolution, population genomics, functional genomics, and adaptation in this and other pest insects. We also identified regions in two chromosomes, including the ancestral-X portion of the neo-X chromosome, with elevated differentiation between northern and southern Canadian populations.

Rapid Adaptation to Temperature via a Potential Genomic Island of Divergence in the Invasive Green Crab, Carcinus maenas

Widespread species often adapt easily to novel conditions – both those found in new habitats and those generated by climate change. However, rapid adaptation may be hindered in the marine realm, where long-distance dispersal and consequently high gene flow are predicted to limit potential for local adaptation. Here, we use a highly dispersive invasive marine crab to test the nature and speed of adaptation to temperature in the sea. Using single nucleotide polymorphisms (SNPs) generated from cardiac transcriptome sequencing, we characterized six populations of the European green crab (Carcinus maenas) located across parallel thermal gradients in their native and invasive ranges. We compared SNP frequencies with local temperatures and previously generated data on cardiac heat and cold tolerance to identify candidate markers associated with population-level differences in thermal physiology. Of 10,790 SNPs, 104 were identified as frequency outliers, a signal that was strongly driven by association with temperature and/or cold tolerance. Seventy-two of these outlier markers, representing 28 different genes, were in a cluster of SNPs identified as a potential inversion polymorphism using linkage disequilibrium network analysis. This SNP cluster was unique in the data set, which was otherwise characterized by low levels of linkage disequilibrium, and markers in this cluster showed a significant enrichment of coding substitutions relative to the full SNP set. These 72 outlier SNPs appear to be transmitted as a unit, and represent a putative genomic island of divergence which varied in frequency with organismal cold tolerance. This relationship was strikingly similar across both native and invasive populations, all of which showed a very strong correlation with cold tolerance (R2 = 0.96 over all six populations). Notably, three of these populations have diverged recently (<100 years) and show little to no neutral divergence, suggesting that this genomic region may be responding to temperature on a relatively short time scale. This relationship indicates adaptation to temperature based on the action of a putative genomic island of divergence, perhaps partially explaining the extraordinary invasive ability of this species.

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.