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.

Genomic inference of contemporary effective population size in a large island population of collared flycatchers (Ficedula albicollis).

Due to its central importance to many aspects of evolutionary biology and population genetics, the long-term effective population size (Ne ) has been estimated for numerous species and populations. However, estimating contemporary Ne is difficult and in practice this parameter is often unknown. In principle, contemporary Ne can be estimated using either analyses of temporal changes in allele frequencies, or the extent of linkage disequilibrium (LD) between unlinked markers. We applied these approaches to estimate contemporary Ne of a relatively recently founded island population of collared flycatchers (Ficedula albicollis). We sequenced the genomes of 85 birds sampled in 1993 and 2015, and applied several temporal methods to estimate Ne at a few thousand (4000-7000). The approach based on LD provided higher estimates of Ne (20,000-32,000) and was associated with high variance, often resulting in infinite Ne . We conclude that whole-genome sequencing data offers new possibilities to estimate high (>1000) contemporary Ne , but also note that such estimates remain challenging, in particular for LD-based methods for contemporary Ne estimation.

Footprints of adaptive evolution revealed by whole Z chromosomes haplotypes in flycatchers.

Detecting positive selection using genomic data is critical to understanding the role of adaptive evolution. Of particular interest in this context is sex chromosomes since they are thought to play a special role in local adaptation and speciation. We sought to circumvent the challenges associated with statistical phasing when using haplotype-based statistics in sweep scans by benefitting from that whole chromosome haplotypes of the sex chromosomes can be obtained by resequencing of individuals of the hemizygous sex. We analyzed whole Z chromosome haplotypes from 100 females from several populations of four black and white flycatcher species (in birds, females are ZW and males ZZ). Based on integrated haplotype score (iHS) and number of segregating sites by length (nSL) statistics, we found strong and frequent haplotype structure in several regions of the Z chromosome in each species. Most of these sweep signals were population-specific, with essentially no evidence for regions under selection shared among species. Some completed sweeps were revealed by the cross-population extended haplotype homozygosity (XP-EHH) statistic. Importantly, by using statistically phased Z chromosome data from resequencing of males, we failed to recover the signals of selection detected in analyses based on whole chromosome haplotypes from females; instead, what likely represent false signals of selection were frequently seen. This highlights the power issues in statistical phasing and cautions against conclusions from selection scans using such data. The detection of frequent selective sweeps on the avian Z chromosome supports a large role of sex chromosomes in adaptive evolution.

Sex-biased gene expression, sexual antagonism and levels of genetic diversity in the collared flycatcher (Ficedula albicollis) genome.

Theoretical work suggests that sexual conflict should promote the maintenance of genetic diversity by the opposing directions of selection on males and females. If such conflict is pervasive, it could potentially lead to genomic heterogeneity in levels of genetic diversity an idea that so far has not been empirically tested on a genomewide scale. We used large-scale population genomic and transcriptomic data from the collared flycatcher (Ficedula albicollis) to analyse how sexual conflict, for which we use sex-biased gene expression as a proxy, relates to genetic variability. Here, we demonstrate that the extent of sex-biased gene expression of both male-biased and female-biased genes is significantly correlated with levels of nucleotide diversity in gene sequences and that this correlation extends to diversity levels also in intergenic DNA and introns. We find signatures of balancing selection in sex-biased genes but also note that relaxed purifying selection could potentially explain part of the observed patterns. The finding of significant genetic differentiation between males and females for male-biased (and gonad-specific) genes indicates ongoing sexual conflict and sex-specific viability selection, potentially driven by sexual selection. Our results thus indicate that sexual antagonism could potentially be considered as one viable explanation to the long-standing question in evolutionary biology of how genomes can remain so genetically variable in face of strong natural and sexual selection.

The genomic landscape of species divergence in Ficedula flycatchers.

Unravelling the genomic landscape of divergence between lineages is key to understanding speciation. The naturally hybridizing collared flycatcher and pied flycatcher are important avian speciation models that show pre- as well as postzygotic isolation. We sequenced and assembled the 1.1-Gb flycatcher genome, physically mapped the assembly to chromosomes using a low-density linkage map and re-sequenced population samples of each species. Here we show that the genomic landscape of species differentiation is highly heterogeneous with approximately 50 'divergence islands' showing up to 50-fold higher sequence divergence than the genomic background. These non-randomly distributed islands, with between one and three regions of elevated divergence per chromosome irrespective of chromosome size, are characterized by reduced levels of nucleotide diversity, skewed allele-frequency spectra, elevated levels of linkage disequilibrium and reduced proportions of shared polymorphisms in both species, indicative of parallel episodes of selection. Proximity of divergence peaks to genomic regions resistant to sequence assembly, potentially including centromeres and telomeres, indicate that complex repeat structures may drive species divergence. A much higher background level of species divergence of the Z chromosome, and a lower proportion of shared polymorphisms, indicate that sex chromosomes and autosomes are at different stages of speciation. This study provides a roadmap to the emerging field of speciation genomics.

Life History Traits, Protein Evolution, and the Nearly Neutral Theory in Amniotes.

The nearly neutral theory of molecular evolution predicts that small populations should accumulate deleterious mutations at a faster rate than large populations. The analysis of nonsynonymous (dN) versus synonymous (dS) substitution rates in birds versus mammals, however, has provided contradictory results, questioning the generality of the nearly neutral theory. Here we analyzed the impact of life history traits, taken as proxies of the effective population size, on molecular evolutionary and population genetic processes in amniotes, including the so far neglected reptiles. We report a strong effect of species body mass, longevity, and age of sexual maturity on genome-wide patterns of polymorphism and divergence across the major groups of amniotes, in agreement with the nearly neutral theory. Our results indicate that the rate of protein evolution in amniotes is determined in the first place by the efficiency of purifying selection against deleterious mutations-and this is true of both radical and conservative amino acid changes. Interestingly, the among-species distribution of dN/dS in birds did not follow this general trend: dN/dS was not higher in large, long-lived than in small, short-lived species of birds. We show that this unexpected pattern is not due to a more narrow range of life history traits, a lack of correlation between traits and Ne, or a peculiar distribution of fitness effects of mutations in birds. Our analysis therefore highlights the bird dN/dS ratio as a molecular evolutionary paradox and a challenge for future research.

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 (Ne ) 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 Ne 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. Ne curves of Atlas, collared, pied and semicollared flycatcher converged 1-2 million years ago (Ma) at an Ne of ≈ 200 000, likely reflecting the time when all four species last shared a common ancestor. Subsequent separate Ne 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 Ne differed among species, and for pied flycatcher, the temporal dynamics of Ne 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.

Demographic divergence history of pied flycatcher and collared flycatcher inferred from whole-genome re-sequencing data.

Profound knowledge of demographic history is a prerequisite for the understanding and inference of processes involved in the evolution of population differentiation and speciation. Together with new coalescent-based methods, the recent availability of genome-wide data enables investigation of differentiation and divergence processes at unprecedented depth. We combined two powerful approaches, full Approximate Bayesian Computation analysis (ABC) and pairwise sequentially Markovian coalescent modeling (PSMC), to reconstruct the demographic history of the split between two avian speciation model species, the pied flycatcher and collared flycatcher. Using whole-genome re-sequencing data from 20 individuals, we investigated 15 demographic models including different levels and patterns of gene flow, and changes in effective population size over time. ABC provided high support for recent (mode 0.3 my, range <0.7 my) species divergence, declines in effective population size of both species since their initial divergence, and unidirectional recent gene flow from pied flycatcher into collared flycatcher. The estimated divergence time and population size changes, supported by PSMC results, suggest that the ancestral species persisted through one of the glacial periods of middle Pleistocene and then split into two large populations that first increased in size before going through severe bottlenecks and expanding into their current ranges. Secondary contact appears to have been established after the last glacial maximum. The severity of the bottlenecks at the last glacial maximum is indicated by the discrepancy between current effective population sizes (20,000-80,000) and census sizes (5-50 million birds) of the two species. The recent divergence time challenges the supposition that avian speciation is a relatively slow process with extended times for intrinsic postzygotic reproductive barriers to evolve. Our study emphasizes the importance of using genome-wide data to unravel tangled demographic histories. Moreover, it constitutes one of the first examples of the inference of divergence history from genome-wide data in non-model species.

Interspecific hybridization increases MHC class II diversity in two sister species of newts.

Our understanding of the evolutionary mechanisms generating variation within the highly polymorphic major histocompatibility complex (MHC) genes remains incomplete. Assessing MHC variation across multiple populations, of recent and ancient divergence, may facilitate understanding of geographical and temporal aspects of variation. Here, we applied 454 sequencing to perform a large-scale, comprehensive analysis of MHC class II in the closely related, hybridizing newts, Lissotriton vulgaris (Lv) and Lissotriton montandoni (Lm). Our study revealed an extensive (299 alleles) geographically structured polymorphism. Populations at the southern margin of the Lv distribution, inhabited by old and distinct lineages (southern Lv), exhibited moderate MHC variation and strong population structure, indicating little gene flow or extensive local adaptation. Lissotriton vulgaris in central Europe and the northern Balkans (northern Lv) and almost all Lm populations had a high MHC variation. A much higher proportion of MHC alleles was shared between Lm and northern Lv than between Lm and southern Lv. Strikingly, the average pairwise F(ST) between northern Lv and Lm was significantly lower than between northern and southern Lv for MHC, but not for microsatellites. Thus, high MHC variation in Lm and northern Lv may result from gene flow between species. We hypothesize that the interspecific exchange of MHC genes may be facilitated by frequency-dependent selection. A marginally significant correlation between the MHC and microsatellite allelic richness indicates that demographic factors may have contributed to the present-day pattern of MHC variation, but unequivocal signatures of adaptive evolution in MHC class II sequences emphasize the role of selection on a longer timescale.

Development and characterization of microsatellite loci in the Centricnemus leucogrammus weevil.

Centricnemus leucogrammus is a weevil characteristic of European xerothermic habitats and steppes. The species was probably more widespread during the Pleistocene glaciations, while its current distribution is limited to "warm-stage refugia." It may be regarded as a typical representative of flightless xerothermophilous beetles. Previous studies concentrated on its genetic variation using mitochondrial genes. Here, we identified, tested and characterized 24 polymorphic microsatellite loci with the use of 454 sequencing of microsatellite enriched genomic libraries. The new set of loci will be used in studies on the population structure of this weevil and may provide valuable information for its conservation.

Temporal Dynamics of Avian Populations during Pleistocene Revealed by Whole-Genome Sequences.

Global climate fluctuations have significantly influenced the distribution and abundance of biodiversity. During unfavorable glacial periods, many species experienced range contraction and fragmentation, expanding again during interglacials. An understanding of the evolutionary consequences of both historical and ongoing climate changes requires knowledge of the temporal dynamics of population numbers during such climate cycles. Variation in abundance should have left clear signatures in the patterns of intraspecific genetic variation in extant species, from which historical effective population sizes (N(e)) can be estimated. We analyzed whole-genome sequences of 38 avian species in a pairwise sequentially Markovian coalescent (PSMC, [5]) framework to quantitatively reveal changes in N(e) from approximately 10 million to 10 thousand years ago. Significant fluctuations in N(e) over time were evident for most species. The most pronounced pattern observed in many species was a severe reduction in N(e) coinciding with the beginning of the last glacial period (LGP). Among species, N(e) varied by at least three orders of magnitude, exceeding 1 million in the most abundant species. Several species on the IUCN Red List of Threatened Species showed long-term reduction in population size, predating recent declines. We conclude that cycles of population expansions and contractions have been a common feature of many bird species during the Quaternary period, likely coinciding with climate cycles. Population size reduction should have increased the risk of extinction but may also have promoted speciation. Species that have experienced long-term declines may be especially vulnerable to recent anthropogenic threats.