Identifying refugia from climate change using coupled ecological and genetic data in a transitional Mediterranean‐temperate tree species

Populations occurring in areas of overlap between the current and future distribution of a species are particularly important because they can represent “refugia from climate change”. We coupled ecological and range‐wide genetic variation data to detect such areas and to evaluate the impacts of habitat suitability changes on the genetic diversity of the transitional Mediterranean‐temperate tree Fraxinus angustifolia. We sampled and genotyped 38 natural populations comprising 1006 individuals from across Europe. We found the highest genetic diversity in western and northern Mediterranean populations, as well as a significant west to east decline in genetic diversity. Areas of potential refugia that correspond to approximately 70% of the suitable habitat may support the persistence of more than 90% of the total number of alleles in the future. Moreover, based on correlations between Bayesian genetic assignment and climate, climate change may favour the westward spread of the Black Sea gene pool in the long term. Overall, our results suggest that the northerly core areas of the current distribution contain the most important part of the genetic variation for this species and may serve as in situ macrorefugia from ongoing climate change. However, rear‐edge populations of the southern Mediterranean may be exposed to a potential loss of unique genetic diversity owing to habitat suitability changes unless populations can persist in microrefugia that have facilitated such persistence in the past.     

Metagenomic insights into microbial metabolism affecting arsenic dispersion in Mediterranean marine sediments

Microorganisms dwelling in sediments have a crucial role in biogeochemical cycles and are expected to have a strong influence on the cycle of arsenic, a metalloid responsible for severe water pollution and presenting major health risks for human populations. We present here a metagenomic study of the sediment from two harbours on the Mediterranean French coast, l'Estaque and St Mandrier. The first site is highly polluted with arsenic and heavy metals, while the arsenic concentration in the second site is below toxicity levels. The goal of this study was to elucidate the potential impact of the microbial community on the chemical parameters observed in complementary geochemical studies performed on the same sites. The metagenomic sequences, along with those from four publicly available metagenomes used as control data sets, were analysed with the RAMMCAP workflow. The resulting functional profiles were compared to determine the over‐represented Gene Ontology categories in the metagenomes of interest. Categories related to arsenic resistance and dissimilatory sulphate reduction were over‐represented in l'Estaque. More importantly, despite very similar profiles, the identification of specific sequence markers for sulphate‐reducing bacteria and sulphur‐oxidizing bacteria showed that sulphate reduction was significantly more associated with l'Estaque than with St Mandrier. We propose that biotic sulphate reduction, arsenate reduction and fermentation may together explain the higher mobility of arsenic observed in l'Estaque in previous physico‐chemical studies of this site. This study also demonstrates that it is possible to draw sound conclusions from comparing complex and similar unassembled metagenomes at the functional level, even with very low sequence coverage.     

Partitioning of genetic variation across the genome using multimarker methods in a wild bird population

The underlying basis of genetic variation in quantitative traits, in terms of the number of causal variants and the size of their effects, is largely unknown in natural populations. The expectation is that complex quantitative trait variation is attributable to many, possibly interacting, causal variants, whose effects may depend upon the sex, age and the environment in which they are expressed. A recently developed methodology in animal breeding derives a value of relatedness among individuals from high‐density genomic marker data, to estimate additive genetic variance within livestock populations. Here, we adapt and test the effectiveness of these methods to partition genetic variation for complex traits across genomic regions within ecological study populations where individuals have varying degrees of relatedness. We then apply this approach for the first time to a natural population and demonstrate that genetic variation in wing length in the great tit (Parus major) reflects contributions from multiple genomic regions. We show that a polygenic additive mode of gene action best describes the patterns observed, and we find no evidence of dosage compensation for the sex chromosome. Our results suggest that most of the genomic regions that influence wing length have the same effects in both sexes. We found a limited amount of genetic variance in males that is attributed to regions that have no effects in females, which could facilitate the sexual dimorphism observed for this trait. Although this exploratory work focuses on one complex trait, the methodology is generally applicable to any trait for any laboratory or wild population, paving the way for investigating sex‐, age‐ and environment‐specific genetic effects and thus the underlying genetic architecture of phenotype in biological study systems.     

Nectar bacteria,but not yeast,weaken a plant–pollinator mutualism

Mutualistic interactions are often subject to exploitation by species that are not directly involved in the mutualism. Understanding which organisms act as such ‘third-party’ species and how they do so is a major challenge in the current study of mutualistic interactions. Here, we show that even species that appear ecologically similar can have contrasting effects as third-party species. We experimentally compared the effects of nectar-inhabiting bacteria and yeasts on the strength of a mutualism between a hummingbird-pollinated shrub, Mimulus aurantiacus, and its pollinators. We found that the common bacterium Gluconobacter sp., but not the common yeast Metschnikowia reukaufii, reduced pollination success, seed set and nectar consumption by pollinators, thereby weakening the plant–pollinator mutualism. We also found that the bacteria reduced nectar pH and total sugar concentration more greatly than the yeasts did and that the bacteria decreased glucose concentration and increased fructose concentration whereas the yeasts affected neither. These distinct changes to nectar chemistry may underlie the microbes'' contrasting effects on the mutualism. Our results suggest that it is necessary to understand the determinants of microbial species composition in nectar and their differential modification of floral rewards to explain the mutual benefits that plants and pollinators gain from each other.     

The genetic rescue of two bottlenecked South Island robin populations using translocations of inbred donors

Populations forced through bottlenecks typically lose genetic variation and exhibit inbreeding depression. ‘Genetic rescue’ techniques that introduce individuals from outbred populations can be highly effective in reversing the deleterious effects of inbreeding, but have limited application for the majority of endangered species, which survive only in a few bottlenecked populations. We tested the effectiveness of using highly inbred populations as donors to rescue two isolated and bottlenecked populations of the South Island robin (Petroica australis). Reciprocal translocations significantly increased heterozygosity and allelic diversity. Increased genetic diversity was accompanied by increased juvenile survival and recruitment, sperm quality, and immunocompetence of hybrid individuals (crosses between the two populations) compared with inbred control individuals (crosses within each population). Our results confirm that the implementation of ‘genetic rescue’ using bottlenecked populations as donors provides a way of preserving endangered species and restoring their viability when outbred donor populations no longer exist.     

Tracking climate change in a dispersal-limited species: reduced spatial and genetic connectivity in a montane salamander

Tropical montane taxa are often locally adapted to very specific climatic conditions, contributing to their lower dispersal potential across complex landscapes. Climate and landscape features in montane regions affect population genetic structure in predictable ways, yet few empirical studies quantify the effects of both factors in shaping genetic structure of montane-adapted taxa. Here, we considered temporal and spatial variability in climate to explain contemporary genetic differentiation between populations of the montane salamander, Pseudoeurycea leprosa. Specifically, we used ecological niche modelling (ENM) and measured spatial connectivity and gene flow (using both mtDNA and microsatellite markers) across extant populations of P. leprosa in the Trans-Mexican Volcanic Belt (TVB). Our results indicate significant spatial and genetic isolation among populations, but we cannot distinguish between isolation by distance over time or current landscape barriers as mechanisms shaping population genetic divergences. Combining ecological niche modelling, spatial connectivity analyses, and historical and contemporary genetic signatures from different classes of genetic markers allows for inference of historical evolutionary processes and predictions of the impacts future climate change will have on the genetic diversity of montane taxa with low dispersal rates. Pseudoeurycea leprosa is one montane species among many endemic to this region and thus is a case study for the continued persistence of spatially and genetically isolated populations in the highly biodiverse TVB of central Mexico.     

A comparison of pedigree‐ and DNA‐based measures for identifying inbreeding depression in the critically endangered Attwater's Prairie‐chicken

The primary goal of captive breeding programmes for endangered species is to prevent extinction, a component of which includes the preservation of genetic diversity and avoidance of inbreeding. This is typically accomplished by minimizing mean kinship in the population, thereby maintaining equal representation of the genetic founders used to initiate the captive population. If errors in the pedigree do exist, such an approach becomes less effective for minimizing inbreeding depression. In this study, both pedigree‐ and DNA‐based methods were used to assess whether inbreeding depression existed in the captive population of the critically endangered Attwater's Prairie‐chicken (Tympanuchus cupido attwateri), a subspecies of prairie grouse that has experienced a significant decline in abundance and concurrent reduction in neutral genetic diversity. When examining the captive population for signs of inbreeding, variation in pedigree‐based inbreeding coefficients (f_pedigree) was less than that obtained from DNA‐based methods (f_DNA). Mortality of chicks and adults in captivity were also positively correlated with parental relatedness (r_DNA) and f_DNA, respectively, while no correlation was observed with pedigree‐based measures when controlling for additional variables such as age, breeding facility, gender and captive/release status. Further, individual homozygosity by loci (HL) and parental r_DNA values were positively correlated with adult mortality in captivity and the occurrence of a lethal congenital defect in chicks, respectively, suggesting that inbreeding may be a contributing factor increasing the frequency of this condition among Attwater's Prairie‐chickens. This study highlights the importance of using DNA‐based methods to better inform management decisions when pedigrees are incomplete or errors may exist due to uncertainty in pairings.     

Population genetics of the American eel (Anguilla rostrata): FST = 0 and North Atlantic Oscillation effects on demographic fluctuations of a panmictic species

We performed population genetic analyses on the American eel (Anguilla rostrata) with three main objectives. First, we conducted the most comprehensive analysis of neutral genetic population structure to date to revisit the null hypothesis of panmixia in this species. Second, we used this data to provide the first estimates of contemporary effective population size (N_e) and to document temporal variation in effective number of breeders (N_b) in American eel. Third, we tested for statistical associations between temporal variation in the North Atlantic Oscillation (NAO), the effective number of breeders and two indices of recruit abundance. A total of 2142 eels from 32 sampling locations were genotyped with 18 microsatellite loci. All measures of differentiation were essentially zero, and no evidence for significant spatial or temporal genetic differentiation was found. The panmixia hypothesis should thus be accepted for this species. N_b estimates varied by a factor of 23 among 12 cohorts, from 473 to 10 999. The effective population size N_e was estimated at 10 532 (95% CI, 9312–11 752). This study also showed that genetically based demographic indices, namely N_b and allelic richness (Ar), can be used as surrogates for the abundance of breeders and recruits, which were both shown to be positively influenced by variation during high (positive) NAO phases. Thus, long‐term genetic monitoring of American glass eels at several sites along the North American Atlantic coast would represent a powerful and efficient complement to census monitoring to track demographic fluctuations and better understand their causes.     

Population genomics from pool sequencing

Next generation sequencing of pooled samples is an effective approach for studies of variability and differentiation in populations. In this paper we provide a comprehensive set of estimators of the most common statistics in population genetics based on the frequency spectrum, namely the Watterson estimator , nucleotide pairwise diversity Π, Tajima's D, Fu and Li's D and F, Fay and Wu's H, McDonald‐Kreitman and HKA tests and , corrected for sequencing errors and ascertainment bias. In a simulation study, we show that pool and individual θ estimates are highly correlated and discuss how the performance of the statistics vary with read depth and sample size in different evolutionary scenarios. As an application, we reanalyse sequences from Drosophila mauritiana and from an evolution experiment in Drosophila melanogaster. These methods are useful for population genetic projects with limited budget, study of communities of individuals that are hard to isolate, or autopolyploid species.