Linkage mapping reveals sex-dimorphic map distances in a passerine bird

Male-biased recombination at chromosome ends in a songbird revealed by precisely mapping crossover positions

Recombination plays a crucial role in evolution by generating novel haplotypes and disrupting linkage between genes, thereby enhancing the efficiency of selection. Here, we analyze the genomes of 12 great reed warblers (Acrocephalus arundinaceus) in a 3-generation pedigree to identify precise crossover positions along the chromosomes. We located more than 200 crossovers and found that these were highly concentrated toward the telomeric ends of the chromosomes. Apart from this major pattern in the recombination landscape, we found significantly higher frequencies of crossovers in genic compared with intergenic regions, and in exons compared with introns. Moreover, while the number of recombination events was similar between the sexes, the crossovers were located significantly closer to the ends of paternal compared with maternal chromosomes. In conclusion, our study of the great reed warbler revealed substantial variation in crossover frequencies within chromosomes, with a distinct bias toward the sub-telomeric regions, particularly on the paternal side. These findings emphasize the importance of thoroughly screening the entire length of chromosomes to characterize the recombination landscape and uncover potential sex-biases in recombination.

The Genetic Architecture of Recombination Rates is Polygenic and Differs Between the Sexes in Wild House Sparrows (Passer domesticus)

Meiotic recombination through chromosomal crossing-over is a fundamental feature of sex and an important driver of genomic diversity. It ensures proper disjunction, allows increased selection responses, and prevents mutation accumulation; however, it is also mutagenic and can break up favorable haplotypes. This cost-benefit dynamic is likely to vary depending on mechanistic and evolutionary contexts, and indeed, recombination rates show huge variation in nature. Identifying the genetic architecture of this variation is key to understanding its causes and consequences. Here, we investigate individual recombination rate variation in wild house sparrows (Passer domesticus). We integrate genomic and pedigree data to identify autosomal crossover counts (ACCs) and intrachromosomal allelic shuffling (r¯intra) in 13,056 gametes transmitted from 2,653 individuals to their offspring. Females had 1.37 times higher ACC, and 1.55 times higher r¯intra than males. ACC and r¯intra were heritable in females and males (ACC h2 = 0.23 and 0.11; r¯intra  h2 = 0.12 and 0.14), but cross-sex additive genetic correlations were low (rA = 0.29 and 0.32 for ACC and r¯intra). Conditional bivariate analyses showed that all measures remained heritable after accounting for genetic values in the opposite sex, indicating that sex-specific ACC and r¯intra can evolve somewhat independently. Genome-wide models showed that ACC and r¯intra are polygenic and driven by many small-effect loci, many of which are likely to act in trans as global recombination modifiers. Our findings show that recombination rates of females and males can have different evolutionary potential in wild birds, providing a compelling mechanism for the evolution of sexual dimorphism in recombination.

Contrasting mtDNA and microsatellite data of great reed warbler Acrocephalus arundinaceus breeding populations on a small geographic scale

The great reed warbler has two genetically distinguishable haplogroups: "Clade A" occurs in higher proportions in Western Europe and Kazakhstan, and colonised Europe and Asia from a refugium in South-West Europe; and "Clade B", which is more common in Eastern Europe, and colonised parts of Europe from a refugium in the Middle East. Our aims were (i) to analyse the rate of differentiation in Hungarian breeding populations in order to see whether European-scale pattern is visible or not on as a small scale as the territory of Hungary and (ii) to compare the results obtained with mtDNA sequencing and microsatellite markers. To analyse the genetic differentiation, the mtDNA control region II was sequenced in 68 adult breeding birds, and 51 were fingerprinted at 11 microsatellite loci, while both analyses were performed on 36 birds (a total of 83 birds). The microsatellite data gave a better resolution and represented the fine-scale pattern of the suspected recolonisation. The lack of genetic differentiation among the breeding populations based on mitochondrial data seems to support this finding, because the admixture of the clades in this particular geographic region obliterates differentiation. Accordingly, the F_st values from different branches are significantly based on microsatellite data only. The mtDNA methods only give reliable results when a geographic and ecological factor plays a role in the population subdivision, but in the case of an intermixing population larger-scale studies are needed.

Extreme variation in recombination rate and genetic diversity along the Sylvioidea neo-sex chromosome

Recombination strongly impacts sequence evolution by affecting the extent of linkage and the efficiency of selection. Here, we study recombination over the Z chromosome in great reed warblers (Acrocephalus arundinaceus) using pedigree-based linkage mapping. This species has extended Z and W chromosomes ("neo-sex chromosomes") formed by a fusion between a part of chromosome 4A and the ancestral sex chromosomes, which provides a unique opportunity to assess recombination and sequence evolution in sex-linked regions of different ages. We assembled an 87.54 Mbp and 90.19 cM large Z with a small pseudoautosomal region (0.89 Mbp) at one end and the fused Chr4A-part at the other end of the chromosome. A prominent feature in our data was an extreme variation in male recombination rate along Z with high values at both chromosome ends, but an apparent lack of recombination over a substantial central section, covering 78% of the chromosome. The nonrecombining region showed a drastic loss of genetic diversity and accumulation of repeats compared to the recombining parts. Thus, our data emphasize a key role of recombination in affecting local levels of polymorphism. Nonetheless, the evolutionary rate of genes (dN/dS) did not differ between high and low recombining regions, suggesting that the efficiency of selection on protein-coding sequences can be maintained also at very low levels of recombination. Finally, the Chr4A-derived part showed a similar recombination rate as the part of the ancestral Z that did recombine, but its sequence characteristics reflected both its previous autosomal, and current Z-linked, recombination patterns.

MHCtools - an R package for MHC high-throughput sequencing data: genotyping, haplotype and supertype inference, and downstream genetic analyses in non-model organisms

The major histocompatibility complex (MHC) plays a central role in the vertebrate adaptive immune system and has been of long-term interest in evolutionary biology. While several protocols have been developed for MHC genotyping, there is a lack of transparent and standardized tools for downstream analysis of MHC data. Here, we present the R package MHCtools and demonstrate the use of its functions to (i) assist accurate MHC genotyping from high-throughput amplicon-sequencing data, (ii) infer functional MHC supertypes using bootstrapped clustering analysis, (iii) identify segregating MHC haplotypes from family data, and (iv) analyse functional and genetic distances between MHC sequences. We employed MHCtools to analyse MHC class I (MHC-I) amplicon data of 559 great reed warblers (Acrocephalus arundinaceus). We identified 390 MHC-I alleles which clustered into 14 functional supertypes. A phylogenetic analysis and analyses of positive selection suggested that the MHC-I alleles belonged to several distinct functional groups. We furthermore identified 107 segregating haplotypes among 116 families, and found substantial variation in diversity with 4-21 MHC-I alleles and 3-13 MHC-I supertypes per haplotype. Finally, we show that the great reed warbler haplotypes harboured combinations of MHC-I supertypes with greater functional divergence than observed in simulated populations of possible haplotypes, a result that is in accordance with the divergent allele advantage hypothesis. Our study demonstrates the power of MHCtools to support genotyping and analysis of MHC in non-model species, which we hope will encourage broad implementation among researchers in MHC genetics and evolution.

A High-Density SNP Genetic Linkage Map and QTL Analysis of Growth-Related Traits in a Hybrid Family of Oysters (Crassostrea gigas × Crassostrea angulata) Using Genotyping-by-Sequencing

Oysters are among the most important species in global aquaculture. Crassostrea gigas, and its subspecies C. angulata, are the major cultured species. To determine the genetic basis of growth-related traits in oysters, we constructed a second-generation linkage map from 3367 single-nucleotide polymorphisms (SNPs) based on genotyping-by-sequencing, genotyped from a C. gigas × C. angulata hybrid family. These 3367 SNPs were distributed on 1695 markers, which were assigned to 10 linkage groups. The genetic linkage map had a total length of 1084.3 cM, with an average of 0.8 cM between markers; it thus represents the densest genetic map constructed for oysters to date. Twenty-seven quantitative trait loci (QTL) for five growth-related traits were detected. These QTL could explain 4.2-7.7% (mean = 5.4%) of the phenotypic variation. In total, 50.8% of phenotypic variance for shell width, 7.7% for mass weight, and 34.1% for soft tissue weight were explained. The detected QTL were distributed among eight linkage groups, and more than half (16) were concentrated within narrow regions in their respective linkage groups. Thirty-eight annotated genes were identified within the QTL regions, two of which are key genes for carbohydrate metabolism. Other genes were found to participate in assembly and regulation of the actin cytoskeleton, signal transduction, and regulation of cell differentiation and development. The newly developed high-density genetic map, and the QTL and candidate genes identified provide a valuable genetic resource and a basis for marker-assisted selection for C. gigas and C. angulata.

High-Resolution Mapping of Crossover and Non-crossover Recombination Events by Whole-Genome Re-sequencing of an Avian Pedigree

Recombination is an engine of genetic diversity and therefore constitutes a key process in evolutionary biology and genetics. While the outcome of crossover recombination can readily be detected as shuffled alleles by following the inheritance of markers in pedigreed families, the more precise location of both crossover and non-crossover recombination events has been difficult to pinpoint. As a consequence, we lack a detailed portrait of the recombination landscape for most organisms and knowledge on how this landscape impacts on sequence evolution at a local scale. To localize recombination events with high resolution in an avian system, we performed whole-genome re-sequencing at high coverage of a complete three-generation collared flycatcher pedigree. We identified 325 crossovers at a median resolution of 1.4 kb, with 86% of the events localized to <10 kb intervals. Observed crossover rates were in excellent agreement with data from linkage mapping, were 52% higher in male (3.56 cM/Mb) than in female meiosis (2.28 cM/Mb), and increased towards chromosome ends in male but not female meiosis. Crossover events were non-randomly distributed in the genome with several distinct hot-spots and a concentration to genic regions, with the highest density in promoters and CpG islands. We further identified 267 non-crossovers, whose location was significantly associated with crossover locations. We detected a significant transmission bias (0.18) in favour of 'strong' (G, C) over 'weak' (A, T) alleles at non-crossover events, providing direct evidence for the process of GC-biased gene conversion in an avian system. The approach taken in this study should be applicable to any species and would thereby help to provide a more comprehensive portray of the recombination landscape across organism groups.

A microsatellite-based linkage map for song sparrows (Melospiza melodia)

Although linkage maps are important tools in evolutionary biology, their availability for wild populations is limited. The population of song sparrows (Melospiza melodia) on Mandarte Island, Canada, is among the more intensively studied wild animal populations. Its long-term pedigree data, together with extensive genetic sampling, have allowed the study of a range of questions in evolutionary biology and ecology. However, the availability of genetic markers has been limited. We here describe 191 new microsatellite loci, including 160 high-quality polymorphic autosomal, 7 Z-linked and 1 W-linked markers. We used these markers to construct a linkage map for song sparrows with a total sex-averaged map length of 1731 cM and covering 35 linkage groups, and hence, these markers cover most of the 38-40 chromosomes. Female and male map lengths did not differ significantly. We then bioinformatically mapped these loci to the zebra finch (Taeniopygia guttata) genome and found that linkage groups were conserved between song sparrows and zebra finches. Compared to the zebra finch, marker order within small linkage groups was well conserved, whereas the larger linkage groups showed some intrachromosomal rearrangements. Finally, we show that as expected, recombination frequency between linked loci explained the majority of variation in gametic phase disequilibrium. Yet, there was substantial overlap in gametic phase disequilibrium between pairs of linked and unlinked loci. Given that the microsatellites described here lie on 35 of the 38-40 chromosomes, these markers will be useful for studies in this species, as well as for comparative genomics studies with other species.

Sex-specific recombination maps for individual macrochromosomes in the Japanese quail (Coturnix japonica)

Meiotic recombination in the Japanese quail was directly studied by immunolocalization of mutL homolog 1 (MLH1), a mismatch repair protein of mature recombination nodules. In total, 15,862 crossovers were scored along the autosomal synaptonemal complexes in 308 meiotic nuclei from males and females. Crossover frequencies calculated from MLH1 foci show wide similitude between males and females with slightly higher number of foci in females. From this analysis, we predict that the sex-averaged map length of the Japanese quail is 2580 cM, with a genome-wide recombination rate of 1.9 cM/Mb. MLH1 focus mapping along the six largest bivalents showed few intersex differences in the distribution of crossovers along with variant patterns in metacentric and acrocentric macrobivalents. These results provide valuable information to complement linkage map analysis in the species while providing insight into our understanding of the mechanisms of crossover distribution along chromosome arms.

Strong artificial selection in domestic mammals did not result in an increased recombination rate

Recombination rates vary in intensity and location at the species, individual, sex and chromosome levels. Despite the fundamental biological importance of this process, the selective forces that operate to shape recombination rate and patterns are unclear. Domestication offers a unique opportunity to study the interplay between recombination and selection. In domesticates, intense selection for particular traits is imposed on small populations over many generations, resulting in organisms that differ, sometimes dramatically, in morphology and physiology from their wild ancestor. Although earlier studies suggested increased recombination rate in domesticates, a formal comparison of recombination rates between domestic mammals and their wild congeners was missing. In order to determine broad-scale recombination rate, we used immunolabeling detection of MLH1 foci as crossover markers in spermatocytes in three pairs of closely related wild and domestic species (dog and wolf, goat and ibex, and sheep and mouflon). In the three pairs, and contrary to previous suggestions, our data show that contemporary recombination rate is higher in the wild species. Subsequently, we inferred recombination breakpoints in sequence data for 16 genomic regions in dogs and wolves, each containing a locus associated with a dog phenotype potentially under selection during domestication. No difference in the number and distribution of recombination breakpoints was found between dogs and wolves. We conclude that our data indicate that strong directional selection did not result in changes in recombination in domestic mammals, and that both upper and lower bounds for crossover rates may be tightly regulated.

A high-density linkage map enables a second-generation collared flycatcher genome assembly and reveals the patterns of avian recombination rate variation and chromosomal evolution

Detailed linkage and recombination rate maps are necessary to use the full potential of genome sequencing and population genomic analyses. We used a custom collared flycatcher 50 K SNP array to develop a high-density linkage map with 37 262 markers assigned to 34 linkage groups in 33 autosomes and the Z chromosome. The best-order map contained 4215 markers, with a total distance of 3132 cm and a mean genetic distance between markers of 0.12 cm. Facilitated by the array being designed to include markers from most scaffolds, we obtained a second-generation assembly of the flycatcher genome that approaches full chromosome sequences (N50 super-scaffold size 20.2 Mb and with 1.042 Gb (of 1.116 Gb) anchored to and mostly ordered and oriented along chromosomes). We found that flycatcher and zebra finch chromosomes are entirely syntenic but that inversions at mean rates of 1.5–2.0 event (6.6–7.5 Mb) per My have changed the organization within chromosomes, rates high enough for inversions to potentially have been involved with many speciation events during avian evolution. The mean recombination rate was 3.1 cm/Mb and correlated closely with chromosome size, from 2 cm/Mb for chromosomes >100 Mb to >10 cm/Mb for chromosomes <10 Mb. This size dependence seemed entirely due to an obligate recombination event per chromosome; if 50 cm was subtracted from the genetic lengths of chromosomes, the rate per physical unit DNA was constant across chromosomes. Flycatcher recombination rate showed similar variation along chromosomes as chicken but lacked the large interior recombination deserts characteristic of zebra finch chromosomes.

Host genotype and age have no effect on rejection of parasitic eggs

Egg rejection belongs to a widely used host tactic to prevent the costs incurred by avian brood parasitism. However, the genetic basis of this behaviour and the effect of host age on the probability of rejecting the parasitic egg remain largely unknown. Here, we used a set of 15 polymorphic microsatellite loci, including a previously detected candidate locus (Ase64), to link genotypes of female great reed warblers (Acrocephalus arundinaceus), a known rejecter, with their egg rejection responses in two host populations. We also tested whether host female age, as a measure of the experience with own eggs, plays a role in rejection of common cuckoo (Cuculus canorus) eggs. We failed to find any consistent association of egg rejection responses with host female genotypes or age. It seems that host decisions on egg rejection show high levels of phenotypic plasticity and are likely to depend on the spatiotemporal variation in the parasitism pressure. Future studies exploring the repeatability of host responses towards parasitic eggs and the role of host individual experience with parasitic eggs would greatly improve our understanding of the variations in host behaviours considering the persistence of brood parasitism in host populations with rejecter phenotypes.

Replicated high-density genetic maps of two great tit populations reveal fine-scale genomic departures from sex-equal recombination rates

Linking variation in quantitative traits to variation in the genome is an important, but challenging task in the study of life-history evolution. Linkage maps provide a valuable tool for the unravelling of such trait-gene associations. Moreover, they give insight into recombination landscapes and between-species karyotype evolution. Here we used genotype data, generated from a 10k single-nucleotide polymorphism (SNP) chip, of over 2000 individuals to produce high-density linkage maps of the great tit (Parus major), a passerine bird that serves as a model species for ecological and evolutionary questions. We created independent maps from two distinct populations: a captive F2-cross from The Netherlands (NL) and a wild population from the United Kingdom (UK). The two maps contained 6554 SNPs in 32 linkage groups, spanning 2010 cM and 1917 cM for the NL and UK populations, respectively, and were similar in size and marker order. Subtle levels of heterochiasmy within and between chromosomes were remarkably consistent between the populations, suggesting that the local departures from sex-equal recombination rates have evolved. This key and surprising result would have been impossible to detect if only one population was mapped. A comparison with zebra finch Taeniopygia guttata, chicken Gallus gallus and the green anole lizard Anolis carolinensis genomes provided further insight into the evolution of avian karyotypes.

A high-density SNP genetic linkage map for the silver-lipped pearl oyster, Pinctada maxima: a valuable resource for gene localisation and marker-assisted selection

Background

The silver-lipped pearl oyster, Pinctada maxima, is an important tropical aquaculture species extensively farmed for the highly sought "South Sea" pearls. Traditional breeding programs have been initiated for this species in order to select for improved pearl quality, but many economic traits under selection are complex, polygenic and confounded with environmental factors, limiting the accuracy of selection. The incorporation of a marker-assisted selection (MAS) breeding approach would greatly benefit pearl breeding programs by allowing the direct selection of genes responsible for pearl quality. However, before MAS can be incorporated, substantial genomic resources such as genetic linkage maps need to be generated. The construction of a high-density genetic linkage map for P. maxima is not only essential for unravelling the genomic architecture of complex pearl quality traits, but also provides indispensable information on the genome structure of pearl oysters.

Results

A total of 1,189 informative genome-wide single nucleotide polymorphisms (SNPs) were incorporated into linkage map construction. The final linkage map consisted of 887 SNPs in 14 linkage groups, spans a total genetic distance of 831.7 centimorgans (cM), and covers an estimated 96% of the P. maxima genome. Assessment of sex-specific recombination across all linkage groups revealed limited overall heterochiasmy between the sexes (i.e. 1.15:1 F/M map length ratio). However, there were pronounced localised differences throughout the linkage groups, whereby male recombination was suppressed near the centromeres compared to female recombination, but inflated towards telomeric regions. Mean values of LD for adjacent SNP pairs suggest that a higher density of markers will be required for powerful genome-wide association studies. Finally, numerous nacre biomineralization genes were localised providing novel positional information for these genes.

Conclusions

This high-density SNP genetic map is the first comprehensive linkage map for any pearl oyster species. It provides an essential genomic tool facilitating studies investigating the genomic architecture of complex trait variation and identifying quantitative trait loci for economically important traits useful in genetic selection programs within the P. maxima pearling industry. Furthermore, this map provides a foundation for further research aiming to improve our understanding of the dynamic process of biomineralization, and pearl oyster evolution and synteny.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-810) contains supplementary material, which is available to authorized users.

Speciation: clash of the genomes

Complete genomes of hybridizing bird species demonstrate the importance of the sex chromosomes, telomeres and centromeres to the initial stages of speciation.

Variation in MHC genotypes in two populations of house sparrow (Passer domesticus) with different population histories

Small populations are likely to have a low genetic ability for disease resistance due to loss of genetic variation through inbreeding and genetic drift. In vertebrates, the highest genetic diversity of the immune system is located at genes within the major histocompatibility complex (MHC). Interestingly, parasite-mediated selection is thought to potentially maintain variation at MHC loci even in populations that are monomorphic at other loci. Therefore, general loss of genetic variation in the genome may not necessarily be associated with low variation at MHC loci. We evaluated inter- and intrapopulation variation in MHC genotypes between an inbred (Aldra) and a relatively outbred population (Hestmannøy) of house sparrows (Passer domesticus) in a metapopulation at Helgeland, Norway. Genomic (gDNA) and transcribed (cDNA) alleles of functional MHC class I and IIB loci, along with neutral noncoding microsatellite markers, were analyzed to obtain relevant estimates of genetic variation. We found lower allelic richness in microsatellites in the inbred population, but high genetic variation in MHC class I and IIB loci in both populations. This suggests that also the inbred population could be under balancing selection to maintain genetic variation for pathogen resistance.

Evidence of a neo-sex chromosome in birds

Neo-sex chromosomes often originate from sex chromosome-autosome fusions and constitute an important basis for the study of gene degeneration and expression in a sex chromosomal context. Neo-sex chromosomes are known from many animal and plant lineages, but have not been reported in birds, a group in which genome organization seems particularly stable. Following indications of sex linkage and unexpected sex-biased gene expression in warblers (Sylvioidea; Passeriformes), we have conducted an extensive marker analysis targeting 31 orthologues of loci on zebra finch chromosome 4a in five species, representative of independent branches of Passerida. We identified a region of sex linkage covering approximately the first half (10 Mb) of chromosome 4a, and associated to both Z and W chromosomes, in three Sylvioidea passerine species. Linkage analysis in an extended pedigree of one species additionally confirmed the association between this part of chromosome 4a and the Z chromosome. Markers located between 10 and 21 Mb of chromosome 4a showed no signs of sex linkage, suggesting that only half of the chromosome was involved in this transition. No sex linkage was observed in non-Sylvioidea passerines, indicating that the neo-sex chromosome arose at the base of the Sylvioidea branch of the avian phylogeny, at 47.4-37.6 millions years ago (MYA), substantially later than the ancestral sex chromosomes (150 MYA). We hypothesize that the gene content of chromosome 4a might be relevant in its transition to a sex chromosome, based on the presence of genes (for example, the androgen receptor) that could offer a selective advantage when associated to Z-linked sex determination loci.

Comparative BAC-based mapping in the white-throated sparrow, a novel behavioral genomics model, using interspecies overgo hybridization

Background

The genomics era has produced an arsenal of resources from sequenced organisms allowing researchers to target species that do not have comparable mapping and sequence information. These new "non-model" organisms offer unique opportunities to examine environmental effects on genomic patterns and processes. Here we use comparative mapping as a first step in characterizing the genome organization of a novel animal model, the white-throated sparrow (Zonotrichia albicollis), which occurs as white or tan morphs that exhibit alternative behaviors and physiology. Morph is determined by the presence or absence of a complex chromosomal rearrangement. This species is an ideal model for behavioral genomics because the association between genotype and phenotype is absolute, making it possible to identify the genomic bases of phenotypic variation.

Findings

We initiated a genomic study in this species by characterizing the white-throated sparrow BAC library via filter hybridization with overgo probes designed for the chicken, turkey, and zebra finch. Cross-species hybridization resulted in 640 positive sparrow BACs assigned to 77 chicken loci across almost all macro-and microchromosomes, with a focus on the chromosomes associated with morph. Out of 216 overgos, 36% of the probes hybridized successfully, with an average number of 3.0 positive sparrow BACs per overgo.

Conclusions

These data will be utilized for determining chromosomal architecture and for fine-scale mapping of candidate genes associated with phenotypic differences. Our research confirms the utility of interspecies hybridization for developing comparative maps in other non-model organisms.

Sex linkage, sex-specific selection, and the role of recombination in the evolution of sexually dimorphic gene expression

Sex-biased genes--genes that are differentially expressed within males and females--are nonrandomly distributed across animal genomes, with sex chromosomes and autosomes often carrying markedly different concentrations of male- and female-biased genes. These linkage patterns are often gene- and lineage-dependent, differing between functional genetic categories and between species. Although sex-specific selection is often hypothesized to shape the evolution of sex-linked and autosomal gene content, population genetics theory has yet to account for many of the gene- and lineage-specific idiosyncrasies emerging from the empirical literature. With the goal of improving the connection between evolutionary theory and a rapidly growing body of genome-wide empirical studies, we extend previous population genetics theory of sex-specific selection by developing and analyzing a biologically informed model that incorporates sex linkage, pleiotropy, recombination, and epistasis, factors that are likely to vary between genes and between species. Our results demonstrate that sex-specific selection and sex-specific recombination rates can generate, and are compatible with, the gene- and species-specific linkage patterns reported in the genomics literature. The theory suggests that sexual selection may strongly influence the architectures of animal genomes, as well as the chromosomal distribution of fixed substitutions underlying sexually dimorphic traits.

Genome-wide SNP detection in the great tit Parus major using high throughput sequencing

Identifying genes that underlie ecological traits will open exiting possibilities to study gene-environment interactions in shaping phenotypes and in measuring natural selection on genes. Evolutionary ecology has been pursuing these objectives for decades, but they come into reach now that next generation sequencing technologies have dramatically lowered the costs to obtain the genomic sequence information that is currently lacking for most ecologically important species. Here we describe how we generated over 2 billion basepairs of novel sequence information for an ecological model species, the great tit Parus major. We used over 16 million short sequence reads for the de novo assembly of a reference sequence consisting of 550 000 contigs, covering 2.5% of the genome of the great tit. This reference sequence was used as the scaffold for mapping of the sequence reads, which allowed for the detection of over 20 000 novel single nucleotide polymorphisms. Contigs harbouring 4272 of the single nucleotide polymorphisms could be mapped to a unique location on the recently sequenced zebra finch genome. Of all the great tit contigs, significantly more were mapped to the microchromosomes than to the intermediate and the macrochromosomes of the zebra finch, indicating a higher overall level of sequence conservation on the microchromosomes than on the other types of chromosomes. The large number of great tit contigs that can be aligned to the zebra finch genome shows that this genome provides a valuable framework for large scale genetics, e.g. QTL mapping or whole genome association studies, in passerines.

The recombination landscape of the zebra finch Taeniopygia guttata genome

Understanding the causes and consequences of variation in the rate of recombination is essential since this parameter is considered to affect levels of genetic variability, the efficacy of selection, and the design of association and linkage mapping studies. However, there is limited knowledge about the factors governing recombination rate variation. We genotyped 1920 single nucleotide polymorphisms in a multigeneration pedigree of more than 1000 zebra finches (Taeniopygia guttata) to develop a genetic linkage map, and then we used these map data together with the recently available draft genome sequence of the zebra finch to estimate recombination rates in 1 Mb intervals across the genome. The average zebra finch recombination rate (1.5 cM/Mb) is higher than in humans, but significantly lower than in chicken. The local rates of recombination in chicken and zebra finch were only weakly correlated, demonstrating evolutionary turnover of the recombination landscape in birds. The distribution of recombination events was heavily biased toward ends of chromosomes, with a stronger telomere effect than so far seen in any organism. In fact, the recombination rate was as low as 0.1 cM/Mb in intervals up to 100 Mb long in the middle of the larger chromosomes. We found a positive correlation between recombination rate and GC content, as well as GC-rich sequence motifs. Levels of linkage disequilibrium (LD) were significantly higher in regions of low recombination, showing that heterogeneity in recombination rates have left a footprint on the genomic landscape of LD in zebra finch populations.

A comparison of SNPs and microsatellites as linkage mapping markers: lessons from the zebra finch (Taeniopygia guttata)

BACKGROUND

Genetic linkage maps are essential tools when searching for quantitative trait loci (QTL). To maximize genome coverage and provide an evenly spaced marker distribution a combination of different types of genetic marker are sometimes used. In this study we created linkage maps of four zebra finch (Taeniopygia guttata) chromosomes (1, 1A, 2 and 9) using two types of marker, Single Nucleotide Polymorphisms (SNPs) and microsatellites. To assess the effectiveness and accuracy of each kind of marker we compared maps built with each marker type separately and with both types of marker combined. Linkage map marker order was validated by making comparisons to the assembled zebra finch genome sequence.

RESULTS

We showed that marker order was less reliable and linkage map lengths were inflated for microsatellite maps relative to SNP maps, apparently due to differing error rates between the two types of marker. Guidelines on how to minimise the effects of error are provided. In particular, we show that when combining both types of marker the conventional process of building linkage maps, whereby the most informative markers are added to the map first, has to be modified in order to improve map accuracy.

CONCLUSIONS

When using multiple types and large numbers of markers to create dense linkage maps, the least error prone loci (SNPs) rather than the most informative should be used to create framework maps before the addition of other potentially more error prone markers (microsatellites). This raises questions about the accuracy of marker order and predicted recombination rates in previous microsatellite linkage maps which were created using the conventional building process, however, provided suitable error detection strategies are followed microsatellite-based maps can continue to be regarded as reasonably reliable.

A strong quantitative trait locus for wing length on chromosome 2 in a wild population of great reed warblers

Wing length is a key character for essential behaviours related to bird flight such as migration and foraging. In the present study, we initiate the search for the genes underlying wing length in birds by studying a long-distance migrant, the great reed warbler (Acrocephalus arundinaceus). In this species wing length is an evolutionary interesting trait with pronounced latitudinal gradient and sex-specific selection regimes in local populations. We performed a quantitative trait locus (QTL) scan for wing length in great reed warblers using phenotypic, genotypic, pedigree and linkage map data from our long-term study population in Sweden. We applied the linkage analysis mapping method implemented in GridQTL (a new web-based software) and detected a genome-wide significant QTL for wing length on chromosome 2, to our knowledge, the first detected QTL in wild birds. The QTL extended over 25 cM and accounted for a substantial part (37%) of the phenotypic variance of the trait. A genome scan for tarsus length (a body-size-related trait) did not show any signal, implying that the wing-length QTL on chromosome 2 was not associated with body size. Our results provide a first important step into understanding the genetic architecture of avian wing length, and give opportunities to study the evolutionary dynamics of wing length at the locus level.

The evolution of heterochiasmy: the role of sexual selection and sperm competition in determining sex-specific recombination rates in eutherian mammals

Early karyotypic work revealed that female and male recombination rates in many species show pronounced differences, and this pattern of heterochiasmy has also been observed in modern linkage mapping studies. Several hypotheses to explain this phenomenon have been offered, ranging from strictly biological mechanisms related to the gametic differences between the sexes, to more evolutionary models based on sexually antagonistic selection. However, despite the long history of interest in heterochiasmy, empirical data has failed to support any theory or pattern consistently. Here I test two alternative evolutionary hypotheses regarding heterochiasmy across the eutherian mammals, and show that sexual dimorphism, but not sperm competition, is strongly correlated with recombination rate, suggesting that sexual antagonism is an important influence. However, the observed relationship between heterochiasmy and sexual dimorphism runs counter to theoretical predictions, with male recombination higher in species with high levels of sexual dimorphism. This may be the response to male-biased dispersal, which, rather than the static male fitness landscape envisioned in the models tested here, could radically shift optimal male fitness parameters among generations.

A genome-wide set of 106 microsatellite markers for the blue tit (Cyanistes caeruleus)

We have characterized a set of 106 microsatellite markers in 26-127 individual blue tits (Cyanistes caeruleus), and assigned their location on the zebra finch (Taeniopygia guttata) and on the chicken (Gallus gallus) genome on the basis of sequence homology. Thirty-one markers are newly designed from zebra finch EST (expressed sequence tags) sequences, 22 markers were developed by others from EST sequences using different methods and the remaining 53 loci were previously designed or modified passerine markers. The 106 microsatellite markers are distributed over 26 and 24 chromosomes in the zebra finch and in the chicken genome respectively and the number of alleles varies between 2 and 49. Eight loci deviate significantly from Hardy-Weinberg equilibrium and show a high frequency of null alleles, and three pairs of markers located in the same chromosome appear to be in linkage disequilibrium. With the exception of these few loci, the polymorphic microsatellite markers presented here provide a useful genome-wide resource for population and evolutionary genetic studies of the blue tit, in addition to their potential utility in other passerine birds.

Genomic linkage map of the human blood fluke Schistosoma mansoni

BACKGROUND

Schistosoma mansoni is a blood fluke that infects approximately 90 million people. The complete life cycle of this parasite can be maintained in the laboratory, making this one of the few experimentally tractable human helminth infections, and a rich literature reveals heritable variation in important biomedical traits such as virulence, host-specificity, transmission and drug resistance. However, there is a current lack of tools needed to study S. mansoni's molecular, quantitative, and population genetics. Our goal was to construct a genetic linkage map for S. mansoni, and thus provide a new resource that will help stimulate research on this neglected pathogen.

RESULTS

We genotyped grandparents, parents and 88 progeny to construct a 5.6 cM linkage map containing 243 microsatellites positioned on 203 of the largest scaffolds in the genome sequence. The map allows 70% of the estimated 300 Mb genome to be ordered on chromosomes, and highlights where scaffolds have been incorrectly assembled. The markers fall into eight main linkage groups, consistent with seven pairs of autosomes and one pair of sex chromosomes, and we were able to anchor linkage groups to chromosomes using fluorescent in situ hybridization. The genome measures 1,228.6 cM. Marker segregation reveals higher female recombination, confirms ZW inheritance patterns, and identifies recombination hotspots and regions of segregation distortion.

CONCLUSIONS

The genetic linkage map presented here is the first for S. mansoni and the first for a species in the phylum Platyhelminthes. The map provides the critical tool necessary for quantitative genetic analysis, aids genome assembly, and furnishes a framework for comparative flatworm genomics and field-based molecular epidemiological studies.

Characterization of 38 microsatellite loci in the European blackbird, Turdus merula (Turdidae, AVES)

We characterized 38 microsatellite loci in the European blackbird, Turdus merula. Thirty-seven loci were identified by testing 242 loci that had been originally isolated in other avian species. One additional locus was isolated from a European blackbird genomic library. All loci were characterized in 20-29 blackbirds from a population in the Czech Republic and displayed between two and 16 alleles, with observed heterozygosity ranging from 0.04 to 1.00. Thirty-seven loci could be assigned a chromosome location in the zebra finch (Taeniopygia guttata) genome based on sequence homology.

Search for informative polymorphisms in candidate genes: clock genes and circadian behaviour in blue tits

The identification of functional polymorphisms in genes that underlie behavioural trait variation is a challenging but intriguing task in evolutionary biology. Given the wealth of genomic data and the increasing number of genotype-phenotype association studies in model organisms, one can ask whether and how this information can be used for non-model organisms. Here we describe two strategies to search for likely functional polymorphisms in candidate genes in a bird species that has been intensively studied by behavioural and population ecologists, the blue tit Cyanistes caeruleus. In the first approach we searched for repeating elements in coding regions of the genome using information about repeats in Gallus gallus genes. The rationale is that tandem-repeat elements have a high potential to be polymorphic and functional. The second strategy aimed to replicate reported genotype-phenotype association studies by extrapolating results from model organisms to our study species. Both strategies showed high success rates with respect to finding homologous gene regions and potentially informative genetic variants in the genes AANAT, ADCYAP1, CKIepsilon, CLOCK, CREB1, NPAS2 and PERIOD2.

A first-generation microsatellite-based genetic linkage map of the Siberian jay (Perisoreus infaustus): insights into avian genome evolution

BACKGROUND

Genomic resources for the majority of free-living vertebrates of ecological and evolutionary importance are scarce. Therefore, linkage maps with high-density genome coverage are needed for progress in genomics of wild species. The Siberian jay (Perisoreus infaustus; Corvidae) is a passerine bird which has been subject to lots of research in the areas of ecology and evolutionary biology. Knowledge of its genome structure and organization is required to advance our understanding of the genetic basis of ecologically important traits in this species, as well as to provide insights into avian genome evolution.

RESULTS

We describe the first genetic linkage map of Siberian jay constructed using 117 microsatellites and a mapping pedigree of 349 animals representing five families from a natural population breeding in western Finland from the years 1975 to 2006. Markers were resolved into nine autosomal and a Z-chromosome-specific linkage group, 10 markers remaining unlinked. The best-position map with the most likely positions of all significantly linked loci had a total sex-average size of 862.8 cM, with an average interval distance of 9.69 cM. The female map covered 988.4 cM, whereas the male map covered only 774 cM. The Z-chromosome linkage group comprised six markers, three pseudoautosomal and three sex-specific loci, and spanned 10.6 cM in females and 48.9 cM in males. Eighty-one of the mapped loci could be ordered on a framework map with odds of >1000:1 covering a total size of 809.6 cM in females and 694.2 cM in males. Significant sex specific distortions towards reduced male recombination rates were revealed in the entire best-position map as well as within two autosomal linkage groups. Comparative mapping between Siberian jay and chicken anchored 22 homologous loci on 6 different linkage groups corresponding to chicken chromosomes Gga1, 2, 3, 4, 5, and Z. Quite a few cases of intra-chromosomal rearrangements within the autosomes and three cases of inter-chromosomal rearrangement between the Siberian jay autosomal linkage groups (LG1, LG2 and LG3) and the chicken sex chromosome GgaZ were observed, suggesting a conserved synteny, but changes in marker order, within autosomes during about 100 million years of avian evolution.

CONCLUSION

The constructed linkage map represents a valuable resource for intraspecific genomics of Siberian jay, as well as for avian comparative genomic studies. Apart from providing novel insights into sex-specific recombination rates and patterns, the described maps - from a previously genomically uncharacterized superfamily (Corvidae) of passerine birds - provide new insights into avian genome evolution. In combination with high-resolution data on quantitative trait variability from the study population, they also provide a foundation for QTL-mapping studies.

Enhanced cross-species utility of conserved microsatellite markers in shorebirds

Background

Microsatellite markers are popular genetic markers frequently used in forensic biology. Despite their popularity, the characterisation of polymorphic microsatellite loci and development of suitable markers takes considerable effort. Newly-available genomic databases make it feasible to identify conserved genetic markers. We examined the utility and characteristics of conserved microsatellite markers in Charadriiformes (plovers, sandpipers, gulls and auks). This order harbours many species with diverse breeding systems, life histories and extraordinary migration biology whose genetics warrant investigation. However, research has been largely restrained by the limited availability of genetic markers. To examine the utility of conserved microsatellite loci as genetic markers we collated a database of Charadriiformes microsatellites, searched for homologues in the chicken genome and tested conserved markers for amplification and polymorphism in a range of charadriiform species.

Results

Sixty-eight (42%) of 161 charadriiform microsatellite loci were assigned to a single location in the chicken genome based on their E-value. Fifty-five primers designed from conserved microsatellite loci with an E-value of E-10 or lower amplified across a wider range of charadriiform species than a control group of primers from ten anonymous microsatellite loci. Twenty-three of 24 examined conserved markers were polymorphic, each in on average 3 of 12 species tested.

Conclusion

Genomic sequence databases are useful tools to identify conserved genetic markers including those located in non-coding regions. By maximising primer sequence similarity between source species and database species, markers can be further improved and provide additional markers to study the molecular ecology of populations of non-model organisms.

A linkage map of the zebra finch Taeniopygia guttata provides new insights into avian genome evolution

Passeriformes are the largest order of birds and one of the most widely studied groups in evolutionary biology and ecology. Until recently genomic tools in passerines relied on chicken genomic resources. Here we report the construction and analysis of a whole-genome linkage map for the zebra finch (Taeniopygia guttata) using a 354-bird pedigree. The map contains 876 SNPs dispersed across 45 linkage groups and we found only a few instances of interchromosomal rearrangement between the zebra finch and the chicken genomes. Interestingly, there was a greater than expected degree of intrachromosomal rearrangements compared to the chicken, suggesting that gene order is not conserved within avian chromosomes. At 1068 cM the map is approximately only one quarter the length of the chicken linkage map, providing further evidence that the chicken has an unusually high recombination rate. Male and female linkage-map lengths were similar, suggesting no heterochiasmy in the zebra finch. This whole-genome map is the first for any passerine and a valuable tool for the zebra finch genome sequence project and for studies of quantitative trait loci.

Gene mapping in the wild with SNPs: guidelines and future directions

One of the biggest challenges facing evolutionary biologists is to identify and understand loci that explain fitness variation in natural populations. This review describes how genetic (linkage) mapping with single nucleotide polymorphism (SNP) markers can lead to great progress in this area. Strategies for SNP discovery and SNP genotyping are described and an overview of how to model SNP genotype information in mapping studies is presented. Finally, the opportunity afforded by new generation sequencing and typing technologies to map fitness genes by genome-wide association studies is discussed.

Extreme heterochiasmy and nascent sex chromosomes in European tree frogs

We investigated sex-specific recombination rates in Hyla arborea, a species with nascent sex chromosomes and male heterogamety. Twenty microsatellites were clustered into six linkage groups, all showing suppressed or very low recombination in males. Seven markers were sex linked, none of them showing any sign of recombination in males (r=0.00 versus 0.43 on average in females). This opposes classical models of sex chromosome evolution, which envision an initially small differential segment that progressively expands as structural changes accumulate on the Y chromosome. For autosomes, maps were more than 14 times longer in females than in males, which seems the highest ratio documented so far in vertebrates. These results support the pleiotropic model of Haldane and Huxley, according to which recombination is reduced in the heterogametic sex by general modifiers that affect recombination on the whole genome.

High degree of transferability of 86 newly developed zebra finch EST-linked microsatellite markers in 8 bird species

High-resolution analysis for population genetic and functional studies requires the use of large numbers of polymorphic markers. The recent increase of available genetic tools is facilitated by the use of publicly available expressed sequence tag (EST) sequence databases that are a valuable resource for identifying gene-linked markers. In the present study, we applied bioinformatics analyses to identify microsatellite markers present in EST sequences from a zebra finch (Taeniopgia guttata) EST database and we explore the success of cross-species amplification of EST-linked microsatellite markers in 7 passerine and 1 nonpasserine species. Eighty-six zebra finch EST-linked microsatellite loci were screened for polymorphism revealing a high amplification success rate and adequate levels of polymorphism (33.3-51%) for relatively closely related species, whereas success decreased in the most distantly related species to zebra finch. EST-linked microsatellites appear to be more highly transferable between taxa than anonymous microsatellites as they revealed higher amplification and polymorphism success between different families indicating that they will be a useful source of gene-linked polymorphic markers in a broad range of avian species.

A gene-based genetic linkage map of the collared flycatcher (Ficedula albicollis) reveals extensive synteny and gene-order conservation during 100 million years of avian evolution

By taking advantage of a recently developed reference marker set for avian genome analysis we have constructed a gene-based genetic map of the collared flycatcher, an important "ecological model" for studies of life-history evolution, sexual selection, speciation, and quantitative genetics. A pedigree of 322 birds from a natural population was genotyped for 384 single nucleotide polymorphisms (SNPs) from 170 protein-coding genes and 71 microsatellites. Altogether, 147 gene markers and 64 microsatellites form 33 linkage groups with a total genetic distance of 1787 cM. Male recombination rates are, on average, 22% higher than female rates (total distance 1982 vs. 1627 cM). The ability to anchor the collared flycatcher map with the chicken genome via the gene-based SNPs revealed an extraordinary degree of both synteny and gene-order conservation during avian evolution. The great majority of chicken chromosomes correspond to a single linkage group in collared flycatchers, with only a few cases of inter- and intrachromosomal rearrangements. The rate of chromosomal diversification, fissions/fusions, and inversions combined is thus considerably lower in birds (0.05/MY) than in mammals (0.6-2.0/MY). A dearth of repeat elements, known to promote chromosomal breakage, in avian genomes may contribute to their stability. The degree of genome stability is likely to have important consequences for general evolutionary patterns and may explain, for example, the comparatively slow rate by which genetic incompatibility among lineages of birds evolves.

Conservation of synteny between guppy and Xiphophorus genomes

The guppy and fish in the genus Xiphophorus have both been important model systems for the study of natural and sexual selection for over 50 years. Whereas the guppy is unique in the degree to which the environmental variables shaping phenotypic variation are known, Xiphophorus has the advantage that genomic resources have been developed due to the utility of this taxon for the study of melanoma. If linkage maps for the guppy and Xiphophorus are similar, genomic resources developed in Xiphophorus will be useful in the guppy. The authors used an F2 mapping cross of divergent populations of the guppy to construct partial female and male genetic linkage maps incorporating microsatellite markers derived from Xiphophorus mapping efforts. Flanking regions for a sample of microsatellites occurring in maps for both taxa were sequenced in the guppy and compared to published sequences from Xiphophorus. This confirmed that these loci were homologous and estimated the divergence in neutral nuclear DNA to be 0.21 substitutions per site. The female map comprises 16 linked markers on six linkage groups, and the male map comprises 24 markers on nine linkage groups. Linkage relationships among loci homologous in the guppy and Xiphophorus primarily show conservation of genetic architecture between species, but several major changes were detected.

Robustness of linkage maps in natural populations: a simulation study

In a number of long-term individual-based studies of vertebrate populations, the genealogical relationships between individuals have been established with molecular markers. As a result, it is possible to construct genetic linkage maps of these study populations by examining the co-segregation of markers through the pedigree. There are now four free-living vertebrate study populations for whom linkage maps have been built. In this study, simulation was used to investigate whether these linkage maps are likely to be accurate. In all four populations, the probability of assigning markers to the correct chromosome is high and framework maps are generally inferred correctly. However, genotyping error can result in incorrect maps being built with very strong statistical support over the correct order. Future applications of linkage maps of natural populations are discussed.