Palaeognathae represent one of the two basal lineages in modern birds, and comprise the volant (flighted) tinamous and the flightless ratites. Resolving palaeognath phylogenetic relationships has historically proved difficult, and short internal branches separating major palaeognath lineages in previous molecular phylogenies suggest that extensive incomplete lineage sorting (ILS) might have accompanied a rapid ancient divergence. Here, we investigate palaeognath relationships using genome-wide data sets of three types of noncoding nuclear markers, together totalling 20,850 loci and over 41 million base pairs of aligned sequence data. We recover a fully resolved topology placing rheas as the sister to kiwi and emu + cassowary that is congruent across marker types for two species tree methods (MP-EST and ASTRAL-II). This topology is corroborated by patterns of insertions for 4,274 CR1 retroelements identified from multi-species whole genome screening, and is robustly supported by phylogenomic subsampling analyses, with MP-EST demonstrating particularly consistent performance across subsampling replicates as compared to ASTRAL. In contrast, analyses of concatenated data supermatrices recover rheas as the sister to all other non-ostrich palaeognaths, an alternative that lacks retroelement support and shows inconsistent behavior under subsampling approaches. While statistically supporting the species tree topology, conflicting patterns of retroelement insertions also occur and imply high amounts of ILS across short successive internal branches, consistent with observed patterns of gene tree heterogeneity. Coalescent simulations indicate that the majority of observed topological incongruence among gene trees is consistent with coalescent variation rather than arising from gene tree estimation error alone, and estimated branch lengths for short successive internodes in the inferred species tree fall within the theoretical range encompassing the anomaly zone. Distributions of empirical gene trees confirm that the most common gene tree topology for each marker type differs from the species tree, signifying the existence of an empirical anomaly zone in palaeognaths.
Supplementary data: species codes
species_codes.txt
Supplementary data: sequence alignments
alignments.tar.gz
Supplementary data: phylogenomic subsampling
phylogenomic_subsampling.tar.gz
Supplementary data: RAxML bestML gene trees
RAxML_bestML_gene_trees.tar.gz
Supplementary data: RAxML CNEE bootstraps part1
RAxML_CNEE_bootstraps_part1.tar.gz
Supplementary data: RAxML CNEE bootstraps part2
RAxML_CNEE_bootstraps_part2.tar.gz
Supplementary data: RAxML CNEE bootstraps part3
RAxML_CNEE_bootstraps_part3.tar.gz
Supplementary data: RAxML intron bootstraps
RAxML_intron_bootstraps.tar.gz
Supplementary data:RAxML UCE bootstraps
RAxML_UCE_bootstraps.tar.gz
Supplementary data: coalescent simulations
coalescent_simulations.tar.gz
Supplementary data: CR1 retroelement insertions final dataset
CR1_retroelement_insertions_final_dataset.tar.gz
Supplementary data: rooted_triplets
rooted_triplets.tar.gz
Supplementary data: species_trees
species_trees.tar.gz
Supplementary data: CR1 retroelement insertions bioinformatic pipeline part1
CR1_retroelement_insertions_bioinformatic_pipeline_part1.tar.gz
Supplementary data: CR1 retroelement insertions bioinformatic pipeline part2
CR1_retroelement_insertions_bioinformatic_pipeline_part2.tar.gz
Supplementary data: CR1 retroelement insertions bioinformatic pipeline part3
CR1_retroelement_insertions_bioinformatic_pipeline_part3.tar.gz
Supplementary data: CR1 retroelement insertions bioinformatic pipeline part4
CR1_retroelement_insertions_bioinformatic_pipeline_part4.tar.gz
Supplementary data: CR1 retroelement insertions bioinformatic pipeline part5
CR1_retroelement_insertions_bioinformatic_pipeline_part5.tar.gz
Supplementary data: CR1 retroelement insertions bioinformatic pipeline part6
CR1_retroelement_insertions_bioinformatic_pipeline_part6.tar.gz
Supplementary data: CR1 retroelement insertions bioinformatic pipeline part7
CR1_retroelement_insertions_bioinformatic_pipeline_part7.tar.gz
Supplementary data: CR1 retroelement insertions bioinformatic pipeline part8
CR1_retroelement_insertions_bioinformatic_pipeline_part8.tar.gz
Supplementary data: CR1 retroelement insertions bioinformatic pipeline part9
CR1_retroelement_insertions_bioinformatic_pipeline_part9.tar.gz
Supplementary data: CR1 retroelement insertions bioinformatic pipeline part10
CR1_retroelement_insertions_bioinformatic_pipeline_part10.tar.gz
Supplementary data: CR1 retroelement insertions bioinformatic pipeline part11
CR1_retroelement_insertions_bioinformatic_pipeline_part11.tar.gz
Supplementary data: CR1 retroelement insertions bioinformatic pipeline part12
CR1_retroelement_insertions_bioinformatic_pipeline_part12.tar.gz
Supplementary data: CR1 retroelement insertions bioinformatic pipeline part13
CR1_retroelement_insertions_bioinformatic_pipeline_part13.tar.gz
Supplementary data: CR1 retroelement insertions bioinformatic pipeline part14
CR1_retroelement_insertions_bioinformatic_pipeline_part14.tar.gz
Supplementary data: CR1 retroelement insertions bioinformatic pipeline part15
CR1_retroelement_insertions_bioinformatic_pipeline_part15.tar.gz
Supplementary Figure S1
Supplementary Figure S1 Palaeognath phylogenetic relationships inferred with MP-EST (a-d) and ASTRAL (e-h) coalescent-based species tree methods and maximum-likelihood inference of fully partitioned concatenated alignments with ExaML (i-l). Topologies were inferred from 12,676 conserved non-exonic elements (CNEEs; parts a,e,i), 5,016 introns (parts b,f,j), 3,158 ultraconserved elements (UCEs; parts c,g,k), and the total evidence data set of 20,850 loci combined across marker types (TENT; parts d,h,l). Bootstrap supports are drawn for clades with < 100% support.
Suppl_FigS1.pdf
Supplementary Figure S2
Supplementary Figure S2 Support for alternative hypotheses of the sister group to emu + cassowary from phylogenomic subsampling using MP-EST (a-c), ASTRAL (d-f), and ExaML (g-i). Plots display the mean bootstrap support for each hypothesis from 10 replicates of randomly sampled loci within each data set size category (e.g. 50–3000 loci, shown on x-axis).
Suppl_FigS2.pdf
Supplementary Figure S3
Supplementary Figure S3 Heatmap of the number of replicates from phylogenomic subsampling that support each alternative hypothesis for the sister group to emu + cassowary at a minimum bootstrap support of 70% (a-i) and 90% (j-r), with results from each analysis method displayed in separate panels (e.g. MP-EST a-c & j-l, ASTRAL d-f & m-o, ExaML g-I & p-r). Rows labelled H1-H4 within each panel correspond to the four alternative hypotheses outlined in the legend. Columns within each panel labelled 50-3000 indicate the number of loci subsampled at random with replacement from all loci for each marker type. Coloring of cells indicates the number of replicates (of 10 in total for each data set size category) that support each hypothesis at the given bootstrap cut off, corresponding to the colouring scheme outlined in the legend.
Suppl_FigS3.pdf
Supplementary Figure S4
Supplementary Figure S4 Support for the four most common gene tree topologies for each marker type. Median bootstrap support (a-c) and the median number of substitutions under a parsimony criterion (d-f) for branches in gene tree topologies that conflict with the species tree are shown for each marker type. Letters below each group refer to topologies shown in (j), with numbers referring to branches labelled on those topologies. Reference lines in (a-c) are drawn at 50% bootstrap support. (g-i) Violin plots of ΔAIC for loci belonging to each AGT category for CNEEs (part g), introns (part h), and UCEs (part i). Values show the difference in AIC for sequence alignments given the majority rule extended consensus gene tree topology, relative to the AIC when sequence alignments are constrained to the species tree topology. A reference line is drawn at ΔAIC(gene tree–species tree)= -2, with values beneath this cut off indicating substantially stronger support in favor of the gene tree topology (Burnham and Anderson 2002).
Suppl_FigS4.pdf
Supplementary Figure S5
Supplementary Figure S5 Summary measures of aligned sequence data for each marker type (CNEEs, introns, UCEs), showing the total aligned sequence length per locus (a-c), variable alignment columns as a percent of the total alignment length (d-f), and the percent of parsimony informative alignment columns per locus (g-i).
Suppl_FigS5.pdf
Supplementary Figure S6
Supplementary Figure S6 Pairs of short successive internal branches are consistent with expectations for the anomaly zone. Species tree topologies for CNEEs (a), introns (b), and UCEs (c) are shown with internal branch lengths in coalescent units estimated from best maximum likelihood gene trees with ASTRAL. Terminal branch lengths are uninformative and are drawn as a constant value across taxa. Coalescent branch lengths for all pairs of branches (x and y) are given in (d), with a(x) calculated following Equation 4 from Degnan and Rosenberg (2006). Anomalous gene trees are expected when y < a(x). Clades fulfilling this anomaly zone criterion are shaded in (d), with the corresponding branches indicated in a-c.
Suppl_FigS6.pdf
Supplementary Tables
Supplementary_Tables.pdf