Projection of current and future distribution of adaptive genetic units in an alpine ungulate.

Climate projections predict major changes in alpine environments by the end of the 21st century. To avoid climate-induced maladaptation and extinction, many animal populations will either need to move to more suitable habitats or adapt in situ to novel conditions. Since populations of a species exhibit genetic variation related to local adaptation, it is important to incorporate this variation into predictive models to help assess the ability of the species to survive climate change. Here, we evaluate how the adaptive genetic variation of a mountain ungulate-the Northern chamois (Rupicapra rupicapra)-could be impacted by future global warming. Based on genotype-environment association analyses of 429 chamois using a ddRAD sequencing approach, we identified genetic variation associated with climatic gradients across the European Alps. We then delineated adaptive genetic units and projected the optimal distribution of these adaptive groups in the future. Our results suggest the presence of local adaptation to climate in Northern chamois with similar genetic adaptive responses in geographically distant but climatically similar populations. Furthermore, our results predict that future climatic changes will modify the Northern chamois adaptive landscape considerably, with various degrees of maladaptation risk.

Inter-Specific Genetic Exchange Despite Strong Divergence in Deep-Sea Hydrothermal Vent Gastropods of the Genus Alviniconcha.

Deep hydrothermal vents are highly fragmented and unstable habitats at all temporal and spatial scales. Such environmental dynamics likely play a non-negligible role in speciation. Little is, however, known about the evolutionary processes that drive population-level differentiation and vent species isolation and, more specifically, how geography and habitat specialisation interplay in the species history of divergence. In this study, the species range and divergence of Alviniconcha snails that occupy active Western Pacific vent fields was assessed by using sequence variation data of the mitochondrial Cox1 gene, RNAseq, and ddRAD-seq. Combining morphological description and sequence datasets of the three species across five basins, we confirmed that A. kojimai, A. boucheti, and A. strummeri, while partially overlapping over their range, display high levels of divergence in the three genomic compartments analysed that usually encompass values retrieved for reproductively isolated species with divergences rang from 9% to 12.5% (mtDNA) and from 2% to 3.1% (nuDNA). Moreover, the three species can be distinguished on the basis of their external morphology by observing the distribution of bristles and the shape of the columella. According to this sampling, A. boucheti and A. kojimai form an east-to-west species abundance gradient, whereas A. strummeri is restricted to the Futuna Arc/Lau and North Fiji Basins. Surprisingly, population models with both gene flow and population size heterogeneities among genomes indicated that these three species are still able to exchange genes due to secondary contacts at some localities after a long period of isolation.

Searching for genetic evidence of demographic decline in an arctic seabird: beware of overlapping generations.

Genetic data are useful for detecting sudden population declines in species that are difficult to study in the field. Yet this indirect approach has its own drawbacks, including population structure, mutation patterns, and generation overlap. The ivory gull (Pagophila eburnea), a long-lived Arctic seabird, is currently suffering from rapid alteration of its primary habitat (i.e., sea ice), and dramatic climatic events affecting reproduction and recruitment. However, ivory gulls live in remote areas, and it is difficult to assess the population trend of the species across its distribution. Here we present complementary microsatellite- and SNP-based genetic analyses to test a recent bottleneck genetic signal in ivory gulls over a large portion of their distribution. With attention to the potential effects of population structure, mutation patterns, and sample size, we found no significant signatures of population decline worldwide. At a finer scale, we found a significant bottleneck signal at one location in Canada. These results were compared with predictions from simulations showing how generation time and generation overlap can delay and reduce the bottleneck microsatellite heterozygosity excess signal. The consistency of the results obtained with independent methods strongly indicates that the species shows no genetic evidence of an overall decline in population size. However, drawing conclusions related to the species' population trends will require a better understanding of the effect of age structure in long-lived species. In addition, estimates of the effective global population size of ivory gulls were surprisingly low (~1000 ind.), suggesting that the evolutionary potential of the species is not assured.

Subtle limits to connectivity revealed by outlier loci within two divergent metapopulations of the deep-sea hydrothermal gastropod Ifremeria nautilei.

Hydrothermal vents form archipelagos of ephemeral deep-sea habitats that raise interesting questions about the evolution and dynamics of the associated endemic fauna, constantly subject to extinction-recolonization processes. These metal-rich environments are coveted for the mineral resources they harbour, thus raising recent conservation concerns. The evolutionary fate and demographic resilience of hydrothermal species strongly depend on the degree of connectivity among and within their fragmented metapopulations. In the deep sea, however, assessing connectivity is difficult and usually requires indirect genetic approaches. Improved detection of fine-scale genetic connectivity is now possible based on genome-wide screening for genetic differentiation. Here, we explored population connectivity in the hydrothermal vent snail Ifremeria nautilei across its species range encompassing five distinct back-arc basins in the Southwest Pacific. The global analysis, based on 10,570 single nucleotide polymorphism (SNP) markers derived from double digest restriction-site associated DNA sequencing (ddRAD-seq), depicted two semi-isolated and homogeneous genetic clusters. Demogenetic modeling suggests that these two groups began to diverge about 70,000 generations ago, but continue to exhibit weak and slightly asymmetrical gene flow. Furthermore, a careful analysis of outlier loci showed subtle limitations to connectivity between neighbouring basins within both groups. This finding indicates that migration is not strong enough to totally counterbalance drift or local selection, hence questioning the potential for demographic resilience at this latter geographical scale. These results illustrate the potential of large genomic data sets to understand fine-scale connectivity patterns in hydrothermal vents and the deep sea.

Age-specific survivorship and fecundity shape genetic diversity in marine fishes.

Genetic diversity varies among species due to a range of eco-evolutionary processes that are not fully understood. The neutral theory predicts that the amount of variation in the genome sequence between different individuals of the same species should increase with its effective population size ( N e ). In real populations, multiple factors that modulate the variance in reproductive success among individuals cause N e to differ from the total number of individuals ( N ). Among these, age-specific mortality and fecundity rates are known to have a direct impact on the N e / N ratio. However, the extent to which vital rates account for differences in genetic diversity among species remains unknown. Here, we addressed this question by comparing genome-wide genetic diversity across 16 marine fish species with similar geographic distributions but contrasted lifespan and age-specific survivorship and fecundity curves. We sequenced the whole genome of 300 individuals to high coverage and assessed their genome-wide heterozygosity with a reference-free approach. Genetic diversity varied from 0.2% to 1.4% among species, and showed a negative correlation with adult lifespan, with a large negative effect ( s l o p e = - 0.089 per additional year of lifespan) that was further increased when brooding species providing intense parental care were removed from the dataset ( s l o p e = - 0.129 per additional year of lifespan). Using published vital rates for each species, we showed that the N e / N ratio resulting simply from life tables parameters can predict the observed differences in genetic diversity among species. Using simulations, we further found that the extent of reduction in N e / N with increasing adult lifespan is particularly strong under Type III survivorship curves (high juvenile and low adult mortality) and increasing fecundity with age, a typical characteristic of marine fishes. Our study highlights the importance of vital rates as key determinants of species genetic diversity levels in nature.

Global 16S rRNA diversity of provannid snail endosymbionts from Indo-Pacific deep-sea hydrothermal vents.

Symbioses between invertebrate animals and chemosynthetic bacteria build the foundation of deep-sea hydrothermal ecosystems worldwide. Despite the importance of these symbioses for ecosystem functioning, the diversity of symbionts within and between host organisms and geographic regions is still poorly understood. In this study we used 16S rRNA amplicon sequencing to determine the diversity of gill endosymbionts in provannid snails of the genera Alviniconcha and Ifremeria, which are key species at deep-sea hydrothermal vents in the Indo-Pacific Ocean. Our analysis of 761 snail samples across the distributional range of these species confirms previous findings that symbiont lineages are strongly partitioned by host species and broad-scale geography. Less structuring was observed within geographic regions, probably due to insufficient strain resolution of the 16S rRNA gene. Symbiont richness in individual hosts appeared to be unrelated to host size, suggesting that provannid snails might acquire their symbionts only during a permissive time window in early developmental stages in contrast to other vent molluscs that obtain their symbionts throughout their lifetime. Despite the extent of our dataset, symbiont accumulation curves did not reach saturation, highlighting the need for increased sampling efforts to uncover the full diversity of symbionts within these and other hydrothermal vent species.

Sexual isolation with and without ecological isolation in marine isopods Jaera albifrons and J. praehirsuta.

Sexual barriers associated with mate choice are often found to be associated with some level of ecological isolation between species. The independence and relative strength of sexual isolation are thus difficult to assess. Here, we take advantage of a pair of marine isopod species (Jaera albifrons and J. praehirsuta) that show sexual isolation and coexist in populations where they share the same microhabitat or not (i.e. without or with ecological isolation). We estimated the strength of sexual isolation between J. albifrons and J. praehirsuta using no-choice trials and a multiple-choice experimental population. We found that sexual isolation is strong in both the presence and the absence of ecological isolation, but that it is asymmetric and fails to prevent interspecific gene flow entirely. First-generation intrinsic post-zygotic barriers were low, and there was no sexual isolation within J. praehirsuta across habitats. The J. albifrons/J. praehirsuta species pair thus provides an example where the role of sexual isolation as a barrier to gene flow (a) does not depend upon current ecological isolation, (b) seems to have evolved independently of local ecological conditions, but (c) is insufficient to complete speciation entirely on its own.

The evolution of strong reproductive isolation between sympatric intertidal snails.

The evolution of strong reproductive isolation (RI) is fundamental to the origins and maintenance of biological diversity, especially in situations where geographical distributions of taxa broadly overlap. But what is the history behind strong barriers currently acting in sympatry? Using whole-genome sequencing and single nucleotide polymorphism genotyping, we inferred (i) the evolutionary relationships, (ii) the strength of RI, and (iii) the demographic history of divergence between two broadly sympatric taxa of intertidal snail. Despite being cryptic, based on external morphology, Littorina arcana and Littorina saxatilis differ in their mode of female reproduction (egg-laying versus brooding), which may generate a strong post-zygotic barrier. We show that egg-laying and brooding snails are closely related, but genetically distinct. Genotyping of 3092 snails from three locations failed to recover any recent hybrid or backcrossed individuals, confirming that RI is strong. There was, however, evidence for a very low level of asymmetrical introgression, suggesting that isolation remains incomplete. The presence of strong, asymmetrical RI was further supported by demographic analysis of these populations. Although the taxa are currently broadly sympatric, demographic modelling suggests that they initially diverged during a short period of geographical separation involving very low gene flow. Our study suggests that some geographical separation may kick-start the evolution of strong RI, facilitating subsequent coexistence of taxa in sympatry. The strength of RI needed to achieve sympatry and the subsequent effect of sympatry on RI remain open questions. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Geographically distinct patterns of reproductive isolation and hybridization in two sympatric species of the Jaera albifrons complex (marine isopods).

Sympatric species that in some populations hybridize and in other populations remain reproductively isolated open interesting research possibilities for the study of hybridization and speciation. Here, we test for such a situation in two littoral isopods (Jaera albifrons and J. praehirsuta) that occur in mixed populations and where past morphological descriptions suggested that the two species are generally reproductively isolated except in rare populations where hybridization may be happening. Using field surveys and microsatellite genetic structure analyses in two regions from France (Normandy and Brittany), we confirmed that introgressive hybridization occurs in a subset of mixed J. albifrons/J. praehirsuta populations (region Normandy) where the two species are found in the same habitat (pebbles on the shore). Moreover, we found that introgression in these populations is differential, 21 of 23 microsatellite markers showing little genetic divergence between species (hierarchical analysis of molecular variance FCT = 0.017) while the remaining two loci were strongly differentiated (FCT = 0.428). By contrast, J. albifrons and J. praehirsuta in mixed populations from region Brittany occupied distinct habitats (pebbles and seaweeds, respectively) with little overlap and showed stronger genetic divergence (FCT = 0.132). In hybridizing populations, the majority of individuals show morphological traits that are characteristic of one or the other species. This raises the question of the forces that act to maintain this polymorphism, noting that hybridizing populations seem to be geographically isolated from potential source parental populations and show no detectable habitat divergence between species.

Unexpected collective larval dispersal but little support for sweepstakes reproductive success in the highly dispersive brooding mollusc Crepidula fornicata.

In many marine invertebrates, long-distance dispersal is achieved during an extended pelagic larval phase. Although such dispersal should result in high gene flow over broad spatial scales, fine-scale genetic structure has often been reported, a pattern attributed to interfamilial variance in reproductive success and limited homogenization during dispersal. To examine this hypothesis, the genetic diversity of dispersing larvae must be compared with the postdispersal stages, that is benthic recruits and adults. Such data remain, however, scarce due to the difficulty to sample and analyse larvae of minute size. Here, we carried out such an investigation using the marine gastropod Crepidula fornicata. Field sampling of three to four larval pools was conducted over the reproductive season and repeated over 3 years. The genetic composition of larval pools, obtained with 16 microsatellite loci, was compared with that of recruits and adults sampled from the same site and years. In contrast to samples of juveniles and adults, large genetic temporal variations between larval pools produced at different times of the same reproductive season were observed. In addition, full- and half-sibs were detected in early larvae and postdispersal juveniles, pointing to correlated dispersal paths between several pairs of individuals. Inbred larvae were also identified. Such collective larval dispersal was unexpected given the long larval duration of the study species. Our results suggest that each larval pool is produced by a small effective number of reproducers but that, over a reproductive season, the whole larval pool is produced by large numbers of reproducers across space and time.

Current hypotheses to explain genetic chaos under the sea.

Chaotic genetic patchiness (CGP) refers to surprising patterns of spatial and temporal genetic structure observed in some marine species at a scale where genetic variation should be efficiently homogenized by gene flow via larval dispersal. Here we review and discuss 4 mechanisms that could generate such unexpected patterns: selection, sweepstakes reproductive success, collective dispersal, and temporal shifts in local population dynamics. First, we review examples where genetic differentiation at specific loci was driven by diversifying selection, which was historically the first process invoked to explain CGP. Second, we turn to neutral demographic processes that may drive genome-wide effects, and whose effects on CGP may be enhanced when they act together. We discuss how sweepstakes reproductive success accelerates genetic drift and can thus generate genetic structure, provided that gene flow is not too strong. Collective dispersal is another mechanism whereby genetic structure can be maintained regardless of dispersal intensity, because it may prevent larval cohorts from becoming entirely mixed. Theoretical analyses of both the sweepstakes and the collective dispersal ideas are presented. Finally, we discuss an idea that has received less attention than the other ones just mentioned, namely temporal shifts in local population dynamics.

The Size Advantage Model of Sex Allocation in the Protandrous Sex-Changer Crepidula fornicata: Role of the Mating System, Sperm Storage, and Male Mobility.

Sequential hermaphroditism is adaptive when the reproductive value of an individual varies with size or age, and this relationship differs between males and females. In this case, theory shows that the lifetime reproductive output of an individual is increased by changing sex (a hypothesis referred to as the size-advantage model). Sex-linked differences in size-fitness curves can stem from differential costs of reproduction, the mating system, and differences in growth and mortality between sexes. Detailed empirical data is required to disentangle the relative roles of each of these factors within the theory. Quantitative data are also needed to explore the role of sperm storage, which has not yet been considered with sequential hermaphrodites. Using experimental rearing and paternity assignment, we report relationships between size and reproductive success of Crepidula fornicata, a protandrous (male-first) gastropod. Male reproductive success increased with size due to the polygamous system and stacking behavior of the species, but females nonetheless had greater reproductive success than males of the same size, in agreement with the size-advantage theory. Sperm storage appeared to be a critical determinant of success for both sexes, and modeling the effect of sperm storage showed that it could potentially accelerate sex change in protandrous species.

Microsatellite markers for the Jaera albifrons species complex (marine isopods).

BACKGROUND: The Jaera albifrons complex contains five species of marine isopods (J. albifrons, J. praehirsuta, J. ischiosetosa, J. forsmani, and J. posthirsuta). These species, occurring on the shores of the North-Atlantic Ocean, are partially reproductively isolated by barriers due to sexual isolation (mate choice), genetic incompatibilities, and ecological specialization. Microsatellite loci would be useful for parentage-based analyses of sexual selection and studies of genetic structure in the context of speciation. FINDINGS: Twenty-four microsatellite markers were developed for J. albifrons using pyrosequencing of enriched libraries. Patterns of polymorphisms were analyzed in 49 J. albifrons adult males sampled in two populations from Brittany (Western France). The average number of alleles per locus was 4.73 ± 2.45 and the average gene diversity was 0.55 ± 0.23. Most markers also successfully amplified in the three sibling species J. praehirsuta, J. ischiosetosa, and J. forsmani. CONCLUSIONS: These polymorphic and cross-amplifiable markers will be useful for population genetics and parentage studies in the J albifrons complex.

Using neutral, selected, and hitchhiker loci to assess connectivity of marine populations in the genomic era.

Estimating the rate of exchange of individuals among populations is a central concern to evolutionary ecology and its applications to conservation and management. For instance, the efficiency of protected areas in sustaining locally endangered populations and ecosystems depends on reserve network connectivity. The population genetics theory offers a powerful framework for estimating dispersal distances and migration rates from molecular data. In the marine realm, however, decades of molecular studies have met limited success in inferring genetic connectivity, due to the frequent lack of spatial genetic structure in species exhibiting high fecundity and dispersal capabilities. This is especially true within biogeographic regions bounded by well-known hotspots of genetic differentiation. Here, we provide an overview of the current methods for estimating genetic connectivity using molecular markers and propose several directions for improving existing approaches using large population genomic datasets. We highlight several issues that limit the effectiveness of methods based on neutral markers when there is virtually no genetic differentiation among samples. We then focus on alternative methods based on markers influenced by selection. Although some of these methodologies are still underexplored, our aim was to stimulate new research to test how broadly they are applicable to nonmodel marine species. We argue that the increased ability to apply the concepts of cline analyses will improve dispersal inferences across physical and ecological barriers that reduce connectivity locally. We finally present how neutral markers hitchhiking with selected loci can also provide information about connectivity patterns within apparently well-mixed biogeographic regions. We contend that one of the most promising applications of population genomics is the use of outlier loci to delineate relevant conservation units and related eco-geographic features across which connectivity can be measured.

The effect of collective dispersal on the genetic structure of a subdivided population.

Correlated dispersal paths between two or more individuals are widespread across many taxa. The population genetic implications of this collective dispersal have received relatively little attention. Here we develop two-sample coalescent theory that incorporates collective dispersal in a finite island model to predict expected coalescence times, genetic diversities, and F-statistics. We show that collective dispersal reduces mixing in the system, which decreases expected coalescence times and increases FST . The effects are strongest in systems with high migration rates. Collective dispersal breaks the invariance of within-deme coalescence times to migration rate, whatever the deme size. It can also cause FST to increase with migration rate because the ratio of within- to between-deme coalescence times can decrease as migration rate approaches unity. This effect is most biologically relevant when deme size is small. We find qualitatively similar results for diploid and gametic dispersal. We also demonstrate with simulations and analytical theory the strong similarity between the effects of collective dispersal and anisotropic dispersal. These findings have implications for our understanding of the balance between drift-migration-mutation in models of neutral evolution. This has applied consequences for the interpretation of genetic structure (e.g., chaotic genetic patchiness) and estimation of migration rates from genetic data.

Genetic drift and collective dispersal can result in chaotic genetic patchiness.

Chaotic genetic patchiness denotes unexpected patterns of genetic differentiation that are observed at a fine scale and are not stable in time. These patterns have been described in marine species with free-living larvae, but are unexpected because they occur at a scale below the dispersal range of pelagic larvae. At the scale where most larvae are immigrants, theory predicts spatially homogeneous, temporally stable genetic variation. Empirical studies have suggested that genetic drift interacts with complex dispersal patterns to create chaotic genetic patchiness. Here we use a co-ancestry model and individual-based simulations to test this idea. We found that chaotic genetic patterns (qualified by global FST and spatio-temporal variation in FST's between pairs of samples) arise from the combined effects of (1) genetic drift created by the small local effective population sizes of the sessile phase and variance in contribution among breeding groups and (2) collective dispersal of related individuals in the larval phase. Simulations show that patchiness levels qualitatively comparable to empirical results can be produced by a combination of strong variance in reproductive success and mild collective dispersal. These results call for empirical studies of the effective number of breeders producing larval cohorts, and population genetics at the larval stage.

Scale-specific sex-biased dispersal in the Valais shrew unveiled by genetic variation on the Y chromosome, autosomes, and mitochondrial DNA.

We investigated sex specificities in the evolutionary processes shaping Y chromosome, autosomes, and mitochondrial DNA patterns of genetic structure in the Valais shrew (Sorex antinorii), a mountain dwelling species with a hierarchical distribution. Both hierarchical analyses of variance and isolation-by-distance analyses revealed patterns of population structure that were not consistent across maternal, paternal, and biparentally inherited markers. Differentiation on a Y microsatellite was lower than expected from the comparison with autosomal microsatellites and mtDNA, and it was mostly due to genetic variance among populations within valleys, whereas the opposite was observed on other markers. In addition, there was no pattern of isolation by distance for the Y, whereas there was strong isolation by distance on mtDNA and autosomes. We use a hierarchical island model of coancestry dynamics to discuss the relative roles of the microevolutionary forces that may induce such patterns. We conclude that sex-biased dispersal is the most important driver of the observed genetic structure, but with an intriguing twist: it seems that dispersal is strongly male biased at large spatial scale, whereas it is mildly biased in favor of females at local scale. These results add to recent reports of scale-specific sex-biased dispersal patterns, and emphasize the usefulness of the Y chromosome in conjunction with mtDNA and autosomes to infer sex specificities.

Genetic bottlenecks driven by population disconnection.

Connectivity among populations plays a crucial role in maintaining genetic variation at a local scale, especially in small populations affected strongly by genetic drift. The negative consequences of population disconnection on allelic richness and gene diversity (heterozygosity) are well recognized and empirically established. It is not well recognized, however, that a sudden drop in local effective population size induced by such disconnection produces a temporary disequilibrium in allelic frequency distributions that is akin to the genetic signature of a demographic bottleneck. To document this effect, we used individual-based simulations and empirical data on allelic richness and gene diversity in six pairs of isolated versus well-connected (core) populations of European tree frogs. In our simulations, population disconnection depressed allelic richness more than heterozygosity and thus resulted in a temporary excess in gene diversity relative to mutation drift equilibrium (i.e., signature of a genetic bottleneck). We observed a similar excess in gene diversity in isolated populations of tree frogs. Our results show that population disconnection can create a genetic bottleneck in the absence of demographic collapse.

Good genes drive female choice for mating partners in the lek-breeding European treefrog.

Investigating the mechanisms underlying female mate choice is important for sexual-selection theory, but also for population-genetic studies, because distinctive breeding strategies affect differently the dynamics of gene diversity within populations. Using field-monitoring, genetic-assignment, and laboratory-rearing methods, we investigated chorus attendance, mating success and offspring fitness in a population of lek-breeding tree-frogs (Hyla arborea) to test whether female choice is driven by good genes or complementary genes. Chorus attendance explained approximately 50% of the variance in male mating success, but did not correlate with offspring fitness. By contrast, offspring body mass and growth rate correlated with male attractiveness, measured as the number of matings obtained per night of calling. Genetic similarity between mating partners did not depart from random, and did not affect offspring fitness. We conclude that females are able to choose good partners under natural settings and obtain benefits from the good genes, rather than compatible genes, their offspring inherit. This heritability of fitness is likely to reduce effective population sizes below values previously estimated.

Inferring recent migration rates from individual genotypes.

We present a novel and straightforward method for estimating recent migration rates between discrete populations using multilocus genotype data. The approach builds upon a two-step sampling design, where individual genotypes are sampled before and after dispersal. We develop a model that estimates all pairwise backwards migration rates (m(ij), the probability that an individual sampled in population i is a migrant from population j) between a set of populations. The method is validated with simulated data and compared with the methods of BayesAss and Structure. First, we use data for an island model and then we consider more realistic data simulations for a metapopulation of the greater white-toothed shrew (Crocidura russula). We show that the precision and bias of estimates primarily depend upon the proportion of individuals sampled in each population. Weak sampling designs may particularly affect the quality of the coverage provided by 95% highest posterior density intervals. We further show that it is relatively insensitive to the number of loci sampled and the overall strength of genetic structure. The method can easily be extended and makes fewer assumptions about the underlying demographic and genetic processes than currently available methods. It allows backwards migration rates to be estimated across a wide range of realistic conditions.