Evolutionary genetics. Clonal inheritance of avian mitochondrial DNA

Evolution, Inheritance, and Strata Formation of the W Chromosome in Duck (Anas platyrhynchos)

The nonrecombining female-limited W chromosome is predicted to experience unique evolutionary processes. Difficulties in assembling W chromosome sequences have hindered the identification of duck W-linked sequences and their evolutionary footprint. To address this, we conducted three initial contig-level genome assemblies and developed a rigorous pipeline by which to successfully expand the W-linked data set, including 11 known genes and 24 newly identified genes. Our results indicate that the W chromosome expression may not be subject to female-specific selection; a significant convergent pattern of upregulation associated with increased female-specific selection was not detected. The genetic stability of the W chromosome is also reflected in the strong evolutionary correlation between it and the mitochondria; the complete consistency of the cladogram topology constructed from their gene sequences proves the shared maternal coevolution. By detecting the evolutionary trajectories of W-linked sequences, we have found that recombination suppression started in four distinct strata, of which three were conserved across Neognathae. Taken together, our results have revealed a unique evolutionary pattern and an independent stratum evolutionary pattern for sex chromosomes.

Phylogeny and evolution of the genus Cervus (Cervidae, Mammalia) as revealed by complete mitochondrial genomes

Mitochondrial DNA (mtDNA) lineages are recognized as important components of intra- and interspecific biodiversity, and allow to reveal colonization routes and phylogeographic structure of many taxa. Among these is the genus Cervus that is widely distributed across the Holarctic. We obtained sequences of complete mitochondrial genomes from 13 Cervus taxa and included them in global phylogenetic analyses of 71 Cervinae mitogenomes. The well-resolved phylogenetic trees confirmed Cervus to be monophyletic. Molecular dating based on several fossil calibration points revealed that ca. 2.6 Mya two main mitochondrial lineages of Cervus separated in Central Asia, the Western (including C. hanglu and C. elaphus) and the Eastern (comprising C. albirostris, C. canadensis and C. nippon). We also observed convergent changes in the composition of some mitochondrial genes in C. hanglu of the Western lineage and representatives of the Eastern lineage. Several subspecies of C. nippon and C. hanglu have accumulated a large portion of deleterious substitutions in their mitochondrial protein-coding genes, probably due to drift in the wake of decreasing population size. In contrast to previous studies, we found that the relic haplogroup B of C. elaphus was sister to all other red deer lineages and that the Middle-Eastern haplogroup E shared a common ancestor with the Balkan haplogroup C. Comparison of the mtDNA phylogenetic tree with a published nuclear genome tree may imply ancient introgressions of mtDNA between different Cervus species as well as from the common ancestor of South Asian deer, Rusa timorensis and R. unicolor, to the Cervus clade.

Genetic Diversity, Heteroplasmy, and Recombination in Mitochondrial Genomes of Daphnia pulex, Daphnia pulicaria, and Daphnia obtusa

Genetic variants of mitochondrial DNA at the individual (heteroplasmy) and population (polymorphism) levels provide insight into their roles in multiple cellular and evolutionary processes. However, owing to the paucity of genome-wide data at the within-individual and population levels, the broad patterns of these two forms of variation remain poorly understood. Here, we analyze 1,804 complete mitochondrial genome sequences from Daphnia pulex, Daphnia pulicaria, and Daphnia obtusa. Extensive heteroplasmy is observed in D. obtusa, where the high level of intraclonal divergence must have resulted from a biparental-inheritance event, and recombination in the mitochondrial genome is apparent, although perhaps not widespread. Global samples of D. pulex reveal remarkably low mitochondrial effective population sizes, <3% of those for the nuclear genome. In addition, levels of population diversity in mitochondrial and nuclear genomes are uncorrelated across populations, suggesting an idiosyncratic evolutionary history of mitochondria in D. pulex. These population-genetic features appear to be a consequence of background selection associated with highly deleterious mutations arising in the strongly linked mitochondrial genome, which is consistent with polymorphism and divergence data suggesting a predominance of strong purifying selection. Nonetheless, the fixation of mildly deleterious mutations in the mitochondrial genome also appears to be driving positive selection on genes encoded in the nuclear genome whose products are deployed in the mitochondrion.

Protein Structure, Models of Sequence Evolution, and Data Type Effects in Phylogenetic Analyses of Mitochondrial Data: A Case Study in Birds

Phylogenomic analyses have revolutionized the study of biodiversity, but they have revealed that estimated tree topologies can depend, at least in part, on the subset of the genome that is analyzed. For example, estimates of trees for avian orders differ if protein-coding or non-coding data are analyzed. The bird tree is a good study system because the historical signal for relationships among orders is very weak, which should permit subtle non-historical signals to be identified, while monophyly of orders is strongly corroborated, allowing identification of strong non-historical signals. Hydrophobic amino acids in mitochondrially-encoded proteins, which are expected to be found in transmembrane helices, have been hypothesized to be associated with non-historical signals. We tested this hypothesis by comparing the evolution of transmembrane helices and extramembrane segments of mitochondrial proteins from 420 bird species, sampled from most avian orders. We estimated amino acid exchangeabilities for both structural environments and assessed the performance of phylogenetic analysis using each data type. We compared those relative exchangeabilities with values calculated using a substitution matrix for transmembrane helices estimated using a variety of nuclear- and mitochondrially-encoded proteins, allowing us to compare the bird-specific mitochondrial models with a general model of transmembrane protein evolution. To complement our amino acid analyses, we examined the impact of protein structure on patterns of nucleotide evolution. Models of transmembrane and extramembrane sequence evolution for amino acids and nucleotides exhibited striking differences, but there was no evidence for strong topological data type effects. However, incorporating protein structure into analyses of mitochondrially-encoded proteins improved model fit. Thus, we believe that considering protein structure will improve analyses of mitogenomic data, both in birds and in other taxa.

When good mitochondria go bad: Cyto-nuclear discordance in landfowl (Aves: Galliformes)

Mitochondrial sequences were among the first molecular data collected for phylogenetic studies and they are plentiful in DNA sequence archives. However, the future value of mitogenomic data in phylogenetics is uncertain, because its phylogenetic signal sometimes conflicts with that of the nuclear genome. A thorough understanding of the causes and prevalence of cyto-nuclear discordance would aid in reconciling different results owing to sequence data type, and provide a framework for interpreting megaphylogenies when taxa which lack substantial nuclear data are placed using mitochondrial data. Here, we examine the prevalence and possible causes of cyto-nuclear discordance in the landfowl (Aves: Galliformes), leveraging 47 new mitogenomes assembled from off-target reads recovered as part of a target-capture study. We evaluated two hypotheses, that cyto-nuclear discordance is "genuine" and a result of biological processes such as incomplete lineage sorting or introgression, and that cyto-nuclear discordance is an artifact of inaccurate mitochondrial tree estimation (the "inaccurate estimation" hypothesis). We identified seven well-supported topological differences between the mitogenomic tree and trees based on nuclear data. These well-supported topological differences were robust to model selection. An examination of sites suggests these differences were driven by small number of sites, particularly from third-codon positions, suggesting that they were not confounded by convergent directional selection. Hence, the hypothesis of genuine discordance was supported.

Evaluating evidence of mitonuclear incompatibilities with the sex chromosomes in an avian hybrid zone

The exploration of hybrid zones and the intergenomic conflicts exposed through hybridization provide windows into the processes of divergence and speciation. Sex chromosomes and mitonuclear incompatibilities have strong associations with the genetics of hybrid dysfunction. In ZW sex-determining systems, maternal co-inheritance of the mitochondrial and W chromosomes immediately exposes incompatibilities between these maternal contributions of one species and the Z chromosome of another. We analyze mitochondrial and Z chromosome admixture in the long-tailed finch (Poephila acuticauda) of Australia, where hybridizing subspecies differ prominently in Z chromosome genotype and in bill color, yet the respective centers of geographic admixture for these two traits are offset by 350 km. We report two well-defined mitochondrial clades that diverged ∼0.5 million years ago. Mitochondrial contact is geographically co-located within a hybrid zone of Z chromosome admixture and is displaced from bill color admixture by nearly 400 km. Consistent with Haldane's rule expectations, hybrid zone females are significantly less likely than males to carry an admixed Z chromosome or have mismatched Z-mitochondrial genotypes. Furthermore, there are significantly fewer than expected mitonuclear mismatches in hybrid zone females and paternal backcross males. Results suggest a potential for mitonuclear/sex chromosome incompatibilities in the emergence of reproductive isolation in this system.

New view on the organization and evolution of Palaeognathae mitogenomes poses the question on the ancestral gene rearrangement in Aves

BACKGROUND

Bird mitogenomes differ from other vertebrates in gene rearrangement. The most common avian gene order, identified first in Gallus gallus, is considered ancestral for all Aves. However, other rearrangements including a duplicated control region and neighboring genes have been reported in many representatives of avian orders. The repeated regions can be easily overlooked due to inappropriate DNA amplification or genome sequencing. This raises a question about the actual prevalence of mitogenomic duplications and the validity of the current view on the avian mitogenome evolution. In this context, Palaeognathae is especially interesting because is sister to all other living birds, i.e. Neognathae. So far, a unique duplicated region has been found in one palaeognath mitogenome, that of Eudromia elegans.

RESULTS

Therefore, we applied an appropriate PCR strategy to look for omitted duplications in other palaeognaths. The analyses revealed the duplicated control regions with adjacent genes in Crypturellus, Rhea and Struthio as well as ND6 pseudogene in three moas. The copies are very similar and were subjected to concerted evolution. Mapping the presence and absence of duplication onto the Palaeognathae phylogeny indicates that the duplication was an ancestral state for this avian group. This feature was inherited by early diverged lineages and lost two times in others. Comparison of incongruent phylogenetic trees based on mitochondrial and nuclear sequences showed that two variants of mitogenomes could exist in the evolution of palaeognaths. Data collected for other avian mitogenomes revealed that the last common ancestor of all birds and early diverging lineages of Neoaves could also possess the mitogenomic duplication.

CONCLUSIONS

The duplicated control regions with adjacent genes are more common in avian mitochondrial genomes than it was previously thought. These two regions could increase effectiveness of replication and transcription as well as the number of replicating mitogenomes per organelle. In consequence, energy production by mitochondria may be also more efficient. However, further physiological and molecular analyses are necessary to assess the potential selective advantages of the mitogenome duplications.

Resolving Phylogenetic Relationships within Passeriformes Based on Mitochondrial Genes and Inferring the Evolution of Their Mitogenomes in Terms of Duplications

Mitochondrial genes are placed on one molecule, which implies that they should carry consistent phylogenetic information. Following this advantage, we present a well-supported phylogeny based on mitochondrial genomes from almost 300 representatives of Passeriformes, the most numerous and differentiated Aves order. The analyses resolved the phylogenetic position of paraphyletic Basal and Transitional Oscines. Passerida occurred divided into two groups, one containing Paroidea and Sylvioidea, whereas the other, Passeroidea and Muscicapoidea. Analyses of mitogenomes showed four types of rearrangements including a duplicated control region (CR) with adjacent genes. Mapping the presence and absence of duplications onto the phylogenetic tree revealed that the duplication was the ancestral state for passerines and was maintained in early diverged lineages. Next, the duplication could be lost and occurred independently at least four times according to the most parsimonious scenario. In some lineages, two CR copies have been inherited from an ancient duplication and highly diverged, whereas in others, the second copy became similar to the first one due to concerted evolution. The second CR copies accumulated over twice as many substitutions as the first ones. However, the second CRs were not completely eliminated and were retained for a long time, which suggests that both regions can fulfill an important role in mitogenomes. Phylogenetic analyses based on CR sequences subjected to the complex evolution can produce tree topologies inconsistent with real evolutionary relationships between species. Passerines with two CRs showed a higher metabolic rate in relation to their body mass.

Sex chromosomes and speciation in birds and other ZW systems

Theory and empirical patterns suggest a disproportionate role for sex chromosomes in evolution and speciation. Focusing on ZW sex determination (females ZW, males ZZ; the system in birds, many snakes, and lepidopterans), I review how evolutionary dynamics are expected to differ between the Z, W and the autosomes, discuss how these differences may lead to a greater role of the sex chromosomes in speciation and use data from birds to compare relative evolutionary rates of sex chromosomes and autosomes. Neutral mutations, partially or completely recessive beneficial mutations, and deleterious mutations under many conditions are expected to accumulate faster on the Z than on autosomes. Sexually antagonistic polymorphisms are expected to arise on the Z, raising the possibility of the spread of preference alleles. The faster accumulation of many types of mutations and the potential for complex evolutionary dynamics of sexually antagonistic traits and preferences contribute to a role for the Z chromosome in speciation. A quantitative comparison among a wide variety of bird species shows that the Z tends to have less within-population diversity and greater between-species differentiation than the autosomes, likely due to both adaptive evolution and a greater rate of fixation of deleterious alleles. The W chromosome also shows strong potential to be involved in speciation, in part because of its co-inheritance with the mitochondrial genome. While theory and empirical evidence suggest a disproportionate role for sex chromosomes in speciation, the importance of sex chromosomes is moderated by their small size compared to the whole genome.

Evolutionary analysis of the female-specific avian W chromosome

The typically repetitive nature of the sex-limited chromosome means that it is often excluded from or poorly covered in genome assemblies, hindering studies of evolutionary and population genomic processes in non-recombining chromosomes. Here, we present a draft assembly of the non-recombining region of the collared flycatcher W chromosome, containing 46 genes without evidence of female-specific functional differentiation. Survival of genes during W chromosome degeneration has been highly non-random and expression data suggest that this can be attributed to selection for maintaining gene dose and ancestral expression levels of essential genes. Re-sequencing of large population samples revealed dramatically reduced levels of within-species diversity and elevated rates of between-species differentiation (lineage sorting), consistent with low effective population size. Concordance between W chromosome and mitochondrial DNA phylogenetic trees demonstrates evolutionary stable matrilineal inheritance of this nuclear–cytonuclear pair of chromosomes. Our results show both commonalities and differences between W chromosome and Y chromosome evolution.

The evolution of non-recombining chromosomes is poorly understood. Here, the authors sequence the collared flycatcher female-specific W chromosome and show nonrandom survival of genes during W chromosome degeneration which is due to selection for maintaining gene dose and expression levels of essential genes.

Does more sequence data improve estimates of galliform phylogeny? Analyses of a rapid radiation using a complete data matrix

The resolution of rapid evolutionary radiations or "bushes" in the tree of life has been one of the most difficult and interesting problems in phylogenetics. The avian order Galliformes appears to have undergone several rapid radiations that have limited the resolution of prior studies and obscured the position of taxa important both agriculturally and as model systems (chicken, turkey, Japanese quail). Here we present analyses of a multi-locus data matrix comprising over 15,000 sites, primarily from nuclear introns but also including three mitochondrial regions, from 46 galliform taxa with all gene regions sampled for all taxa. The increased sampling of unlinked nuclear genes provided strong bootstrap support for all but a small number of relationships. Coalescent-based methods to combine individual gene trees and analyses of datasets that are independent of published data indicated that this well-supported topology is likely to reflect the galliform species tree. The inclusion or exclusion of mitochondrial data had a limited impact upon analyses upon analyses using either concatenated data or multispecies coalescent methods. Some of the key phylogenetic findings include support for a second major clade within the core phasianids that includes the chicken and Japanese quail and clarification of the phylogenetic relationships of turkey. Jackknifed datasets suggested that there is an advantage to sampling many independent regions across the genome rather than obtaining long sequences for a small number of loci, possibly reflecting the differences among gene trees that differ due to incomplete lineage sorting. Despite the novel insights we obtained using this increased sampling of gene regions, some nodes remain unresolved, likely due to periods of rapid diversification. Resolving these remaining groups will likely require sequencing a very large number of gene regions, but our analyses now appear to support a robust backbone for this order.

The erratic mitochondrial clock: variations of mutation rate, not population size, affect mtDNA diversity across birds and mammals

BACKGROUND

During the last ten years, major advances have been made in characterizing and understanding the evolution of mitochondrial DNA, the most popular marker of molecular biodiversity. Several important results were recently reported using mammals as model organisms, including (i) the absence of relationship between mitochondrial DNA diversity and life-history or ecological variables, (ii) the absence of prominent adaptive selection, contrary to what was found in invertebrates, and (iii) the unexpectedly large variation in neutral substitution rate among lineages, revealing a possible link with species maximal longevity. We propose to challenge these results thanks to the bird/mammal comparison. Direct estimates of population size are available in birds, and this group presents striking life-history trait differences with mammals (higher mass-specific metabolic rate and longevity). These properties make birds the ideal model to directly test for population size effects, and to discriminate between competing hypotheses about the causes of substitution rate variation.

RESULTS

A phylogenetic analysis of cytochrome b third-codon position confirms that the mitochondrial DNA mutation rate is quite variable in birds, passerines being the fastest evolving order. On average, mitochondrial DNA evolves slower in birds than in mammals of similar body size. This result is in agreement with the longevity hypothesis, and contradicts the hypothesis of a metabolic rate-dependent mutation rate. Birds show no footprint of adaptive selection on cytochrome b evolutionary patterns, but no link between direct estimates of population size and cytochrome b diversity. The mutation rate is the best predictor we have of within-species mitochondrial diversity in birds. It partly explains the differences in mitochondrial DNA diversity patterns observed between mammals and birds, previously interpreted as reflecting Hill-Robertson interferences with the W chromosome.

CONCLUSION

Mitochondrial DNA diversity patterns in birds are strongly influenced by the wide, unexpected variation of mutation rate across species. From a fundamental point of view, these results are strongly consistent with a relationship between species maximal longevity and mitochondrial mutation rate, in agreement with the mitochondrial theory of ageing. Form an applied point of view, this study reinforces and extends the message of caution previously expressed for mammals: mitochondrial data tell nothing about species population sizes, and strongly depart the molecular clock assumption.

Genetic evidence from Indian red jungle fowl corroborates multiple domestication of modern day chicken

BACKGROUND

Domestication of chicken is believed to have occurred in Southeast Asia, especially in Indus valley. However, non-inclusion of Indian red jungle fowl (RJF), Gallus gallus murghi in previous studies has left a big gap in understanding the relationship of this major group of birds. In the present study, we addressed this issue by analyzing 76 Indian birds that included 56 G. g. murghi (RJF), 16 G. g. domesticus (domestic chicken) and 4 G. sonneratii (Grey JF) using both microsatellite markers and mitochondrial D-loop sequences. We also compared the D-loop sequences of Indian birds with those of 779 birds obtained from GenBank.

RESULTS

Microsatellite marker analyses of Indian birds indicated an average FST of 0.126 within G. g. murghi, and 0.154 within G. g. domesticus while it was more than 0.2 between the two groups. The microsatellite-based phylogenetic trees showed a clear separation of G. g. domesticus from G. g. murghi, and G. sonneratii. Mitochondrial DNA based mismatch distribution analyses showed a lower Harpending's raggedness index in both G. g. murghi (0.001515) and in Indian G. g. domesticus (0.0149) birds indicating population expansion. When meta analysis of global populations of 855 birds was carried out using median joining haplotype network, 43 Indian birds of G. g. domesticus (19 haplotypes) were distributed throughout the network sharing haplotypes with the RJFs of different origins.

CONCLUSION

Our results suggest that the domestication of chicken has occurred independently in different locations of Asia including India. We found evidence for domestication of Indian birds from G. g. spadiceus and G. g. gallus as well as from G. g. murghi, corroborating multiple domestication of Indian and other domestic chicken. In contrast to the commonly held view that RJF and domestic birds hybridize in nature, the present study shows that G. g. murghi is relatively pure. Further, the study also suggested that the chicken populations have undergone population expansion, especially in the Indus valley.

Genetic evidence from Indian red jungle fowl corroborates multiple domestication of modern day chicken

Background

Domestication of chicken is believed to have occurred in Southeast Asia, especially in Indus valley. However, non-inclusion of Indian red jungle fowl (RJF), Gallus gallus murghi in previous studies has left a big gap in understanding the relationship of this major group of birds. In the present study, we addressed this issue by analyzing 76 Indian birds that included 56 G. g. murghi (RJF), 16 G. g. domesticus (domestic chicken) and 4 G. sonneratii (Grey JF) using both microsatellite markers and mitochondrial D-loop sequences. We also compared the D-loop sequences of Indian birds with those of 779 birds obtained from GenBank.

Results

Microsatellite marker analyses of Indian birds indicated an average F_ST of 0.126 within G. g. murghi, and 0.154 within G. g. domesticus while it was more than 0.2 between the two groups. The microsatellite-based phylogenetic trees showed a clear separation of G. g. domesticus from G. g. murghi, and G. sonneratii. Mitochondrial DNA based mismatch distribution analyses showed a lower Harpending's raggedness index in both G. g. murghi (0.001515) and in Indian G. g. domesticus (0.0149) birds indicating population expansion. When meta analysis of global populations of 855 birds was carried out using median joining haplotype network, 43 Indian birds of G. g. domesticus (19 haplotypes) were distributed throughout the network sharing haplotypes with the RJFs of different origins.

Conclusion

Our results suggest that the domestication of chicken has occurred independently in different locations of Asia including India. We found evidence for domestication of Indian birds from G. g. spadiceus and G. g. gallus as well as from G. g. murghi, corroborating multiple domestication of Indian and other domestic chicken. In contrast to the commonly held view that RJF and domestic birds hybridize in nature, the present study shows that G. g. murghi is relatively pure. Further, the study also suggested that the chicken populations have undergone population expansion, especially in the Indus valley.

Association of MT-ND5 gene variation with mitochondrial respiratory control ratio and NADH dehydrogenase activity in Tibet chicken embryos

NADH dehydrogenase (complex I) couples the oxidation of NADH for the reduction of ubiquinone with the generation of a proton gradient that can be used for the synthesis of ATP. We have found a missense mutation in the MT-ND5 subunit of NADH dehydrogenase in the Tibet chicken breed. In the present study, the mitochondrial respiratory control ratio (RCR) and NADH dehydrogenase activity in Tibet chicken embryonic brain with different genotypes were measured. Significant differences between animals carrying mitochondria with the EF493865.1:m.1627A vs. EF493865.1:m.1627C alleles were observed for RCR and enzyme activity.

Low mitochondrial variability in birds may indicate Hill–Robertson effects on the W chromosome

Interference among loci subject to selection (the Hill-Robertson effect) may considerably reduce levels of adaptation and variability in genomic regions that lack recombination. Y- or W chromosomes are particularly vulnerable to such effects, since they represent large, non-recombining blocks of genetic material. In birds, the W chromosome and mitochondrial genomes are both maternally transmitted, and hence fail to recombine with each other, whereas in mammals the Y chromosome is paternally transmitted. We show here that mitochondrial DNA sequence diversity is reduced in non-ratite birds compared with mammals. After considering possible confounding factors, such as differences in generation times, mutation rates and demography, we conclude that Hill-Robertson effects associated with the W chromosome provide the most likely explanation for this difference.

Conflicting patterns of mitochondrial and nuclear DNA diversity in Phylloscopus warblers

Molecular variation is often used to infer the demographic history of species, but sometimes the complexity of species history can make such inference difficult. The willow warbler, Phylloscopus trochilus, shows substantially less geographical variation than the chiffchaff, Phylloscopus collybita, both in morphology and in mitochondrial DNA (mtDNA) divergence. We therefore predicted that the willow warbler should harbour less nuclear DNA diversity than the chiffchaff. We analysed sequence data obtained from multiple samples of willow warblers and chiffchaffs for the mtDNA cytochrome b gene and four nuclear genes. We confirmed that the mtDNA diversity among willow warblers is low (pi = 0.0021). Sequence data from three nuclear genes (CHD-Z, AFLP-WW1 and MC1R) not linked to the mitochondria demonstrated unexpectedly high nucleotide diversity (pi values of 0.0172, 0.0141 and 0.0038) in the willow warbler, on average higher than the nucleotide diversity for the chiffchaff (pi values of 0.0025, 0.0017 and 0.0139). In willow warblers, Tajima's D analyses showed that the mtDNA diversity, but not the nuclear DNA diversity, has been reduced relative to the neutral expectation of molecular evolution, suggesting the action of a selective sweep affecting the maternally inherited genes. The large nuclear diversity seen within willow warblers is not compatible with processes of neutral evolution occurring in a population with a constant population size, unless the long-term effective population size has been very large (N(e) > 10(6)). We suggest that the contrasting patterns of genetic diversity in the willow warbler may reflect a more complex evolutionary history, possibly including historical demographic fluctuations or historical male-biased introgression of nuclear genes from a differentiated population of Phylloscopus warblers.

Mutation pressure and the evolution of organelle genomic architecture

The nuclear genomes of multicellular animals and plants contain large amounts of noncoding DNA, the disadvantages of which can be too weak to be effectively countered by selection in lineages with reduced effective population sizes. In contrast, the organelle genomes of these two lineages evolved to opposite ends of the spectrum of genomic complexity, despite similar effective population sizes. This pattern and other puzzling aspects of organelle evolution appear to be consequences of differences in organelle mutation rates. These observations provide support for the hypothesis that the fundamental features of genome evolution are largely defined by the relative power of two nonadaptive forces: random genetic drift and mutation pressure.

Fast accumulation of nonsynonymous mutations on the female-specific W chromosome in birds

Following cessation of recombination during sex chromosome evolution, the nonrecombining sex chromosome is affected by a number of degenerative forces, possibly resulting in the fixation of deleterious mutations. This might take place because of weak selection against recessive or partly recessive deleterious mutations due to permanent heterozygosity of nonrecombining chromosomes. Furthermore, population genetic processes, such as selective sweeps, background selection, and Muller's ratchet, result in a reduction in Ne, which increase the likelihood of fixation of deleterious mutations. Theory thus predicts that nonrecombining genes should show increased levels of nonsynonymous (dN) to synonymous substitutions (dS). We tested this in an avian system by estimating the ratio between dN and dS in six gametologous gene pairs located on the Z chromosome and the nonrecombining, female-specific W chromosome. In comparisons, we found a significantly higher dN/dS ratio for the W-linked than the Z-linked copy in three of the investigated genes. In a concatenated alignment of all six genes, the dN/dS ratio was six times higher for W-linked than Z-linked genes. By using human and mouse as outgroup in maximum likelihood analyses, W-linked genes were found to evolve differently compared with their Z-linked gametologues and outgroup sequences. This seems not to be a consequence of functional diversification because d(N)/d(S) ratios between gametologous gene copies were consistently low. We conclude that deleterious mutations are accumulating at a high rate on the avian W chromosome, probably as a result of the lack of recombination in this female-specific chromosome.

Chicken W: a genetically uniform chromosome in a highly variable genome

The Y chromosome of organisms with male heterogamety is expected to show reduced levels of genetic diversity, because the effective population size is one-fourth that of autosomes. However, studies in mammals, flies, and plants show that Y chromosome diversity is lower than expected even when differences in effective population size are taken into account. This may be explained by skewed reproductive success among males, leading to low male effective population size, or by a strong role of selection in shaping levels of nucleotide diversity in nonrecombining chromosomes. We tested these hypotheses in a system with female heterogamety by estimating nucleotide diversity in the female-specific W chromosome of the domestic chicken by resequencing of 7,643 base pairs in 47 birds from 10 highly divergent breeds. The screening revealed only one single segregating site, which is in sharp contrast to our previous observation, using a similar panel of birds of, on average, one segregating site every 39 base pairs in autosomal sequence. When taking sex-specific mutation rates and differences in effective population size into account, the observed degree of W chromosome polymorphism is 28-fold lower than expected for the frequency of segregating sites and 13-fold lower than expected for estimates of nucleotide diversity (autosomes, 6.5 x 10(-3); W, 7.0 x 10(-5)). We note that selection is the only factor that can explain the reduced diversity in the sex-limited chromosome irrespective of mode of reproduction or whether there is male or female heterogamety. Reduced variability in female-specific W chromosomes is not easily explained by sexual selection.

Do avian mitochondria recombine?

The dogma of strict maternal inheritance of mitochondria is now being tested with population genetics methods on sequence data from many species. In this study we investigated whether recombination occurs in the mitochondria of the blue tit ( Parus caeruleus) by studying polymorphisms in the mitochondrial control region and in a recently identified (A)(n) microsatellite on the W chromosome. The female heterogamety of avian sex chromosomes allows a test of whether mitochondrial recombination affects genealogical inference by comparison of mitochondrial and W-linked sequence variation. There is no discrepancy between mitochondrial and W-linked genealogies in blue tits, consistent with no recombination. We also analyzed mitochondrial sequence variation in both blue tits and peregrine falcons ( Falco peregrinus) using a coalescent-based approach which accounts for recurrent mutation; in neither bird species did we find evidence of recombination. We conclude that it is unlikely that mitochondrial recombination has large effects on mitochondrial genetic variability in birds.

No association between mitochondrial DNA haplotypes and a female-limited mimicry phenotype in Papilio glaucus

Alternative alleles at a locus on the W chromosome of Papilio glaucus (causing dark or yellow wing colors, respectively) underlie a female-limited mimicry polymorphism thought to be maintained by balancing selection. In species with heterogametic females (i.e., the ZZ-male/ZW-female sex chromosome system), the mitochondrial DNA and the W chromosome are genetically linked because they are both maternally transmitted. We investigate the association of COI and COII mitochondrial DNA haplotypes with alternative W-linked phenotypes. Surprisingly, we find no congruence between mitochondrial DNA genealogies and inferred W-linked color alleles in P. glaucus. Using a maximum-likelihood phylogenetic approach, we reject the hypothesis of monophyly for dark-morph mitochondrial DNA lineages, even in the presence of putative low-frequency mimicry suppressor alleles or alternative melanizing factors. The most likely genealogical tree topologies assume more than one exchange event between mitochondrial DNA cytotype and the W-linked color morph. These results suggest that there is either paternal leakage of mitochondrial DNA or that more than two W-linked alleles underlie the alternative color morphs. Using data from an additional mitochondrial DNA locus, ND5, we show that pairwise linkage disequilibrium decays with physical distance between polymorphic sites. This finding suggests that genetic exchanges between maternal and paternal mitochondrial DNAs may have contributed to the lack of association we observe between phenotype and genotype.

Sex chromosome evolution and speciation in Ficedula flycatchers

Speciation is the combination of evolutionary processes that leads to the reproductive isolation of different populations. We investigate the significance of sex-chromosome evolution on the development of post- and prezygotic isolation in two naturally hybridizing Ficedula flycatcher species. Applying a tag-array-based mini-sequencing assay to genotype single nucleotide polymorphisms (SNPs) and interspecific substitutions, we demonstrate rather extensive hybridization and backcrossing in sympatry. However, gene flow across the partial postzygotic barrier (introgression) is almost exclusively restricted to autosomal loci, suggesting strong selection against introgression of sex-linked genes. In addition to this partial postzygotic barrier, character displacement of male plumage characteristics has previously been shown to reinforce prezygotic isolation in these birds. We show that male plumage traits involved in reinforcing prezygotic isolation are sex linked. These results suggest a major role of sex-chromosome evolution in mediating post- and prezygotic barriers to gene flow and point to a causal link in the development of the two forms of reproductive isolation.