Population genomic footprints of fine-scale differentiation between habitats in Mediterranean blue tits

Fine-scale genetic structure and phenotypic divergence of a passerine bird population inhabiting a continuous Mediterranean woodland

Genetic differentiation between populations inhabiting ecologically different habitats might appear because of limited dispersal and gene flow, which may lead to patterns of phenotypic divergence and local adaptation. In this study, we use dispersal, genotypic (24 microsatellite loci) and phenotypic (body size and clutch size) data to analyse patterns of genetic structuring and phenotypic divergence in a blue tit (Cyanistes caeruleus) population inhabiting a continuous and heterogeneous woodland along a valley. The two slopes of the valley differ in their forest formations and environmental conditions. Findings showed that most blue tits reproduced within their natal slope. Accordingly, microsatellite analyses revealed that populations of blue tits established in the two slopes show subtle genetic differentiation. The two genetic populations diverged in clutch size, exceeding the level of differentiation expected based on genetic drift, hence suggesting divergent selection (or other processes promoting divergence) on this life-history trait. Our findings reveal that restricted dispersal and spatial heterogeneity may lead to genetic differentiation among bird populations at a surprisingly small scale. In this respect, it is worth highlighting that such differentiation occurs for an organism with high dispersal capacity and within a continuous woodland. Moreover, we show that small-scale ecological differences, together with limited gene flow, can result in selection favouring different phenotypes even within the same continuum population.

Small-Mammal Genomics Highlights Viaducts as Potential Dispersal Conduits for Fragmented Populations

Wildlife crossings are implemented in many countries to facilitate the dispersal of animals among habitats fragmented by roads. However, the efficacy of different types of habitat corridors remains poorly understood. We used a comprehensive sampling regime in two lowland dipterocarp forest areas in peninsular Malaysia to sample pairs of small mammal individuals in three treatment types: (1) viaduct sites, at which sampling locations were separated by a highway but connected by a vegetated viaduct; (2) non-viaduct sites, at which sampling locations were separated by a highway and not connected by a viaduct; and (3) control sites, at which there was no highway fragmenting the forest. For four small mammal species, the common tree shrew Tupaia glis, Rajah's spiny rat Maxomys rajah, Whitehead's spiny rat Maxomys whiteheadi and dark-tailed tree rat Niviventer cremoriventer, we used genome-wide markers to assess genetic diversity, gene flow and genetic structure. The differences in genetic distance across sampling settings among the four species indicate that they respond differently to the presence of highways and viaducts. Viaducts connecting forests separated by highways appear to maintain higher population connectivity than forest fragments without viaducts, at least in M. whiteheadi, but apparently not in the other species.

Adaptive phenotypic and genomic divergence in the common chaffinch (Fringilla coelebs) following niche expansion within a small oceanic island

According to models of ecological speciation, adaptation to adjacent, contrasting habitat types can lead to population divergence given strong enough environment-driven selection to counteract the homogenizing effect of gene flow. We tested this hypothesis in the common chaffinch (Fringilla coelebs) on the small island of La Palma, Canary Islands, where it occupies two markedly different habitats. Isotopic (δ13 C, δ15 N) analysis of feathers indicated that birds in the two habitats differed in ecosystem and/or diet, and analysis of phenotypic traits revealed significant differences in morphology and plumage colouration that are consistent with ecomorphological and ecogeographical predictions respectively. A genome-wide survey of single-nucleotide polymorphism revealed marked neutral structure that was consistent with geography and isolation by distance, suggesting low dispersal. In contrast, loci putatively under selection identified through genome-wide association and genotype-environment association analyses, revealed amarked adaptive divergence between birds in both habitats. Loci associated with phenotypic and environmental differences among habitats were distributed across the genome, as expected for polygenic traits involved in local adaptation. Our results suggest a strong role for habitat-driven local adaptation in population divergence in the chaffinches of La Palma, a process that appears to be facilitated by a strong reduction in effective dispersal distances despite the birds' high dispersal capacity.

Genetic diversity and spatial genetic structure support the specialist-generalist variation hypothesis in two sympatric woodpecker species

Species are often arranged along a continuum from “specialists” to “generalists”. Specialists typically use fewer resources, occur in more patchily distributed habitats and have overall smaller population sizes than generalists. Accordingly, the specialist-generalist variation hypothesis (SGVH) proposes that populations of habitat specialists have lower genetic diversity and are genetically more differentiated due to reduced gene flow compared to populations of generalists. Here, expectations of the SGVH were tested by examining genetic diversity, spatial genetic structure and contemporary gene flow in two sympatric woodpecker species differing in habitat specialization. Compared to the generalist great spotted woodpecker (Dendrocopos major), lower genetic diversity was found in the specialist middle spotted woodpecker (Dendrocoptes medius). Evidence for recent bottlenecks was revealed in some populations of the middle spotted woodpecker, but in none of the great spotted woodpecker. Substantial spatial genetic structure and a significant correlation between genetic and geographic distances were found in the middle spotted woodpecker, but only weak spatial genetic structure and no significant correlation between genetic and geographic distances in the great spotted woodpecker. Finally, estimated levels of contemporary gene flow did not differ between the two species. Results are consistent with all but one expectations of the SGVH. This study adds to the relatively few investigations addressing the SGVH in terrestrial vertebrates.

Interactions between microenvironment, selection and genetic architecture drive multiscale adaptation in a simulation experiment

When environmental conditions differ both within and among populations, multiscale adaptation results from processes at both scales and interference across scales. We hypothesize that within-population environmental heterogeneity influences the chance of success of migration events, both within and among populations, and maintains within-population adaptive differentiation. We used a simulation approach to analyze the joint effects of environmental heterogeneity patterns, selection intensity and number of QTL controlling a selected trait on local adaptation in a hierarchical metapopulation design. We show the general effects of within-population environmental heterogeneity: (i) it increases occupancy rate at the margins of distribution ranges, under extreme environments and high levels of selection; (ii) it increases the adaptation lag in all environments; (iii) it impacts the genetic variance in each environment, depending on the ratio of within- to between-populations environmental heterogeneity; (iv) it reduces the selection-induced erosion of adaptive gene diversity. Most often, the smaller the number of QTL involved, the stronger are these effects. We also show that both within- and between-populations phenotypic differentiation (Q_ST ) mainly results from covariance of QTL effects rather than QTL differentiation (F_STq ), that within-population QTL differentiation is negligible, and that stronger divergent selection is required to produce adaptive differentiation within populations than among populations. With a high number of QTL, when the difference between environments within populations exceeds the smallest difference between environments across populations, high levels of within-population differentiation can be reached, reducing differentiation among populations. Our study stresses the need to account for within-population environmental heterogeneity when investigating local adaptation.

Divergence of a genomic island leads to the evolution of melanization in a halophyte root fungus

Understanding how organisms adapt to extreme living conditions is central to evolutionary biology. Dark septate endophytes (DSEs) constitute an important component of the root mycobiome and they are often able to alleviate host abiotic stresses. Here, we investigated the molecular mechanisms underlying the beneficial association between the DSE Laburnicola rhizohalophila and its host, the native halophyte Suaeda salsa, using population genomics. Based on genome-wide Fst (pairwise fixation index) and Vst analyses, which compared the variance in allele frequencies of single-nucleotide polymorphisms (SNPs) and copy number variants (CNVs), respectively, we found a high level of genetic differentiation between two populations. CNV patterns revealed population-specific expansions and contractions. Interestingly, we identified a ~20 kbp genomic island of high divergence with a strong sign of positive selection. This region contains a melanin-biosynthetic polyketide synthase gene cluster linked to six additional genes likely involved in biosynthesis, membrane trafficking, regulation, and localization of melanin. Differences in growth yield and melanin biosynthesis between the two populations grown under 2% NaCl stress suggested that this genomic island contributes to the observed differences in melanin accumulation. Our findings provide a better understanding of the genetic and evolutionary mechanisms underlying the adaptation to saline conditions of the L. rhizohalophila-S. salsa symbiosis.

Genomic Characterization of the Barnacle Balanus improvisus Reveals Extreme Nucleotide Diversity in Coding Regions

Barnacles are key marine crustaceans in several habitats, and they constitute a common practical problem by causing biofouling on man-made marine constructions and ships. Despite causing considerable ecological and economic impacts, there is a surprising void of basic genomic knowledge, and a barnacle reference genome is lacking. We here set out to characterize the genome of the bay barnacle Balanus improvisus (= Amphibalanus improvisus) based on short-read whole-genome sequencing and experimental genome size estimation. We show both experimentally (DNA staining and flow cytometry) and computationally (k-mer analysis) that B. improvisus has a haploid genome size of ~ 740 Mbp. A pilot genome assembly rendered a total assembly size of ~ 600 Mbp and was highly fragmented with an N50 of only 2.2 kbp. Further assembly-based and assembly-free analyses revealed that the very limited assembly contiguity is due to the B. improvisus genome having an extremely high nucleotide diversity (π) in coding regions (average π ≈ 5% and average π in fourfold degenerate sites ≈ 20%), and an overall high repeat content (at least 40%). We also report on high variation in the α-octopamine receptor OctA (average π = 3.6%), which might increase the risk that barnacle populations evolve resistance toward antifouling agents. The genomic features described here can help in planning for a future high-quality reference genome, which is urgently needed to properly explore and understand proteins of interest in barnacle biology and marine biotechnology and for developing better antifouling strategies.

Genetic structure of urban and non-urban populations differs between two common parid species

Landscape conversions induced by human activities can affect dispersal patterns of various bird species and, as a result, affect genetic structure of their populations. Genetic differentiation of bird populations may be enhanced by habitat variation, especially in urban-non-urban systems. The majority of population genetic studies focus on single species, which inflicts limitations for direct comparisons of genetic responses of avian populations to urbanization. Here, we used a set of microsatellite markers to examine genetic diversity, gene flow and population structure in two common parid species, great tits Parus major and blue tits Cyanistes caeruleus occupying three sites in habitats with contrasting urbanization level in central Poland. We found low but significant divergence of urban park population with both suburban and non-urban forest great tit populations, while no differentiation was found between suburban forest and non-urban forest populations. In contrast, no evidence for genetic differentiation was found between blue tit populations from the urban park, suburban forest and non-urban forest sites. We conclude that great tits and blue tits respond to urbanization-related changes in a different way, which may be a result of different rates of migration and/or dispersal, likely higher in blue tits. Some impact may be also induced by interspecific competition. We suggest that changing the focus of urban genetic research from single to multiple species may provide novel insights into how natural populations respond to the processes of urbanization.

Genomic variation in the American pika: signatures of geographic isolation and implications for conservation

Background

Distributional responses by alpine taxa to repeated, glacial-interglacial cycles throughout the last two million years have significantly influenced the spatial genetic structure of populations. These effects have been exacerbated for the American pika (Ochotona princeps), a small alpine lagomorph constrained by thermal sensitivity and a limited dispersal capacity. As a species of conservation concern, long-term lack of gene flow has important consequences for landscape genetic structure and levels of diversity within populations. Here, we use reduced representation sequencing (ddRADseq) to provide a genome-wide perspective on patterns of genetic variation across pika populations representing distinct subspecies. To investigate how landscape and environmental features shape genetic variation, we collected genetic samples from distinct geographic regions as well as across finer spatial scales in two geographically proximate mountain ranges of eastern Nevada.

Results

Our genome-wide analyses corroborate range-wide, mitochondrial subspecific designations and reveal pronounced fine-scale population structure between the Ruby Mountains and East Humboldt Range of eastern Nevada. Populations in Nevada were characterized by low genetic diversity (π = 0.0006–0.0009; θ_W = 0.0005–0.0007) relative to populations in California (π = 0.0014–0.0019; θ_W = 0.0011–0.0017) and the Rocky Mountains (π = 0.0025–0.0027; θ_W = 0.0021–0.0024), indicating substantial genetic drift in these isolated populations. Tajima’s D was positive for all sites (D = 0.240–0.811), consistent with recent contraction in population sizes range-wide.

Conclusions

Substantial influences of geography, elevation and climate variables on genetic differentiation were also detected and may interact with the regional effects of anthropogenic climate change to force the loss of unique genetic lineages through continued population extirpations in the Great Basin and Sierra Nevada.

Demographic history and genomics of local adaptation in blue tit populations

Understanding the genomic processes underlying local adaptation is a central aim of modern evolutionary biology. This task requires identifying footprints of local selection but also estimating spatio‐temporal variations in population demography and variations in recombination rate and in diversity along the genome. Here, we investigated these parameters in blue tit populations inhabiting deciduous versus evergreen forests, and insular versus mainland areas, in the context of a previously described strong phenotypic differentiation. Neighboring population pairs of deciduous and evergreen habitats were weakly genetically differentiated (F_ST = 0.003 on average), nevertheless with a statistically significant effect of habitat type on the overall genetic structure. This low differentiation was consistent with the strong and long‐lasting gene flow between populations inferred by demographic modeling. In turn, insular and mainland populations were moderately differentiated (F_ST = 0.08 on average), in line with the inference of moderate ancestral migration, followed by isolation since the end of the last glaciation. Effective population sizes were large, yet smaller on the island than on the mainland. Weak and nonparallel footprints of divergent selection between deciduous and evergreen populations were consistent with their high connectivity and the probable polygenic nature of local adaptation in these habitats. In turn, stronger footprints of divergent selection were identified between long isolated insular versus mainland birds and were more often found in regions of low recombination, as expected from theory. Lastly, we identified a genomic inversion on the mainland, spanning 2.8 Mb. These results provide insights into the demographic history and genetic architecture of local adaptation in blue tit populations at multiple geographic scales.

Fine-scale barriers to connectivity across a fragmented South-East Asian landscape in six songbird species

Habitat fragmentation is a major extinction driver. Despite dramatically increasing fragmentation across the globe, its specific impacts on population connectivity across species with differing life histories remain difficult to characterize, let alone quantify. Here, we investigate patterns of population connectivity in six songbird species from Singapore, a highly fragmented tropical rainforest island. Using massive panels of genome-wide single nucleotide polymorphisms across dozens of samples per species, we examined population genetic diversity, inbreeding, gene flow and connectivity among species along a spectrum of ecological specificities. We found a higher resilience to habitat fragmentation in edge-tolerant and forest-canopy species as compared to forest-dependent understorey insectivores. The latter exhibited levels of genetic diversity up to three times lower in Singapore than in populations from contiguous forest elsewhere. Using dense genomic and geographic sampling, we identified individual barriers such as reservoirs that effectively minimize gene flow in sensitive understorey birds, revealing that terrestrial forest species may exhibit levels of sensitivity to fragmentation far greater than previously expected. This study provides a blueprint for conservation genomics at small scales with a view to identifying preferred locations for habitat corridors, flagging candidate populations for restocking with translocated individuals and improving the design of future reserves.

Phenotypic plasticity in response to climate change: the importance of cue variation

Phenotypic plasticity is a major mechanism of response to global change. However, current plastic responses will only remain adaptive under future conditions if informative environmental cues are still available. We briefly summarize current knowledge of the evolutionary origin and mechanistic underpinnings of environmental cues for phenotypic plasticity, before highlighting the potentially complex effects of global change on cue availability and reliability. We then illustrate some of these aspects with a case study, comparing plasticity of blue tit breeding phenology in two contrasted habitats: evergreen and deciduous forests. Using long-term datasets, we investigate the climatic factors linked to the breeding phenology of the birds and their main food source. Blue tits occupying different habitats differ extensively in the cues affecting laying date plasticity, as well as in the reliability of these cues as predictors of the putative driver of selective pressure, the date of caterpillar peak. The temporal trend for earlier laying date, detected only in the evergreen populations, is explained by increased temperature during their cue windows. Our results highlight the importance of integrating ecological mechanisms shaping variation in plasticity if we are to understand how global change will affect plasticity and its consequences for population biology. This article is part of the theme issue 'The role of plasticity in phenotypic adaptation to rapid environmental change'.

Population genomics of Bombus terrestris reveals high but unstructured genetic diversity in a potential glacial refugium

Ongoing climate change is expected to cause an increase in temperature and a reduction of precipitation levels in the Mediterranean region, which might cause changes in many species distributions. These effects negatively influence species gene pools, decreasing genetic variability and adaptive potential. Here, we use mitochondrial DNA and RADseq to analyse population genetic structure and genetic diversity of the bumblebee species Bombus terrestris (subspecies Bombus terrestris lusitanicus), in the Iberian Peninsula. Although this subspecies shows a panmictic pattern of population structure across Iberia and beyond, we found differentiation between subspecies B. t. lusitanicus and B. t. africanus, probably caused by the existence of barriers to gene flow between Iberia and North Africa. Furthermore, the results revealed that the Iberian Peninsula harbours a large fraction of B. terrestris intraspecific genetic variation, with the highest number of mitochondrial haplotypes found when compared with any other region in Europe studied so far, suggesting a potential role for the Iberian Peninsula as a glacial refugium. Our findings strengthen the idea that Iberia is a very important source of diversity for the global genetic pool of this species, because rare alleles might play a role in population resilience against human- or climate-mediated changes.

New insights from RADseq data on differentiation in the Hottentot golden mole species complex from South Africa

Golden moles (Family Chrysochloridae) are small subterranean mammals, endemic to sub-Saharan Africa, and many of the 21 species are listed as threatened on the IUCN Red List. Most species have highly restricted ranges; however two species, the Hottentot golden mole (Amblysomus hottentotus) and the Cape golden mole (Chrysochloris asiatica) have relatively wide ranges. We recently uncovered cryptic diversity within A. hottentotus, through a phylogeographic analysis of this taxon using two mitochondrial gene regions and a nuclear intron. To further investigate this cryptic diversity, we generated nuclear SNP data from across the genome of A. hottentotus, by means of double-digest restriction-site associated DNA sequencing (ddRADSeq), and mapped reads to the Cape golden mole genome. We conducted a phylogenetic analysis and investigated population differentiation. Our results support the distinctiveness of A. h. meesteri. Furthermore, we provide evidence from nuclear SNPs in support of our previous finding that Central coastal samples represent a unique cryptic lineage that is highly divergent from A. h. pondoliae farther south. Although mtDNA suggests that Umtata may represent a unique lineage sister to A. h. longiceps, mito-nuclear discordance from our RADseq data indicate that these samples may instead be closer to A. h. pondoliae, and therefore may not represent a distinct lineage. We stress the importance of recognizing that understudied populations, such as that of Umtata, may represent populations or ESUs under threat and in need of conservation attention. We present a high-quality filtered SNP dataset, comprising thousands of SNPs, which may serve as a useful resource for future golden mole studies. We have thus added to the growing body of research demonstrating the power and utility of RADseq to investigate population differentiation.

Evaluation of a Chicken 600K SNP genotyping array in non-model species of grouse

The use of single nucleotide polymorphism (SNP) arrays to generate large SNP datasets for comparison purposes have recently become an attractive alternative to other genotyping methods. Although most SNP arrays were originally developed for domestic organisms, they can be effectively applied to wild relatives to obtain large panels of SNPs. In this study, we tested the cross-species application of the Affymetrix 600K Chicken SNP array in five species of North American prairie grouse (Centrocercus and Tympanuchus genera). Two individuals were genotyped per species for a total of ten samples. A high proportion (91%) of the total 580 961 SNPs were genotyped in at least one individual (73–76% SNPs genotyped per species). Principal component analysis with autosomal SNPs separated the two genera, but failed to clearly distinguish species within genera. Gene ontology analysis identified a set of genes related to morphogenesis and development (including genes involved in feather development), which may be primarily responsible for large phenotypic differences between Centrocercus and Tympanuchus grouse. Our study provided evidence for successful cross-species application of the chicken SNP array in grouse which diverged ca. 37 mya from the chicken lineage. As far as we are aware, this is the first reported application of a SNP array in non-passerine birds, and it demonstrates the feasibility of using commercial SNP arrays in research on non-model bird species.

Genomic insights into the vulnerability of sympatric whitefish species flocks

The erosion of habitat heterogeneity can reduce species diversity directly but can also lead to the loss of distinctiveness of sympatric species through speciation reversal. We know little about changes in genomic differentiation during the early stages of these processes, which can be mediated by anthropogenic perturbation. Here, we analyse three sympatric whitefish species (Coregonus spp) sampled across two neighbouring and connected Swiss pre-alpine lakes, which have been differentially affected by anthropogenic eutrophication. Our data set comprises 16,173 loci genotyped across 138 whitefish using restriction-site associated DNA sequencing (RADseq). Our analysis suggests that in each of the two lakes, the population of a different, but ecologically similar, whitefish species declined following a recent period of eutrophication. Genomic signatures consistent with hybridization are more pronounced in the more severely impacted lake. Comparisons between sympatric pairs of whitefish species with contrasting ecology, where one is shallow benthic and the other one more profundal pelagic, reveal genomic differentiation that is largely correlated along the genome, while differentiation is uncorrelated between pairs of allopatric provenance with similar ecology. We identify four genomic loci that provide evidence of parallel divergent adaptation between the shallow benthic species and the two different more profundal species. Functional annotations available for two of those loci are consistent with divergent ecological adaptation. Our genomic analysis indicates the action of divergent natural selection between sympatric whitefish species in pre-alpine lakes and reveals the vulnerability of these species to anthropogenic alterations of the environment and associated adaptive landscape.

Genome-wide signals of drift and local adaptation during rapid lineage divergence in a songbird

The formation of independent evolutionary lineages involves neutral and selective factors, and understanding their relative roles in population divergence is a fundamental goal of speciation research. Correlations between allele frequencies and environmental variability can reveal the role of selection, yet the relative contribution of drift can be difficult to establish. Recently diversified taxa like the Oregon junco (Aves, Passerellidae, Junco hyemalis oreganus) of western North America provide ideal scenarios to apply genetic-environment association analyses (GEA) while controlling for population structure. Analysis of genome-wide SNP loci revealed marked genetic structure consisting of differentiated populations in isolated, dry southern mountain ranges, and less divergent, recently expanded populations in humid northern latitudes. We used correlations between genomic and environmental variance to test for three specific modes of evolutionary divergence: (a) drift in geographic isolation, (b) differentiation along continuous selective gradients and (c) isolation-by-adaptation. We found evidence of strong drift in southern mountains, but also signals of local adaptation driven by temperature, precipitation, elevation and vegetation, especially when controlling for population history. We identified numerous variants under selection scattered across the genome, suggesting that local adaptation can promote rapid differentiation when acting over multiple independent loci.

Unravelling the Genetic Diversity among Cassava Bemisia tabaci Whiteflies Using NextRAD Sequencing

Bemisia tabaci threatens production of cassava in Africa through vectoring viruses that cause cassava mosaic disease (CMD) and cassava brown streak disease (CBSD). B. tabaci sampled from cassava in eight countries in Africa were genotyped using NextRAD sequencing, and their phylogeny and population genetics were investigated using the resultant single nucleotide polymorphism (SNP) markers. SNP marker data and short sequences of mitochondrial DNA cytochrome oxidase I (mtCOI) obtained from the same insect were compared. Eight genetically distinct groups were identified based on mtCOI, whereas phylogenetic analysis using SNPs identified six major groups, which were further confirmed by PCA and multidimensional analyses. STRUCTURE analysis identified four ancestral B. tabaci populations that have contributed alleles to the six SNP-based groups. Significant gene flows were detected between several of the six SNP-based groups. Evidence of gene flow was strongest for SNP-based groups occurring in central Africa. Comparison of the mtCOI and SNP identities of sampled insects provided a strong indication that hybrid populations are emerging in parts of Africa recently affected by the severe CMD pandemic. This study reveals that mtCOI is not an effective marker at distinguishing cassava-colonizing B. tabaci haplogroups, and that more robust SNP-based multilocus markers should be developed. Significant gene flows between populations could lead to the emergence of haplogroups that might alter the dynamics of cassava virus spread and disease severity in Africa. Continuous monitoring of genetic compositions of whitefly populations should be an essential component in efforts to combat cassava viruses in Africa.

Landscape genetics identifies streams and drainage infrastructure as dispersal corridors for an endangered wetland bird

Anthropogenic alterations to landscape structure and composition can have significant impacts on biodiversity, potentially leading to species extinctions. Population‐level impacts of landscape change are mediated by animal behaviors, in particular dispersal behavior. Little is known about the dispersal habits of rails (Rallidae) due to their cryptic behavior and tendency to occupy densely vegetated habitats. The effects of landscape structure on the movement behavior of waterbirds in general are poorly studied due to their reputation for having high dispersal abilities. We used a landscape genetic approach to test hypotheses of landscape effects on dispersal behavior of the Hawaiian gallinule (Gallinula galeata sandvicensis), an endangered subspecies endemic to the Hawaiian Islands. We created a suite of alternative resistance surfaces representing biologically plausible a priori hypotheses of how gallinules might navigate the landscape matrix and ranked these surfaces by their ability to explain observed patterns in genetic distance among 12 populations on the island of O`ahu. We modeled effective distance among wetland locations on all surfaces using both cumulative least‐cost‐path and resistance‐distance approaches and evaluated relative model performance using Mantel tests, a causal modeling approach, and the mixed‐model maximum‐likelihood population‐effects framework. Across all genetic markers, simulation methods, and model comparison metrics, surfaces that treated linear water features like streams, ditches, and canals as corridors for gallinule movement outperformed all other models. This is the first landscape genetic study on the movement behavior of any waterbird species to our knowledge. Our results indicate that lotic water features, including drainage infrastructure previously thought to be of minimal habitat value, contribute to habitat connectivity in this listed subspecies.

Gene flow does not prevent personality and morphological differentiation between two blue tit populations

Understanding the causes and consequences of population phenotypic divergence is a central goal in ecology and evolution. Phenotypic divergence among populations can result from genetic divergence, phenotypic plasticity or a combination of the two. However, few studies have deciphered these mechanisms for populations geographically close and connected by gene flow, especially in the case of personality traits. In this study, we used a common garden experiment to explore the genetic basis of the phenotypic divergence observed between two blue tit (Cyanistes caeruleus) populations inhabiting contrasting habitats separated by 25 km, for two personality traits (exploration speed and handling aggression), one physiological trait (heart rate during restraint) and two morphological traits (tarsus length and body mass). Blue tit nestlings were removed from their population and raised in a common garden for up to 5 years. We then compared adult phenotypes between the two populations, as well as trait-specific Q_st and F_st . Our results revealed differences between populations similar to those found in the wild, suggesting a genetic divergence for all traits. Q_st -F_st comparisons revealed that the trait divergences likely result from dissimilar selection patterns rather than from genetic drift. Our study is one of the first to report a Q_st -F_st comparison for personality traits and adds to the growing body of evidence that population genetic divergence is possible at a small scale for a variety of traits including behavioural traits.

Great tits and the city: Distribution of genomic diversity and gene-environment associations along an urbanization gradient

Urbanization is a growing concern challenging the evolutionary potential of wild populations by reducing genetic diversity and imposing new selection regimes affecting many key fitness traits. However, genomic footprints of urbanization have received little attention so far. Using RAD sequencing, we investigated the genomewide effects of urbanization on neutral and adaptive genomic diversity in 140 adult great tits Parus major collected in locations with contrasted urbanization levels (from a natural forest to highly urbanized areas of a city; Montpellier, France). Heterozygosity was slightly lower in the more urbanized sites compared to the more rural ones. Low but significant effect of urbanization on genetic differentiation was found, at the site level but not at the nest level, indicative of the geographic scale of urbanization impact and of the potential for local adaptation despite gene flow. Gene-environment association tests identified numerous SNPs with small association scores to urbanization, distributed across the genome, from which a subset of 97 SNPs explained up to 81% of the variance in urbanization, overall suggesting a polygenic response to selection in the urban environment. These findings open stimulating perspectives for broader applications of high-resolution genomic tools on other cities and larger sample sizes to investigate the consistency of the effects of urbanization on the spatial distribution of genetic diversity and the polygenic nature of gene-urbanization association.

Inferring spatial patterns and drivers of population divergence of Neolitsea sericea (Lauraceae), based on molecular phylogeography and landscape genomics

The relative roles of geography, climate and ecology in driving population divergence and (incipient) speciation has so far been largely neglected in studies addressing the evolution of East Asia's island flora. Here, we employed chloroplast and ribosomal DNA sequences and restriction site-associated DNA sequencing (RADseq) loci to investigate the phylogeography and drivers of population divergence of Neolitsea sericea. These data sets support the subdivision of N. sericea populations into the Southern and Northern lineages across the 'Tokara gap'. Two distinct sublineages were further identified for the Northern lineage of N. sericea from the RADseq data. RADseq was also used along with approximate Bayesian computation to show that the current distribution and differentiation of N. sericea populations resulted from a combination of relatively ancient migration and successive vicariant events that likely occurred during the mid to late Pleistocene. Landscape genomic analyses showed that, apart from geographic barriers, barrier, potentially local adaptation to different climatic conditions appears to be one of the major drivers for lineage diversification of N. sericea.

Genomic Signature of Adaptive Divergence despite Strong Nonadaptive Forces on Edaphic Islands: A Case Study of Primulina juliae

Both genetic drift and divergent selection are expected to be strong evolutionary forces driving population differentiation on edaphic habitat islands. However, the relative contribution of genetic drift and divergent selection to population divergence has rarely been tested simultaneously. In this study, restriction-site associated DNA-based population genomic analyses were applied to assess the relative importance of drift and divergent selection on population divergence of Primulina juliae, an edaphic specialist from southern China. All populations were found with low standing genetic variation, small effective population size (NE), and signatures of bottlenecks. Populations with the lowest genetic variation were most genetically differentiated from other populations and the extent of genetic drift increased with geographic distance from other populations. Together with evidence of isolation by distance, these results support neutral drift as a critical evolutionary driver. Nonetheless, redundancy analysis revealed that genomic variation is significantly associated with both edaphic habitats and climatic factors independently of spatial effects. Moreover, more genomic variation was explained by environmental factors than by geographic variables, suggesting that local adaptation might have played an important role in driving population divergence. Finally, outlier tests and environment association analyses identified 31 single-nucleotide polymorphisms as candidates for adaptive divergence. Among these candidates, 26 single-nucleotide polymorphisms occur in/near genes that potentially play a role in adaptation to edaphic specialization. This study has important implications that improve our understanding of the joint roles of genetic drift and adaptation in generating population divergence and diversity of edaphic specialists.

Genetic Diversity and Connectivity in Maurolicus muelleri in the Bay of Biscay Inferred from Thousands of SNP Markers

Mesopelagic fish are largely abundant poorly studied fish that are still intact, but which, due to their potentially great added value, will be imminently exploited by humans. Therefore, studies that provide information to anticipate the anthropogenic impact on this important resource are urgently needed. In particular, knowledge about their connectivity, potential adaptation and resilience are needed. This information can be obtained through the analysis of genome-wide markers which are now relatively easily and cost-efficiently discovered thanks to high-throughput sequencing technologies. Here, we have generated thousands of SNP markers in Maurolicus muelleri, based on the restriction-site associated DNA sequencing method, and preformed population connectivity and genetic diversity analyses in a subset of samples collected from the Bay of Biscay. Our study proves the method valid for obtaining genome-wide markers in this species and provides the first insights into the population genomics of M. muelleri. Importantly, the genomic resources developed here are made available for future studies and set the basics for additional endeavors on this issue.

Using NextRAD sequencing to infer movement of herbivores among host plants

Herbivores often move among spatially interspersed host plants, tracking high-quality resources through space and time. This dispersal is of particular interest for vectors of plant pathogens. Existing molecular tools to track such movement have yielded important insights, but often provide insufficient genetic resolution to infer spread at finer spatiotemporal scales. Here, we explore the use of Nextera-tagmented reductively-amplified DNA (NextRAD) sequencing to infer movement of a highly-mobile winged insect, the potato psyllid (Bactericera cockerelli), among host plants. The psyllid vectors the pathogen that causes zebra chip disease in potato (Solanum tuberosum), but understanding and managing the spread of this pathogen is limited by uncertainty about the insect's host plant(s) outside of the growing season. We identified 1,978 polymorphic loci among psyllids separated spatiotemporally on potato or in patches of bittersweet nightshade (S. dulcumara), a weedy plant proposed to be the source of potato-colonizing psyllids. A subset of the psyllids on potato exhibited genetic similarity to insects on nightshade, consistent with regular movement between these two host plants. However, a second subset of potato-collected psyllids was genetically distinct from those collected on bittersweet nightshade; this suggests that a currently unrecognized source, i.e., other nightshade patches or a third host-plant species, could be contributing to psyllid populations in potato. Oftentimes, dispersal of vectors of pathogens must be tracked at a fine scale in order to understand, predict, and manage disease spread. We demonstrate that emerging sequencing technologies that detect genome-wide SNPs of a vector can be used to infer such localized movement.

Conserved G-matrices of morphological and life-history traits among continental and island blue tit populations

The genetic variance-covariance matrix (G-matrix) summarizes the genetic architecture of multiple traits. It has a central role in the understanding of phenotypic divergence and the quantification of the evolutionary potential of populations. Laboratory experiments have shown that G-matrices can vary rapidly under divergent selective pressures. However, because of the demanding nature of G-matrix estimation and comparison in wild populations, the extent of its spatial variability remains largely unknown. In this study, we investigate spatial variation in G-matrices for morphological and life-history traits using long-term data sets from one continental and three island populations of blue tit (Cyanistes caeruleus) that have experienced contrasting population history and selective environment. We found no evidence for differences in G-matrices among populations. Interestingly, the phenotypic variance-covariance matrices (P) were divergent across populations, suggesting that using P as a substitute for G may be inadequate. These analyses also provide the first evidence in wild populations for additive genetic variation in the incubation period (that is, the period between last egg laid and hatching) in all four populations. Altogether, our results suggest that G-matrices may be stable across populations inhabiting contrasted environments, therefore challenging the results of previous simulation studies and laboratory experiments.

Absence of population structure across elevational gradients despite large phenotypic variation in mountain chickadees (Poecile gambeli)

Montane habitats are characterized by predictably rapid heterogeneity along elevational gradients and are useful for investigating the consequences of environmental heterogeneity for local adaptation and population genetic structure. Food-caching mountain chickadees inhabit a continuous elevation gradient in the Sierra Nevada, and birds living at harsher, high elevations have better spatial memory ability and exhibit differences in male song structure and female mate preference compared to birds inhabiting milder, low elevations. While high elevation birds breed, on average, two weeks later than low elevation birds, the extent of gene flow between elevations is unknown. Despite phenotypic variation and indirect evidence for local adaptation, population genetic analyses based on 18 073 single nucleotide polymorphisms across three transects of high and low elevation populations provided no evidence for genetic differentiation. Analyses based on individual genotypes revealed no patterns of clustering, pairwise estimates of genetic differentiation (F_ST, Nei's D) were very low, and AMOVA revealed no evidence for genetic variation structured by transect or by low and high elevation sites within transects. In addition, we found no consistent evidence for strong parallel allele frequency divergence between low and high elevation sites within the three transects. Large elevation-related phenotypic variation may be maintained by strong selection despite gene flow and future work should focus on the mechanisms underlying such variation.

Assortative mating by colored ornaments in blue tits: space and time matter

Assortative mating is a potential outcome of sexual selection, and estimating its level is important to better understand local adaptation and underlying trait evolution. However, assortative mating studies frequently base their conclusions on small numbers of individuals sampled over short periods of time and limited spatial scales even though spatiotemporal variation is common. Here, we characterized assortative mating patterns over 10 years in four populations of the blue tit (Cyanistes caeruleus), a passerine bird. We focused on two plumage ornaments-the blue crown and the yellow breast patch. Based on data for 1,657 pairs of birds, we found large interannual variation: assortative mating varied from positive to negative. To determine whether there was nonetheless a general trend in the data, we ran a within-study meta-analysis. It revealed that assortative mating was moderately positive for both ornaments. It also showed that mating patterns differed among populations and especially between two neighboring populations that displayed phenotypic divergence. Our results therefore underscore that long-term studies are needed to draw broad conclusions about mating patterns in natural populations. They also call for studying the potential role of assortative mating in local adaptation and evolution of ornaments in both sexes.

RAD sequencing resolves fine-scale population structure in a benthic invertebrate: implications for understanding phenotypic plasticity

The field of molecular ecology is transitioning from the use of small panels of classical genetic markers such as microsatellites to much larger panels of single nucleotide polymorphisms (SNPs) generated by approaches like RAD sequencing. However, few empirical studies have directly compared the ability of these methods to resolve population structure. This could have implications for understanding phenotypic plasticity, as many previous studies of natural populations may have lacked the power to detect genetic differences, especially over micro-geographic scales. We therefore compared the ability of microsatellites and RAD sequencing to resolve fine-scale population structure in a commercially important benthic invertebrate by genotyping great scallops (Pecten maximus) from nine populations around Northern Ireland at 13 microsatellites and 10 539 SNPs. The shells were then subjected to morphometric and colour analysis in order to compare patterns of phenotypic and genetic variation. We found that RAD sequencing was superior at resolving population structure, yielding higher Fst values and support for two distinct genetic clusters, whereas only one cluster could be detected in a Bayesian analysis of the microsatellite dataset. Furthermore, appreciable phenotypic variation was observed in size-independent shell shape and coloration, including among localities that could not be distinguished from one another genetically, providing support for the notion that these traits are phenotypically plastic. Taken together, our results suggest that RAD sequencing is a powerful approach for studying population structure and phenotypic plasticity in natural populations.

RAD sequencing resolves fine-scale population structure in a benthic invertebrate: implications for understanding phenotypic plasticity

The field of molecular ecology is transitioning from the use of small panels of classical genetic markers such as microsatellites to much larger panels of single nucleotide polymorphisms (SNPs) generated by approaches like RAD sequencing. However, few empirical studies have directly compared the ability of these methods to resolve population structure. This could have implications for understanding phenotypic plasticity, as many previous studies of natural populations may have lacked the power to detect genetic differences, especially over micro-geographic scales. We therefore compared the ability of microsatellites and RAD sequencing to resolve fine-scale population structure in a commercially important benthic invertebrate by genotyping great scallops (Pecten maximus) from nine populations around Northern Ireland at 13 microsatellites and 10 539 SNPs. The shells were then subjected to morphometric and colour analysis in order to compare patterns of phenotypic and genetic variation. We found that RAD sequencing was superior at resolving population structure, yielding higher F_st values and support for two distinct genetic clusters, whereas only one cluster could be detected in a Bayesian analysis of the microsatellite dataset. Furthermore, appreciable phenotypic variation was observed in size-independent shell shape and coloration, including among localities that could not be distinguished from one another genetically, providing support for the notion that these traits are phenotypically plastic. Taken together, our results suggest that RAD sequencing is a powerful approach for studying population structure and phenotypic plasticity in natural populations.

Environmental heterogeneity and population differences in blue tits personality traits

Environmental heterogeneity can result in spatial variation in selection pressures that can produce local adaptations. The pace-of-life syndrome hypothesis predicts that habitat-specific selective pressures will favor the coevolution of personality, physiological, and life-history phenotypes. Few studies so far have compared these traits simultaneously across different ecological conditions. In this study, we compared 3 personality traits (handling aggression, exploration speed in a novel environment, and nest defense behavior) and 1 physiological trait (heart rate during manual restraint) across 3 Corsican blue tit (Cyanistes caeruleus) populations. These populations are located in contrasting habitats (evergreen vs. deciduous) and are situated in 2 different valleys 25 km apart. Birds from these populations are known to differ in life-history characteristics, with birds from the evergreen habitat displaying a slow pace-of-life, and birds from the deciduous habitat a comparatively faster pace-of-life. We expected personality to differ across populations, in line with the differences in pace-of-life documented for life-history traits. As expected, we found behavioral differences among populations. Despite considerable temporal variation, birds exhibited lower handling aggression in the evergreen populations. Exploration speed and male heart rate also differed across populations, although our results for exploration speed were more consistent with a phenotypic difference between the 2 valleys than between habitats. There were no clear differences in nest defense intensity among populations. Our study emphasizes the role of environmental heterogeneity in shaping population divergence in personality traits at a small spatial scale.

An ecological model organism flies into the genomics era

Despite the very rapid 'genomicization' of the field of Molecular Ecology in recent years, there have been relatively few annotated whole-genome assemblies of nonmodel organisms published. Instead, molecular ecologists have more frequently utilized next-generation sequencing technologies to develop genome-wide markers or to generate transcriptome data. Whole-genome assemblies are more expensive and require considerable computational resources and bioinformatic expertise. However, the availability of an annotated genome offers exciting opportunities to address fundamental questions in ecology and evolution that are difficult to address with moderate sets of markers or by transcriptome sequencing. Such questions include elucidating the roles of natural and sexual selection in shaping diversity, determining the roles of regulatory and protein-coding change in the evolution of traits, and determining the genomic architecture of sex-specific trait variation. Arguably, these questions are most tractable--and most interesting--in well-characterized species for which there is already some knowledge of natural and sexual selection, and of the traits that are most likely to link to fitness. In this issue, Mueller et al. (2016) present the assembly and annotation of the genome of the blue tit (Cyanistes caeruleus), a model ecological species. In addition, by sequencing the transcriptome of male and female blue tits, the authors identify and annotate sex-biased gene expression and conclude that noncoding RNA genes are likely to play a significant role in sex-biased expression. By making their assembly and annotation publically available and accessible via a genome browser, Mueller et al. (2016) offer exciting possibilities for further research into the genomic basis of adaptation, and investigation of the roles of natural and sexual selection, in this well-studied ecological model species.

Fine-scale genetic structure among greater sage-grouse leks in central Nevada

Background

Mating systems that reduce dispersal and lead to non-random mating might increase the potential for genetic structure to arise at fine geographic scales. Greater sage-grouse (Centrocercus urophasianus) have a lek-based mating system and exhibit high site fidelity and skewed mating ratios. We quantified population structure by analyzing variation at 27,866 single-nucleotide polymorphisms in 140 males from ten leks (within five lek complexes) occurring in a small geographic region in central Nevada.

Results

Lek complexes, and to a lesser extent individual leks, formed statistically identifiable clusters in ordination analyses, providing evidence for fine-scale geographic genetic differentiation. Lek geography predicted genetic differentiation even at a small geographic scale, which could be sharpened by strong site fidelity. Relatedness was also higher among individuals within lek complexes (and leks), suggesting that reproductive skew, where few males participate in most of the successful matings, could also potentially contribute to genetic differentiation. Models incorporating a habitat resistance surface as a proxy for potentially reduced movement due to landscape features indicated that both geographic distance and habitat suitability (i.e. preferred habitat) predicted genetic structure, with no significant effect of man-made barriers to movement (i.e. power lines and roads). Finally, we illustrate how data sets containing fewer loci (<4000) had less statistical precision and failed to detect the full degree of genetic structure.

Conclusion

Our results suggest that habitat features and lek site geography of sage-grouse shape fine scale genetic structure, and highlight how larger data sets can have increased precision and accuracy for quantifying ecologically relevant genetic structure over small geographic scales.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0702-4) contains supplementary material, which is available to authorized users.

Mediterranean blue tits as a case study of local adaptation

While the study of the origins of biological diversity across species has provided numerous examples of adaptive divergence, the realization that it can occur at microgeographic scales despite gene flow is recent, and scarcely illustrated. We review here evidence suggesting that the striking phenotypic differentiation in ecologically relevant traits exhibited by blue tits Cyanistes caeruleus in their southern range‐edge putatively reflects adaptation to the heterogeneity of the Mediterranean habitats. We first summarize the phenotypic divergence for a series of life history, morphological, behavioural, acoustic and colour ornament traits in blue tit populations of evergreen and deciduous forests. For each divergent trait, we review the evidence obtained from common garden experiments regarding a possible genetic origin of the observed phenotypic differentiation as well as evidence for heterogeneous selection. Second, we argue that most phenotypically differentiated traits display heritable variation, a fundamental requirement for evolution to occur. Third, we discuss nonrandom dispersal, selective barriers and assortative mating as processes that could reinforce local adaptation. Finally, we show how population genomics supports isolation – by – environment across landscapes. Overall, the combination of approaches converges to the conclusion that the strong phenotypic differentiation observed in Mediterranean blue tits is a fascinating case of local adaptation.