Whole-Genome Analyses Resolve the Phylogeny of Flightless Birds (Palaeognathae) in the Presence of an Empirical Anomaly Zone.

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.

A nuclear genome assembly of an extinct flightless bird, the little bush moa.

We present a draft genome of the little bush moa (Anomalopteryx didiformis)-one of approximately nine species of extinct flightless birds from Aotearoa, New Zealand-using ancient DNA recovered from a fossil bone from the South Island. We recover a complete mitochondrial genome at 249.9× depth of coverage and almost 900 megabases of a male moa nuclear genome at ~4 to 5× coverage, with sequence contiguity sufficient to identify more than 85% of avian universal single-copy orthologs. We describe a diverse landscape of transposable elements and satellite repeats, estimate a long-term effective population size of ~240,000, identify a diverse suite of olfactory receptor genes and an opsin repertoire with sensitivity in the ultraviolet range, show that the wingless moa phenotype is likely not attributable to gene loss or pseudogenization, and identify potential function-altering coding sequence variants in moa that could be synthesized for future functional assays. This genomic resource should support further studies of avian evolution and morphological divergence.

Positive selection in the adhesion domain of Mus sperm Adam genes through gene duplications and function-driven gene complex formations.

BACKGROUND: Sperm and testes-expressed Adam genes have been shown to undergo bouts of positive selection in mammals. Despite the pervasiveness of positive selection signals, it is unclear what has driven such selective bouts. The fact that only sperm surface Adam genes show signals of positive selection within their adhesion domain has led to speculation that selection might be driven by species-specific adaptations to fertilization or sperm competition. Alternatively, duplications and neofunctionalization of Adam sperm surface genes, particularly as it is now understood in rodents, might have contributed to an acceleration of evolutionary rates and possibly adaptive diversification. RESULTS: Here we sequenced and conducted tests of selection within the adhesion domain of sixteen known sperm-surface Adam genes among five species of the Mus genus. We find evidence of positive selection associated with all six Adam genes known to interact to form functional complexes on Mus sperm. A subset of these complex-forming sperm genes also displayed accelerated branch evolution with Adam5 evolving under positive selection. In contrast to our previous findings in primates, selective bouts within Mus sperm Adams showed no associations to proxies of sperm competition. Expanded phylogenetic analysis including sequence data from other placental mammals allowed us to uncover ancient and recent episodes of adaptive evolution. CONCLUSIONS: The prevailing signals of rapid divergence and positive selection detected within the adhesion domain of interacting sperm Adams is driven by duplications and potential neofunctionalizations that are in some cases ancient (Adams 2, 3 and 5) or more recent (Adams 1b, 4b and 6).

Improved cytosine base editors generated from TadA variants.

Cytosine base editors (CBEs) enable programmable genomic C·G-to-T·A transition mutations and typically comprise a modified CRISPR-Cas enzyme, a naturally occurring cytidine deaminase, and an inhibitor of uracil repair. Previous studies have shown that CBEs utilizing naturally occurring cytidine deaminases may cause unguided, genome-wide cytosine deamination. While improved CBEs that decrease stochastic genome-wide off-targets have subsequently been reported, these editors can suffer from suboptimal on-target performance. Here, we report the generation and characterization of CBEs that use engineered variants of TadA (CBE-T) that enable high on-target C·G to T·A across a sequence-diverse set of genomic loci, demonstrate robust activity in primary cells and cause no detectable elevation in genome-wide mutation. Additionally, we report cytosine and adenine base editors (CABEs) catalyzing both A-to-I and C-to-U editing (CABE-Ts). Together with ABEs, CBE-Ts and CABE-Ts enable the programmable installation of all transition mutations using laboratory-evolved TadA variants with improved properties relative to previously reported CBEs.

Pervasive male-biased expression throughout the germline-specific regions of the sea lamprey genome supports key roles in sex differentiation and spermatogenesis.

Sea lamprey undergo programmed genome rearrangement (PGR) in which ∼20% of the genome is jettisoned from somatic cells during embryogenesis. Although the role of PGR in embryonic development has been studied, the role of the germline-specific region (GSR) in gonad development is unknown. We analysed RNA-sequence data from 28 sea lamprey gonads sampled across life-history stages, generated a genome-guided de novo superTranscriptome with annotations, and identified germline-specific genes (GSGs). Overall, we identified 638 GSGs that are enriched for reproductive processes and exhibit 36x greater odds of being expressed in testes than ovaries. Next, while 55% of the GSGs have putative somatic paralogs, the somatic paralogs are not differentially expressed between sexes. Further, putative orthologs of some the male-biased GSGs have known functions in sex determination or differentiation in other vertebrates. We conclude that the GSR of sea lamprey plays an important role in testicular differentiation and potentially sex determination.

Developmental Biology: Hox Timing Determines Limb Placement.

Hox genes are known to determine vertebral identity along with being required for normal limb patterning. A new study now finds that differential expression timing of Hox genes in the lateral plate mesoderm determines limb placement as well.

Evolutionary Dynamics of Sex Chromosomes of Paleognathous Birds.

Standard models of sex chromosome evolution propose that recombination suppression leads to the degeneration of the heterogametic chromosome, as is seen for the Y chromosome in mammals and the W chromosome in most birds. Unlike other birds, paleognaths (ratites and tinamous) possess large nondegenerate regions on their sex chromosomes (PARs or pseudoautosomal regions). It remains unclear why these large PARs are retained over >100 Myr, and how this retention impacts the evolution of sex chromosomes within this system. To address this puzzle, we analyzed Z chromosome evolution and gene expression across 12 paleognaths, several of whose genomes have recently been sequenced. We confirm at the genomic level that most paleognaths retain large PARs. As in other birds, we find that all paleognaths have incomplete dosage compensation on the regions of the Z chromosome homologous to degenerated portions of the W (differentiated regions), but we find no evidence for enrichments of male-biased genes in PARs. We find limited evidence for increased evolutionary rates (faster-Z) either across the chromosome or in differentiated regions for most paleognaths with large PARs, but do recover signals of faster-Z evolution in tinamou species with mostly degenerated W chromosomes, similar to the pattern seen in neognaths. Unexpectedly, in some species, PAR-linked genes evolve faster on average than genes on autosomes, suggested by diverse genomic features to be due to reduced efficacy of selection in paleognath PARs. Our analysis shows that paleognath Z chromosomes are atypical at the genomic level, but the evolutionary forces maintaining largely homomorphic sex chromosomes in these species remain elusive.

Attenuated Fgf Signaling Underlies the Forelimb Heterochrony in the Emu Dromaius novaehollandiae.

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.

Convergent regulatory evolution and loss of flight in paleognathous birds.

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.

Integrating natural history collections and comparative genomics to study the genetic architecture of convergent evolution.

Evolutionary convergence has been long considered primary evidence of adaptation driven by natural selection and provides opportunities to explore evolutionary repeatability and predictability. In recent years, there has been increased interest in exploring the genetic mechanisms underlying convergent evolution, in part, owing to the advent of genomic techniques. However, the current 'genomics gold rush' in studies of convergence has overshadowed the reality that most trait classifications are quite broadly defined, resulting in incomplete or potentially biased interpretations of results. Genomic studies of convergence would be greatly improved by integrating deep 'vertical', natural history knowledge with 'horizontal' knowledge focusing on the breadth of taxonomic diversity. Natural history collections have and continue to be best positioned for increasing our comprehensive understanding of phenotypic diversity, with modern practices of digitization and databasing of morphological traits providing exciting improvements in our ability to evaluate the degree of morphological convergence. Combining more detailed phenotypic data with the well-established field of genomics will enable scientists to make progress on an important goal in biology: to understand the degree to which genetic or molecular convergence is associated with phenotypic convergence. Although the fields of comparative biology or comparative genomics alone can separately reveal important insights into convergent evolution, here we suggest that the synergistic and complementary roles of natural history collection-derived phenomic data and comparative genomics methods can be particularly powerful in together elucidating the genomic basis of convergent evolution among higher taxa. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Comparative Genomics as a Foundation for Evo-Devo Studies in Birds.

Developmental genomics is a rapidly growing field, and high-quality genomes are a useful foundation for comparative developmental studies. A high-quality genome forms an essential reference onto which the data from numerous assays and experiments, including ChIP-seq, ATAC-seq, and RNA-seq, can be mapped. A genome also streamlines and simplifies the development of primers used to amplify putative regulatory regions for enhancer screens, cDNA probes for in situ hybridization, microRNAs (miRNAs) or short hairpin RNAs (shRNA) for RNA interference (RNAi) knockdowns, mRNAs for misexpression studies, and even guide RNAs (gRNAs) for CRISPR knockouts. Finally, much can be gleaned from comparative genomics alone, including the identification of highly conserved putative regulatory regions. This chapter provides an overview of laboratory and bioinformatics protocols for DNA extraction, library preparation, library quantification, and genome assembly, from fresh or frozen tissue to a draft avian genome. Generating a high-quality draft genome can provide a developmental research group with excellent resources for their study organism, opening the doors to many additional assays and experiments.

Izumo1 and Juno: the evolutionary origins and coevolution of essential sperm-egg binding partners.

Reproductive proteins are among the most rapidly evolving classes of proteins. For a subset of these, rapid evolution is driven by positive Darwinian selection despite vital, well-conserved, reproductive functions. Izumo1 is the only essential sperm-egg fusion protein currently known on mammalian sperm, and its egg receptor (Juno; formerly Folr4) was recently discovered. Male knockout mice for Izumo1 and female knockout mice for Juno are both healthy but sterile. Here, both sperm-egg binding proteins are shown to be evolving under positive selection. Within mammals, coevolution of Izumo1 and Juno is also uncovered, suggesting that similar forces have shaped the evolutionary histories of these binding partners within Mammalia. Additionally, genomic analyses reveal an ancient origin for the Izumo gene family, initially reported as conserved exclusively in mammals. Newly identified Izumo1 orthologues could serve reproductive functions in birds, fish and reptiles. Surprisingly, these same analyses support Juno's presence in mammals alone, suggesting a recent mammalian-specific duplication and neofunctionalization of the ancestral folate receptor. Despite the indispensability of their reproductive interaction, and their apparent coevolution within Mammalia, this binding pair arose through strikingly different evolutionary forces.

Positive Selection and the Evolution of izumo Genes in Mammals.

Most genes linked to male reproductive function have been known to evolve rapidly among species and to show signatures of positive selection. Different male species-specific reproductive strategies have been proposed to underlie positive selection, such as sperm competitive advantage and control over females postmating physiology. However, an underexplored aspect potentially affecting male reproductive gene evolution in mammals is the effect of gene duplications. Here we analyze the molecular evolution of members of the izumo gene family in mammals, a family of four genes mostly expressed in the sperm with known and potential roles in sperm-egg fusion. We confirm a previously reported bout of selection for izumo1 and establish that the bout of selection is restricted to the diversification of species of the superorder Laurasiatheria. None of the izumo genes showed evidence of positive selection in Glires (Rodentia and Lagomorpha), and in the case of the non-testes-specific izumo4, rapid evolution was driven by relaxed selection. We detected evidence of positive selection for izumo3 among Primates. Interestingly, positively selected sites include several serine residues suggesting modifications in protein function and/or localization among Primates. Our results suggest that positive selection is driven by aspects related to species-specific adaptations to fertilization rather than sexual selection.