John Gatesy

Adaptation to ant and/or termite consumption (myrmecophagy) in mammals constitutes a textbook example of convergent evolution, being independently derived in several mammalian lineages. Myrmecophagous species are characterized by striking convergent morphological adaptations such as skull elongation, enlargement of salivary glands, and long claws to dig into ant and termite nests. These evolutionary modifications also include anatomical regression, such as dental simplification or loss, reduction of masticatory muscles, and possessing a reduced set of taste buds. To gain insights into the molecular changes underlying the regression of these morpho-anatomical traits, we investigated the functionality of candidate genes related to dentition, gustation, and mastication in nine convergent myrmecophagous mammalian lineages. We examined these genes in a comparative phylogenetic context, paired with molecular evolutionary analyses, to estimate the relative timing of loss of gene function over the evolutionary history of each convergent lineage. We found that gustatory reduction and pseudogenization of masticatory myosin were often associated with the regression of dental genes. Evidence of pseudogenization events linked to oral anatomy dates to as early as the Cretaceous/Paleogene boundary, and is an ongoing process including examples of incipient gene inactivations. Whereas we found evidence for gene inactivations across all three functional categories occurring during distinct temporal intervals, there was variation in the sets of genes lost and the relative timing of inactivation events. The combined evidence suggests that the convergent evolution of myrmecophagy has occurred as a protracted process with distinct phases of anatomical evolution, over timescales as long as 60 Myr.

Objectives: The new data provide an important genomic resource for the Critically Endangered Cuban crocodile (Crocodylus rhombifer). Cuban crocodiles are restricted to the Zapata Swamp in southern Matanzas Province, Cuba, and readily hybridize with the widespread American crocodile (Crocodylus acutus) in areas of sympatry. The reported de novo assembly will contribute to studies of crocodylian evolutionary history and provide a resource for informing Cuban crocodile conservation.Data description: The final 2.2 Gb draft genome forC. rhombiferconsists of 41,387 scaffolds (contigs: N50 = 104.67 Kb; scaffold: N50-518.55 Kb). Benchmarking Universal Single-Copy Orthologs (BUSCO) identified 92.3% of the 3,354 genes in the vertebrata_odb10 database. Approximately 42% of the genome (960Mbp) comprises repeat elements. We predicted 30,138 unique protein-coding sequences (17,737 unique genes) in the genome assembly. Functional annotation found the top Gene Ontology annotations for Biological Processes, Molecular Function, and Cellular Component were regulation, protein, and intracellular, respectively. This assembly will support future macroevolutionary, conservation, and molecular studies of the Cuban crocodile.

The precise pattern and timing of speciation events that gave rise to all living placental mammals remain controversial. We provide a comprehensive phylogenetic analysis of genetic variation across an alignment of 241 placental mammal genome assemblies, addressing prior concerns regarding limited genomic sampling across species. We compared neutral genome-wide phylogenomic signal using concatenation and coalescent-based approaches, interrogated phylogenetic variation across chromosomes and analyzed extensive catalogs of structural variants. Interordinal relationships exhibit relatively low rates of phylogenomic conflict across diverse datasets and analytical methods. Conversely, X-chromosome versus autosome conflicts characterize multiple independent clades that radiated during the Cenozoic. Genomic timetrees reveal an accumulation of cladogenic events before and immediately following the KPg boundary implying important roles for Cretaceous continental vicariance and the KPg extinction in the placental radiation. One-Sentence Summary A comprehensive whole genome phylogeny of extant placental mammals reveals timing and patterns of ordinal diversification.

Retroelement insertions (RIs) are low-homoplasy characters that are ideal data for addressing deep evolutionary radiations, where gene tree reconstruction errors can severely hinder phylogenetic inference with DNA and protein sequence data. Phylogenomic studies of Neoaves, a large clade of birds (>9000 species) that first diversified near the Cretaceous–Paleogene boundary, have yielded an array of robustly supported, contradictory relationships among deep lineages. Here, we reanalyzed a large RI matrix for birds using recently proposed quartet-based coalescent methods that enable inference of large species trees including branch lengths in coalescent units, clade-support, statistical tests for gene flow, and combined analysis with DNA-sequence-based gene trees. Genome-scale coalescent analyses revealed extremely short branches at the base of Neoaves, meager branch support, and limited congruence with previous work at the most challenging nodes. Despite widespread topological conflicts with DNA-sequence-based trees, combined analyses of RIs with thousands of gene trees show emergent support for multiple higher-level clades (Columbea, Passerea, Columbimorphae, Otidimorphae, Phaethoquornithes). RIs express asymmetrical support for deep relationships within the subclade Afroaves that hints at ancient gene flow involving the owl lineage (Strigiformes). Because DNA-sequence data are challenged by gene tree-reconstruction error, analysis of RIs represents one approach for improving gene tree-based methods when divergences are deep, internodes are short, terminal branches are long, and introgressive hybridization further confounds species–tree inference.

The transition in Mysticeti (Cetacea) from capture of individual prey using teeth to bulk filtering batches of small prey using baleen ranks among the most dramatic evolutionary transformations in mammalian history. We review phylogenetic work on the homology of mysticete feeding structures from anatomical, ontogenetic, and genomic perspectives. Six characters with key functional significance for filter-feeding behavior are mapped to cladograms based on 11 morphological datasets to reconstruct evolutionary change across the teeth-to-baleen transition. This comparative summary within a common parsimony framework reveals extensive conflicts among independent systematic efforts but also broad support for the newly named clade Kinetomenta (Aetiocetidae + Chaeomysticeti). Complementary anatomical studies using CTscans and ontogenetic series hint at commonalities between the developmental programs for teeth and baleen, lending further support for a ‘transitional chimaeric feeder’ scenario that best explains current knowledge on the transition to filter feeding. For some extant mysticetes, the ontogenetic sequence in fetal specimens recapitulates the inferred evolutionary transformation: from teeth, to teeth and baleen, to just baleen. Phylogenetic mapping of inactivating mutations reveals mutational decay of ‘dental genes’ related to enamel formation before the emergence of crown Mysticeti, while ‘baleen genes’ that were repurposed or newly derived during the evolutionary elaboration of baleen currently are poorly characterized. Review and meta-analysis of available data suggest that the teeth-to-baleen transition in Mysticeti ranks among the best characterized macroevolutionary shifts due to the diversity of data from the genome, the fossil record, comparative anatomy, and ontogeny that directly bears on this remarkable evolutionary transformation.

Summary coalescent methods offer an alternative to the concatenation (supermatrix) approach for inferring phylogenetic relationships from genome-scale datasets. Given huge datasets, broad congruence between contrasting phylogenomic paradigms is often obtained, but empirical studies commonly show some well supported conflicts between concatenation and coalescence results and also between species trees estimated from alternative coalescent methods. Partitioned support indices can help arbitrate these discrepancies by pinpointing outlier loci that are unjustifiably influential at conflicting nodes. Partitioned coalescence support (PCS) recently was developed for summary coalescent methods, such as ASTRAL and MP-EST, that use the summed fits of individual gene trees to estimate the species tree. However, PCS cannot be implemented when distance-based coalescent methods (e.g., STAR, NJst, ASTRID, STEAC) are applied. Here, this deficiency is addressed by automating computation of ‘partitioned coalescent branch length’ (PCBL), a novel index that uses iterative removal of individual gene trees to assess the impact of each gene on every clade in a distance-based coalescent tree. Reanalyses of five phylogenomic datasets show that PCBL for STAR and NJst trees helps quantify the overall stability/instability of clades and clarifies disagreements with results from optimality-based coalescent analyses. PCBL scores reveal severe ‘missing taxa’, ‘apical nesting’, ‘misrooting’, and ‘basal dragdown’ biases. Contrived examples demonstrate the gross overweighting of outlier gene trees that drives these biases. Because of interrelated biases revealed by PCBL scores, caution should be exercised when using STAR and NJst, in particular when many taxa are analyzed, missing data are non-randomly distributed, and widespread gene-tree reconstruction error is suspected. Similar biases in the optimality-based coalescent method MP-EST indicate that congruence among species trees estimated via STAR, NJst, and MP-EST should not be interpreted as independent corroboration for phylogenetic relationships. Such agreements among methods instead might be due to the common defects of all three summary coalescent methods.

The loss of teeth and evolution of baleen racks in Mysticeti was a profound transformation that permitted baleen whales to radiate and diversify into a previously underutilized ecological niche of bulk filter-feeding on zooplankton and other small prey. Ancestral state reconstructions suggest that teeth were lost in the common ancestor of crown Mysticeti. Genomic studies provide some support for this hypothesis and suggest that the genetic toolkit for enamel production was inactivated in the common ancestor of living baleen whales. However, molecular studies to date have not provided direct evidence for the complete loss of teeth, including their dentin component, on the stem mysticete branch. Given these results, several questions remain unanswered: (1) Were teeth lost in a single step or did enamel loss precede dentin loss? (2) Was enamel lost early or late on the stem mysticete branch? (3) If enamel and dentin/tooth loss were decoupled in the ancestry of baleen whales, did dentin loss occur on the stem mysticete branch or independently in different crown mysticete lineages? To address these outstanding questions, we compiled and analyzed complete protein-coding sequences for nine tooth-related genes from cetaceans with available genome data. Seven of these genes are associated with enamel formation ( ACP4, AMBN, AMELX, AMTN, ENAM, KLK4, MMP20 ) whereas two other genes are either dentin-specific ( DSPP ) or tooth-specific ( ODAPH ) but not enamel-specific. Molecular evolutionary analyses indicate that all seven enamel-specific genes have inactivating mutations that are scattered across branches of the mysticete tree. Three of the enamel genes ( ACP4, KLK4, MMP20 ) have inactivating mutations that are shared by all mysticetes. The two genes that are dentin-specific ( DSPP ) or tooth-specific ( ODAPH ) do not have any inactivating mutations that are shared by all mysticetes, but there are shared mutations in Balaenidae as well as in Plicogulae (Neobalaenidae + Balaenopteroidea). These shared mutations suggest that teeth were lost at most two times. Shared inactivating mutations and dN/dS analyses, in combination with cetacean divergence times, were used to estimate inactivation times of genes and by proxy enamel and tooth phenotypes. The results of these analyses are most compatible with a two-step model for the loss of teeth in the ancestry of living baleen whales: enamel was lost very early on the stem Mysticeti branch followed by the independent loss of dentin (and teeth) in the common ancestors of Balaenidae and Plicogulae, respectively. These results imply that some stem mysticetes, and even early crown mysticetes, may have had vestigial teeth comprised of dentin with no enamel. Our results also demonstrate that all odontocete species (in our study) with absent or degenerative enamel have inactivating mutations in one or more of their enamel genes.