Mammalian evolution and the interphotoreceptor retinoid binding protein (IRBP) gene: Convincing evidence for several superordinal clades

Phylogenetic relationships of 25 mammalian species representing 17 of the 18 eutherian orders were examined using DNA sequences from a 1.2-kb region of the 5′ end of exon 1 of the single-copy nuclear gene known as interphotoreceptor retinoid binding protein (IRBP). A wide variety of methods of analysis of the DNA sequence, and of the translated products, all supported a five-order clade consisting of elephant shrew (Macroscelidea)/aardvark (Tubulidentata)/and the paenungulates (hyracoids, sirenians, and elephants), with bootstrap support in all cases of 100%. The Paenungulata was also strongly supported by these IRBP data. In the majority of analyses this monophyletic five-order grouping was the first branch off the tree after the Edentata. These results are highly congruent with two other recent sources of molecular data. Another superordinal grouping, with similar 100% bootstrap support in all of the same wide-ranging types of analyses, was Artiodactyla/Cetacea. Other superordinal affinities, suggested by the analyses, but with less convincing support, included a Perissodactyla/Artiodactyla/Cetacea clade, an Insectivora/Chiroptera clade, and Glires (an association of rodents and lagomorphs). Correspondence to: M. J. Stanhope     

Genomic inferences from Afrotheria and the evolution of elephants

Recent genetic studies have established that African forest and savanna elephants are distinct species with dissociated cytonuclear genomic patterns, and have identified Asian elephants from Borneo and Sumatra as conservation priorities. Representative of Afrotheria, a superordinal clade encompassing six eutherian orders, the African savanna elephant was among the first mammals chosen for whole-genome sequencing to provide a comparative understanding of the human genome. Elephants have large and complex brains and display advanced levels of social structure, communication, learning and intelligence. The elephant genome sequence might prove useful for comparative genomic studies of these advanced traits, which have appeared independently in only three mammalian orders: primates, cetaceans and proboscideans.     

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.     

Towards the delineation of the ancestral eutherian genome organization: comparative genome maps of human and the African elephant (Loxodonta africana) generated by chromosome painting

This study presents a whole-genome comparison of human and a representative of the Afrotherian clade, the African elephant, generated by reciprocal Zoo-FISH. An analysis of Afrotheria genomes is of special interest, because recent DNA sequence comparisons identify them as the oldest placental mammalian clade. Complete sets of whole-chromosome specific painting probes for the African elephant and human were constructed by degenerate oligonucleotide-primed PCR amplification of flow-sorted chromosomes. Comparative genome maps are presented based on their hybridization patterns. These maps show that the elephant has a moderately rearranged chromosome complement when compared to humans. The human paint probes identified 53 evolutionary conserved segments on the 27 autosomal elephant chromosomes and the X chromosome. Reciprocal experiments with elephant probes delineated 68 conserved segments in the human genome. The comparison with a recent aardvark and elephant Zoo-FISH study delineates new chromosomal traits which link the two Afrotherian species phylogenetically. In the absence of any morphological evidence the chromosome painting data offer the first non-DNA sequence support for an Afrotherian clade. The comparative human and elephant genome maps provide new insights into the karyotype organization of the proto-afrotherian, the ancestor of extant placental mammals, which most probably consisted of 2n=46 chromosomes.     

Identification of two major conformational aquaporin-4 epitopes for neuromyelitis optica autoantibody binding

Neuromyelitis optica (NMO) is an autoimmune demyelinating disease characterized by the presence of anti-aquaporin-4 (AQP4) antibodies in the patient sera. We recently reported that these autoantibodies are able to bind AQP4 when organized in the supramolecular structure called the orthogonal array of particles (OAP). To map the antigenic determinants, we produced a series of AQP4 mutants based on multiple alignment sequence analysis between AQP4 and other OAP-forming AQPs. Mutations were introduced in the three extracellular loops (A, C, and E), and the binding capacity of NMO sera was tested on AQP4 mutants. Results indicate that one group of sera was able to recognize a limited portion of loop C containing the amino acid sequence (146)GVT(T/M)V(150). A second group of sera was characterized by a predominant role of loop A. Deletion of four AQP4-specific amino acids ((61)G(S/T)E(N/K)(64)) in loop A substantially affected the binding of this group of sera. However, the binding capacity was further reduced when amino acids in loop A were mutated together with those in loop E or when those in loop C were mutated in combination with loop E. Finally, a series of AQP0 mutants were produced in which the extracellular loops were progressively changed to make them identical to AQP4. Results showed that none of the mutants was able to reproduce in AQP0 the NMO-IgG epitopes, indicating that the extracellular loop sequence by itself was not sufficient to determine the rearrangement required to create the epitopes. Although our data highlight the complexity of the disease, this study identifies key immunodominant epitopes and provides direct evidence that the transition from AQP4 tetramers to AQP4-OAPs involves conformational changes of the extracellular loops.     

Highly Congruent Molecular Support for a Diverse Superordinal Clade of Endemic African Mammals

A solution to higher level mammalian phylogeny is going to depend on the congruent establishment of superordinal groupings followed by a linking together of these clades. We present congruent and convincing evidence from four disparate nuclear protein coding genes and from a tandem alignment of the 12S–16S mitochondrial region, for a superordinal clade of endemic African mammals that includes elephant shrews, aardvarks, golden mole, elephants, sirenians, and hyraxes. Because of strong support for golden mole as part of this clade, the Insectivora are rendered paraphyletic or polyphyletic, with constrained monophyly of the insectivores judged significantly worse in the vast majority of tests. Branching arrangement within this clade remains highly uncertain; however, a tandem alignment of the protein coding genes suggests elephant shrew is the earliest African lineage. None of the individual data sets or combinations of data sets support the widely held view of a mirorder Tethytheria (Sirenia/Proboscidea), although only a tandem alignment of protein coding and mitochondrial loci significantly rejects this association. The majority of the data sets and analyses provide strong support for Caviomorpha as part of a monophyletic Rodentia.     

A molecular perspective on mammalian evolution from the gene encoding interphotoreceptor retinoid binding protein, with convincing evidence for bat monophyly.

The evolutionary relationships of the various orders of placental mammals remain an issue of uncertainty and controversy. Molecular studies of mammalian phylogeny at the DNA level that include more than just a few orders are still relatively meager. Here we report results on mammalian phylogeny deduced from the coding sequence of the single-copy nuclear gene for the interphotoreceptor retinoid binding protein (IRBP). Analysis of 13 species representing eight eutherian orders and one marsupial yielded results that falsify the hypothesis that megachiropteran bats are "flying primates," only convergently resembling microchiropteran bats. Instead, in agreement with more traditional views, as well as those from other recent molecular studies, the results strongly support a monophyletic Chiroptera (micro- and megabats grouped together). The IRBP results also offer some rare molecular support for the Glires concept, in which rodents and lagomorphs form a superordinal grouping. Also in congruence with other recent molecular evidence, IRBP sequences do not support the view of a superorder Archonta that includes Chiroptera along with Dermoptera (flying lemur), Scandentia (tree shrew), and Primates. IRBP was not however, without its shortcomings as a molecular phylogenetic system: high levels of homoplasy, evident in the marsupial outgroup, did not allow us to properly root the tree, and several of the higher level eutherian clades were only weakly supported (e.g., a Carnivora/Chiroptera clade and an Artiodactyla/Carnivora/Chiroptera clade). We suggest that these shortcomings may be diminished as the phylogenetic density of the data set is increased.     

Interordinal relationships and timescale of eutherian evolution as inferred from mitochondrial genome data

Extensive phylogenetic analyses of the updated sequence data of mammalian mitochondrial genomes were carried out using the maximum likelihood method in order to resolve deep branchings in eutherian evolution. The divergence times in the mammalian tree were estimated by a relaxed molecular clock of the mitochondrial proteins calibrated with multiple references. A Chiroptera/Eulipotyphla (i.e. bat/mole) clade and a close relationship of this clade to Fereuungulata (Carnivora+Perissodactyla+Cetartiodactyla) were reconfirmed with high statistical significance. However, a support for a monophyly of Fereuungulata relative to the Chiroptera/Eulipotyphla clade was fragile, and we suggest that the three branchings among Carnivora, Perissodactyla, Cetartiodactyla and Chiroptera/Eulipotyphla occurred successively in a short time period, estimated to be approximately 77Myr BP. The Chiroptera/Eulipotyphla divergence was estimated to roughly coincide with the Cretaceous-Tertiary boundary (65Myr BP). The monophyly of Rodentia, the Lagomorpha/Rodentia clade (traditionally called Glires), and the Afrotheria/Xenarthra clade were preferred over alternative relationships, but the supports of these clades were not strong enough to exclude other possibilities. Although several super-order taxa of eutherians were strongly supported by the analyses of the mitochondrial genome data, the branching order in the deepest part of the eutherian tree remained ambiguous from the data presently available.     

Housekeeping genes for phylogenetic analysis of eutherian relationships

The molecular relationship of placental mammals has attracted great interest in recent years. However, 2 crucial and conflicting hypotheses remain, one with respect to the position of the root of the eutherian tree and the other the relationship between the orders Rodentia, Lagomorpha (rabbits, hares), and Primates. Although most mitochondrial (mt) analyses have suggested that rodents have a basal position in the eutherian tree, some nuclear data in combination with mt-rRNA genes have placed the root on the so-called African clade or on a branch that includes this clade and the Xenarthra (e.g., anteater and armadillo). In order to generate a new and independent set of molecular data for phylogenetic analysis, we have established cDNA sequences from different tissues of various mammalian species. With this in mind, we have identified and sequenced 8 housekeeping genes with moderately fast rate of evolution from 22 placental mammals, representing 11 orders. In order to determine the root of the eutherian tree, the same genes were also sequenced for 3 marsupial species, which were used as outgroup. Inconsistent with the analyses of nuclear + mt-rRNA gene data, the current data set did not favor a basal position of the African clade or Xenarthra in the eutherian tree. Similarly, by joining rodents and lagomorphs on the same basal branch (Glires hypothesis), the data set is also inconsistent with the tree commonly favored in mtDNA analyses. The analyses of the currently established sequences have helped examination of problematic parts in the eutherian tree at the same time as they caution against suggestions that have claimed that basal eutherian relationships have been conclusively settled.     

Transmembrane helix 5 is critical for the high water permeability of aquaporin

Aquaporin-2 (AQP2), a vasopressin-regulated water channel, plays a major role in urinary concentration. AQP2 and the major intrinsic protein (MIP) of lens fiber are highly homologous (58% amino acid identity) and share a topology of six transmembrane helices connected by five loops (loops A-E). Despite the similarities of these proteins, however, the water channel activity of AQP2 is much higher than that of MIP. To determine the site responsible for this gain of activity in AQP2, several parts of MIP were replaced with the corresponding parts of AQP2. When expressed in Xenopus oocytes, the osmotic water permeability (P(f)) of MIP and AQP2 was 48 and 245 x 10(-)(4) cm/s, respectively. Substitutions in loops B-D failed to increase P(f), whereas substitution of loop E significantly increased P(f) 1.5-fold. A similar increase in P(f) was observed with the substitution of the front half of loop E. P(f) measurements taken in a yeast vesicle expression system also confirmed that loop E had a complementary effect, whereas loops B-D did not. However, P(f) values of the loop E chimeras were only approximately 30% of that of AQP2. Simultaneous exchanges of loop E and a distal half of transmembrane helix 5 just proximal to loop E increased P(f) to the level of that of AQP2. Replacement of helix 5 alone stimulated P(f) 2.7-fold. Conversely, P(f) was decreased by 73% when helix 5 of AQP2 was replaced with that of MIP. Moreover, P(f) was stimulated 2.6- and 3.3-fold after helix 5 of AQP1 and AQP4 was spliced into MIP, respectively. Our findings suggested that the distal half of helix 5 is necessary for maximum water channel activity in AQP. We speculate that this portion contributes to the formation of the aqueous pore and the determination of the flux rate.     

Structure and evolution of opossum,guinea pig,and porcupine cytochrome b genes

Summary We have sequenced the mitochondrial cytochrome b gene from the guinea pig, the African porcupine, and a South American opossum. A phylogenetic analysis, which includes 22 eutherian and four other vertebrate cytochrome b sequences, indicates that the guinea pig and the porcupine constitute a natural clade (Hystricomorpha) that is not a sister group to the clade of mice and rats (Myomorpha). Therefore, the hypothesis that the Rodentia is paraphyletic receives additional support. The artiodactyls, the perissodactyls, and the cetaceans form a group that is separated from the primates and the rodents. The 26 sequences are used to study the structure/function relationships in cytochrome b, whose function is electron transport. Most of the amino acid residues involved in the two reaction centers are well conserved in evolution. The four histidines that are believed to ligate the two hemes are invariant among the 26 sequences, but their nearby residues are not well conserved in evolution. The eight transmembrane domains represent some of the most divergent regions in the cytochrome b sequence. The rate of nonsynonymous substitution is considerably faster in the human and elephant lineages than in other eutherian lineages; the faster rate might be due to coevolution between cytochrome b and cytochrome c. Offprint requests to: W.-H. Li     

Predicting the magnitude of the reflex response to insertions in ubiquitin

The ability to predict the structural response of a protein to an insertion would be a significant advance for the fields of homology modeling and protein design. However, the effects of insertions on protein conformation are not well understood. Previous work has demonstrated that for two loops in ubiquitin, the primary determinant of the structural adaptation to insertions is the insertion site rather than the sequence of the insertion; this phenomenon was termed the reflex response of loops to insertions. We report herein the analysis of ubiquitin mutants with insertions in two other loops. This study demonstrates that the insertion site is the primary determinant of the response to insertions for these two new loops as well, which further supports the reflex response hypothesis. We also attempted to predict the relative magnitudes of the responses at each site but were unsuccessful. Using the additional data collected in this work, we have refined our predictive hypothesis.     

A Functional Aqp1 Gene Product Localizes on The Contractile Vacuole Complex in Paramecium multimicronucleatum

In a ciliate Paramecium, the presence of water channels on the membrane of contractile vacuole has long been predicted by both morphological and physiological data, however, to date either the biochemical or the molecular biological data have not been provided. In the present study, to examine the presence of aquaporin in Paramecium, we carried out RT-PCR with degenerated primers designed based on the ParameciumDB, and an aquaporin cDNA (aquaporin 1, aqp1) with a full-length ORF encoding 251 amino acids was obtained from Paramecium multimicronucleatum by using RACE. The deduced amino acid sequence of AQP1 had NPA-NPG motifs, and the prediction of protein secondary structure by CNR5000 and hydropathy plot showed the presence of six putative transmembrane domains and five connecting loops. Phylogenetic analysis results showed that the amino acid sequence of AQP1 was close to that of the Super-aquaporin group. The AQP1-GFP fusion protein clearly demonstrated the subcellular localization of AQP1 on the contractile vacuole complex, except for the decorated spongiome membrane. The functional analyses of aqp1 were done by RNA interference-based gene silencing, using an established feeding method. The aqp1 was found to be crucial for the total fluid output of the cell, the function of contractile vacuole membranes.     

Genome Size: A Novel Genomic Signature in Support of Afrotheria

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

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

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

Evidence of positive selection suggests possible role of aquaporins in the water-to-land transition of mudskippers

Aquaporins are integral membrane proteins that exchange water and small solutes. They played an important role in the colonisation of terrestrial environments by tetrapod ancestors via the appearance of three exclusive paralogs. Like early tetrapods, mudskippers represent an independent case of amphibious lifestyle evolution that is unparalleled by other extant fish groups. Given this lifestyle parallelism and that aquaporins were relevant for tetrapod terrestrialisation, this study examines the aquaporins in mudskippers to investigate whether similar changes in aquaporins could have possibly occurred during their water-to-land transition. We have catalogued aquaporin genes in four mudskipper genomes and studied their diversity and molecular evolution (including detection of positive selection) in a broad phylogenetic context of vertebrates. Our genomic screening returned 55 aquaporin genes for mudskippers (none of them constituting novel paralogs) that can be assigned to 10 different known classes. We detected signatures of positive selection in AQP10a and AQP11b in mudskippers (both the entire clade and the clade containing the most terrestrial species, implying different evolutionary times). This suggests possible alteration of the molecular function of such paralogs caused by changes at specific protein sequence positions, some of them located in relatively close proximity to parts of the molecule involved in pore formation and substrate selectivity. Given the importance of aquaporins for osmotic regulation in fishes, it might be possible that these selective changes (perhaps allowing permeability to new solutes) could have played a role during the adaptation of mudskippers to an amphibious lifestyle.     

Generation and Characterization of Antibodies against Asian Elephant (Elephas maximus) IgG,IgM, and IgA

Asian elephant (Elephas maximus) immunity is poorly characterized and understood. This gap in knowledge is particularly concerning as Asian elephants are an endangered species threatened by a newly discovered herpesvirus known as elephant endotheliotropic herpesvirus (EEHV), which is the leading cause of death for captive Asian elephants born after 1980 in North America. While reliable diagnostic assays have been developed to detect EEHV DNA, serological assays to evaluate elephant anti-EEHV antibody responses are lacking and will be needed for surveillance and epidemiological studies and also for evaluating potential treatments or vaccines against lethal EEHV infection. Previous studies have shown that Asian elephants produce IgG in serum, but they failed to detect IgM and IgA, further hampering development of informative serological assays for this species. To begin to address this issue, we determined the constant region genomic sequence of Asian elephant IgM and obtained some limited protein sequence information for putative serum IgA. The information was used to generate or identify specific commercial antisera reactive against IgM and IgA isotypes. In addition, we generated a monoclonal antibody against Asian elephant IgG. These three reagents were used to demonstrate that all three immunoglobulin isotypes are found in Asian elephant serum and milk and to detect antibody responses following tetanus toxoid booster vaccination or antibodies against a putative EEHV structural protein. The results indicate that these new reagents will be useful for developing sensitive and specific assays to detect and characterize elephant antibody responses for any pathogen or vaccine, including EEHV.     

A Molecular Phylogeny of Plesiorycteropus Reassigns the Extinct Mammalian Order ‘Bibymalagasia’

Madagascar is well known for its diverse fauna and flora, being home to many species not found anywhere else in the world. However, its biodiversity in the recent past included a range of extinct enigmatic fauna, such as elephant birds, giant lemurs and dwarfed hippopotami. The ‘Malagasy aardvark’ (Plesiorycteropus) has remained one of Madagascar’s least well-understood extinct species since its discovery in the 19^th century. Initially considered a close relative of the aardvark (Orycteropus) within the order Tubulidentata, more recent morphological analyses challenged this placement on the grounds that the identifiably derived traits supporting this allocation were adaptations to digging rather than shared ancestry. Because the skeletal evidence showed many morphological traits diagnostic of different eutherian mammal orders, they could not be used to resolve its closest relatives. As a result, the genus was tentatively assigned its own taxonomic order ‘Bibymalagasia’, yet how this order relates to other eutherian mammal orders remains unclear despite numerous morphological investigations. This research presents the first known molecular sequence data for Plesiorycteropus, obtained from the bone protein collagen (I), which places the ‘Malagasy aardvark’ as more closely related to tenrecs than aardvarks. More specifically, Plesiorycteropus was recovered within the order Tenrecoidea (golden moles and tenrecs) within Afrotheria, suggesting that the taxonomic order ‘Bibymalagasia’ is obsolete. This research highlights the potential for collagen sequencing in investigating the phylogeny of extinct species as a viable alternative to ancient DNA (aDNA) sequencing, particularly in cases where aDNA cannot be recovered.