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.

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.

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.

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.

Life history trade-offs at a single locus maintain sexually selected genetic variation.

Sexual selection, through intra-male competition or female choice, is assumed to be a source of strong and sustained directional selection in the wild. In the presence of such strong directional selection, alleles enhancing a particular trait are predicted to become fixed within a population, leading to a decrease in the underlying genetic variation. However, there is often considerable genetic variation underlying sexually selected traits in wild populations, and consequently, this phenomenon has become a long-discussed issue in the field of evolutionary biology. In wild Soay sheep, large horns confer an advantage in strong intra-sexual competition, yet males show an inherited polymorphism for horn type and have substantial genetic variation in their horn size. Here we show that most genetic variation in this trait is maintained by a trade-off between natural and sexual selection at a single gene, relaxin-like receptor 2 (RXFP2). We found that an allele conferring larger horns, Ho(+), is associated with higher reproductive success, whereas a smaller horn allele, Ho(P), confers increased survival, resulting in a net effect of overdominance (that is, heterozygote advantage) for fitness at RXFP2. The nature of this trade-off is simple relative to commonly proposed explanations for the maintenance of sexually selected traits, such as genic capture ('good genes') and sexually antagonistic selection. Our results demonstrate that by identifying the genetic architecture of trait variation, we can determine the principal mechanisms maintaining genetic variation in traits under strong selection and explain apparently counter-evolutionary observations.

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.

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.

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).

Long sperm fertilize more eggs in a bird.

Sperm competition, in which the ejaculates of multiple males compete to fertilize a female's ova, results in strong selection on sperm traits. Although sperm size and swimming velocity are known to independently affect fertilization success in certain species, exploring the relationship between sperm length, swimming velocity and fertilization success still remains a challenge. Here, we use the zebra finch (Taeniopygia guttata), where sperm size influences sperm swimming velocity, to determine the effect of sperm total length on fertilization success. Sperm competition experiments, in which pairs of males whose sperm differed only in length and swimming speed, revealed that males producing long sperm were more successful in terms of (i) the number of sperm reaching the ova and (ii) fertilizing those ova. Our results reveal that although sperm length is the main factor determining the outcome of sperm competition, complex interactions between male and female reproductive traits may also be important. The mechanisms underlying these interactions are poorly understood, but we suggest that differences in sperm storage and utilization by females may contribute to the outcome of sperm competition.

Localization of quantitative trait loci for diapause and other photoperiodically regulated life history traits important in adaptation to seasonally varying environments.

Seasonally changing environments at high latitudes present great challenges for the reproduction and survival of insects, and photoperiodic cues play an important role in helping them to synchronize their life cycle with prevalent and forthcoming conditions. We have mapped quantitative trait loci (QTL) responsible for the photoperiodic regulation of four life history traits, female reproductive diapause, cold tolerance, egg-to-eclosion development time and juvenile body weight in Drosophila montana strains from different latitudes in Canada and Finland. The F2 progeny of the cross was reared under a single photoperiod (LD cycle 16:8), which the flies from the Canadian population interpret as early summer and the flies from the Finnish population as late summer. The analysis revealed a unique QTL for diapause induction on the X chromosome and several QTL for this and the other measured traits on the 4th chromosome. Flies' cold tolerance, egg-to-eclosion development time and juvenile body weight had several QTL also on the 2nd, 3rd and 5th chromosome, some of the peaks overlapping with each other. These results suggest that while the downstream output of females' photoperiodic diapause response is partly under a different genetic control from that of the other traits in the given day length, all traits also share some QTL, possibly involving genes with pleiotropic effects and/or multiple tightly linked genes. Nonoverlapping QTL detected for some of the traits also suggest that the traits are potentially capable of independent evolution, even though this may be restricted by epistatic interactions and/or correlations and trade-offs between the traits.

Replicated analysis of the genetic architecture of quantitative traits in two wild great tit populations.

Currently, there is much debate on the genetic architecture of quantitative traits in wild populations. Is trait variation influenced by many genes of small effect or by a few genes of major effect? Where is additive genetic variation located in the genome? Do the same loci cause similar phenotypic variation in different populations? Great tits (Parus major) have been studied extensively in long-term studies across Europe and consequently are considered an ecological 'model organism'. Recently, genomic resources have been developed for the great tit, including a custom SNP chip and genetic linkage map. In this study, we used a suite of approaches to investigate the genetic architecture of eight quantitative traits in two long-term study populations of great tits--one in the Netherlands and the other in the United Kingdom. Overall, we found little evidence for the presence of genes of large effects in either population. Instead, traits appeared to be influenced by many genes of small effect, with conservative estimates of the number of contributing loci ranging from 31 to 310. Despite concordance between population-specific heritabilities, we found no evidence for the presence of loci having similar effects in both populations. While population-specific genetic architectures are possible, an undetected shared architecture cannot be rejected because of limited power to map loci of small and moderate effects. This study is one of few examples of genetic architecture analysis in replicated wild populations and highlights some of the challenges and limitations researchers will face when attempting similar molecular quantitative genetic studies in free-living populations.

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.

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.

High-utility conserved avian microsatellite markers enable parentage and population studies across a wide range of species.

BACKGROUND: Microsatellites are widely used for many genetic studies. In contrast to single nucleotide polymorphism (SNP) and genotyping-by-sequencing methods, they are readily typed in samples of low DNA quality/concentration (e.g. museum/non-invasive samples), and enable the quick, cheap identification of species, hybrids, clones and ploidy. Microsatellites also have the highest cross-species utility of all types of markers used for genotyping, but, despite this, when isolated from a single species, only a relatively small proportion will be of utility. Marker development of any type requires skill and time. The availability of sufficient "off-the-shelf" markers that are suitable for genotyping a wide range of species would not only save resources but also uniquely enable new comparisons of diversity among taxa at the same set of loci. No other marker types are capable of enabling this. We therefore developed a set of avian microsatellite markers with enhanced cross-species utility. RESULTS: We selected highly-conserved sequences with a high number of repeat units in both of two genetically distant species. Twenty-four primer sets were designed from homologous sequences that possessed at least eight repeat units in both the zebra finch (Taeniopygia guttata) and chicken (Gallus gallus). Each primer sequence was a complete match to zebra finch and, after accounting for degenerate bases, at least 86% similar to chicken. We assessed primer-set utility by genotyping individuals belonging to eight passerine and four non-passerine species. The majority of the new Conserved Avian Microsatellite (CAM) markers amplified in all 12 species tested (on average, 94% in passerines and 95% in non-passerines). This new marker set is of especially high utility in passerines, with a mean 68% of loci polymorphic per species, compared with 42% in non-passerine species. CONCLUSIONS: When combined with previously described conserved loci, this new set of conserved markers will not only reduce the necessity and expense of microsatellite isolation for a wide range of genetic studies, including avian parentage and population analyses, but will also now enable comparisons of genetic diversity among different species (and populations) at the same set of loci, with no or reduced bias. Finally, the approach used here can be applied to other taxa in which appropriate genome sequences are available.

Genome mapping in intensively studied wild vertebrate populations.

Over the past decade, long-term studies of vertebrate populations have been the focus of many quantitative genetic studies. As a result, we have a clearer understanding of why some fitness-related traits are heritable and under selection, but are apparently not evolving. An exciting extension of this work is to identify the genes underlying phenotypic variation in natural populations. The advent of next-generation sequencing and high-throughput single nucleotide polymorphism (SNP) genotyping platforms means that mapping studies are set to become widespread in those wild populations for whom appropriate phenotypic data and DNA samples are available. Here, we highlight the progress made in this area and define evolutionary genetic questions that have become tractable with the arrival of these new genomics technologies.

Pronounced inter- and intrachromosomal variation in linkage disequilibrium across the zebra finch genome.

The extent of nonrandom association of alleles at two or more loci, termed linkage disequilibrium (LD), can reveal much about population demography, selection, and recombination rate, and is a key consideration when designing association mapping studies. Here, we describe a genome-wide analysis of LD in the zebra finch (Taeniopygia guttata) using 838 single nucleotide polymorphisms and present LD maps for all assembled chromosomes. We found that LD declined with physical distance approximately five times faster on the microchromosomes compared to macrochromosomes. The distribution of LD across individual macrochromosomes also varied in a distinct pattern. In the center of the macrochromosomes there were large blocks of markers, sometimes spanning tens of mega bases, in strong LD whereas on the ends of macrochromosomes LD declined more rapidly. Regions of high LD were not simply the result of suppressed recombination around the centromere and this pattern has not been observed previously in other taxa. We also found evidence that this pattern of LD has remained stable across many generations. The variability in LD between and within chromosomes has important implications for genome wide association studies in birds and for our understanding of the distribution of recombination events and the processes that govern them.

Partitioning of genetic variation across the genome using multimarker methods in a wild bird population.

The underlying basis of genetic variation in quantitative traits, in terms of the number of causal variants and the size of their effects, is largely unknown in natural populations. The expectation is that complex quantitative trait variation is attributable to many, possibly interacting, causal variants, whose effects may depend upon the sex, age and the environment in which they are expressed. A recently developed methodology in animal breeding derives a value of relatedness among individuals from high-density genomic marker data, to estimate additive genetic variance within livestock populations. Here, we adapt and test the effectiveness of these methods to partition genetic variation for complex traits across genomic regions within ecological study populations where individuals have varying degrees of relatedness. We then apply this approach for the first time to a natural population and demonstrate that genetic variation in wing length in the great tit (Parus major) reflects contributions from multiple genomic regions. We show that a polygenic additive mode of gene action best describes the patterns observed, and we find no evidence of dosage compensation for the sex chromosome. Our results suggest that most of the genomic regions that influence wing length have the same effects in both sexes. We found a limited amount of genetic variance in males that is attributed to regions that have no effects in females, which could facilitate the sexual dimorphism observed for this trait. Although this exploratory work focuses on one complex trait, the methodology is generally applicable to any trait for any laboratory or wild population, paving the way for investigating sex-, age- and environment-specific genetic effects and thus the underlying genetic architecture of phenotype in biological study systems.

Genomic dissection of variation in clutch size and egg mass in a wild great tit (Parus major) population.

Clutch size and egg mass are life history traits that have been extensively studied in wild bird populations, as life history theory predicts a negative trade-off between them, either at the phenotypic or at the genetic level. Here, we analyse the genomic architecture of these heritable traits in a wild great tit (Parus major) population, using three marker-based approaches - chromosome partitioning, quantitative trait locus (QTL) mapping and a genome-wide association study (GWAS). The variance explained by each great tit chromosome scales with predicted chromosome size, no location in the genome contains genome-wide significant QTL, and no individual SNPs are associated with a large proportion of phenotypic variation, all of which may suggest that variation in both traits is due to many loci of small effect, located across the genome. There is no evidence that any regions of the genome contribute significantly to both traits, which combined with a small, nonsignificant negative genetic covariance between the traits, suggests the absence of genetic constraints on the independent evolution of these traits. Our findings support the hypothesis that variation in life history traits in natural populations is likely to be determined by many loci of small effect spread throughout the genome, which are subject to continued input of variation by mutation and migration, although we cannot exclude the possibility of an additional input of major effect genes influencing either trait.

Introgression and the fate of domesticated genes in a wild mammal population.

When domesticated species are not reproductively isolated from their wild relatives, the opportunity arises for artificially selected variants to be re-introduced into the wild. However, the evolutionary consequences of introgression of domesticated genes back into the wild are poorly understood. By combining high-throughput genotyping with 25 years of long-term ecological field data, we describe the occurrence and consequences of admixture between a primitive sheep breed, the free-living Soay sheep of St Kilda, and more modern breeds. Utilizing data from a 50 K ovine SNP chip, together with forward simulations of demographic scenarios, we show that admixture occurred between Soay sheep and a more modern breed, consistent with historical accounts, approximately 150 years ago. Haplotype-sharing analyses with other breeds revealed that polymorphisms in coat colour and pattern in Soay sheep arose as a result of introgression of genetic variants favoured by artificial selection. Because the haplotypes carrying the causative mutations are known to be under natural selection in free-living Soay sheep, the admixture event created an opportunity to observe the outcome of a 'natural laboratory' experiment where ancestral and domesticated genes competed with each other. The haplotype carrying the domesticated light coat colour allele was favoured by natural selection, while the haplotype associated with the domesticated self coat pattern allele was associated with decreased survival. Therefore, we demonstrate that introgression of domesticated alleles into wild populations can provide a novel source of variation capable of generating rapid evolutionary changes.