Parallel adaptive radiations in two major clades of placental mammals

Molecular cloning and evolutionary analysis of the GJA1 (connexin43) gene from bats (Chiroptera)

Gap junction protein connexin43 (Cx43), encoded by the GJA1 gene, is the most abundant connexin in the cardiovascular system and was reported as a crucial factor maintaining cardiac electrical conduction, as well as having a very important function in facilitating the recycling of potassium ions from hair cells in the cochlea back into the cochlear endolymph during auditory transduction processes. In mammals, bats are the only taxon possessing powered flight, placing exceptional demand on many organismal processes. To meet the demands of flying, the hearts of bats show many specialties. Moreover, ultrasonic echolocation allows bat species to orientate and often detect and locate food in darkness. In this study, we cloned the full-length coding region of GJA1 gene from 12 different species of bats and obtained orthologous sequences from other mammals. We used the maximum likelihood method to analyse the evolution of GJA1 gene in mammals and the lineage of bats. Our results showed this gene is much conserved in mammals, as well as in bats' lineage. Compared with other mammals, we found one private amino acid substitution shared by bats, which is located on the inner loop domain, as well as some species-specific amino acid substitutions. The evolution rate analyses showed the signature of purifying selection on not only different classification level lineages but also the different domains and amino acid residue sites of this gene. Also, we suggested that GJA1 gene could be used as a good molecular marker to do the phylogenetic reconstruction.

Loss of teeth and enamel in tetrapods: fossil record, genetic data and morphological adaptations

Since their recruitment in the oral cavity, approximately 450 million years ago, teeth have been subjected to strong selective constraints due to the crucial role that they play in species survival. It is therefore quite surprising that the ability to develop functional teeth has subsequently been lost several times, independently, in various lineages. In this review, we concentrate our attention on tetrapods, the only vertebrate lineage in which several clades lack functional teeth from birth to adulthood. Indeed, in other lineages, teeth can be absent in adults but be functionally present in larvae and juveniles, can be absent in the oral cavity but exist in the pharyngeal region, or can develop on the upper jaw but be absent on the lower jaw. Here, we analyse the current data on toothless (edentate) tetrapod taxa, including information available on enamel-less species. Firstly, we provide an analysis of the dispersed and fragmentary morphological data published on the various living taxa concerned (and their extinct relatives) with the aim of tracing the origin of tooth or enamel loss, i.e. toads in Lissamphibia, turtles and birds in Sauropsida, and baleen whales, pangolins, anteaters, sloths, armadillos and aardvark in Mammalia. Secondly, we present current hypotheses on the genetic basis of tooth loss in the chicken and thirdly, we try to answer the question of how these taxa have survived tooth loss given the crucial importance of this tool. The loss of teeth (or only enamel) in all of these taxa was not lethal because it was always preceded in evolution by the pre-adaptation of a secondary tool (beak, baleens, elongated adhesive tongues or hypselodonty) useful for improving efficiency in food uptake. The positive selection of such secondary tools would have led to relaxed functional constraints on teeth and would have later compensated for the loss of teeth. These hypotheses raise numerous questions that will hopefully be answered in the near future.

Extreme convergence in stick insect evolution: phylogenetic placement of the Lord Howe Island tree lobster

The 'tree lobsters' are an enigmatic group of robust, ground-dwelling stick insects (order Phasmatodea) from the subfamily Eurycanthinae, distributed in New Guinea, New Caledonia and associated islands. Its most famous member is the Lord Howe Island stick insect Dryococelus australis (Montrouzier), which was believed to have become extinct but was rediscovered in 2001 and is considered to be one of the rarest insects in the world. To resolve the evolutionary position of Dryococelus, we constructed a phylogeny from approximately 2.4 kb of mitochondrial and nuclear sequence data from representatives of all major phasmatodean lineages. Our data placed Dryococelus and the New Caledonian tree lobsters outside the New Guinean Eurycanthinae as members of an unrelated Australasian stick insect clade, the Lanceocercata. These results suggest a convergent origin of the 'tree lobster' body form. Our reanalysis of tree lobster characters provides additional support for our hypothesis of convergent evolution. We conclude that the phenotypic traits leading to the traditional classification are convergent adaptations to ground-living behaviour. Our molecular dating analyses indicate an ancient divergence (more than 22 Myr ago) between Dryococelus and its Australian relatives. Hence, Dryococelus represents a long-standing separate evolutionary lineage within the stick insects and must be regarded as a key taxon to protect with respect to phasmatodean diversity.

Retroposon analysis and recent geological data suggest near-simultaneous divergence of the three superorders of mammals

As a consequence of recent developments in molecular phylogenomics, all extant orders of placental mammals have been grouped into 3 lineages: Afrotheria, Xenarthra, and Boreotheria, which originated in Africa, South America, and Laurasia, respectively. Despite this advancement, the order of divergence of these 3 lineages remains unresolved. Here, we performed extensive retroposon analysis with mammalian genomic data. Surprisingly, we identified a similar number of informative retroposon loci that support each of 3 possible phylogenetic hypotheses: the basal position for Afrotheria (22 loci), Xenarthra (25 loci), and Boreotheria (21 loci). This result indicates that the divergence of the placental common ancestor into the 3 lineages occurred nearly simultaneously. Thus, we examined whether these molecular data could be integrated into the geological context by incorporating recent geological data. We obtained firm evidence that complete separation of Gondwana into Africa and South America occurred 120 +/- 10 Ma. Accordingly, the previous reported time frame (division of Pangea into Gondwana and Laurasia at 148-138 Ma and division of Gondwana at 105 Ma) cannot be used to validate mammalian divergence order. Instead, we use our retroposon results and the recent geological data to propose that near-simultaneous divisions of continents leading to isolated Africa, South America, and Laurasia caused nearly concomitant divergence of the ancient placental ancestor into 3 lineages, Afrotheria, Xenarthra, and Boreotheria, approximately 120 Ma.

Morphological diversity and the roles of contingency, chance and determinism in african cichlid radiations

Background

Deterministic evolution, phylogenetic contingency and evolutionary chance each can influence patterns of morphological diversification during adaptive radiation. In comparative studies of replicate radiations, convergence in a common morphospace implicates determinism, whereas non-convergence suggests the importance of contingency or chance.

Methodology/Principal Findings

The endemic cichlid fish assemblages of the three African great lakes have evolved similar sets of ecomorphs but show evidence of non-convergence when compared in a common morphospace, suggesting the importance of contingency and/or chance. We then analyzed the morphological diversity of each assemblage independently and compared their axes of diversification in the unconstrained global morphospace. We find that despite differences in phylogenetic composition, invasion history, and ecological setting, the three assemblages are diversifying along parallel axes through morphospace and have nearly identical variance-covariance structures among morphological elements.

Conclusions/Significance

By demonstrating that replicate adaptive radiations are diverging along parallel axes, we have shown that non-convergence in the common morphospace is associated with convergence in the global morphospace. Applying these complimentary analyses to future comparative studies will improve our understanding of the relationship between morphological convergence and non-convergence, and the roles of contingency, chance and determinism in driving morphological diversification.

Mosaic retroposon insertion patterns in placental mammals

One and a half centuries after Charles Darwin and Alfred Russel Wallace outlined our current understanding of evolution, a new scientific era is dawning that enables direct observations of genetic variation. However, pure sequence-based molecular attempts to resolve the basal origin of placental mammals have so far resulted only in apparently conflicting hypotheses. By contrast, in the mammalian genomes where they were highly active, the insertion of retroelements and their comparative insertion patterns constitute a neutral, virtually homoplasy-free archive of evolutionary histories. The "presence" of a retroelement at an orthologous genomic position in two species indicates their common ancestry in contrast to its "absence" in more distant species. To resolve the placental origin controversy we extracted approximately 2 million potentially phylogenetically informative, retroposon-containing loci from representatives of the major placental mammalian lineages and found highly significant evidence challenging all current single hypotheses of their basal origin. The Exafroplacentalia hypothesis (Afrotheria as the sister group to all remaining placentals) is significantly supported by five retroposon insertions, the Epitheria hypothesis (Xenarthra as the sister group to all remaining placentals) by nine insertion patterns, and the Atlantogenata hypothesis (a monophyletic clade comprising Xenarthra and Afrotheria as the sister group to Boreotheria comprising all remaining placentals) by eight insertion patterns. These findings provide significant support for a "soft" polytomy of the major mammalian clades. Ancestral successive hybridization events and/or incomplete lineage sorting associated with short speciation intervals are viable explanations for the mosaic retroposon insertion patterns of recent placental mammals and for the futile search for a clear root dichotomy.

Support patterns from different outgroups provide a strong phylogenetic signal

It is known that the accuracy of phylogenetic reconstruction decreases when more distant outgroups are used. We quantify this phenomenon with a novel scoring method, the outgroup score pOG. This score expresses if the support for a particular branch of a tree decreases with increasingly distant outgroups. Large-scale simulations confirmed that the outgroup support follows this expectation and that the pOG score captures this pattern. The score often identifies the correct topology even when the primary reconstruction methods fail, particularly in the presence of model violations. In simulations of problematic phylogenetic scenarios such as rate variation among lineages (which can lead to long-branch attraction artifacts) and quartet-based reconstruction, the pOG analysis outperformed the primary reconstruction methods. Because the pOG method does not make any assumptions about the evolutionary model (besides the decreasing support from increasingly distant outgroups), it can detect cases of violations not treated by a specific model or too strong to be fully corrected. When used as an optimization criterion in the construction of a tree of 23 mammals, the outgroup signal confirmed many well-accepted mammalian orders and superorders. It supports Atlantogenata, a clade of Afrotheria and Xenarthra, and suggests an Artiodactyla-Chiroptera clade.

Resolving arthropod phylogeny: exploring phylogenetic signal within 41 kb of protein-coding nuclear gene sequence

This study attempts to resolve relationships among and within the four basal arthropod lineages (Pancrustacea, Myriapoda, Euchelicerata, Pycnogonida) and to assess the widespread expectation that remaining phylogenetic problems will yield to increasing amounts of sequence data. Sixty-eight regions of 62 protein-coding nuclear genes (approximately 41 kilobases (kb)/taxon) were sequenced for 12 taxonomically diverse arthropod taxa and a tardigrade outgroup. Parsimony, likelihood, and Bayesian analyses of total nucleotide data generally strongly supported the monophyly of each of the basal lineages represented by more than one species. Other relationships within the Arthropoda were also supported, with support levels depending on method of analysis and inclusion/exclusion of synonymous changes. Removing third codon positions, where the assumption of base compositional homogeneity was rejected, altered the results. Removing the final class of synonymous mutations--first codon positions encoding leucine and arginine, which were also compositionally heterogeneous--yielded a data set that was consistent with a hypothesis of base compositional homogeneity. Furthermore, under such a data-exclusion regime, all 68 gene regions individually were consistent with base compositional homogeneity. Restricting likelihood analyses to nonsynonymous change recovered trees with strong support for the basal lineages but not for other groups that were variably supported with more inclusive data sets. In a further effort to increase phylogenetic signal, three types of data exploration were undertaken. (1) Individual genes were ranked by their average rate of nonsynonymous change, and three rate categories were assigned--fast, intermediate, and slow. Then, bootstrap analysis of each gene was performed separately to see which taxonomic groups received strong support. Five taxonomic groups were strongly supported independently by two or more genes, and these genes mostly belonged to the slow or intermediate categories, whereas groups supported only by a single gene region tended to be from genes of the fast category, arguing that fast genes provide a less consistent signal. (2) A sensitivity analysis was performed in which increasing numbers of genes were excluded, beginning with the fastest. The number of strongly supported nodes increased up to a point and then decreased slightly. Recovery of Hexapoda required removal of fast genes. Support for Mandibulata (Pancrustacea + Myriapoda) also increased, at times to "strong" levels, with removal of the fastest genes. (3) Concordance selection was evaluated by clustering genes according to their ability to recover Pancrustacea, Euchelicerata, or Myriapoda and analyzing the three clusters separately. All clusters of genes recovered the three concordance clades but were at times inconsistent in the relationships recovered among and within these clades, a result that indicates that the a priori concordance criteria may bias phylogenetic signal in unexpected ways. In a further attempt to increase support of taxonomic relationships, sequence data from 49 additional taxa for three slow genes (i.e., EF-1 alpha, EF-2, and Pol II) were combined with the various 13-taxon data sets. The 62-taxon analyses supported the results of the 13-taxon analyses and provided increased support for additional pancrustacean clades found in an earlier analysis including only EF-1 alpha, EF-2, and Pol II.

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.

Amniote phylogenomics: testing evolutionary hypotheses with BAC library scanning and targeted clone analysis of large-scale DNA sequences from reptiles

Phylogenomics research integrating established principles of systematic biology and taking advantage of the wealth of DNA sequences being generated by genome science holds promise for answering long-standing evolutionary questions with orders of magnitude more primary data than in the past. Although it is unrealistic to expect whole-genome initiatives to proceed rapidly for commercially unimportant species such as reptiles, practical approaches utilizing genomic libraries of large-insert clones pave the way for a phylogenomics of species that are nevertheless essential for testing evolutionary hypotheses within a phylogenetic framework. This chapter reviews the case for adopting genome-enabled approaches to evolutionary studies and outlines a program for using bacterial artificial chromosome (BAC) libraries or plasmid libraries as a basis for completing "genome scans" of reptiles. We have used BACs to close a critical gap in the genome database for Reptilia, the sister group of mammals, and present the methodological approaches taken to achieve this as a guideline for designing similar comparative studies. In addition, we provide a detailed step-by-step protocol for BAC-library screening and shotgun sequencing of specific clones containing target genes of evolutionary interest. Taken together, the genome scanning and shotgun sequencing techniques offer complementary diagnostic potential and can substantially increase the scale and power of analyses aimed at testing evolutionary hypotheses for nonmodel species.

Affect cues in vocalizations of the bat, Megaderma lyra, during agonistic interactions

Some features of emotional prosody in human speech may be traced back to affect cues in mammalian vocalizations. The present study addresses the question whether affect intensity, as expressed by the intensity of behavioral displays, is encoded in vocal cues, i.e., changes in the structure of associated calls, in bats, a group evolutionarily remote from primates. A frame-by-frame video analysis of 109 dyadic agonistic interactions recorded in approach situations was performed to categorize displays into two intensity levels based on a cost-benefit estimate. M. lyra showed graded visual displays accompanied by specific calls and response calls of the second bat. A sound analysis revealed systematic changes of call sequence parameters with display level. At the high intensity level, total call duration, number of syllables within a call, and the number of calls within a sequence were increased, while intervals between call syllables were decreased for both call types. In addition, the latency of the response call was shorter, and its main syllable-type durations and fundamental frequency were increased. These systematic changes of vocal parameters with affect intensity correspond to prosodic changes in human speech, suggesting that emotion-related acoustic cues are a common feature of vocal communication in mammals.

Rooting the eutherian tree: the power and pitfalls of phylogenomics

In an attempt to root the eutherian tree using genome-scale data with the maximum likelihood method, a concatenate analysis supports a putatively wrong tree, whereas separate analyses of different genes reduced the bias.

Background

Ongoing genome sequencing projects have led to a phylogenetic approach based on genome-scale data (phylogenomics), which is beginning to shed light on longstanding unresolved phylogenetic issues. The use of large datasets in phylogenomic analysis results in a global increase in resolution due to a decrease in sampling error. However, a fully resolved tree can still be wrong if the phylogenetic inference is biased.

Results

Here, in an attempt to root the eutherian tree using genome-scale data with the maximum likelihood method, we demonstrate a case in which a concatenate analysis strongly supports a putatively wrong tree, whereas the total evaluation of separate analyses of different genes grossly reduced the bias of the phylogenetic inference. A conventional method of concatenate analysis of nucleotide sequences from our dataset, which includes a more than 1 megabase alignment of 2,789 nuclear genes, suggests a misled monophyly of Afrotheria (for example, elephant) and Xenarthra (for example, armadillo) with 100% bootstrap probability. However, this tree is not supported by our 'separate method', which takes into account the different tempos and modes of evolution among genes, and instead the basal Afrotheria tree is favored.

Conclusion

Our analysis demonstrates that in cases in which there is great variation in evolutionary features among different genes, the separate model, rather than the concatenate model, should be used for phylogenetic inference, especially in genome-scale data.

A quantitative analysis of the Eutherian orbit: correlations with masticatory apparatus

The mammalian orbit, or eye-socket, is a highly plastic region of the skull. It comprises between seven and nine bones, all of which vary widely in their contribution to this region among the different mammalian orders and families. It is hypothesised that the structure of the mammalian orbit is principally influenced by the forces generated by the jaw-closing musculature. In order to quantify the orbit, fourteen linear, angular and area measurements were taken from 84 species of placental mammals using a Microscribe-3D digitiser. The results were then analysed using principal components analysis. The results of the multivariate analysis on untransformed data showed a clear division of the mammalian taxa into temporalis-dominant forms and masseter-dominant forms. This correlation between orbital structure and masticatory musculature was reinforced by results from the size-corrected data, which showed a separation of the taxa into the three specialised feeding types proposed by Turnbull (1970): i.e. 'carnivore-shear', 'ungulate-grinding' and 'rodent-gnawing'. Moreover, within the rodents there was a clear distinction between species in which the masseter is highly developed and those in which the temporalis has more prominence. These results were reinforced by analysis of variance which showed significant differences in the relative orbital areas of certain bones between temporalis-dominant and masseter-dominant taxa. Subsequent cluster analysis suggested that most of the variables could be grouped into three assemblages: those associated with the length of the rostrum; those associated with the width of the skull; and those associated with the relative size of the orbit and the shape of the face. However, the relative area of the palatine bone showed weak correlations with the other variables and did not fit into any group. Overall the relative area of the palatine was most closely correlated with feeding type, and this measure that appeared to be most strongly associated with the arrangement of the masticatory musculature. These results give a strong indication that, although orbital structure is in part determined by the relative size and orientation of the orbits, the forces generated by the muscles of mastication also have a large effect.

Imprinted expression of the canine IGF2R, in the absence of an anti-sense transcript or promoter methylation

Imprinted genes are epigenetically modified in a parent of origin-dependent manner, and as a consequence, are differentially expressed. Although the evolution of genomic imprinting is a subject of intense debate, imprinted genes have been studied primarily in mice and humans and in a small number of marsupial mammals. Comparative studies involving rodents and primates are of limited value, as they belong to the same superordinal group of eutherian mammals (Euarchontoglires). On the other hand, comparisons involving marsupials may not be informative, due to phylogenetic distance. Canis familiaris belongs to Laurasiatheria, a sister-group of Euarchontoglires, and should prove useful in comparative studies of imprinted genes. Using RT-PCR we demonstrate monoallelic expression of the canine IGF2R in several tissues, including uterus and umbilical cord. In the case of umbilical cord, we identify the expressed allele as maternally derived. The canine IGF2R is thus an imprinted gene. Using bisulfite sequencing, we show that the canine IGF2R resembles the imprinted mouse Igf2r in having a CpG island in intron 2 that is hemi-methylated. However, it differs from the mouse gene in that maintenance of the monoallelic expression of canine IGF2R does not require expression of an anti-sense transcript from the paternally derived allele, or methylation of the repressed IGF2R promoter. In these two important features, the imprinted canine gene resembles the imprinted opossum IGF2R. Our data suggest that these features were properties of the ancestral imprinted IGF2R and that the anti-sense transcript (Air) and promoter methylation observed in mouse are derived features of the mouse Igf2r locus.

Computer simulation of feeding behaviour in the thylacine and dingo as a novel test for convergence and niche overlap

The extinct marsupial thylacine (Thylacinus cynocephalus) and placental grey wolf (Canis lupus) are commonly presented as an iconic example of convergence. However, various analyses suggest distinctly different behaviours and specialization towards either relatively small or large prey in the thylacine, bringing the degree of apparent convergence into question. Here we apply a powerful engineering tool, three-dimensional finite element analysis incorporating multiple material properties for bone, to examine mechanical similarity and niche overlap in the thylacine and the wolf subspecies implicated in its extinction from mainland Australia, Canis lupus dingo. Comparisons of stress distributions not only reveal considerable similarity, but also informative differences. The thylacine's mandible performs relatively poorly where only the actions of the jaw muscles are considered, although this must be considered in the light of relatively high bite forces. Stresses are high in the posterior of the thylacine's cranium under loads that simulate struggling prey. We conclude that relative prey size may have been comparable where both species acted as solitary predators, but that the dingo is better adapted to withstand the high extrinsic loads likely to accompany social hunting of relatively large prey. It is probable that there was considerable ecological overlap. As a large mammalian hypercarnivore adapted to taking small-medium sized prey, the thylacine may have been particularly vulnerable to disturbance.

The isolator piglet: a model for studying the development of adaptive immunity

The period from late gestation to weaning in neonatal mammals is a critical window when the adaptive immune system develops and replaces the protection temporarily provided by passive immunity and pre-adaptive antibodies. It is also when oral tolerance to dietary antigen and the distinction between commensal and pathogenic gut bacteria becomes established resulting in immune homeostasis. The reproductive biology of swine provides a unique model for distinguishing the effects of different factors on immune development during this critical period because all extrinsic factors are controlled by the experimenter. This chapter reviews this early stage of development and the usefulness of the piglet model for understanding events during this transitional stage. The review also describes the major features of the porcine immune system and the immune stimulatory and dysregulatory factors that act during this period. The value of the model to medical science in such areas as food allergy, organ transplantation, cystic fibrosis and the production of humanized antibodies for immuno-therapy is discussed.

Recruitment of mitochondrial tRNA genes as auxiliary variability markers for both intra- and inter-species analysis: The paradigm of brown hare (Lepus europaeus)

We sequenced and analyzed the mitochondrial tRNA(Thr) and tRNA(Pro) genes from brown hare (Lepus europaeus) individuals of different geographic distribution and we investigated the role of various nucleotide substitutions that were detected. We compared these tRNAs with the respective available mitochondrial tRNA genes sequences within Lepus species and among mammals. The mutations that were detected represent specific and conserved polymorphisms that do not seem to affect the structural and functional features that are required for participation of tRNA molecules in mitochondrial protein synthesis. These changes however, possibly reflect on the evolutionary background of the species, which is based on the high intra-genomic variability and the evolutionary dynamic of the mitochondrial DNA. In an attempt to compare the phylogeny that is based on these specific tRNA genes with the phylogeny that is produced from sequencing data of the mitochondrial variable loop, we came up with results that indicate similar phylogeographic clusters. This observation implies that the tRNA mutations that were used for the present study have been well tolerated during evolution and they define an additional genetic and biochemical tag that can be used for such studies. Based on this notion and according to our results, we propose that mitochondrial tRNA genes can be used as valuable auxiliary molecular markers for contemporaneous and linked biochemical and genetic analyses.

Evolution of the placenta during the early radiation of placental mammals

The chorioallantoic placenta is an organ of gaseous exchange that exhibits a high degree of structural diversity. One factor determining oxygen transfer across the placenta, the diffusion distance, is in part dependent on the number of cell layers separating maternal from fetal blood. This interhaemal barrier occurs in three principal variants. The focus of this review is on determining how the barrier evolved in placental mammals. The analysis was based on current knowledge of placental structure, as far as possible using ultrastructural data, and on current views about the evolution of placental mammals, derived from molecular phylogenetics. We show that epitheliochorial placentation, the least invasive type, is a derived state and discuss factors that may have determined its evolution with reference to conflict theory, as applied to the allocation of resources between mother and fetus. It is not yet possible to determine which of the two more invasive types of placentation occurred in the last common ancestor of crown placentals. Depending on tree topology and taxon sampling, the result achieved is either endotheliochorial, haemochorial or unresolved. Finally we discuss other factors important to placental gas exchange and point to physiological variables that might become amenable to phylogenetic analysis.

OrthoMaM: a database of orthologous genomic markers for placental mammal phylogenetics

Background

Molecular sequence data have become the standard in modern day phylogenetics. In particular, several long-standing questions of mammalian evolutionary history have been recently resolved thanks to the use of molecular characters. Yet, most studies have focused on only a handful of standard markers. The availability of an ever increasing number of whole genome sequences is a golden mine for modern systematics. Genomic data now provide the opportunity to select new markers that are potentially relevant for further resolving branches of the mammalian phylogenetic tree at various taxonomic levels.

Description

The EnsEMBL database was used to determine a set of orthologous genes from 12 available complete mammalian genomes. As targets for possible amplification and sequencing in additional taxa, more than 3,000 exons of length > 400 bp have been selected, among which 118, 368, 608, and 674 are respectively retrieved for 12, 11, 10, and 9 species. A bioinformatic pipeline has been developed to provide evolutionary descriptors for these candidate markers in order to assess their potential phylogenetic utility. The resulting OrthoMaM (Orthologous Mammalian Markers) database can be queried and alignments can be downloaded through a dedicated web interface .

Conclusion

The importance of marker choice in phylogenetic studies has long been stressed. Our database centered on complete genome information now makes possible to select promising markers to a given phylogenetic question or a systematic framework by querying a number of evolutionary descriptors. The usefulness of the database is illustrated with two biological examples. First, two potentially useful markers were identified for rodent systematics based on relevant evolutionary parameters and sequenced in additional species. Second, a complete, gapless 94 kb supermatrix of 118 orthologous exons was assembled for 12 mammals. Phylogenetic analyses using probabilistic methods unambiguously supported the new placental phylogeny by retrieving the monophyly of Glires, Euarchontoglires, Laurasiatheria, and Boreoeutheria. Muroid rodents thus do not represent a basal placental lineage as it was mistakenly reasserted in some recent phylogenomic analyses based on fewer taxa. We expect the OrthoMaM database to be useful for further resolving the phylogenetic tree of placental mammals and for better understanding the evolutionary dynamics of their genomes, i.e., the forces that shaped coding sequences in terms of selective constraints.

Molecular and genomic data identify the closest living relative of primates

A full understanding of primate morphological and genomic evolution requires the identification of their closest living relative. In order to resolve the ancestral relationships among primates and their closest relatives, we searched multispecies genome alignments for phylogenetically informative rare genomic changes within the superordinal group Euarchonta, which includes the orders Primates, Dermoptera (colugos), and Scandentia (treeshrews). We also constructed phylogenetic trees from 14 kilobases of nuclear genes for representatives from most major primate lineages, both extant colugos, and multiple treeshrews, including the pentail treeshrew, Ptilocercus lowii, the only living member of the family Ptilocercidae. A relaxed molecular clock analysis including Ptilocercus suggests that treeshrews arose approximately 63 million years ago. Our data show that colugos are the closest living relatives of primates and indicate that their divergence occurred in the Cretaceous.

A new estimate of afrotherian phylogeny based on simultaneous analysis of genomic, morphological, and fossil evidence

BACKGROUND

The placental mammalian clade Afrotheria is now supported by diverse forms of genomic data, but interordinal relationships within, and morphological support for, the group remains elusive. As a means for addressing these outstanding problems, competing hypotheses of afrotherian interordinal relationships were tested through simultaneous parsimony analysis of a large data set (> 4,590 parsimony informative characters) containing genomic data (> 17 kb of nucleotide data, chromosomal associations, and retroposons) and 400 morphological characters scored across 16 extant and 35 extinct afrotherians.

RESULTS

Parsimony analysis of extant taxa alone recovered the interordinal topology (Afrosoricida, ((Macroscelidea, Tubulidentata), (Hyracoidea, (Proboscidea, Sirenia)))). Analysis following addition of extinct taxa instead supported Afroinsectivora (Afrosoricida + Macroscelidea) and Pseudoungulata (Tubulidentata + Paenungulata), as well as Tethytheria (Proboscidea + Sirenia). This latter topology is, however, sensitive to taxon deletion and different placements of the placental root, and numerous alternative interordinal arrangements within Afrotheria could not be statistically rejected. Relationships among extinct stem members of each afrotherian clade were more stable, but one alleged stem macroscelidean (Herodotius) never grouped with that clade and instead consistently joined pseudoungulates or paenungulates. When character transformations were optimized onto a less resolved afrotherian tree that reflects uncertainty about the group's interordinal phylogeny, a total of 21 morphological features were identified as possible synapomorphies of crown Afrotheria, 9 of which optimized unambiguously across all character treatments and optimization methods.

CONCLUSION

Instability in afrotherian interordinal phylogeny presumably reflects rapid divergences during two pulses of cladogenesis - the first in the Late Cretaceous, at and just after the origin of crown Afrotheria, and the second in the early Cenozoic, with the origin of crown Paenungulata. Morphological evidence for divergences during these two pulses either never existed or has largely been "erased" by subsequent evolution along long ordinal branches. There may, nevertheless, be more morphological character support for crown Afrotheria than is currently assumed; the features identified here as possible afrotherian synapomorphies can be further scrutinized through future phylogenetic analyses with broader taxon sampling, as well as recovery of primitive fossil afrotherians from the Afro-Arabian landmass, where the group is likely to have first diversified.

Chromosomal instability in Afrotheria: fragile sites, evolutionary breakpoints and phylogenetic inference from genome sequence assemblies

BACKGROUND

Extant placental mammals are divided into four major clades (Laurasiatheria, Supraprimates, Xenarthra and Afrotheria). Given that Afrotheria is generally thought to root the eutherian tree in phylogenetic analysis of large nuclear gene data sets, the study of the organization of the genomes of afrotherian species provides new insights into the dynamics of mammalian chromosomal evolution. Here we test if there are chromosomal bands with a high tendency to break and reorganize in Afrotheria, and by analyzing the expression of aphidicolin-induced common fragile sites in three afrotherian species, whether these are coincidental with recognized evolutionary breakpoints.

RESULTS

We described 29 fragile sites in the aardvark (OAF) genome, 27 in the golden mole (CAS), and 35 in the elephant-shrew (EED) genome. We show that fragile sites are conserved among afrotherian species and these are correlated with evolutionary breakpoints when compared to the human (HSA) genome. Inddition, by computationally scanning the newly released opossum (Monodelphis domestica) and chicken sequence assemblies for use as outgroups to Placentalia, we validate the HSA 3/21/5 chromosomal synteny as a rare genomic change that defines the monophyly of this ancient African clade of mammals. On the other hand, support for HSA 1/19p, which is also thought to underpin Afrotheria, is currently ambiguous.

CONCLUSION

We provide evidence that (i) the evolutionary breakpoints that characterise human syntenies detected in the basal Afrotheria correspond at the chromosomal band level with fragile sites, (ii) that HSA 3p/21 was in the amniote ancestor (i.e., common to turtles, lepidosaurs, crocodilians, birds and mammals) and was subsequently disrupted in the lineage leading to marsupials. Its expansion to include HSA 5 in Afrotheria is unique and (iii) that its fragmentation to HSA 3p/21 + HSA 5/21 in elephant and manatee was due to a fission within HSA 21 that is probably shared by all Paenungulata.

Species-specific evolution of immune receptor tyrosine based activation motif-containing CEACAM1-related immune receptors in the dog

BACKGROUND

Although the impact of pathogens on the evolution of the mammalian immune system is still under debate, proteins, which both regulate immune responses and serve as cellular receptors for pathogens should be at the forefront of pathogen-driven host evolution. The CEA (carcinoembryonic antigen) gene family codes for such proteins and indeed shows tremendous species-specific variation between human and rodents. Since little is known about the CEA gene family in other lineages of placental mammals, we expected to gain new insights into the evolution of the rapidly diverging CEA family by analyzing the CEA family of the dog.

RESULTS

Here we describe the complete CEA gene family in the dog. We found that the gene coding for the ITIM-bearing immunoregulatory molecule CEACAM1 gave rise to a recent expansion of the canine CEA gene family by gene duplication, similar to that previously found in humans and mice. However, while the murine and human CEACAMs (carcinoembryonic antigen-related cell adhesion molecules) are predominantly secreted and GPI-anchored, respectively, in the dog, most of the CEACAMs represent ITAM-bearing transmembrane proteins. One of these proteins, CEACAM28, exhibits nearly complete sequence identity with the ligand-binding N domain of CEACAM1, but antagonizing signaling motifs in the cytoplasmic tail. Comparison of nonsynonymous and synonymous substitutions indicates that the CEACAM28 N domain is under the strongest purifying selection of all canine CEACAM1-related CEACAMs. In addition, CEACAM28 shows a similar expression pattern in resting immune cells and tissues as CEACAM1. However, upon activation CEACAM28 mRNA and CEACAM1 mRNA are differentially regulated.

CONCLUSION

Thus, CEACAM1 and CEACAM28 are the first paired immune receptors identified within the CEA gene family, which are expressed on T cells and are most likely involved in the fine-tuning of T cell responses. The direction of gene conversion accompanied by purifying selection and expression in immune cells suggests the possibility that CEACAM28 evolved in response to selective pressure imposed by species-specific pathogens.

THE LIKELIHOOD NODE DENSITY EFFECT AND CONSEQUENCES FOR EVOLUTIONARY STUDIES OF MOLECULAR RATES

Many molecular phylogenies show longer root-to-tip path lengths in species-rich groups, encouraging hypotheses linking cladogenesis with accelerated molecular evolution. However, the pattern can also be caused by an artifact called the node density effect (NDE): this effect occurs when the method used to reconstruct a tree underestimates multiple hits that would have been revealed by extra nodes, leading to longer root-to-tip path lengths in clades with more terminal taxa. Here we use a twofold approach to demonstrate that maximum likelihood and Bayesian methods also suffer from the NDE known to affect parsimony. First, simulations deliberately mismatching the simulation and reconstruction models show that the greater the model disparity, the greater the gap between actual and reconstructed tree lengths, and the greater the NDE. Second, taxon sampling manipulation with empirical data shows that NDE can still be present when using optimized models: across 12 datasets, 70 out of 109 sister path comparisons showed significant evidence of NDE. Unless the model fairly accurately reconstructs the real tree length-and given the complexity of real sequence evolution this may be uncommon -- it will consistently produce a node density artifact. At commonly encountered divergence levels, a 10% underestimation of tree length results in > or = 80% of simulated phylogenies showing a positive NDE. Bayesian trees have a slight but consistently stronger effect. This pervasive methodological artifact increases apparent rate heterogeneity, and can compromise investigations of factors influencing molecular evolutionary rate that use path lengths in topologically asymmetric trees.

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.

Detecting and overcoming systematic errors in genome-scale phylogenies

Genome-scale data sets result in an enhanced resolution of the phylogenetic inference by reducing stochastic errors. However, there is also an increase of systematic errors due to model violations, which can lead to erroneous phylogenies. Here, we explore the impact of systematic errors on the resolution of the eukaryotic phylogeny using a data set of 143 nuclear-encoded proteins from 37 species. The initial observation was that, despite the impressive amount of data, some branches had no significant statistical support. To demonstrate that this lack of resolution is due to a mutual annihilation of phylogenetic and nonphylogenetic signals, we created a series of data sets with slightly different taxon sampling. As expected, these data sets yielded strongly supported but mutually exclusive trees, thus confirming the presence of conflicting phylogenetic and nonphylogenetic signals in the original data set. To decide on the correct tree, we applied several methods expected to reduce the impact of some kinds of systematic error. Briefly, we show that (i) removing fast-evolving positions, (ii) recoding amino acids into functional categories, and (iii) using a site-heterogeneous mixture model (CAT) are three effective means of increasing the ratio of phylogenetic to nonphylogenetic signal. Finally, our results allow us to formulate guidelines for detecting and overcoming phylogenetic artefacts in genome-scale phylogenetic analyses.

Primate Origins: Implications of a Cretaceous Ancestry

It has long been accepted that the adaptive radiation of modern placental mammals, like that of modern birds, did not begin until after the Cretaceous/Tertiary (K/T) boundary 65 million years (Ma) ago, following the extinction of the dinosaurs. The first undoubted fossil relatives of modern primates appear in the record 55 Ma ago. However, in agreement with evidence from molecular phylogenies calibrated with dates from denser parts of the fossil record, a statistical analysis of the primate record allowing for major gaps now indicates a Cretaceous origin of euprimates 80-90 Ma ago. If this interpretation is correct, primates overlapped with dinosaurs by some 20 Ma prior to the K/T boundary, and the initial radiation of primates was probably truncated as part of the major extinction event that occurred at the end of the Cretaceous. Following a review of evidence for an early origin of primates, implications of this are discussed with respect to the likely ancestral condition for primates, including a southern continental area of origin and moderately large body size. The known early Tertiary primates are re-interpreted as northern continental offshoots of a 'second wave' of primate evolution.

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.

Early Tertiary mammals from North Africa reinforce the molecular Afrotheria clade

The phylogenetic pattern and timing of the radiation of mammals, especially the geographical origins of major crown clades, are areas of controversy among molecular biologists, morphologists and palaeontologists. Molecular phylogeneticists have identified an Afrotheria clade, which includes several taxa as different as tenrecs (Tenrecidae), golden moles (Chrysochloridae), elephant-shrews (Macroscelididae), aardvarks (Tubulidentata) and paenungulates (elephants, sea cows and hyracoids). Molecular data also suggest a Cretaceous African origin for Afrotheria within Placentalia followed by a long period of endemic evolution on the Afro-Arabian continent after the mid-Cretaceous Gondwanan breakup (approx. 105-25 Myr ago). However, there was no morphological support for such a natural grouping so far. Here, we report new dental and postcranial evidence of Eocene stem hyrax and macroscelidid from North Africa that, for the first time, provides a congruent phylogenetic view with the molecular Afrotheria clade. These new fossils imply, however, substantial changes regarding the historical biogeography of afrotheres. Their long period of isolation in Africa, as assumed by molecular inferences, is now to be reconsidered inasmuch as Eocene paenungulates and elephant-shrews are here found to be related to some Early Tertiary Euramerican 'hyopsodontid condylarths' (archaic hoofed mammals). As a result, stem members of afrotherian clades are not strictly African but also include some Early Paleogene Holarctic mammals.

Congruence of molecules and morphology using a narrow allometric approach

There are many cases of incongruence between phylogenetic hypotheses produced from morphological data and those produced from molecular data. In such instances, many researchers prefer to accept the results of molecular phylogenies. For example, in a recent analysis of primate phylogenies based on craniodental characters, Collard and Wood [Collard M, Wood BA (2000) Proc Natl Acad Sci USA 97:5003-5006] argued that, because craniodental data do not yield relationships concordant with molecular studies, the results of studies that employ such characters must be considered suspect. As most of our knowledge of mammalian evolution and phylogenetic history comes from craniodental fossils, these results have dramatic implications. However, the aforementioned analysis did not take into account the potentially confounding effects of allometry on quantitative craniodental characters. In this article, we employ a previously undescribed narrow allometric coding method that accounts for such confounding influences in phylogenetic analyses of craniodental morphology. By using essentially the same raw data set as Collard and Wood [Collard M, Wood BA (2000) Proc Natl Acad Sci USA 97:5003-5006], 65 quantitative craniodental characters were adjusted in a parsimony analysis for the primate tribe Papionini, a group of monkeys argued to display extensive homoplasy. The resulting phylogenetic tree was congruent with the phylogenetic tree based on molecular data for these species, thereby meeting the "criterion of congruence." These results suggest that morphological data, when treated properly, can be considered as reliable as molecular data in phylogenetic reconstruction. Rather than accepting phylogenetic hypotheses from one data source over another, cases of incongruence should be examined with greater scrutiny.

Accounting for gene rate heterogeneity in phylogenetic inference

Traditionally, phylogenetic analyses over many genes combine data into a contiguous block. Under this concatenated model, all genes are assumed to evolve at the same rate. However, it is clear that genes evolve at very different rates and that accounting for this rate heterogeneity is important if we are to accurately infer phylogenies from heterogeneous multigene data sets. There remain open questions regarding how best to incorporate gene rate parameters into phylogenetic models and which properties of real data correlate with improved fit over the concatenated model. In this study, two methods of accounting for gene rate heterogeneity are compared: the n-parameter method, which allows for each of the n gene partitions to have a gene rate parameter, and the alpha-parameter method, which fits a distribution to the gene rates. Results demonstrate that the n-parameter method is both computationally faster and in general provides a better fit over the concatenated model than the alpha-parameter method. Furthermore, improved model fit over the concatenated model is highly correlated with the presence of a gene with a slow relative rate of evolution.

The Origin and Evolution of Enamel Mineralization Genes

BACKGROUND/AIMS

Enamel and enameloid were identified in early jawless vertebrates, about 500 million years ago (MYA). This suggests that enamel matrix proteins (EMPs) have at least the same age. We review the current data on the origin, evolution and relationships of enamel mineralization genes.

METHODS AND RESULTS

Three EMPs are secreted by ameloblasts during enamel formation: amelogenin (AMEL), ameloblastin (AMBN) and enamelin (ENAM). Recently, two new genes, amelotin (AMTN) and odontogenic ameloblast associated (ODAM), were found to be expressed by ameloblasts during maturation, increasing the group of ameloblast-secreted proteins to five members. The evolutionary analysis of these five genes indicates that they are related: AMEL is derived from AMBN, AMTN and ODAM are sister genes, and all are derived from ENAM. Using molecular dating, we showed that AMBN/AMEL duplication occurred >600 MYA. The large sequence dataset available for mammals and reptiles was used to study AMEL evolution. In the N- and C-terminal regions, numerous residues were unchanged during >200 million years, suggesting that they are important for the proper function of the protein.

CONCLUSION

The evolutionary analysis of AMEL led to propose a dataset that will be useful to validate AMEL mutations leading to X- linked AI.

CONVERGENCE AND REMARKABLY CONSISTENT CONSTRAINT IN THE EVOLUTION OF CARNIVORE SKULL SHAPE

Phenotypic similarities between distantly related marsupials and placentals are commonly presented as examples of convergence and support for the role of adaptive evolution in shaping morphological and ecological diversity. Here we compare skull shape in a wide range of carnivoran placentals (Carnivora) and nonherbivorous marsupials using a three-dimensional (3-D) geometric morphometric approach. Morphological and ecological diversity among extant carnivorans is considerably greater than is evident in the marsupial order Dasyuromorphia with which they have most commonly been compared. To examine convergence across a wider, but broadly comparable range of feeding ecologies, a dataset inclusive of nondasyuromorphian marsupials and extinct taxa representing morphotypes no longer present was assembled. We found support for the adaptive paradigm, with correlations between morphology, feeding behavior, and bite force, although skull shape better predicted feeding ecology in the phylogenetically diverse marsupial sample than in carnivorans. However, we also show that remarkably consistent but differing constraints have influenced the evolution of cranial shape in both groups. These differences between carnivorans and marsupials, which correlate with brain size and bite force, are maintained across the full gamut of morphologies and feeding categories, from small insectivores and omnivores to large meat-specialists.

Validation of amelogenesis imperfecta inferred from amelogenin evolution

We used the evolutionary analysis of amelogenin (AMEL) in 80 amniotes (52 mammalian and 28 reptilian sequences) to aid in the genetic diagnosis of X-linked amelogenesis imperfecta (AIH1). Out of 191 residues, 77 were found to be unchanged in mammals, and only 34 in amniotes. The latter are considered crucial residues for enamel formation, while the 43 residues conserved only in mammals could indicate that they play new, important roles for enamel formation in this lineage. The 5 substitutions leading to AIH1 were validated when the mammalian dataset was used, and 4 of them with the amniote dataset. These 2 sequence datasets will facilitate the validation of any human AMEL mutation suspected of involvement in AIH1. This evolutionary analysis also revealed numerous residues that appeared to be important for correct AMEL function, but their role remains to be elucidated.

Using genomic data to unravel the root of the placental mammal phylogeny

The phylogeny of placental mammals is a critical framework for choosing future genome sequencing targets and for resolving the ancestral mammalian genome at the nucleotide level. Despite considerable recent progress defining superordinal relationships, several branches remain poorly resolved, including the root of the placental tree. Here we analyzed the genome sequence assemblies of human, armadillo, elephant, and opossum to identify informative coding indels that would serve as rare genomic changes to infer early events in placental mammal phylogeny. We also expanded our species sampling by including sequence data from >30 ongoing genome projects, followed by PCR and sequencing validation of each indel in additional taxa. Our data provide support for a sister-group relationship between Afrotheria and Xenarthra (the Atlantogenata hypothesis), which is in turn the sister-taxon to Boreoeutheria. We failed to recover any indels in support of a basal position for Xenarthra (Epitheria), which is suggested by morphology and a recent retroposon analysis, or a hypothesis with Afrotheria basal (Exafricoplacentalia), which is favored by phylogenetic analysis of large nuclear gene data sets. In addition, we identified two retroposon insertions that also support Atlantogenata and none for the alternative hypotheses. A revised molecular timescale based on these phylogenetic inferences suggests Afrotheria and Xenarthra diverged from other placental mammals approximately 103 (95-114) million years ago. We discuss the impacts of this topology on earlier phylogenetic reconstructions and repeat-based inferences of phylogeny.

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.

Members of a novel gene family, Gsdm, are expressed exclusively in the epithelium of the skin and gastrointestinal tract in a highly tissue-specific manner

Gasdermin (Gsdm) was originally identified as a candidate causative gene for several mouse skin mutants. Several Gsdm-related genes sharing a protein domain with DFNA5, the causative gene of human nonsyndromic hearing loss, have been found in the mouse and human genomes, and this group is referred to as the DFNA5-Gasdermin domain family. However, our current comparative genomic analysis identified several novel motifs distinct from the previously reported domain in the Gsdm-related genes. We also identified three new Gsdm genes clustered on mouse chromosome 15. We named these genes collectively the Gsdm family. Extensive expression analysis revealed exclusive expression of Gsdm family genes in the epithelium of the skin and gastrointestinal tract in a highly tissue-specific manner. Further database searching revealed the presence of other related genes with a similar N-terminal motif. These results suggest that the Gsdm family and related genes have evolved divergent epithelial expression profiles.