Genetic mapping in a natural population of collared flycatchers (Ficedula albicollis): conserved synteny but gene order rearrangements on the avian Z chromosome

Dynamic Patterns of Sex Chromosome Evolution in Neognath Birds: Many Independent Barriers to Recombination at the ATP5F1A Locus

Avian sex chromosomes evolved after the divergence of birds and crocodilians from their common ancestor, so they are younger than the better-studied chromosomes of mammals. It has long been recognized that there may have been several stages to the evolution of avian sex chromosomes. For example, the CHD1 undergoes recombination in paleognaths but not neognaths. Genome assemblies have suggested that there may be variation in the timing of barriers to recombination among Neognathae, but there remains little understanding of the extent of this variability. Here, we look at partial sequences of ATP5F1A, which is on the avian Z and W chromosomes. It is known that recombination of this gene has independently ceased in Galliformes, Anseriformes, and at least five neoavian orders, but whether there are other independent cessations of recombination among Neoaves is not understood. We analyzed a combination of data extracted from published chromosomal-level genomes with data collected using PCR and cloning to identify Z and W copies in 22 orders. Our results suggest that there may be at least 19 independent cessations of recombination within Neognathae, and 3 clades that may still be undergoing recombination (or have only recently ceased recombination). Analyses of ATP5F1A protein sequences revealed an increased amino acid substitution rate for W chromosome gametologs, suggesting relaxed purifying selection on the W chromosome. Supporting this hypothesis, we found that the increased substitution rate was particularly pronounced for buried residues, which are expected to be more strongly constrained by purifying selection. This highlights the dynamic nature of avian sex chromosomes, and that this level of variation among clades means they should be a good system to understand sex chromosome evolution.

Quaternary land bridges have not been universal conduits of gene flow

Quaternary climate oscillations are a well-known driver of animal diversification, but their effects are most well studied in areas where glaciations lead to habitat fragmentation. In large areas of the planet, however, glaciations have had the opposite effect, but here their impacts are much less well understood. This is especially true in Southeast Asia, where cyclical changes in land distribution have generated enormous land expansions during glacial periods. In this study, we selected a panel of five songbird species complexes covering a range of ecological specificities to investigate the effects Quaternary land bridges have had on the connectivity of Southeast Asian forest biota. Specifically, we combined morphological and bioacoustic analysis with an arsenal of population genomic and modelling approaches applied to thousands of genome-wide DNA markers across a total of more than 100 individuals. Our analyses show that species dependent on forest understorey exhibit deep differentiation between Borneo and western Sundaland, with no evidence of gene flow during the land bridges accompanying the last 1-2 ice ages. In contrast, dispersive canopy species and habitat generalists have experienced more recent gene flow. Our results argue that there remains much cryptic species-level diversity to be discovered in Southeast Asia even in well-known animal groups such as birds, especially in nondispersive forest understorey inhabitants. We also demonstrate that Quaternary land bridges have not been equally suitable conduits of gene flow for all species complexes and that life history is a major factor in predicting relative population divergence time across Quaternary climate fluctuations.

Emerging patterns of genome organization in Notopteridae species (Teleostei, Osteoglossiformes) as revealed by Zoo-FISH and Comparative Genomic Hybridization (CGH)

Notopteridae (Teleostei, Osteoglossiformes) represents an old fish lineage with ten currently recognized species distributed in African and Southeastern Asian rivers. Their karyotype structures and diploid numbers remained conserved over long evolutionary periods, since African and Asian lineages diverged approximately 120 Mya. However, a significant genetic diversity was already identified for these species using molecular data. Thus, why the evolutionary relationships within Notopteridae are so diverse at the genomic level but so conserved in terms of their karyotypes? In an attempt to develop a more comprehensive picture of the karyotype and genome evolution in Notopteridae, we performed comparative genomic hybridization (CGH) and cross-species (Zoo-FISH) whole chromosome painting experiments to explore chromosome-scale intergenomic divergence among seven notopterid species, collected in different African and Southeast Asian river basins. CGH demonstrated an advanced stage of sequence divergence among the species and Zoo-FISH experiments showed diffuse and limited homology on inter-generic level, showing a temporal reduction of evolutionarily conserved syntenic regions. The sharing of a conserved chromosomal region revealed by Zoo-FISH in these species provides perspectives that several other homologous syntenic regions have remained conserved among their genomes despite long temporal isolation. In summary, Notopteridae is an interesting model for tracking the chromosome evolution as it is (i) ancestral vertebrate group with Gondwanan distribution and (ii) an example of animal group exhibiting karyotype stasis. The present study brings new insights into degree of genome divergence vs. conservation at chromosomal and sub-chromosomal level in representative sampling of this group.

The niches of nuthatches affect their lineage evolution differently across latitude

Ecological niche evolution can promote or hinder the differentiation of taxa and determine their distribution. Niche-mediated evolution may differ among climatic regimes, and thus, species that occur across a wide latitudinal range offer a chance to test these heterogeneous evolutionary processes. In this study, we examine (a) how many lineages have evolved across the continent-wide range of the Eurasian nuthatch (Sitta europaea), (b) whether the lineages' niches are significantly divergent or conserved and (c) how their niche evolution explains their geographic distribution. Phylogenetic reconstruction and ecological niche models (ENMs) showed that the Eurasian nuthatch contained six parapatric lineages that diverged within 2 Myr and did not share identical climatic niches. However, the niche discrepancy between these distinct lineages was relatively conserved compared with the environmental differences between their ranges and thus was unlikely to drive lineage divergence. The ENMs of southern lineages tended to cross-predict with their neighbouring lineages whereas those of northern lineages generally matched with their abutting ranges. The coalescence-based analyses revealed more stable populations for the southern lineages than the northern ones during the last glaciation cycle. In contrast to the overlapping ENMs, the smaller parapatric distribution suggests that the southern lineages might have experienced competitive exclusion to prevent them from becoming sympatric. On the other hand, the northern lineages have expanded their ranges and their current abutting distribution might have resulted from lineages adapting to different climatic conditions in allopatry. This study suggests that niche evolution may affect lineage distribution in different ways across latitude.

Breeding consequences of flavivirus infection in the collared flycatcher

Background

The breeding consequences of virus infections have rarely been studied in avian natural breeding populations. In this paper we investigated the links between humoral immunity following a natural flavivirus infection and reproduction in a wild bird population of collared flycatcher (Ficedula albicollis). We analyzed plasma from 744 birds for antibodies and correlated these results to a number of reproductive components.

Results

Nearly one third (27.8%) of the sampled collared flycatchers were found seropositive for flavivirus. Males had significantly more frequently flavivirus antibodies (32.3%) than females (25.1%). Seropositive females differed significantly from seronegative females in four traits: they had earlier lay date, higher body weight, higher survival rate and were older than seronegative females. The females did not differ in clutch size, number of fledged young or number of recruited young. Seropositive males had female partners with earlier lay date, i.e. the males bred earlier and they also produced more fledged young than seronegative males. In contrast, the males did not differ in clutch size, number of recruited young, male weight, age or survival. Interestingly, seropositive males had larger ornament, forehead badge size, than seronegative males.

Conclusions

Collared flycatchers with an antibody response against flavivirus were more successful than birds with no antibody response, for any of the measured life history traits. The positive link between flavivirus antibody presence and life-history trait levels suggest that it is condition dependent in the collared flycatcher.

Erbin and ErbB2 play roles in the sexual differentiation of the song system nucleus HVC in bengalese finches (Lonchura Striata var. domestica)

Song control nuclei have distinct sexual differences in songbirds. However, the mechanism that underlies the sexual differentiation of song nuclei is still not well understood. Using a combination of anatomical, pharmacological, genetic, and behavioral approaches, the present study investigated the role of erbb2 (a homolog of the avian erythroblastic leukemia viral oncogene homolog 2) and the erbb2-interacting gene, erbin, in the sexual differentiation of the song nucleus HVC in the Bengalese finch. We first found that both erbin and erbb2 were expressed in the developing HVC at posthatch day (PHD) 15 in a male-biased fashion using qRT-PCR and in situ hybridization. Following the addition of a pharmaceutical inhibitor of the ErbB2 signaling pathway to the culture medium, cell proliferation in the cultured ventricle zone (VZ) that overlies the developing HVC decreased significantly. After the injection of erbin- or erbb2-interfering lentiviruses into the HVC and its overlying VZ at PHD 15, the cell proliferation in the VZ at PHD 24, the number of the differentiated neurons (Hu⁺ /BrdU⁺ or NeuN⁺ /BrdU⁺ ) in the HVC at PHD 31 or PHD 130, and the number of RA-projecting cells at PHD 130 all decreased significantly. Additionally, the adult songs displayed serious abnormalities. Finally, 173 male-biased genes were expressed in the developing HVC at PHD 15 using cDNA microarrays, of which 27.2% were Z-linked genes and approximately 20 genes were involved in the Erbin- or ErbB2-related signaling pathways. Our results provide some specific genetic factors that contribute to neurogenesis and sex differentiation in a song nucleus of songbirds. © 2017 Wiley Periodicals, Inc. Develop Neurobiol 78: 15-38, 2018.

Matching loci surveyed to questions asked in phylogeography

Although mitochondrial DNA (mtDNA) has long been used for assessing genetic variation within and between populations, its workhorse role in phylogeography has been criticized owing to its single-locus nature. The only choice for testing mtDNA results is to survey nuclear loci, which brings into contrast the difference in locus effective size and coalescence times. Thus, it remains unclear how erroneous mtDNA-based estimates of species history might be, especially for evolutionary events in the recent past. To test the robustness of mtDNA and nuclear sequences in phylogeography, we provide one of the largest paired comparisons of summary statistics and demographic parameters estimated from mitochondrial, five Z-linked and 10 autosomal genes of 30 avian species co-distributed in the Caucasus and Europe. The results suggest that mtDNA is robust in estimating inter-population divergence but not in intra-population diversity, which is sensitive to population size change. Here, we provide empirical evidence showing that mtDNA was more likely to detect population divergence than any other single locus owing to its smaller Ne and thus faster coalescent time. Therefore, at least in birds, numerous studies that have based their inferences of phylogeographic patterns solely on mtDNA should not be readily dismissed.

DLGP: A database for lineage-conserved and lineage-specific gene pairs in animal and plant genomes

The conservation of gene organization in the genome with lineage-specificity is an invaluable resource to decipher their potential functionality with diverse selective constraints, especially in higher animals and plants. Gene pairs appear to be the minimal structure for such kind of gene clusters that tend to reside in their preferred locations, representing the distinctive genomic characteristics in single species or a given lineage. Despite gene families having been investigated in a widespread manner, the definition of gene pair families in various taxa still lacks adequate attention. To address this issue, we report DLGP (http://lcgbase.big.ac.cn/DLGP/) that stores the pre-calculated lineage-based gene pairs in currently available 134 animal and plant genomes and inspect them under the same analytical framework, bringing out a set of innovational features. First, the taxonomy or lineage has been classified into four levels such as Kingdom, Phylum, Class and Order. It adopts all-to-all comparison strategy to identify the possible conserved gene pairs in all species for each gene pair in certain species and reckon those that are conserved in over a significant proportion of species in a given lineage (e.g. Primates, Diptera or Poales) as the lineage-conserved gene pairs. Furthermore, it predicts the lineage-specific gene pairs by retaining the above-mentioned lineage-conserved gene pairs that are not conserved in any other lineages. Second, it carries out pairwise comparison for the gene pairs between two compared species and creates the table including all the conserved gene pairs and the image elucidating the conservation degree of gene pairs in chromosomal level. Third, it supplies gene order browser to extend gene pairs to gene clusters, allowing users to view the evolution dynamics in the gene context in an intuitive manner. This database will be able to facilitate the particular comparison between animals and plants, between vertebrates and arthropods, and between monocots and eudicots, accounting for the significant contribution of gene pairs to speciation and diversification in specific lineages.

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.

Differential introgression and effective size of marker type influence phylogenetic inference of a recently divergent avian group (Phasianidae: Tympanuchus)

Life history strategies can influence the effective population size (Ne) of loci differently based on their mode of inheritance. Recognizing how this may affect the rate of lineage sorting among marker types is important for studies focused on resolving phylogenetic relationships among recently divergent taxa. In this study, we use gene tree, coalescent-based species tree, and isolation-with-migration analyses to explore the differences between marker types (autosomal, Z-linked, and mitochondrial) in resolving phylogenetic relationships among North American prairie grouse (Tympanuchus). We found that Z-linked loci were more likely to identify monophyletic relationships among prairie grouse species compared to autosomal and mtDNA loci in both species and gene tree analyses, with species tree analyses outperforming gene trees. These results were further supported with isolation-with-migration analyses, where Z-linked loci largely followed a strict isolation model while autosomal loci were more likely to fit a model with gene flow between species following population divergence. While accounting for differences in inheritance pattern (or Ne) for marker type, results suggest that additional factors, such as strong sexual selection and sex-biased introgression (i.e., male-biased postzygotic hybrid behavioral isolation or "unsexy son"), may further explain the decreased diversity levels and increased rate of lineage sorting observed with the Z-linked loci relative to autosomal and mtDNA loci. In fact, to our knowledge no hybrid male prairie grouse have been observed breeding in the wild, yet hybrid females along with backcross females are known to produce viable offspring. Overall, this study highlights that more work is needed to determine how complex models of gene flow (i.e., sex biased introgression) and differences in the effective size among marker types based on differing life history strategies influence divergence date estimation and species delimitation.

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.

Genomic structures of the kW1 loci on the Z and W chromosomes in ratite birds: structural changes at an early stage of W chromosome differentiation

The W chromosome of ratite birds shows minimal morphological differentiation and retains homology of genetic linkage and gene order with a substantial stretch of the Z chromosome; however, the molecular structure in the differentiated region is still not well known. The kW1 sequence was isolated from the kiwi as a W-specific DNA marker for PCR-based molecular sexing of ratite birds. In ratite W chromosomes, this sequence commonly contains a ∼200-bp deletion. To characterize the very early event of avian sex chromosome differentiation, we performed molecular cytogenetic analyses of kW1 and its flanking sequences in paleognathous and neognathous birds and reptiles. Female-specific repeats were found in the kW1-flanking sequence of the cassowary (Casuarius casuarius), and the repeats have been amplified in the pericentromeric region of the W chromosomes of ratites, which may have resulted from the cessation of meiotic recombination between the Z and W chromosomes at an early stage of sex chromosome differentiation. The presence of the kW1 sequence in neognathous birds and a crocodilian species suggests that the kW1 sequence was present in the ancestral genome of Archosauria; however, it disappeared in other reptilian taxa and several lineages of neognathous birds.

Deep phylogeographic divergence of a migratory passerine in Sino-Himalayan and Siberian forests: the Red-flanked Bluetail (Tarsiger cyanurus) complex

Enormous mountainous forests in Sino-Himalayans and Siberia harbor important avian biodiversity in the Northern Hemisphere. Numerous studies in last two decades have been contributed to systematics and taxonomy of passerines birds in these regions and have revealed various and complex phylogeographic patterns. A passerine species Red-flanked Bluetail Tarsiger cyanurus provided a good system to manifest such evolutionary complexity. The subspecies T. c. cyanurus and T. c. rufilatus (or/and T. c. pallidior), divergent in morphology, acoustics, and migratory strategies are allopatric in Siberia and Sino-Himalayan forests, respectively. The two taxa most likely deserve full species status but rigorous genetic analysis is missing. In this study, multilocus phylogeography based on mitochondrial DNA and Z-linked DNA reveals that T. c. cyanurus and T. c. rufilatus are reciprocally monophyletic with significant statistical support and differ with a large number of diagnostic nucleotide sites resulting substantial genetic divergence. Our finding supports the proposed split of Tarsiger cyanurus s.l. that T. cyanurus and T. rufilatus should be treated as two full species. Whether "pallidior" is a subspecies or geographical form of T. rufilatus is still uncertain. Additionally, these two forest passerine species may have diverged 1.88 (3.25-1.30) Mya, which might be shaped by geographical vicariance due to grassland and desert steppe on the central Loess Plateau during the Pliocene. Taken together, this study and further suggests another independent example of North Palearctic-Sino-Himalayan phylogeographic pattern in Palearctic birds.

Distinguishing the effects of selection from demographic history in the genetic variation of two sister passerines based on mitochondrial-nuclear comparison

Determining the mechanisms responsible for the distribution of genetic diversity in natural populations has occupied a central role in molecular evolution. Our study was motivated by the unprecedented observation that a widespread Eurasian flycatcher, Ficedula albicilla, exhibited no variation at the mitochondrial DNA (mtDNA) ND2 gene in 75 individuals sampled over a 5000-km distance. In contrast, its sister species, F. parva, had low but considerably higher levels of mtDNA variation. We assessed whether natural selection or demographic factors could explain the absence of mtDNA variation in F. albicilla. Eighteen nuclear genes were sequenced to estimate the two species' phylogeographic histories, and for comparison to the mtDNA data. Multilocus coalescence analyses suggested that F. albicilla experienced a population expansion perhaps following a population bottleneck. Simulations based on this demographic history, however, did not replicate the extremely low level of mtDNA variation. Historical range changes based on ecological niche models also failed to explain the observed mtDNA patterns. Neutrality tests (DHEW and ML-HKA) suggested a non-neutral pattern in the mtDNA of F. albicilla. We found a transmembrane-skewed distribution of nonsynonymous substitutions between the two species, three of which caused functional change; the results implied that positive selection could have targeted mtDNA. Several lines of evidence support selection rather than demographic history as the main force influencing the patterns of mtDNA variation. Despite the influence of natural selection, many of the phylogeographic inferences derived from mtDNA were robust, including species limits and a high level of gene flow among populations within species.

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.

Demographic divergence history of pied flycatcher and collared flycatcher inferred from whole-genome re-sequencing data

Profound knowledge of demographic history is a prerequisite for the understanding and inference of processes involved in the evolution of population differentiation and speciation. Together with new coalescent-based methods, the recent availability of genome-wide data enables investigation of differentiation and divergence processes at unprecedented depth. We combined two powerful approaches, full Approximate Bayesian Computation analysis (ABC) and pairwise sequentially Markovian coalescent modeling (PSMC), to reconstruct the demographic history of the split between two avian speciation model species, the pied flycatcher and collared flycatcher. Using whole-genome re-sequencing data from 20 individuals, we investigated 15 demographic models including different levels and patterns of gene flow, and changes in effective population size over time. ABC provided high support for recent (mode 0.3 my, range <0.7 my) species divergence, declines in effective population size of both species since their initial divergence, and unidirectional recent gene flow from pied flycatcher into collared flycatcher. The estimated divergence time and population size changes, supported by PSMC results, suggest that the ancestral species persisted through one of the glacial periods of middle Pleistocene and then split into two large populations that first increased in size before going through severe bottlenecks and expanding into their current ranges. Secondary contact appears to have been established after the last glacial maximum. The severity of the bottlenecks at the last glacial maximum is indicated by the discrepancy between current effective population sizes (20,000-80,000) and census sizes (5-50 million birds) of the two species. The recent divergence time challenges the supposition that avian speciation is a relatively slow process with extended times for intrinsic postzygotic reproductive barriers to evolve. Our study emphasizes the importance of using genome-wide data to unravel tangled demographic histories. Moreover, it constitutes one of the first examples of the inference of divergence history from genome-wide data in non-model species.

TECHNICAL ADVANCES: A microarray for large-scale genomic and transcriptional analyses of the zebra finch (Taeniopygia guttata) and other passerines

The microarray technology has revolutionized biological research in the last decade. By monitoring the expression of many genes simultaneously, microarrays can elucidate gene function, as well as scan entire genomes for candidate genes encoding complex traits. However, because of high costs of sequencing and design, microarrays have largely been restricted to a few model species. Cross-species microarray (CSM) analyses, where microarrays are used for other species than the one they were designed for, have had varied success. We have conducted a CSM analysis by hybridizing genomic DNA from the common whitethroat (Sylvia communis) on a newly developed Affymetrix array designed for the zebra finch (Taeniopygia guttata), the Lund-zf array. The results indicate a very high potential for the zebra finch array to act as a CSM utility in other passerine birds. When hybridizing zebra finch genomic DNA, 98% of the gene representatives had higher signal intensities than the background cut-off, and for the common whitethroat, we found the equivalent proportion to be as high as 96%. This was surprising given the fact that finches and warblers diverged 25-50 million years ago, but may be explained by a relatively low sequence divergence between passerines (89-93%). Passerine birds are widely used in studies of ecology and evolution, and a zebra finch array that can be used for many species may have a large impact on future research directions.

Evidence for asymmetrical divergence-gene flow of nuclear loci, but not mitochondrial loci, between seabird sister species: blue-footed (Sula nebouxii) and Peruvian (S. variegata) boobies

Understanding the process of speciation requires understanding how gene flow influences divergence. Recent analyses indicate that divergence can take place despite gene flow and that the sex chromosomes can exhibit different levels of gene flow than autosomes and mitochondrial DNA. Using an eight marker dataset including autosomal, z-linked, and mitochondrial loci we tested the hypothesis that blue-footed (Sula nebouxii) and Peruvian (S. variegata) boobies diverged from their common ancestor with gene flow, paying specific attention to the differences in gene flow estimates from nuclear and mitochondrial markers. We found no gene flow at mitochondrial markers, but found evidence from the combined autosomal and z-linked dataset that blue-footed and Peruvian boobies experienced asymmetrical gene flow during or after their initial divergence, predominantly from Peruvian boobies into blue-footed boobies. This gene exchange may have occurred either sporadically between periods of allopatry, or regularly throughout the divergence process. Our results add to growing evidence that diverging species can remain distinct but exchange genes.

The causes and evolutionary consequences of mixed singing in two hybridizing songbird species (Luscinia spp.)

Bird song plays an important role in the establishment and maintenance of prezygotic reproductive barriers. When two closely related species come into secondary contact, song convergence caused by acquisition of heterospecific songs into the birds’ repertoires is often observed. The proximate mechanisms responsible for such mixed singing, and its effect on the speciation process, are poorly understood. We used a combination of genetic and bioacoustic analyses to test whether mixed singing observed in the secondary contact zone of two passerine birds, the Thrush Nightingale (Luscinia luscinia) and the Common Nightingale (L. megarhynchos), is caused by introgressive hybridization. We analysed song recordings of both species from allopatric and sympatric populations together with genotype data from one mitochondrial and seven nuclear loci. Semi-automated comparisons of our recordings with an extensive catalogue of Common Nightingale song types confirmed that most of the analysed sympatric Thrush Nightingale males were ‘mixed singers’ that use heterospecific song types in their repertoires. None of these ‘mixed singers’ possessed any alleles introgressed from the Common Nightingale, suggesting that they were not backcross hybrids. We also analysed songs of five individuals with intermediate phenotype, which were identified as F₁ hybrids between the Thrush Nightingale female and the Common Nightingale male by genetic analysis. Songs of three of these hybrids corresponded to the paternal species (Common Nightingale) but the remaining two sung a mixed song. Our results suggest that although hybridization might increase the tendency for learning songs from both parental species, interspecific cultural transmission is the major proximate mechanism explaining the occurrence of mixed singers among the sympatric Thrush Nightingales. We also provide evidence that mixed singing does not substantially increase the rate of interspecific hybridization and discuss the possible adaptive value of this phenomenon in nightingales.

Inferring the demographic history of European Ficedula flycatcher populations

Background

Inference of population and species histories and population stratification using genetic data is important for discriminating between different speciation scenarios and for correct interpretation of genome scans for signs of adaptive evolution and trait association. Here we use data from 24 intronic loci re-sequenced in population samples of two closely related species, the pied flycatcher and the collared flycatcher.

Results

We applied Isolation-Migration models, assignment analyses and estimated the genetic differentiation and diversity between species and between populations within species. The data indicate a divergence time between the species of <1 million years, significantly shorter than previous estimates using mtDNA, point to a scenario with unidirectional gene-flow from the pied flycatcher into the collared flycatcher and imply that barriers to hybridisation are still permeable in a recently established hybrid zone. Furthermore, we detect significant population stratification, predominantly between the Spanish population and other pied flycatcher populations.

Conclusions

Our results provide further evidence for a divergence process where different genomic regions may be at different stages of speciation. We also conclude that forthcoming analyses of genotype-phenotype relations in these ecological model species should be designed to take population stratification into account.

Recent allopatric divergence and niche evolution in a widespread Palearctic bird, the common rosefinch (Carpodacus erythrinus)

A previously published phylogeographic analysis of mtDNA sequences from the widespread Palearctic common rosefinch (Carpodacus erythrinus) suggested the existence of three recently diverged groups, corresponding to the Caucasus, central-western Eurasia, and northeastern Eurasia. We re-evaluated the mtDNA data using coalescence methods and added sequence data from a sex-linked gene. The mtDNA gene tree and SAMOVA supported the distinctiveness of the Caucasian group but not the other two groups. However, UPGMA clustering of mtDNA Φ(ST)-values among populations recovered the three groups. The sex-linked gene tree recovered no phylogeographic signal, which was attributed to recent divergence and insufficient time for sorting of alleles. Overall, coalescence methods indicated a lack of gene flow among the three groups, and population expansion in the central-western and northeastern Eurasia groups. These three groups corresponded to named subspecies, further supporting their validity. A species distribution model revealed potential refugia at the Last Glacial Maximum. These three groups, which we hypothesized are in the early stages of speciation, provided an opportunity for testing tenets of ecological speciation. We showed that the early stages of speciation were not accompanied by ecological niche divergence, consistent with other avian studies.

Extensive introgressive hybridization within the northern oriole group (Genus Icterus) revealed by three-species isolation with migration analysis

Until recently, studies of divergence and gene flow among closely-related taxa were generally limited to pairs of sister taxa. However, organisms frequently exchange genes with other non-sister taxa. The “northern oriole” group within genus Icterus exemplifies this problem. This group involves the extensively studied hybrid zone between Baltimore oriole (Icterus galbula) and Bullock's oriole (I. bullockii), an alleged hybrid zone between I. bullockii and black-backed oriole (I. abeillei), and likely mtDNA introgression between I. galbula and I. abeillei. Here, we examine the divergence population genetics of the entire northern oriole group using a multipopulation Isolation-with-Migration (IM) model. In accordance with Haldane's rule, nuclear loci introgress extensively beyond the I. galbula–I. bullockii hybrid zone, while mtDNA does not. We found no evidence of introgression between I. bullockii and I. abeillei or between I. galbula and I. abeillei when all three species were analyzed together in a three-population model. However, traditional pairwise analysis suggested some nuclear introgression from I. abeillei into I. galbula, probably reflecting genetic contributions from I. bullockii unaccounted for in a two-population model. Thus, only by including all members of this group in the analysis was it possible to rigorously estimate the level of gene flow among these three closely related species.

A genetic polymorphism in the sex-linked ATP5A1 gene is associated with individual fitness in Ovenbirds (Seiurus aurocapilla)

While testing genetic sexing techniques in Ovenbirds (Seiurus aurocapilla), we found a genetic polymorphism in the ATP5A1 gene in 38% of individuals. The Z′ allele included changes in both intronic and exonic portions of the sequenced region, but there was no evidence that this changed the resulting ATP synthase product. Males that had one or more copies of this allele had higher relative body mass (mass corrected for size) than other genotypes. This allele was unrelated to stable isotope signatures, and so was not a useful predictor of latitude within the eastern portion of the Ovenbird breeding range. Future studies are needed to determine whether this polymorphism may be a useful geographic marker. This study is the first to link polymorphisms in the sex-linked ATP5A1 gene with fitness effects.

Multilocus coalescence analyses support a mtDNA-based phylogeographic history for a widespread Palearctic passerine bird, Sitta europaea

Our understanding of species phylogeography in much of the Palearctic is incomplete. In addition, many existing studies based solely on mitochondrial DNA (mtDNA) can provide a biased view of phylogeographic history because of the effects of lineage sorting, natural selection, or hybridization. We analyzed 13 introns to assess a mtDNA study of the Eurasian nuthatch (Sitta europaea) that suggested a seemingly contemporaneous origin of distinct taxa in the Caucasus, Europe, and Asia. Neutrality tests showed no evidence of selection on either the mtDNA or nuclear sequences. Most nuclear gene trees, except for Z-linked ones, did not recover the three lineages, which we attribute to recent splitting. Analyses of the 13 introns combined revealed the same three groups as did the mtDNA and suggested that nuthatches experienced a trichotomous (or two indistinguishable) split(s) 1-2 million years ago (Mya) and have remained isolated with trifling if not zero gene flow since then, and the Asian group increased in population size. This result demonstrates the usefulness of mtDNA in discovering phylogeographic patterns. The use of multiple nuclear loci facilitated detection of an introgressed individual and improved estimates of process parameters such as divergence time and population expansion. We recommend that phylogeographic studies should be based on both mtDNA and nuclear genes.

Hybrid speciation in sparrows II: a role for sex chromosomes?

Homoploid hybrid speciation in animals is poorly understood, mainly because of the scarcity of well-documented cases. Here, we present the results of a multilocus sequence analysis on the house sparrow (Passer domesticus), Spanish sparrow (P. hispaniolensis) and their proposed hybrid descendant, the Italian sparrow (P. italiae). The Italian sparrow is shown to be genetically intermediate between the house sparrow and Spanish sparrow, exhibiting genealogical discordance and a mosaic pattern of alleles derived from either of the putative parental species. The average variation on the Z chromosome was significantly reduced compared with autosomal variation in the putative parental species, the house sparrow and Spanish sparrow. Additionally, divergence between the two species was elevated on the Z chromosome relative to the autosomes. This pattern of variation and divergence is consistent with reduced introgression of Z-linked genes and/or a faster-Z effect (increased rate of adaptive divergence on the Z). F(ST) -outlier tests were consistent with the faster-Z hypothesis: two of five Z-linked loci (CHD1Z and PLAA) were identified as candidates for being subject to positive, divergent selection in the putative parental species. Interestingly, the two latter genes showed a mosaic pattern in the (hybrid) Italian sparrow; that is, the Italian sparrow was found to be fixed for Spanish sparrow alleles at CHD1Z and to mainly have house sparrow alleles at PLAA. Preliminary evidence presented in this study thus suggests that sex chromosomes may play a significant role in this case of homoploid hybrid speciation.

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.

Beyond a morphological paradox: complicated phylogenetic relationships of the parrotbills (Paradoxornithidae, Aves)

The parrotbills (Paradoxornithidae, meaning "birds of paradox," Aves) are a group of Old World passerines with perplexing taxonomic histories due to substantial morphological and ecological variation at various levels. In this study, phylogenetic relationships of the parrotbills were reconstructed based on sequences of two mitochondrial segments and three nuclear coding regions. Three major clades with characteristic body size and plumage coloration were found in both mtDNA and nuclear gene trees. However, mtDNA phylogeny suggested that the Paradoxornithidae is paraphyletic and relationships among three major parrotbill clades were poorly resolved. On the contrary, apparent and well-supported monophyletic relationships among the three major clades of Paradoxornithidae were revealed by concatenated nuclear dataset. Since paraphyly based on mtDNA data has commonly been found within avian taxa, the conflicting phylogenetic signal between mtDNA and nuclear loci revealed in this study indicates that results obtained from mtDNA dataset alone need to be evaluated with caution. Taxonomic implications of our phylogenetic findings are discussed. These phylogenies also point out areas for future investigation regarding the rapid diversification, morphological evolution and environmental adaptation of various parrotbill species or species complexes.

Samia cynthia versus Bombyx mori: comparative gene mapping between a species with a low-number karyotype and the model species of Lepidoptera

We performed gene-based comparative FISH mapping between a wild silkmoth, Samia cynthia ssp. with a low number of chromosomes (2n=25-28) and the model species, Bombyx mori (2n=56), in order to identify the genomic components that make up the chromosomes in a low-number karyotype. Mapping of 64 fosmid probes containing orthologs of B. mori genes revealed that the homologues of either two or four B. mori chromosomes constitute the S. c. ricini (Vietnam population, 2n=27♀/28♂, Z0/ZZ) autosomes. Where tested, even the gene order was conserved between S. c. ricini and B. mori. This was also true for the originally autosomal parts of the neo-sex chromosomes in S. c. walkeri (Sapporo population, 2n=26♀/26♂, neo-Wneo-Z/neo-Zneo-Z) and S. cynthia subsp. indet. (Nagano population, 2n=25♀/26♂, neo-WZ₁Z₂/Z₁Z₁Z₂Z₂). The results are evidence for an internal stability of lepidopteran chromosomes even when all autosomes had undergone fusion processes to form a low-number karyotype.

Sex- and species-biased gene flow in a spotted eagle hybrid zone

Background

Recent theoretical and empirical work points toward a significant role for sex-chromosome linked genes in the evolution of traits that induce reproductive isolation and for traits that evolve under influence of sexual selection. Empirical studies including recently diverged (Pleistocene), short-lived avian species pairs with short generation times have found that introgression occurs on the autosomes but not on the Z-chromosome. Here we study genetic differentiation and gene flow in the long-lived greater spotted eagle (Aquila clanga) and lesser spotted eagle (A. pomarina), two species with comparatively long generation times.

Results

Our data suggest that there is a directional bias in migration rates between hybridizing spotted eagles in eastern Europe. We find that a model including post divergence gene flow fits our data best for both autosomal and Z-chromosome linked loci but, for the Z-chromosome, the rate is reduced in the direction from A. pomarina to A. clanga.

Conclusions

The fact that some introgression still occurs on the Z-chromosome between these species suggests that the differentiation process is in a more premature phase in our study system than in previously studied avian species pairs and that could be explained by a shorter divergence time and/or a longer average generation time in the spotted eagles. The results are in agreement with field observations and provide further insight into the role of sex-linked loci for the build-up of barriers to gene flow among diverging populations and 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.

Genome evolution in Reptilia, the sister group of mammals

The genomes of birds and nonavian reptiles (Reptilia) are critical for understanding genome evolution in mammals and amniotes generally. Despite decades of study at the chromosomal and single-gene levels, and the evidence for great diversity in genome size, karyotype, and sex chromosome diversity, reptile genomes are virtually unknown in the comparative genomics era. The recent sequencing of the chicken and zebra finch genomes, in conjunction with genome scans and the online publication of the Anolis lizard genome, has begun to clarify the events leading from an ancestral amniote genome--predicted to be large and to possess a diverse repeat landscape on par with mammals and a birdlike sex chromosome system--to the small and highly streamlined genomes of birds. Reptilia exhibit a wide range of evolutionary rates of different subgenomes and, from isochores to mitochondrial DNA, provide a critical contrast to the genomic paradigms established in mammals.

No evidence for Z-chromosome rearrangements between the pied flycatcher and the collared flycatcher as judged by gene-based comparative genetic maps

Revealing the genetic basis of reproductive isolation is fundamental for understanding the speciation process. Chromosome speciation models propose a role for chromosomal rearrangements in promoting the build up of reproductive isolation between diverging populations and empirical data from several animal and plant taxa support these models. The pied flycatcher and the collared flycatcher are two closely related species that probably evolved reproductive isolation during geographical separation in Pleistocene glaciation refugia. Despite the short divergence time and current hybridization, these two species demonstrate a high degree of intrinsic post-zygotic isolation and previous studies have shown that traits involved in mate choice and hybrid viability map to the Z-chromosome. Could rearrangements of the Z-chromosome between the species explain their reproductive isolation? We developed high coverage Z-chromosome linkage maps for both species, using gene-based markers and large-scale SNP genotyping. Best order maps contained 57-62 gene markers with an estimated average density of one every 1-1.5 Mb. We estimated the recombination rates in flycatcher Z-chromosomes to 1.1-1.3 cM/Mb. A comparison of the maps of the two species revealed extensive co-linearity with no strong evidence for chromosomal rearrangements. This study does therefore not provide support the idea that sex chromosome rearrangements have caused the relatively strong post-zygotic reproductive isolation between these two Ficedula species.

Speciation in Ficedula flycatchers

Speciation in animals often requires that population divergence goes through three major evolutionary stages, i.e. ecological divergence, development of sexual isolation and the build-up of genetic incompatibility. There is theoretical consensus regarding favourable conditions required for speciation to reach its final and irreversible stage, but empirical tests remain rare. Here, we review recent research on processes of speciation, based on studies in hybrid zones between collared (Ficedula albicollis) and pied flycatchers (Ficedula hypoleuca). A major advantage of this study system is that questions concerning all three major sources of reproductive isolation and their interconnections can be addressed. We conclude that (i) ecological divergence is caused by divergence in life-history traits, (ii) females prefer mates of their own species based on differences in both plumage and song characteristics, (iii) male plumage characteristics have diverged but their song has converged in sympatry, (iv) there is genetic incompatibility in accordance with Haldane's rule, and (v) the Z-chromosome appears to be a hotspot for genes involved in sexual isolation and genetic incompatibility. We discuss how identification of the genes underlying the three major sources of reproductive isolation can be used to draw conclusions about links between the processes driving their evolution.

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.

Contrasting patterns of phenotypic plasticity in reproductive traits in two great tit (Parus major) populations

Phenotypic plasticity is an important mechanism via which populations can respond to changing environmental conditions, but we know very little about how natural populations vary with respect to plasticity. Here we use random-regression animal models to understand the multivariate phenotypic and genetic patterns of plasticity variation in two key life-history traits, laying date and clutch size, using data from long-term studies of great tits in The Netherlands (Hoge Veluwe [HV]) and UK (Wytham Woods [WW]). We show that, while population-level responses of laying date and clutch size to temperature were similar in the two populations, between-individual variation in plasticity differed markedly. Both populations showed significant variation in phenotypic plasticity (IxE) for laying date, but IxE was significantly higher in HV than in WW. There were no significant genotype-by-environment interactions (GxE) for laying date, yet differences in GxE were marginally nonsignificant between HV and WW. For clutch size, we only found significant IxE and GxE in WW but no significant difference between populations. From a multivariate perspective, plasticity in laying date was not correlated with plasticity in clutch size in either population. Our results suggest that generalizations about the form and cause of any response to changing environmental conditions across populations may be difficult.

Female heterogamety and speciation: reduced introgression of the Z chromosome between two species of nightingales

Several lines of evidence suggest that the X chromosome plays a large role in intrinsic postzygotic isolation. The role of the Z chromosome in speciation is much less understood. To explore the role of the Z chromosome in reproductive isolation, we studied nucleotide variation in two closely related bird species, the Thrush Nightingale (Luscinia luscinia) and the Common Nightingale (L. megarhynchos). These species are isolated by incomplete prezygotic isolation and female hybrid sterility. We sequenced introns of four Z-linked and eight autosomal loci and analyzed patterns of polymorphism and divergence using a divergence-with-gene flow framework. Our results suggest that the nightingale species diverged approximately 1.8 Mya. We found strong evidence of gene flow after divergence in both directions, although more introgression occurred from L. megarhynchos into L. luscinia. Gene flow was significantly higher on the autosomes than on the Z chromosome. Our results support the idea that the Z chromosome plays an important role in intrinsic postzygotic isolation in birds, although it may also contribute to the evolution of prezygotic isolation through sexual selection. This highlights the similarities in the genetic basis of reproductive isolation between organisms with heterogametic males and organisms with heterogametic females during the early stages of speciation.

Molecular cytogenetics of the california condor: evolutionary and conservation implications

Evolutionary cytogenetic comparisons involved 5 species of birds (California condor, chicken, zebra finch, collared flycatcher and black stork) belonging to divergent taxonomic orders. Seventy-four clones from a condor BAC library containing 80 genes were mapped to condor chromosomes using FISH, and 15 clones containing 16 genes were mapped to the stork Z chromosome. Maps for chicken and finch were derived from genome sequence databases, and that for flycatcher from the published literature. Gene content and gene order were highly conserved when individual condor, chicken, and zebra finch autosomes were compared, confirming that these species largely retain karyotypes close to the ancestral condition for neognathous birds. However, several differences were noted: zebra finch chromosomes 1 and 1A are homologous to condor and chicken chromosomes 1, the CHUNK1 gene appears to have transposed on condor chromosome 1, condor chromosomes 4 and 9 and zebra finch chromosomes 4 and 4A are homologous to chicken chromosome arms 4q and 4p, and novel inversions on chromosomes 4, 12 and 13 were found. Condor and stork Z chromosome gene orders are collinear and differentiated by a series of inversions/transpositions when compared to chicken, zebra finch, or flycatcher; phylogenetic analyses suggest independent rearrangement along the chicken, finch, and flycatcher lineages.

Genomic evidence for a large-Z effect

The 'large-X effect' suggests that sex chromosomes play a disproportionate role in adaptive evolution. Theoretical work indicates that this effect may be most pronounced in genetic systems with female heterogamety under both good-genes and Fisher's runaway models of sexual selection (males ZZ, females ZW). Here, I use a comparative genomic approach (alignments of several thousands of chicken-zebra finch-human-mouse-opossum orthologues) to show that avian Z-linked genes are highly overrepresented among those bird-mammalian orthologues that show evidence of accelerated rate of functional evolution in birds relative to mammals; the data suggest a twofold excess of such genes on the Z chromosome. A reciprocal analysis of genes accelerated in mammals found no evidence for an excess of X-linkage. This would be compatible with theoretical expectations for differential selection on sex-linked genes under male and female heterogamety, although the power in this case was not sufficient to statistically show that 'large-Z' was more pronounced than 'large-X'. Accelerated Z-linked genes include a variety of functional categories and are characterized by higher non-synonymous to synonymous substitution rate ratios than both accelerated autosomal and non-accelerated genes. This points at a genomic 'large-Z effect', which is widespread and of general significance for adaptive divergence in birds.