Investigating Difficult Nodes in the Placental Mammal Tree with Expanded Taxon Sampling and Thousands of Ultraconserved Elements

Phylogenomics of Afrotherian mammals and improved resolution of extant Paenungulata

Molecular investigations have gathered a diverse set of mammals-predominantly African natives like elephants, hyraxes, and aardvarks-into a clade known as Afrotheria. Nevertheless, the precise phylogenetic relationships among these species remain contentious. Here, we sourced orthologous markers and ultraconserved elements to discern the interordinal connections among Afrotherian mammals. Our phylogenetic analyses bolster the common origin of Afroinsectiphilia and Paenungulata, and propose Afrosoricida as the closer relative to Macroscelidea rather than Tubulidentata, while also challenging the notion of Sirenia and Hyracoidea as sister taxa. The approximately unbiased test and the gene concordance factor uniformly recognized the alliance of Proboscidea with Hyracoidea as the dominant topology within Paenungulata. Investigation into sites with extremly high phylogenetic signal unveiled their potential to intensify conflicts in the Paenungulata topology. Subsequent exploration suggested that incomplete lineage sorting was predominantly responsible for the observed contentious relationships, whereas introgression exerted a subsidiary influence. The divergence times estimated in our study hint at the Cretaceous-Paleogene (K-Pg) extinction event as a catalyst for Afrotherian diversification. Overall, our findings deliver a tentative but insightful overview of Afrotheria phylogeny and divergence, elucidating these relationships through the lens of phylogenomics.

Evolution and conservation genetics of pangolins

Pangolins (Pholidota, Manidae) are classified as an evolutionarily distinct and globally endangered mammal due to their unique morphology (nail-like scales and a myrmecophagous diet) and being the victim of heavy poaching and worldwide trafficking. As such, pangolins serve as a textbook example for studying the special phenotypic evolutionary adaptations and conservation genetics of an endangered species. Recent years have demonstrated significant advancements in the fields of molecular genetics and genomics, which have translated to a series of important research achievements and breakthroughs concerning the evolution and conservation genetics of pangolins. This review comprehensively presents the hitherto advances in phylogeny, adaptive evolution, conservation genetics, and conservation genomics that are related to pangolins, which will provide an ample understanding of their diversity, molecular adaptation mechanisms, and evolutionary potentials. In addition, we highlight the priority of investigating species/population diversity among pangolins and suggest several avenues of research that are highly relevant for future pangolin conservation.

Metazoa-level USCOs as markers in species delimitation and classification

Metazoa-level universal single-copy orthologs (mzl-USCOs) are universally applicable markers for DNA taxonomy in animals that can replace or supplement single-gene barcodes. Previously, mzl-USCOs from target enrichment data were shown to reliably distinguish species. Here, we tested whether USCOs are an evenly distributed, representative sample of a given metazoan genome and therefore able to cope with past hybridization events and incomplete lineage sorting. This is relevant for coalescent-based species delimitation approaches, which critically depend on the assumption that the investigated loci do not exhibit autocorrelation due to physical linkage. Based on 239 chromosome-level assembled genomes, we confirmed that mzl-USCOs are genetically unlinked for practical purposes and a representative sample of a genome in terms of reciprocal distances between USCOs on a chromosome and of distribution across chromosomes. We tested the suitability of mzl-USCOs extracted from genomes for species delimitation and phylogeny in four case studies: Anopheles mosquitos, Drosophila fruit flies, Heliconius butterflies and Darwin's finches. In almost all instances, USCOs allowed delineating species and yielded phylogenies that corresponded to those generated from whole genome data. Our phylogenetic analyses demonstrate that USCOs may complement single-gene DNA barcodes and provide more accurate taxonomic inferences. Combining USCOs from sources that used different versions of ortholog reference libraries to infer marker orthology may be challenging and, at times, impact taxonomic conclusions. However, we expect this problem to become less severe as the rapidly growing number of reference genomes provides a better representation of the number and diversity of organismal lineages.

Basicranial evidence suggests picrodontid mammals are not stem primates

The Picrodontidae from the middle Palaeocene of North America are enigmatic placental mammals that were allied with various mammalian groups but are generally now considered to have close affinities to paromomyid and palaechthonid plesiadapiforms based on proposed dental synapomorphies. The picrodontid fossil record consists entirely of dental and gnathic remains except for one partial cranium of Zanycteris paleocenus (AMNH 17180). Here, we use µCT technology to unveil previously undocumented morphology in AMNH 17180, describe and compare the basicranial morphology of a picrodontid for the first time, and incorporate these new data into cladistic analyses. The basicranial morphology of Z. paleocenus is distinct from plesiadapiforms and shares similarities with the Palaeogene Apatemyidae and Nyctitheriidae. Results of cladistic analyses incorporating these novel data suggest picrodontids are not stem primates nor euarchontan mammals and that the proposed dental synapomorphies uniting picrodontids with plesiadapiforms and, by extension, primates evolved independently. Results highlight the need to scrutinize proposed synapomorphies of highly autapomorphic taxa with limited fossil records.

Improved mammalian family phylogeny using gap-rare multiple sequence alignment: A timetree of extant placentals and marsupials

The timing of mammalian diversification in relation to the Cretaceous-Paleogene (KPg) mass extinction continues to be a subject of substantial debate. Previous studies have either focused on limited taxonomic samples with available whole-genome data or relied on short sequence alignments coupled with extensive species samples. In the present study, we improved an existing dataset from the landmark study of Meredith et al. (2011) by filling in missing fragments and further generated another dataset containing 120 taxa and 98 exonic markers. Using these two datasets, we then constructed phylogenies for extant mammalian families, providing improved resolution of many conflicting relationships. Moreover, the timetrees generated, which were calibrated using appropriate molecular clock models and multiple fossil records, indicated that the interordinal diversification of placental mammals initiated before the Late Cretaceous period. Additionally, intraordinal diversification of both extant placental and marsupial lineages accelerated after the KPg boundary, supporting the hypothesis that the availability of numerous vacant ecological niches subsequent to the mass extinction event facilitated rapid diversification. Thus, our results support a scenario of placental radiation characterized by both basal cladogenesis and active interordinal divergences spanning from the Late Cretaceous into the Paleogene.

Quantification of congruence among gene trees with polytomies using overall success of resolution for phylogenomic coalescent analyses

Gene-tree-inference error can cause species-tree-inference artefacts in summary phylogenomic coalescent analyses. Here we integrate two ways of accommodating these inference errors: collapsing arbitrarily or dubiously resolved gene-tree branches, and subsampling gene trees based on their pairwise congruence. We tested the effect of collapsing gene-tree branches with 0% approximate-likelihood-ratio-test (SH-like aLRT) support in likelihood analyses and strict consensus trees for parsimony, and then subsampled those partially resolved trees based on congruence measures that do not penalize polytomies. For this purpose we developed a new TNT script for congruence sorting (congsort), and used it to calculate topological incongruence for eight phylogenomic datasets using three distance measures: standard Robinson-Foulds (RF) distances; overall success of resolution (OSR), which is based on counting both matching and contradicting clades; and RF contradictions, which only counts contradictory clades. As expected, we found that gene-tree incongruence was often concentrated in clades that are arbitrarily or dubiously resolved and that there was greater congruence between the partially collapsed gene trees and the coalescent and concatenation topologies inferred from those genes. Coalescent branch lengths typically increased as the most incongruent gene trees were excluded, although branch supports typically did not. We investigated two successful and complementary approaches to prioritizing genes for investigation of alignment or homology errors. Coalescent-tree clades that contradicted concatenation-tree clades were generally less robust to gene-tree subsampling than congruent clades. Our preferred approach to collapsing likelihood gene-tree clades (0% SH-like aLRT support) and subsampling those trees (OSR) generally outperformed competing approaches for a large fungal dataset with respect to branch lengths, support and congruence. We recommend widespread application of this approach (and strict consensus trees for parsimony-based analyses) for improving quantification of gene-tree congruence/conflict, estimating coalescent branch lengths, testing robustness of coalescent analyses to gene-tree-estimation error, and improving topological robustness of summary coalescent analyses. This approach is quick and easy to implement, even for huge datasets.

The evolution of antimicrobial peptides in Chiroptera

High viral tolerance coupled with an extraordinary regulation of the immune response makes bats a great model to study host-pathogen evolution. Although many immune-related gene gains and losses have been previously reported in bats, important gene families such as antimicrobial peptides (AMPs) remain understudied. We built an exhaustive bioinformatic pipeline targeting the major gene families of defensins and cathelicidins to explore AMP diversity and analyze their evolution and distribution across six bat families. A combination of manual and automated procedures identified 29 AMP families across queried species, with α-, β-defensins, and cathelicidins representing around 10% of AMP diversity. Gene duplications were inferred in both α-defensins, which were absent in five species, and three β-defensin gene subfamilies, but cathelicidins did not show significant shifts in gene family size and were absent in Anoura caudifer and the pteropodids. Based on lineage-specific gains and losses, we propose diet and diet-related microbiome evolution may determine the evolution of α- and β-defensins gene families and subfamilies. These results highlight the importance of building species-specific libraries for genome annotation in non-model organisms and shed light on possible drivers responsible for the rapid evolution of AMPs. By focusing on these understudied defenses, we provide a robust framework for explaining bat responses to pathogens.

An atlas of rabbit development as a model for single-cell comparative genomics

Traditionally, the mouse has been the favoured vertebrate model for biomedical research, due to its experimental and genetic tractability. However, non-rodent embryological studies highlight that many aspects of early mouse development, such as its egg-cylinder gastrulation and method of implantation, diverge from other mammals, thus complicating inferences about human development. Like the human embryo, rabbits develop as a flat-bilaminar disc. Here we constructed a morphological and molecular atlas of rabbit development. We report transcriptional and chromatin accessibility profiles for over 180,000 single cells and high-resolution histology sections from embryos spanning gastrulation, implantation, amniogenesis and early organogenesis. Using a neighbourhood comparison pipeline, we compare the transcriptional landscape of rabbit and mouse at the scale of the entire organism. We characterize the gene regulatory programmes underlying trophoblast differentiation and identify signalling interactions involving the yolk sac mesothelium during haematopoiesis. We demonstrate how the combination of both rabbit and mouse atlases can be leveraged to extract new biological insights from sparse macaque and human data. The datasets and computational pipelines reported here set a framework for a broader cross-species approach to decipher early mammalian development, and are readily adaptable to deploy single-cell comparative genomics more broadly across biomedical research.

A global catalog of whole-genome diversity from 233 primate species

The rich diversity of morphology and behavior displayed across primate species provides an informative context in which to study the impact of genomic diversity on fundamental biological processes. Analysis of that diversity provides insight into long-standing questions in evolutionary and conservation biology and is urgent given severe threats these species are facing. Here, we present high-coverage whole-genome data from 233 primate species representing 86% of genera and all 16 families. This dataset was used, together with fossil calibration, to create a nuclear DNA phylogeny and to reassess evolutionary divergence times among primate clades. We found within-species genetic diversity across families and geographic regions to be associated with climate and sociality, but not with extinction risk. Furthermore, mutation rates differ across species, potentially influenced by effective population sizes. Lastly, we identified extensive recurrence of missense mutations previously thought to be human specific. This study will open a wide range of research avenues for future primate genomic research.

Reduction of Paraoxonase Expression Followed by Inactivation across Independent Semiaquatic Mammals Suggests Stepwise Path to Pseudogenization

Convergent adaptation to the same environment by multiple lineages frequently involves rapid evolutionary change at the same genes, implicating these genes as important for environmental adaptation. Such adaptive molecular changes may yield either change or loss of protein function; loss of function can eliminate newly deleterious proteins or reduce energy necessary for protein production. We previously found a striking case of recurrent pseudogenization of the Paraoxonase 1 (Pon1) gene among aquatic mammal lineages-Pon1 became a pseudogene with genetic lesions, such as stop codons and frameshifts, at least four times independently in aquatic and semiaquatic mammals. Here, we assess the landscape and pace of pseudogenization by studying Pon1 sequences, expression levels, and enzymatic activity across four aquatic and semiaquatic mammal lineages: pinnipeds, cetaceans, otters, and beavers. We observe in beavers and pinnipeds an unexpected reduction in expression of Pon3, a paralog with similar expression patterns but different substrate preferences. Ultimately, in all lineages with aquatic/semiaquatic members, we find that preceding any coding-level pseudogenization events in Pon1, there is a drastic decrease in expression, followed by relaxed selection, thus allowing accumulation of disrupting mutations. The recurrent loss of Pon1 function in aquatic/semiaquatic lineages is consistent with a benefit to Pon1 functional loss in aquatic environments. Accordingly, we examine diving and dietary traits across pinniped species as potential driving forces of Pon1 functional loss. We find that loss is best associated with diving activity and likely results from changes in selective pressures associated with hypoxia and hypoxia-induced inflammation.

Evolutionary constraint and innovation across hundreds of placental mammals

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

A genomic timescale for placental mammal evolution

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 signals 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 time trees reveal an accumulation of cladogenic events before and immediately after the Cretaceous-Paleogene (K-Pg) boundary, implying important roles for Cretaceous continental vicariance and the K-Pg extinction in the placental radiation.

Ultraconserved elements improve resolution of marmot phylogeny and offer insights into biogeographic history

Marmots (Marmota spp.) comprise a lineage of large-bodied ground squirrels that diversified rapidly in the Pleistocene, when the planet quickly transitioned to a drier, colder, and highly seasonal climate-particularly at high latitudes. Fossil evidence indicates the genus spread from North America, across Beringia, and into the European Alps over the course of only a few million years, beginning in the late Pliocene. Marmots are highly adapted to survive long and severely cold winters, and this likely favored their expansion and diversification over this time period. Previous phylogenetic studies have identified two major subgenera of marmots, but the timing of important speciation events and some species relationships have been difficult to resolve. Here we use ultraconserved elements and mitogenomes, with samples from all 15 extant species, to more precisely retrace how and when marmots came to inhabit a vast Holarctic range. Our results indicate marmots arose in North America in the mid Miocene (∼16.3 Mya) and dispersed across the Bering Land Bridge in the late Pliocene (∼3-4 Mya); in addition, our fossil-calibrated timeline is suggestive of the rise and spread of open grasslands as being particularly important to marmot diversification. The woodchuck (M. monax) and the Alaska marmot (M. broweri) are found to be more closely related to the Eurasian species than to the other North American species. Paraphyly is evident in the bobak marmot (M. bobak) and the hoary marmot (M. caligata), and in the case of the latter the data are highly suggestive of a second, cryptic species in the Cascade Mountains of Washington.

Molecular data from the holotype of the enigmatic Bornean Black Shrew, Suncusater Medway, 1965 (Soricidae, Crocidurinae), place it in the genus Palawanosorex

Although Borneo has received more attention from biologists than most other islands in the Malay Archipelago, many questions regarding the systematic relationships of Bornean mammals remain. Using next-generation sequencing technology, we obtained mitochondrial DNA sequences from the holotype of Suncusater, the only known specimen of this shrew. Several shrews collected recently in Sarawak are closely aligned, both morphologically and mitochondrially, with the holotype of S.ater. Phylogenetic analyses of mitochondrial sequences indicate that the S.ater holotype and new Sarawak specimens do not belong to the genus Suncus, but instead are most closely related to Palawanosorexmuscorum. Until now Palawanosorex has been known only from the neighboring Philippine island of Palawan. Additional sequences from nuclear ultra-conserved elements from the new Sarawak specimens strongly support a sister relationship to P.muscorum. We therefore transfer ater to Palawanosorex. The new specimens demonstrate that P.ater is more widespread in northern Borneo than previously recorded. Continued sampling of Bornean mammal diversity and reexamination of type material are critical in understanding the evolutionary history of the biologically rich Malay Archipelago.

Evolutionary analyses of polymeric immunoglobulin receptor (pIgR) in the mammals reveals an outstanding mutation rate in the lagomorphs

BACKGROUND

The transcytosis of polymeric immunoglobulins, IgA and IgM, across the epithelial barrier to the luminal side of mucosal tissues is mediated by the polymeric immunoglobulin receptor (pIgR). At the luminal side the extracellular ligand binding region of pIgR, the secretory component (SC), is cleaved and released bound to dimeric IgA (dIgA), protecting it from proteolytic degradation, or in free form, protecting the mucosa form pathogens attacks. The pIgR was first cloned for rabbit in early 1980's and since then has been described for all vertebrates, from fish to mammals. The existence of more than one functional pIgR alternative-spliced variant in the European rabbit, the complete pIgR as other mammals and a shorter pIgR lacking two SC exons, raised the question whether other lagomorphs share the same characteristics and how has the PIGR gene evolved in these mammals.

RESULTS

To investigate these questions, we sequenced expressed pIgR genes for other leporid genus, Lepus spp., and obtained and aligned pIgR sequences from representative species of all mammalian orders. The obtained mammalian phylogeny, as well as the Bayesian inference of evolutionary rates and genetic distances, show that Lagomorpha pIgR is evolving at a higher substitution rate. Codon-based analyses of positive selection show that mammalian pIgR is evolving under strong positive selection, with strong incidence in the domains excised from the rabbit short pIgR isoform. We further confirmed that the hares also express the two rabbit pIgR isoforms.

CONCLUSIONS

The Lagomorpha pIgR unique evolutionary pattern may reflect a group specific adaptation. The pIgR evolution may be linked to the unusual expansion of IgA genes observed in lagomorphs, or to neofunctionalization in this group. Further studies are necessary to clarify the driving forces behind the unique lagomorph pIgR evolution.

Hearing in African pygmy hedgehogs (Atelerix albiventris): audiogram, sound localization, and ear anatomy

The behavioral audiogram and sound localization performance, together with the middle and inner ear anatomy, were examined in African pygmy hedgehogs Atelerix albiventris. Their auditory sensitivity at 60 dB SPL extended from 2 to 46 kHz, revealing a relatively narrow hearing range of 4.6 octaves, with a best sensitivity of 21 dB at 8 kHz. Their noise-localization acuity around the midline (minimum audible angle) was 14°, matching the mean of terrestrial mammals. The African pygmy hedgehog was not able to localize low-frequency pure tones or a 3-kHz amplitude-modulated tone when forced to rely on the interaural phase-difference cue, a trait shared by at least nine other mammals. The middle ear of Atelerix has a configuration including an ectotympanic which is not fused to the surrounding bones, a substantial pars flaccida, a synostosed malleo-ectotympanic articulation and a 'microtype' malleus. The hearing and sound localization of A. albiventris is compared to that of a broad range of other mammals. It is shown that a malleus morphology like that of Atelerix, including a stiff articulation with the ectotympanic, is a consistent feature of other mammals that do not hear frequencies below 400 Hz.

Taxon-specific ultraconserved element probe design for phylogenetic analyses of scale insects (Hemiptera: Sternorrhyncha: Coccoidea)

Scale insects (Coccoidea) are morphologically specialized members of the order Hemiptera, with 56 families recognized to date. However, the phylogenetic relationships within and among families are poorly resolved. In this study, to further characterize the phylogenetic relationships among scale insects, an ultraconserved element (UCE) probe set was designed specifically for Coccoidea based on three low-coverage whole genome sequences along with three publicly available genomes. An in silico test including eight additional genomes was performed to evaluate the effectiveness of the probe set. Most scale insect lineages were recovered by the phylogenetic analysis. This study recovered the monophyly of neococcoids. The newly developed UCE probe set has the potential to reshape and improve our understanding of the phylogenetic relationships within and among families of scale insects at the genome level.

Total evidence time-scaled phylogenetic and biogeographic models for the evolution of sea cows (Sirenia, Afrotheria)

Molecular phylogenetic studies that have included sirenians from the genera Trichechus, Dugong, and Hydrodamalis have resolved their interrelationships but have yielded divergence age estimates that are problematically discordant. The ages of these lineage splits have profound implications for how to interpret the sirenian fossil record-including clade membership, biogeographic patterns, and correlations with Earth history events. In an effort to address these issues, here we present a total evidence phylogenetic analysis of Sirenia that includes living and fossil species and applies Bayesian tip-dating methods to estimate their interrelationships and divergence times. In addition to extant sirenians, our dataset includes 56 fossil species from 106 dated localities and numerous afrotherian outgroup taxa. Genetic, morphological, temporal, and biogeographic data are assessed simultaneously to bring all available evidence to bear on sirenian phylogeny. The resulting time-tree is then used for Bayesian geocoordinates reconstruction analysis, which models ancestral geographic areas at splits throughout the phylogeny, thereby allowing us to infer the direction and timing of dispersals. Our results suggest that Pan-Sirenia arose in North Africa during the latest Paleocene and that the Eocene evolution of stem sirenians was primarily situated in the Tethyan realm. In the late Eocene, some lineages moved into more northern European latitudes, an area that became the source region for a key trans-Atlantic dispersal towards the Caribbean and northern-adjacent west Atlantic. This event led to the phylogenetic and biogeographic founding of crown Sirenia with the Dugongidae-Trichechidae split occurring at the Eocene-Oligocene boundary (~33.9 Ma), temporally coincident with the onset of dropping global sea levels and temperatures. This region became the nexus of sirenian diversification and supported taxonomically-rich dugongid communities until the earliest Pliocene. The Dugonginae-Hydrodamalinae split occurred near Florida during the early Miocene (~21.2 Ma) and was followed by a west-bound dispersal that gave rise to the Pacific hydrodamalines. The late middle Miocene (~12.2 Ma) split of Dugong from all other dugongines also occurred near Florida and our analyses suggest that the Indo-Pacific distribution of modern dugongs is the result of a trans-Pacific dispersal. From at least the early Miocene, trichechid evolution was based entirely in South America, presumably within the Pebas Wetlands System. We infer that the eventual establishment of Amazon drainage into the South Atlantic allowed the dispersal of Trichechus out of South America no earlier than the mid-Pliocene. Our analyses provide a new temporal and biogeographic framework for understanding major events in sirenian evolution and their possible relationships to oceanographic and climatic changes. These hypotheses can be further tested with the recovery and integration of new fossil evidence.

Predicting Species Boundaries and Assessing Undescribed Diversity in Pneumocystis, an Obligate Lung Symbiont

Far more biodiversity exists in Fungi than has been described, or could be described in several lifetimes, given current rates of species discovery. Although this problem is widespread taxonomically, our knowledge of animal-associated fungi is especially lacking. Fungi in the genus Pneumocystis are obligate inhabitants of mammal lungs, and they have been detected in a phylogenetically diverse array of species representing many major mammal lineages. The hypothesis that Pneumocystis cospeciate with their mammalian hosts suggests that thousands of Pneumocystis species may exist, potentially equal to the number of mammal species. However, only six species have been described, and the true correspondence of Pneumocystis diversity to host species boundaries is unclear. Here, we use molecular species delimitation to estimate the boundaries of Pneumocystis species sampled from 55 mammal species representing eight orders. Our results suggest that Pneumocystis species often colonize several closely related mammals, especially those in the same genus. Using the newly estimated ratio of fungal to host diversity, we estimate ≈4600 to 6250 Pneumocystis species inhabit the 6495 currently recognized extant mammal species. Additionally, we review the literature and find that only 240 (~3.7%) mammal species have been screened for Pneumocystis, and many detected Pneumocystis lineages are not represented by any genetic data. Although crude, our findings challenge the dominant perspective of strict specificity of Pneumocystis to their mammal hosts and highlight an abundance of undescribed diversity.

Systematics of the New World bats Eptesicus and Histiotus suggest trans-marine dispersal followed by Neotropical cryptic diversification

Biodiversity can be boosted by colonization of new habitats such as remote islands and separated continents. Molecular studies have suggested that recently evolved organisms probably colonized already separated continents by dispersal, either via land bridge connections or crossing the ocean. Here we test the on-land and trans-marine dispersal hypotheses by evaluating possibilities of colonization routes over the Bering land bridge and across the Atlantic Ocean in the cosmopolitan bat genus Eptesicus (Chiroptera, Vespertilionidae). Previous molecular studies have found New World Eptesicus more closely related to Histiotus, a Neotropical endemic lineage with enlarged ears, than to Old World Eptesicus. However, phylogenetic relationships within the New World group remained unresolved and their evolutionary history was unclear. Here we studied the systematics of New World Eptesicus and Histiotus using extensive taxonomic and geographic sampling, and genomic data from thousands of ultra-conserved elements (UCEs). We estimated phylogenetic trees using concatenation and multispecies coalescent. All analyses supported four major New World clades and a novel topology where E. fuscus and Histiotus are sister clades that together diverged from two sister clades of Neotropical Eptesicus. Intra-clade divergence suggested cryptic diversity that has been concealed by morphological features, especially in the Neotropics where taxonomic re-evaluations are warranted. Molecular dating estimated that Old World and New World clades diverged around 17 million years ago followed by radiation of major New World clades in the mid-Miocene, when climatic changes might have facilitated global dispersal and radiation events. Biogeographic ancestral reconstruction supported the Neotropical origin of the New World clades, suggesting a trans-Atlantic colonization route from North Africa to the northern Neotropics. We highlight that trans-marine dispersal may be more prevalent than currently acknowledged and may be an important first step to global biodiversification.

Old specimens for old branches: Assessing effects of sample age in resolving a rapid Neotropical radiation of squirrels

Ultraconserved Elements (UCEs) have been useful to resolve challenging phylogenies of non-model clades, unpuzzling long-conflicted relationships in key branches of the Tree of Life at both deep and shallow levels. UCEs are often reliably recovered from historical samples, unlocking a vast number of preserved natural history specimens for analysis. However, the extent to which sample age and preservation method impact UCE recovery as well as downstream inferences remains unclear. Furthermore, there is an ongoing debate on how to curate, filter, and properly analyze UCE data when locus recovery is uneven across sample age and quality. In the present study we address these questions with an empirical dataset composed of over 3800 UCE loci from 219 historical and modern samples of Sciuridae, a globally distributed and ecologically important family of rodents. We provide a genome-scale phylogeny of two squirrel subfamilies (Sciurillinae and Sciurinae: Sciurini) and investigate their placement within Sciuridae. For historical specimens, recovery of UCE loci and mean length per locus were inversely related to sample age; deeper sequencing improved the number of UCE loci recovered but not locus length. Most of our phylogenetic inferences-performed on six datasets with alternative data-filtering strategies, and using three distinct optimality criteria-resulted in distinct topologies. Datasets containing more loci (40% and 50% taxa representativeness matrices) yielded more concordant topologies and higher support values than strictly filtered datasets (60% matrices) particularly with IQ-Tree and SVDquartets, while filtering based on information content provided better topological resolution for inferences with the coalescent gene-tree based approach in ASTRAL-III. We resolved deep relationships in Sciuridae (including among the five currently recognized subfamilies) and relationships among the deepest branches of Sciurini, but conflicting relationships remain at both genus- and species-levels for the rapid Neotropical tree squirrel radiation. Our results suggest that phylogenomic consensus can be difficult and heavily influenced by the age of available samples and the filtering steps used to optimize dataset properties.

Phylogenomics and Diversification of the Schistosomatidae Based on Targeted Sequence Capture of Ultra-Conserved Elements

Schistosomatidae Stiles and Hassall 1898 is a medically significant family of digenetic trematodes (Trematoda: Digenea), members of which infect mammals or birds as definitive hosts and aquatic or amphibious gastropods as intermediate hosts. Currently, there are 17 named genera, for many of which evolutionary interrelationships remain unresolved. The lack of a resolved phylogeny has encumbered our understanding of schistosomatid evolution, specifically patterns of host-use and the role of host-switching in diversification. Here, we used targeted sequence capture of ultra-conserved elements (UCEs) from representatives of 13 of the 17 named genera and 11 undescribed lineages that are presumed to represent either novel genera or species to generate a phylogenomic dataset for the estimation of schistosomatid interrelationships. This study represents the largest phylogenetic effort within the Schistosomatidae in both the number of loci and breadth of taxon sampling. We present a near-comprehensive family-level phylogeny providing resolution to several clades of long-standing uncertainty within Schistosomatidae, including resolution for the placement of the North American mammalian schistosomes, implying a second separate capture of mammalian hosts. Additionally, we present evidence for the placement of Macrobilharzia at the base of the Schistosoma + Bivitellobilharzia radiation. Patterns of definitive and intermediate host use and a strong role for intermediate host-switching are discussed relative to schistosomatid diversification.

Contradictory Phylogenetic Signals in the Laurasiatheria Anomaly Zone

Relationships among laurasiatherian clades represent one of the most highly disputed topics in mammalian phylogeny. In this study, we attempt to disentangle laurasiatherian interordinal relationships using two independent genome-level approaches: (1) quantifying retrotransposon presence/absence patterns, and (2) comparisons of exon datasets at the levels of nucleotides and amino acids. The two approaches revealed contradictory phylogenetic signals, possibly due to a high level of ancestral incomplete lineage sorting. The positions of Eulipotyphla and Chiroptera as the first and second earliest divergences were consistent across the approaches. However, the phylogenetic relationships of Perissodactyla, Cetartiodactyla, and Ferae, were contradictory. While retrotransposon insertion analyses suggest a clade with Cetartiodactyla and Ferae, the exon dataset favoured Cetartiodactyla and Perissodactyla. Future analyses of hitherto unsampled laurasiatherian lineages and synergistic analyses of retrotransposon insertions, exon and conserved intron/intergenic sequences might unravel the conflicting patterns of relationships in this major mammalian clade.

Testing Phylogenetic Stability with Variable Taxon Sampling

Over the past three decades, computational capabilities have grown at such a rapid rate that they have given rise to many computationally heavy science fields such as phylogenomics. As increasingly more genomes are sequenced in the three domains of life, larger and more species-complete phylogenetic tree reconstructions are leading to a better understanding of the tree of life and the evolutionary histories in deep times. However, these large datasets pose unique challenges from a modeling and computational perspective: accurately describing the evolutionary process of thousands of species is still beyond the capability of current models, while the computational burden limits our ability to test multiple hypotheses. Thus, it is common practice to reduce the size of a dataset by selecting species to represent a clade (taxon sampling). Unfortunately, this process is subjective, and comparisons of large tree of life studies show that choice and number of species used in a dataset can alter the topology obtained. Thus, taxon sampling is, in itself, a process that needs to be fully investigated to determine its effect on phylogenetic stability. Here, we present the theory and practical application of an automated pipeline that can be easily implemented to explore the effect of taxon sampling on phylogenetic reconstructions. The application of this approach was recently discussed in a study of Terrabacteria and shows its power in investigating the accuracy of deep nodes of a phylogeny.

Genome-wide ultraconserved elements resolve phylogenetic relationships and biogeographic history among Neotropical leaf-nosed bats in the genus Anoura (Phyllostomidae)

AnouraGray, 1838 are Neotropical nectarivorous bats and the most speciose genus within the phyllostomid subfamily Glossophaginae. However, Anoura species limits remain debated, and phylogenetic relationships remain poorly known, because previous studies used limited Anoura taxon sampling or focused primarily on higher-level relationships. Here, we conduct the first phylogenomic study of Anoura by analyzing 2039 genome-wide ultraconserved elements (UCEs) sequenced for 42 individuals from 8 Anoura species/lineages plus two outgroups. Overall, our results based on UCEs resolved relationships in the genus and supported (1) the monophyly of small-bodied Anoura species (previously genus Lonchoglossa); (2) monotypic status of A. caudifer; and (3) nested positions of "A. carishina", A. caudifer aequatoris, and A. geoffroyi peruana specimens within A. latidens, A. caudifer and A. geoffroyi, respectively (suggesting that these taxa are not distinct species). Additionally, (4) phylogenetic networks allowing reticulate edges did not explain gene tree discordance better than the species tree (without introgression), indicating that a coalescent model accounting for discordance solely through incomplete lineage sorting fit our data well. Sensitivity analyses indicated that our species tree results were not adversely affected by varying taxon sampling across loci. Tree calibration and Bayesian coalescent analyses dated the onset of diversification within Anoura to around ∼ 6-9 million years ago in the Miocene, with extant species diverging mainly within the past ∼ 4 million years. We inferred a historical biogeographical scenario for Anoura of parapatric speciation fragmenting the range of a wide-ranging ancestral lineage centered in the Central to Northern Andes, along with Pliocene-Pleistocene dispersal or founder event speciation in Amazonia and the Brazilian Atlantic forest during the last ∼ 2.5 million years.

The complete mitochondrial genome of Tylonycteris fulvida (Chiroptera: Vespertilionidae) in South China

The complete mitochondrial genome of Tylonycteris fulvida (Peters, 1872) was obtained using high-throughput sequencing technology. The genome is a circular molecule of 16,621 bp length, containing 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a control region. A phylogenetic tree of 13 protein-coding genes was constructed using IQ-TREE. Our result suggests that T. fulvida cluster within Chiroptera and Fereuungulata. The complete mitochondrial genome sequence of T. fulvida will be helpful for future taxonomic and phylogenetic studies on Chiroptera.

Odorant-binding proteins in canine anal sac glands indicate an evolutionarily conserved role in mammalian chemical communication

Background

Chemical communication is an important aspect of the behavioural ecology of a wide range of mammals. In dogs and other carnivores, anal sac glands are thought to convey information to conspecifics by secreting a pallet of small volatile molecules produced by symbiotic bacteria. Because these glands are unique to carnivores, it is unclear how their secretions relate to those of other placental mammals that make use of different tissues and secretions for chemical communication. Here we analyse the anal sac glands of domestic dogs to verify the secretion of proteins and infer their evolutionary relationship to those involved in the chemical communication of non-carnivoran mammals.

Results

Proteomic analysis of anal sac gland secretions of 17 dogs revealed the consistently abundant presence of three related proteins. Homology searches against online databases indicate that these proteins are evolutionary related to ‘odorant binding proteins’ (OBPs) found in a wide range of mammalian secretions and known to contribute to chemical communication. Screening of the dog’s genome sequence show that the newly discovered OBPs are encoded by a single cluster of three genes in the pseudoautosomal region of the X-chromosome. Comparative genomic screening indicates that the same locus is shared by a wide range of placental mammals and that it originated at least before the radiation of extant placental orders. Phylogenetic analyses suggest a dynamic evolution of gene duplication and loss, resulting in large gene clusters in some placental taxa and recurrent loss of this locus in others. The homology of OBPs in canid anal sac glands and those found in other mammalian secretions implies that these proteins maintained a function in chemical communication throughout mammalian evolutionary history by multiple shifts in expression between secretory tissues involved in signal release and nasal mucosa involved in signal reception.

Conclusions

Our study elucidates a poorly understood part of the biology of a species that lives in close association with humans. In addition, it shows that the protein repertoire underlying chemical communication in mammals is more evolutionarily stable than the variation of involved glands and tissues would suggest.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01910-w.

The effects of fossil taxa, hypothetical predicted ancestors, and a molecular scaffold on pseudoextinction analyses of extant placental orders

Pseudoextinction analyses, which simulate extinction in extant taxa, use molecular phylogenetics to assess the accuracy of morphological phylogenetics. Previous pseudoextinction analyses have shown a failure of morphological phylogenetics to place some individual placental orders in the correct superordinal clade. Recent work suggests that the inclusion of hypothetical ancestors of extant placental clades, estimated by ancestral state reconstructions of morphological characters, may increase the accuracy of morphological phylogenetic analyses. However, these studies reconstructed direct hypothetical ancestors for each extant taxon based on a well-corroborated molecular phylogeny, which is not possible for extinct taxa that lack molecular data. It remains to be determined if pseudoextinct taxa, and by proxy extinct taxa, can be accurately placed when their immediate hypothetical ancestors are unknown. To investigate this, we employed molecular scaffolds with the largest available morphological data set for placental mammals. Each placental order was sequentially treated as pseudoextinct by exempting it from the molecular scaffold and recoding soft morphological characters as missing for all its constituent species. For each pseudoextinct data set, we omitted the pseudoextinct taxon and performed a parsimony ancestral state reconstruction to obtain hypothetical predicted ancestors. Each pseudoextinct order was then evaluated in seven parsimony analyses that employed combinations of fossil taxa, hypothetical predicted ancestors, and a molecular scaffold. In treatments that included fossils, hypothetical predicted ancestors, and a molecular scaffold, only 8 of 19 pseudoextinct placental orders (42%) retained the same interordinal placement as on the molecular scaffold. In treatments that included hypothetical predicted ancestors but not fossils or a scaffold, only four placental orders (21%) were recovered in positions that are congruent with the scaffold. These results indicate that hypothetical predicted ancestors do not increase the accuracy of pseudoextinct taxon placement when the immediate hypothetical ancestor of the taxon is unknown. Hypothetical predicted ancestors are not a panacea for morphological phylogenetics.

Phylogenetic Signal and Bias in Paleontology

An unprecedented amount of evidence now illuminates the phylogeny of living mammals and birds on the Tree of Life. We use this tree to measure phylogenetic value of data typically used in paleontology (bones and teeth) from six datasets derived from five published studies. We ask three interrelated questions: 1) Can these data adequately reconstruct known parts of the Tree of Life? 2) Is accuracy generally similar for studies using morphology, or do some morphological datasets perform better than others? 3) Does the loss of non-fossilizable data cause taxa to occur in misleadingly basal positions? Adding morphology to DNA datasets usually increases congruence of resulting topologies to the well corroborated tree, but this varies among morphological datasets. Extant taxa with a high proportion of missing morphological characters can greatly reduce phylogenetic resolution when analyzed together with fossils. Attempts to ameliorate this by deleting extant taxa missing morphology are prone to decreased accuracy due to long-branch artefacts. We find no evidence that fossilization causes extinct taxa to incorrectly appear at or near topologically basal branches. Morphology comprises the evidence held in common by living taxa and fossils, and phylogenetic analysis of fossils greatly benefits from inclusion of molecular and morphological data sampled for living taxa, whatever methods are used for phylogeny estimation.

The whale shark genome reveals patterns of vertebrate gene family evolution

Chondrichthyes (cartilaginous fishes) are fundamental for understanding vertebrate evolution, yet their genomes are understudied. We report long-read sequencing of the whale shark genome to generate the best gapless chondrichthyan genome assembly yet with higher contig contiguity than all other cartilaginous fish genomes, and studied vertebrate genomic evolution of ancestral gene families, immunity, and gigantism. We found a major increase in gene families at the origin of gnathostomes (jawed vertebrates) independent of their genome duplication. We studied vertebrate pathogen recognition receptors (PRRs), which are key in initiating innate immune defense, and found diverse patterns of gene family evolution, demonstrating that adaptive immunity in gnathostomes did not fully displace germline-encoded PRR innovation. We also discovered a new toll-like receptor (TLR29) and three NOD1 copies in the whale shark. We found chondrichthyan and giant vertebrate genomes had decreased substitution rates compared to other vertebrates, but gene family expansion rates varied among vertebrate giants, suggesting substitution and expansion rates of gene families are decoupled in vertebrate genomes. Finally, we found gene families that shifted in expansion rate in vertebrate giants were enriched for human cancer-related genes, consistent with gigantism requiring adaptations to suppress cancer.

Uneven Missing Data Skew Phylogenomic Relationships within the Lories and Lorikeets

The resolution of the Tree of Life has accelerated with advances in DNA sequencing technology. To achieve dense taxon sampling, it is often necessary to obtain DNA from historical museum specimens to supplement modern genetic samples. However, DNA from historical material is generally degraded, which presents various challenges. In this study, we evaluated how the coverage at variant sites and missing data among historical and modern samples impacts phylogenomic inference. We explored these patterns in the brush-tongued parrots (lories and lorikeets) of Australasia by sampling ultraconserved elements in 105 taxa. Trees estimated with low coverage characters had several clades where relationships appeared to be influenced by whether the sample came from historical or modern specimens, which were not observed when more stringent filtering was applied. To assess if the topologies were affected by missing data, we performed an outlier analysis of sites and loci, and a data reduction approach where we excluded sites based on data completeness. Depending on the outlier test, 0.15% of total sites or 38% of loci were driving the topological differences among trees, and at these sites, historical samples had 10.9× more missing data than modern ones. In contrast, 70% data completeness was necessary to avoid spurious relationships. Predictive modeling found that outlier analysis scores were correlated with parsimony informative sites in the clades whose topologies changed the most by filtering. After accounting for biased loci and understanding the stability of relationships, we inferred a more robust phylogenetic hypothesis for lories and lorikeets.

Lagomorpha as a Model Morphological System

Due to their global distribution, invasive history, and unique characteristics, European rabbits are recognizable almost anywhere on our planet. Although they are members of a much larger group of living and extinct mammals [Mammalia, Lagomorpha (rabbits, hares, and pikas)], the group is often characterized by several well-known genera (e.g., Oryctolagus, Sylvilagus, Lepus, and Ochotona). This representation does not capture the extraordinary diversity of behavior and form found throughout the order. Model organisms are commonly used as exemplars for biological research, but there are a limited number of model clades or lineages that have been used to study evolutionary morphology in a more explicitly comparative way. We present this review paper to show that lagomorphs are a strong system in which to study macro- and micro-scale patterns of morphological change within a clade that offers underappreciated levels of diversity. To this end, we offer a summary of the status of relevant aspects of lagomorph biology.

Myoglobin primary structure reveals multiple convergent transitions to semi-aquatic life in the world's smallest mammalian divers

The speciose mammalian order Eulipotyphla (moles, shrews, hedgehogs, solenodons) combines an unusual diversity of semi-aquatic, semi-fossorial, and fossorial forms that arose from terrestrial forbearers. However, our understanding of the ecomorphological pathways leading to these lifestyles has been confounded by a fragmentary fossil record, unresolved phylogenetic relationships, and potential morphological convergence, calling for novel approaches. The net surface charge of the oxygen-storing muscle protein myoglobin (ZMb), which can be readily determined from its primary structure, provides an objective target to address this question due to mechanistic linkages with myoglobin concentration. Here, we generate a comprehensive 71 species molecular phylogeny that resolves previously intractable intra-family relationships and then ancestrally reconstruct ZMb evolution to identify ancient lifestyle transitions based on protein sequence alone. Our phylogenetically informed analyses confidently resolve fossorial habits having evolved twice in talpid moles and reveal five independent secondary aquatic transitions in the order housing the world's smallest endothermic divers.

Ultraconserved Elements Improve the Resolution of Difficult Nodes within the Rapid Radiation of Neotropical Sigmodontine Rodents (Cricetidae: Sigmodontinae)

Sigmodontine rodents (Cricetidae, Sigmodontinae) represent the second largest muroid subfamily and the most species-rich group of New World mammals, encompassing above 410 living species and ca. 87 genera. Even with advances on the clarification of sigmodontine phylogenetic relationships that have been made recently, the phylogenetic relationships among the 12 main group of genera (i.e., tribes) remain poorly resolved, in particular among those forming the large clade Oryzomyalia. This pattern has been interpreted as consequence of a rapid radiation upon the group entrance into South America. Here, we attempted to resolve phylogenetic relationships within Sigmodontinae using target capture and high-throughput sequencing of ultraconserved elements (UCEs). We enriched and sequenced UCEs for 56 individuals and collected data from four already available genomes. Analyses of distinct data sets, based on the capture of 4,634 loci, resulted in a highly resolved phylogeny consistent across different methods. Coalescent species-tree based approaches, concatenated matrices, and Bayesian analyses recovered similar topologies that were congruent at the resolution of difficult nodes. We recovered good support for the intertribal relationships within Oryzomyalia; for instance, the tribe Oryzomyini appears as the sister taxa of the remaining oryzomyalid tribes. The estimates of divergence times agree with results of previous studies. We inferred the crown age of the sigmodontine rodents at the end of Middle Miocene, while the main lineages of Oryzomyalia appear to have radiated in a short interval during the Late Miocene. Thus, the collection of a genomic scale data set with a wide taxonomic sampling, provided resolution for the first time of the relationships among the main lineages of Sigmodontinae. We expect the phylogeny presented here will become the backbone for future systematic and evolutionary studies of the group.

Assessing topological congruence among concatenation-based phylogenomic approaches in empirical datasets

Assessing the effect of methodological decisions on the resulting hypotheses is critical in phylogenetics. Recent studies have focused on evaluating how model selection, orthology definition and confounding factors affect phylogenomic results. Here, we compare the results of three concatenated phylogenetic methods (Maximum Likelihood, ML; Bayesian Inference, BI; Maximum Parsimony, MP) in 157 empirical phylogenomic datasets. The resulting trees were very similar, with 96.7% of all nodes shared between BI and ML (90.6% for ML-MP and 89.1% for BI-MP). Differing nodes were predominantly those of lower support. The main conclusions of most of the studies agreed for the three phylogenetic methods and the discordance involved nodes considered as recalcitrant problems in systematics. The differences between methods were proportionally larger in datasets that analyze the relationships at higher taxonomic levels (particularly phyla and kingdoms), and independent of the number of characters included in the datasets. Note: a spanish version of this article is available in the Supplementary material (Supplementary material online).

Collapsing dubiously resolved gene-tree branches in phylogenomic coalescent analyses

In two-step coalescent analyses of phylogenomic data, gene-tree topologies are treated as fixed prior to species-tree inference. Although all gene-tree conflict is assumed to be caused by lineage sorting when applying these methods, in empirical datasets much of the conflict can be caused by estimation error. Weakly supported and even arbitrarily resolved clades are important sources of this estimation error for gene trees inferred from few informative characters relative to the number of sampled terminals, and the resulting extraneous conflict among gene trees can negatively impact species-tree inference. In this study, we quantified the relative severity of alternative methods for collapsing gene-tree branches for seven empirical datasets and quantified their effects on species-tree inference. The branch-collapsing methods that we employed were based on the strict consensus of optimal topologies, various bootstrap thresholds, and 0% approximate likelihood ratio test (SH-like aLRT) support. Up to 86% of internal gene-tree branches are dubiously or arbitrarily resolved in reanalyses of these published phylogenomic datasets, and collapsing these branches increased inferred species-tree coalescent branch lengths by up to 455%. For two datasets, the longer inferred branch lengths sometimes impacted inference of anomaly-zone conditions. Although branch-collapsing methods did not consistently affect the species-tree topology, they often increased branch support. The more severe and clearly justified gene-tree branch-collapsing methods, which we recommend be broadly applied for two-step coalescent analyses, are use of the strict consensus in parsimony analyses and the collapse clades with 0% SH-like aLRT support in likelihood analyses. Collapsing dubiously or arbitrarily resolved branches in gene trees sometimes improved congruence between coalescent-based results and concatenation trees. In such cases, we contend that the resolution provided by concatenation should be preferred and that incomplete lineage sorting is a poor explanation for the initial conflict between phylogenetic approaches.

How to Tackle Phylogenetic Discordance in Recent and Rapidly Radiating Groups? Developing a Workflow Using Loricaria (Asteraceae) as an Example

A major challenge in phylogenetics and -genomics is to resolve young rapidly radiating groups. The fast succession of species increases the probability of incomplete lineage sorting (ILS), and different topologies of the gene trees are expected, leading to gene tree discordance, i.e., not all gene trees represent the species tree. Phylogenetic discordance is common in phylogenomic datasets, and apart from ILS, additional sources include hybridization, whole-genome duplication, and methodological artifacts. Despite a high degree of gene tree discordance, species trees are often well supported and the sources of discordance are not further addressed in phylogenomic studies, which can eventually lead to incorrect phylogenetic hypotheses, especially in rapidly radiating groups. We chose the high-Andean Asteraceae genus Loricaria to shed light on the potential sources of phylogenetic discordance and generated a phylogenetic hypothesis. By accounting for paralogy during gene tree inference, we generated a species tree based on hundreds of nuclear loci, using Hyb-Seq, and a plastome phylogeny obtained from off-target reads during target enrichment. We observed a high degree of gene tree discordance, which we found implausible at first sight, because the genus did not show evidence of hybridization in previous studies. We used various phylogenomic analyses (trees and networks) as well as the D-statistics to test for ILS and hybridization, which we developed into a workflow on how to tackle phylogenetic discordance in recent radiations. We found strong evidence for ILS and hybridization within the genus Loricaria. Low genetic differentiation was evident between species located in different Andean cordilleras, which could be indicative of substantial introgression between populations, promoted during Pleistocene glaciations, when alpine habitats shifted creating opportunities for secondary contact and hybridization.

Diverse phylogenomic datasets uncover a concordant scenario of laurasiatherian interordinal relationships

Resolving the interordinal relationships in the mammalian superorder Laurasiatheria has been among the most intractable problems in higher-level mammalian systematics, with many conflicting hypotheses having been proposed. The present study collected three different sources of genome-scale data with comprehensive taxon sampling of laurasiatherian species, including two protein-coding datasets (4,186 protein-coding genes for an amino acid dataset comprising 2,761,247 amino acid residues and a nucleotide dataset comprising 5,516,340 nucleotides from 1st and 2nd codon positions), an intronic dataset (1,210 introns comprising 1,162,723 nucleotides) and an ultraconserved elements (UCEs) dataset (1,246 UCEs comprising 1,946,472 nucleotides) from 40 species representing all six laurasiatherian orders and 7 non-laurasiatherian outgroups. Remarkably, phylogenetic trees reconstructed with the four datasets using different tree-building methods (RAxML, FastTree, ASTRAL and MP-EST) all supported the relationship (Eulipotyphla, (Chiroptera, ((Carnivora, Pholidota), (Cetartiodactyla, Perissodactyla)))). We find a resolution of interordinal relationships of Laurasiatheria among all types of markers used in the present study, and the likelihood ratio tests for tree comparisons confirmed that the present tree topology is the optimal hypothesis compared to other examined hypotheses. Jackknifing subsampling analyses demonstrate that the results of laurasiatherian tree reconstruction varied with the number of loci and ordinal representatives used, which are likely the two main contributors to phylogenetic disagreements of Laurasiatheria seen in previous studies. Our study provides significant insight into laurasiatherian evolution, and moreover, an important methodological strategy and reference for resolving phylogenies of adaptive radiation, which have been a long-standing challenge in the field of phylogenetics.

Phylogenomic Discordance in the Eared Seals is best explained by Incomplete Lineage Sorting following Explosive Radiation in the Southern Hemisphere

The phylogeny and systematics of fur seals and sea lions (Otariidae) have long been studied with diverse data types, including an increasing amount of molecular data. However, only a few phylogenetic relationships have reached acceptance because of strong gene-tree species tree discordance. Divergence times estimates in the group also vary largely between studies. These uncertainties impeded the understanding of the biogeographical history of the group, such as when and how trans-equatorial dispersal and subsequent speciation events occurred. Here, we used high-coverage genome-wide sequencing for 14 of the 15 species of Otariidae to elucidate the phylogeny of the family and its bearing on the taxonomy and biogeographical history. Despite extreme topological discordance among gene trees, we found a fully supported species tree that agrees with the few well-accepted relationships and establishes monophyly of the genus Arctocephalus. Our data support a relatively recent trans-hemispheric dispersal at the base of a southern clade, which rapidly diversified into six major lineages between 3 and 2.5 Ma. Otaria diverged first, followed by Phocarctos and then four major lineages within Arctocephalus. However, we found Zalophus to be nonmonophyletic, with California (Zalophus californianus) and Steller sea lions (Eumetopias jubatus) grouping closer than the Galapagos sea lion (Zalophus wollebaeki) with evidence for introgression between the two genera. Overall, the high degree of genealogical discordance was best explained by incomplete lineage sorting resulting from quasi-simultaneous speciation within the southern clade with introgresssion playing a subordinate role in explaining the incongruence among and within prior phylogenetic studies of the family. [Hybridization; ILS; phylogenomics; Pleistocene; Pliocene; monophyly.].

Phylogenomics and the Genetic Architecture of the Placental Mammal Radiation

The genomes of placental mammals are being sequenced at an unprecedented rate. Alignments of hundreds, and one day thousands, of genomes spanning the rich living and extinct diversity of species offer unparalleled power to resolve phylogenetic controversies, identify genomic innovations of adaptation, and dissect the genetic architecture of reproductive isolation. We highlight outstanding questions about the earliest phases of placental mammal diversification and the promise of newer methods, as well as remaining challenges, toward using whole genome data to resolve placental mammal phylogeny. The next phase of mammalian comparative genomics will see the completion and application of finished-quality, gapless genome assemblies from many ordinal lineages and closely related species. Interspecific comparisons between the most hypervariable genomic loci will likely reveal large, but heretofore mostly underappreciated, effects on population divergence, morphological innovation, and the origin of new species.

Phylogenetic Reconstruction Based on Synteny Block and Gene Adjacencies

Gene order can be used as an informative character to reconstruct phylogenetic relationships between species independently from the local information present in gene/protein sequences. PhyChro is a reconstruction method based on chromosomal rearrangements, applicable to a wide range of eukaryotic genomes with different gene contents and levels of synteny conservation. For each synteny breakpoint issued from pairwise genome comparisons, the algorithm defines two disjoint sets of genomes, named partial splits, respectively, supporting the two block adjacencies defining the breakpoint. Considering all partial splits issued from all pairwise comparisons, a distance between two genomes is computed from the number of partial splits separating them. Tree reconstruction is achieved through a bottom-up approach by iteratively grouping sister genomes minimizing genome distances. PhyChro estimates branch lengths based on the number of synteny breakpoints and provides confidence scores for the branches. PhyChro performance is evaluated on two data sets of 13 vertebrates and 21 yeast genomes by using up to 130,000 and 179,000 breakpoints, respectively, a scale of genomic markers that has been out of reach until now. PhyChro reconstructs very accurate tree topologies even at known problematic branching positions. Its robustness has been benchmarked for different synteny block reconstruction methods. On simulated data PhyChro reconstructs phylogenies perfectly in almost all cases, and shows the highest accuracy compared with other existing tools. PhyChro is very fast, reconstructing the vertebrate and yeast phylogenies in <15 min.

Genomic Characterization and Curation of UCEs Improves Species Tree Reconstruction

Ultraconserved genomic elements (UCEs) are generally treated as independent loci in phylogenetic analyses. The identification pipeline for UCE probes does not require prior knowledge of genetic identity, only selecting loci that are highly conserved, single copy, without repeats, and of a particular length. Here, we characterized UCEs from 11 phylogenomic studies across the animal tree of life, from birds to marine invertebrates. We found that within vertebrate lineages, UCEs are mostly intronic and intergenic, while in invertebrates, the majority are in exons. We then curated four different sets of UCE markers by genomic category from five different studies including: birds, mammals, fish, Hymenoptera (ants, wasps, and bees), and Coleoptera (beetles). Of genes captured by UCEs, we find that many are represented by two or more UCEs, corresponding to nonoverlapping segments of a single gene. We considered these UCEs to be nonindependent, merged all UCEs that belonged to a particular gene, constructed gene and species trees, and then evaluated the subsequent effect of merging cogenic UCEs on gene and species tree reconstruction. Average bootstrap support for merged UCE gene trees was significantly improved across all data sets apparently driven by the increase in loci length. Additionally, we conducted simulations and found that gene trees generated from merged UCEs were more accurate than those generated by unmerged UCEs. As loci length improves gene tree accuracy, this modest degree of UCE characterization and curation impacts downstream analyses and demonstrates the advantages of incorporating basic genomic characterizations into phylogenomic analyses. [Anchored hybrid enrichment; ants; ASTRAL; bait capture; carangimorph; Coleoptera; conserved nonexonic elements; exon capture; gene tree; Hymenoptera; mammal; phylogenomic markers; songbird; species tree; ultraconserved elements; weevils.]

Comparative and functional morphology of chevron bones among mammals

Tail morphology and function vary considerably across mammals. While studies of the mammalian tail have paid increasing attention to the caudal vertebrae, the chevron bones, ventrally positioned elements that articulate with the caudal vertebrae of most species and that serve to protect blood vessels and provide attachment sites for tail flexor musculature, have largely been ignored. Here, morphological variation in chevron bones is documented systematically among mammals possessing different tail locomotor functions, including prehensility, terrestrial propulsion (use for pentapedal locomotion), and postural prop, during which chevron bones are presumably under different mechanical stresses or serve different mechanical roles. Several chevron bone morphotypes were identified along the tail, varying both within and between tail regions. While some morphotypes were present across many or all clades surveyed, other morphotypes were clade-specific. Chevron bone dorsoventral height was examined at key vertebral levels among closely related species with different tail locomotor functions to assess whether variation followed any functional patterns. Dorsoventral height of chevron bones differed between prehensile- and nonprehensile-tailed, prop-tailed and nonprop-tailed, and pentapedal and nonpentapedal mammals. However, small sample sizes precluded rigorous statistical analyses. Distinctions were also qualified among species (not grouped by tail locomotor function), and the utility of metrics for quantifying specific aspects of chevron bone anatomy is discussed. This study offers information about the functional morphology of mammalian tails and has implications for reconstructing tail function in the fossil record.

Developmental origins and homologies of the hyracoid dentition

Understanding the origins of morphological specializations in mammals is a key goal in evolutionary biology. It can be accomplished by studying dental homology, which is at the core of most evolutionary and developmental studies. Here, we focused on the evolution and development of the specialized dentition of hyraxes for which dental homologies have long been debated, and could have implications on early placental evolution. Specifically, we analysed dental mineralization sequences of the three living genera of hyraxes and 17 fossil species using X-ray computed microtomography. Our results point out the labile position of vestigial upper teeth on jaw bones in extant species, associated with the frequently unusual premolar shape of deciduous canines over 50 Ma of hyracoid evolution. We proposed two evolutionary and developmental hypotheses to explain these original hyracoid dental characteristics. (a) The presence of a vestigial teeth on the maxilla in front of a complex deciduous canine could be interpreted as extra-teeth reminiscent of early placental evolution or sirenians, an order phylogenetically close to hyracoids and showing five premolars. (b) These vestigial teeth could also correspond to third incisors with a position unusually shifted on the maxilla, which could be explained by the dual developmental origin of these most posterior incisors and their degenerated condition. This integrative study allows discussion on the current evolutionary and developmental paradigms associated with the mammalian dentition. It also highlights the importance of nonmodel species to understand dental homologies.

A phylogenomic rodent tree reveals the repeated evolution of masseter architectures

Understanding the number of times a trait has evolved is a necessary foundation for comprehending its potential relationships with selective regimes, developmental constraints and evolutionary diversification. Rodents make up over 40% of extant mammalian species, and their ecological and evolutionary success has been partially attributed to the increase in biting efficiency that resulted from a forward shift of one or two portions of the masseter muscle from the zygomatic arch onto the rostrum. This forward shift has occurred in three discrete ways, but the number of times it has occurred has never been explicitly quantified. We estimated an ultrametric phylogeny, the first to include all rodent families, using thousands of ultraconserved elements. We examined support for evolutionary relationships among the five rodent suborders and then incorporated relevant fossils, fitted models of character evolution, and used stochastic character mapping to determine that a portion of the masseter muscle has moved forward onto the rostrum at least seven times (with one reversal) during the approximately 70 Myr history of rodents. Combined, the repeated evolution of this key innovation, its increasing prevalence through time, and the species diversity of clades with this character underscores the adaptive value of improved biting efficiency and the relative ease with which some advantageous traits arise.