Avian sex chromosomes evolved after the divergence birds and crocodilians from their common ancestor, so they are much 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 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 used 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 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.

The phylogeny of Neoaves, the largest clade of extant birds, has remained unclear despite intense study. The difficulty associated with resolving the early branches in Neoaves is likely driven by the rapid radiation of this group. However, conflicts among studies may be exacerbated by the hypothesis that relationships are sensitive to the data type analyzed. For example, analyses of coding exons typically yield trees that place Strisores (nightjars and allies) sister to the remaining Neoaves, while analyses of non-coding data typically yield trees where Mirandornites (flamingos and grebes) is the sister of the remaining Neoaves. Our understanding of data type effects is hampered by the fact that previous analyses have used different taxa, loci, and types of non-coding data. Herein, we provide strong corroboration of the data type effects hypothesis for Neoaves by comparing trees based on coding and non-coding data derived from the same taxa and gene regions. A simple analytical method known to minimize biases due to base composition (coding nucleotides as purines and pyrimidines) resulted in coding exon data with increased congruence to the non-coding topology using concatenated analyses. These results improve our understanding of the resolution of neoavian phylogeny and point to a challenge - data type effects - that is likely to be an important factor in phylogenetic analyses of birds (and many other taxonomic groups). Using our results, we provide a summary phylogeny that identifies well-corroborated relationships and highlights specific nodes where future efforts should focus.

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 acids 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 dataset for transmembrane helices from a variety of sampled set 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.