Genomic consequences of intensive inbreeding in an isolated wolf population.

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 ROH). 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 P) was strongly correlated with realized inbreeding measured with the entire genome (r 2 = 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 ROH estimated with as few as 10,000 single nucleotide polymorphisms, than with F P. These results document in fine detail the genomic consequences of intensive inbreeding in a population of conservation concern.

A Linkage Map and QTL Analysis for Pyrethroid Resistance in the Bed Bug Cimex lectularius.

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.

Spatial and temporal genetic structure at the fourth trophic level in a fragmented landscape.

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.

Oxygen and energy availability interact to determine flight performance in the Glanville fritillary butterfly.

Flying insects have the highest known mass-specific demand for oxygen, which makes it likely that reduced availability of oxygen might limit sustained flight, either instead of or in addition to the limitation due to metabolite resources. The Glanville fritillary butterfly (Melitaea cinxia) occurs as a large metapopulation in which adult butterflies frequently disperse between small local populations. Here, we examine how the interaction between oxygen availability and fuel use affects flight performance in the Glanville fritillary. Individuals were flown under either normoxic (21 kPa O2) or hypoxic (10 kPa O2) conditions and their flight metabolism was measured. To determine resource use, levels of circulating glucose, trehalose and whole-body triglyceride were recorded after flight. Flight performance was significantly reduced in hypoxic conditions. When flown under normoxic conditions, we observed a positive correlation among individuals between post-flight circulating trehalose levels and flight metabolic rate, suggesting that low levels of circulating trehalose constrains flight metabolism. To test this hypothesis experimentally, we measured the flight metabolic rate of individuals injected with a trehalase inhibitor. In support of the hypothesis, experimental butterflies showed significantly reduced flight metabolic rate, but not resting metabolic rate, in comparison to control individuals. By contrast, under hypoxia there was no relationship between trehalose and flight metabolic rate. Additionally, in this case, flight metabolic rate was reduced in spite of circulating trehalose levels that were high enough to support high flight metabolic rate under normoxic conditions. These results demonstrate a significant interaction between oxygen and energy availability for the control of flight performance.

Predictable allele frequency changes due to habitat fragmentation in the Glanville fritillary butterfly.

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.

Flight-induced changes in gene expression in the Glanville fritillary butterfly.

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.

Outbreeding effects in an inbreeding insect, Cimex lectularius.

In some species, populations with few founding individuals can be resilient to extreme inbreeding. Inbreeding seems to be the norm in the common bed bug, Cimex lectularius, a flightless insect that, nevertheless, can reach large deme sizes and persist successfully. However, bed bugs can also be dispersed passively by humans, exposing inbred populations to gene flow from genetically distant populations. The introduction of genetic variation through this outbreeding could lead to increased fitness (heterosis) or be costly by causing a loss of local adaptation or exposing genetic incompatibility between populations (outbreeding depression). Here, we addressed how inbreeding within demes and outbreeding between distant populations impact fitness over two generations in this re-emerging public health pest. We compared fitness traits of families that were inbred (mimicking reproduction following a founder event) or outbred (mimicking reproduction following a gene flow event). We found that outbreeding led to increased starvation resistance compared to inbred families, but this benefit was lost after two generations of outbreeding. No other fitness benefits of outbreeding were observed in either generation, including no differences in fecundity between the two treatments. Resilience to inbreeding is likely to result from the history of small founder events in the bed bug. Outbreeding benefits may only be detectable under stress and when heterozygosity is maximized without disruption of coadaptation. We discuss the consequences of these results both in terms of inbreeding and outbreeding in populations with genetic and spatial structuring, as well as for the recent resurgence of bed bug populations.

Human-facilitated metapopulation dynamics in an emerging pest species, Cimex lectularius.

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.