Molecules remodel the mammalian tree

Characterization of Common Minke Whale (Balaenoptera Acutorostrata) Cell Lines Immortalized with the Expression of Cell Cycle Regulators

Primary cultured cells cannot proliferate infinite. The overcoming of this limit can be classified as immortalization. Bypass of p16 senescence protein induces efficient immortalization various types of mammalians is previously reported. However, the Cetacea species is not known. Here, that common minke whale-derived cells can be immortalized with a combination of human genes, mutant cyclin-dependent kinase 4 (CDK4R24C ), cyclin D1, and Telomerase Reverse Transcriptase (TERT) is reported. These results indicate that the function of cell cycle regulators in premature senescence is evolutionarily conserved. This study describes the conserved roles of cell cycle regulators in the immortalization of cells from humans to Cetacea species. Furthermore, using RNA-seq based on next-generation sequencing, the gene expression profiles of immortalized cells are compared with parental cells as well as those immortalized with SV40 large T antigen, which is once a popular method for cellular immortalization. The profiling results show that newly established common minke-whale-derived immortaliozed cells have completely different profiles from SV40 cells. This result indicates that the expression of mutant CDK4, cyclin D1, and TERT enables to establish immortalized cell lines with different biological nature from SV40 expressing cells.

Aging: What We Can Learn From Elephants

Elephants are large-brained, social mammals with a long lifespan. Studies of elephants can provide insight into the aging process, which may be relevant to understanding diseases that affect elderly humans because of their shared characteristics that have arisen through independent evolution. Elephants become sexually mature at 12 to 14 years of age and are known to live into, and past, their 7^th decade of life. Because of their relatively long lifespans, elephants may have evolved mechanisms to counter age-associated morbidities, such as cancer and cognitive decline. Elephants rely heavily on their memory, and engage in multiple levels of competitive and collaborative relationships because they live in a fission-fusion system. Female matrilineal relatives and dependent offspring form tight family units led by an older-aged matriarch, who serves as the primary repository for social and ecological knowledge in the herd. Similar to humans, elephants demonstrate a dependence on social bonds, memory, and cognition to navigate their environment, behaviors that might be associated with specializations of brain anatomy. Compared with other mammals, the elephant hippocampus is proportionally smaller, whereas the temporal lobe is disproportionately large and expands laterally. The elephant cerebellum is also relatively enlarged, and the cerebral cortex is highly convoluted with numerous gyral folds, more than in humans. Last, an interesting characteristic unique to elephants is the presence of at least 20 copies of the TP53 tumor suppressor gene. Humans have only a single copy. TP53 encodes for the p53 protein, which is known to orchestrate cellular response to DNA damage. The effects of these multiple copies of TP53 are still being investigated, but it may be to protect elephants against multiple age-related diseases. For these reasons, among others, studies of elephants would be highly informative for aging research. Elephants present an underappreciated opportunity to explore further common principles of aging in a large-brained mammal with extended longevity. Such research can contribute to contextualizing our knowledge of age-associated morbidities in humans.

Elaboration and Innervation of the Vibrissal System in the Rock Hyrax (Procavia capensis)

Mammalian tactile hairs are commonly found on specific, restricted regions of the body, but Florida manatees represent a unique exception, exhibiting follicle-sinus complexes (FSCs, also known as vibrissae or tactile hairs) on their entire body. The orders Sirenia (including manatees and dugongs) and Hyracoidea (hyraxes) are thought to have diverged approximately 60 million years ago, yet hyraxes are among the closest relatives to sirenians. We investigated the possibility that hyraxes, like manatees, are tactile specialists with vibrissae that cover the entire postfacial body. Previous studies suggested that rock hyraxes possess postfacial vibrissae in addition to pelage hair, but this observation was not verified through histological examination. Using a detailed immunohistochemical analysis, we characterized the gross morphology, innervation and mechanoreceptors present in FSCs sampled from facial and postfacial vibrissae body regions to determine that the long postfacial hairs on the hyrax body are in fact true vibrissae. The types and relative densities of mechanoreceptors associated with each FSC also appeared to be relatively consistent between facial and postfacial FSCs. The presence of vibrissae covering the hyrax body presumably facilitates navigation in the dark caves and rocky crevices of the hyrax's environment where visual cues are limited, and may alert the animal to predatory or conspecific threats approaching the body. Furthermore, the presence of vibrissae on the postfacial body in both manatees and hyraxes indicates that this distribution may represent the ancestral condition for the supraorder Paenungulata.

Feature selection methods for identifying genetic determinants of host species in RNA viruses

Despite environmental, social and ecological dependencies, emergence of zoonotic viruses in human populations is clearly also affected by genetic factors which determine cross-species transmission potential. RNA viruses pose an interesting case study given their mutation rates are orders of magnitude higher than any other pathogen – as reflected by the recent emergence of SARS and Influenza for example. Here, we show how feature selection techniques can be used to reliably classify viral sequences by host species, and to identify the crucial minority of host-specific sites in pathogen genomic data. The variability in alleles at those sites can be translated into prediction probabilities that a particular pathogen isolate is adapted to a given host. We illustrate the power of these methods by: 1) identifying the sites explaining SARS coronavirus differences between human, bat and palm civet samples; 2) showing how cross species jumps of rabies virus among bat populations can be readily identified; and 3) de novo identification of likely functional influenza host discriminant markers.

Moving away from genome scan methods used for human GWAS (ultimately inappropriate for the short highly polymorphic genomes of RNA viruses), our work shows the power and potential of multi-class machine learning algorithms in inferring the functional genetic changes associated with phenotypic change (e.g. crossing a species barrier). We show that even distantly related viruses within a viral family share highly conserved genetic signatures of host specificity; reinforce how fitness landscapes of host adaptation are shaped by host phylogeny; and highlight the evolutionary trajectories of RNA viruses in rapid expansion and under great evolutionary pressure. We do so by (for each dataset) unveiling a set of phenotype characteristic mutations which are shown to be functionally relevant, thus providing new insights into phenotypic relationships between RNA viruses. These methods also provide a solid statistical framework with which the degree of host adaptation can be inferred, thus serving as a valuable tool for studying host transition events with particular relevance for emerging infectious diseases. These methods can then serve as rigorous tools of emergence potential assessment, specifically in scenarios where rapid host classification of newly emerging viruses can be more important than identifying putative functional sites.

Cellular scaling rules of insectivore brains

Insectivores represent extremes in mammalian body size and brain size, retaining various “primitive” morphological characteristics, and some species of Insectivora are thought to share similarities with small-bodied ancestral eutherians. This raises the possibility that insectivore brains differ from other taxa, including rodents and primates, in cellular scaling properties. Here we examine the cellular scaling rules for insectivore brains and demonstrate that insectivore scaling rules overlap somewhat with those for rodents and primates such that the insectivore cortex shares scaling rules with rodents (increasing faster in size than in numbers of neurons), but the insectivore cerebellum shares scaling rules with primates (increasing isometrically). Brain structures pooled as “remaining areas” appear to scale similarly across all three mammalian orders with respect to numbers of neurons, and the numbers of non-neurons appear to scale similarly across all brain structures for all three orders. Therefore, common scaling rules exist, to different extents, between insectivore, rodent, and primate brain regions, and it is hypothesized that insectivores represent the common aspects of each order. The olfactory bulbs of insectivores, however, offer a noteworthy exception in that neuronal density increases linearly with increasing structure mass. This implies that the average neuronal cell size decreases with increasing olfactory bulb mass in order to accommodate greater neuronal density, and represents the first documentation of a brain structure gaining neurons at a greater rate than mass. This might allow insectivore brains to concentrate more neurons within the olfactory bulbs without a prohibitively large and metabolically costly increase in structure mass.

Large brains and cognition: Where do elephants fit in?

Among terrestrial mammals, elephants share the unique status, along with humans and great apes, of having large brains, being long-lived and having offspring that require long periods of dependency. Elephants have the largest brains of all terrestrial mammals, including the greatest volume of cerebral cortex. In contrast to what one might expect from such a large-brained species, the performance of elephants in cognitive feats, such as tool use, visual discrimination learning and tests of "insight" behavior, is unimpressive in comparison to the performance by chimpanzees and, of course, humans. Where elephants do seem to excel is in long-term, extensive spatial-temporal and social memory. In addition, elephants appear to be somewhat unique among non-human species in their reactions to disabled and deceased conspecifics, exhibiting behaviors that are mindful of "theory-of-mind" phenomena. Information gleaned from studies on the neural cytoarchitecture of large brains reveals that the neurons of the cerebral cortex of elephants are much less densely populated than in large-brained primates. The interactions between cortical neurons would appear to be more global and less compartmentalized into local areas, and cortical information processing slower, than in great apes and humans. Although focused neural cytoarchitecture studies on the elephant are needed, this comparative perspective on the cortical neural cytoarchitecture appears to relate to differences in behavior between elephants and their primate counterparts.

Site specific rates of mitochondrial genomes and the phylogeny of eutheria

BACKGROUND

Traditionally, most studies employing data from whole mitochondrial genomes to diagnose relationships among the major lineages of mammals have attempted to exclude regions that potentially complicate phylogenetic analysis. Components generally excluded are 3rd codon positions of protein-encoding genes, the control region, rRNAs, tRNAs, and the ND6 gene (encoded on the opposite strand). We present an approach that includes all the data, with the exception of the control region. This approach is based on a site-specific rate model that accommodates excessive homoplasy and that utilizes secondary structure as a reference for proper alignment of rRNAs and tRNAs.

RESULTS

Mitochondrial genomic data for 78 eutherian mammals, 8 metatherians, and 3 monotremes were analyzed with a Bayesian analysis and our site specific rate model. The resultant phylogeny revealed strong support for most nodes and was highly congruent with more recent phylogenies based on nuclear DNA sequences. In addition, many of the conflicting relationships observed by earlier mitochondrial-based analyses were resolved without need for the exclusion of large subsets of the data.

CONCLUSION

Rather than exclusion of data to minimize presumed noise associated with non-protein encoding genes in the mitochondrial genome, our results indicate that selection of an appropriate model that accommodates rate heterogeneity across data partitions and proper treatment of RNA genes can result in a mitochondrial genome-based phylogeny of eutherian mammals that is reasonably congruent with recent phylogenies derived from nuclear genes.

The evolutionary radiation of Arvicolinae rodents (voles and lemmings): relative contribution of nuclear and mitochondrial DNA phylogenies

BACKGROUND

Mitochondrial and nuclear genes have generally been employed for different purposes in molecular systematics, the former to resolve relationships within recently evolved groups and the latter to investigate phylogenies at a deeper level. In the case of rapid and recent evolutionary radiations, mitochondrial genes like cytochrome b (CYB) are often inefficient for resolving phylogenetic relationships. One of the best examples is illustrated by Arvicolinae rodents (Rodentia; Muridae), the most impressive mammalian radiation of the Northern Hemisphere which produced voles, lemmings and muskrats. Here, we compare the relative contribution of a nuclear marker--the exon 10 of the growth hormone receptor (GHR) gene--to the one of the mitochondrial CYB for inferring phylogenetic relationships among the major lineages of arvicoline rodents.

RESULTS

The analysis of GHR sequences improves the overall resolution of the Arvicolinae phylogeny. Our results show that the Caucasian long-clawed vole (Prometheomys schaposnikowi) is one of the basalmost arvicolines, and confirm that true lemmings (Lemmus) and collared lemmings (Dicrostonyx) are not closely related as suggested by morphology. Red-backed voles (Myodini) are found as the sister-group of a clade encompassing water vole (Arvicola), snow vole (Chionomys), and meadow voles (Microtus and allies). Within the latter, no support is recovered for the generic recognition of Blanfordimys, Lasiopodomys, Neodon, and Phaiomys as suggested by morphology. Comparisons of parameter estimates for branch lengths, base composition, among sites rate heterogeneity, and GTR relative substitution rates indicate that CYB sequences consistently exhibit more heterogeneity among codon positions than GHR. By analyzing the contribution of each codon position to node resolution, we show that the apparent higher efficiency of GHR is due to their third positions. Although we focus on speciation events spanning the last 10 million years (Myr), CYB sequences display highly saturated codon positions contrary to the nuclear exon. Lastly, variable length bootstrap predicts a significant increase in resolution of arvicoline phylogeny through the sequencing of nuclear data in an order of magnitude three to five times greater than the size of GHR exon 10.

CONCLUSION

Our survey provides a first resolved gene tree for Arvicolinae. The comparison of CYB and GHR phylogenetic efficiency supports recent assertions that nuclear genes are useful for resolving relationships of recently evolved animals. The superiority of nuclear exons may reside both in (i) less heterogeneity among sites, and (ii) the presence of highly informative sites in third codon positions, that evolve rapidly enough to accumulate synapomorphies, but slow enough to avoid substitutional saturation.

The mid-domain effect revisited

We revisit the proposition that boundary constraints on species' ranges cause species richness gradients (the mid-domain effect [MDE] hypothesis). In the absence of environmental gradients, species should not retain their observed range sizes as assumed by MDE models. Debate remains regarding the definition of domain limits, valid predictions for testing the models, and their statistical assessment. Empirical support for the MDE is varied but often weak, suggesting that geometric constraints on species' ranges do not provide a general explanation for richness gradients. Criticism of MDE model assumptions does not, however, imply opposition to the use of null models in ecology.

Visualizing conflicting evolutionary hypotheses in large collections of trees: using consensus networks to study the origins of placentals and hexapods

Many phylogenetic methods produce large collections of trees as opposed to a single tree, which allows the exploration of support for various evolutionary hypotheses. However, to be useful, the information contained in large collections of trees should be summarized; frequently this is achieved by constructing a consensus tree. Consensus trees display only those signals that are present in a large proportion of the trees. However, by their very nature consensus trees require that any conflicts between the trees are necessarily disregarded. We present a method that extends the notion of consensus trees to allow the visualization of conflicting hypotheses in a consensus network. We demonstrate the utility of this method in highlighting differences amongst maximum likelihood bootstrap values and Bayesian posterior probabilities in the placental mammal phylogeny, and also in comparing the phylogenetic signal contained in amino acid versus nucleotide characters for hexapod monophyly.

Primate phylogeny, evolutionary rate variations, and divergence times: a contribution from the nuclear gene IRBP

The first third (ca. 1200 bp) of exon 1 of the nuclear gene encoding the interstitial retinoid-binding protein (IRBP) has been sequenced for 12 representative primates belonging to Lemuriformes, Lorisiformes, Tarsiiformes, Platyrrhini, and Catarrhini, and combined with available data (13 other primates, 11 nonprimate placentals, and 2 marsupials). Phylogenetic analyses using maximum likelihood on nucleotides and amino acids robustly support the monophyly of primates, Strepsirrhini, Lemuriformes, Lorisiformes, Anthropoidea, Catarrhini, and Platyrrhini. It is interesting to note that 1) Tarsiidae grouped with Anthropoidea, and the support for this node depends on the molecular characters considered; 2) Cheirogaleidae grouped within Lemuriformes; and 3) Daubentonia was the sister group of all other Lemuriformes. Study of the IRBP evolutionary rate shows a high heterogeneity within placentals and also within primates. Maximum likelihood local molecular clocks were assigned to three clades displaying significantly contrasted evolutionary rates. Paenungulata were shown to evolve 2.5-3 times faster than Perissodactyla and Lemuriformes. Six independent calibration points were used to estimate splitting ages of the main primate clades, and their compatibility was evaluated. Divergence ages were obtained for the following crown groups: 13.8-14.2 MY for Lorisiformes, 26.5-27.2 MY for Lemuroidea, 39.6-40.7 MY for Lemuriformes, 45.4-46.7 MY for Strepsirrhini, and 56.7-58.4 MY for Haplorrhini. The incompatibility between some paleontological and molecular estimates may reflect the incompleteness of the placental fossil record, and/or indicate that the variable IRBP evolutionary rates are not fully accommodated by local molecular clocks.

Characterization of a cetacean aromatase (CYP19) and the phylogeny and functional conservation of vertebrate aromatase

Aromatase (cytochrome P450 19, CYP19, P450arom) is the enzyme responsible for the production of estrogens, hormones critical for development and reproduction. Aromatase was sequenced from a white-sided dolphin (Lagenorhynchus acutus) ovary, transiently transfected into HEK 293 cells, and the expressed protein was characterized for aromatase activity in the presence of androstenedione and testosterone and after exposure to the aromatase inhibitor letrazole. The Kms for androstenedione and testosterone were 63.5 and 75 nM, respectively, values that are very similar to those reported for other mammalian aromatases. A Bayesian phylogenetic analysis of the vertebrate aromatases was performed on the amino acid sequences of aromatases from fish, amphibians, reptiles, birds, and mammals. Based on known species phylogeny, the cetacean aromatase showed an expected grouping with artiodactyls (cow, sheep, and goat). An analysis of functional divergence showed strong conservation of aromatase across the entire protein, which indicates that the observed sequence divergence is functionally neutral.

Fetal lung development in the elephant reflects the adaptations required for snorkeling in adult life

The adult elephant is unique among mammals in that the pleural membranes are thickened and the pleural cavity is obliterated by connective tissue. It has been suggested that this peculiar anatomy developed because the animal can snorkel at depth, and this behavior subjects the microvessels in the parietal pleura to a very large transmural pressure. To investigate the development of the parietal pleura, the thickness of the endothoracic fascia (ET) was measured in four fetal African elephants of approximate gestational age 111-130 days, and the appearances were compared with those in human, rabbit, rat and mouse fetuses of approximately the same stage of lung organogenesis. The mean thicknesses of ET in the elephant, human, rabbit, rat and mouse were 403, 53, 29, 27 and 37 microm, respectively. This very early development of a thick parietal pleura in the elephant fetus is consistent with the hypothesis of a long history of snorkeling in the elephant's putative aquatic ancestors.

Evaluating placental inter-ordinal phylogenies with novel sequences including RAG1, gamma-fibrinogen, ND6, and mt-tRNA, plus MCMC-driven nucleotide, amino acid, and codon models

It is essential to test a priori scientific hypotheses with independent data, not least to partly negate factors such as gene-specific base composition biases misleading our models. Seven new gene segments and sequences plus Bayesian likelihood phylogenetic methods were used to compare and test five recent placental phylogenies. These five phylogenies are similar to each other, yet quite different from Fthose of previously proposed trees, and span Waddell et al. [Syst. Biol. 48 (1999) 1] to Murphy et al. [Science 294 (2001b) 2348]. Trees for RAG1, gamma-fibrinogen, ND6, mt-tRNA, mt-RNA, c-MYC, epsilon -globin, and GHR are significantly congruent with the four main groups of mammals common to the five phylogenies, i.e., Afrotheria, Laurasiatheria, Euarchontoglires, Xenarthra plus Boreoeutheria (Laurasiatheria plus Euarchontoglires). Where these five a priori phylogenies differ, remain areas generally hard to resolve with the new sequences. The root remains ambiguous and does not reject a basal Afrotheria (the Exafroplacentalia hypothesis), Afrotheria plus Xenarthra together with basal (Atlantogenata), or Epitheria (Xenarthra basal) convincingly. Good evidence is found that Eulipotyphla is monophyletic and is located at the base of Laurasiatheria. The shrew mole, Uropsilus, is found to cluster consistently with other moles, while Solenodon may be the sister taxa to all other eulipotyphlans. Support is found for a probable sister pairing of just hedgehogs/gymnures and shrews. Relationships within Afrotheria, except the Paenungulata clade, remain hard to resolve, although there is congruent support for Afroinsectiphillia (aardvark, elephant shrews, golden moles, and tenrecs). A first-time use is made of MCMC enacted general time-reversible (GTR) amino acid and codon-based models for general tree selection. Even with ND6, a GTR amino acid model provided resolution of fine features, such as the sister group relationship of walrus to Otatriidae, and with BRCA a more reasonable rooting. An extensive analysis of GHR sequences reveals strong congruence with prior phylogenies, including strong support for Eulipotyphla, and good resolution within Rodentia. A codon model gives a worse likelihood than a nucleotide model and sometimes switches support, e.g., with RAG1+gamma-fibrinogen from a hyrax-sirenian association to support for Tethytheria. An analysis of the concatenated data is in accordance with well-resolved features of the gene trees. Taken all together, this work suggests that we are on the right path finding strong confirmation of prior phylogenies. However, with the use of robust criteria for assessing trees (i.e., not Bayesian posteriors), it is apparent parts of the tree remain hard to resolve. Since our current models are far from fitting the sequence data, we should continue with our exploratory analyses to arrive at a refined set of hypotheses for future testing using more model independent characters (e.g., rare indels, gene rearrangement, and SINE data).

Recreating a functional ancestral archosaur visual pigment

The ancestors of the archosaurs, a major branch of the diapsid reptiles, originated more than 240 MYA near the dawn of the Triassic Period. We used maximum likelihood phylogenetic ancestral reconstruction methods and explored different models of evolution for inferring the amino acid sequence of a putative ancestral archosaur visual pigment. Three different types of maximum likelihood models were used: nucleotide-based, amino acid-based, and codon-based models. Where possible, within each type of model, likelihood ratio tests were used to determine which model best fit the data. Ancestral reconstructions of the ancestral archosaur node using the best-fitting models of each type were found to be in agreement, except for three amino acid residues at which one reconstruction differed from the other two. To determine if these ancestral pigments would be functionally active, the corresponding genes were chemically synthesized and then expressed in a mammalian cell line in tissue culture. The expressed artificial genes were all found to bind to 11-cis-retinal to yield stable photoactive pigments with lambda(max) values of about 508 nm, which is slightly redshifted relative to that of extant vertebrate pigments. The ancestral archosaur pigments also activated the retinal G protein transducin, as measured in a fluorescence assay. Our results show that ancestral genes from ancient organisms can be reconstructed de novo and tested for function using a combination of phylogenetic and biochemical methods.

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.

The evolution of armadillos, anteaters and sloths depicted by nuclear and mitochondrial phylogenies: implications for the status of the enigmatic fossil Eurotamandua

The mammalian order Xenarthra (armadillos, anteaters and sloths) is one of the four major clades of placentals, but it remains poorly studied from the molecular phylogenetics perspective. We present here a study encompassing most of the order's diversity in order to establish xenarthrans' intra-ordinal relationships, discuss the evolution of their morphological characters, search for their extant sister group and specify the timing of their radiation with special emphasis on the status of the controversial fossil Eurotamandua. Sequences of three genes (nuclear exon 28 of the Von Willebrand factor and mitochondrial 12S and 16S rRNAs) are compared for eight of the 13 living genera. Phylogenetic analyses confirm the order's monophyly and that of its three major lineages: armadillos (Cingulata), anteaters (Vermilingua) and sloths ('Tardigrada', renamed in 'Folivora'), and our results strongly support the grouping of hairy xenarthrans (anteaters and sloths) into Pilosa. Within placentals, Afrotheria might be the first lineage to branch off, followed by Xenarthra. The morphological adaptative convergence between New World xenarthrans and Old World pangolins is confirmed. Molecular datings place the early emergence of armadillos around the Cretaceous/Tertiary boundary, followed by the divergence between anteaters and sloths in the Early Eocene era. These Tertiary dates contradict the concept of a very ancient origin of modern xenarthran lineages. They also question the placement of the purported fossil anteater (Eurotamandua) from the Middle Eocene period of Europe with the Vermilingua and instead suggest the independent and convergent evolution of this enigmatic taxon.

The molecular genetics and evolution of red and green color vision in vertebrates

To better understand the evolution of red-green color vision in vertebrates, we inferred the amino acid sequences of the ancestral pigments of 11 selected visual pigments: the LWS pigments of cave fish (Astyanax fasciatus), frog (Xenopus laevis), chicken (Gallus gallus), chameleon (Anolis carolinensis), goat (Capra hircus), and human (Homo sapiens);and the MWS pigments of cave fish, gecko (Gekko gekko), mouse (Mus musculus), squirrel (Sciurus carolinensis), and human. We constructed these ancestral pigments by introducing the necessary mutations into contemporary pigments and evaluated their absorption spectra using an in vitro assay. The results show that the common ancestor of vertebrates and most other ancestors had LWS pigments. Multiple regression analyses of ancestral and contemporary MWS and LWS pigments show that single mutations S180A, H197Y, Y277F, T285A, A308S, and double mutations S180A/H197Y shift the lambda(max) of the pigments by -7, -28, -8, -15, -27, and 11 nm, respectively. It is most likely that this "five-sites" rule is the molecular basis of spectral tuning in the MWS and LWS pigments during vertebrate evolution.

Snorkel breathing in the elephant explains the unique anatomy of its pleura

It has been known for over 300 years that the anatomy of the elephant lung is unique among mammals in that the pleural cavity is obliterated by connective tissue. However no satisfactory explanation has been advanced. Recent studies suggest that the elephant has an aquatic ancestry and the trunk may have developed for snorkeling. In addition, the modern day elephant is the only mammal that can remain submerged far below the surface of the water while snorkeling. The resulting differences of pressures within the thorax mean that the small blood vessels of the pleura are in great danger of rupturing or causing severe edema. The same distribution of pressures occurs when the animal raises water inside its trunk prior to drinking although in this case the pressure differences are relatively short-lived. Evolution has provided a remarkable solution to this problem by replacing the normally delicate parietal and visceral pleurae by dense connective tissue, and separating the two pleurae by loose connective tissue to allow some sliding movement.

Molecular and morphological supertrees for eutherian (placental) mammals

A large body of diverse comparative data now exists for a major phylogenetic synthesis of the higher-level relationships among eutherian (placental) mammals. We present such a phylogenetic synthesis using the composite trees or supertrees from the combined and separate analyses of their published molecular and morphological source phylogenies. Our combined and separate supertrees largely support the same suprafamilial taxa and orders, but different interordinal clades. These similarities and differences reinforce the continuing contributions of morphological studies, while highlighting the growing influence of molecular information on the field. As current summaries of past research, our supertrees emphasize opportunities for future work, while providing a step toward the eventual integration of the data and characters themselves.

Parallel adaptive radiations in two major clades of placental mammals

Higher level relationships among placental mammals, as well as the historical biogeography and morphological diversification of this group, remain unclear. Here we analyse independent molecular data sets, having aligned lengths of DNA of 5,708 and 2,947 base pairs, respectively, for all orders of placental mammals. Phylogenetic analyses resolve placental orders into four groups: Xenarthra, Afrotheria, Laurasiatheria, and Euarchonta plus Glires. The first three groups are consistently monophyletic with different methods of analysis. Euarchonta plus Glires is monophyletic or paraphyletic depending on the phylogenetic method. A unique nine-base-pair deletion in exon 11 of the BRCA1 gene provides additional support for the monophyly of Afrotheria, which includes proboscideans, sirenians, hyracoids, tubulidentates, macroscelideans, chrysochlorids and tenrecids. Laurasiatheria contains cetartiodactyls, perissodactyls, carnivores, pangolins, bats and eulipotyphlan insectivores. Parallel adaptive radiations have occurred within Laurasiatheria and Afrotheria. In each group, there are aquatic, ungulate and insectivore-like forms.

Mitochondrial versus nuclear gene sequences in deep-level mammalian phylogeny reconstruction

Both mitochondrial and nuclear gene sequences have been employed in efforts to reconstruct deep-level phylogenetic relationships. A fundamental question in molecular systematics concerns the efficacy of different types of sequences in recovering clades at different taxonomic levels. We compared the performance of four mitochondrial data sets (cytochrome b, cytochrome oxidase II, NADH dehydrogenase subunit I, 12S rRNA-tRNA-16S rRNA) and eight nuclear data sets (exonic regions of alpha-2B adrenergic receptor, aquaporin, ss-casein, gamma-fibrinogen, interphotoreceptor retinoid binding protein, kappa-casein, protamine, von Willebrand Factor) in recovering deep-level mammalian clades. We employed parsimony and minimum-evolution with a variety of distance corrections for superimposed substitutions. In 32 different pairwise comparisons between these mitochondrial and nuclear data sets, we used the maximum set of overlapping taxa. In each case, the variable-length bootstrap was used to resample at the size of the smaller data set. The nuclear exons consistently performed better than mitochondrial protein and rRNA-tRNA coding genes on a per-residue basis in recovering benchmark clades. We also concatenated nuclear genes for overlapping taxa and made comparisons with concatenated mitochondrial protein-coding genes from complete mitochondrial genomes. The variable-length bootstrap was used to score the recovery of benchmark clades as a function of the number of resampled base pairs. In every case, the nuclear concatenations were more efficient than the mitochondrial concatenations in recovering benchmark clades. Among genes included in our study, the nuclear genes were much less affected by superimposed substitutions. Nuclear genes having appropriate rates of substitution should receive strong consideration in efforts to reconstruct deep-level phylogenetic relationships.

An efficient cis-element discovery method using multiple sequence comparisons based on evolutionary relationships

The discovery of cis-element control motifs in noncoding DNA poses a difficult problem in genome analysis. Functional analysis by means of reporter constructs expressed in transgenic organisms is the most reliable method, but is by itself time-consuming and expensive. Searching noncoding DNA for known control motifs by sequence analysis is problematic, since protein binding motifs are short, in the range of 8-10 bp, and occur frequently by chance. Heretofore, the most reliable sequence analysis method has been the comparison of homologous sequence domains in related but moderately evolutionarily divergent species such as, for example, mouse and human. In such pairwise combinations, control regions are conserved because they serve a vital function and can be identified by their similar sequences. Single pairwise comparisons, however, allow the discovery of conserved sequence strings only at low resolution and without specific identity. We have investigated the possibility of using multiple sequence comparisons to correct these shortcomings. We applied this method to the Hoxc8 early enhancer region that has been previously analyzed in depth by functional methods and through its application successfully identified known protein binding cis-element motifs. Candidate protein binding sites could also be identified. This method, based on evolutionarily related sequence comparisons, should be quite useful as a prescreening step prior to functional analysis with corresponding savings in time and resources.

Protein sequence signatures support the African clade of mammals

DNA sequence evidence supports a superordinal clade of mammals that comprises elephants, sea cows, hyraxes, aardvarks, elephant shrews, golden moles, and tenrecs, which all have their origins in Africa, and therefore are dubbed Afrotheria. Morphologically, this appears an unlikely assemblage, which challenges-by including golden moles and tenrecs-the monophyly of the order Lipotyphla (Insectivora). We here identify in three proteins unique combinations of apomorphous amino acid replacements that support this clade. The statistical support for such "sequence signatures" as unambiguous synapomorphic evidence for the naturalness of the Afrotherian clade is reported. Using likelihood, combinatorial, and Bayesian methods we show that the posterior probability of the mammalian tree containing the Afrotherian clade is effectively 1.0, based on conservative assumptions. Presenting sequence data for another African insectivore, the otter shrew Micropotamogale lamottei, we demonstrate that such signatures are diagnostic for including newly investigated species in the Afrotheria. Sequence signatures provide "protein-morphological" synapomorphies that may aid in visualizing monophyletic groupings.

Protein sequence signatures support the African clade of mammals

DNA sequence evidence supports a superordinal clade of mammals that comprises elephants, sea cows, hyraxes, aardvarks, elephant shrews, golden moles, and tenrecs, which all have their origins in Africa, and therefore are dubbed Afrotheria. Morphologically, this appears an unlikely assemblage, which challenges-by including golden moles and tenrecs-the monophyly of the order Lipotyphla (Insectivora). We here identify in three proteins unique combinations of apomorphous amino acid replacements that support this clade. The statistical support for such "sequence signatures" as unambiguous synapomorphic evidence for the naturalness of the Afrotherian clade is reported. Using likelihood, combinatorial, and Bayesian methods we show that the posterior probability of the mammalian tree containing the Afrotherian clade is effectively 1.0, based on conservative assumptions. Presenting sequence data for another African insectivore, the otter shrew Micropotamogale lamottei, we demonstrate that such signatures are diagnostic for including newly investigated species in the Afrotheria. Sequence signatures provide "protein-morphological" synapomorphies that may aid in visualizing monophyletic groupings.

Bias in phylogenetic reconstruction of vertebrate rhodopsin sequences

Two spurious nodes were found in phylogenetic analyses of vertebrate rhodopsin sequences in comparison with well-established vertebrate relationships. These spurious reconstructions were well supported in bootstrap analyses and occurred independently of the method of phylogenetic analysis used (parsimony, distance, or likelihood). Use of this data set of vertebrate rhodopsin sequences allowed us to exploit established vertebrate relationships, as well as the considerable amount known about the molecular evolution of this gene, in order to identify important factors contributing to the spurious reconstructions. Simulation studies using parametric bootstrapping indicate that it is unlikely that the spurious nodes in the parsimony analyses are due to long branches or other topological effects. Rather, they appear to be due to base compositional bias at third positions, codon bias, and convergent evolution at nucleotide positions encoding the hydrophobic residues isoleucine, leucine, and valine. LogDet distance methods, as well as maximum-likelihood methods which allow for nonstationary changes in base composition, reduce but do not entirely eliminate support for the spurious resolutions. Inclusion of five additional rhodopsin sequences in the phylogenetic analyses largely corrected one of the spurious reconstructions while leaving the other unaffected. The additional sequences not only were more proximal to the corrected node, but were also found to have intermediate levels of base composition and codon bias as compared with neighboring sequences on the tree. This study shows that the spurious reconstructions can be corrected either by excluding third positions, as well as those encoding the amino acids Ile, Val, and Leu (which may not be ideal, as these sites can contain useful phylogenetic signal for other parts of the tree), or by the addition of sequences that reduce problems associated with convergent evolution.

Factors influencing phylogenetic inference: a case study using the mammalian carnivores

Phylogenetic reconstruction has undergone numerous developments in tree selection criteria (e.g., phenetics, cladistics, maximum-likelihood), available data sources (morphology versus molecules, and subsets of the latter), and practical limits on study size. Together with study age, I examined the effects of these variables on inferences of phylogeny for the mammalian order Carnivora. The raw data comprised 274 source trees spread among 13 carnivore taxa (generally families), which I divided into categories for each variable and combined using the supertree technique matrix representation with parsimony analysis. Incongruence between the resultant tree topologies or the underlying data was assessed using four comparison measures, each with slightly different properties: the triplet measures "do not conflict" and "explicitly agree," the partition metric, and the incongruence length difference metric. Except for a few cases reflecting historical problem areas in carnivore systematics, no significant differences in incongruence levels were found among the different categories within each variable, between the variables themselves, or between the taxa. Thus, most estimates of carnivore phylogeny cannot be distinguished from one another (and may even point toward the same solution) regardless of the methodology or data source employed. This conclusion held regardless of the comparison measure used.

Lineage-specific evolutionary rate in mammalian mtDNA

The existence of a lineage-specific nucleotide substitution rate in mammalian mtDNA has been investigated by analyzing the mtDNA of all available species, that is, 35 complete mitochondrial genomes from 14 mammalian orders. A detailed study of their evolutionary dynamics has been carried out on both ribosomal RNA and first and second codon positions (P12) of H-strand protein-coding genes by using two different types of relative-rate tests. Results are quite congruent between ribosomal and P12 sites. Significant rate variations have been observed among orders and among species of the same order. However, rate variation does not exceed 1.8-fold between the fastest (Proboscidea and Primates) and the slowest (Perissodactyla) evolving orders. Thus, the observed mitochondrial rate variations among taxa do not invalidate the suitability of mtDNA for drawing mammalian phylogeny. Dependence of evolutionary rate differences on variations in mutation and/or fixation rates was examined. Body size, generation time, and metabolic rate were tested, and no significant correlation was observed between them and the taxon-specific evolutionary rates, most likely because the latter might be influenced by multiple overlapping variable constraints.

Extracting species trees from complex gene trees: reconciled trees and vertebrate phylogeny

Paralogy is a pervasive problem in trying to use nuclear gene sequences to infer species phylogenies. One strategy for dealing with this problem is to infer species phylogenies from gene trees using reconciled trees, rather than directly from the sequences themselves. In this approach, the optimal species tree is the tree that requires the fewest gene duplications to be invoked. Because reconciled trees can identify orthologous from paralogous sequences, there is no need to do this prior to the analysis. Multiple gene trees can be analyzed simultaneously; however, the problem of nonuniform gene sampling raises practical problems which are discussed. In this paper the technique is applied to phylogenies for nine vertebrate genes (aldolase, alpha-fetoprotein, lactate dehydrogenase, prolactin, rhodopsin, trypsinogen, tyrosinase, vassopressin, and Wnt-7). The resulting species tree shows much similarity with currently accepted vertebrate relationships.

The dynamics of chromosome evolution in birds and mammals

Comparative mapping, which compares the location of homologous genes in different species, is a powerful tool for studying genome evolution. Comparative maps suggest that rates of chromosomal change in mammals can vary from one to ten rearrangements per million years. On the basis of these rates we would expect 84 to 600 conserved segments in a chicken comparison with human or mouse. Here we build comparative maps between these species and estimate that numbers of conserved segments are in the lower part of this range. We conclude that the organization of the human genome is closer to that of the chicken than the mouse and by adding comparative mapping results from a range of vertebrates, we identify three possible phases of chromosome evolution. The relative stability of genomes such as those of the chicken and human will enable the reconstruction of maps of ancestral vertebrates.

Model dependence of the phylogenetic inference: relationship among carnivores, Perissodactyls and cetartiodactyls as inferred from mitochondrial genome sequences

Some previous analysis of mitochondrial proteins strongly support the Carnivora/Perissodactyla grouping excluding Cetartiodactyla (Artiodactyla + Cetacea) as an outgroup, but the support of the hypothesis remains equivocal from the analysis of several nuclear-encoded proteins. In order to evaluate the strength of the support by mitochondrial proteins, phylogenetic relationship among Carnivora, Perissodactyla, and Cetartiodactyla was estimated with the ML method by using the updated data set of the 12 mitochondrial proteins with several alternative models. The analyses demonstrate that the phylogenetic inference depends on the model used in the ML analysis; i.e., whether the site-heterogeneity is taken into account and whether the rate parameters are estimated for each individual proteins or for the concatenated sequences. Although the analysis of concatenated sequences strongly supports the Carnivora/Perissodactyla grouping, the total evaluation of the separate analyses of individual proteins, which approximates the data better than the concatenated analysis, gives only ambiguous results, and therefore it is concluded that more data are needed to resolve this trichotomy.

The developing renal, reproductive, and respiratory systems of the African elephant suggest an aquatic ancestry

The early embryology of the elephant has never been studied before. We have obtained a rare series of African elephant (Loxodonta africana) embryos and fetuses ranging in weight from 0.04 to 18.5 g, estimated gestational ages 58-166 days (duration of gestation is approximately 660 days). Nephrostomes, a feature of aquatic vertebrates, were found in the mesonephric kidneys at all stages of development whereas they have never been recorded in the mesonephric kidneys of other viviparous mammals. The trunk was well developed even in the earliest fetus. The testes were intra-abdominal, and there was no evidence of a gubernaculum, pampiniform plexus, processus vaginalis, or a scrotum, confirming that the elephant, like the dugong, is one of the few primary testicond mammals. The palaeontological evidence suggests that the elephant's ancestors were aquatic, and recent immunological and molecular evidence shows an extremely close affinity between present-day elephants and the aquatic Sirenia (dugong and manatees). The evidence from our embryological study of the elephant also suggests that it evolved from an aquatic mammal.

Growing up with dinosaurs: molecular dates and the mammalian radiation

Dates of divergence derived from molecular data have been used to place the beginning of the radiation of modern mammalian orders in the Cretaceous, long before the final extinction of the dinosaurs. These molecular dates have been used to challenge the idea that the ordinal diversification of mammals was triggered by the availability of 'empty niches' left vacant by the disappearance of the dinosaurs. However, the broad discrepancies between molecular date estimates from different studies warn that molecular dates should not be accepted uncritically. Consideration of the wide confidence intervals around molecular date estimates, and the potential for geographic bias in the fossil record, could lessen the discrepancy between molecular and palaeontological data but might still prompt a re-evaluation of the timing and causes of the mammalian radiation.

Evolutionary and preservational constraints on origins of biologic groups: divergence times of eutherian mammals

Some molecular clock estimates of divergence times of taxonomic groups undergoing evolutionary radiation are much older than the groups' first observed fossil record. Mathematical models of branching evolution are used to estimate the maximal rate of fossil preservation consistent with a postulated missing history, given the sum of species durations implied by early origins under a range of species origination and extinction rates. The plausibility of postulated divergence times depends on origination, extinction, and preservation rates estimated from the fossil record. For eutherian mammals, this approach suggests that it is unlikely that many modern orders arose much earlier than their oldest fossil records.

Comparative studies on mammalian Hoxc8 early enhancer sequence reveal a baleen whale-specific deletion of a cis-acting element

Variations in regulatory regions of developmental control genes have been implicated in the divergence of axial morphologies. To find potentially significant changes in cis-regulatory regions, we compared nucleotide sequences and activities of mammalian Hoxc8 early enhancers. The nucleotide sequence of the early enhancer region is extremely conserved among mammalian clades, with five previously described cis-acting elements, A-E, being invariant. However, a 4-bp deletion within element C of the Hoxc8 early enhancer sequence is observed in baleen whales. When assayed in transgenic mouse embryos, a baleen whale enhancer (unlike other mammalian enhancers) directs expression of the reporter gene to more posterior regions of the neural tube but fails to direct expression to posterior mesoderm. We suggest that regulation of Hoxc8 in baleen whales differs from other mammalian species and may be associated with variation in axial morphology.