Parallel adaptive radiations in two major clades of placental mammals

Bootstrapping phylogenies inferred from rearrangement data

BACKGROUND

Large-scale sequencing of genomes has enabled the inference of phylogenies based on the evolution of genomic architecture, under such events as rearrangements, duplications, and losses. Many evolutionary models and associated algorithms have been designed over the last few years and have found use in comparative genomics and phylogenetic inference. However, the assessment of phylogenies built from such data has not been properly addressed to date. The standard method used in sequence-based phylogenetic inference is the bootstrap, but it relies on a large number of homologous characters that can be resampled; yet in the case of rearrangements, the entire genome is a single character. Alternatives such as the jackknife suffer from the same problem, while likelihood tests cannot be applied in the absence of well established probabilistic models.

RESULTS

We present a new approach to the assessment of distance-based phylogenetic inference from whole-genome data; our approach combines features of the jackknife and the bootstrap and remains nonparametric. For each feature of our method, we give an equivalent feature in the sequence-based framework; we also present the results of extensive experimental testing, in both sequence-based and genome-based frameworks. Through the feature-by-feature comparison and the experimental results, we show that our bootstrapping approach is on par with the classic phylogenetic bootstrap used in sequence-based reconstruction, and we establish the clear superiority of the classic bootstrap for sequence data and of our corresponding new approach for rearrangement data over proposed variants. Finally, we test our approach on a small dataset of mammalian genomes, verifying that the support values match current thinking about the respective branches.

CONCLUSIONS

Our method is the first to provide a standard of assessment to match that of the classic phylogenetic bootstrap for aligned sequences. Its support values follow a similar scale and its receiver-operating characteristics are nearly identical, indicating that it provides similar levels of sensitivity and specificity. Thus our assessment method makes it possible to conduct phylogenetic analyses on whole genomes with the same degree of confidence as for analyses on aligned sequences. Extensions to search-based inference methods such as maximum parsimony and maximum likelihood are possible, but remain to be thoroughly tested.

Impact of missing data on phylogenies inferred from empirical phylogenomic data sets

Progress in sequencing technology allows researchers to assemble ever-larger supermatrices for phylogenomic inference. However, current phylogenomic studies often rest on patchy data sets, with some having 80% missing (or ambiguous) data or more. Though early simulations had suggested that missing data per se do not harm phylogenetic inference when using sufficiently large data sets, Lemmon et al. (Lemmon AR, Brown JM, Stanger-Hall K, Lemmon EM. 2009. The effect of ambiguous data on phylogenetic estimates obtained by maximum likelihood and Bayesian inference. Syst Biol. 58:130-145.) have recently cast doubt on this consensus in a study based on the introduction of parsimony-uninformative incomplete characters. In this work, we empirically reassess the issue of missing data in phylogenomics while exploring possible interactions with the model of sequence evolution. First, we note that parsimony-uninformative incomplete characters are actually informative in a probabilistic framework. A reanalysis of Lemmon's data set with this in mind gives a very different interpretation of their results and shows that some of their conclusions may be unfounded. Second, we investigate the effect of the progressive introduction of missing data in a complete supermatrix (126 genes × 39 species) capable of resolving animal relationships. These analyses demonstrate that missing data perturb phylogenetic inference slightly beyond the expected decrease in resolving power. In particular, they exacerbate systematic errors by reducing the number of species effectively available for the detection of multiple substitutions. Consequently, large sparse supermatrices are more sensitive to phylogenetic artifacts than smaller but less incomplete data sets, which argue for experimental designs aimed at collecting a modest number (~50) of highly covered genes. Our results further confirm that including incomplete yet short-branch taxa (i.e., slowly evolving species or close outgroups) can help to eschew artifacts, as predicted by simulations. Finally, it appears that selecting an adequate model of sequence evolution (e.g., the site-heterogeneous CAT model instead of the site-homogeneous WAG model) is more beneficial to phylogenetic accuracy than reducing the level of missing data.

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.

The immune gene repertoire of an important viral reservoir, the Australian black flying fox

Background

Bats are the natural reservoir host for a range of emerging and re-emerging viruses, including SARS-like coronaviruses, Ebola viruses, henipaviruses and Rabies viruses. However, the mechanisms responsible for the control of viral replication in bats are not understood and there is little information available on any aspect of antiviral immunity in bats. Massively parallel sequencing of the bat transcriptome provides the opportunity for rapid gene discovery. Although the genomes of one megabat and one microbat have now been sequenced to low coverage, no transcriptomic datasets have been reported from any bat species. In this study, we describe the immune transcriptome of the Australian flying fox, Pteropus alecto, providing an important resource for identification of genes involved in a range of activities including antiviral immunity.

Results

Towards understanding the adaptations that have allowed bats to coexist with viruses, we have de novo assembled transcriptome sequence from immune tissues and stimulated cells from P. alecto. We identified about 18,600 genes involved in a broad range of activities with the most highly expressed genes involved in cell growth and maintenance, enzyme activity, cellular components and metabolism and energy pathways. 3.5% of the bat transcribed genes corresponded to immune genes and a total of about 500 immune genes were identified, providing an overview of both innate and adaptive immunity. A small proportion of transcripts found no match with annotated sequences in any of the public databases and may represent bat-specific transcripts.

Conclusions

This study represents the first reported bat transcriptome dataset and provides a survey of expressed bat genes that complement existing bat genomic data. In addition, these data provide insight into genes relevant to the antiviral responses of bats, and form a basis for examining the roles of these molecules in immune response to viral infection.

Morphology and ultrastructure of the chorioallantoic placenta of the Iberian mole (Talpa occidentalis) with special reference to heterophagous areolas and the nature of interhaemal barrier

This study provides a contribution to the reconstruction of the eulipotyphlan placental morphotype and also may help resolving a long-standing conflict about the interhaemal barrier in moles. As detailed descriptions of talpid placentation, only available for Talpa europaea and Scalopus aquaticus, led to a controversial debate on the nature of interhaemal barrier, the collection of more placental data of further mole species was strongly desired. Hence, the placentas of six gestational stages of Talpa occidentalis have been studied concerning their morphogenesis and ultrastructure with special focus on the structure of the interhaemal barrier and heterophagous regions. Generally, the mode of placentation in T. occidentalis resembles that of T. europaea, including a broad, discoid, antimesometrial, definitive chorioallantoic placenta of labyrinthine type being still villous in earlier stages. Within the labyrinth, the zona intima shows an endotheliochorial interhaemal barrier with a two-layered trophoblast. This clearly contradicts former statements on the S. aquaticus placenta made by Prasad et al. (1979), although their findings cannot exclude a totally different interpretation. Regardless, the placenta of moles represents the least invasive mode of placentation among Eulipotyphla, which otherwise have highly invasive placentas. Although the phagocytic areolas situated above uterine gland openings are heterophagous, they mainly seem to serve fetal histiotrophic nutrition, at least early in pregnancy. In later stages the number of glands and areolas decreases. This special type of additional phagocytic region is usually most common in species with noninvasive, epitheliochorial placentation, which suggests a correlation between placental invasiveness and the occurrence and type of phagocytic placental structures. The compact and invasive mode of placentation of Talpidae and all other Eulipotyphla seems to be plesiomorphic within Laurasiatheria and is always correlated with an altricial neonate.

Chromosome evolution in Xenarthra: new insights from an ancient group

The Magnaorder Xenarthra is one of the four main supraordinal eutherian clades, together with Afrotheria, Euarchontoglires and Laurasiatheria. Xenarthra is an eminently Central and South American group of special interest in phylogenetic studies due to its possible position at the base of the eutherian tree. The use of modern cytogenetic techniques in some species of Xenarthra has provided important insights into the karyotypic evolution of mammals. Nevertheless, chromosome analyses in the group are still restricted, with only a few individuals of each species studied and karyotype descriptions mostly without banding patterns. In addition, it is likely that still unknown species exist and that the chromosome variability in the group is underestimated. We present a review of the currently available data on Xenarthra chromosomes and genomes and on the impact that their study has had in the understanding of mammalian genome evolution. It is clear that further cytogenetic analyses in Xenarthra, including banding patterns and molecular approaches, are likely to help in the identification of new species, reveal still undetected chromosome variations, provide information to support conservation strategies planning, and greatly contribute to a better understanding of mammalian genome evolution.

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.

Rooting gene trees without outgroups: EP rooting

Gene sequences are routinely used to determine the topologies of unrooted phylogenetic trees, but many of the most important questions in evolution require knowing both the topologies and the roots of trees. However, general algorithms for calculating rooted trees from gene and genomic sequences in the absence of gene paralogs are few. Using the principles of evolutionary parsimony (EP) (Lake JA. 1987a. A rate-independent technique for analysis of nucleic acid sequences: evolutionary parsimony. Mol Biol Evol. 4:167–181) and its extensions (Cavender, J. 1989. Mechanized derivation of linear invariants. Mol Biol Evol. 6:301–316; Nguyen T, Speed TP. 1992. A derivation of all linear invariants for a nonbalanced transversion model. J Mol Evol. 35:60–76), we explicitly enumerate all linear invariants that solely contain rooting information and derive algorithms for rooting gene trees directly from gene and genomic sequences. These new EP linear rooting invariants allow one to determine rooted trees, even in the complete absence of outgroups and gene paralogs. EP rooting invariants are explicitly derived for three taxon trees, and rules for their extension to four or more taxa are provided. The method is demonstrated using 18S ribosomal DNA to illustrate how the new animal phylogeny (Aguinaldo AMA et al. 1997. Evidence for a clade of nematodes, arthropods, and other moulting animals. Nature 387:489–493; Lake JA. 1990. Origin of the metazoa. Proc Natl Acad Sci USA 87:763–766) may be rooted directly from sequences, even when they are short and paralogs are unavailable. These results are consistent with the current root (Philippe H et al. 2011. Acoelomorph flatworms are deuterostomes related to Xenoturbella. Nature 470:255–260).

Chromosome painting in three-toed sloths: a cytogenetic signature and ancestral karyotype for Xenarthra

BACKGROUND

Xenarthra (sloths, armadillos and anteaters) represent one of four currently recognized Eutherian mammal supraorders. Some phylogenomic studies point to the possibility of Xenarthra being at the base of the Eutherian tree, together or not with the supraorder Afrotheria. We performed painting with human autosomes and X-chromosome specific probes on metaphases of two three-toed sloths: Bradypus torquatus and B. variegatus. These species represent the fourth of the five extant Xenarthra families to be studied with this approach.

RESULTS

Eleven human chromosomes were conserved as one block in both B. torquatus and B. variegatus: (HSA 5, 6, 9, 11, 13, 14, 15, 17, 18, 20, 21 and the X chromosome). B. torquatus, three additional human chromosomes were conserved intact (HSA 1, 3 and 4). The remaining human chromosomes were represented by two or three segments on each sloth. Seven associations between human chromosomes were detected in the karyotypes of both B. torquatus and B. variegatus: HSA 3/21, 4/8, 7/10, 7/16, 12/22, 14/15 and 17/19. The ancestral Eutherian association 16/19 was not detected in the Bradypus species.

CONCLUSIONS

Our results together with previous reports enabled us to propose a hypothetical ancestral Xenarthran karyotype with 48 chromosomes that would differ from the proposed ancestral Eutherian karyotype by the presence of the association HSA 7/10 and by the split of HSA 8 into three blocks, instead of the two found in the Eutherian ancestor. These same chromosome features point to the monophyly of Xenarthra, making this the second supraorder of placental mammals to have a chromosome signature supporting its monophyly.

A comparative study of mammalian diversification pattern

Although mammals have long been regarded as a successful radiation, the diversification pattern among the clades is still poorly known. Higher-level phylogenies are conflicting and comprehensive comparative analyses are still lacking. Using a recently published supermatrix encompassing nearly all extant mammalian families and a novel comparative likelihood approach (MEDUSA), the diversification pattern of mammalian groups was examined. Both order- and family-level phylogenetic analyses revealed the rapid radiation of Boreoeutheria and Euaustralidelphia in the early mammalian history. The observation of a diversification burst within Boreoeutheria at approximately 100 My supports the Long Fuse model in elucidating placental diversification progress, and the rapid radiation of Euaustralidelphia suggests an important role of biogeographic dispersal events in triggering early Australian marsupial rapid radiation. Diversification analyses based on family-level diversity tree revealed seven additional clades with exceptional diversification rate shifts, six of which represent accelerations in net diversification rate as compared to the background pattern. The shifts gave origin to the clades Muridae+Cricetidae, Bovidae+Moschidae+Cervidae, Simiiformes, Echimyidae, Odontoceti (excluding Physeteridae+Kogiidae+Platanistidae), Macropodidae, and Vespertilionidae. Moderate to high extinction rates from background and boreoeutherian diversification patterns indicate the important role of turnovers in shaping the heterogeneous taxonomic richness observed among extant mammalian groups. Furthermore, the present results emphasize the key role of extinction on erasing unusual diversification signals, and suggest that further studies are needed to clarify the historical radiation of some mammalian groups for which MEDUSA did not detect exceptional diversification rates.

Joint reconstruction of divergence times and life-history evolution in placental mammals using a phylogenetic covariance model

Violation of the molecular clock has been amply documented, and is now routinely taken into account by molecular dating methods. Comparative analyses have revealed a systematic component in rate variation, relating it to the evolution of life-history traits, such as body size or generation time. Life-history evolution can be reconstructed using Brownian models. However, the resulting estimates are typically uncertain, and potentially sensitive to the underlying assumptions. As a way of obtaining more accurate ancestral trait and divergence time reconstructions, correlations between life-history traits and substitution rates could be used as an additional source of information. In this direction, a Bayesian framework for jointly reconstructing rates, traits, and dates was previously introduced. Here, we apply this model to a 17 protein-coding gene alignment for 73 placental taxa. Our analysis indicates that the coupling between molecules and life history can lead to a reevaluation of ancestral life-history profiles, in particular for groups displaying convergent evolution in body size. However, reconstructions are sensitive to fossil calibrations and to the Brownian assumption. Altogether, our analysis suggests that further integrating inference of rates and traits might be particularly useful for neontological macroevolutionary comparative studies.

Fast and accurate phylogenetic reconstruction from high-resolution whole-genome data and a novel robustness estimator

The rapid accumulation of whole-genome data has renewed interest in the study of genomic rearrangements. Comparative genomics, evolutionary biology, and cancer research all require models and algorithms to elucidate the mechanisms, history, and consequences of these rearrangements. However, even simple models lead to NP-hard problems, particularly in the area of phylogenetic analysis. Current approaches are limited to small collections of genomes and low-resolution data (typically a few hundred syntenic blocks). Moreover, whereas phylogenetic analyses from sequence data are deemed incomplete unless bootstrapping scores (a measure of confidence) are given for each tree edge, no equivalent to bootstrapping exists for rearrangement-based phylogenetic analysis. We describe a fast and accurate algorithm for rearrangement analysis that scales up, in both time and accuracy, to modern high-resolution genomic data. We also describe a novel approach to estimate the robustness of results-an equivalent to the bootstrapping analysis used in sequence-based phylogenetic reconstruction. We present the results of extensive testing on both simulated and real data showing that our algorithm returns very accurate results, while scaling linearly with the size of the genomes and cubically with their number. We also present extensive experimental results showing that our approach to robustness testing provides excellent estimates of confidence, which, moreover, can be tuned to trade off thresholds between false positives and false negatives. Together, these two novel approaches enable us to attack heretofore intractable problems, such as phylogenetic inference for high-resolution vertebrate genomes, as we demonstrate on a set of six vertebrate genomes with 8,380 syntenic blocks. A copy of the software is available on demand.

Phylogenomic analysis resolves the interordinal relationships and rapid diversification of the laurasiatherian mammals

Although great progress has been made in resolving the relationships of placental mammals, the position of several clades in Laurasiatheria remain controversial. In this study, we performed a phylogenetic analysis of 97 orthologs (46,152 bp) for 15 taxa, representing all laurasiatherian orders. Additionally, phylogenetic trees of laurasiatherian mammals with draft genome sequences were reconstructed based on 1608 exons (2,175,102 bp). Our reconstructions resolve the interordinal relationships within Laurasiatheria and corroborate the clades Scrotifera, Fereuungulata, and Cetartiodactyla. Furthermore, we tested alternative topologies within Laurasiatheria, and among alternatives for the phylogenetic position of Perissodactyla, a sister-group relationship with Cetartiodactyla receives the highest support. Thus, Pegasoferae (Perissodactyla + Carnivora + Pholidota + Chiroptera) does not appear to be a natural group. Divergence time estimates from these genes were compared with published estimates for splits within Laurasiatheria. Our estimates were similar to those of several studies and suggest that the divergences among these orders occurred within just a few million years.

Ultraconserved elements are novel phylogenomic markers that resolve placental mammal phylogeny when combined with species-tree analysis

Phylogenomics offers the potential to fully resolve the Tree of Life, but increasing genomic coverage also reveals conflicting evolutionary histories among genes, demanding new analytical strategies for elucidating a single history of life. Here, we outline a phylogenomic approach using a novel class of phylogenetic markers derived from ultraconserved elements and flanking DNA. Using species-tree analysis that accounts for discord among hundreds of independent loci, we show that this class of marker is useful for recovering deep-level phylogeny in placental mammals. In broad outline, our phylogeny agrees with recent phylogenomic studies of mammals, including several formerly controversial relationships. Our results also inform two outstanding questions in placental mammal phylogeny involving rapid speciation, where species-tree methods are particularly needed. Contrary to most phylogenomic studies, our study supports a first-diverging placental mammal lineage that includes elephants and tenrecs (Afrotheria). The level of conflict among gene histories is consistent with this basal divergence occurring in or near a phylogenetic "anomaly zone" where a failure to account for coalescent stochasticity will mislead phylogenetic inference. Addressing a long-standing phylogenetic mystery, we find some support from a high genomic coverage data set for a traditional placement of bats (Chiroptera) sister to a clade containing Perissodactyla, Cetartiodactyla, and Carnivora, and not nested within the latter clade, as has been suggested recently, although other results were conflicting. One of the most remarkable findings of our study is that ultraconserved elements and their flanking DNA are a rich source of phylogenetic information with strong potential for application across Amniotes.

Sequencing of 21 varicella-zoster virus genomes reveals two novel genotypes and evidence of recombination

Genotyping of 21 varicella-zoster virus (VZV) strains using a scattered single nucleotide polymorphism (SNP) method revealed ambiguous SNPs and two nontypeable isolates. For a further genetic characterization, the genomes of all strains were sequenced using the 454 technology. Almost-complete genome sequences were assembled, and most remaining gaps were closed with Sanger sequencing. Phylogenetic analysis of 42 genomes revealed five established and two novel VZV genotypes, provisionally termed VIII and IX. Genotypes VIII and IX are distinct from the previously reported provisional genotypes VI and VII as judged from the SNP pattern. The alignments showed evidence of ancient recombination events in the phylogeny of clade 4 and recent recombinations within single strains: 3/2005 (clade 1), 11 and 405/2007 (clade 3), 8 and DR (clade 4), CA123 and 413/2000 (clade 5), and strains of the novel genotypes VIII and IX. Bayesian tree inference of the thymidine kinase and the polymerase genes of the VZV clades and other varicelloviruses revealed that VZV radiation began some 110,000 years ago, which correlates with the out-of-Africa dispersal of modern humans. The split of ancestral clades 2/4 and 1/3/5/VIII/IX shows the greatest node height.

The genome as a life-history character: why rate of molecular evolution varies between mammal species

DNA sequences evolve at different rates in different species. This rate variation has been most closely examined in mammals, revealing a large number of characteristics that can shape the rate of molecular evolution. Many of these traits are part of the mammalian life-history continuum: species with small body size, rapid generation turnover, high fecundity and short lifespans tend to have faster rates of molecular evolution. In addition, rate of molecular evolution in mammals might be influenced by behaviour (such as mating system), ecological factors (such as range restriction) and evolutionary history (such as diversification rate). I discuss the evidence for these patterns of rate variation, and the possible explanations of these correlations. I also consider the impact of these systematic patterns of rate variation on the reliability of the molecular date estimates that have been used to suggest a Cretaceous radiation of modern mammals, before the final extinction of the dinosaurs.

The historical biogeography of Mammalia

Palaeobiogeographic reconstructions are underpinned by phylogenies, divergence times and ancestral area reconstructions, which together yield ancestral area chronograms that provide a basis for proposing and testing hypotheses of dispersal and vicariance. Methods for area coding include multi-state coding with a single character, binary coding with multiple characters and string coding. Ancestral reconstruction methods are divided into parsimony versus Bayesian/likelihood approaches. We compared nine methods for reconstructing ancestral areas for placental mammals. Ambiguous reconstructions were a problem for all methods. Important differences resulted from coding areas based on the geographical ranges of extant species versus the geographical provenance of the oldest fossil for each lineage. Africa and South America were reconstructed as the ancestral areas for Afrotheria and Xenarthra, respectively. Most methods reconstructed Eurasia as the ancestral area for Boreoeutheria, Euarchontoglires and Laurasiatheria. The coincidence of molecular dates for the separation of Afrotheria and Xenarthra at approximately 100 Ma with the plate tectonic sundering of Africa and South America hints at the importance of vicariance in the early history of Placentalia. Dispersal has also been important including the origins of Madagascar's endemic mammal fauna. Further studies will benefit from increased taxon sampling and the application of new ancestral area reconstruction methods.

Convergent evolution of a symbiotic duet: the case of the lichen genus Polychidium (Peltigerales, Ascomycota)

PREMISE OF THE STUDY

Thallus architecture has long been a powerful guide for classifying lichens and has often trumped photobiont association and ascomatal type, but the reliability of these characters to predict phylogenetic affinity has seldom been tested. The cyanolichen genus Polychidium unites species that have strikingly similar gross morphology but consort with different photobiont genera. If Polychidium were found to be monophyletic, photobiont switching among closely related species would be suggested. If, however, species were found to arise in different lineages, a convergent body plan and ascomatal type evolution would be inferred.

METHODS

We tested the monophyly of Polychidium with a multilocus phylogeny based on nuclear and mitochondrial sequence data from all known Peltigeralean families and reconstructed ancestral states for specific thallus architecture and ascomatal ontogeny types relative to Polychidium and other clades.

KEY RESULTS

We found that Polychidium consists of two species groups that arose independently in different suborders within the Peltigerales, associated with Nostoc and Scytonema photobionts, respectively. We infer from ancestral character state reconstruction that dendroid thallus architecture evolved independently in these two lineages.

CONCLUSIONS

The independent development of similar dendroid thallus architecture in different fungal suborders with different photobionts represents a clear and previously overlooked example of convergent evolution in lichens. Our results also suggest a pattern of character state conservation, loss, and reversion in ascomatal ontogeny types, hitherto considered conserved traits useful for higher level ascomycete systematics.

Potential pitfalls of modelling ribosomal RNA data in phylogenetic tree reconstruction: evidence from case studies in the Metazoa

Background

Failure to account for covariation patterns in helical regions of ribosomal RNA (rRNA) genes has the potential to misdirect the estimation of the phylogenetic signal of the data. Furthermore, the extremes of length variation among taxa, combined with regional substitution rate variation can mislead the alignment of rRNA sequences and thus distort subsequent tree reconstructions. However, recent developments in phylogenetic methodology now allow a comprehensive integration of secondary structures in alignment and tree reconstruction analyses based on rRNA sequences, which has been shown to correct some of these problems. Here, we explore the potentials of RNA substitution models and the interactions of specific model setups with the inherent pattern of covariation in rRNA stems and substitution rate variation among loop regions.

Results

We found an explicit impact of RNA substitution models on tree reconstruction analyses. The application of specific RNA models in tree reconstructions is hampered by interaction between the appropriate modelling of covarying sites in stem regions, and excessive homoplasy in some loop regions. RNA models often failed to recover reasonable trees when single-stranded regions are excessively homoplastic, because these regions contribute a greater proportion of the data when covarying sites are essentially downweighted. In this context, the RNA6A model outperformed all other models, including the more parametrized RNA7 and RNA16 models.

Conclusions

Our results depict a trade-off between increased accuracy in estimation of interdependencies in helical regions with the risk of magnifying positions lacking phylogenetic signal. We can therefore conclude that caution is warranted when applying rRNA covariation models, and suggest that loop regions be independently screened for phylogenetic signal, and eliminated when they are indistinguishable from random noise. In addition to covariation and homoplasy, other factors, like non-stationarity of substitution rates and base compositional heterogeneity, can disrupt the signal of ribosomal RNA data. All these factors dictate sophisticated estimation of evolutionary pattern in rRNA data, just as other molecular data require similarly complicated (but different) corrections.

Morphological cranial diversity contributes to phylogeny in soft-furred sengis (Afrotheria, Macroscelidea)

Despite the well-supported Macroscelidea phylogeny proposed at the end of the 1960s, several systematic arrangements have been suggested in the last 20 years, raising doubts about the phylogeny of the Macroscelidinae; sengi inter-specific relationships are still debated to this day. The main issue of concern involves the supposed Elephantulus diphyly. To solve this persisting debate about sengi phylogeny, we examined the cranium ventral surface of 13 species using geometric morphometric techniques and neighbour-joining algorithms. This study supported the idea that the ventral side of the sengi cranium has the potential to provide important signals for reconstructing the Macroscelidea phylogeny. The phylogenetic signals seemed to differentiate between two major clades in the sengi radiation. In the first clade, the two monospecific genera (Petrodromus and Macroscelides), the two African Horn species (Elephantulus revoilii and E. rufescens), and the only North African species (E. rozeti) were clustered together. The second clade includes the remnant south-central African Elephantulus species. Our results were in agreement with both mitochondrial and nuclear data, confirmed that there is no Elephantulus monophyly and highlighted the close relationship between Petrodromus and E. rozeti. It appears that all the soft-furred sengi species are organised in two evolutionary lines: an old monophyletic clade, comprising only Elephantulus species, and a new polyphyletic clade, including P. tetradactylus, M. proboscideus, and E. rozeti. This requires a taxonomic and nomenclatural rearrangement within Macroscelidinae, where the phylogenetic position of the remnant 4 (of 12) Elephantulus species has yet to be fully defined.

Site-specific time heterogeneity of the substitution process and its impact on phylogenetic inference

Background

Model violations constitute the major limitation in inferring accurate phylogenies. Characterizing properties of the data that are not being correctly handled by current models is therefore of prime importance. One of the properties of protein evolution is the variation of the relative rate of substitutions across sites and over time, the latter is the phenomenon called heterotachy. Its effect on phylogenetic inference has recently obtained considerable attention, which led to the development of new models of sequence evolution. However, thus far focus has been on the quantitative heterogeneity of the evolutionary process, thereby overlooking more qualitative variations.

Results

We studied the importance of variation of the site-specific amino-acid substitution process over time and its possible impact on phylogenetic inference. We used the CAT model to define an infinite mixture of substitution processes characterized by equilibrium frequencies over the twenty amino acids, a useful proxy for qualitatively estimating the evolutionary process. Using two large datasets, we show that qualitative changes in site-specific substitution properties over time occurred significantly. To test whether this unaccounted qualitative variation can lead to an erroneous phylogenetic tree, we analyzed a concatenation of mitochondrial proteins in which Cnidaria and Porifera were erroneously grouped. The progressive removal of the sites with the most heterogeneous CAT profiles across clades led to the recovery of the monophyly of Eumetazoa (Cnidaria+Bilateria), suggesting that this heterogeneity can negatively influence phylogenetic inference.

Conclusion

The time-heterogeneity of the amino-acid replacement process is therefore an important evolutionary aspect that should be incorporated in future models of sequence change.

Molecular characterisation of Toll-like receptors in the black flying fox Pteropus alecto

Bats are believed to be reservoir hosts for a number of emerging and re-emerging viruses, many of which are responsible for illness and mortality in humans, livestock and other animals. In other vertebrates, early responses to viral infection involve engagement of Toll-like receptors (TLRs), which induce changes in gene expression collectively leading to an "antiviral state". In this study we report the cloning and bioinformatic analysis of a complete set of TLRs from the black flying fox Pteropus alecto, and perform quantitative tissue expression analysis of the nucleic acid-sensing TLRs 3, 7, 8 and 9. Full-length mRNA transcripts from TLRs homologous to human TLRs 1-10 were sequenced, as well as a nearly intact TLR13 pseudogene that was spliced and polyadenylated. This prototype data can now be used to design functional studies of the bat innate immune system.

Conservation and evolution in and among SRF- and MEF2-type MADS domains and their binding sites

Serum response factor (SRF) and myocyte enhancer factor 2 (MEF2) represent two types of members of the MCM1, AGAMOUS, DEFICIENS, and SRF (MADS)-box transcription factor family present in animals and fungi. Each type has distinct biological functions, which are reflected by the distinct specificities of the proteins bound to their cognate DNA-binding sites and activated by their respective cofactors. However, little is known about the evolution of MADS domains and their DNA-binding sites. Here, we report on the conservation and evolution of the two types of MADS domains with their cognate DNA-binding sites by using phylogenetic analyses. First, there are great similarities between the two types of proteins with amino acid positions highly conserved, which are critical for binding to the DNA sequence and for the maintenance of the 3D structure. Second, in contrast to MEF2-type MADS domains, distinct conserved residues are present at some positions in SRF-type MADS domains, determining specificity and the configuration of the MADS domain bound to DNA sequences. Furthermore, the ancestor sequence of SRF- and MEF2-type MADS domains is more similar to MEF2-type MADS domains than to SRF-type MADS domains. In the case of DNA-binding sites, the MEF2 site has a T-rich core in one DNA sequence and an A-rich core in the reverse sequence as compared with the SRF site, no matter whether where either A or T is present in the two complementary sequences. In addition, comparing SRF sites in the human and the mouse genomes reveals that the evolution rate of CArG-boxes is faster in mouse than in human. Moreover, interestingly, a CArG-like sequence, which is probably functionless, could potentially mutate to a functional CArG-box that can be bound by SRF and vice versa. Together, these results significantly improve our knowledge on the conservation and evolution of the MADS domains and their binding sites to date and provide new insights to investigate the MADS family, which is not only on evolution of MADS factors but also on evolution of their binding sites and even on coevolution of MADS factors with their binding sites.

Molecular analysis of neocortical layer structure in the ferret

Molecular markers that distinguish specific layers of rodent neocortex are increasingly employed to study cortical development and the physiology of cortical circuits. The extent to which these markers represent general features of neocortical cell type identity across mammals, however, is unknown. To assess the conservation of layer markers more broadly, we isolated orthologs for 15 layer-enriched genes in the ferret, a carnivore with a large, gyrencephalic brain, and analyzed their patterns of neocortical gene expression. Our major findings are: 1) Many but not all layer markers tested show similar patterns of layer-specific gene expression between mouse and ferret cortex, supporting the view that layer-specific cell type identity is conserved at a molecular level across mammalian superorders; 2) Our panel of deep layer markers (ER81/ETV1, SULF2, PCP4, FEZF2/ZNF312, CACNA1H, KCNN2/SK2, SYT6, FOXP2, CTGF) provides molecular evidence that the specific stratifications of layers 5 and 6 into 5a, 5b, 6a, and 6b are also conserved between rodents and carnivores; 3) Variations in layer-specific gene expression are more pronounced across areas of ferret cortex than between homologous areas of mouse and ferret cortex; 4) This variation of area gene expression was clearest with the superficial layer markers studied (SERPINE2, MDGA1, CUX1, UNC5D, RORB/NR1F2, EAG2/KCNH5). Most dramatically, the layer 4 markers RORB and EAG2 disclosed a molecular sublamination to ferret visual cortex and demonstrated a molecular dissociation among the so-called agranular areas of the neocortex. Our findings establish molecular markers as a powerful complement to cytoarchitecture for neocortical layer and cell-type comparisons across mammals.

Divergence time estimates of mammals from molecular clocks and fossils: relevance of new fossil finds from India

This paper presents a brief review of recent advances in the classification of mammals at higher levels using fossils and molecular clocks. It also discusses latest fossil discoveries from the Cretaceous - Eocene (66-55 m.y.) rocks of India and their relevance to our current understanding of placental mammal origins and diversifications.

The new framework for understanding placental mammal evolution

An unprecedented level of confidence has recently crystallized around a new hypothesis of how living placental mammals share a pattern of common descent. The major groups are afrotheres (e.g., aardvarks, elephants), xenarthrans (e.g., anteaters, sloths), laurasiatheres (e.g., horses, shrews), and euarchontoglires (e.g., humans, rodents). Compared with previous hypotheses this tree is remarkably stable; however, some uncertainty persists about the location of the placental root, and (for example) the position of bats within laurasiatheres, of sea cows and aardvarks within afrotheres, and of dermopterans within euarchontoglires. A variety of names for sub-clades within the new placental mammal tree have been proposed, not all of which follow conventions regarding priority and stability. More importantly, the new phylogenetic framework enables the formulation of new hypotheses and testing thereof, for example regarding the possible developmental dichotomy that seems to distinguish members of the newly identified southern and northern radiations of living placental mammals.

Deep evolutionary roots of strepsirrhine primate labyrinthine morphology

The cavity system of the inner ear of mammals is a complex three-dimensional structure that houses the organs of equilibrium and hearing. Morphological variation of the inner ear across mammals reflects differences in locomotor behaviour and hearing performance, and the good preservation of this structure in many fossil specimens permits analogous inferences. However, it is less well known to what extent the morphology of the bony labyrinth conveys information about the evolutionary history of primate taxa. We studied this question in strepsirrhine primates with the aim to assess the potential and limitations of using the inner ear as a phylogenetic marker. Geometric morphometric analysis showed that the labyrinthine morphology of extant strepsirrhines contains a mixed locomotor, allometric and phylogenetic signal. Discriminant analysis at the family level confirmed that labyrinthine shape is a good taxonomic marker. Our results support the hypothesis that evolutionary change in labyrinthine morphology is adequately described with a random walk model, i.e. random phenotypic dispersal in morphospace. Under this hypothesis, average shapes calculated for each node of the phylogenetic tree give an estimate of inner ear shapes of the respective last common ancestors (LCAs), and this information can be used to infer character state polarity. The labyrinthine morphology of the fossil Adapinae is close to the inferred basal morphology of the strepsirrhines. The inner ear of Daubentonia, one of the most derived extant strepsirrhines, is autapomorphic in many respects, but also presents unique similarities with adapine labyrinths.

Phylogenetic evidence for early hemochorial placentation in eutheria

The eutherian placenta is remarkable for its structural and functional variability. In order to construct and test comparative hypotheses relating ecological, behavioral and physiological traits to placental characteristics it is first necessary to reconstruct the historical course of placental evolution. Previous attempts to do so have yielded inconsistent results, particularly with respect to the early evolution of structural relationships between fetal and maternal circulatory systems. Here, we bring a battery of phylogenetic methods - including parsimony, likelihood and Bayesian approaches - to bear on the question of placental evolution. All of these approaches are consistent in indicating that highly invasive hemochorial placentation, as found in human beings and numerous other taxa, was an early evolutionary innovation present in the most ancient ancestors of the living placental mammals.

Whole-genome phylogeny of mammals: evolutionary information in genic and nongenic regions

Ten complete mammalian genome sequences were compared by using the "feature frequency profile" (FFP) method of alignment-free comparison. This comparison technique reveals that the whole nongenic portion of mammalian genomes contains evolutionary information that is similar to their genic counterparts--the intron and exon regions. We partitioned the complete genomes of mammals (such as human, chimp, horse, and mouse) into their constituent nongenic, intronic, and exonic components. Phylogenic species trees were constructed for each individual component class of genome sequence data as well as the whole genomes by using standard tree-building algorithms with FFP distances. The phylogenies of the whole genomes and each of the component classes (exonic, intronic, and nongenic regions) have similar topologies, within the optimal feature length range, and all agree well with the evolutionary phylogeny based on a recent large dataset, multispecies, and multigene-based alignment. In the strictest sense, the FFP-based trees are genome phylogenies, not species phylogenies. However, the species phylogeny is highly related to the whole-genome phylogeny. Furthermore, our results reveal that the footprints of evolutionary history are spread throughout the entire length of the whole genome of an organism and are not limited to genes, introns, or short, highly conserved, nongenic sequences that can be adversely affected by factors (such as a choice of sequences, homoplasy, and different mutation rates) resulting in inconsistent species phylogenies.

Genomic imprinting in Canis familiaris

For the vast majority of mammalian genes, maternally- and paternally-derived alleles behave identically and are either expressed or repressed, regardless of whether they were inherited from egg or sperm. For imprinted genes, however, this is not the case. The alleles of imprinted genes are epigenetically modified in a parent-of-origin-specific manner and, as a consequence, maternally- and paternally-derived alleles behave differently. Typically one allele is expressed while the other is silent. Although relatively few in number, imprinted genes are the focus of intensive study, as they have important roles in embryonic development. Abnormal expression of imprinted genes results in growth disorders and is implicated in several clinical conditions. Most studies of imprinted genes have been performed in rodents or primates, with limited studies in other mammals such as bovine and opossum. We have recently demonstrated the existence of imprinted genes in the canine, by showing that the canine insulin-like growth factor 2 receptor gene (IGF2R) is monoallelically expressed, with predominant expression of the maternally-derived allele and repression of the paternally-inherited allele. Our ultimate goal is to characterize all imprinted genes in the canine, and to understand how they contribute to canine reproduction, development and disease. Such knowledge will be vital for optimizing the success of most reproductive strategies in the canine.

Assessing the genomic evidence for conserved transcribed pseudogenes under selection

BACKGROUND

Transcribed pseudogenes are copies of protein-coding genes that have accumulated indicators of coding-sequence decay (such as frameshifts and premature stop codons), but nonetheless remain transcribed. Recent experimental evidence indicates that transcribed pseudogenes may regulate the expression of homologous genes, through antisense interference, or generation of small interfering RNAs (siRNAs). Here, we assessed the genomic evidence for such transcribed pseudogenes of potential functional importance, in the human genome. The most obvious indicators of such functional importance are significant evidence of conservation and selection pressure.

RESULTS

A variety of pseudogene annotations from multiple sources were pooled and filtered to obtain a subset of sequences that have significant mid-sequence disablements (frameshifts and premature stop codons), and that have clear evidence of full-length mRNA transcription. We found 1750 such transcribed pseudogene annotations (TPAs) in the human genome (corresponding to approximately 11.5% of human pseudogene annotations). We checked for syntenic conservation of TPAs in other mammals (rhesus monkey, mouse, rat, dog and cow). About half of the human TPAs are conserved in rhesus monkey, but strikingly, very few in mouse (approximately 3%). The TPAs conserved in rhesus monkey show evidence of selection pressure (relative to surrounding intergenic DNA) on: (i) their GC content, and (ii) their rate of nucleotide substitution. This is in spite of distributions of Ka/Ks (ratios of non-synonymous to synonymous substitution rates), congruent with a lack of protein-coding ability. Furthermore, we have identified 68 human TPAs that are syntenically conserved in at least two other mammals. Interestingly, we observe three TPA sequences conserved in dog that have intermediate character (i.e., evidence of both protein-coding ability and pseudogenicity), and discuss the implications of this.

CONCLUSION

Through evolutionary analysis, we have identified candidate sequences for functional human transcribed pseudogenes, and have pinpointed 68 strong candidates for further investigation as potentially functional transcribed pseudogenes across multiple mammal species.

Paleocene emergence of elephant relatives and the rapid radiation of African ungulates

Elephants are the only living representatives of the Proboscidea, a formerly diverse mammalian order whose history began with the 55-million years (mys) old Phosphatherium. Reported here is the discovery from the early late Paleocene of Morocco, ca. 60 mys, of the oldest and most primitive elephant relative, Eritherium azzouzorum n.g., n.sp., which is one of the earliest known representatives of modern placental orders. This well supported stem proboscidean is extraordinarily primitive and condylarth-like. It provides the first dental evidence of a resemblance between the proboscideans and African ungulates (paenungulates) on the one hand and the louisinines and early macroscelideans on the other. Eritherium illustrates the origin of the elephant order at a previously unknown primitive stage among paenungulates and "ungulates." The primitive morphology of Eritherium suggests a recent and rapid paenungulate radiation after the Cretaceous-Tertiary boundary, probably favoured by early endemic African paleoecosystems. At a broader scale, Eritherium provides a new old calibration point of the placental tree and supports an explosive placental radiation. The Ouled Abdoun basin, which yields the oldest known African placentals, is a key locality for elucidating phylogeny and early evolution of paenungulates and other related endemic African lineages.

Breakpoint graphs and ancestral genome reconstructions

Recently completed whole-genome sequencing projects marked the transition from gene-based phylogenetic studies to phylogenomics analysis of entire genomes. We developed an algorithm MGRA for reconstructing ancestral genomes and used it to study the rearrangement history of seven mammalian genomes: human, chimpanzee, macaque, mouse, rat, dog, and opossum. MGRA relies on the notion of the multiple breakpoint graphs to overcome some limitations of the existing approaches to ancestral genome reconstructions. MGRA also generates the rearrangement-based characters guiding the phylogenetic tree reconstruction when the phylogeny is unknown.