Genomic landscapes of divergence among island bird populations: Evidence of parallel adaptation but at different loci?

When populations colonise new environments, they may be exposed to novel selection pressures but also suffer from extensive genetic drift due to founder effects, small population sizes and limited interpopulation gene flow. Genomic approaches enable us to study how these factors drive divergence, and disentangle neutral effects from differentiation at specific loci due to selection. Here, we investigate patterns of genetic diversity and divergence using whole-genome resequencing (>22× coverage) in Berthelot's pipit (Anthus berthelotii), a passerine endemic to the islands of three north Atlantic archipelagos. Strong environmental gradients, including in pathogen pressure, across populations in the species range, make it an excellent system in which to explore traits important in adaptation and/or incipient speciation. First, we quantify how genomic divergence accumulates across the speciation continuum, that is, among Berthelot's pipit populations, between sub species across archipelagos, and between Berthelot's pipit and its mainland ancestor, the tawny pipit (Anthus campestris). Across these colonisation timeframes (2.1 million-ca. 8000 years ago), we identify highly differentiated loci within genomic islands of divergence and conclude that the observed distributions align with expectations for non-neutral divergence. Characteristic signatures of selection are identified in loci associated with craniofacial/bone and eye development, metabolism and immune response between population comparisons. Interestingly, we find limited evidence for repeated divergence of the same loci across the colonisation range but do identify different loci putatively associated with the same biological traits in different populations, likely due to parallel adaptation. Incipient speciation across these island populations, in which founder effects and selective pressures are strong, may therefore be repeatedly associated with morphology, metabolism and immune defence.

Genotype-environment associations reveal genes potentially linked to avian malaria infection in populations of an endemic island bird.

Patterns of pathogen prevalence are, at least partially, the result of coevolutionary host-pathogen interactions. Thus, exploring the distribution of host genetic variation in relation to infection by a pathogen within and across populations can provide important insights into mechanisms of host defence and adaptation. Here, we use a landscape genomics approach (Bayenv) in conjunction with genome-wide data (ddRADseq) to test for associations between avian malaria (Plasmodium) prevalence and host genetic variation across 13 populations of the island endemic Berthelot's pipit (Anthus berthelotii). Considerable and consistent spatial heterogeneity in malaria prevalence was observed among populations over a period of 15 years. The prevalence of malaria infection was also strongly positively correlated with pox (Avipoxvirus) prevalence. Multiple host loci showed significant associations with malaria prevalence after controlling for genome-wide neutral genetic structure. These sites were located near to or within genes linked to metabolism, stress response, transcriptional regulation, complement activity and the inflammatory response, many previously implicated in vertebrate responses to malarial infection. Our findings identify diverse genes - not just limited to the immune system - that may be involved in host protection against malaria and suggest that spatially variable pathogen pressure may be an important evolutionary driver of genetic divergence among wild animal populations, such as Berthelot's pipit. Furthermore, our data indicate that spatio-temporal variation in multiple different pathogens (e.g. malaria and pox in this case) may have to be studied together to develop a more holistic understanding of host pathogen-mediated evolution.

Life history correlations and trade-offs resulting from selection for dispersal in Tribolium castaneum.

Dispersal is an important facet of the life history of many organisms and is, therefore, subject to selective pressure but does not evolve in isolation. Across nature, there are examples of dispersal syndromes and life history strategies in which suites of traits coevolve and covary with dispersal in combinations that serve to maximize fitness in a given ecological context. The red rust flour beetle, Tribolium castaneum, is a model organism and globally significant post-harvest pest that relies on dispersal to reach new patches of ephemeral habitat. Dispersal behaviour in Tribolium has a strong genetic basis. However, a robust understanding of the relationship between dispersal and other life-history components, which could elucidate evolutionary processes and allow pest managers to control their spread and reduce the impact of infestation, is currently lacking. Here, we use highly replicated lines of T. castaneum previously artificially selected for divergent small-scale dispersal propensity to robustly test several important life history components: reproductive strategy, development time, and longevity. As predicted, we find that a suite of important changes as a result of our selection on dispersal: high dispersal propensity is associated with a lower number of longer mating attempts by males, lower investment in early life reproduction by females, slower development of later-laid offspring, and longer female life span. These findings indicate that correlated intraspecific variation in dispersal and related traits may represent alternative life history strategies in T. castaneum. We therefore suggest that pest management efforts to mitigate the species' agro-economic impact should consider the eco-evolutionary dynamics within multiple life histories. The benefits of doing so could be felt both through improved targeting of efforts to reduce spread and also in forecasting how the selection pressures applied through pest management are likely to affect pest evolution.

Runs of homozygosity reveal past bottlenecks and contemporary inbreeding across diverging populations of an island-colonizing bird.

Genomes retain evidence of the demographic history and evolutionary forces that have shaped populations and drive speciation. Across island systems, contemporary patterns of genetic diversity reflect population demography, including colonization events, bottlenecks, gene flow and genetic drift. Here, we investigate genome-wide diversity and the distribution of runs of homozygosity (ROH) using whole-genome resequencing of individuals (>22× coverage) from six populations across three archipelagos of Berthelot's pipit (Anthus berthelotii)-a passerine that has recently undergone island speciation. We show the most dramatic reduction in diversity occurs between the mainland sister species (the tawny pipit) and Berthelot's pipit and is lowest in the populations that have experienced sequential bottlenecks (i.e., the Madeiran and Selvagens populations). Pairwise sequential Markovian coalescent (PSMC) analyses estimated that Berthelot's pipit diverged from its sister species ~2 million years ago, with the Madeiran archipelago founded 50,000 years ago, and the Selvagens colonized 8000 years ago. We identify many long ROH (>1 Mb) in these most recently colonized populations. Population expansion within the last 100 years may have eroded long ROH in the Madeiran archipelago, resulting in a prevalence of short ROH (<1 Mb). However, the extensive long and short ROH detected in the Selvagens suggest strong recent inbreeding and bottleneck effects, with as much as 38% of the autosomes consisting of ROH >250 kb. These findings highlight the importance of demographic history, as well as selection and genetic drift, in shaping contemporary patterns of genomic diversity across diverging populations.

Genetic architecture of dispersal behaviour in the post-harvest pest and model organism Tribolium castaneum.

Dispersal behaviour is an important aspect of the life-history of animals. However, the genetic architecture of dispersal-related traits is often obscure or unknown, even in well studied species. Tribolium castaneum is a globally significant post-harvest pest and established model organism, yet studies of its dispersal have shown ambiguous results and the genetic basis of this behaviour remains unresolved. We combine experimental evolution and agent-based modelling to investigate the number of loci underlying dispersal in T. castaneum, and whether the trait is sex-linked. Our findings demonstrate rapid evolution of dispersal behaviour under selection. We find no evidence of sex-biases in the dispersal behaviour of the offspring of crosses, supporting an autosomal genetic basis of the trait. Moreover, simulated data approximates experimental data under simulated scenarios where the dispersal trait is controlled by one or few loci, but not many loci. Levels of dispersal in experimentally inbred lines, compared with simulations, indicate that a single locus model is not well supported. Taken together, these lines of evidence support an oligogenic architecture underlying dispersal in Tribolium castaneum. These results have implications for applied pest management and for our understanding of the evolution of dispersal in the coleoptera, the world's most species-rich order.

Causes and consequences of telomere lengthening in a wild vertebrate population.

Telomeres have been advocated to be important markers of biological age in evolutionary and ecological studies. Telomeres usually shorten with age and shortening is frequently associated with environmental stressors and increased subsequent mortality. Telomere lengthening - an apparent increase in telomere length between repeated samples from the same individual - also occurs. However, the exact circumstances, and consequences, of telomere lengthening are poorly understood. Using longitudinal data from the Seychelles warbler (Acrocephalus sechellensis), we tested whether telomere lengthening - which occurs in adults of this species - is associated with specific stressors (reproductive effort, food availability, malarial infection and cooperative breeding) and predicts subsequent survival. In females, telomere shortening was observed under greater stress (i.e., low food availability, malaria infection), while telomere lengthening was observed in females experiencing lower stress (i.e., high food availability, assisted by helpers, without malaria). The telomere dynamics of males were not associated with the key stressors tested. These results indicate that, at least for females, telomere lengthening occurs in circumstances more conducive to self-maintenance. Importantly, both females and males with lengthened telomeres had improved subsequent survival relative to individuals that displayed unchanged, or shortened, telomeres - indicating that telomere lengthening is associated with individual fitness. These results demonstrate that telomere dynamics are bidirectionally responsive to the level of stress that an individual faces, but may poorly reflect the accumulation of stress over an individuals lifetime.

Genomic associations with poxvirus across divergent island populations in Berthelot's pipit.

Understanding the mechanisms and genes that enable animal populations to adapt to pathogens is important from an evolutionary, health and conservation perspective. Berthelot's pipit (Anthus berthelotii) experiences extensive and consistent spatial heterogeneity in avian pox infection pressure across its range of island populations, thus providing an excellent system with which to examine how pathogen-mediated selection drives spatial variation in immunogenetic diversity. Here, we test for evidence of genetic variation associated with avian pox at both an individual and population-level. At the individual level, we find no evidence that variation in MHC class I and TLR4 (both known to be important in recognising viral infection) was associated with pox infection within two separate populations. However, using genotype-environment association (Bayenv) in conjunction with genome-wide (ddRAD-seq) data, we detected strong associations between population-level avian pox prevalence and allele frequencies of single nucleotide polymorphisms (SNPs) at a number of sites across the genome. These sites were located within genes involved in cellular stress signalling and immune responses, many of which have previously been associated with responses to viral infection in humans and other animals. Consequently, our analyses indicate that pathogen-mediated selection may play a role in shaping genomic variation among relatively recently colonised island bird populations and highlight the utility of genotype-environment associations for identifying candidate genes potentially involved in host-pathogen interactions.

Telomere heritability and parental age at conception effects in a wild avian population.

Individual variation in telomere length is predictive of health and mortality risk across a range of species. However, the relative influence of environmental and genetic variation on individual telomere length in wild populations remains poorly understood. Heritability of telomere length has primarily been calculated using parent-offspring regression which can be confounded by shared environments. To control for confounding variables, quantitative genetic "animal models" can be used, but few studies have applied animal models in wild populations. Furthermore, parental age at conception may also influence offspring telomere length, but most studies have been cross-sectional. We investigated within- and between-parental age at conception effects and heritability of telomere length in the Seychelles warbler using measures from birds caught over 20 years and a multigenerational pedigree. We found a weak negative within-paternal age at conception effect (as fathers aged, their offspring had shorter telomeres) and a weak positive between-maternal age at conception effect (females that survived to older ages had offspring with longer telomeres). Animal models provided evidence that heritability and evolvability of telomere length were low in this population, and that variation in telomere length was not driven by early-life effects of hatch period or parental identities. Quantitative polymerase chain reaction plate had a large influence on telomere length variation and not accounting for it in the models would have underestimated heritability. Our study illustrates the need to include and account for technical variation in order to accurately estimate heritability, as well as other environmental effects, on telomere length in natural populations.

Tribolium beetles as a model system in evolution and ecology.

Flour beetles of the genus Tribolium have been utilised as informative study systems for over a century and contributed to major advances across many fields. This review serves to highlight the significant historical contribution that Tribolium study systems have made to the fields of ecology and evolution, and to promote their use as contemporary research models. We review the broad range of studies employing Tribolium to make significant advances in ecology and evolution. We show that research using Tribolium beetles has contributed a substantial amount to evolutionary and ecological understanding, especially in the fields of population dynamics, reproduction and sexual selection, population and quantitative genetics, and behaviour, physiology and life history. We propose a number of future research opportunities using Tribolium, with particular focus on how their amenability to forward and reverse genetic manipulation may provide a valuable complement to other insect models.

Sexual selection protects against extinction.

Reproduction through sex carries substantial costs, mainly because only half of sexual adults produce offspring. It has been theorized that these costs could be countered if sex allows sexual selection to clear the universal fitness constraint of mutation load. Under sexual selection, competition between (usually) males and mate choice by (usually) females create important intraspecific filters for reproductive success, so that only a subset of males gains paternity. If reproductive success under sexual selection is dependent on individual condition, which is contingent to mutation load, then sexually selected filtering through 'genic capture' could offset the costs of sex because it provides genetic benefits to populations. Here we test this theory experimentally by comparing whether populations with histories of strong versus weak sexual selection purge mutation load and resist extinction differently. After evolving replicate populations of the flour beetle Tribolium castaneum for 6 to 7 years under conditions that differed solely in the strengths of sexual selection, we revealed mutation load using inbreeding. Lineages from populations that had previously experienced strong sexual selection were resilient to extinction and maintained fitness under inbreeding, with some families continuing to survive after 20 generations of sib × sib mating. By contrast, lineages derived from populations that experienced weak or non-existent sexual selection showed rapid fitness declines under inbreeding, and all were extinct after generation 10. Multiple mutations across the genome with individually small effects can be difficult to clear, yet sum to a significant fitness load; our findings reveal that sexual selection reduces this load, improving population viability in the face of genetic stress.

Kinship and familiarity mitigate costs of social conflict between Seychelles warbler neighbors.

Because virtually all organisms compete with others in their social environment, mechanisms that reduce conflict between interacting individuals are crucial for the evolution of stable families, groups, and societies. Here, we tested whether costs of social conflict over territorial space between Seychelles warblers (Acrocephalus sechellensis) are mitigated by kin-selected (genetic relatedness) or mutualistic (social familiarity) mechanisms. By measuring longitudinal changes in individuals' body mass and telomere length, we demonstrated that the fitness costs of territoriality are driven by a complex interplay between relatedness, familiarity, local density, and sex. Physical fights were less common at territory boundaries shared between related or familiar males. In line with this, male territory owners gained mass when living next to related or familiar males and also showed less telomere attrition when living next to male kin. Importantly, these relationships were strongest in high-density areas of the population. Males also had more rapid telomere attrition when living next to unfamiliar male neighbors, but mainly when relatedness to those neighbors was also low. In contrast, neither kinship nor familiarity was linked to body mass or telomere loss in female territory owners. Our results indicate that resolving conflict over territorial space through kin-selected or mutualistic pathways can reduce both immediate energetic costs and permanent somatic damage, thus providing an important mechanism to explain fine-scale population structure and cooperation between different social units across a broad range of taxa.

Differential female sociality is linked with the fine-scale structure of sexual interactions in replicate groups of red junglefowl, Gallus gallus.

Recent work indicates that social structure has extensive implications for patterns of sexual selection and sexual conflict. However, little is known about the individual variation in social behaviours linking social structure to sexual interactions. Here, we use network analysis of replicate polygynandrous groups of red junglefowl (Gallus gallus) to show that the association between social structure and sexual interactions is underpinned by differential female sociality. Sexual dynamics are largely explained by a core group of highly social, younger females, which are more fecund and more polyandrous, and thus associated with more intense postcopulatory competition for males. By contrast, less fecund females from older cohorts, which tend to be socially dominant, avoid male sexual attention by clustering together and perching on branches, and preferentially reproduce with dominant males by more exclusively associating and mating with them. Collectively, these results indicate that individual females occupy subtly different social niches and demonstrate that female sociality can be an important factor underpinning the landscape of intrasexual competition and the emergent structure of animal societies.

Spatio-temporal variation in lifelong telomere dynamics in a long-term ecological study.

Understanding individual-level variation in response to the environment is fundamental to understanding life-history evolution and population dynamics. Telomeres, the protective caps at the ends of chromosomes, shorten in response to oxidative stress, and telomere shortening is correlated with reduced survival and life span. Investigating telomere dynamics may help us quantify individual variation in the costs experienced from social and ecological factors, and enhance our understanding of the dynamics of natural populations. Here, we study spatio-temporal variation in lifelong telomere dynamics in the Seychelles warbler, Acrocephalus sechellensis. We combine long-term life history and ecological data with a large longitudinal dataset of mean telomere lengths, consisting of 1,808 samples from 22 cohorts born between 1993 and 2014. We provide a detailed analysis of how telomere dynamics vary over individual life spans and cohorts, and with spatio-temporal variation in the social and ecological environment. We found that telomere length decreases with cross-sectional and longitudinal measures of age, and most rapidly very early in life. However, both cross-sectional and longitudinal data suggested that against this overall pattern of shortening, bouts of telomere length increase occur in some individuals. Using a large number of repeated measurements we show statistically that these increases are unlikely to be explained solely by qPCR measurement error. Telomere length varied markedly among cohorts. Telomere length was positively associated with temporal variation in island-wide insect abundance-a key resource for the insectivorous Seychelles warbler-suggesting that the costs associated with living in harsher environments can be studied by investigating telomere dynamics. We also found evidence for sex-specific relationships between telomeres and tarsus length, potentially reflecting differential costs of growth. Our long-term data show that in a natural population, telomere dynamics vary in a complex manner over individual life spans, and across space and time. Variance in telomere dynamics among individuals is the product of a wide array of genetic, parental and environmental factors. Explaining this variation more fully will require the integration of comprehensive long-term ecological and genetic data from multiple populations and species.

Telomere length reveals cumulative individual and transgenerational inbreeding effects in a passerine bird.

Inbreeding results in more homozygous offspring that should suffer reduced fitness, but it can be difficult to quantify these costs for several reasons. First, inbreeding depression may vary with ecological or physiological stress and only be detectable over long time periods. Second, parental homozygosity may indirectly affect offspring fitness, thus confounding analyses that consider offspring homozygosity alone. Finally, measurement of inbreeding coefficients, survival and reproductive success may often be too crude to detect inbreeding costs in wild populations. Telomere length provides a more precise measure of somatic costs, predicts survival in many species and should reflect differences in somatic condition that result from varying ability to cope with environmental stressors. We studied relative telomere length in a wild population of Seychelles warblers (Acrocephalus sechellensis) to assess the lifelong relationship between individual homozygosity, which reflects genome-wide inbreeding in this species, and telomere length. In juveniles, individual homozygosity was negatively associated with telomere length in poor seasons. In adults, individual homozygosity was consistently negatively related to telomere length, suggesting the accumulation of inbreeding depression during life. Maternal homozygosity also negatively predicted offspring telomere length. Our results show that somatic inbreeding costs are environmentally dependent at certain life stages but may accumulate throughout life.

Landscape-scale variation in an anthropogenic factor shapes immune gene variation within a wild population.

Understanding the spatial scale at which selection acts upon adaptive genetic variation in natural populations is fundamental to our understanding of evolutionary ecology, and has important ramifications for conservation. The environmental factors to which individuals of a population are exposed can vary at fine spatial scales, potentially generating localized patterns of adaptation. Here, we compared patterns of neutral and major histocompatibility complex (MHC) variation within an island population of Berthelot's pipit (Anthus berthelotii) to assess whether landscape-level differences in pathogen-mediated selection generate fine-scale spatial structuring in these immune genes. Specifically, we tested for spatial associations between the distribution of avian malaria, and the factors previously shown to influence that distribution, and MHC variation within resident individuals. Although we found no overall genetic structure across the population for either neutral or MHC loci, we did find localized associations between environmental factors and MHC variation. One MHC class I allele (ANBE48) was directly associated with malaria infection risk, while the presence of the ANBE48 and ANBE38 alleles within individuals correlated (positively and negatively, respectively) with distance to the nearest poultry farm, an anthropogenic factor previously shown to be an important determinant of disease distribution in the study population. Our findings highlight the importance of considering small spatial scales when studying the patterns and processes involved in evolution at adaptive loci.

454 screening of individual MHC variation in an endemic island passerine.

Genes of the major histocompatibility complex (MHC) code for receptors that are central to the adaptive immune response of vertebrates. These genes are therefore important genetic markers with which to study adaptive genetic variation in the wild. Next-generation sequencing (NGS) has increasingly been used in the last decade to genotype the MHC. However, NGS methods are highly prone to sequencing errors, and although several methodologies have been proposed to deal with this, until recently there have been no standard guidelines for the validation of putative MHC alleles. In this study, we used the 454 NGS platform to screen MHC class I exon 3 variation in a population of the island endemic Berthelot's pipit (Anthus berthelotii). We were able to characterise MHC genotypes across 309 individuals with high levels of repeatability. We were also able to determine alleles that had low amplification efficiencies, whose identification within individuals may thus be less reliable. At the population level we found lower levels of MHC diversity in Berthelot's pipit than in its widespread continental sister species the tawny pipit (Anthus campestris), and observed trans-species polymorphism. Using the sequence data, we identified signatures of gene conversion and evidence of maintenance of functionally divergent alleles in Berthelot's pipit. We also detected positive selection at 10 codons. The present study therefore shows that we have an efficient method for screening individual MHC variation across large datasets in Berthelot's pipit, and provides data that can be used in future studies investigating spatio-temporal patterns and scales of selection on the MHC.

Drift, not selection, shapes toll-like receptor variation among oceanic island populations.

Understanding the relative role of different evolutionary forces in shaping the level and distribution of functional genetic diversity among natural populations is a key issue in evolutionary and conservation biology. To do so accurately genetic data must be analysed in conjunction with an unambiguous understanding of the historical processes that have acted upon the populations. Here, we focused on diversity at toll-like receptor (TLR) loci, which play a key role in the vertebrate innate immune system and, therefore, are expected to be under pathogen-mediated selection. We assessed TLR variation within and among 13 island populations (grouped into three archipelagos) of Berthelot's pipit, Anthus berthelotii, for which detailed population history has previously been ascertained. We also compared the variation observed with that found in its widespread sister species, the tawny pipit, Anthus campestris. We found strong evidence for positive selection at specific codons in TLR1LA, TLR3 and TLR4. Despite this, we found that at the allele frequency level, demographic history has played the major role in shaping patterns of TLR variation in Berthelot's pipit. Levels of diversity and differentiation within and across archipelagos at all TLR loci corresponded very closely with neutral microsatellite variation and with the severity of the bottlenecks that occurred during colonization. Our study shows that despite the importance of TLRs in combating pathogens, demography can be the main driver of immune gene variation within and across populations, resulting in patterns of functional variation that can persist over evolutionary timescales.

Genetic and phenotypic divergence in an island bird: isolation by distance, by colonization or by adaptation?

Discerning the relative roles of adaptive and nonadaptive processes in generating differences among populations and species, as well as how these processes interact, is a fundamental aim in biology. Both genetic and phenotypic divergence across populations can be the product of limited dispersal and gradual genetic drift across populations (isolation by distance), of colonization history and founder effects (isolation by colonization) or of adaptation to different environments preventing migration between populations (isolation by adaptation). Here, we attempt to differentiate between these processes using island populations of Berthelot's pipit (Anthus berthelotii), a passerine bird endemic to three Atlantic archipelagos. Using microsatellite markers and approximate Bayesian computation, we reveal that the northward colonization of this species ca. 8500 years ago resulted in genetic bottlenecks in the colonized archipelagos. We then show that high levels of genetic structure exist across archipelagos and that these are consistent with a pattern of isolation by colonization, but not with isolation by distance or adaptation. Finally, we show that substantial morphological divergence also exists and that this is strongly concordant with patterns of genetic structure and bottleneck history, but not with environmental differences or geographic distance. Overall, our data suggest that founder effects are responsible for both genetic and phenotypic changes across archipelagos. Our findings provide a rare example of how founder effects can persist over evolutionary timescales and suggest that they may play an important role in the early stages of speciation.

The impact of translocations on neutral and functional genetic diversity within and among populations of the Seychelles warbler.

Translocations are an increasingly common tool in conservation. The maintenance of genetic diversity through translocation is critical for both the short- and long-term persistence of populations and species. However, the relative spatio-temporal impacts of translocations on neutral and functional genetic diversity, and how this affects genetic structure among the conserved populations overall, have received little investigation. We compared the impact of translocating different numbers of founders on both microsatellite and major histocompatibility complex (MHC) class I diversity over a 23-year period in the Seychelles warbler (Acrocephalus sechellensis). We found low and stable microsatellite and MHC diversity in the source population and evidence for only a limited loss of either type of diversity in the four new populations. However, we found evidence of significant, but low to moderate, genetic differentiation between populations, with those populations established with fewer founders clustering separately. Stochastic genetic capture (as opposed to subsequent drift) was the main determinant of translocated population diversity. Furthermore, a strong correlation between microsatellite and MHC differentiation suggested that neutral processes outweighed selection in shaping MHC diversity in the new populations. These data provide important insights into how to optimize the use of translocation as a conservation tool.