Insights into the genetic architecture of morphological traits in two passerine bird species.

Knowledge about the underlying genetic architecture of phenotypic traits is needed to understand and predict evolutionary dynamics. The number of causal loci, magnitude of the effects and location in the genome are, however, still largely unknown. Here, we use genome-wide single-nucleotide polymorphism (SNP) data from two large-scale data sets on house sparrows and collared flycatchers to examine the genetic architecture of different morphological traits (tarsus length, wing length, body mass, bill depth, bill length, total and visible badge size and white wing patches). Genomic heritabilities were estimated using relatedness calculated from SNPs. The proportion of variance captured by the SNPs (SNP-based heritability) was lower in house sparrows compared with collared flycatchers, as expected given marker density (6348 SNPs in house sparrows versus 38 689 SNPs in collared flycatchers). Indeed, after downsampling to similar SNP density and sample size, this estimate was no longer markedly different between species. Chromosome-partitioning analyses demonstrated that the proportion of variance explained by each chromosome was significantly positively related to the chromosome size for some traits and, generally, that larger chromosomes tended to explain proportionally more variation than smaller chromosomes. Finally, we found two genome-wide significant associations with very small-effect sizes. One SNP on chromosome 20 was associated with bill length in house sparrows and explained 1.2% of phenotypic variation (VP), and one SNP on chromosome 4 was associated with tarsus length in collared flycatchers (3% of VP). Although we cannot exclude the possibility of undetected large-effect loci, our results indicate a polygenic basis for morphological traits.

An extensive candidate gene approach to speciation: diversity, divergence and linkage disequilibrium in candidate pigmentation genes across the European crow hybrid zone.

Colouration patterns have an important role in adaptation and speciation. The European crow system, in which all-black carrion crows and grey-coated hooded crows meet in a narrow hybrid zone, is a prominent example. The marked phenotypic difference is maintained by assortative mating in the absence of neutral genetic divergence, suggesting the presence of few pigmentation genes of major effect. We made use of the rich phenotypic and genetic resources in mammals and identified a comprehensive panel of 95 candidate pigmentation genes for birds. Based on functional annotation, we chose a subset of the most promising 37 candidates, for which we developed a marker system that demonstrably works across the avian phylogeny. In total, we sequenced 107 amplicons (∼3 loci per gene, totalling 60 kb) in population samples of crows (n=23 for each taxon). Tajima's D, Fu's FS, DHEW and HKA (Hudson-Kreitman-Aguade) statistics revealed several amplicons that deviated from neutrality; however, none of these showed significantly elevated differentiation between the two taxa. Hence, colour divergence in this system may be mediated by uncharacterized pigmentation genes or regulatory regions outside genes. Alternatively, the observed high population recombination rate (4Ner∼0.03), with overall linkage disequilibrium dropping rapidly within the order of few 100 bp, may compromise the power to detect causal loci with nearby markers. Our results add to the debate as to the utility of candidate gene approaches in relation to genomic features and the genetic architecture of the phenotypic trait in question.

Conservation of neutral substitution rate and substitutional asymmetries in mammalian genes.

Local variation in neutral substitution rate across mammalian genomes is governed by several factors, including sequence context variables and structural variables. In addition, the interplay of replication and transcription, known to induce a strand bias in mutation rate, gives rise to variation in substitutional strand asymmetries. Here, we address the conservation of variation in mutation rate and substitutional strand asymmetries using primate- and rodent-specific repeat elements located within the introns of protein-coding genes. We find significant but weak conservation of local mutation rates between human and mouse orthologs. Likewise, substitutional strand asymmetries are conserved between human and mouse, where substitution rate asymmetries show a higher degree of conservation than mutation rate. Moreover, we provide evidence that replication and transcription are correlated to the strength of substitutional asymmetries. The effect of transcription is particularly visible for genes with highly conserved gene expression. In comparison with replication and transcription, mutation rate influences the strength of substitutional asymmetries only marginally.

Segregation distortion in chicken and the evolutionary consequences of female meiotic drive in birds.

As all four meiotic products give rise to sperm in males, female meiosis result in a single egg in most eukaryotes. Any genetic element with the potential to influence chromosome segregation, so that it is preferentially included in the egg, should therefore gain a transmission advantage; a process termed female meiotic drive. We are aware of two chromosomal components, centromeres and telomeres, which share the potential to influence chromosome movement during meioses and make the following predictions based on the presence of female meiotic drive: (1) centromere-binding proteins should experience rapid evolution as a result of a conflict between driving centromeres and the rest of the genome; and (2) segregation patterns should be skewed near centromeres and telomeres. To test these predictions, we first analyze the molecular evolution of seven centromere-binding proteins in nine divergent bird species. We find strong evidence for positive selection in two genes, lending support to the genomic conflict hypothesis. Then, to directly test for the presence of segregation distortion, we also investigate the transmission of approximately 9000 single-nucleotide polymorphisms in 197 chicken families. By simulating fair Mendelian meioses, we locate chromosomal regions with statistically significant transmission ratio distortion. One region is located near the centromere on chromosome 1 and a second region is located near the telomere on the p-arm of chromosome 1. Although these observations do not provide conclusive evidence in favour of the meiotic drive/genome conflict hypothesis, they do lend support to the hypothesis that centromeres and telomeres drive during female meioses in chicken.

Avian genome evolution: insights from a linkage map of the blue tit (Cyanistes caeruleus).

We provide a first-generation linkage map of the blue tit (Cyanistes caeruleus), a passerine within the previously genetically uncharacterized family Paridae, which includes 91 orthologous loci with a single anchored position in the chicken (Gallus gallus) sequence assembly. The map consists of 18 linkage groups and covers 935 cM. There was highly conserved synteny between blue tit and chicken with the exception of a split on chromosome 1, potential splits on chromosome 4 and the translocation of two markers from chromosome 2 and 3, respectively, to chromosome 5. Gene order was very well conserved for the majority of chromosomes, an exception being chromosome 1 where multiple rearrangements were detected. Similar results were obtained in a comparison to the zebra finch (Taeniopygia guttata) genome assembly. The recombination rate in females was slightly higher than in males, implying a moderate degree of heterochiasmy in the blue tit. The map distance of the blue tit was approximately 78% of that of the Wageningen chicken broiler population, and very similar to the Uppsala chicken mapping population, over homologous genome regions. Apart from providing insights into avian recombination and genome evolution, our blue tit linkage map forms a valuable genetic resource for ecological and evolutionary research in Paridae.

All dosage compensation is local: gene-by-gene regulation of sex-biased expression on the chicken Z chromosome.

Recent reports have suggested that birds lack a mechanism of wholesale dosage compensation for the Z sex chromosome. This discovery was rather unexpected, as all other animals investigated with chromosomal mechanisms of sex determination have some method to counteract the effects of gene dosage of the dominant sex chromosome in males and females. Despite the lack of a global mechanism of avian dosage compensation, the pattern of gene expression difference between males and females varies a great deal for individual Z-linked genes. This suggests that some genes may be individually dosage compensated, and that some less-than-global pattern of dosage compensation, such as local or temporal, exists on the avian Z chromosome. We used global gene expression profiling in males and females for both somatic and gonadal tissue at several time points in the life cycle of the chicken to assess the pattern of sex-biased gene expression on the Z chromosome. Average fold-change between males and females varied somewhat among tissue time-point combinations, with embryonic brain samples having the smallest gene dosage effects, and adult gonadal tissue having the largest degree of male bias. Overall, there were no neighborhoods of overall dosage compensation along the Z. Taken together, this suggests that dosage compensation is regulated on the Z chromosome entirely on a gene-by-gene level, and can vary during the life cycle and by tissue type. This regulation may be an indication of how critical a given gene's functionality is, as the expression level for essential genes will be tightly regulated in order to avoid perturbing important pathways and networks with differential expression levels in males and females.

Molecular evolutionary genomics of birds.

Insight into the molecular evolution of birds has been offered by the steady accumulation of avian DNA sequence data, recently culminating in the first draft sequence of an avian genome, that of chicken. By studying avian molecular evolution we can learn about adaptations and phenotypic evolution in birds, and also gain an understanding of the similarities and differences between mammalian and avian genomes. In both these lineages, there is pronounced isochore structure with highly variable GC content. However, while mammalian isochores are decaying, they are maintained in the chicken lineage, which is consistent with a biased gene conversion model where the high and variable recombination rate of birds reinforces heterogeneity in GC. In Galliformes, GC is positively correlated with the rate of nucleotide substitution; the mean neutral mutation rate is 0.12-0.15% at each site per million years but this estimate comes with significant local variation in the rate of mutation. Comparative genomics reveals lower d(N)/d(S) ratios on micro- compared to macrochromosomes, possibly due to population genetic effects or a non-random distribution of genes with respect to chromosome size. A non-random genomic distribution is shown by genes with sex-biased expression, with male-biased genes over-represented and female-biased genes under-represented on the Z chromosome. A strong effect of selection is evident on the non-recombining W chromosome with high d(N)/d(S) ratios and limited polymorphism. Nucleotide diversity in chicken is estimated at 4-5 x 10(-3) which might be seen as surprisingly high given presumed bottlenecks during domestication, but is lower than that recently observed in several natural populations of other species. Several important aspects of the molecular evolutionary process of birds remain to be understood and it can be anticipated that the upcoming genome sequence of a second bird species, the zebra finch, as well as the integration of data on gene expression, shall further advance our knowledge of avian evolution.

Unequal contribution of sexes in the origin of dog breeds.

Dogs (Canis familiaris) were domesticated from the gray wolf (Canis lupus) at least 14,000 years ago, and there is evidence of dogs with phenotypes similar to those in modern breeds 4000 years ago. However, recent genetic analyses have suggested that modern dog breeds have a much more recent origin, probably <200 years ago. To study the origin of contemporaneous breeds we combined the analysis of paternally inherited Y chromosome markers with maternally inherited mitochondrial DNA and biparentally inherited autosomal microsatellite markers in both domestic dogs and their wild ancestor, the gray wolf. Our results show a sex bias in the origin of breeds, with fewer males than females contributing genetically, which clearly differs from the breeding patterns in wild gray wolf populations where both sexes have similar contributions. Furthermore, a comparison of mitochondrial DNA and Y chromosome diversity in dog groups recognized by the World Canine Organization, as well as in groups defined by the breeds' genetic composition, shows that paternal lineages are more differentiated among groups than maternal lineages. This demonstrates a lower exchange of males than of females between breeds belonging to different groups, which illustrates how breed founders may have been chosen.

Understanding the phylogeographic patterns of European hedgehogs, Erinaceus concolor and E. europaeus using the MHC.

The genome of the European hedgehog, Erinaceus concolor and E. europaeus, shows a strong signal of cycles of restriction to glacial refugia and postglacial expansion. Patterns of expansion, however, differ for mitochondrial DNA (mtDNA) and preliminary analysis of nuclear markers. In this study, we determine phylogeographic patterns in the hedgehog using two loci of the major histocompatibility complex (MHC), isolated for the first time in hedgehogs. These genes show long persistence times and high polymorphism in many species because of the actions of balancing selection. Among 84 individuals screened for variation, only two DQA alleles were identified in each species, but 10 DQB alleles were found in E. concolor and six in E. europaeus. A strong effect of demography on patterns of DQB variability is observed, with only weak evidence of balancing selection. While data from mtDNA clearly subdivide both species into monophyletic subgroups, the MHC data delineate only E. concolor into distinct subgroups, supporting the preliminary findings of other nuclear markers. Together with differences in variability, this suggests that the refugia history and/or expansion patterns of E. concolor and E. europaeus differ.

SNPs in ecological and conservation studies: a test in the Scandinavian wolf population.

Single nucleotide polymorphisms (SNPs) have the potential to become the genetic marker of choice in studies of the ecology and conservation of natural populations because of their capacity to access variability across the genome. In this study, we provide one of the first demonstrations of SNP discovery in a wild population in order to address typical issues of importance in ecology and conservation in the recolonized Scandinavian and neighbouring Finnish wolf Canis lupus populations. Using end sequence from BAC (bacterial artificial chromosome) clones specific for dogs, we designed assays for 24 SNP loci, 20 sites of which had previously been shown to be polymorphic in domestic dogs and four sites were newly identified as polymorphic in wolves. Of the 24 assayed loci, 22 SNPs were found to be variable within the Scandinavian population and, importantly, these were able to distinguish individual wolves from one another (unbiased probability of identity of 4.33 x 10(-8)), providing equivalent results to that derived from 12 variable microsatellites genotyped in the same population. An assignment test shows differentiation between the Scandinavian and neighbouring Finnish wolf populations, although not all known immigrants are accurately identified. An exploration of the misclassification rates in the identification of relationships shows that neither 22 SNP nor 20 microsatellite loci are able to discriminate across single order relationships. Despite the remaining obstacle of SNP discovery in nonmodel organisms, the use of SNPs in ecological and conservation studies is encouraged by the advent of large scale screening methods. Furthermore, the ability to amplify extremely small fragments makes SNPs of particular use for population monitoring, where faecal and other noninvasive samples are routinely used.

A temporal analysis shows major histocompatibility complex loci in the Scandinavian wolf population are consistent with neutral evolution.

The major histocompatibility complex (MHC) has an integral role in the immune system, and hence diversity at its genes may be of particular importance for the health of populations. In large populations, balancing selection maintains diversity in MHC genes, but theoretical expectations indicate that this form of selection is absent or inefficient in small populations. We examine the level of diversity at three MHC class II loci in the wolf population of Scandinavia, a population naturally recolonized with a genetic contribution from as few as three founders, and in four neighbouring wolf populations. In the Scandinavian wolf population, two alleles were found for each locus and the distribution of alleles is compatible with their linkage into two haplotypes. Changes in the level of heterozygosity over time since recolonization demonstrate the effects of the proposed arrival of an immigrant wolf. The maintenance of diversity is shown to be compatible with a neutral, random allocation of alleles, in conjunction with crossing between packs. A total of 15 DRB1, seven DQA and 10 DQB1 alleles are found in four neighbouring wolf populations, with substantial sharing across populations. Even in these larger populations, bottlenecks and fragmentation with consequent genetic drift are likely to have resulted in few indicators for balancing selection and significant differentiation of populations.

Sexual variation in heritability and genetic correlations of morphological traits in house sparrow (Passer domesticus).

Estimates of genetic components are important for our understanding of how individual characteristics are transferred between generations. We show that the level of heritability varies between 0.12 and 0.68 in six morphological traits in house sparrows (Passer domesticus L.) in northern Norway. Positive and negative genetic correlations were present among traits, suggesting evolutionary constraints on the evolution of some of these characters. A sexual difference in the amount of heritable genetic variation was found in tarsus length, wing length, bill depth and body condition index, with generally higher heritability in females. In addition, the structure of the genetic variance-covariance matrix for the traits differed between the sexes. Genetic correlations between males and females for the morphological traits were however large and not significantly different from one, indicating that sex-specific responses to selection will be influenced by intersexual differences in selection differentials. Despite this, some traits had heritability above 0.1 in females, even after conditioning on the additive genetic covariance between sexes and the additive genetic variances in males. Moreover, a meta-analysis indicated that higher heritability in females than in males may be common in birds. Thus, this indicates sexual differences in the genetic architecture of birds. Consequently, as in house sparrows, the evolutionary responses to selection will often be larger in females than males. Hence, our results suggest that sex-specific additive genetic variances and covariances, although ignored in most studies, should be included when making predictions of evolutionary changes from standard quantitative genetic models.

Cryptic population structure in a large, mobile mammalian predator: the Scandinavian lynx.

The Eurasian lynx (Lynx lynx) is an example of a species that has gone through a severe bottleneck, leading to near extinction in Scandinavia around 1930-- a pattern shared with several other large carnivorous mammals. Here we extend previous genetic analyses of northern European lynx, confirming that lynx from the Scandinavian Peninsula represent a distinct clade differing clearly from European conspecifics. Furthermore, and despite a recent bottleneck and subsequent range expansion, we detect marked genetic differentiation within Scandinavia. This differentiation is largely manifested as a north-south gradient, with a linear increase in the quantity FST/(1 - FST). Aided by computer simulations we find that this pattern is unlikely to have arisen by random genetic drift in the short time since lynx started to expand in the 1950s, suggesting that the spatial structure may predate the bottleneck. Individual-based analyses indicate that, instead of a continuous gradient, Scandinavian lynx may be structured into three more or less distinct groups, possibly corresponding to northern, central and southern subpopulations. The presence of such structuring was unknown previously and was unexpected from general considerations on the mobility of the species, historical data and the absence of geographical barriers. Our study demonstrates how molecular markers may be used to detect cryptic population structure, invisible using traditional methods.

Two centuries of the Scandinavian wolf population: patterns of genetic variability and migration during an era of dramatic decline.

The grey wolf (Canis lupus) was numerous on the Scandinavian peninsula in the early 19th century. However, as a result of intense persecution, the population declined dramatically and was virtually extinct from the peninsula by the 1960s. We examined historical patterns of genetic variability throughout the period of decline, from 1829 to 1979. Contemporary Finnish wolves, considered to be representative of a large eastern wolf population, were used for comparison. Mitochondrial DNA (mtDNA) variability among historical Scandinavian wolves was significantly lower than in Finland while Y chromosome variability was comparable between the two populations. This may suggest that long-distance migration from the east has been male-biased. Importantly though, as the historical population was significantly differentiated from contemporary Finnish wolves, the overall immigration rate to the Scandinavian peninsula appears to have been low. Levels of variability at autosomal microsatellite loci were high by the early 1800s but declined considerably towards the mid-20th century. At this time, approximately 40% of the allelic diversity and 30% of the heterozygosity had been lost. After 1940, however, there is evidence of several immigration events, coinciding with episodes of marked population increase in Russian Karelia and subsequent westwards migration.

Combined use of maternal, paternal and bi-parental genetic markers for the identification of wolf-dog hybrids.

The identification of hybrids is often a subject of primary concern for the development of conservation and management strategies, but can be difficult when the hybridizing species are closely related and do not possess diagnostic genetic markers. However, the combined use of mitochondrial DNA (mtDNA), autosomal and Y chromosome genetic markers may allow the identification of hybrids and of the direction of hybridization. We used these three types of markers to genetically characterize one possible wolf-dog hybrid in the endangered Scandinavian wolf population. We first characterized the variability of mtDNA and Y chromosome markers in Scandinavian wolves as well as in neighboring wolf populations and in dogs. While the mtDNA data suggested that the target sample could correspond to a wolf, its Y chromosome type had not been observed before in Scandinavian wolves. We compared the genotype of the target sample at 18 autosomal microsatellite markers with those expected in pure specimens and in hybrids using assignment tests. The combined results led to the conclusion that the animal was a hybrid between a Scandinavian female wolf and a male dog. This finding confirms that inter-specific hybridization between wolves and dogs can occur in natural wolf populations. A possible correlation between hybridization and wolf population density and disturbance deserves further research.

Contrasting levels of nucleotide diversity on the avian Z and W sex chromosomes.

Sex chromosomes may provide a context for studying the local effects of mutation rate on molecular evolution, since the two types of sex chromosomes are generally exposed to different mutational environments in male and female germ lines. Importantly, recent studies of some vertebrates have provided evidence for a higher mutation rate among males than among females. Thus, in birds, the Z chromosome, which spends two thirds of its time in the male germ line, is exposed to more mutations than the female-specific W chromosome. We show here that levels of nucleotide diversity are drastically higher on the avian Z chromosome than in paralogous sequences on the W chromosome. In fact, no intraspecific polymorphism whatsoever was seen in about 3.4 kb of CHD1W intron sequence from a total of >150 W chromosome copies of seven different bird species. In contrast, the amount of genetic variability in paralogous sequences on the Z chromosome was significant, with an average pairwise nucleotide diversity (d) of 0.0020 between CHD1Z introns and with 37 segregating sites in a total of 3.8 kb of Z sequence. The contrasting levels of genetic variability on the avian sex chromosomes are thus in a direction predicted from a male-biased mutation rate. However, although a low gene number, as well as some other factors, argues against background selection and/or selective sweeps shaping the genetic variability of the avian W chromosome, we cannot completely exclude selection as a contributor to the low levels of variation on the W chromosome.

Evolutionary genetics. Clonal inheritance of avian mitochondrial DNA.

We have taken a new approach to test the commonly accepted, but recently questioned, principle of clonal inheritance of vertebrate mitochondrial DNA (mtDNA) by relating its inheritance to a female-specific marker of nuclear DNA. Whereas this is impossible in organisms with male heterogamy (such as mammals), we show here that genealogies of mtDNA and the female-specific W chromosome of a bird species are completely concordant. Our results indicate that inheritance of mtDNA is free of detectable recombination effects over an evolutionary timescale.