Phylogenomic data analyses provide evidence that Xenarthra and Afrotheria are sister groups

The phylogenetic positions of the 4 clades, Euarchontoglires, Laurasiatheria, Afrotheria, and Xenarthra, have been major issues in the recent discussion of basal relationships among placental mammals. However, despite considerable efforts these relationships, crucial to the understanding of eutherian evolution and biogeography, have remained essentially unresolved. Euarchontoglires and Laurasiatheria are generally joined into a common clade (Boreoeutheria), whereas the position of Afrotheria and Xenarthra relative to Boreoeutheria has been equivocal in spite of the use of comprehensive amounts of nuclear encoded sequences or the application of genome-level characters such as retroposons. The probable reason for this uncertainty is that the divergences took place long time ago and within a narrow temporal window, leaving only short common branches. With the aim of further examining basal eutherian relationships, we have collected conserved protein-coding sequences from 11 placental mammals, a marsupial and a bird, whose nuclear genomes have been largely sequenced. The length of the alignment of homologous sequences representing each individual species is 2,168,859 nt. This number of sites, representing 2840 protein-coding genes, exceeds by a considerable margin that of any previous study. The phylogenetic analysis joined Xenarthra and Afrotheria on a common branch, Atlantogenata. This topology was found to fit the data significantly better than the alternative trees.     

Congruent mammalian trees from mitochondrial and nuclear genes using Bayesian methods

Analyses of mitochondrial and nuclear gene sequences have often produced different mammalian tree topologies, undermining confidence in the merit of molecular approaches with respect to "traditional" morphological classification. The recent sequencing of the complete mitochondrial genomes of two additional rodents (Spalax judaei and Jaculus jaculus) and one lagomorph (Ochotona princeps) has prompted us to reinvestigate the issue. Using Bayesian phylogenetics, we found phylogenetic relationships between mammalian species highly congruent with previous results based on nuclear genes. Our results show the existence of four primary lineages of placental mammals: Xenarthra, Afrotheria, Laurasiatheria, and Euarchontoglires. Relationships between and within these lineages strongly suggest that the gene trees may also be congruent with the underlying species phylogeny.     

Chromosomal instability in Afrotheria: fragile sites,evolutionary breakpoints and phylogenetic inference from genome sequence assemblies

Background  

Extant placental mammals are divided into four major clades (Laurasiatheria, Supraprimates, Xenarthra and Afrotheria). Given that Afrotheria is generally thought to root the eutherian tree in phylogenetic analysis of large nuclear gene data sets, the study of the organization of the genomes of afrotherian species provides new insights into the dynamics of mammalian chromosomal evolution. Here we test if there are chromosomal bands with a high tendency to break and reorganize in Afrotheria, and by analyzing the expression of aphidicolin-induced common fragile sites in three afrotherian species, whether these are coincidental with recognized evolutionary breakpoints.     

PATTERNS OF OSSIFICATION IN SOUTHERN VERSUS NORTHERN PLACENTAL MAMMALS

Consensus on placental mammal phylogeny is fairly recent compared to that for vertebrates as a whole. A stable phylogenetic hypothesis enables investigation into the possibility that placental clades differ from one another in terms of their development. Here, we focus on the sequence of skeletal ossification as a possible source of developmental distinctiveness in “northern” (Laurasiatheria and Euarchontoglires) versus “southern” (Afrotheria and Xenarthra) placental clades. We contribute data on cranial and postcranial ossification events during growth in Afrotheria, including elephants, hyraxes, golden moles, tenrecs, sengis, and aardvarks. We use three different techniques to quantify sequence heterochrony: continuous method, sequence‐ANOVA (analysis of variance) and event‐paring/Parsimov. We show that afrotherians significantly differ from other placentals by an early ossification of the orbitosphenoid and caudal vertebrae. Our analysis also suggests that both southern placental groups show a greater degree of developmental variability; however, they rarely seem to vary in the same direction, especially regarding the shifts that differ statistically. The latter observation is inconsistent with the Atlantogenata hypothesis in which afrotherians are considered as the sister clade of xenarthrans. Interestingly, ancestral nodes for Laurasiatheria and Euarchontoglires show very similar trends and our results suggest that developmental homogeneity in some ossification sequences may be restricted to northern placental mammals (Boreoeutheria).     

Recapitulating the evolution of Afrotheria: 57 genes and rare genomic changes (RGCs) consolidate their history

The decipherment of higher level relationships among the orders of Afrotheria – an extraordinary assumption in mammalian evolution – constitutes one of the major disputes in the evolutionary history of mammals. Recent comprehensive studies of various genomic data, including mitochondrial and nuclear DNA sequences, chromosomal syntenic associations and retroposon insertions support strongly the monophyly of Afrotheria. However, the relationships within Afrotheria have remained ambiguous and there is a necessity for a more sophisticated analysis (i.e. combination of gene phylogeny and Rare Genomic Changes (RGCs)), which could aid in the comprehension of the evolutionary history of this old group of mammals. The present study investigated the phylogenetic relationships within Afrotheria by analysing a data set of coding and non-coding sequences (～32 000 bp) comprising 57 orthologous genes and 31 RGCs, such as chromosomal associations and retroposon insertions, and re-evaluated a molecular timescale for afrotherian mammals using a Bayesian relaxed clock approach. The interordinal afrotherians phylogeny presented here contributed to the elucidation of the evolutionary history of this ancient clade of mammals, which is one of the most unorthodox proposals in mammalian biology. This is critical not only for understanding how Afrotheria evolved in Africa, but also to comprehend the early biogeographical history of placental mammals.     

Molecular phylogeny of living xenarthrans and the impact of character and taxon sampling on the placental tree rooting

Extant xenarthrans (armadillos, anteaters and sloths) are among the most derived placental mammals ever evolved. South America was the cradle of their evolutionary history. During the Tertiary, xenarthrans experienced an extraordinary radiation, whereas South America remained isolated from other continents. The 13 living genera are relics of this earlier diversification and represent one of the four major clades of placental mammals. Sequences of the three independent protein-coding nuclear markers alpha2B adrenergic receptor (ADRA2B), breast cancer susceptibility (BRCA1), and von Willebrand Factor (VWF) were determined for 12 of the 13 living xenarthran genera. Comparative evolutionary dynamics of these nuclear exons using a likelihood framework revealed contrasting patterns of molecular evolution. All codon positions of BRCA1 were shown to evolve in a strikingly similar manner, and third codon positions appeared less saturated within placentals than those of ADRA2B and VWF. Maximum likelihood and Bayesian phylogenetic analyses of a 47 placental taxa data set rooted by three marsupial outgroups resolved the phylogeny of Xenarthra with some evidence for two radiation events in armadillos and provided a strongly supported picture of placental interordinal relationships. This topology was fully compatible with recent studies, dividing placentals into the Southern Hemisphere clades Afrotheria and Xenarthra and a monophyletic Northern Hemisphere clade (Boreoeutheria) composed of Laurasiatheria and Euarchontoglires. Partitioned likelihood statistical tests of the position of the root, under different character partition schemes, identified three almost equally likely hypotheses for early placental divergences: a basal Afrotheria, an Afrotheria + Xenarthra clade, or a basal Xenarthra (Epitheria hypothesis). We took advantage of the extensive sampling realized within Xenarthra to assess its impact on the location of the root on the placental tree. By resampling taxa within Xenarthra, the conservative Shimodaira-Hasegawa likelihood-based test of alternative topologies was shown to be sensitive to both character and taxon sampling.     

A new phylogenetic marker,apolipoprotein B,provides compelling evidence for eutherian relationships

Higher-level relationships within, and the root of Placentalia, remain contentious issues. Resolution of the placental tree is important to the choice of mammalian genome projects and model organisms, as well as for understanding the biogeography of the eutherian radiation. We present phylogenetic analyses of 63 species representing all extant eutherian mammal orders for a new molecular phylogenetic marker, a 1.3kb portion of exon 26 of the apolipoprotein B (APOB) gene. In addition, we analyzed a multigene concatenation that included APOB sequences and a previously published data set (Murphy et al., 2001b) of three mitochondrial and 19 nuclear genes, resulting in an alignment of over 17kb for 42 placentals and two marsupials. Due to computational difficulties, previous maximum likelihood analyses of large, multigene concatenations for placental mammals have used quartet puzzling, less complex models of sequence evolution, or phylogenetic constraints to approximate a full maximum likelihood bootstrap. Here, we utilize a Unix load sharing facility to perform maximum likelihood bootstrap analyses for both the APOB and concatenated data sets with a GTR+Gamma+I model of sequence evolution, tree-bisection and reconnection branch-swapping, and no phylogenetic constraints. Maximum likelihood and Bayesian analyses of both data sets provide support for the superordinal clades Boreoeutheria, Euarchontoglires, Laurasiatheria, Xenarthra, Afrotheria, and Ostentoria (pangolins+carnivores), as well as for the monophyly of the orders Eulipotyphla, Primates, and Rodentia, all of which have recently been questioned. Both data sets recovered an association of Hippopotamidae and Cetacea within Cetartiodactyla, as well as hedgehog and shrew within Eulipotyphla. APOB showed strong support for an association of tarsier and Anthropoidea within Primates. Parsimony, maximum likelihood and Bayesian analyses with both data sets placed Afrotheria at the base of the placental radiation. Statistical tests that employed APOB to examine a priori hypotheses for the root of the placental tree rejected rooting on myomorphs and hedgehog, but did not discriminate between rooting at the base of Afrotheria, at the base of Xenarthra, or between Atlantogenata (Xenarthra+Afrotheria) and Boreoeutheria. An orthologous deletion of 363bp in the aligned APOB sequences proved phylogenetically informative for the grouping of the order Carnivora with the order Pholidota into the superordinal clade Ostentoria. A smaller deletion of 237-246bp was diagnostic of the superordinal clade Afrotheria.     

The chromosomes of afrotheria and their bearing on Mammalian genome evolution

Afrotheria is the clade of placental mammals that, together with Xenarthra, Euarchontoglires and Laurasiatheria, represents 1 of the 4 main recognized supraordinal eutherian clades. It reunites 6 orders of African origin: Proboscidea, Sirenia, Hyracoidea, Macroscelidea, Afrosoricida and Tubulidentata. The apparently unlikely relationship among such disparate morphological taxa and their possible basal position at the base of the eutherian phylogenetic tree led to a great deal of attention and research on the group. The use of biomolecular data was pivotal in Afrotheria studies, as they were the basis for the recognition of this clade. Although morphological evidence is still scarce, a plethora of molecular data firmly attests to the phylogenetic relationship among these mammals of African origin. Modern cytogenetic techniques also gave a significant contribution to the study of Afrotheria, revealing chromosome signatures for the group as a whole, as well as for some of its internal relationships. The associations of human chromosomes HSA1/19 and 5/21 were found to be chromosome signatures for the group and provided further support for Afrotheria. Additional chromosome synapomorphies were also identified linking elephants and manatees in Tethytheria (the associations HSA2/3, 3/13, 8/22, 18/19 and the lack of HSA4/8) and elephant shrews with the aardvark (HSA2/8, 3/20 and 10/17). Herein, we review the current knowledge on Afrotheria chromosomes and genome evolution. The already available data on the group suggests that further work on this apparently bizarre assemblage of mammals will provide important data to a better understanding on mammalian genome evolution.     

Timing of cranial suture closure in placental mammals: Phylogenetic patterns,intraspecific variation,and comparison with marsupials

Used as markers of postnatal growth closure sequences of 22 ectocranial sutures and synchondroses were recorded in a sample of 1161 skulls belonging to 38 species from all major placental clades: Afrotheria, Xenarthra, Laurasiatheria and Euarchontoglires (Boreoeutheria). The maximum closure level, which is not significantly correlated to body mass, is higher in Afrotheria and Xenarthra than in Boreoeutheria. Only the basioccipito‐basisphenoid and the basioccipito‐exoccipital synchondroses close in all species sampled, the supraoccipito‐exoccipital and the inter‐parietal sutures do in most species. Parsimov retrieved more heterochronic shifts for Afrotheria and Xenarthra than for Boreoeutheria. The amount of intraspecific variation differs among the species sampled being high among xenarthran species and low among afrotherians. Specimens (162) representing 12 marsupial genera were also analysed. Placentals exhibit a larger number of suture closures than marsupials and in both groups the sutures at the base of the skull are the first to fuse starting with the basioccipito‐exoccipital. J. Morphol. 275:125–140, 2014. © 2013 Wiley Periodicals, Inc.     

Genome size evolution: sizing Mammalian genomes

The study of genome size (GS) and its variation is so fascinating to the scientific community because it constitutes the link between the present-day analytical and molecular studies of the genome and the old trunk of the holistic and synthetic view of the genome. The GS of several taxa vary over a broad range and do not correlate with the complexity of the organisms (the C-value paradox). However, the biology of transposable elements has let us reach a satisfactory view of the molecular mechanisms that give rise to GS variation and novelties, providing a less perplexing view of the significance of the GS (C-enigma). The knowledge of the composition and structure of a genome is a pre-requisite for trying to understand the evolution of the main genome signature: its size. The radiation of mammals provides an approximately 180-million-year test case for theories of how GS evolves. It has been found from data-mining GS databases that GS is a useful cyto-taxonomical instrument at the level of orders/superorders, providing genomic signatures characterizing Monotremata, Marsupialia, Afrotheria, Xenarthra, Laurasiatheria, and Euarchontoglires. A hypothetical ancestral mammalian-like GS of 2.9-3.7 pg has been suggested. This value appears compatible with the average values calculated for the high systematic levels of the extant Monotremata (～2.97 pg) and Marsupialia (～4.07 pg), suggesting invasion of mobile DNA elements concurrently with the separation of the older clades of Afrotheria (～5.5 pg) and Xenarthra (～4.5 pg) with larger GS, leaving the Euarchontoglires (～3.4 pg) and Laurasiatheria (～2.8 pg) genomes with fewer transposable elements. However, the paucity of GS data (546 mammalian species sized from 5,488 living species) for species, genera, and families calls for caution. Considering that mammalian species may be vanished even before they are known, GS data are sorely needed to phenotype the effects brought about by their variation and to validate any hypotheses on GS evolution in mammals.     

Understanding mammalian evolution using Bayesian phylogenetic inference

1. Phylogenetic trees are critical in addressing evolutionary hypotheses; however, the reconstruction of a phylogeny is no easy task. This process has recently been made less arduous by using a Bayesian statistical approach. This method offers the advantage that one can determine the probability of some hypothesis (i.e. a phylogeny), conditional on the observed data (i.e. nucleotide sequences). 2. By reconstructing phylogenies using Bayes’ theorem in combination with Markov chain Monte Carlo, phylogeneticists are able to test hypotheses more quickly compared with using standard methods such as neighbour-joining, maximum likelihood or parsimony. Critics of the Bayesian approach suggest that it is not a panacea, and argue that the prior probability is too subjective and the resulting posterior probability is too liberal compared with maximum likelihood. 3. These issues are currently debated in the arena of mammalian evolution. Recently, proponents and opponents of the Bayesian approach have constructed the mammalian phylogeny using different methods under different conditions and with a variety of parameters. These analyses showed the robustness (or lack of) of the Bayesian approach. In the end, consensus suggests that Bayesian methods are robust and give essentially the same answer as maximum likelihood methods but in less time. 4. Approaches based on fossils and molecules typically agree on ordinal-level relationships among mammals but not on higher-level relationships, as Bayesian analyses recognize the African radiation, Afrotheria, and the two Laurasian radiations, Laurasiatheria and Euarchontoglires, whereas fossils did not predict Afrotheria.     

Indel evolution of mammalian introns and the utility of non-coding nuclear markers in eutherian phylogenetics

Nuclear DNA intron sequences are increasingly used to investigate evolutionary relationships among closely related organisms. The phylogenetic usefulness of intron sequences at higher taxonomic levels has, however, not been firmly established and very few studies have used these markers to address evolutionary questions above the family level. In addition, the mechanisms driving intron evolution are not well understood. We compared DNA sequence data derived from three presumably independently segregating introns (THY, PRKC I and MGF) across 158 mammalian species. All currently recognized extant eutherian mammalian orders were included with the exception of Cingulata, Dermoptera and Scandentia. The total aligned length of the data was 6366 base pairs (bp); after the exclusion of autapomorphic insertions, 1511 bp were analyzed. In many instances the Bayesian and parsimony analyses were complementary and gave significant posterior probability and bootstrap support (>80) for the monophyly of Afrotheria, Euarchontoglires, Laurasiatheria and Boreoeutheria. Apart from finding congruent support when using these methods, the intron data also provided several indels longer than 3 bp that support, among others, the monophyly of Afrotheria, Paenungulata, Ferae and Boreoeutheria. A quantitative analysis of insertions and deletions suggested that there was a 75% bias towards deletions. The average insertion size in the mammalian data set was 16.49 bp +/- 57.70 while the average deletion was much smaller (4.47 bp +/- 14.17). The tendency towards large insertions and small deletions is highlighted by the observation that out of a total of 17 indels larger than 100 bp, 15 were insertions. The majority of indels (>60% of all events) were 1 or 2 bp changes. Although the average overall indel substitution rate of 0.00559 per site is comparable to that previously reported for rodents and primates, individual analyses among different evolutionary lineages provide evidence for differences in the formation rate of indels among the different mammalian groups.     

Genome sizes in afrotheria, xenarthra, euarchontoglires, and laurasiatheria

Topical literature and Web site databases provide genome sizes for approximately 4,000 animal species, invertebrates and vertebrates, 330 of which are mammals. We provide the genome size for 67 mammalian species, including 51 never reported before. Knowledge of genome size facilitates sequencing projects. The data presented here encompassed 5 Metatheria (order Didelphimorphia) and 62 Eutheria: 15 Xenarthra, 24 Euarchontoglires (Rodentia), as well as 23 Laurasiatheria (22 Chiroptera and 1 species from Perissodactyla). Already available karyotypes supplement the haploid nuclear DNA contents of the respective species. Thus, we established the first comprehensive set of genome size measurements for 15 Xenarthra species (armadillos) and for 12 house-mouse species; each group was previously represented by only one species. The Xenarthra exhibited much larger genomes than the modal 3 pg DNA known for mammals. Within the genus Mus, genome sizes varied between 2.98 pg and 3.68 pg. The 22 bat species we measured support the low 2.63 pg modal value for Chiroptera. In general, the genomes of Euarchontoglires and Laurasiatheria were found being smaller than those of (Afrotheria and) Xenarthra. Interspecific variation in genome sizes is discussed with particular attention to repetitive elements, which probably promoted the adaptation of extant mammals to their environment.     

Confirming the phylogeny of mammals by use of large comparative sequence data sets

The ongoing generation of prodigious amounts of genomic sequence data from myriad vertebrates is providing unparalleled opportunities for establishing definitive phylogenetic relationships among species. The size and complexities of such comparative sequence data sets not only allow smaller and more difficult branches to be resolved but also present unique challenges, including large computational requirements and the negative consequences of systematic biases. To explore these issues and to clarify the phylogenetic relationships among mammals, we have analyzed a large data set of over 60 megabase pairs (Mb) of high-quality genomic sequence, which we generated from 41 mammals and 3 other vertebrates. All sequences are orthologous to a 1.9-Mb region of the human genome that encompasses the cystic fibrosis transmembrane conductance regulator gene (CFTR). To understand the characteristics and challenges associated with phylogenetic analyses of such a large data set, we partitioned the sequence data in several ways and utilized maximum likelihood, maximum parsimony, and Neighbor-Joining algorithms, implemented in parallel on Linux clusters. These studies yielded well-supported phylogenetic trees, largely confirming other recent molecular phylogenetic analyses. Our results provide support for rooting the placental mammal tree between Atlantogenata (Xenarthra and Afrotheria) and Boreoeutheria (Euarchontoglires and Laurasiatheria), illustrate the difficulty in resolving some branches even with large amounts of data (e.g., in the case of Laurasiatheria), and demonstrate the valuable role that very large comparative sequence data sets can play in refining our understanding of the evolutionary relationships of vertebrates.     

Genome Size: A Novel Genomic Signature in Support of Afrotheria

Molecular phylogenetic analyses suggest an emerging phylogeny for the extant Placentalia (eutherian) that radically departs from morphologically based constructions of the past. Placental mammals are partitioned into four supraordinal clades: Afrotheria, Xenarthra, Laurasiatheria, and Euarchontoglires. Afrotheria form an endemic African clade that includes elephant shrews, golden moles, tenrecs, aardvarks, hyraxes, elephants, dugongs, and manatees. Datamining databases of genome size (GS) shows that till today just one afrotherian GS has been evaluated, that of the aardvark Orycteropus afer. We show that the GSs of six selected representatives across the Afrotheria supraordinal group are among the highest for the extant Placentalia, providing a novel genomic signature of this enigmatic group. The mean GS value of Afrotheria, 5.3 ± 0.7 pg, is the highest reported for the extant Placentalia. This should assist in planning new genome sequencing initiatives. [Reviewing Editor: Dmitri Petrov]     

Calcitonin receptor-stimulating peptide: Its evolutionary and functional relationship with calcitonin/calcitonin gene-related peptide based on gene structure

This review focuses on the evolutionary and functional relationship of calcitonin receptor-stimulating peptide (CRSP) with calcitonin (CT)/calcitonin gene-related peptide (CGRP) in mammals. CRSP shows high sequence identity with CGRP, but distinct biological properties. CRSP genes (CRSPs) have been identified in mammals such as pigs and dogs of the Laurasiatheria, but not in primates and rodents of the Euarchontoglires or in non-placental mammals. CRSPs have genomic organizations highly similar to those of CT/CGRP genes (CT/CGRPs), which are located along with CGRPs in a locus between CYP2R1 and INSC, while the other members of the CGRP superfamily, adrenomedullin and amylin, show genomic organizations and locations distinct from CT, CGRP, and CRSP. Thus, we categorized these three peptides into the CT/CGRP/CRSP family. Non-placental mammals having one and placental mammals having multiple CT/CGRP/CRSP family genes suggests that multiplicity of CT/CGRP started at an early stage of mammalian evolution. In the placental mammals, Laurasiatheria generally possesses multiple CRSPs and only one CT/CGRP, while Euarchontoglires possesses CT/CGRP and CGRPβ but no CRSP, indicating an increase in the diversity and multiplicity of this family of genes in mammalian evolution. Phylogenetic analysis suggests that some CRSPs have been generated very recently in mammalian evolution. Taken together, the increase in the number and complexity of the CT/CGRP/CRSP family genes may have due to evolutionary pressure to facilitate adaptation during mammalian evolution. In this regard, it is important to elucidate the physiological roles of CT, CGRP and CRSP from the viewpoint of the CT/CGRP/CRSP family even in Euarchontoglires.     

Early history of mammals is elucidated with the ENCODE multiple species sequencing data

Understanding the early evolution of placental mammals is one of the most challenging issues in mammalian phylogeny. Here, we addressed this question by using the sequence data of the ENCODE consortium, which include 1% of mammalian genomes in 18 species belonging to all main mammalian lineages. Phylogenetic reconstructions based on an unprecedented amount of coding sequences taken from 218 genes resulted in a highly supported tree placing the root of Placentalia between Afrotheria and Exafroplacentalia (Afrotheria hypothesis). This topology was validated by the phylogenetic analysis of a new class of genomic phylogenetic markers, the conserved noncoding sequences. Applying the tests of alternative topologies on the coding sequence dataset resulted in the rejection of the Atlantogenata hypothesis (Xenarthra grouping with Afrotheria), while this test rejected the second alternative scenario, the Epitheria hypothesis (Xenarthra at the base), when using the noncoding sequence dataset. Thus, the two datasets support the Afrotheria hypothesis; however, none can reject both of the remaining topological alternatives.     

LINE-1 distribution in Afrotheria and Xenarthra: implications for understanding the evolution of LINE-1 in eutherian genomes

Long interspersed nuclear elements (LINEs) comprise about 21% of the human genome (of which L1 is most abundant) and are preferentially accumulated in AT-rich regions, as well as the X and Y chromosomes. Most knowledge of L1 distribution in mammals is restricted to human and mouse. Here we report the first investigation of L1 distribution in the genomes of a wide variety of eutherian mammals, including species in the two basal clades, Afrotheria and Xenarthra. Our results show L1 accumulation on the X of all eutherian mammals, an observation consistent with an ancestral involvement of these elements in the X-inactivation process (the Lyon repeat hypothesis). Surprisingly, conspicuous accumulation of L1 in AT-rich regions of the genome was not observed in any species outside of Euarchontoglires (represented by human, mouse and rabbit). Although several features were common to most species investigated, our comprehensive survey shows that the patterns observed in human and mouse are, in many aspects, far from typical for all mammals. We discuss these findings with reference to models that have previously been proposed to explain the AT distribution bias of L1 in human and mouse, and how this relates to the evolution of these elements in other eutherian genomes.Paul D. Waters and Gauthier Dobigny contributed equally to this work