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.     

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.     

Molecular Evidence for the Major Clades of Placental Mammals

Higher-level relationships among placental mammals, as well as the historical biogeography of this group against the backdrop of continental fragmentation and reassembly, remain poorly understood. Here, we analyze two independent molecular data sets that represent all placental orders. The first data set includes six genes (A2AB, IRBP, vWF, 12S rRNA, tRNA valine, 16S rRNA; total = 5.71 kb) for 26 placental taxa and two marsupials; the second data set includes 2.95 kb of exon 11 of the BRCA1 gene for 51 placental taxa and four marsupials. We also analyzed a concatenation of these data sets (8.66 kb) for 26 placentals and one marsupial. Unrooted and rooted analyses were performed with parsimony, distance methods, maximum likelihood, and a Bayesian approach. Unrooted analyses provide convincing support for a fundamental separation of placental orders into groups with southern and northern hemispheric origins according to the current fossil record. On rooted trees, one or both of these groups are monophyletic depending on the position of the root. Maximum likelihood and Bayesian analyses with the BRCA1 and combined 8.66 kb data sets provide strong support for the monophyly of the northern hemisphere group (Boreoeutheria). Boreoeutheria is divided into Laurasiatheria (Carnivora + Cetartiodactyla + Chiroptera + Eulipotyphla + Perissodactyla + Pholidota) and Euarchonta (Dermoptera + Primates + Scandentia) + Glires (Lagomorpha + Rodentia). The southern hemisphere group is either monophyletic or paraphyletic, depending on the method of analysis used. Within this group, Afrotheria (Proboscidea + Sirenia + Hyracoidea + Tubulidentata + Macroscelidea + Afrosoricida) is monophyletic. A unique nine base-pair deletion in exon 11 of the BRCA1 gene also supports Afrotheria monophyly. Given molecular dates that suggest that the southern hemisphere group and Boreoeutheria diverged in the Early Cretaceous, a single trans-hemispheric dispersal event may have been of fundamental importance in the early history of crown-group Eutheria. Parallel adaptive radiations have subsequently occurred in the four major groups: Laurasiatheria, Euarchonta + Glires, Afrotheria, and Xenarthra.     

A retroposon analysis of Afrotherian phylogeny

Recent comprehensive studies of DNA sequences support the monophyly of Afrotheria, comprising elephants, sirenians (dugongs and manatees), hyraxes, tenrecs, golden moles, aardvarks, and elephant shrews, as well as that of Paenungulata, comprising elephants, sirenians, and hyraxes. However, phylogenetic relationships among paenungulates, as well as among nonpaenungulates, have remained ambiguous. Here we applied an extensive retroposon analysis to these problems to support the monophyly of aardvarks, tenrecs, and golden moles, with elephant shrews as their sister group. Regarding phylogenetic relationships in Paenungulata, we could characterize only one informative locus, although we could isolate many insertions specific to each of three lineages, namely, Proboscidea, Sirenia, and Hyracoidea. These data prompted us to reexamine phylogenetic relationships among Paenungulata using 19 nuclear gene sequences resulting in three different analyses, namely, short interspersed element (SINE) insertions, nuclear sequence analyses, and morphological cladistics, supporting different respective phylogenies. We concluded that these three lineages diverged very rapidly in a very short evolutionary period, with the consequence that ancestral polymorphism present in the last common ancestor of Paenungulata results in such incongruence. Our results suggest the rapid fixation of many large-scale morphological synapomorphies for Tethytheria; implications of this in relation to the morphological evolution in Paenungulata are discussed.     

Patterns of Insertion and Deletion in Mammalian Genomes

Nucleotide insertions and deletions (indels) are responsible for gaps in the sequence alignments. Indel is one of the major sources of evolutionary change at the molecular level. We have examined the patterns of insertions and deletions in the 19 mammalian genomes, and found that deletion events are more common than insertions in the mammalian genomes. Both the number of insertions and deletions decrease rapidly when the gap length increases and single nucleotide indel is the most frequent in all indel events. The frequencies of both insertions and deletions can be described well by power law.Key Words: Insertion, deletion, gap, indel, mammalian genome.     

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.     

Pattern and rate of indel evolution inferred from whole chloroplast intergenic regions in sugarcane, maize and rice.

Microstructural changes such as insertions and deletions (=indels) are a major driving force in the evolution of non-coding DNA sequences. To better understand the mechanisms by which indel mutations arise, as well as the molecular evolution of non-coding regions, the number and pattern of indels and nucleotide substitutions were compared in the whole chloroplast genomes. Comparisons were made for a total of over 38 kb non-coding DNA sequences from 126 intergenic regions in two data sets representing species with different divergence times: sugarcane and maize and Oryza sativa var. indica and japonica. The main findings of this study are: (i) Approximately half of all indels are single nucleotide indels. This observation agrees with previous studies in various organisms. (ii) The distribution and number of indels was different between two data sets, and different patterns were observed for tandem repeat and non-repeat indels. (iii) Distribution pattern of tandem repeat indels showed statistically significant bias towards A/T-rich. (iv) The rate of indel mutation was estimated to be approximately 0.8 +/- 0.04 x 10(-9) per site per year, which was similar to previous estimates in other organisms. (v) The frequencies of nucleotide substitutions and indels were significantly lower in inverted repeat (IR).     

Heterogeneous yet similar introns reside in identical positions of the rRNA genes in natural isolates of the archaeon Aeropyrum pernix

Some archaeal ribosomal DNA (rDNA) introns carry homing endonuclease-like genes and are therefore assumed to propagate by "intron homing". A previous study demonstrated that three introns are located within the rRNA operon (arnSL) of Aeropyrum pernix strain K1, two of which, Ialpha and Igamma, harbor open reading frames (ORFs) encoding putative LAGLIDADG-type endonucleases. In an effort to understand further the rDNA intron distribution in natural A. pernix populations, 11 A. pernix strains were isolated from marine hydrothermal biotopes, and comparative nucleotide sequence analysis of the arnSL alleles was performed. Of the 11 isolates, eight contained multiple introns, and three patterns of intron insertion were found. Three novel introns, Idelta (62 bp in length), Ivarepsilon (122 bp) and Izeta (57 bp) were identified. They were all ORF-less, but their predicted RNA secondary structure at the exon-intron junctions was consistent with the bulge-helix-bulge motif. The insertion positions and the terminal inverted repeat sequences of Idelta and Izeta were in agreement with those of Ialpha and Igamma, respectively. This suggests that these intron variants were generated by large indels (insertions/deletions) during their evolution.     

Molecular phylogenetics of the lizard genus Microlophus (squamata:tropiduridae): aligning and retrieving indel signal from nuclear introns

We use a multigene data set (the mitochondrial locus and nine nuclear gene regions) to test phylogenetic relationships in the South American "lava lizards" (genus Microlophus) and describe a strategy for aligning noncoding sequences that accounts for differences in tempo and class of mutational events. We focus on seven nuclear introns that vary in size and frequency of multibase length mutations (i.e., indels) and present a manual alignment strategy that incorporates insertions and deletions (indels) for each intron. Our method is based on mechanistic explanations of intron evolution that does not require a guide tree. We also use a progressive alignment algorithm (Probabilistic Alignment Kit; PRANK) and distinguishes insertions from deletions and avoids the "gapcost" conundrum. We describe an approach to selecting a guide tree purged of ambiguously aligned regions and use this to refine PRANK performance. We show that although manual alignment is successful in finding repeat motifs and the most obvious indels, some regions can only be subjectively aligned, and there are limits to the size and complexity of a data matrix for which this approach can be taken. PRANK alignments identified more parsimony-informative indels while simultaneously increasing nucleotide identity in conserved sequence blocks flanking the indel regions. When comparing manual and PRANK with two widely used methods (CLUSTAL, MUSCLE) for the alignment of the most length-variable intron, only PRANK recovered a tree congruent at deeper nodes with the combined data tree inferred from all nuclear gene regions. We take this concordance as an objective function of alignment quality and present a strongly supported phylogenetic hypothesis for Microlophus relationships. From this hypothesis we show that (1) a coded indel data partition derived from the PRANK alignment contributed significantly to nodal support and (2) the indel data set permitted detection of significant conflict between mitochondrial and nuclear data partitions, which we hypothesize arose from secondary contact of distantly related taxa, followed by hybridization and mtDNA introgression.     

Selective constraints on intron evolution in Drosophila

Intron sizes show an asymmetrical distribution in a number of organisms, with a large number of "short" introns clustered around a minimal intron length and a much broader distribution of longer introns. In Drosophila melanogaster, the short intron class is centered around 61 bp. The narrow length distribution suggests that natural selection may play a role in maintaining intron size. A comparison of 15 orthologous introns among species of the D. melanogaster subgroup indicates that, in general, short introns are not under greater DNA sequence or length constraints than long introns. There is a bias toward deletions in all introns (deletion/insertion ratio is 1.66), and the vast majority of indels are of short length (<10 bp). Indels occurring on the internal branches of the phylogenetic tree are significantly longer than those occurring on the terminal branches. These results are consistent with a compensatory model of intron length evolution in which slightly deleterious short deletions are frequently fixed within species by genetic drift, and relatively rare larger insertions that restore intron length are fixed by positive selection. A comparison of paralogous introns shared among duplicated genes suggests that length constraints differ between introns within the same gene. The janusA, janusB, and ocnus genes share two short introns derived from a common ancestor. The first of these introns shows significantly fewer indels than the second intron, although the two introns show a comparable number of substitutions. This indicates that intron-specific selective constraints have been maintained following gene duplication, which preceded the divergence of the D. melanogaster species subgroup.     

Characterization and phylogenetic utility of the mammalian protamine p1 gene

We sequenced the protamine P1 gene (ca. 450 bp) from 20 bats (order Chiroptera) and the flying lemur (order Dermoptera). We compared these sequences with published sequences from 19 other mammals representing seven orders (Artiodactyla, Carnivora, Cetacea, Perissodactyla, Primates, Proboscidea, and Rodentia) to assess structure, base compositional bias, and phylogenetic utility. Approximately 80% of second codon positions were guanine, resulting in protamine proteins containing a high frequency of arginine residues. Our data indicate that codon usage for arginine differs among higher mammalian taxa. Parsimony analysis of 40 species representing nine orders produced a well-resolved tree in which most nodes were supported strongly, except at the lowest taxonomic levels (e.g., within Artiodactyla and Vespertilionidae). These data support monophyly of several taxa proposed by morphologic and molecular studies (all nine orders: Laurasiatheria, Cetartiodactytla, Yangochiroptera, Noctilionoidea, Rhinolophoidea, Vespertilionoidea, Phyllostomidae, Natalidae, and Vespertilionidae) and, in agreement with recent molecular studies, reject monophyly of Archonta, Volitantia, and Microchiroptera. Bats were sister to a clade containing Perissodactyla, Carnivora, and Cetartiodactyla, and, although not unequivocally, rhinolophoid bats (traditional microchiropterans) were sister to megachiropterans. Sequences of the protamine P1 gene are useful for resolving relationships at and above the familial level in bats, and generally within and among mammalian orders, but with some drawbacks. The coding and intervening sequences are small, producing few phylogenetically informative characters, and aligning the intron is difficult, even among closely related families. Given these caveats, the protamine P1 gene may be important to future systematic studies because its functional and evolutionary constraints differ from other genes currently used in systematic studies.     

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.     

The genomic landscape of short insertion and deletion polymorphisms in the chicken (Gallus gallus) Genome: a high frequency of deletions in tandem duplicates

It is increasingly recognized that insertions and deletions (indels) are an important source of genetic as well as phenotypic divergence and diversity. We analyzed length polymorphisms identified through partial (0.25x) shotgun sequencing of three breeds of domestic chicken made by the International Chicken Polymorphism Map Consortium. A data set of 140,484 short indel polymorphisms in unique DNA was identified after filtering for microsatellite structures. There was a significant excess of tandem duplicates at indel sites, with deletions of a duplicate motif outnumbering the generation of duplicates through insertion. Indel density was lower in microchromosomes than in macrochromosomes, in the Z chromosome than in autosomes, and in 100 bp of upstream sequence, 5'-UTR, and first introns than in intergenic DNA and in other introns. Indel density was highly correlated with single nucleotide polymorphism (SNP) density. The mean density of indels in pairwise sequence comparisons was 1.9 x 10(-4) indel events/bp, approximately 5% the density of SNPs segregating in the chicken genome. The great majority of indels involved a limited number of nucleotides (median 1 bp), with A-rich motifs being overrepresented at indel sites. The overrepresentation of deletions at tandem duplicates indicates that replication slippage in duplicate sequences is a common mechanism behind indel mutation. The correlation between indel and SNP density indicates common effects of mutation and/or selection on the occurrence of indels and point mutations.     

Polymorphism and structural variation of rps16 group-II intron in the Solanum species

Nucleotide sequences of the self-splicing group-II intron of rps16 have first been determined in nine species of the Solanum genus. It was found that the observed variations in the intron length (855–864 bp) was associated with indels of 1 to 9 bp. Altogether, five indels and 50 nucleotide substitutions were detected, which were used to identify six Solanum haplotypes. Although the intron sequence was in general fairly well conserved, the distribution of the described mutations among its structural elements corresponding to six pre-RNA domains was qualitatively and quantitatively nonuniform. The highest polymorphism levels were observed in domains I, II, and IV. The sequence of domain V was absolutely invariable, which is in agreement with its functional significance. The chloroplast rpS16 intron sequences have been characterized in nine Solanum species. The intron length ranged from 855 bp to 864 bp, which is associated with 1–9-nucleotide indels. In total five indels and 50 nucleotide substitutions have been detected and six Solanum haplotypes have been revealed. Solanum rpS16 introns has been characterized by mutation rate heterogeneity between structure regions of all six domains its pre-RNA. Intron domains I, II, IV are shown to be more variable. Sequences of the domain V are invariant, that agrees with its functional significance.     

Is morphology telling the truth about the evolution of the species rich genus Peperomia (Piperaceae)?

Peperomia is with approximately 1,600 species one of the species rich angiosperm genera. Several characters on which current infrageneric classifications are based are influenced by parallel evolution. A well-resolved molecular backbone phylogeny of the genus is needed to address evolutionary questions about morphological traits. Based on separate and combined analyses of a morphological data set and three molecular data sets, phylogenetic relationships within Peperomia are investigated with respect to character evolution. The resulting trees from different datasets are highly congruent. Morphological characters are mapped on a combined molecular tree, visualizing the contrast between previously used homoplastic characters and some newly observed characters, that can be used to delimit monophyletic groups. Length mutational events of the chloroplast dataset are coded and plotted on the respective tree, to test if indels support alternative hypothesis of relationships found in the nuclear datasets as well as the overall performance of indels compared with substitutional mutations. Our findings indicate that length distribution of indels is highest among five and six bp events. Autapomorphic and synapomorphic length mutations are most frequent in both insertions and deletions and are also more frequent independent of the length of the mutation. Concluding, independent of the length, mutations are of phylogenetic importance and should not be disregarded. None of the homoplastic indels turn into synapomorphic indels, supporting the different topology of the nrDNA tree but indicate areas of molecular evolution in favour of length mutations resulting in independent events.     

Growth and decline of introns

To gauge the processes that might direct the length of introns, I studied the balance of indels (insertions or deletions, determined using Alu and LINE1 retroposon repeats) and the density of these repeats in the introns of the human genome. The indel balance is biased in favour of deletions and correlated with the divergence of repeats. At fixed repeat divergence, the indel bias correlated with the intron size: the shorter the intron, the more deletions were favoured over insertions. This correlation with the intron size was stronger than with the gene-wide or isochore-wide parameters. The density of repeats (the number of repeats in a unit of intron length) correlated positively with the intron size. Thus, quite different mechanisms, the indel bias and the integration and/or persistence of retroposons, act in the same direction in regards to intron size, which suggests selection for the size of individual introns.     

A macaque's-eye view of human insertions and deletions: differences in mechanisms

Insertions and deletions (indels) cause numerous genetic diseases and lead to pronounced evolutionary differences among genomes. The macaque sequences provide an opportunity to gain insights into the mechanisms generating these mutations on a genome-wide scale by establishing the polarity of indels occurring in the human lineage since its divergence from the chimpanzee. Here we apply novel regression techniques and multiscale analyses to demonstrate an extensive regional indel rate variation stemming from local fluctuations in divergence, GC content, male and female recombination rates, proximity to telomeres, and other genomic factors. We find that both replication and, surprisingly, recombination are significantly associated with the occurrence of small indels. Intriguingly, the relative inputs of replication versus recombination differ between insertions and deletions, thus the two types of mutations are likely guided in part by distinct mechanisms. Namely, insertions are more strongly associated with factors linked to recombination, while deletions are mostly associated with replication-related features. Indel as a term misleadingly groups the two types of mutations together by their effect on a sequence alignment. However, here we establish that the correct identification of a small gap as an insertion or a deletion (by use of an outgroup) is crucial to determining its mechanism of origin. In addition to providing novel insights into insertion and deletion mutagenesis, these results will assist in gap penalty modeling and eventually lead to more reliable genomic alignments.     

Weighting indels as phylogenetic markers of 18S rDNA sequences in Diptera and Strepsiptera

Opinions split when it comes to the significance and thus the weighting of indel characters as phylogenetic markers. This paper attempts to test the phylogenetic information content of indels and nucleotide substitutions by proposing an a priori weighting system of non-protein-coding genes. Theoretically, the system rests on a weighting scheme which is based on a falsificationist approach to cladistic inference. It provides insertions, deletions and nucleotide substitutions weights according to their specific number of identical classes of potential falsifiers, resulting in the following system: nucleotide substitutions weight = 3, deletions of n nucleotides weight = (2ⁿ–1), and insertions of n nucleotides weight = (5ⁿ–1). This weighting system and the utility of indels as phylogenetic markers are tested against a suitable data set of 18S rDNA sequences of Diptera and Strepsiptera taxa together with other Metazoa species. The indels support the same clades as the nucleotide substitution data, and the application of the weighting system increases the corresponding consistency indices of the differentially weighted character types. As a consequence, applying the weighting system seems to be reasonable, and indels appear to be good phylogenetic markers.