Inbreeding (mating between relatives) is a major concern for conservation as it decreases individual fitness and can increase the risk of population extinction. We used whole-genome resequencing of 97 grey wolves (Canis lupus) from the highly inbred Scandinavian wolf population to identify 'identical-by-descent' (IBD) chromosome segments as runs of homozygosity (ROH). This gave the high resolution required to precisely measure realized inbreeding as the IBD fraction of the genome in ROH (F ). We found a striking pattern of complete or near-complete homozygosity of entire chromosomes in many individuals. The majority of individual inbreeding was due to long IBD segments (>5 cM) originating from ancestors ≤10 generations ago, with 10 genomic regions showing very few ROH and forming candidate regions for containing loci contributing strongly to inbreeding depression. Inbreeding estimated with an extensive pedigree (F ) was strongly correlated with realized inbreeding measured with the entire genome (r = 0.86). However, inbreeding measured with the whole genome was more strongly correlated with multi-locus heterozygosity estimated with as few as 500 single nucleotide polymorphisms, and with F estimated with as few as 10,000 single nucleotide polymorphisms, than with F. These results document in fine detail the genomic consequences of intensive inbreeding in a population of conservation concern.

Describing the evolutionary dynamics of now extinct populations is challenging, as their genetic composition before extinction is generally unknown. The Glanville fritillary butterfly has a large extant metapopulation in the Åland Islands in Finland, but declined to extinction in the nearby fragmented southwestern (SW) Finnish archipelago in the 20th century. We genotyped museum samples for 222 SNPs across the genome, including SNPs from candidate genes and neutral regions. SW Finnish populations had significantly reduced genetic diversity before extinction, and their allele frequencies gradually diverged from those in contemporary Åland populations over 80 y. We identified 15 outlier loci among candidate SNPs, mostly related to flight, in which allele frequencies have changed more than the neutral expectation. At outlier loci, allele frequencies in SW Finland shifted in the same direction as newly established populations deviated from old local populations in contemporary Åland. Moreover, outlier allele frequencies in SW Finland resemble those in fragmented landscapes as opposed to continuous landscapes in the Baltic region. These results indicate selection for genotypes associated with good colonization capacity in the highly fragmented landscape before the extinction of the populations. Evolutionary response to habitat fragmentation may have enhanced the viability of the populations, but it did not save the species from regional extinction in the face of severe habitat loss and fragmentation. These results highlight a potentially common situation in changing environments: evolutionary changes are not strong enough to fully compensate for the direct adverse effects of environmental change and thereby rescue populations from extinction.

The number and demographic history of colonists can have dramatic consequences for the way in which genetic diversity is distributed and maintained in a metapopulation. The bed bug (Cimex lectularius) is a re-emerging pest species whose close association with humans has led to frequent local extinction and colonization, that is, to metapopulation dynamics. Pest control limits the lifespan of subpopulations, causing frequent local extinctions, and human-facilitated dispersal allows the colonization of empty patches. Founder events often result in drastic reductions in diversity and an increased influence of genetic drift. Coupled with restricted migration, this can lead to rapid population differentiation. We therefore predicted strong population structuring. Here, using 21 newly characterized microsatellite markers and approximate Bayesian computation (ABC), we investigate simplified versions of two classical models of metapopulation dynamics, in a coalescent framework, to estimate the number and genetic composition of founders in the common bed bug. We found very limited diversity within infestations but high degrees of structuring across the city of London, with extreme levels of genetic differentiation between infestations (FST = 0.59). ABC results suggest a common origin of all founders of a given subpopulation and that the numbers of colonists were low, implying that even a single mated female is enough to found a new infestation successfully. These patterns of colonization are close to the predictions of the propagule pool model, where all founders originate from the same parental infestation. These results show that aspects of metapopulation dynamics can be captured in simple models and provide insights that are valuable for the future targeted control of bed bug infestations.

Dispersal is important for determining both species ecological processes, such as population viability, and its evolutionary processes, like gene flow and local adaptation. Yet obtaining accurate estimates in the wild through direct observation can be challenging or even impossible, particularly over large spatial and temporal scales. Genotyping many individuals from wild populations can provide detailed inferences about dispersal. We therefore utilized genomewide marker data to estimate dispersal in the classic metapopulation of the Glanville fritillary butterfly (Melitaea cinxia L.), in the Åland Islands in SW Finland. This is an ideal system to test the effectiveness of this approach due to the wealth of information already available covering dispersal across small spatial and temporal scales, but lack of information at larger spatial and temporal scales. We sampled three larvae per larval family group from 3732 groups over a six‐year period and genotyped for 272 SNPs across the genome. We used this empirical data set to reconstruct cases where full‐sibs were detected in different local populations to infer female effective dispersal distance, that is, dispersal events directly contributing to gene flow. On average this was one kilometre, closely matching previous dispersal estimates made using direct observation. To evaluate our power to detect full‐sib families, we performed forward simulations using an individual‐based model constructed and parameterized for the Glanville fritillary metapopulation. Using these simulations, 100% of predicted full‐sibs were correct and over 98% of all true full‐sib pairs were detected. We therefore demonstrate that even in a highly dynamic system with a relatively small number of markers, we can accurately reconstruct full‐sib families and for the first time make inferences on female effective dispersal. This highlights the utility of this approach in systems where it has previously been impossible to obtain accurate estimates of dispersal over both ecological and evolutionary scales.

A fragmented habitat becomes increasingly fragmented for species at higher trophic levels, such as parasitoids. To persist, these species are expected to possess life-history traits, such as high dispersal, that facilitate their ability to use resources that become scarce in fragmented landscapes. If a specialized parasitoid disperses widely to take advantage of a sparse host, then the parasitoid population should have lower genetic structure than the host. We investigated the temporal and spatial genetic structure of a hyperparasitoid (fourth trophic level) in a fragmented landscape over 50 Â 70 km, using microsatellite markers, and compared it with the known structures of its host parasitoid, and the butterfly host which lives as a classic metapopulation. We found that population genetic structure decreases with increasing trophic level. The hyperparasitoid has fewer genetic clusters (K ¼ 4), than its host parasitoid (K ¼ 15), which in turn is less structured than the host butterfly (K ¼ 27). The genetic structure of the hyperparasitoid also shows temporal variation, with genetic differentiation increasing due to reduction of the population size, which reduces the effective population size. Overall, our study confirms the idea that specialized species must be dispersive to use a fragmented host resource, but that this adaptation has limits.

Insect flight is one of the most energetically demanding activities in the animal kingdom, yet for many insects flight is necessary for reproduction and foraging. Moreover, dispersal by flight is essential for the viability of species living in fragmented landscapes. Here, working on the Glanville fritillary butterfly (Melitaea cinxia), we use transcriptome sequencing to investigate gene expression changes caused by 15 min of flight in two contrasting populations and the two sexes. Male butterflies and individuals from a large metapopulation had significantly higher peak flight metabolic rate (FMR) than female butterflies and those from a small inbred population. In the pooled data, FMR was significantly positively correlated with genome-wide heterozygosity, a surrogate of individual inbreeding. The flight experiment changed the expression level of 1513 genes, including genes related to major energy metabolism pathways, ribosome biogenesis and RNA processing, and stress and immune responses. Males and butterflies from the population with high FMR had higher basal expression of genes related to energy metabolism, whereas females and butterflies from the small population with low FMR had higher expression of genes related to ribosome/RNA processing and immune response. Following the flight treatment, genes related to energy metabolism were generally down-regulated, while genes related to ribosome/RNA processing and immune response were up-regulated. These results suggest that common molecular mechanisms respond to flight and can influence differences in flight metabolic capacity between populations and sexes.

The rapid evolution of insecticide resistance remains one of the biggest challenges in the control of medically and economically important pests. Insects have evolved a diverse range of mechanisms to reduce the efficacy of the commonly used classes of insecticides, and finding the genetic basis of resistance is a major aid to management. In a previously unstudied population, we performed an F2 resistance mapping cross for the common bed bug, Cimex lectularius, for which insecticide resistance is increasingly widespread. Using 334 SNP markers obtained through RAD-sequencing, we constructed the first linkage map for the species, consisting of 14 putative linkage groups (LG), with a length of 407 cM and an average marker spacing of 1.3 cM. The linkage map was used to reassemble the recently published reference genome, facilitating refinement and validation of the current genome assembly. We detected a major QTL on LG12 associated with insecticide resistance, occurring in close proximity (1.2 Mb) to a carboxylesterase encoding candidate gene for pyrethroid resistance. This provides another example of this candidate gene playing a major role in determining survival in a bed bug population following pesticide resistance evolution. The recent availability of the bed bug genome, complete with a full list of potential candidate genes related to insecticide resistance, in addition to the linkage map generated here, provides an excellent resource for future research on the development and spread of insecticide resistance in this resurging pest species.

For conservation genetic studies using non-invasively collected samples, genome-wide data may be hard to acquire. Until now, such studies have instead mostly relied on analyses of traditional genetic markers such as microsatellites (SSRs). Recently, high throughput genotyping of single nucleotide polymorphisms (SNPs) has become available, expanding the use of genomic methods to include non-model species of conservation concern. We have developed a 96-marker SNP array for use in applied conservation monitoring of the Scandinavian wolverine (Gulo gulo) population. By genotyping more than a thousand non-invasively collected samples, we were able to obtain precise estimates of different types of genotyping errors and sample dropout rates. The SNP panel significantly outperforms the SSR markers (and DBY intron markers for sexing) both in terms of precision in genotyping, sex assignment and individual identification, as well as in the proportion of samples successfully genotyped. Furthermore, SNP genotyping offers a simplified laboratory and analysis pipeline with fewer samples needed to be repeatedly genotyped in order to obtain reliable consensus data. In addition, we utilised a unique opportunity to successfully demonstrate the application of SNP genotype data for reconstructing pedigrees in wild populations, by validating the method with samples from wild individuals with known relatedness. By offering a simplified workflow with improved performance, we anticipate this methodology will facilitate the use of non-invasive samples to improve genetic management of many different types of populations that have previously been challenging to survey.