A core question in evolutionary biology is whether convergent phenotypic evolution is driven by convergent molecular changes in proteins or regulatory regions. We combined phylogenomic, developmental, and epigenomic analysis of 11 new genomes of paleognathous birds, including an extinct moa, to show that convergent evolution of regulatory regions, more so than protein-coding genes, is prevalent among developmental pathways associated with independent losses of flight. A Bayesian analysis of 284,001 conserved noncoding elements, 60,665 of which are corroborated as enhancers by open chromatin states during development, identified 2355 independent accelerations along lineages of flightless paleognaths, with functional consequences for driving gene expression in the developing forelimb. Our results suggest that the genomic landscape associated with morphological convergence in ratites has a substantial shared regulatory component.

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 totaling 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 4274 CR1 retroelements identified from multispecies 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 nonostrich 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 and topology tests 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.

24Palaeognathae represent one of the two basal lineages in modern birds, and comprise the 25 volant (flighted) tinamous and the flightless ratites. Resolving palaeognath phylogenetic 26 relationships has historically proved difficult, and short internal branches separating major 27 palaeognath lineages in previous molecular phylogenies suggest that extensive incomplete 28 lineage sorting (ILS) might have accompanied a rapid ancient divergence. Here, we investigate 29 palaeognath relationships using genome-wide data sets of three types of noncoding nuclear 30 markers, together totalling 20,850 loci and over 41 million base pairs of aligned sequence data. 31We recover a fully resolved topology placing rheas as the sister to kiwi and emu + cassowary 32 that is congruent across marker types for two species tree methods (MP-EST and ASTRAL-II). 33This topology is corroborated by patterns of insertions for 4,274 CR1 retroelements identified 34 from multi-species whole genome screening, and is robustly supported by phylogenomic 35 subsampling analyses, with MP-EST demonstrating particularly consistent performance across 36 subsampling replicates as compared to ASTRAL. In contrast, analyses of concatenated data 37 supermatrices recover rheas as the sister to all other non-ostrich palaeognaths, an alternative that 38 lacks retroelement support and shows inconsistent behavior under subsampling approaches. 39 While statistically supporting the species tree topology, conflicting patterns of retroelement 40 insertions also occur and imply high amounts of ILS across short successive internal branches, 41 consistent with observed patterns of gene tree heterogeneity. Coalescent simulations indicate 42 that the majority of observed topological incongruence among gene trees is consistent with 43 coalescent variation rather than arising from gene tree estimation error alone, and estimated 44 branch lengths for short successive internodes in the inferred species tree fall within the 45 theoretical range encompassing the anomaly zone. Distributions of empirical gene trees confirm 46 3 that the most common gene tree topology for each marker type differs from the species tree, 47 signifying the existence of an empirical anomaly zone in palaeognaths. 48 49

Powered flight was fundamental to the establishment and radiation of birds. However, flight has been lost multiple times throughout avian evolution. Convergent losses of flight within the ratites (flightless paleognaths, including the emu and ostrich) often coincide with reduced wings. Although there is a wealth of anatomical knowledge for several ratites, the genetic mechanisms causing these changes remain debated. Here, we use a multidisciplinary approach employing embryological, genetic, and genomic techniques to interrogate the mechanisms underlying forelimb heterochrony in emu embryos. We show that the initiation of limb formation, an epithelial to mesenchymal transition (EMT) in the lateral plate mesoderm (LPM) and myoblast migration into the LPM, occur at equivalent stages in the emu and chick. However, the emu forelimb fails to subsequently proliferate. The unique emu forelimb expression of Nkx2.5, previously associated with diminished wing development, initiates after this stage (concomitant with myoblast migration into the LPM) and is therefore unlikely to cause this developmental delay. In contrast, RNA sequencing of limb tissue reveals significantly lower Fgf10 expression in the emu forelimb. Artificially increasing Fgf10 expression in the emu LPM induces ectodermal Fgf8 expression and a limb bud. Analyzing open chromatin reveals differentially active regulatory elements near Fgf10 and Sall-1 in the emu wing, and the Sall-1 enhancer activity is dependent on a likely Fgf-mediated Ets transcription factor-binding site. Taken together, our results suggest that regulatory changes result in lower expression of Fgf10 and a concomitant failure to express genes required for limb proliferation in the early emu wing bud.

The relative roles of regulatory and protein evolution in the origin and loss of convergent phenotypic traits is a core question in evolutionary biology. Here we combine phylogenomic, epigenomic and developmental data to show that convergent evolution of regulatory regions, but not protein-coding genes, is associated with flightlessness in palaeognathous birds, a classic example of a convergent phenotype. Eleven new genomes, including a draft genome from an extinct moa, resolve palaeognath phylogeny and show that the incidence of independent, convergent accelerations among 284,000 conserved non-exonic elements is significantly more frequent in ratites than other bird lineages. Ratite-specific acceleration of conserved regions and measures of open chromatin across eight tissues in the developing chick identify candidate regulatory regions that may have modified or lost function in ratites. Enhancer activity assays conducted in the early developing chicken forelimb confirm that volant versions of a conserved element in the first intron of the TEAD1 gene display conserved enhancer activity, whereas an accelerated flightless version fails to drive reporter gene expression. Our results show that convergent molecular changes associated with loss of flight are largely regulatory in nature.