Climate change leads to differential shifts in the timing of annual cycle stages in a migratory bird

Shifts in reproductive phenology due to climate change have been well documented in many species but how, within the same species, other annual cycle stages (e.g. moult, migration) shift relative to the timing of breeding has rarely been studied. When stages shift at different rates, the interval between stages may change resulting in overlaps, and as each stage is energetically demanding, these overlaps may have negative fitness consequences. We used long‐term data of a population of European pied flycatchers (Ficedula hypoleuca) to investigate phenological shifts in three annual cycle stages: spring migration (arrival dates), breeding (egg‐laying and hatching dates) and the onset of postbreeding moult. We found different advancements in the timing of breeding compared with moult (moult advances faster) and no advancement in arrival dates. To understand these differential shifts, we explored which temperatures best explain the year‐to‐year variation in the timing of these stages, and show that they respond differently to temperature increases in the Netherlands, causing the intervals between arrival and breeding and between breeding and moult to decrease. Next, we tested the fitness consequences of these shortened intervals. We found no effect on clutch size, but the probability of a fledged chick to recruit increased with a shorter arrival‐breeding interval (earlier breeding). Finally, mark–recapture analyses did not detect an effect of shortened intervals on adult survival. Our results suggest that the advancement of breeding allows more time for fledgling development, increasing their probability to recruit. This may incur costs to other parts of the annual cycle, but, despite the shorter intervals, there was no effect on adult survival. Our results show that to fully understand the consequences of climate change, it is necessary to look carefully at different annual cycle stages, especially for organisms with complex cycles, such as migratory birds.     

Relaxed predation results in reduced phenotypic integration in a suite of dragonflies

Although changes in magnitude of single traits responding to selective agents have been studied intensively, little is known about selection shaping networks of traits and their patterns of covariation. However, this is central for our understanding of phenotypic evolution as traits are embedded in a multivariate environment with selection affecting a multitude of traits simultaneously rather than individually. Here, we investigate inter‐ and intraspecific patterns of trait integration (trait correlations) in the larval abdomen of dragonflies as a response to a change in predator selection. Species of the dragonfly genus Leucorrhinia underwent a larval habitat shift from predatory fish to predatory dragonfly‐dominated lakes with an associated relaxation in selection pressure from fish predation. Our results indicate that the habitat‐shift‐induced relaxed selection pressure caused phenotypic integration of abdominal traits to be reduced. Intraspecific findings matched patterns comparing species from both habitats with higher abdominal integration in response to predatory fish. This higher integration is probably a result of faster burst swimming speed. The abdomen holds the necessary morphological machinery to successfully evade predatory fish via burst swimming. Hence, abdominal traits have to function in a tight coordinated manner, as maladaptive variation and consequently nonoptimal burst swimming would cause increased mortality. In predatory dragonfly‐dominated lakes, no such strong link between burst swimming and mortality is present. Our findings highlight the importance of studying multivariate trait relationships as a response to selection for understanding patterns of phenotypic diversification.     

PSMC analysis of effective population sizes in molecular ecology and its application to black‐and‐white Ficedula flycatchers

Climatic fluctuations during the Quaternary period governed the demography of species and contributed to population differentiation and ultimately speciation. Studies of these past processes have previously been hindered by a lack of means and genetic data to model changes in effective population size (N_e) through time. However, based on diploid genome sequences of high quality, the recently developed pairwise sequentially Markovian coalescent (PSMC) can estimate trajectories of changes in N_e over considerable time periods. We applied this approach to resequencing data from nearly 200 genomes of four species and several populations of the Ficedula species complex of black‐and‐white flycatchers. N_e curves of Atlas, collared, pied and semicollared flycatcher converged 1–2 million years ago (Ma) at an N_e of ≈ 200 000, likely reflecting the time when all four species last shared a common ancestor. Subsequent separate N_e trajectories are consistent with lineage splitting and speciation. All species showed evidence of population growth up until 100–200 thousand years ago (kya), followed by decline and then start of a new phase of population expansion. However, timing and amplitude of changes in N_e differed among species, and for pied flycatcher, the temporal dynamics of N_e differed between Spanish birds and central/northern European populations. This cautions against extrapolation of demographic inference between lineages and calls for adequate sampling to provide representative pictures of the coalescence process in different species or populations. We also empirically evaluate criteria for proper inference of demographic histories using PSMC and arrive at recommendations of using sequencing data with a mean genome coverage of ≥18X, a per‐site filter of ≥10 reads and no more than 25% of missing data.     

Quantifying the links between land use and population growth rate in a declining farmland bird

Land use is likely to be a key driver of population dynamics of species inhabiting anthropogenic landscapes, such as farmlands. Understanding the relationships between land use and variation in population growth rates is therefore critical for the management of many farmland species. Using 24 years of data of a declining farmland bird in an integrated population model, we examined how spatiotemporal variation in land use (defined as habitats with “Short” and “Tall” ground vegetation during the breeding season) and habitat‐specific demographic parameters relates to variation in population growth taking into account individual movements between habitats. We also evaluated contributions to population growth using transient life table response experiments which gives information on contribution of past variation of parameters and real‐time elasticities which suggests future scenarios to change growth rates. LTRE analyses revealed a clear contribution of Short habitats to the annual variation in population growth rate that was mostly due to fledgling recruitment, whereas there was no evidence for a contribution of Tall habitats. Only 18% of the variation in population growth was explained by the modeled local demography, the remaining variation being explained by apparent immigration (i.e., the residual variation). We discuss potential biological and methodological reasons for high contributions of apparent immigration in open populations. In line with LTRE analysis, real‐time elasticity analysis revealed that demographic parameters linked to Short habitats had a stronger potential to influence population growth rate than those of Tall habitats. Most particularly, an increase of the proportion of Short sites occupied by Old breeders could have a distinct positive impact on population growth. High‐quality Short habitats such as grazed pastures have been declining in southern Sweden. Converting low‐quality to high‐quality habitats could therefore change the present negative population trend of this, and other species with similar habitat requirements.     

Long‐term stasis and short‐term divergence in the phenotypes of microsnails on oceanic islands

Phenotypic divergence is often unrelated to genotypic divergence. An extreme example is rapid phenotypic differentiation despite genetic similarity. Another extreme is morphological stasis despite substantial genetic divergence. These opposite patterns have been viewed as reflecting opposite properties of the lineages. In this study, phenotypic radiation accompanied by both rapid divergence and long‐term conservatism is documented in the inferred molecular phylogeny of the micro land snails Cavernacmella (Assimineidae) on the Ogasawara Islands. The populations of Cavernacmella on the Sekimon limestone outcrop of Hahajima Island showed marked divergence in shell morphology. Within this area, one lineage diversified into types with elongated turret shells, conical shells and flat disc‐like shells without substantial genetic differentiation. Additionally, a co‐occurring species with these types developed a much larger shell size. Moreover, a lineage adapted to live inside caves in this area. In contrast, populations in the other areas exhibited no morphological differences despite high genetic divergence among populations. Accordingly, the phenotypic evolution of Cavernacmella in Ogasawara is characterized by a pattern of long‐term stasis and periodic bursts of change. This pattern suggests that even lineages with phenotypic conservatism could shift to an alternative state allowing rapid phenotypic divergence.     

Phenotypic divergence among threespine stickleback that differ in nuptial coloration

By studying systems in their earliest stages of differentiation, we can learn about the evolutionary forces acting within and among populations and how those forces could contribute to reproductive isolation. Such an understanding would help us to better discern and predict how selection leads to the maintenance of multiple morphs within a species, rather than speciation. The postglacial adaptive radiation of the threespine stickleback (Gasterosteus aculeatus) is one of the best‐studied cases of evolutionary diversification and rapid, repeated speciation. Following deglaciation, marine stickleback have continually invaded freshwater habitats across the northern hemisphere and established resident populations that diverged innumerable times from their oceanic ancestors. Independent freshwater colonization events have yielded broadly parallel patterns of morphological differences in freshwater and marine stickleback. However, there is also much phenotypic diversity within and among freshwater populations. We studied a lesser‐known freshwater “species pair” found in southwest Washington, where male stickleback in numerous locations have lost the ancestral red sexual signal and instead develop black nuptial coloration. We measured phenotypic variation in a suite of traits across sites where red and black stickleback do not overlap in distribution and at one site where they historically co‐occurred. We found substantial phenotypic divergence between red and black morphs in noncolor traits including shape and lateral plating, and additionally find evidence that supports the hypothesis of sensory drive as the mechanism responsible for the evolutionary switch in color from red to black. A newly described third “mixed” morph in Connor Creek, Washington, differs in head shape and size from the red and black morphs, and we suggest that their characteristics are most consistent with hybridization between anadromous and freshwater stickleback. These results lay the foundation for future investigation of the underlying genetic basis of this phenotypic divergence as well as the evolutionary processes that may drive, maintain, or limit divergence among morphs.     

A population genomic scan in Chorthippus grasshoppers unveils previously unknown phenotypic divergence

Understanding the genetics of speciation and the processes that drive it is a central goal of evolutionary biology. Grasshoppers of the Chorthippus species group differ strongly in calling song (and corresponding female preferences) but are exceedingly similar in other characteristics such as morphology. Here, we performed a population genomic scan on three Chorthippus species (Chorthippus biguttulus, C. mollis and C. brunneus) to gain insight into the genes and processes involved in divergence and speciation in this group. Using an RNA‐seq approach, we examined functional variation between the species by calling SNPs for each of the three species pairs and using F_ST‐based approaches to identify outliers. We found approximately 1% of SNPs in each comparison to be outliers. Between 37% and 40% of these outliers were nonsynonymous SNPs (as opposed to a global level of 17%) indicating that we recovered loci under selection. Among the outliers were several genes that may be involved in song production and hearing as well as genes involved in other traits such as food preferences and metabolism. Differences in food preferences between species were confirmed with a behavioural experiment. This indicates that multiple phenotypic differences implicating multiple evolutionary processes (sexual selection and natural selection) are present between the species.     

The role of immigration and local adaptation on fine‐scale genotypic and phenotypic population divergence in a less mobile passerine

Dispersal and local patterns of adaptation play a major role on the ecological and evolutionary trajectory of natural populations. In this study, we employ a combination of genetic (25 microsatellite markers) and field‐based information (seven study years) to analyse the impact of immigration and local patterns of adaptation in two nearby (< 7 km) blue tit (Cyanistes caeruleus) populations. We used genetic assignment analyses to identify immigrant individuals and found that dispersal rate is female‐biased (72%). Data on lifetime reproductive success indicated that immigrant females produced fewer local recruits than their philopatric counterparts whereas immigrant males recruited more offspring than those that remained in their natal location. In spite of the considerably higher immigration rates of females, our results indicate that, in absolute terms, their demographic and genetic impact in the receiving populations is lower than that in immigrant males. Immigrants often brought novel alleles into the studied populations and a high proportion of them were transmitted to their recruits, indicating that the genetic impact of immigrants is not ephemeral. Although only a few kilometres apart, the two study populations were genetically differentiated and showed strong divergence in different phenotypic and life‐history traits. An almost absent inter‐population dispersal, together with the fact that both populations receive immigrants from different source populations, is probably the main cause of the observed pattern of genetic differentiation. However, phenotypic differentiation (P_ST) for all the studied traits greatly exceeded neutral genetic differentiation (F_ST), indicating that divergent natural selection is the prevailing factor determining the evolutionary trajectory of these populations. Our study highlights the importance of integrating individual‐ and population‐based approaches to obtain a comprehensive view about the role of dispersal and natural selection on structuring the genotypic and phenotypic characteristics of natural populations.     

Genetic and phenotypic population divergence on a microgeographic scale in brown trout

Salmonid populations of many rivers are rapidly declining. One possible explanation is that habitat fragmentation increases genetic drift and reduces the populations' potential to adapt to changing environmental conditions. We measured the genetic and eco-morphological diversity of brown trout (Salmo trutta) in a Swiss stream system, using multivariate statistics and Bayesian clustering. We found large genetic and phenotypic variation within only 40 km of stream length. Eighty-eight percent of all pairwise F(ST) comparisons and 50% of the population comparisons in body shape were significant. High success rates of population assignment tests confirmed the distinctiveness of populations in both genotype and phenotype. Spatial analysis revealed that divergence increased with waterway distance, the number of weirs, and stretches of poor habitat between sampling locations, but effects of isolation-by-distance and habitat fragmentation could not be fully disentangled. Stocking intensity varied between streams but did not appear to erode genetic diversity within populations. A lack of association between phenotypic and genetic divergence points to a role of local adaptation or phenotypically plastic responses to habitat heterogeneity. Indeed, body shape could be largely explained by topographic stream slope, and variation in overall phenotype matched the flow regimes of the respective habitats.     

Genome divergence and diversification within a geographic mosaic of coevolution

Despite substantial interest in coevolution's role in diversification, examples of coevolution contributing to speciation have been elusive. Here, we build upon past studies that have shown both coevolution between South Hills crossbills and lodgepole pine (Pinus contorta), and high levels of reproductive isolation between South Hills crossbills and other ecotypes in the North American red crossbill (Loxia curvirostra) complex. We used genotyping by sequencing to generate population genomic data and applied phylogenetic and population genetic analyses to characterize the genetic structure within and among nine of the ecotypes. Although genome‐wide divergence was slight between ecotypes (F_ST = 0.011–0.035), we found evidence of relative genetic differentiation (as measured by F_ST) between and genetic cohesiveness within many of them. As expected for nomadic and opportunistic breeders, we detected no evidence of isolation by distance. The one sedentary ecotype, the South Hills crossbill, was genetically most distinct because of elevated divergence at a small number of loci rather than pronounced overall genome‐wide divergence. These findings suggest that mechanisms related to recent local coevolution between South Hills crossbills and lodgepole pine (e.g. strong resource‐based density dependence limiting gene flow) have been associated with genome divergence in the face of gene flow. Our results further characterize a striking example of coevolution driving speciation within perhaps as little as 6000 years.     

Estimation of linkage disequilibrium and interspecific gene flow in Ficedula flycatchers by a newly developed 50k single‐nucleotide polymorphism array

With the access to draft genome sequence assemblies and whole‐genome resequencing data from population samples, molecular ecology studies will be able to take truly genome‐wide approaches. This now applies to an avian model system in ecological and evolutionary research: Old World flycatchers of the genus Ficedula, for which we recently obtained a 1.1 Gb collared flycatcher genome assembly and identified 13 million single‐nucleotide polymorphism (SNP)s in population resequencing of this species and its sister species, pied flycatcher. Here, we developed a custom 50K Illumina iSelect flycatcher SNP array with markers covering 30 autosomes and the Z chromosome. Using a number of selection criteria for inclusion in the array, both genotyping success rate and polymorphism information content (mean marker heterozygosity = 0.41) were high. We used the array to assess linkage disequilibrium (LD) and hybridization in flycatchers. Linkage disequilibrium declined quickly to the background level at an average distance of 17 kb, but the extent of LD varied markedly within the genome and was more than 10‐fold higher in ‘genomic islands’ of differentiation than in the rest of the genome. Genetic ancestry analysis identified 33 F₁ hybrids but no later‐generation hybrids from sympatric populations of collared flycatchers and pied flycatchers, contradicting earlier reports of backcrosses identified from much fewer number of markers. With an estimated divergence time as recently as <1 Ma, this suggests strong selection against F₁ hybrids and unusually rapid evolution of reproductive incompatibility in an avian system.     

Phylogeography and population genetics of pine butterflies: Sky islands increase genetic divergence

The sky islands of southeastern Arizona (AZ) mark a major transition zone between tropical and temperate biota and are considered a neglected biodiversity hotspot. Dispersal ability and host plant specificity are thought to impact the history and diversity of insect populations across the sky islands. We aimed to investigate the population structure and phylogeography of two pine‐feeding pierid butterflies, the pine white (Neophasia menapia) and the Mexican pine white (Neophasia terlooii), restricted to these “islands” at this transition zone. Given their dependence on pines as the larval hosts, we hypothesized that habitat connectivity affects population structure and is at least in part responsible for their allopatry. We sampled DNA from freshly collected butterflies from 17 sites in the sky islands and adjacent high‐elevation habitats and sequenced these samples using ddRADSeq. Up to 15,399 SNPs were discovered and analyzed in population genetic and phylogenetic contexts with Stacks and pyRAD pipelines. Low genetic differentiation in N. menapia suggests that it is panmictic. Conversely, there is strong evidence for population structure within N. terlooii. Each sky island likely contains a population of N. terlooii, and clustering is hierarchical, with populations on proximal mountains being more related to each other. The N. menapia habitat, which is largely contiguous, facilitates panmixia, while the N. terlooii habitat, restricted to the higher elevations on each sky island, creates distinct population structure. Phylogenetic results corroborate those from population genetic analyses. The historical climate‐driven fluxes in forest habitat connectivity have implications for understanding the biodiversity of fragmented habitats.     

Qualitative population divergence in proximate determination of a sexually selected trait in the collared flycatcher

We examined proximate determination of sexually selected forehead patch size in a Central‐European population of Ficedula albicollis, the collared flycatcher, using a 9‐year database, and compared our results with those obtained in other populations of the same and the sister species. Between‐individual variation of forehead patch size was large, its repeatability larger than, and heritability similar to the Swedish population. Unlike in the other populations, the trait proved unaffected by body condition, and only very slightly influenced by age. There was no relationship between forehead patch size and breeding lifespan, and a marginal negative association with survivorship in adult males. Our results suggest that additive genetic variance of the trait in this population is large, but genes act independently of body condition, and there is no viability indicator value of the trait. This is the first report of a qualitative intraspecific difference in proximate determination of a sexually selected trait.     

The extent and genetic basis of phenotypic divergence in life history traits in Mimulus guttatus

Differential natural selection acting on populations in contrasting environments often results in adaptive divergence in multivariate phenotypes. Multivariate trait divergence across populations could be caused by selection on pleiotropic alleles or through many independent loci with trait‐specific effects. Here, we assess patterns of association between a suite of traits contributing to life history divergence in the common monkey flower, Mimulus guttatus, and examine the genetic architecture underlying these correlations. A common garden survey of 74 populations representing annual and perennial strategies from across the native range revealed strong correlations between vegetative and reproductive traits. To determine whether these multitrait patterns arise from pleiotropic or independent loci, we mapped QTLs using an approach combining high‐throughput sequencing with bulk segregant analysis on a cross between populations with divergent life histories. We find extensive pleiotropy for QTLs related to flowering time and stolon production, a key feature of the perennial strategy. Candidate genes related to axillary meristem development colocalize with the QTLs in a manner consistent with either pleiotropic or independent QTL effects. Further, these results are analogous to previous work showing pleiotropy‐mediated genetic correlations within a single population of M. guttatus experiencing heterogeneous selection. Our findings of strong multivariate trait associations and pleiotropic QTLs suggest that patterns of genetic variation may determine the trajectory of adaptive divergence.     

Strong dispersal in a parasitoid wasp overwhelms habitat fragmentation and host population dynamics

The population dynamics of a parasite depend on species traits, host dynamics and the environment. Those dynamics are reflected in the genetic structure of the population. Habitat fragmentation has a greater impact on parasites than on their hosts because resource distribution is increasingly fragmented for species at higher trophic levels. This could lead to either more or less genetic structure than the host, depending on the relative dispersal rates of species. We examined the spatial genetic structure of the parasitoid wasp Hyposoter horticola, and how it was influenced by dispersal, host population dynamics and habitat fragmentation. The host, the Glanville fritillary butterfly, lives as a metapopulation in a fragmented landscape in the Åland Islands, Finland. We collected wasps throughout the 50 by 70 km archipelago and determined the genetic diversity, spatial population structure and genetic differentiation using 14 neutral DNA microsatellite loci. We compared the genetic structure of the wasp with that of the host butterfly using published genetic data collected over the shared landscape. Using maternity assignment, we also identified full‐siblings among the sampled parasitoids to estimate the dispersal range of individual females. We found that because the parasitoid is dispersive, it has low genetic structure, is not very sensitive to habitat fragmentation and has less spatial genetic structure than its butterfly host. The wasp is sensitive to regional rather than local host dynamics, and there is a geographic mosaic landscape for antagonistic co‐evolution of host resistance and parasite virulence.     

Gene flow and population structure of a solitary top carnivore in a human‐dominated landscape

While African leopard populations are considered to be continuous as demonstrated by their high genetic variation, the southernmost leopard population exists in the Eastern and Western Cape, South Africa, where anthropogenic activities may be affecting this population's structure. Little is known about the elusive, last free‐roaming top predator in the region and this study is the first to report on leopard population structuring using nuclear DNA. By analyzing 14 microsatellite markers from 40 leopard tissue samples, we aimed to understand the populations' structure, genetic distance, and gene flow (Nm). Our results, based on spatially explicit analysis with Bayesian methods, indicate that leopards in the region exist in a fragmented population structure with lower than expected genetic diversity. Three population groups were identified, between which low to moderate levels of gene flow were observed (Nm 0.5 to 3.6). One subpopulation exhibited low genetic differentiation, suggesting a continuous population structure, while the remaining two appear to be less connected, with low emigration and immigration between these populations. Therefore, genetic barriers are present between the subpopulations, and while leopards in the study region may function as a metapopulation, anthropogenic activities threaten to decrease habitat and movement further. Our results indicate that the leopard population may become isolated within a few generations and suggest that management actions should aim to increase habitat connectivity and reduce human–carnivore conflict. Understanding genetic diversity and connectivity of populations has important conservation implications that can highlight management of priority populations to reverse the effects of human‐caused extinctions.     

Change in sexual signalling traits outruns morphological divergence across an ecological gradient in the post‐glacial radiation of the songbird genus Junco

The relative roles of natural and sexual selection in promoting evolutionary lineage divergence remains controversial and difficult to assess in natural systems. Local adaptation through natural selection is known to play a central role in promoting evolutionary divergence, yet secondary sexual traits can vary widely among species in recent radiations, suggesting that sexual selection may also be important in the early stages of speciation. Here, we compare rates of divergence in ecologically relevant traits (morphology) and sexually selected signalling traits (coloration) relative to neutral structure in genome‐wide molecular markers and examine patterns of variation in sexual dichromatism to explore the roles of natural and sexual selection in the diversification of the songbird genus Junco (Aves: Passerellidae). Juncos include divergent lineages in Central America and several dark‐eyed junco (J. hyemalis) lineages that diversified recently as the group recolonized North America following the last glacial maximum (ca. 18,000 years ago). We found an accelerated rate of divergence in sexually selected characters relative to ecologically relevant traits. Moreover, sexual dichromatism measurements suggested a positive relationship between the degree of colour divergence and the strength of sexual selection when controlling for neutral genetic distance. We also found a positive correlation between dichromatism and latitude, which coincides with the geographic axis of decreasing lineage age in juncos but also with a steep ecological gradient. Finally, we found significant associations between genome‐wide variants linked to functional genes and proxies of both sexual and natural selection. These results suggest that the joint effects of sexual and ecological selection have played a prominent role in the junco radiation.     

Sampling related individuals within ponds biases estimates of population structure in a pond‐breeding amphibian

Effective conservation and management of pond‐breeding amphibians depends on the accurate estimation of population structure, demographic parameters, and the influence of landscape features on breeding‐site connectivity. Population‐level studies of pond‐breeding amphibians typically sample larval life stages because they are easily captured and can be sampled nondestructively. These studies often identify high levels of relatedness between individuals from the same pond, which can be exacerbated by sampling the larval stage. Yet, the effect of these related individuals on population genetic studies using genomic data is not yet fully understood. Here, we assess the effect of within‐pond relatedness on population and landscape genetic analyses by focusing on the barred tiger salamanders (Ambystoma mavortium) from the Nebraska Sandhills. Utilizing genome‐wide SNPs generated using a double‐digest RADseq approach, we conducted standard population and landscape genetic analyses using datasets with and without siblings. We found that reduced sample sizes influenced parameter estimates more than the inclusion of siblings, but that within‐pond relatedness led to the inference of spurious population structure when analyses depended on allele frequencies. Our landscape genetic analyses also supported different models across datasets depending on the spatial resolution analyzed. We recommend that future studies not only test for relatedness among larval samples but also remove siblings before conducting population or landscape genetic analyses. We also recommend alternative sampling strategies to reduce sampling siblings before sequencing takes place. Biases introduced by unknowingly including siblings can have significant implications for population and landscape genetic analyses, and in turn, for species conservation strategies and outcomes.