The great tit HapMap project: A continental-scale analysis of genomic variation in a songbird.

A major aim of evolutionary biology is to understand why patterns of genomic diversity vary within taxa and space. Large-scale genomic studies of widespread species are useful for studying how environment and demography shape patterns of genomic divergence. Here, we describe one of the most geographically comprehensive surveys of genomic variation in a wild vertebrate to date; the great tit (Parus major) HapMap project. We screened ca 500,000 SNP markers across 647 individuals from 29 populations, spanning ~30 degrees of latitude and 40 degrees of longitude - almost the entire geographical range of the European subspecies. Genome-wide variation was consistent with a recent colonisation across Europe from a South-East European refugium, with bottlenecks and reduced genetic diversity in island populations. Differentiation across the genome was highly heterogeneous, with clear 'islands of differentiation', even among populations with very low levels of genome-wide differentiation. Low local recombination rates were a strong predictor of high local genomic differentiation (FST), especially in island and peripheral mainland populations, suggesting that the interplay between genetic drift and recombination causes highly heterogeneous differentiation landscapes. We also detected genomic outlier regions that were confined to one or more peripheral great tit populations, probably as a result of recent directional selection at the species' range edges. Haplotype-based measures of selection were related to recombination rate, albeit less strongly, and highlighted population-specific sweeps that likely resulted from positive selection. Our study highlights how comprehensive screens of genomic variation in wild organisms can provide unique insights into spatio-temporal evolutionary dynamics.

A sex-linked supergene with large effects on sperm traits has little impact on reproductive traits in female zebra finches.

Despite constituting an essential component of fitness, reproductive success can vary remarkably between individuals and the causes of such variation are not well understood across taxa. In the zebra finch-a model songbird, almost all the variation in sperm morphology and swimming speed is maintained by a large polymorphic inversion (commonly known as a supergene) on the Z chromosome. The relationship between this polymorphism and reproductive success is not fully understood, particularly for females. Here, we explore the effects of female haplotype, and the combination of male and female genotype, on several primary reproductive traits in a captive population of zebra finches. Despite the inversion polymorphism's known effects on sperm traits, we find no evidence that inversion haplotype influences egg production by females or survival of embryos through to hatching. However, our findings do reinforce existing evidence that the inversion polymorphism is maintained by a heterozygote advantage for male fitness. This work provides an important step in understanding the causes of variation in reproductive success in this model species.

The genomics of adaptation to climate in European great tit (Parus major) populations.

The recognition that climate change is occurring at an unprecedented rate means that there is increased urgency in understanding how organisms can adapt to a changing environment. Wild great tit (Parus major) populations represent an attractive ecological model system to understand the genomics of climate adaptation. They are widely distributed across Eurasia and they have been documented to respond to climate change. We performed a Bayesian genome-environment analysis, by combining local climate data with single nucleotide polymorphisms genotype data from 20 European populations (broadly spanning the species' continental range). We found 36 genes putatively linked to adaptation to climate. Following an enrichment analysis of biological process Gene Ontology (GO) terms, we identified over-represented terms and pathways among the candidate genes. Because many different genes and GO terms are associated with climate variables, it seems likely that climate adaptation is polygenic and genetically complex. Our findings also suggest that geographical climate adaptation has been occurring since great tits left their Southern European refugia at the end of the last ice age. Finally, we show that substantial climate-associated genetic variation remains, which will be essential for adaptation to future changes.

DNA methylation differences between stick insect ecotypes.

Epigenetic mechanisms, such as DNA methylation, can influence gene regulation and affect phenotypic variation, raising the possibility that they contribute to ecological adaptation. Beginning to address this issue requires high-resolution sequencing studies of natural populations to pinpoint epigenetic regions of potential ecological and evolutionary significance. However, such studies are still relatively uncommon, especially in insects, and are mainly restricted to a few model organisms. Here, we characterize patterns of DNA methylation for natural populations of Timema cristinae adapted to two host plant species (i.e. ecotypes). By integrating results from sequencing of whole transcriptomes, genomes and methylomes, we investigate whether environmental, host and genetic differences of these stick insects are associated with methylation levels of cytosine nucleotides in the CpG context. We report an overall genome-wide methylation level for T. cristinae of ~14%, with methylation being enriched in gene bodies and impoverished in repetitive elements. Genome-wide DNA methylation variation was strongly positively correlated with genetic distance (relatedness), but also exhibited significant host-plant effects. Using methylome-environment association analysis, we pinpointed specific genomic regions that are differentially methylated between ecotypes, with these regions being enriched for genes with functions in membrane processes. The observed association between methylation variation and genetic relatedness, and with the ecologically important variable of host plant, suggests a potential role for epigenetic modification in T. cristinae adaptation. To substantiate such adaptive significance, future studies could test whether methylation can be transmitted across generations and the extent to which it responds to experimental manipulation in field and laboratory studies.

A polygenic basis for birth weight in a wild population of red deer (Cervus elaphus).

The genetic architecture of traits under selection has important consequences for the response to selection and potentially for population viability. Early QTL mapping studies in wild populations have reported loci with large effect on trait variation. However, these results are contradicted by more recent genome-wide association analyses, which strongly support the idea that most quantitative traits have a polygenic basis. This study aims to re-evaluate the genetic architecture of a key morphological trait, birth weight, in a wild population of red deer (Cervus elaphus), using genomic approaches. A previous study using 93 microsatellite and allozyme markers and linkage mapping on a kindred of 364 deer detected a pronounced QTL on chromosome 21 explaining 29% of the variance in birth weight, suggesting that this trait is partly controlled by genes with large effects. Here, we used data for more than 2,300 calves genotyped at >39,000 SNP markers and two approaches to characterise the genetic architecture of birth weight. First, we performed a genome-wide association (GWA) analysis, using a genomic relatedness matrix to account for population structure. We found no SNPs significantly associated with birth weight. Second, we used genomic prediction to estimate the proportion of variance explained by each SNP and chromosome. This analysis confirmed that most genetic variance in birth weight was explained by loci with very small effect sizes. Third, we found that the proportion of variance explained by each chromosome was slightly positively correlated with its size. These three findings highlight a highly polygenic architecture for birth weight, which contradicts the previous QTL study. These results are probably explained by the differences in how associations are modelled between QTL mapping and GWA. Our study suggests that models of polygenic adaptation are the most appropriate to study the evolutionary trajectory of this trait.

Genetic variance in fitness indicates rapid contemporary adaptive evolution in wild animals.

The rate of adaptive evolution, the contribution of selection to genetic changes that increase mean fitness, is determined by the additive genetic variance in individual relative fitness. To date, there are few robust estimates of this parameter for natural populations, and it is therefore unclear whether adaptive evolution can play a meaningful role in short-term population dynamics. We developed and applied quantitative genetic methods to long-term datasets from 19 wild bird and mammal populations and found that, while estimates vary between populations, additive genetic variance in relative fitness is often substantial and, on average, twice that of previous estimates. We show that these rates of contemporary adaptive evolution can affect population dynamics and hence that natural selection has the potential to partly mitigate effects of current environmental change.

Genomic prediction in the wild: A case study in Soay sheep.

Genomic prediction, the technique whereby an individual's genetic component of their phenotype is estimated from its genome, has revolutionised animal and plant breeding and medical genetics. However, despite being first introduced nearly two decades ago, it has hardly been adopted by the evolutionary genetics community studying wild organisms. Here, genomic prediction is performed on eight traits in a wild population of Soay sheep. The population has been the focus of a >30 year evolutionary ecology study and there is already considerable understanding of the genetic architecture of the focal Mendelian and quantitative traits. We show that the accuracy of genomic prediction is high for all traits, but especially those with loci of large effect segregating. Five different methods are compared, and the two methods that can accommodate zero-effect and large-effect loci in the same model tend to perform best. If the accuracy of genomic prediction is similar in other wild populations, then there is a real opportunity for pedigree-free molecular quantitative genetics research to be enabled in many more wild populations; currently the literature is dominated by studies that have required decades of field data collection to generate sufficiently deep pedigrees. Finally, some of the potential applications of genomic prediction in wild populations are discussed.

Cost, risk, and avoidance of inbreeding in a cooperatively breeding bird.

Inbreeding is often avoided in natural populations by passive processes such as sex-biased dispersal. But, in many social animals, opposite-sexed adult relatives are spatially clustered, generating a risk of incest and hence selection for active inbreeding avoidance. Here we show that, in long-tailed tits (Aegithalos caudatus), a cooperative breeder that risks inbreeding by living alongside opposite-sex relatives, inbreeding carries fitness costs and is avoided by active kin discrimination during mate choice. First, we identified a positive association between heterozygosity and fitness, indicating that inbreeding is costly. We then compared relatedness within breeding pairs to that expected under multiple mate-choice models, finding that pair relatedness is consistent with avoidance of first-order kin as partners. Finally, we show that the similarity of vocal cues offers a plausible mechanism for discrimination against first-order kin during mate choice. Long-tailed tits are known to discriminate between the calls of close kin and nonkin, and they favor first-order kin in cooperative contexts, so we conclude that long-tailed tits use the same kin discrimination rule to avoid inbreeding as they do to direct help toward kin.

Phenotypic Responses to and Genetic Architecture of Sterility Following Exposure to Sub-Lethal Temperature During Development.

Thermal tolerance range, based on temperatures that result in incapacitating effects, influences species' distributions and has been used to predict species' response to increasing temperature. Reproductive performance may also be negatively affected at less extreme temperatures, but such sublethal heat-induced sterility has been relatively ignored in studies addressing the potential effects of, and ability of species' to respond to, predicted climate warming. The few studies examining the link between increased temperature and reproductive performance typically focus on adults, although effects can vary between life history stages. Here we assessed how sublethal heat stress during development impacted subsequent adult fertility and its plasticity, both of which can provide the raw material for evolutionary responses to increased temperature. We quantified phenotypic and genetic variation in fertility of Drosophila melanogaster reared at standardized densities in three temperatures (25, 27, and 29°C) from a set of lines of the Drosophila Genetic Reference Panel (DGRP). We found little phenotypic variation at the two lower temperatures with more variation at the highest temperature and for plasticity. Males were more affected than females. Despite reasonably large broad-sense heritabilities, a genome-wide association study found little evidence for additive genetic variance and no genetic variants were robustly linked with reproductive performance at specific temperatures or for phenotypic plasticity. We compared results on heat-induced male sterility with other DGRP results on relevant fitness traits measured after abiotic stress and found an association between male susceptibility to sterility and male lifespan reduction following oxidative stress. Our results suggest that sublethal stress during development has profound negative consequences on male adult reproduction, but despite phenotypic variation in a population for this response, there is limited evolutionary potential, either through adaptation to a specific developmental temperature or plasticity in response to developmental heat-induced sterility.

A genome-wide linkage map for the house sparrow (Passer domesticus) provides insights into the evolutionary history of the avian genome.

The house sparrow is an important model species for studying physiological, ecological and evolutionary processes in wild populations. Here, we present a medium density, genome wide linkage map for house sparrow (Passer domesticus) that has aided the assembly of the house sparrow reference genome, and that will provide an important resource for ongoing mapping of genes controlling important traits in the ecology and evolution of this species. Using a custom house sparrow 10 K iSelect Illumina SNP chip we have assigned 6,498 SNPs to 29 autosomal linkage groups, based on a mean of 430 informative meioses per SNP. The map was constructed by combining the information from linkage with that of the physical position of SNPs within scaffold sequences in an iterative process. Averaged between the sexes; the linkage map had a total length of 2,004 cM, with a longer map for females (2,240 cM) than males (1,801 cM). Additionally, recombination rates also varied along the chromosomes. Comparison of the linkage map to the reference genomes of zebra finch, collared flycatcher and chicken, showed a chromosome fusion of the two avian chromosomes 8 and 4A in house sparrow. Lastly, information from the linkage map was utilized to conduct analysis of linkage disequilibrium (LD) in eight populations with different effective population sizes (Ne ) in order to quantify the background level LD. Together, these results aid the design of future association studies, facilitate the development of new genomic tools and support the body of research that describes the evolution of the avian genome.

The Genomic Complexity of a Large Inversion in Great Tits.

Chromosome inversions have clear effects on genome evolution and have been associated with speciation, adaptation, and the evolution of the sex chromosomes. In birds, these inversions may play an important role in hybridization of species and disassortative mating. We identified a large (≈64 Mb) inversion polymorphism in the great tit (Parus major) that encompasses almost 1,000 genes and more than 90% of Chromosome 1A. The inversion occurs at a low frequency in a set of over 2,300 genotyped great tits in the Netherlands with only 5% of the birds being heterozygous for the inversion. In an additional analysis of 29 resequenced birds from across Europe, we found two heterozygotes. The likely inversion breakpoints show considerable genomic complexity, including multiple copy number variable segments. We identified different haplotypes for the inversion, which differ in the degree of recombination in the center of the chromosome. Overall, this remarkable genetic variant is widespread among distinct great tit populations and future studies of the inversion haplotype, including how it affects the fitness of carriers, may help to understand the mechanisms that maintain it.

Recent natural selection causes adaptive evolution of an avian polygenic trait.

We used extensive data from a long-term study of great tits (Parus major) in the United Kingdom and Netherlands to better understand how genetic signatures of selection translate into variation in fitness and phenotypes. We found that genomic regions under differential selection contained candidate genes for bill morphology and used genetic architecture analyses to confirm that these genes, especially the collagen gene COL4A5, explained variation in bill length. COL4A5 variation was associated with reproductive success, which, combined with spatiotemporal patterns of bill length, suggested ongoing selection for longer bills in the United Kingdom. Last, bill length and COL4A5 variation were associated with usage of feeders, suggesting that longer bills may have evolved in the United Kingdom as a response to supplementary feeding.

A sex-linked supergene controls sperm morphology and swimming speed in a songbird.

Sperm competition is an important selective force in many organisms. As a result, sperm have evolved to be among the most diverse cells in the animal kingdom. However, the relationship between sperm morphology, sperm motility and fertilization success is only partially understood. The extent to which between-male variation is heritable is largely unknown, and remarkably few studies have investigated the genetic architecture of sperm traits, especially sperm morphology. Here we use high-density genotyping and gene expression profiling to explore the considerable sperm trait variation that exists in the zebra finch Taeniopygia guttata. We show that nearly all of the genetic variation in sperm morphology is caused by an inversion polymorphism on the Z chromosome acting as a 'supergene'. These results provide a striking example of two evolutionary genetic predictions. First, that in species where females are the heterogametic sex, genetic variation affecting sexually dimorphic traits will accumulate on the Z chromosome. Second, recombination suppression at the inversion allows beneficial dominant alleles to become fixed on whichever haplotype they first arise, without being exchanged onto other haplotypes. Finally, we show that the inversion polymorphism will be stably maintained by heterozygote advantage, because heterozygous males have the fastest and most successful sperm.

Determinants of the Efficacy of Natural Selection on Coding and Noncoding Variability in Two Passerine Species.

Population genetic theory predicts that selection should be more effective when the effective population size (Ne) is larger, and that the efficacy of selection should correlate positively with recombination rate. Here, we analyzed the genomes of ten great tits and ten zebra finches. Nucleotide diversity at 4-fold degenerate sites indicates that zebra finches have a 2.83-fold larger Ne. We obtained clear evidence that purifying selection is more effective in zebra finches. The proportion of substitutions at 0-fold degenerate sites fixed by positive selection (α) is high in both species (great tit 48%; zebra finch 64%) and is significantly higher in zebra finches. When α was estimated on GC-conservative changes (i.e., between A and T and between G and C), the estimates reduced in both species (great tit 22%; zebra finch 53%). A theoretical model presented herein suggests that failing to control for the effects of GC-biased gene conversion (gBGC) is potentially a contributor to the overestimation of α, and that this effect cannot be alleviated by first fitting a demographic model to neutral variants. We present the first estimates in birds for α in the untranslated regions, and found evidence for substantial adaptive changes. Finally, although purifying selection is stronger in high-recombination regions, we obtained mixed evidence for α increasing with recombination rate, especially after accounting for gBGC. These results highlight that it is important to consider the potential confounding effects of gBGC when quantifying selection and that our understanding of what determines the efficacy of selection is incomplete.

Social and spatial effects on genetic variation between foraging flocks in a wild bird population.

Social interactions are rarely random. In some instances, animals exhibit homophily or heterophily, the tendency to interact with similar or dissimilar conspecifics, respectively. Genetic homophily and heterophily influence the evolutionary dynamics of populations, because they potentially affect sexual and social selection. Here, we investigate the link between social interactions and allele frequencies in foraging flocks of great tits (Parus major) over three consecutive years. We constructed co-occurrence networks which explicitly described the splitting and merging of 85,602 flocks through time (fission-fusion dynamics), at 60 feeding sites. Of the 1,711 birds in those flocks, we genotyped 962 individuals at 4,701 autosomal single nucleotide polymorphisms (SNPs). By combining genomewide genotyping with repeated field observations of the same individuals, we were able to investigate links between social structure and allele frequencies at a much finer scale than was previously possible. We explicitly accounted for potential spatial effects underlying genetic structure at the population level. We modelled social structure and spatial configuration of great tit fission-fusion dynamics with eigenvector maps. Variance partitioning revealed that allele frequencies were strongly affected by group fidelity (explaining 27%-45% of variance) as individuals tended to maintain associations with the same conspecifics. These conspecifics were genetically more dissimilar than expected, shown by genomewide heterophily for pure social (i.e., space-independent) grouping preferences. Genomewide homophily was linked to spatial configuration, indicating spatial segregation of genotypes. We did not find evidence for homophily or heterophily for putative socially relevant candidate genes or any other SNP markers. Together, these results demonstrate the importance of distinguishing social and spatial processes in determining population structure.

Robust inference of genetic architecture in mapping studies.

The genetic architecture of a trait usually refers to the number and magnitude of loci that explain phenotypic variation. A description of genetic architecture can help us to understand how genetic variation is maintained, how traits have evolved and how phenotypes might respond to selection. However, linkage mapping and association studies can suffer from problems of bias, especially when conducted in natural populations where the opportunity to perform studies with very large sample sizes can be limited. In this issue of Molecular Ecology, Li and colleagues perform an association study of brain traits in ninespine sticklebacks Pungitius pungitius. They use a sophisticated approach that models all of the genotyped markers simultaneously; conventional approaches fit each marker individually. Although the single-marker and multi-marker approaches find similar regions of the genome that explain phenotypic variation, the overall conclusions about trait architecture are somewhat different, depending on the approach used. Single-marker methods identify regions that explain quite large proportions of genetic variation, whereas the multi-marker approach suggests the traits are far more polygenic. Simulations suggest the multi-marker approach is robust. This study highlights how molecular quantitative genetics in wild populations can be used to address hypothesis-driven questions, without making unrealistic assumptions about effect sizes of individual quantitative trait loci.

Sperm morphology, adenosine triphosphate (ATP) concentration and swimming velocity: unexpected relationships in a passerine bird.

The relationship between sperm energetics and sperm function is poorly known, but is central to our understanding of the evolution of sperm traits. The aim of this study was to examine how sperm morphology and ATP content affect sperm swimming velocity in the zebra finch Taeniopygia guttata We exploited the high inter-male variation in this species and created extra experimental power by increasing the number of individuals with very long or short sperm through artificial selection. We found a pronounced quadratic relationship between total sperm length and swimming velocity, with velocity increasing with length up to a point, but declining in the very longest sperm. We also found an unexpected negative association between midpiece length and ATP content: sperm with a short midpiece generally contained the highest concentration of ATP. Low intracellular ATP is therefore unlikely to explain reduced swimming velocity among the very longest sperm (which tend to have a shorter midpiece).