OrthoMaM v12: a database of curated single-copy ortholog alignments and trees to study mammalian evolutionary genomics.

To date, the databases built to gather information on gene orthology do not provide end-users with descriptors of the molecular evolution information and phylogenetic pattern of these orthologues. In this context, we developed OrthoMaM, a database of ORTHOlogous MAmmalian Markers describing the evolutionary dynamics of coding sequences in mammalian genomes. OrthoMaM version 12 includes 15,868 alignments of orthologous coding sequences (CDS) from the 190 complete mammalian genomes currently available. All annotations and 1-to-1 orthology assignments are based on NCBI. Orthologous CDS can be mined for potential informative markers at the different taxonomic levels of the mammalian tree. To this end, several evolutionary descriptors of DNA sequences are provided for querying purposes (e.g. base composition and relative substitution rate). The graphical web interface allows the user to easily browse and sort the results of combined queries. The corresponding multiple sequence alignments and ML trees, inferred using state-of-the art approaches, are available for download both at the nucleotide and amino acid levels. OrthoMaM v12 can be used by researchers interested either in reconstructing the phylogenetic relationships of mammalian taxa or in understanding the evolutionary dynamics of coding sequences in their genomes. OrthoMaM is available for browsing, querying and complete or filtered download at https://orthomam.mbb.cnrs.fr/.

Exon capture museomics deciphers the nine-banded armadillo species complex and identifies a new species endemic to the Guiana Shield.

The nine-banded armadillo (Dasypus novemcinctus) is the most widespread xenarthran species across the Americas. Recent studies have suggested it is composed of four morphologically and genetically distinct lineages of uncertain taxonomic status. To address this issue, we used a museomic approach to sequence 80 complete mitogenomes and capture 997 nuclear loci for 71 Dasypus individuals sampled across the entire distribution. We carefully cleaned up potential genotyping errors and cross contaminations that could blur species boundaries by mimicking gene flow. Our results unambiguously support four distinct lineages within the D. novemcinctus complex. We found cases of mito-nuclear phylogenetic discordance but only limited contemporary gene flow confined to the margins of the lineage distributions. All available evidence including the restricted gene flow, phylogenetic reconstructions based on both mitogenomes and nuclear loci, and phylogenetic delimitation methods consistently supported the four lineages within D. novemcinctus as four distinct species. Comparable genetic differentiation values to other recognized Dasypus species further reinforced their status as valid species. Considering congruent morphological results from previous studies, we provide an integrative taxonomic view to recognise four species within the D. novemcinctus complex: D. novemcinctus, D. fenestratus, D. mexicanus, and D. guianensis sp. nov., a new species endemic of the Guiana Shield that we describe here. The two available individuals of D. mazzai and D. sabanicola were consistently nested within D. novemcinctus lineage and their status remains to be assessed. The present work offers a case study illustrating the power of museomics to reveal cryptic species diversity within a widely distributed and emblematic species of mammals.

PhylteR: Efficient Identification of Outlier Sequences in Phylogenomic Datasets.

In phylogenomics, incongruences between gene trees, resulting from both artifactual and biological reasons, can decrease the signal-to-noise ratio and complicate species tree inference. The amount of data handled today in classical phylogenomic analyses precludes manual error detection and removal. However, a simple and efficient way to automate the identification of outliers from a collection of gene trees is still missing. Here, we present PhylteR, a method that allows rapid and accurate detection of outlier sequences in phylogenomic datasets, i.e. species from individual gene trees that do not follow the general trend. PhylteR relies on DISTATIS, an extension of multidimensional scaling to 3 dimensions to compare multiple distance matrices at once. In PhylteR, these distance matrices extracted from individual gene phylogenies represent evolutionary distances between species according to each gene. On simulated datasets, we show that PhylteR identifies outliers with more sensitivity and precision than a comparable existing method. We also show that PhylteR is not sensitive to ILS-induced incongruences, which is a desirable feature. On a biological dataset of 14,463 genes for 53 species previously assembled for Carnivora phylogenomics, we show (i) that PhylteR identifies as outliers sequences that can be considered as such by other means, and (ii) that the removal of these sequences improves the concordance between the gene trees and the species tree. Thanks to the generation of numerous graphical outputs, PhylteR also allows for the rapid and easy visual characterization of the dataset at hand, thus aiding in the precise identification of errors. PhylteR is distributed as an R package on CRAN and as containerized versions (docker and singularity).

Pangolin Genomes Offer Key Insights and Resources for the World's Most Trafficked Wild Mammals.

Pangolins form a group of scaly mammals that are trafficked at record numbers for their meat and purported medicinal properties. Despite their conservation concern, knowledge of their evolution is limited by a paucity of genomic data. We aim to produce exhaustive genomic resources that include 3,238 orthologous genes and whole-genome polymorphisms to assess the evolution of all eight extant pangolin species. Robust orthologous gene-based phylogenies recovered the monophyly of the three genera and highlighted the existence of an undescribed species closely related to Southeast Asian pangolins. Signatures of middle Miocene admixture between an extinct, possibly European, lineage and the ancestor of Southeast Asian pangolins, provide new insights into the early evolutionary history of the group. Demographic trajectories and genome-wide heterozygosity estimates revealed contrasts between continental versus island populations and species lineages, suggesting that conservation planning should consider intraspecific patterns. With the expected loss of genomic diversity from recent, extensive trafficking not yet realized in pangolins, we recommend that populations be genetically surveyed to anticipate any deleterious impact of the illegal trade. Finally, we produce a complete set of genomic resources that will be integral for future conservation management and forensic endeavors for pangolins, including tracing their illegal trade. These comprise the completion of whole-genomes for pangolins through the hybrid assembly of the first reference genome for the giant pangolin (Smutsia gigantea) and new draft genomes (∼43x-77x) for four additional species, as well as a database of orthologous genes with over 3.4 million polymorphic sites.

Bayesian total-evidence dating revisits sloth phylogeny and biogeography: a cautionary tale on morphological clock analyses.

Combining morphological and molecular characters through Bayesian total-evidence dating allows inferring the phylogenetic and timescale framework of both extant and fossil taxa, while accounting for the stochasticity and incompleteness of the fossil record. Such an integrative approach is particularly needed when dealing with clades such as sloths (Mammalia: Folivora), for which developmental and biomechanical studies have shown high levels of morphological convergence whereas molecular data can only account for a limited percentage of their total species richness. Here, we propose an alternative hypothesis of sloth evolution that emphasizes the pervasiveness of morphological convergence and the importance of considering the fossil record and an adequate taxon sampling in both phylogenetic and biogeographic inferences. Regardless of different clock models and morphological datasets, the extant sloth Bradypus is consistently recovered as a megatherioid, and Choloepus as a mylodontoid, in agreement with molecular-only analyses. The recently extinct Caribbean sloths (Megalocnoidea) are found to be a monophyletic sister-clade of Megatherioidea, in contrast to previous phylogenetic hypotheses. Our results contradict previous morphological analyses and further support the polyphyly of "Megalonychidae", whose members were found in five different clades. Regardless of taxon sampling and clock models, the Caribbean colonization of sloths is compatible with the exhumation of islands along Aves Ridge and its geological time frame. Overall, our total-evidence analysis illustrates the difficulty of positioning highly incomplete fossils, although a robust phylogenetic framework was recovered by an a posteriori removal of taxa with high percentages of missing characters. Elimination of these taxa improved topological resolution by reducing polytomies and increasing node support. However, it introduced a systematic and geographic bias because most of these incomplete specimens are from northern South America. This is evident in biogeographic reconstructions, which suggest Patagonia as the area of origin of many clades when taxa are underrepresented, but Amazonia and/or Central and Southern Andes when all taxa are included. More generally, our analyses demonstrate the instability of topology and divergence time estimates when using different morphological datasets and clock models, and thus caution against making macroevolutionary inferences when node support is weak or when uncertainties in the fossil record are not considered.

From teeth to pad: tooth loss and development of keratinous structures in sirenians.

Sirenians are a well-known example of morphological adaptation to a shallow-water grazing diet characterized by a modified feeding apparatus and orofacial morphology. Such adaptations were accompanied by an anterior tooth reduction associated with the development of keratinized pads, the evolution of which remains elusive. Among sirenians, the recently extinct Steller's sea cow represents a special case for being completely toothless. Here, we used µ-CT scans of sirenian crania to understand how motor-sensor systems associated with tooth innervation responded to innovations such as keratinized pads and continuous dental replacement. In addition, we surveyed nine genes associated with dental reduction for signatures of loss of function. Our results reveal how patterns of innervation changed with modifications of the dental formula, especially continuous replacement in manatees. Both our morphological and genomic data show that dental development was not completely lost in the edentulous Steller's sea cows. By tracing the phylogenetic history of tooth innervation, we illustrate the role of development in promoting the innervation of keratinized pads, similar to the secondary use of dental canals for innervating neomorphic keratinized structures in other tetrapod groups.

OrthoMaM v10: Scaling-Up Orthologous Coding Sequence and Exon Alignments with More than One Hundred Mammalian Genomes.

We present version 10 of OrthoMaM, a database of orthologous mammalian markers. OrthoMaM is already 11 years old and since the outset it has kept on improving, providing alignments and phylogenetic trees of high-quality computed with state-of-the-art methods on up-to-date data. The main contribution of this version is the increase in the number of taxa: 116 mammalian genomes for 14,509 one-to-one orthologous genes. This has been made possible by the combination of genomic data deposited in Ensembl complemented by additional good-quality genomes only available in NCBI. Version 10 users will benefit from pipeline improvements and a completely redesigned web-interface.

ANISEED 2017: extending the integrated ascidian database to the exploration and evolutionary comparison of genome-scale datasets.

ANISEED (www.aniseed.cnrs.fr) is the main model organism database for tunicates, the sister-group of vertebrates. This release gives access to annotated genomes, gene expression patterns, and anatomical descriptions for nine ascidian species. It provides increased integration with external molecular and taxonomy databases, better support for epigenomics datasets, in particular RNA-seq, ChIP-seq and SELEX-seq, and features novel interactive interfaces for existing and novel datatypes. In particular, the cross-species navigation and comparison is enhanced through a novel taxonomy section describing each represented species and through the implementation of interactive phylogenetic gene trees for 60% of tunicate genes. The gene expression section displays the results of RNA-seq experiments for the three major model species of solitary ascidians. Gene expression is controlled by the binding of transcription factors to cis-regulatory sequences. A high-resolution description of the DNA-binding specificity for 131 Ciona robusta (formerly C. intestinalis type A) transcription factors by SELEX-seq is provided and used to map candidate binding sites across the Ciona robusta and Phallusia mammillata genomes. Finally, use of a WashU Epigenome browser enhances genome navigation, while a Genomicus server was set up to explore microsynteny relationships within tunicates and with vertebrates, Amphioxus, echinoderms and hemichordates.

Odontogenic ameloblast-associated (ODAM) is inactivated in toothless/enamelless placental mammals and toothed whales.

BACKGROUND: The gene for odontogenic ameloblast-associated (ODAM) is a member of the secretory calcium-binding phosphoprotein gene family. ODAM is primarily expressed in dental tissues including the enamel organ and the junctional epithelium, and may also have pleiotropic functions that are unrelated to teeth. Here, we leverage the power of natural selection to test competing hypotheses that ODAM is tooth-specific versus pleiotropic. Specifically, we compiled and screened complete protein-coding sequences, plus sequences for flanking intronic regions, for ODAM in 165 placental mammals to determine if this gene contains inactivating mutations in lineages that either lack teeth (baleen whales, pangolins, anteaters) or lack enamel on their teeth (aardvarks, sloths, armadillos), as would be expected if the only essential functions of ODAM are related to tooth development and the adhesion of the gingival junctional epithelium to the enamel tooth surface. RESULTS: We discovered inactivating mutations in all species of placental mammals that either lack teeth or lack enamel on their teeth. A surprising result is that ODAM is also inactivated in a few additional lineages including all toothed whales that were examined. We hypothesize that ODAM inactivation is related to the simplified outer enamel surface of toothed whales. An alternate hypothesis is that ODAM inactivation in toothed whales may be related to altered antimicrobial functions of the junctional epithelium in aquatic habitats. Selection analyses on ODAM sequences revealed that the composite dN/dS value for pseudogenic branches is close to 1.0 as expected for a neutrally evolving pseudogene. DN/dS values on transitional branches were used to estimate ODAM inactivation times. In the case of pangolins, ODAM was inactivated ~ 65 million years ago, which is older than the oldest pangolin fossil (Eomanis, 47 Ma) and suggests an even more ancient loss or simplification of teeth in this lineage. CONCLUSION: Our results validate the hypothesis that the only essential functions of ODAM that are maintained by natural selection are related to tooth development and/or the maintenance of a healthy junctional epithelium that attaches to the enamel surface of teeth.

MACSE v2: Toolkit for the Alignment of Coding Sequences Accounting for Frameshifts and Stop Codons.

Multiple sequence alignment is a prerequisite for many evolutionary analyses. Multiple Alignment of Coding Sequences (MACSE) is a multiple sequence alignment program that explicitly accounts for the underlying codon structure of protein-coding nucleotide sequences. Its unique characteristic allows building reliable codon alignments even in the presence of frameshifts. This facilitates downstream analyses such as selection pressure estimation based on the ratio of nonsynonymous to synonymous substitutions. Here, we present MACSE v2, a major update with an improved version of the initial algorithm enriched with a complete toolkit to handle multiple alignments of protein-coding sequences. A graphical interface now provides user-friendly access to the different subprograms.

Convergent Acquisition of Nonembryonic Development in Styelid Ascidians.

Asexual propagation and whole body regeneration are forms of nonembryonic development (NED) widespread across animal phyla and central in life history and evolutionary diversification of metazoans. Whereas it is challenging to reconstruct the gains or losses of NED at large phylogenetic scale, comparative studies could benefit from being conducted at more restricted taxonomic scale, in groups for which phylogenetic relationships are well established. The ascidian family of Styelidae encompasses strictly sexually reproducing solitary forms as well as colonial species that combine sexual reproduction with different forms of NED. To date, the phylogenetic relationships between colonial and solitary styelids remain controversial and so is the pattern of NED evolution. In this study, we built an original pipeline to combine eight genomes with 18 de novo assembled transcriptomes and constructed data sets of unambiguously orthologous genes. Using a phylogenomic super-matrix of 4,908 genes from these 26 tunicates we provided a robust phylogeny of this family of chordates, which supports two convergent acquisitions of NED. This result prompted us to further describe the budding process in the species Polyandrocarpa zorritensis, leading to the discovery of a novel mechanism of asexual development. Whereas the pipeline and the data sets produced can be used for further phylogenetic reconstructions in tunicates, the phylogeny provided here sets an evolutionary framework for future experimental studies on the emergence and disappearance of complex characters such as asexual propagation and whole body regeneration.

A phylogenomic framework and timescale for comparative studies of tunicates.

BACKGROUND: Tunicates are the closest relatives of vertebrates and are widely used as models to study the evolutionary developmental biology of chordates. Their phylogeny, however, remains poorly understood, and to date, only the 18S rRNA nuclear gene and mitogenomes have been used to delineate the major groups of tunicates. To resolve their evolutionary relationships and provide a first estimate of their divergence times, we used a transcriptomic approach to build a phylogenomic dataset including all major tunicate lineages, consisting of 258 evolutionarily conserved orthologous genes from representative species. RESULTS: Phylogenetic analyses using site-heterogeneous CAT mixture models of amino acid sequence evolution resulted in a strongly supported tree topology resolving the relationships among four major tunicate clades: (1) Appendicularia, (2) Thaliacea + Phlebobranchia + Aplousobranchia, (3) Molgulidae, and (4) Styelidae + Pyuridae. Notably, the morphologically derived Thaliacea are confirmed as the sister group of the clade uniting Phlebobranchia + Aplousobranchia within which the precise position of the model ascidian genus Ciona remains uncertain. Relaxed molecular clock analyses accommodating the accelerated evolutionary rate of tunicates reveal ancient diversification (~ 450-350 million years ago) among the major groups and allow one to compare their evolutionary age with respect to the major vertebrate model lineages. CONCLUSIONS: Our study represents the most comprehensive phylogenomic dataset for the main tunicate lineages. It offers a reference phylogenetic framework and first tentative timescale for tunicates, allowing a direct comparison with vertebrate model species in comparative genomics and evolutionary developmental biology studies.

Resolving the phylogenetic position of Darwin's extinct ground sloth (Mylodon darwinii) using mitogenomic and nuclear exon data.

Mylodon darwinii is the extinct giant ground sloth named after Charles Darwin, who first collected its remains in South America. We have successfully obtained a high-quality mitochondrial genome at 99-fold coverage using an Illumina shotgun sequencing of a 12 880-year-old bone fragment from Mylodon Cave in Chile. Low level of DNA damage showed that this sample was exceptionally well preserved for an ancient subfossil, probably the result of the dry and cold conditions prevailing within the cave. Accordingly, taxonomic assessment of our shotgun metagenomic data showed a very high percentage of endogenous DNA with 22% of the assembled metagenomic contigs assigned to Xenarthra. Additionally, we enriched over 15 kb of sequence data from seven nuclear exons, using target sequence capture designed against a wide xenarthran dataset. Phylogenetic and dating analyses of the mitogenomic dataset including all extant species of xenarthrans and the assembled nuclear supermatrix unambiguously place Mylodon darwinii as the sister-group of modern two-fingered sloths, from which it diverged around 22 million years ago. These congruent results from both the mitochondrial and nuclear data support the diphyly of the two modern sloth lineages, implying the convergent evolution of their unique suspensory behaviour as an adaption to arboreality. Our results offer promising perspectives for whole-genome sequencing of this emblematic extinct taxon.

ANISEED 2015: a digital framework for the comparative developmental biology of ascidians.

Ascidians belong to the tunicates, the sister group of vertebrates and are recognized model organisms in the field of embryonic development, regeneration and stem cells. ANISEED is the main information system in the field of ascidian developmental biology. This article reports the development of the system since its initial publication in 2010. Over the past five years, we refactored the system from an initial custom schema to an extended version of the Chado schema and redesigned all user and back end interfaces. This new architecture was used to improve and enrich the description of Ciona intestinalis embryonic development, based on an improved genome assembly and gene model set, refined functional gene annotation, and anatomical ontologies, and a new collection of full ORF cDNAs. The genomes of nine ascidian species have been sequenced since the release of the C. intestinalis genome. In ANISEED 2015, all nine new ascidian species can be explored via dedicated genome browsers, and searched by Blast. In addition, ANISEED provides full functional gene annotation, anatomical ontologies and some gene expression data for the six species with highest quality genomes. ANISEED is publicly available at: http://www.aniseed.cnrs.fr.

Shotgun Mitogenomics Provides a Reference Phylogenetic Framework and Timescale for Living Xenarthrans.

Xenarthra (armadillos, sloths, and anteaters) constitutes one of the four major clades of placental mammals. Despite their phylogenetic distinctiveness in mammals, a reference phylogeny is still lacking for the 31 described species. Here we used Illumina shotgun sequencing to assemble 33 new complete mitochondrial genomes, establishing Xenarthra as the first major placental clade to be fully sequenced at the species level for mitogenomes. The resulting data set allowed the reconstruction of a robust phylogenetic framework and timescale that are consistent with previous studies conducted at the genus level using nuclear genes. Incorporating the full species diversity of extant xenarthrans points to a number of inconsistencies in xenarthran systematics and species definition. We propose to split armadillos into two distinct families Dasypodidae (dasypodines) and Chlamyphoridae (euphractines, chlamyphorines, and tolypeutines) to better reflect their ancient divergence, estimated around 42 Ma. Species delimitation within long-nosed armadillos (genus Dasypus) appeared more complex than anticipated, with the discovery of a divergent lineage in French Guiana. Diversification analyses showed Xenarthra to be an ancient clade with a constant diversification rate through time with a species turnover driven by high but constant extinction. We also detected a significant negative correlation between speciation rate and past temperature fluctuations with an increase in speciation rate corresponding to the general cooling observed during the last 15 My. Biogeographic reconstructions identified the tropical rainforest biome of Amazonia and the Guiana Shield as the cradle of xenarthran evolutionary history with subsequent dispersions into more open and dry habitats.

Genome-Wide Screening of Retroviral Envelope Genes in the Nine-Banded Armadillo (Dasypus novemcinctus, Xenarthra) Reveals an Unfixed Chimeric Endogenous Betaretrovirus Using the ASCT2 Receptor.

UNLABELLED: Retroviruses enter host cells through the interaction of their envelope (Env) protein with a cell surface receptor, which triggers the fusion of viral and cellular membranes. The sodium-dependent neutral amino acid transporter ASCT2 is the common receptor of the large RD114 retrovirus interference group, whose members display frequent env recombination events. Germ line retrovirus infections have led to numerous inherited endogenous retroviruses (ERVs) in vertebrate genomes, which provide useful insights into the coevolutionary history of retroviruses and their hosts. Rare ERV-derived genes display conserved viral functions, as illustrated by the fusogenic syncytin env genes involved in placentation. Here, we searched for functional env genes in the nine-banded armadillo (Dasypus novemcinctus) genome and identified dasy-env1.1, which clusters with RD114 interference group env genes and with two syncytin genes sharing ASCT2 receptor usage. Using ex vivo pseudotyping and cell-cell fusion assays, we demonstrated that the Dasy-Env1.1 protein is fusogenic and can use both human and armadillo ASCT2s as receptors. This gammaretroviral env gene belongs to a provirus with betaretrovirus-like features, suggesting acquisition through recombination. Provirus insertion was found in several Dasypus species, where it has not reached fixation, whereas related family members integrated before diversification of the genus Dasypus >12 million years ago (Mya). This newly described ERV lineage is potentially useful as a population genetic marker. Our results extend the usage of ASCT2 as a retrovirus receptor to the mammalian clade Xenarthra and suggest that the acquisition of an ASCT2-interacting env gene is a major selective force driving the emergence of numerous chimeric viruses in vertebrates. IMPORTANCE: Retroviral infection is initiated by the binding of the viral envelope glycoprotein to a host cell receptor(s), triggering membrane fusion. Ancient germ line infections have generated numerous endogenous retroviruses (ERVs) in nearly all vertebrate genomes. Here, we report a previously uncharacterized ERV lineage from the genome of a xenarthran species, the nine-banded armadillo (Dasypus novemcinctus). It entered the Dasypus genus >12 Mya, with one element being inserted more recently in some Dasypus species, where it could serve as a useful marker for population genetics. This element exhibits an env gene, acquired by recombination events, with conserved viral fusogenic properties through binding to ASCT2, a receptor used by a wide range of recombinant retroviruses infecting other vertebrate orders. This specifies the ASCT2 transporter as a successful receptor for ERV endogenization and suggests that ASCT2-binding env acquisition events have favored the emergence of numerous chimeric viruses in a wide range of species.

Genetic structuring in a relictual population of screaming hairy armadillo (Chaetophractus vellerosus) in Argentina revealed by a set of novel microsatellite loci.

The screaming hairy armadillo (Chaetophractus vellerosus) is a mammal species containing disjunct and isolated populations. In order to assess the effect of habitat fragmentation and geographic isolation, we developed seven new microsatellite loci isolated from low-coverage genome shotgun sequencing data for this species. Among these loci, six microsatellites were found to be polymorphic with 8-26 alleles per locus detected across 69 samples analyzed from a relictual population of the species located in the northeast of the Buenos Aires Province (Argentina). Mean allelic richness and polymorphic information content were 15 and 0.75, with observed and expected heterozygosities ranging from 0.40 to 0.67 and 0.58 to 0.90, respectively. All loci showed departures from Hardy-Weinberg equilibrium. The analysis of population structure in this relictual population revealed three groups of individuals that are genetically differentiated. These newly developed microsatellites will constitute a very useful tool for the estimation of genetic diversity and structure, population dynamics, social structure, parentage and mating system in this little-studied armadillo species. Such genetic data will be particularly helpful for the development of conservation strategies for this isolated population and also for the endangered Bolivian populations previously recognized as a distinct species (Chaetophractus nationi).

Predicting biotic interactions and their variability in a changing environment.

Global environmental change is altering the patterns of biodiversity worldwide. Observation and theory suggest that species' distributions and abundances depend on a suite of processes, notably abiotic filtering and biotic interactions, both of which are constrained by species' phylogenetic history. Models predicting species distribution have historically mostly considered abiotic filtering and are only starting to integrate biotic interaction. However, using information on present interactions to forecast the future of biodiversity supposes that biotic interactions will not change when species are confronted with new environments. Using bacterial microcosms, we illustrate how biotic interactions can vary along an environmental gradient and how this variability can depend on the phylogenetic distance between interacting species.

Evolutionary analysis of selective constraints identifies ameloblastin (AMBN) as a potential candidate for amelogenesis imperfecta.

BACKGROUND: Ameloblastin (AMBN) is a phosphorylated, proline/glutamine-rich protein secreted during enamel formation. Previous studies have revealed that this enamel matrix protein was present early in vertebrate evolution and certainly plays important roles during enamel formation although its precise functions remain unclear. We performed evolutionary analyses of AMBN in order to (i) identify residues and motifs important for the protein function, (ii) predict mutations responsible for genetic diseases, and (iii) understand its molecular evolution in mammals. RESULTS: In silico searches retrieved 56 complete sequences in public databases that were aligned and analyzed computationally. We showed that AMBN is globally evolving under moderate purifying selection in mammals and contains a strong phylogenetic signal. In addition, our analyses revealed codons evolving under significant positive selection. Evidence for positive selection acting on AMBN was observed in catarrhine primates and the aye-aye. We also found that (i) an additional translation initiation site was recruited in the ancestral placental AMBN, (ii) a short exon was duplicated several times in various species including catarrhine primates, and (iii) several polyadenylation sites are present. CONCLUSIONS: AMBN possesses many positions, which have been subjected to strong selective pressure for 200 million years. These positions correspond to several cleavage sites and hydroxylated, O-glycosylated, and phosphorylated residues. We predict that these conserved positions would be potentially responsible for enamel disorder if substituted. Some motifs that were previously identified as potentially important functionally were confirmed, and we found two, highly conserved, new motifs, the function of which should be tested in the near future. This study illustrates the power of evolutionary analyses for characterizing the functional constraints acting on proteins with yet uncharacterized structure.

OrthoMaM v8: a database of orthologous exons and coding sequences for comparative genomics in mammals.

Comparative genomic studies extensively rely on alignments of orthologous sequences. Yet, selecting, gathering, and aligning orthologous exons and protein-coding sequences (CDS) that are relevant for a given evolutionary analysis can be a difficult and time-consuming task. In this context, we developed OrthoMaM, a database of ORTHOlogous MAmmalian Markers describing the evolutionary dynamics of orthologous genes in mammalian genomes using a phylogenetic framework. Since its first release in 2007, OrthoMaM has regularly evolved, not only to include newly available genomes but also to incorporate up-to-date software in its analytic pipeline. This eighth release integrates the 40 complete mammalian genomes available in Ensembl v73 and provides alignments, phylogenies, evolutionary descriptor information, and functional annotations for 13,404 single-copy orthologous CDS and 6,953 long exons. The graphical interface allows to easily explore OrthoMaM to identify markers with specific characteristics (e.g., taxa availability, alignment size, %G+C, evolutionary rate, chromosome location). It hence provides an efficient solution to sample preprocessed markers adapted to user-specific needs. OrthoMaM has proven to be a valuable resource for researchers interested in mammalian phylogenomics, evolutionary genomics, and has served as a source of benchmark empirical data sets in several methodological studies. OrthoMaM is available for browsing, query and complete or filtered downloads at http://www.orthomam.univ-montp2.fr/.