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

Heterochiasmy and Sexual Dimorphism: The Case of the Barn Swallow (Hirundo rustica, Hirundinidae, Aves)

Heterochiasmy, a sex-based difference in recombination rate, has been detected in many species of animals and plants. Several hypotheses about evolutionary causes of heterochiasmy were proposed. However, there is a shortage of empirical data. In this paper, we compared recombination related traits in females and males of the barn swallow Hirundo rustica (Linnaeus, 1758), the species under strong sexual selection, with those in the pale martin Riparia diluta (Sharpe and Wyatt, 1893), a related and ecologically similar species with the same karyotype (2N = 78), but without obvious sexual dimorphism. Recombination traits were examined in pachytene chromosome spreads prepared from spermatocytes and oocytes. Synaptonemal complexes and mature recombination nodules were visualized with antibodies to SYCP3 and MLH1 proteins, correspondingly. Recombination rate was significantly higher (p = 0.0001) in barn swallow females (55.6 ± 6.3 recombination nodules per autosomal genome), caused by the higher number of nodules at the macrochromosomes, than in males (49.0 ± 4.5). They also showed more even distribution of recombination nodules along the macrochromosomes. At the same time, in the pale martin, sexual differences in recombination rate and distributions were rather small. We speculate that an elevated recombination rate in the female barn swallows might have evolved as a compensatory reaction to runaway sexual selection in males.

The Impact of Natural Selection on Short Insertion and Deletion Variation in the Great Tit Genome

Insertions and deletions (INDELs) remain understudied, despite being the most common form of genetic variation after single nucleotide polymorphisms. This stems partly from the challenge of correctly identifying the ancestral state of an INDEL and thus identifying it as an insertion or a deletion. Erroneously assigned ancestral states can skew the site frequency spectrum, leading to artificial signals of selection. Consequently, the selective pressures acting on INDELs are, at present, poorly resolved. To tackle this issue, we have recently published a maximum likelihood approach to estimate the mutation rate and the distribution of fitness effects for INDELs. Our approach estimates and controls for the rate of ancestral state misidentification, overcoming issues plaguing previous INDEL studies. Here, we apply the method to INDEL polymorphism data from ten high coverage (∼44×) European great tit (Parus major) genomes. We demonstrate that coding INDELs are under strong purifying selection with a small proportion making it into the population (∼4%). However, among fixed coding INDELs, 71% of insertions and 86% of deletions are fixed by positive selection. In noncoding regions, we estimate ∼80% of insertions and ∼52% of deletions are effectively neutral, the remainder show signatures of purifying selection. Additionally, we see evidence of linked selection reducing INDEL diversity below background levels, both in proximity to exons and in areas of low recombination.

Divergent Evolutionary Trajectories of Two Young, Homomorphic, and Closely Related Sex Chromosome Systems

There exists extraordinary variation among species in the degree and nature of sex chromosome divergence. However, much of our knowledge about sex chromosomes is based on comparisons between deeply diverged species with different ancestral sex chromosomes, making it difficult to establish how fast and why sex chromosomes acquire variable levels of divergence. To address this problem, we studied sex chromosome evolution in two species of African clawed frog (Xenopus), both of whom acquired novel systems for sex determination from a recent common ancestor, and both of whom have female (ZW/ZZ) heterogamy. Derived sex chromosomes of one species, X. laevis, have a small region of suppressed recombination that surrounds the sex determining locus, and have remained this way for millions of years. In the other species, X. borealis, a younger sex chromosome system exists on a different pair of chromosomes, but the region of suppressed recombination surrounding an unidentified sex determining gene is vast, spanning almost half of the sex chromosomes. Differences between these sex chromosome systems are also apparent in the extent of nucleotide divergence between the sex chromosomes carried by females. Our analyses also indicate that in autosomes of both of these species, recombination during oogenesis occurs more frequently and in different genomic locations than during spermatogenesis. These results demonstrate that new sex chromosomes can assume radically different evolutionary trajectories, with far-reaching genomic consequences. They also suggest that in some instances the origin of new triggers for sex determination may be coupled with rapid evolution sex chromosomes, including recombination suppression of large genomic regions.

Broad-scale recombination pattern in the primitive bird Rhea americana (Ratites, Palaeognathae)

Birds have genomic and chromosomal features that make them an attractive group to analyze the evolution of recombination rate and the distribution of crossing over. Yet, analyses are biased towards certain species, especially domestic poultry and passerines. Here we analyze for the first time the recombination rate and crossover distribution in the primitive ratite bird, Rhea americana (Rheiformes, Palaeognathae). Using a cytogenetic approach for in situ mapping of crossovers we found that the total genetic map is 3050 cM with a global recombination rate of 2.1 cM/Mb for female rheas. In the five largest macrobivalents there were 3 or more crossovers in most bivalents. Recombination rates for macrobivalents ranges between 1.8-2.1 cM/Mb and the physical length of their synaptonemal complexes is highly predictive of their genetic lengths. The crossover rate at the pseudoautosomal region is 2.1 cM/Mb, similar to those of autosomal pairs 5 and 6 and only slightly higher compared to other macroautosomes. It is suggested that the presence of multiple crossovers on the largest macrobivalents is a feature common to many avian groups, irrespective of their position throughout phylogeny. These data provide new insights to analyze the heterogeneous recombination landscape of birds.

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.

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.

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.

Pronounced fixation, strong population differentiation and complex population history in the Canary Islands blue tit subspecies complex

Evolutionary molecular studies of island radiations may lead to insights in the role of vicariance, founder events, population size and drift in the processes of population differentiation. We evaluate the degree of population genetic differentiation and fixation of the Canary Islands blue tit subspecies complex using microsatellite markers and aim to get insights in the population history using coalescence based methods. The Canary Island populations were strongly genetically differentiated and had reduced diversity with pronounced fixation including many private alleles. In population structure models, the relationship between the central island populations (La Gomera, Tenerife and Gran Canaria) and El Hierro was difficult to disentangle whereas the two European populations showed consistent clustering, the two eastern islands (Fuerteventura and Lanzarote) and Morocco weak clustering, and La Palma a consistent unique lineage. Coalescence based models suggested that the European mainland forms an outgroup to the Afrocanarian population, a split between the western island group (La Palma and El Hierro) and the central island group, and recent splits between the three central islands, and between the two eastern islands and Morocco, respectively. It is clear that strong genetic drift and low level of concurrent gene flow among populations have shaped complex allelic patterns of fixation and skewed frequencies over the archipelago. However, understanding the population history remains challenging; in particular, the pattern of extreme divergence with low genetic diversity and yet unique genetic material in the Canary Island system requires an explanation. A potential scenario is population contractions of a historically large and genetically variable Afrocanarian population, with vicariance and drift following in the wake. The suggestion from sequence-based analyses of a Pleistocene extinction of a substantial part of North Africa and a Pleistocene/Holocene eastward re-colonisation of western North Africa from the Canaries remains possible.

Microsatellite resources for Passeridae species: a predicted microsatellite map of the house sparrow Passer domesticus

We identified microsatellite sequences of potential utility in the house sparrow (Passer domesticus) and assigned their predicted genome locations. These sequences included newly isolated house sparrow loci, which we fully characterized. Many of the newly isolated loci were polymorphic in two other species of Passeridae: Berthelot's pipit Anthus berthelotii and zebra finch Taeniopygia guttata. In total, we identified 179 microsatellite markers that were either isolated directly from, or are of known utility in, the house sparrow. Sixty-seven of these markers were designed from unique sequences that we isolated from a house sparrow genomic library. These new markers were combined with 36 house sparrow markers isolated by other studies and 76 markers isolated from other passerine species but known to be polymorphic in the house sparrow. We utilized sequence homology to assign chromosomal locations for these loci in the assembled zebra finch genome. One hundred and thirty-four loci were assigned to 25 different autosomes and eight loci to the Z chromosome. Examination of the genotypes of known-sex house sparrows for 37 of the new loci revealed a W-linked locus and an additional Z-linked locus. Locus Pdoμ2, previously reported as autosomal, was found to be Z-linked. These loci enable the creation of powerful and cost-effective house sparrow multiplex primer sets for population and parentage studies. They can be used to create a house sparrow linkage map and will aid the identification of quantitative trait loci in passerine species.

Characterization of a divergent chromosome region in the willow warbler Phylloscopus trochilus using avian genomic resources

Genome scans have made it possible to find outlier markers thought to have been influenced by divergent selection in almost any wild population. However, the lack of genomic information in nonmodel species often makes it difficult to associate these markers with certain genes or chromosome regions. Furthermore, the extent of linkage disequilibrium (LD) in the genome will determine the density of markers required to identify the genes under selection. In this study, we investigated a chromosome region in the willow warbler Phylloscopus trochilus surrounding a single marker previously identified in a genome scan. We first located the marker in the assembled genome of another species, the zebra finch Taeniopygia guttata, and amplified surrounding sequences in Fennoscandian willow warblers. Within an investigated chromosome region of 7.3 Mb as mapped to the zebra finch genome, we observed elevated genetic differentiation between a southern and a northern population across a 2.5-Mb interval comprising numerous coding genes. Within the southern and northern populations, higher values of LD were mostly found between SNPs within the same locus, but extended across distantly situated loci when the analyses were restricted to sampling sites showing intermediate allele frequencies of southern and northern alleles. Our study shows that cross-species genome information is a useful resource to obtain candidate sequences adjacent to outlier markers in nonmodel species.

General lack of global dosage compensation in ZZ/ZW systems? Broadening the perspective with RNA-seq

BACKGROUND

Species with heteromorphic sex chromosomes face the challenge of large-scale imbalance in gene dose. Microarray-based studies in several independent male heterogametic XX/XY systems suggest that dosage compensation mechanisms are in place to mitigate the detrimental effects of gene dose differences. However, recent genomic research on female heterogametic ZZ/ZW systems has generated surprising results. In two bird species and one lepidopteran no evidence for a global dosage compensating mechanism has been found. The recent advent of massively parallel RNA sequencing now opens up the possibility to gauge the generality of this observation with a broader phylogenetic sampling. It further allows assessing the validity of microarray-based inference on dosage compensation with a novel technology.

RESULTS

We here exemplify this approach using massively parallel sequencing on barcoded individuals of a bird species, the European crow (Corvus corone), where previously no genetic resources were available. Testing for Z-linkage with quantitative PCR (qPCR,) we first establish that orthology with distantly related species (chicken, zebra finch) can be used as a good predictor for chromosomal affiliation of a gene. We then use a digital measure of gene expression (RNA-seq) on brain transcriptome and confirm a global lack of dosage compensation on the Z chromosome. RNA-seq estimates of male-to-female (m:f) expression difference on the Z compare well to previous microarray-based estimates in birds and lepidopterans. The data further lends support that an up-regulation of female Z-linked genes conveys partial compensation and suggest a relationship between sex-bias and absolute expression level of a gene. Correlation of sex-biased gene expression on the Z chromosome across all three bird species further suggests that the degree of compensation has been partly conserved across 100 million years of avian evolution.

CONCLUSIONS

This work demonstrates that the study of dosage compensation has become amenable to species where previously no genetic resources were available. Massively parallele transcriptome sequencing allows re-assessing the degree of dosage compensation with a novel tool in well-studies species and, in addition, gain valuable insights into the generality of mechanisms across independent taxonomic group for both the XX/XY and ZZ/ZW system.

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.