Marsupials and Eutherians reunited: genetic evidence for the Theria hypothesis of mammalian evolution

Cooperative inflammation: The recruitment of inflammatory signaling in marsupial and eutherian pregnancy

The evolution of viviparity in therian mammals, i.e. marsupials and "placental" mammals, occurred by retention of the conceptus in the female reproductive tract and precocious "hatching" from the shell coat. Both eutherian embryo implantation and the opossum embryo attachment reaction are evolutionarily derived from and homologous to a defensive inflammatory process induced after shell coat hatching. However, both lineages, marsupials and placental mammals, have modified the inflammatory response substantially. We review the induction, maintenance, and effects of inflammation throughout pregnancy, with special attention to the role of prostaglandins and the mucosal inflammatory response, both of which likely had roles in early mammalian viviparity. We propose that the key step was not only suppression of the inflammatory response after implantation in placental mammals, but also the transfer of the inflammatory cell-cell communication network to a different set of cell types than in generic inflammation. To support this conclusion we discuss evidence that pro-inflammatory signal production in the opossum is not limited to maternal cells, as expected in bona fide defensive inflammation, but also includes fetal tissues, in a process we term cooperative inflammation. The ways in which the inflammatory reaction was independently modified in these two lineages helps explain major life history differences between extant marsupials and eutherians.

The Type 3 Deiodinase: Epigenetic Control of Brain Thyroid Hormone Action and Neurological Function

Thyroid hormones (THs) influence multiple processes in the developing and adult central nervous system, and their local availability needs to be maintained at levels that are tailored to the requirements of their biological targets. The local complement of TH transporters, deiodinase enzymes, and receptors is critical to ensure specific levels of TH action in neural cells. The type 3 iodothyronine deiodinase (DIO3) inactivates THs and is highly present in the developing and adult brain, where it limits their availability and action. DIO3 deficiency in mice results in a host of neurodevelopmental and behavioral abnormalities, demonstrating the deleterious effects of TH excess, and revealing the critical role of DIO3 in the regulation of TH action in the brain. The fact the Dio3 is an imprinted gene and that its allelic expression pattern varies across brain regions and during development introduces an additional level of control to deliver specific levels of hormone action in the central nervous system (CNS). The sensitive epigenetic nature of the mechanisms controlling the genomic imprinting of Dio3 renders brain TH action particularly susceptible to disruption due to exogenous treatments and environmental exposures, with potential implications for the etiology of human neurodevelopmental disorders.

Transcriptomic Changes Associated with Pregnancy in a Marsupial, the Gray Short-Tailed Opossum Monodelphis domestica

Live birth has emerged as a reproductive strategy many times across vertebrate evolution; however, mammals account for the majority of viviparous vertebrates. Marsupials are a mammalian lineage that last shared a common ancestor with eutherians (placental mammals) over 148 million years ago. Marsupials are noted for giving birth to highly altricial young after a short gestation, and represent humans’ most distant viviparous mammalian relatives. Here we ask what insight can be gained into the evolution of viviparity in mammals specifically and vertebrates in general by analyzing the global uterine transcriptome in a marsupial. Transcriptome analyses were performed using NextGen sequencing of uterine RNA samples from the gray short-tailed opossum, Monodelphis domestica. Samples were collected from late stage pregnant, virgin, and non-pregnant experienced breeders. Three different algorithms were used to determine differential expression, and results were confirmed by quantitative PCR. Over 900 opossum gene transcripts were found to be significantly more abundant in the pregnant uterus than non-pregnant, and over 1400 less so. Most with increased abundance were genes related to metabolism, immune systems processes, and transport. This is the first study to characterize the transcriptomic differences between pregnant, non-pregnant breeders, and virgin marsupial uteruses and helps to establish a set of pregnancy-associated genes in the opossum. These observations allowed for comparative analyses of the differentially transcribed genes with other mammalian and non-mammalian viviparous species, revealing similarities in pregnancy related gene expression over 300 million years of amniote evolution.

Micro-computed tomography and microdissection of the temporal bone of tree shrews

OBJECTIVE

To understand the morphology and anatomical data of the temporal bone of tree shrews through micro-computed tomography (micro-CT) and microdissection.

METHODS

Skull specimens from 10 tree shrews were scanned using micro-CT examination. The acquired images were used for three-dimensional reconstruction and measurement using the Mimics 10.01 software. Twenty tree shrews were subjected to microdissection and the data were measured.

RESULTS

Micro-CT and three-dimensional reconstruction could clearly define the three-dimensional spatial position of the ear structure. Micro-CT and microdissection showed that the otic vesicles of the tree shrews were located on both sides of posterior-inferior skull bone. The location of the otic vesicles was superficial, and the bone was thin. All of the structures of the middle and inner ear of the tree shrews were well developed. The ossicular chain was differentiated into the malleus, incus and stapes. The location of the three semi-circular canals of the tree shrews was superficial and easy to dissect. In vivo, the three semi-circular canals were easy to localize and the surface bone was thin. The contour and structure of the cochlea and number of cochlear turns were similar to those in humans.

CONCLUSION

This study could provide anatomical data to allow tree shrews to be used as animal models for studying ear diseases.

The Evolution and Structure of Atypical T Cell Receptors

The T cell receptor structure and genetic organization have been thought to have been stable in vertebrate evolution relative to the immunoglobulins. For the most part, this has been true and the content and organization of T cell receptor genes has been fairly conserved over the past 400 million years of gnathostome evolution. Analyses of TCRδ chains in a broad range of vertebrate lineages over the past decade have revealed a remarkable and previously unrealized degree of plasticity. This plasticity can generally be described in two forms. The first is broad use of antibody heavy chain variable genes in place of the conventional Vδ. The second form containing an unusual three extracellular domain structures has evolved independently in both cartilaginous fishes and mammals. Two well-studied vertebrate lineages, the eutherian mammals such as mice and humans and teleost fishes, lack any of these alternative TCR forms, contributing to why they went undiscovered for so long after the initial description of the conventional TCR chains three decades ago. This chapter describes the state of knowledge of these unusual TCR forms, both their structure and genetics, and current ideas on their function.

Molecular evolution of a novel marsupial S100 protein (S100A19) which is expressed at specific stages of mammary gland and gut development

S100 proteins are calcium-binding proteins involved in controlling diverse intracellular and extracellular processes such as cell growth, differentiation, and antimicrobial function. We recently identified a S100-like cDNA from the tammar wallaby (Macropus eugenii) stomach. Phylogentic analysis shows wallaby S100A19 forms a new clade with other marsupial and monotreme S100A19, while this group shows similarity to eutherian S100A7 and S100A15 genes. This is also supported by amino acid and domain comparisons. We show S100A19 is developmentally-regulated in the tammar wallaby gut by demonstrating the gene is expressed in the forestomach of young animals at a time when the diet consists of only milk, but is absent in older animals when the diet is supplemented with herbage. During this transition the forestomach phenotype changes from a gastric stomach into a fermentation sac and intestinal flora changes with diet. We also show that S100A19 is expressed in the mammary gland of the tammar wallaby only during specific stages of lactation; the gene is up-regulated during pregnancy and involution and not expressed during the milk production phase of lactation. Comparison of the tammar wallaby S100A19 protein sequence with S100 protein sequences from eutherian, monotreme and other marsupial species suggest the marsupial S100A19 has two functional EF hand domains, and an extended His tail. An evolutionary analysis of S100 family proteins was carried out to gain a better understanding of the relationship between the S100 family member functions. We propose that S100A19 gene/protein is the ancestor of the eutherian S100A7 gene/protein, which has subsequently modified its original function in eutherians. This modified function may have arisen due to differentiation of evolutionary pressures placed on gut and mammary gland developmental during mammal evolution. The highly regulated differential expression patterns of S100A19 in the tammar wallaby suggests that S100A19 may play a role in gut development, which differs between metatherians and eutherians, and/or include a potential antibacterial role in order to establish the correct flora and protect against spiral bacteria in the immature forestomach. In the mammary gland it may protect the tissue from infection at times of vulnerability during the lactation cycle.

Expression of antisense of insulin-like growth factor-2 receptor RNA non-coding (AIRN) during early gestation in cattle

The insulin-like growth factor type 2 receptor (IGF2R) regulates fetal growth by removing IGF2 from circulation. In mice, expression of the Igf2r gene is only imprinted after implantation and is associated with expression of the antisense non-coding (nc)RNA, Airn. The objectives of this study were, first, to determine if bovine AIRN was expressed during developmentally important stages of gestation, and second, to determine if expression of bAIRN was affected by method of embryo production. Control reactions confirmed that sequence verified bAIRN PCR amplicons resulted from RNA within the sample and not from genomic DNA contamination. IGF2R mRNA was expressed in all fetal liver samples at Days 35-55 and 70 of gestation as well as in 8 of 9 Day 15 conceptuses, 10 of 10 Day 18 conceptuses, and in all day 7 blastocyst pools. bAIRN was expressed in all samples of fetal liver at Days 35-55 and 70 of gestation. The proportion of conceptuses that expressed bAIRN increased from 1 of 9 at Day 15 of gestation to 8 of 10 at Day 18 of gestation. No bAIRN was expressed in any blastocyst pools. The relative level of bAIRN was greater (P<0.05) in fetal liver from embryos produced in vivo compared to that from embryos produced in vitro. In summary bAIRN was not expressed in blastocyst-stage embryos, was expressed in an increasing proportion of embryos around the time of maternal recognition of pregnancy and was expressed following implantation. Furthermore, relative levels of bAIRN in bovine fetal liver can be altered by method of embryo production.

A single glycine-alanine exchange directs ligand specificity of the elephant progestin receptor

The primary gestagen of elephants is 5α-dihydroprogesterone (DHP), which is unlike all other mammals studied until now. The level of DHP in elephants equals that of progesterone in other mammals, and elephants are able to bind DHP with similar affinity to progesterone indicating a unique ligand-binding specificity of the elephant progestin receptor (PR). Using site-directed mutagenesis in combination with in vitro binding studies we here report that this change in specificity is due to a single glycine to alanine exchange at position 722 (G722A) of PR, which specifically increases DHP affinity while not affecting binding of progesterone. By conducting molecular dynamics simulations comparing human and elephant PR ligand-binding domains (LBD), we observed that the alanine methyl group at position 722 is able to push the DHP A-ring into a position similar to progesterone. In the human PR, the DHP A-ring position is twisted towards helix 3 of PR thereby disturbing the hydrogen bond pattern around the C3-keto group, resulting in a lower binding affinity. Furthermore, we observed that the elephant PR ligand-binding pocket is more rigid than the human analogue, which probably explains the higher affinity towards both progesterone and DHP. Interestingly, the G722A substitution is not elephant-specific, rather it is also present in five independent lineages of mammalian evolution, suggesting a special role of the substitution for the development of distinct mammalian gestagen systems.

Convergent and divergent evolution of genomic imprinting in the marsupial Monodelphis domestica

BACKGROUND

Genomic imprinting is an epigenetic phenomenon resulting in parent-of-origin specific monoallelic gene expression. It is postulated to have evolved in placental mammals to modulate intrauterine resource allocation to the offspring. In this study, we determined the imprint status of metatherian orthologues of eutherian imprinted genes.

RESULTS

L3MBTL and HTR2A were shown to be imprinted in Monodelphis domestica (the gray short-tailed opossum). MEST expressed a monoallelic and a biallelic transcript, as in eutherians. In contrast, IMPACT, COPG2, and PLAGL1 were not imprinted in the opossum. Differentially methylated regions (DMRs) involved in regulating imprinting in eutherians were not found at any of the new imprinted loci in the opossum. Interestingly, a novel DMR was identified in intron 11 of the imprinted IGF2R gene, but this was not conserved in eutherians. The promoter regions of the imprinted genes in the opossum were enriched for the activating histone modification H3 Lysine 4 dimethylation.

CONCLUSIONS

The phenomenon of genomic imprinting is conserved in Therians, but the marked difference in the number and location of imprinted genes and DMRs between metatherians and eutherians indicates that imprinting is not fully conserved between the two Therian infra-classes. The identification of a novel DMR at a non-conserved location as well as the first demonstration of histone modifications at imprinted loci in the opossum suggest that genomic imprinting may have evolved in a common ancestor of these two Therian infra-classes with subsequent divergence of regulatory mechanisms in the two lineages.

Evolution of immune functions of the mammary gland and protection of the infant

Abstract The evolution of immunological agents in milk is intertwined with the general aspects of the evolution of the mammary gland. In that respect, mammalian precursors emerged from basal amniotes some 300 million years ago. In contrast to the predominant dinosaurs, proto-mammals possessed a glandular skin. A secondary palate in the roof of the mouth that directed airflow from the nostrils to the oropharynx and thus allowed mammals to ingest and breathe simultaneously first appeared in cynodonts 230 million years ago. This set the stage for mammalian newborns to nurse from the future mammary gland. Interplays between environmental and genetic changes shaped mammalian evolution including the mammary gland from dermal glands some 160 millions of years ago. It is likely that secretions from early mammary glands provided nutrients and immunological agents for the infant. Natural selection culminated in milks uniquely suited to nourish and protect infants of each species. In human milk, antimicrobial, anti-inflammatory, and immunoregulatory agents and living leukocytes are qualitatively or quantitatively different from those in other mammalian milks. Those in human milk compensate for developmental delays in the immunological system of the recipient infant. Consequently, the immune system in human milk provided by evolution is much of the basis for encouraging breastfeeding for human infants.

Genomic imprinting of the type 3 thyroid hormone deiodinase gene: regulation and developmental implications

BACKGROUND

In recent years, findings in a number of animal and human models have ignited renewed interest in the type 3 deiodinase (D3), the main enzyme responsible for the inactivation of thyroid hormones. The induction of D3 in models of illness and injury has raised critical questions about the physiological significance of reduced thyroid hormone availability in those states. Phenotypes in transgenic mice lacking this enzyme also point to important developmental roles for D3. A critical determinant of D3 expression is genomic imprinting, an epigenetic phenomenon that regulates a small number of dosage-critical genes in the mammalian genome. The D3 gene (Dio3) is imprinted and preferentially expressed from one of the alleles in most tissues.

SCOPE OF REVIEW

In the context of the physiological significance of D3 and the characteristics and purported origins of genomic imprinting, we review the current knowledge about the epigenetic mechanisms specifying gene dosage in the Dio3 locus.

MAJOR CONCLUSIONS

Altered Dio3 dosage is detrimental to development, suggesting that the level of thyroid hormone action needs to be exquisitely tailored in a timely fashion to the requirements of particular tissues. An appropriate Dio3 dosage is the result of the coordinated action of certain genomic elements and epigenetic marks in the Dlk1-Dio3 domain.

GENERAL SIGNIFICANCE

The imprinting of Dio3 prompts intriguing questions about why the level of thyroid hormone signaling should be regulated in this rare epigenetic manner, and to what extent altered Dio3 expression due to aberrant imprinting may be implicated in human conditions. This article is part of a Special Issue entitled Thyroid hormone signalling.

Detecting the symplesiomorphy trap: a multigene phylogenetic analysis of terebelliform annelids

BACKGROUND

For phylogenetic reconstructions, conflict in signal is a potential problem for tree reconstruction. For instance, molecular data from different cellular components, such as the mitochondrion and nucleus, may be inconsistent with each other. Mammalian studies provide one such case of conflict where mitochondrial data, which display compositional biases, support the Marsupionta hypothesis, but nuclear data confirm the Theria hypothesis. Most observations of compositional biases in tree reconstruction have focused on lineages with different composition than the majority of the lineages under analysis. However in some situations, the position of taxa that lack compositional bias may be influenced rather than the position of taxa that possess compositional bias. This situation is due to apparent symplesiomorphic characters and known as "the symplesiomorphy trap".

RESULTS

Herein, we report an example of the sympleisomorphy trap and how to detect it. Worms within Terebelliformia (sensu Rouse & Pleijel 2001) are mainly tube-dwelling annelids comprising five 'families': Alvinellidae, Ampharetidae, Terebellidae, Trichobranchidae and Pectinariidae. Using mitochondrial genomic data, as well as data from the nuclear 18S, 28S rDNA and elongation factor-1α genes, we revealed incongruence between mitochondrial and nuclear data regarding the placement of Trichobranchidae. Mitochondrial data favored a sister relationship between Terebellidae and Trichobranchidae, but nuclear data placed Trichobranchidae as sister to an Ampharetidae/Alvinellidae clade. Both positions have been proposed based on morphological data.

CONCLUSIONS

Our investigation revealed that mitochondrial data of Ampharetidae and Alvinellidae exhibited strong compositional biases. However, these biases resulted in a misplacement of Trichobranchidae, rather than Alvinellidae and Ampharetidae. Herein, we document that Trichobranchidae was apparently caught in the symplesiomorphy trap suggesting that in certain situations even homologies can be misleading.

An atlas of the speed of copy number changes in animal gene families and its implications

The notion that gene duplications generating new genes and functions is commonly accepted in evolutionary biology. However, this assumption is more speculative from theory rather than well proven in genome-wide studies. Here, we generated an atlas of the rate of copy number changes (CNCs) in all the gene families of ten animal genomes. We grouped the gene families with similar CNC dynamics into rate pattern groups (RPGs) and annotated their function using a novel bottom-up approach. By comparing CNC rate patterns, we showed that most of the species-specific CNC rates groups are formed by gene duplication rather than gene loss, and most of the changes in rates of CNCs may be the result of adaptive evolution. We also found that the functions of many RPGs match their biological significance well. Our work confirmed the role of gene duplication in generating novel phenotypes, and the results can serve as a guide for researchers to connect the phenotypic features to certain gene duplications.

Burst of young retrogenes and independent retrogene formation in mammals

Retroposition and retrogenes gain increasing attention as recent studies show that they play an important role in human new gene formation. Here we examined the patterns of retrogene distribution in 8 mammalian genomes using 4 non-mammalian genomes as a contrast. There has been a burst of young retrogenes not only in primate lineages as suggested in a recent study, but also in other mammalian lineages. In mammals, most of the retrofamilies (the gene families that have retrogenes) are shared between species. In these shared retrofamilies, 14%–18% of functional retrogenes may have originated independently in multiple mammalian species. Notably, in the independently originated retrogenes, there is an enrichment of ribosome related gene function. In sharp contrast, none of these patterns hold in non-mammals. Our results suggest that the recruitment of the specific L1 retrotransposons in mammals might have been an important evolutionary event for the split of mammals and non-mammals and retroposition continues to be an important active process in shaping the dynamics of mammalian genomes, as compared to being rather inert in non-mammals.

The opossum genome: insights and opportunities from an alternative mammal

The strategic importance of the genome sequence of the gray, short-tailed opossum, Monodelphis domestica, accrues from both the unique phylogenetic position of metatherian (marsupial) mammals and the fundamental biologic characteristics of metatherians that distinguish them from other mammalian species. Metatherian and eutherian (placental) mammals are more closely related to one another than to other vertebrate groups, and owing to this close relationship they share fundamentally similar genetic structures and molecular processes. However, during their long evolutionary separation these alternative mammals have developed distinctive anatomical, physiologic, and genetic features that hold tremendous potential for examining relationships between the molecular structures of mammalian genomes and the functional attributes of their components. Comparative analyses using the opossum genome have already provided a wealth of new evidence regarding the importance of noncoding elements in the evolution of mammalian genomes, the role of transposable elements in driving genomic innovation, and the relationships between recombination rate, nucleotide composition, and the genomic distributions of repetitive elements. The genome sequence is also beginning to enlarge our understanding of the evolution and function of the vertebrate immune system, and it provides an alternative model for investigating mechanisms of genomic imprinting. Equally important, availability of the genome sequence is fostering the development of new research tools for physical and functional genomic analyses of M. domestica that are expanding its versatility as an experimental system for a broad range of research applications in basic biology and biomedically oriented research.

Confirming the phylogeny of mammals by use of large comparative sequence data sets

The ongoing generation of prodigious amounts of genomic sequence data from myriad vertebrates is providing unparalleled opportunities for establishing definitive phylogenetic relationships among species. The size and complexities of such comparative sequence data sets not only allow smaller and more difficult branches to be resolved but also present unique challenges, including large computational requirements and the negative consequences of systematic biases. To explore these issues and to clarify the phylogenetic relationships among mammals, we have analyzed a large data set of over 60 megabase pairs (Mb) of high-quality genomic sequence, which we generated from 41 mammals and 3 other vertebrates. All sequences are orthologous to a 1.9-Mb region of the human genome that encompasses the cystic fibrosis transmembrane conductance regulator gene (CFTR). To understand the characteristics and challenges associated with phylogenetic analyses of such a large data set, we partitioned the sequence data in several ways and utilized maximum likelihood, maximum parsimony, and Neighbor-Joining algorithms, implemented in parallel on Linux clusters. These studies yielded well-supported phylogenetic trees, largely confirming other recent molecular phylogenetic analyses. Our results provide support for rooting the placental mammal tree between Atlantogenata (Xenarthra and Afrotheria) and Boreoeutheria (Euarchontoglires and Laurasiatheria), illustrate the difficulty in resolving some branches even with large amounts of data (e.g., in the case of Laurasiatheria), and demonstrate the valuable role that very large comparative sequence data sets can play in refining our understanding of the evolutionary relationships of vertebrates.

An analysis of IGFBP evolution

OBJECTIVES

To perform a synonymous, non-synonymous codon mutational analysis of the IGFBP gene family and identify mechanisms by which the IGFBP subfamilies diverged.

METHODS

We identified 78 intact nucleotide sequences from 6 IGFBP subfamilies and 12 different species and used them for phylogenetic and synonymous, non-synonymous codon mutational analysis. Deletion and insertion comparisons were performed across subfamilies to determine if this might play a unique role in subfamily genesis.

RESULTS

IGFBP-2 was identified by phylogenetic analysis to be the most related subfamily of the IGFBP progenitor, followed by IGFBP-4. Insertions and deletions within the variable domains were associated with divergence of each subfamily from its progenitor, suggesting a common motif for IGFBP evolution. Insertions unique to mammals were also found within the amino terminus of IGFBP-2.

CONCLUSION

IGFBP subfamily divergence is associated with variable domain insertion or deletion and vigorous non-synonymous codon mutation. Our findings suggest strong selective pressure for IGFBP divergence in terrestrial vertebrates.

Loss of genes implicated in gastric function during platypus evolution

BACKGROUND

The duck-billed platypus (Ornithorhynchus anatinus) belongs to the mammalian subclass Prototheria, which diverged from the Theria line early in mammalian evolution. The platypus genome sequence provides a unique opportunity to illuminate some aspects of the biology and evolution of these animals.

RESULTS

We show that several genes implicated in food digestion in the stomach have been deleted or inactivated in platypus. Comparison with other vertebrate genomes revealed that the main genes implicated in the formation and activity of gastric juice have been lost in platypus. These include the aspartyl proteases pepsinogen A and pepsinogens B/C, the hydrochloric acid secretion stimulatory hormone gastrin, and the alpha subunit of the gastric H+/K+-ATPase. Other genes implicated in gastric functions, such as the beta subunit of the H+/K+-ATPase and the aspartyl protease cathepsin E, have been inactivated because of the acquisition of loss-of-function mutations. All of these genes are highly conserved in vertebrates, reflecting a unique pattern of evolution in the platypus genome not previously seen in other mammalian genomes.

CONCLUSION

The observed loss of genes involved in gastric functions might be responsible for the anatomical and physiological differences in gastrointestinal tract between monotremes and other vertebrates, including small size, lack of glands, and high pH of the monotreme stomach. This study contributes to a better understanding of the mechanisms that underlie the evolution of the platypus genome, might extend the less-is-more evolutionary model to monotremes, and provides novel insights into the importance of gene loss events during mammalian evolution.

Mammalian evolution and biomedicine: new views from phylogeny

Recent progress resolving the phylogenetic relationships of the major lineages of mammals has had a broad impact in evolutionary biology, comparative genomics and the biomedical sciences. Novel insights into the timing and historical biogeography of early mammalian diversification have resulted from a new molecular tree for placental mammals coupled with dating approaches that relax the assumption of the molecular clock. We highlight the numerous applications to come from a well-resolved phylogeny and genomic prospecting in multiple lineages of mammals, from identifying regulatory elements in mammalian genomes to assessing the functional consequences of mutations in human disease loci and those driving adaptive evolution.

Visualizing differences in phylogenetic information content of alignments and distinction of three classes of long-branch effects

BACKGROUND

Published molecular phylogenies are usually based on data whose quality has not been explored prior to tree inference. This leads to errors because trees obtained with conventional methods suppress conflicting evidence, and because support values may be high even if there is no distinct phylogenetic signal. Tools that allow an a priori examination of data quality are rarely applied.

RESULTS

Using data from published molecular analyses on the phylogeny of crustaceans it is shown that tree topologies and popular support values do not show existing differences in data quality. To visualize variations in signal distinctness, we use network analyses based on split decomposition and split support spectra. Both methods show the same differences in data quality and the same clade-supporting patterns. Both methods are useful to discover long-branch effects. We discern three classes of long branch effects. Class I effects consist of attraction of terminal taxa caused by symplesiomorphies, which results in a false monophyly of paraphyletic groups. Addition of carefully selected taxa can fix this effect. Class II effects are caused by drastic signal erosion. Long branches affected by this phenomenon usually slip down the tree to form false clades that in reality are polyphyletic. To recover the correct phylogeny, more conservative genes must be used. Class III effects consist of attraction due to accumulated chance similarities or convergent character states. This sort of noise can be reduced by selecting less variable portions of the data set, avoiding biases, and adding slower genes.

CONCLUSION

To increase confidence in molecular phylogenies an exploratory analysis of the signal to noise ratio can be conducted with split decomposition methods. If long-branch effects are detected, it is necessary to discern between three classes of effects to find the best approach for an improvement of the raw data.

Calibration Choice, Rate Smoothing, and the Pattern of Tetrapod Diversification According to the Long Nuclear Gene RAG-1

A phylogeny of tetrapods is inferred from nearly complete sequences of the nuclear RAG-1 gene sampled across 88 taxa encompassing all major clades, analyzed via parsimony and Bayesian methods. The phylogeny provides support for Lissamphibia, Theria, Lepidosauria, a turtle-archosaur clade, as well as most traditionally accepted groupings. This tree allows simultaneous molecular clock dating for all tetrapod groups using a set of well-corroborated calibrations. Relaxed clock (PLRS) methods, using the amniote = 315 Mya (million years ago) calibration or a set of consistent calibrations, recovers reasonable divergence dates for most groups. However, the analysis systematically underestimates divergence dates within archosaurs. The bird-crocodile split, robustly documented in the fossil record as being around approximately 245 Mya, is estimated at only approximately 190 Mya, and dates for other divergences within archosaurs are similarly underestimated. Archosaurs, and particulary turtles have slow apparent rates possibly confounding rate modeling, and inclusion of calibrations within archosaurs (despite their high deviances) not only improves divergence estimates within archosaurs, but also across other groups. Notably, the monotreme-therian split ( approximately 210 Mya) matches the fossil record; the squamate radiation ( approximately 190 Mya) is younger than suggested by some recent molecular studies and inconsistent with identification of approximately 220 and approximately 165 Myo (million-year-old) fossils as acrodont iguanians and approximately 95 Myo fossils colubroid snakes; the bird-lizard (reptile) split is considerably older than fossil estimates (< or = 285 Mya); and Sphenodon is a remarkable phylogenetic relic, being the sole survivor of a lineage more than a quarter of a billion years old. Comparison with other molecular clock studies of tetrapod divergences suggests that the common practice of enforcing most calibrations as minima, with a single liberal maximal constraint, will systematically overestimate divergence dates. Similarly, saturation of mitochondrial DNA sequences, and the resultant greater compression of basal branches means that using only external deep calibrations will also lead to inflated age estimates within the focal ingroup.

Reliable safety assessment depends on species differences in epigenetic mechanisms of gene regulation

The potential carcinogenic hazard of chemical agents to humans is presently based primarily on the results of long-term animal bioassays. The validity of this toxicologic approach to human risk assessment depends on two fundamental assumptions. First, the results of an animal bioassay are directly applicable to humans (interspecies extrapolation). Second, the doses used in an animal bioassay are relevant for estimating risk at known or expected human exposure levels (dose extrapolation). Although progress has been made over the past four decades in understanding the mode of action of chemical carcinogens, it is increasingly important to determine mechanistically the relevance of these modes of action in humans. There is now evidence that M6P/IGF2R functions as a novel tumor-suppressor gene in a variety of human and rodent cancers. M6p/Igf2r is imprinted in rodents and expressed only from the maternal allele after embryonic implantation. In contrast, both alleles are functional in humans. This marked species difference in M6P/IGF2R imprinting has important implications for human carcinogen risk assessment since only one rather than two alleles needs to be mutated in rodents to completely inactivate the function of this tumor suppressor gene. This striking species difference in the imprint status of M6P/IGF2R clearly demonstrates that we need to understand better variations in epigenetic mechanisms of gene regulation between rodents and humans to perform accurately chemical safety assessments.

Altered epigenetic variance in surviving litters from nutritionally restricted lactating primiparous sows

Feed restriction of primiparous sows during the last week of lactation has been shown to decrease embryonic growth and female embryo survival to Day 30 of gestation. This study sought to determine whether global DNA methylation and epigenetic gene expression of the candidate genes Igf2, Igf2r, and Xist were associated with these treatment effects. Given that these epigenetic traits are expected to be important for embryo viability, changes in variance for these traits at Day 30 were predicted to be reflected in the loss of abnormal embryos at this time. Consistent with this prediction, variance in DNA methylation was reduced (P < 0.001) in Restrict male embryo, and there was a tendency for reduced variance (P < 0.06) in Restrict female embryos. Variation in DNA methylation tended to be correlated (R = 0.42, P < 0.1) with the difference in variance of embryo weights between treatments (P < 0.01), suggesting a relationship between epigenetic changes and embryonic development. Variance in Igf2r expression tended to decrease (P < 0.07) in Restrict female embryos while variance in Xist expression tended to decrease in Restrict male embryos (P < 0.08), suggesting that maternally inherited epigenetic defects may cause female embryonic loss and reduced growth before Day 30 of gestation.

Cross-species clues of an epigenetic imprinting regulatory code for the IGF2R gene

The epigenetic marks on the IGF2R gene that encodes a receptor responsible for IGF-II degradation consist of differentially methylated DNA in association with multiple modifications on the associated histones. We review these epigenetic marks across various species during the evolution of IGF2R imprinting. Both IGF2 and IGF2R genesare imprinted in the mammal lineage that diverged from Monotremata approximately 150 million years ago. While IGF2 is consistently imprinted in all mammals following its divergence, IGF2R imprinting disappears in the Euarchonta lineage, including human species, approximately 75 million years ago. Differential DNA methylation marks on the two parental alleles correlate with imprinting in all imprinted genes including IGF2R. While the DNA methylation marks in the IGF2R promoter region 1 (DMR1) correlate with IGF2R allelic expression, the DNA methylation marks in the intron region 2 (DMR2) fail to correlate with IGF2R imprinting status in a number of species. Human IGF2R and mouse neuronal Igf2r are not imprinted despite the presence of DMR2. We have noted that human IGF2R is not imprinted in more than 100 informative samples including various tumor tissues. Furthermore, opossum (Marsupialia) IGF2R is consistently imprinted despite the absence of DMR2. These lines of evidence indicate that DNA methylation marks in DMR2 are neither necessary nor sufficient for consistent imprinting of IGF2R across species. Histone modification marks, however, correlate more consistently with the tissue-specific and species-specific imprinting status of IGF2R in human and mouse. Acetylated histone H3 and H4 and methylated lysine 4 of H3 (H3-K4Me) associate with transcriptionally active alleles while tri-methylated lysine 9 of H3 (H3-K9Me3) marks the silenced alleles. In the mouse, an antisense non-coding transcript called Air is transcribed from DMR2 on the paternal allele, and this imprinted transcript plays a central role in Igf2r imprinting. Mouse Igf2r imprinting depends on an Air RNA while the existence of AIR in other species is unknown. Overall, DNA methylation, histone acetylation, and histone methylation play a vital role in coordinating IGF2R allelic expression across all species. Rare monoallelic or skewed allelic expression of human IGF2R and their biological importance warrants further rigorous study.

Retroposed elements as archives for the evolutionary history of placental mammals

Reconstruction of the placental mammalian (eutherian) evolutionary tree has undergone diverse revisions, and numerous aspects remain hotly debated. Initial hierarchical divisions based on morphology contained many misgroupings due to features that evolved independently by similar selection processes. Molecular analyses corrected many of these misgroupings and the superordinal hierarchy of placental mammals was recently assembled into four clades. However, long or rapid evolutionary periods, as well as directional mutation pressure, can produce molecular homoplasies, similar characteristics lacking common ancestors. Retroposed elements, by contrast, integrate randomly into genomes with negligible probabilities of the same element integrating independently into orthologous positions in different species. Thus, presence/absence analyses of these elements are a superior strategy for molecular systematics. By computationally scanning more than 160,000 chromosomal loci and judiciously selecting from only phylogenetically informative retroposons for experimental high-throughput PCR applications, we recovered 28 clear, independent monophyly markers that conclusively verify the earliest divergences in placental mammalian evolution. Using tests that take into account ancestral polymorphisms, multiple long interspersed elements and long terminal repeat element insertions provide highly significant evidence for the monophyletic clades Boreotheria (synonymous with Boreoeutheria), Supraprimates (synonymous with Euarchontoglires), and Laurasiatheria. More importantly, two retropositions provide new support for a prior scenario of early mammalian evolution that places the basal placental divergence between Xenarthra and Epitheria, the latter comprising all remaining placentals. Due to its virtually homoplasy-free nature, the analysis of retroposon presence/absence patterns avoids the pitfalls of other molecular methodologies and provides a rapid, unequivocal means for revealing the evolutionary history of organisms.

The authors identified and sequenced retroposons in mammalian genomes. The presence and absence of these retroposons provided evolutionary markers from which the authors reconstructed the phylogenetic history of placental mammals.

Status and applications of genomic resources for the gray, short-tailed opossum, Monodelphis domestica, an American marsupial model for comparative biology

Owing to its small size, favourable reproductive characteristics, and simple husbandry, the gray, short-tailed opossum, Monodelphis domestica, has become the most widely distributed and intensively utilised laboratory-bred research marsupial in the world today. This article provides an overview of the current state and future projections of genomic resources for this species and discusses the potential impact of this growing resource base on active research areas that use M. domestica as a model system. The resources discussed include: fully arrayed, bacterial artificial chromosome (BAC) libraries; an expanding linkage map; developing full-genome BAC-contig and chromosomal fluorescence in situ hybridisation maps; public websites providing access to the M. domestica whole-genome-shotgun sequence trace database and the whole-genome sequence assembly; and a new project underway to create an expressed-sequence database and microchip expression arrays for functional genomics applications. Major research areas discussed span a variety of genetic, evolutionary, physiologic, reproductive, developmental, and behavioural topics, including: comparative immunogenetics; genomic imprinting; reproductive biology; neurobiology; photobiology and carcinogenesis; genetics of lipoprotein metabolism; developmental and behavioural endocrinology; sexual differentiation and development; embryonic and fetal development; meiotic recombination; genome evolution; molecular evolution and phylogenetics; and more.

The Platypus Is in Its Place: Nuclear Genes and Indels Confirm the Sister Group Relation of Monotremes and Therians

Morphological data supports monotremes as the sister group of Theria (extant marsupials + eutherians), but phylogenetic analyses of 12 mitochondrial protein-coding genes have strongly supported the grouping of monotremes with marsupials: the Marsupionta hypothesis. Various nuclear genes tend to support Theria, but a comprehensive study of long concatenated sequences and broad taxon sampling is lacking. We therefore determined sequences from six nuclear genes and obtained additional sequences from the databases to create two large and independent nuclear data sets. One (data set I) emphasized taxon sampling and comprised five genes, with a concatenated length of 2,793 bp, from 21 species (two monotremes, six marsupials, nine placentals, and four outgroups). The other (data set II) emphasized gene sampling and comprised eight genes and three proteins, with a concatenated length of 10,773 bp or 3,669 amino acids, from five taxa (a monotreme, a marsupial, a rodent, human, and chicken). Both data sets were analyzed by parsimony, minimum evolution, maximum likelihood, and Bayesian methods using various models and data partitions. Data set I gave bootstrap support values for Theria between 55% and 100%, while support for Marsupionta was at most 12.3%. Taking base compositional bias into account generally increased the support for Theria. Data set II exclusively supported Theria, with the highest possible values and significantly rejected Marsupionta. Independent phylogenetic evidence in support of Theria was obtained from two single amino acid deletions and one insertion, while no supporting insertions and deletions were found for Marsupionta. On the basis of our data sets, the time of divergence between Monotremata and Theria was estimated at 231-217 MYA and between Marsupialia and Eutheria at 193-186 MYA. The morphological evidence for a basal position of Monotremata, well separated from Theria, is thus fully supported by the available molecular data from nuclear genes.

Comparative sequencing provides insights about the structure and conservation of marsupial and monotreme genomes

Sequencing and comparative analyses of genomes from multiple vertebrates are providing insights about the genetic basis for biological diversity. To date, these efforts largely have focused on eutherian mammals, chicken, and fish. In this article, we describe the generation and study of genomic sequences from noneutherian mammals, a group of species occupying unusual phylogenetic positions. A large sequence data set (totaling >5 Mb) was generated for the same orthologous region in three marsupial (North American opossum, South American opossum, and Australian tammar wallaby) and one monotreme (platypus) genomes. These ancient mammalian genomes are characterized by unusual architectural features with respect to G + C and repeat content, as well as compression relative to human. Approximately 14% and 34% of the human sequence forms alignments with the orthologous sequence from platypus and the marsupials, respectively; these numbers are distinctly lower than that observed with nonprimate eutherian mammals (45-70%). The alignable sequences between human and each marsupial species are not completely overlapping (only 80% common to all three species) nor are the platypus-alignable sequences completely contained within the marsupial-alignable sequences. Phylogenetic analysis of synonymous coding positions reveals that platypus has a notably long branch length, with the human-platypus substitution rate being on average 55% greater than that seen with human-marsupial pairs. Finally, analyses of the major mammalian lineages reveal distinct patterns with respect to the common presence of evolutionarily conserved vertebrate sequences. Our results confirm that genomic sequence from noneutherian mammals can contribute uniquely to unraveling the functional and evolutionary histories of the mammalian genome.

Mammalian phylogenomics comes of age

The relatively new field of phylogenomics is beginning to reveal the potential of genomic data for evolutionary studies. As the cost of whole genome sequencing falls, anticipation of complete genome sequences from divergent species, reflecting the major lineages of modern mammals, is no longer a distant dream. In this article, we describe how comparative genomic data from mammals is progressing to resolve long-standing phylogenetic controversies, to refine dogma on how chromosomes evolve and to guide annotation of human and other vertebrate genomes.

Adaptive evolution of HoxA-11 and HoxA-13 at the origin of the uterus in mammals

The evolution of morphological characters is mediated by the evolution of developmental genes. Evolutionary changes can either affect cis-regulatory elements, leading to differences in their temporal and spatial regulation, or affect the coding region. Although there is ample evidence for the importance of cis-regulatory evolution, it has only recently been shown that transcription factors do not remain functionally equivalent during evolution. These results suggest that the evolution of transcription factors may play an active role in the evolution of development. To test this idea we investigated the molecular evolution of two genes essential for the development and function of the mammalian female reproductive organs, HoxA-11 and HoxA-13. We predicted that if coding-region evolution plays an active role in developmental evolution, then these genes should have experienced adaptive evolution at the origin of the mammalian female reproductive system. We report the sequences of HoxA-11 from basal mammalian and amniote taxa and analyse HoxA-11 and HoxA-13 for signatures of adaptive molecular evolution. The data demonstrate that these genes were under strong positive (directional) selection in the stem lineage of therian and eutherian mammals, coincident with the evolution of the uterus and vagina. These results support the idea that adaptive evolution of transcription factors can be an integral part in the evolution of novel structures.

Evolution of the vertebral formulae in mammals: A perspective on developmental constraints

Developmental constraints refer to biases that limit phenotypic changes during evolution. To examine the contribution of developmental constraints in the evolution of vertebrate morphology, we analyzed the distribution pattern of mammalian vertebral formulae. Data on mammalian vertebral formulae were collected from the Descriptive Catalogue of the Osteological Series Contained in the Museum of the Royal College of Surgeons of England by Richard Owen (1853) and were plotted onto the most reliable mammalian phylogenetic tree based on recent molecular studies. In addition to the number of cervical vertebrae that is almost fixed to 7, we found that the number of thoracolumbar vertebrae tends to be 19 in many groups of mammals. Since fidelity of the number of thoracolumbar vertebrae was also completely maintained in Monotremata and Marsupialia, we presumed that thoracolumbar vertebral number as well as cervical vertebral number might have been fixed in the primitive mammalian lineage. On the basis of primitive vertebral formulae, we could clarify the polarity of evolution and identify several deviations from the primitive states during the mammalian evolution. The changes in the vertebral formulae in eutherian mammals seem to be lineage-specific, such that most species in Carnivora have 20 instead of 19 thoracolumbar vertebrae. Because such lineage-specific vertebral formulae contrast with the estimated distribution pattern on the assumption of evolution only through the selective pressure, we concluded that developmental constraints played an important role in the evolution of mammalian vertebral formulae.

First-generation linkage map of the gray, short-tailed opossum, Monodelphis domestica, reveals genome-wide reduction in female recombination rates

The gray, short-tailed opossum, Monodelphis domestica, is the most extensively used, laboratory-bred marsupial resource for basic biologic and biomedical research worldwide. To enhance the research utility of this species, we are building a linkage map, using both anonymous markers and functional gene loci, that will enable the localization of quantitative trait loci (QTL) and provide comparative information regarding the evolution of mammalian and other vertebrate genomes. The current map is composed of 83 loci distributed among eight autosomal linkage groups and the X chromosome. The autosomal linkage groups appear to encompass a very large portion of the genome, yet span a sex-average distance of only 633.0 cM, making this the most compact linkage map known among vertebrates. Most surprising, the male map is much larger than the female map (884.6 cM vs. 443.1 cM), a pattern contrary to that in eutherian mammals and other vertebrates. The finding of genome-wide reduction in female recombination in M. domestica, coupled with recombination data from two other, distantly related marsupial species, suggests that reduced female recombination might be a widespread metatherian attribute. We discuss possible explanations for reduced female recombination in marsupials as a consequence of the metatherian characteristic of determinate paternal X chromosome inactivation.

The complexity of expressed kappa light chains in egg-laying mammals

Complementary DNAs encoding immunoglobulin light chains were isolated from two monotreme species, Ornithorhynchus anatinus (duckbill platypus) and Tachyglossus aculeatus (echidna). The sequences of both the variable and constant regions of these clones had greater similarity to IGK than to other light chain classes and phylogenetic analyses place them squarely within the mammalian IGK group, establishing them as monotreme IGK homologues. The constant region sequences of all clones were essentially identical within each species and, along with Southern blot results, the data are consistent with a single IGKC in each species. The expressed IGKV repertoires from both platypus and echidna were randomly sampled and there appear to be at least four platypus and at least nine echidna IGKV subgroups. The IGKV subgroups are highly divergent within species, in some cases sharing as little as 57% nucleotide identity. Two of the IGKV subgroups are present in both species, so there is some degree of overlap in the germline repertoires of these two monotremes. Overall the complexity seen in platypus and echidna IGK light chains is comparable with that of other mammals considered to have high levels of germline diversity and is in contrast to what has been found so far for monotreme IGL.

Modeling the impact of DNA methylation on the evolution of BRCA1 in mammals

The modified base 5-methylcytosine ((m)C) plays an important functional role in the biology of mammals as an epigenetic modification and appears to exert a striking impact on the molecular evolution of mammal genomes. The collective epigenetic functions of (m)C revolve around its effect on gene transcription, while the influence of this modified base on the evolution of mammal genomes derives from the greatly elevated spontaneous mutation rate of (m)C to T. In mammals, (m)C occurs at the dinucleotides CpG, CpA, and CpT. As a step toward a comprehensive statistical examination of the role of (m)C in mammal molecular evolution, we have developed novel Markov models of codon substitution that incorporate dinucleotide-level terms relevant to (m)C mutation. We apply these models to two data sets of aligned BRCA1 exon 11 sequences from bats and primates. In all cases, terms specific to mutations that affect the dinucleotides CpG, CpA, and CpT significantly improved model fit. For the CpG-specific terms, both transition and transversion substitution rates were elevated. These rates differed between the data sets. Bats exhibited a lower relative rate of substitutions at CpG-containing codons. Transition substitutions were significantly less than 1 at CpA-containing codons but greater than 1 at CpT-containing codons. The inclusion of interaction terms in the codon models to represent possible confounding with the effect of natural selection were supported for codons that contained CpG and CpT, but not CpA. From the results, we infer that mutation of (m)C is a probable factor that affects BRCA1 codons containing the dinucleotide CpG, a possible factor for CpA-containing codons, and an unlikely factor that affects CpT-containing codons. The confounding of estimated terms with the effect of natural selection indicate this confounding must be addressed for comparisons between different coding and noncoding regions.

Cloning of IgE from the echidna (Tachyglossus aculeatus) and a comparative analysis of epsilon chains from all three extant mammalian lineages

In continuation of our evolutionary studies of immunoglobulin (Ig) expression, we present here the cloning of IgE from a monotreme, the short-beaked echidna (Tachyglossus aculeatus). Including echidna IgE, 15 epsilon chain sequences have been isolated and each of the three mammalian lineages (placentals, marsupials and monotremes) is now represented by at least two sequences. Phylogenetic analyses based on all available epsilon chains and a selection of other mammalian Ig isotypes (IgM, IgA and IgG) were generated using three different algorithms. The resulting trees strongly support the Theria hypothesis, which states that the monotreme lineage was the first of the three extant mammalian lineages to appear in evolution. Furthermore, to increase our understanding of IgE we have done a detailed comparative analysis, with focus on primary structure, potential N-glycosylation, charge distribution and conservation of residues in the putative receptor-binding site. The overall structure of IgE, i.e. four constant domains and the positions of putative disulfide-bridge formations, are conserved, as is an N-glycosylation site in the third constant domain. An increased homology was observed in the putative receptor-binding site, which suggests an important function for the IgE/Fc epsilon RI interaction. IgE has been found exclusively in mammals, but it is present in all extant mammalian lineages. This, together with the overall conservation of structure, indicates that IgE appeared as a separate isotype early in mammalian evolution and that structural maintenance may have a selective advantage.

The monotreme genome: a patchwork of reptile, mammal and unique features?

The first specimen of platypus (Ornithorhynchus anatinus) that reached Britain in the late 18th century was regarded a scientific hoax. Over decades the anatomical characteristics of these unique mammals, such as egg laying and the existence of mammary glands, were hotly debated before they were accepted. Within the last 40 years, more and more details of monotreme physiology, histology, reproduction and genetics have been revealed. Some show similarities with birds or reptiles, some with therian mammals, but many are very specific to monotremes. The genome is no exception to monotreme uniqueness. An early opinion was that the karyotype, composed of a few large chromosomes and many small ones, resembled bird and reptile macro- and micro-chromosomes. However, the platypus genome also features characteristics that are not present in other mammals, such as a complex translocation system. The sex chromosome system is still not resolved. Nothing is known about dosage compensation and, unlike in therian mammals, there seems to be no genomic imprinting. In this article we will recount the mysteries of the monotreme genome and describe how we are using recently developed technology to identify chromosomes in mitosis, meiosis and sperm, to map genes to chromosomes, to unravel the sex chromosome system and the translocation chain and investigate X inactivation and genomic imprinting in monotremes.

Review of the monotreme fossil record and comparison of palaeontological and molecular data

Monotremes have traditionally been considered a remnant group of mammals descended from archaic Mesozoic stock, surviving to the present day on the relatively isolated Australian continent. Challenges to this orthodoxy have been spurred by discoveries of 'advanced' Cretaceous monotremes (Steropodon galmani, Archer, M., et al., 1985. First Mesozoic mammal from Australia-an Early Cretaceous monotreme, Nature. 318, 363-366) as well as by results from molecular data linking monotremes to therian mammals (specifically to marsupials in some studies). This paper reviews the monotreme fossil record and briefly discusses significant new information from additional Cretaceous Australian material. Mesozoic monotremes (including S. galmani) were a diverse group as evidenced by new material from the Early Cretaceous of New South Wales and Victoria currently under study. Although most of these new finds are edentulous jaws (limiting dental comparisons and determination of dietary niches), a range of sizes and forms has been determined. Some of these Cretaceous jaws exhibit archaic features-in particular evidence for the presence of a splenial bone in S. galmani-not seen in therian mammals or in post-Mesozoic (Tertiary and Quaternary) monotreme taxa. Tertiary monotremes were either archaic ornithorhynchids (toothed platypuses in the genera Monotrematum and Obdurodon) or tachyglossids (large echidnas in the genera Megalibgwilia and Zaglossus). Quaternary ornithorhynchid material is referable to the sole living platypus species Ornithorhynchus anatinus. Quaternary echidnas, however, were moderately diverse and several forms are known (Megalibgwilia species; 'Zaglossus' hacketti; Zaglossus species and Tachyglossus aculeatus).

A despecialization step underlying evolution of a family of serine proteases

In the trypsin superfamily of serine proteases, non-trypsin-like primary specificities have arisen in only two monophyletic descendent subbranches. We have recreated an ancestor to one of these subbranches (granzyme) using phylogenetic inference, gene synthesis, and protein expression. This ancestor has two unusual properties. First, it has broad primary specificity encompassing the entire repertoire of novel primary specificities found in its descendents. Second, unlike extant members that have narrow primary specificities, the ancestor exhibits tolerance to mutational changes in primary specificity-conferring residues-that is, structural plasticity. Molecular modeling and mutagenesis studies indicate that these unusual properties are due to a particularly wide substrate binding pocket. These two crucial properties of the ancestor not only distinguish it from its extant descendents but also from the trypsin-like proteases that preceded it. This indicates that a despecialization step, characterized by broad specificity and structural plasticity, underlies evolution of new primary specificities in this protease superfamily.

Complete nucleotide sequence of the mitochondrial genome of a salamander, Mertensiella luschani

The complete nucleotide sequence (16,650 bp) of the mitochondrial genome of the salamander Mertensiella luschani (Caudata, Amphibia) was determined. This molecule conforms to the consensus vertebrate mitochondrial gene order. However, it is characterized by a long non-coding intervening sequence with two 124-bp repeats between the tRNA(Thr) and tRNA(Pro) genes. The new sequence data were used to reconstruct a phylogeny of jawed vertebrates. Phylogenetic analyses of all mitochondrial protein-coding genes at the amino acid level recovered a robust vertebrate tree in which lungfishes are the closest living relatives of tetrapods, salamanders and frogs are grouped together to the exclusion of caecilians (the Batrachia hypothesis) in a monophyletic amphibian clade, turtles show diapsid affinities and are placed as sister group of crocodiles+birds, and the marsupials are grouped together with monotremes and basal to placental mammals. The deduced phylogeny was used to characterize the molecular evolution of vertebrate mitochondrial proteins. Amino acid frequencies were analyzed across the main lineages of jawed vertebrates, and leucine and cysteine were found to be the most and least abundant amino acids in mitochondrial proteins, respectively. Patterns of amino acid replacements were conserved among vertebrates. Overall, cartilaginous fishes showed the least variation in amino acid frequencies and replacements. Constancy of rates of evolution among the main lineages of jawed vertebrates was rejected.

Evolution of the major histocompatibility complex: Isolation of class II beta cDNAs from two monotremes, the platypus and the short-beaked echidna

Extant mammals are composed of three lineages: the eutherians, the marsupials and the monotremes. The majority of the mammalian major histocompatibility complex (MHC) data is based on the eutherian mammals, which generally have three classical MHC class II beta chain gene clusters - DRB, DQB and DPB, as well as the non-classical DMB and DOB. Marsupial DMB, DAB and DBB have been characterised. Confusion still surrounds the relationship of the marsupial DAB and DBB genes with the classical eutherian class II clusters. Here we present the first monotreme MHC class II beta chain sequences. Four MHC class II beta chain sequences were isolated from a spleen cDNA library from the short-beaked echidna, and one from a spleen cDNA library from platypus using a brushtail possum DAB probe. Given the non-orthologous relationship of the monotreme sequences with marsupial and eutherian beta chain clusters, we recommend that the five new monotreme sequences be assigned the nomenclature 'DZB', signifying the description of a new mammalian beta chain cluster. Our analysis suggests that all mammalian beta chain sequences (except DMB) evolved from a common ancestor. Maximum likelihood analysis places the monotreme beta chain sequences at the base of the mammalian clade, indicating their ancestral status. However, within the mammalian clade, monophyletic clades are not robust, and elucidation of the order of gene duplication that gave rise to the present-day gene clusters is not yet possible.

The root of the mammalian tree inferred from whole mitochondrial genomes

Morphological and molecular data are currently contradictory over the position of monotremes with respect to marsupial and placental mammals. As part of a re-evaluation of both forms of data we examine complete mitochondrial genomes in more detail. There is a particularly large discrepancy in the frequencies of thymine and cytosine (T-C) between mitochondrial genomes that appears to affect some deep divergences in the mammalian tree. We report that recoding nucleotides to RY-characters, and partitioning maximum-likelihood analyses among subsets of data reduces such biases, and improves the fit of models to the data, respectively. RY-coding also increases the signal on the internal branches relative to external, and thus increases the phylogenetic signal. In contrast to previous analyses of mitochondrial data, our analyses favor Theria (marsupials plus placentals) over Marsupionta (monotremes plus marsupials). However, a short therian stem lineage is inferred, which is at variance with the traditionally deep placement of monotremes on morphological data.

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).

The evolution of tribospheny and the antiquity of mammalian clades

The evolution of tribosphenic molars is a key innovation in the history of Mammalia. Tribospheny allows for both shearing and grinding occlusal functions. Marsupials and placentals are advanced tribosphenic mammals (i.e., Theria) that show additional modifications of the tribosphenic dentition including loss of the distal metacristid and development of double-rank postvallum/prevallid shear. The recent discovery of Eomaia [Nature 416 (2002) 816], regarded as the oldest eutherian mammal, implies that the marsupial-placental split is at least 125 million years old. The conventional scenario for the evolution of tribosphenic and therian mammals hypothesizes that each group evolved once, in the northern hemisphere, and is based on a predominantly Laurasian fossil record. With the recent discovery of the oldest tribosphenic mammal (Ambondro) from the Mesozoic of Gondwana, Flynn et al. [Nature 401 (1999) 57] suggested that tribospheny evolved in Gondwana rather than in Laurasia. Luo et al. [Nature 409 (2001) 53; Acta Palaeontol. Pol. 47 (2002) 1] argued for independent origins of tribospheny in northern (Boreosphenida) and southern (Australosphenida) hemisphere clades, with the latter including Ambondro, ausktribosphenids, and monotremes. Here, we present cladistic evidence for a single origin of tribosphenic molars. Further, Ambondro may be a stem eutherian, making the split between marsupials and placentals at least 167 m.y. old. To test this hypothesis, we used the relaxed molecular clock approach of Thorne/Kishino with amino acid data sets for BRCA1 [J. Mammal. Evol. 8 (2001) 239] and the IGF2 receptor [Mammal. Genome 12 (2001) 513]. Point estimates for the marsupial-placental split were 182-190 million years based on BRCA1 and 185-187 million years based on the IGF2 receptor. These estimates are fully compatible with the results of our cladistic analyses.

Characterization of beta(2)-microglobulin coding sequence from three non-placental mammals: the duckbill platypus, the short-beaked echidna, and the grey short-tailed opossum

To further characterize genes of immunological importance from non-placental mammals, cDNAs encoding beta(2)-microglobulin (beta(2)m) were isolated from two prototherians, the platypus and an echidna, and one metatherian, a grey short-tailed opossum. In addition, a second allele of beta(2)m was identified in another metatherian species, the brushtail possum. Analysis of the deduced translations revealed conservation of key residues in these molecules over a long evolutionary history. The types of nucleotide substitutions present among the various taxa are also consistent with purifying selection at this conserved locus. An evolutionary tree of beta(2)m was constructed that supports the classic view of evolution with prototherians as the basal mammalian group.