Genomic variation, population history and within-archipelago adaptation between island bird populations

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.

Demographic responses of oceanic island birds to local and regional ecological disruptions revealed by whole-genome sequencing

Disentangling the effects of ecological disruptions operating at different spatial and temporal scales in shaping past species' demography is particularly important in the current context of rapid environmental changes driven by both local and regional factors. We argue that volcanic oceanic islands provide useful settings to study the influence of past ecological disruptions operating at local and regional scales on population demographic histories. We investigate potential drivers of past population dynamics for three closely related species of passerine birds from two volcanic oceanic islands, Reunion and Mauritius (Mascarene archipelago), with distinct volcanic history. Using ABC and PSMC inferences from complete genomes, we reconstructed the demographic history of the Reunion Grey White-eye (Zosterops borbonicus (Pennant, 1781)), the Reunion Olive White-eye (Z. olivaceus (Linnaeus, 1766)) and the Mauritius Grey White-eye (Z. mauritianus (Gmelin, 1789)) and searched for possible causes underlying similarities or differences between species living on the same or different islands. Both demographic inferences strongly support ancient and long-term expansions in all species. They also reveal different trajectories between species inhabiting different islands, but consistent demographic trajectories in species or populations from the same island. Species from Reunion appear to have experienced synchronous reductions in population size during the Last Glacial Maximum, a trend not seen in the Mauritian species. Overall, this study suggests that local events may have played a role in shaping population trajectories of these island species. It also highlights the potential of our conceptual framework to disentangle the effects of local and regional drivers on past species' demography and long-term population processes.

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.

Identification of Signatures of Selection for Litter Size and Pubertal Initiation in Two Sheep Populations

Fecundity is an important economic trait in sheep that directly affects their economic and productive efficiency. Our study aimed to identify SNP loci associated with sheep puberty or litter size which could be used in future breeding programs to improve fertility. Genomic DNA was obtained from Hetian and Cele Black sheep breeds and used for reduced-representation genome sequencing to identify SNP loci associated with pubertal initiation and litter size. Selective signatures analysis was performed based on the fixation index and nucleotide diversity, followed by pathway analysis of the genes contained in the selected regions. The selected SNP loci in the genes associated with pubertal initiation and litter size were validated using both sheep breeds. In total, 384,718 high quality SNPs were obtained and 376 genes were selected. Functional annotation of genes and enrichment analysis identified 12 genes associated with pubertal initiation and 11 genes associated with litter size. SNP locus validation showed that two SNP on PAK1 and four on ADCY1 may be associated with pubertal initiation, and one SNP on GNAQ gene (NC_040253.1: g.62677376G > A) was associated with litter size in Cele Black sheep. Our results provide new theoretical support for sheep breeding.

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 indicates that pathogen-mediated selection may play a role in shaping genomic variation among relatively recently colonised island bird populations and highlights the utility of genotype-environment associations for identifying candidate genes potentially involved in host-pathogen interactions.