Distribution of human endogenous retrovirus HERV-K genomes in humans and different primates

The distribution of the human endogenous retrovirus (HERV)-K genome was investigated by Southern-blot analyses using a HERV-K-env DNA probe. With the exception of one DNA-sample, obtained from a Chinese individual in whom an amplification of HERV-K was detected, Southern-blot analyses yielded identical hybridization patterns with DNA from peripheral blood lymphocytes of 37 normal healthy blood donors, with DNA from six tumor cell lines, or with 23 DNA samples prepared from various carcinoma tissues. To elucidate whether the integration of HERV-K genomes into the primate lineage occurred as a single event or as an integration with later expansion, we further examined the evolutionary history of HERV-K by Southern blot analyses with DNA samples from different primate species. We detected HERV-K genomes in Macaca mulatta and Macaca silenus, which represent Old World monkeys, but not in prosimians (Galago demidovii) and New World monkeys, represented by Saguinus fuscicollis, Saguinus oedipus, and Calliihrix iacchus. Thus, we assume that the infection of the primate lineage with HERV-K had occurred after the divergence of New World and Old World monkeys, but before the evolutionary expansion of large hominoids. In contrast to the apparent lack of HERV-Kenv sequences in DNA from tissue of the New World monkey Saguinus oedipus (cotton-top marmoset), we found HERV-K-DNA in the B95-8 cell-line, which is a Saguinus oedipus leukocyte cell-line, immortalized in vitro by Epstein-Barr virus (EBV) and cultivated in human cells. It may be speculated that HERV-K-DNA or HERV-K-particles were introduced into these cells during in vitro transformation with EBV.     

The T-cell receptor β chain-encoding gene repertoire of a New World primate species,the cotton-top tamarin

 The New World primate, the cotton-top tamarin (Saguinus oedipus), expresses major histocompatibility complex (MHC) class I molecules with limited diversity. The uniqueness of the cotton-top tamarin MHC class I loci may contribute to this species’ unusual susceptibility to viral infections and high incidence of ulcerative colitis. As a prelude to examining the effect of this limited MHC class I diversity on the tamarin CD8⁺ T-cell receptor (TCR) repertoire, we identified expressed tamarin TCR β chain (TCRB) cDNAs by anchored and inverse polymerase chain reaction. Sequence alignments and phylogenetic comparisons with human and rhesus macaque sequences identified homologues of 21 human variable (V) gene families. Only single variable region genes were identified in each of these tamarin VB families, with the exception of the VB 5, 9, and 13 families which were comprised of two or three distinct members. The multiple genes within these three VB families do not appear to have separate human homologues, but rather aligned equally well to a single human gene from their respective VB families. These genes appear to have arisen, therefore, by duplication of certain VB genes in the tamarin ancestors following their divergence from the lineage leading to Old World primates and hominoids. Homologues of 12 of the 13 human joining (J) region genes were also identified in the tamarin. Comparison of the proportion of nonsynonymous (p_N) and synonymous (p_S) substitutions occurring per site within tamarin variable region genes demonstrated a reduction in p_N in the framework regions compared with p_N in the presumed MHC contact regions (CDR1 and CDR2). Taken together, these findings illustrate that the TCR β chain-encoding genes of the cotton-top tamarin are similar in structure and degree of complexity compared with their Old World primate and human counterparts. Received: 19 July 1996 / Revised: 12 August 1996     

Systematics of the Saguinus oedipus group of the bare-face tamarins: Evidence from facial morphology

The systematics of the Saguinus oedipus group within the bare-face tamarins remains open to question. Hershkovitz (Living New World Monkeys (Platyrrhini), Vol. 1. Chicago: University of Chicago Press, 1977) places the cotton-top and rufus-naped tamarins as subspecies of Saguinus oedipus (S. o. oedipus and S. o. geoffroyi, respectively). In contrast, several other authors have argued that these two taxa should be considered separate species (S. oedipus and S. geoffroyi). Phylogenetic relationships within the group are also disputed. Resolving these different interpretations has been difficult in part because no study of this group has included an objective measure of expected levels of specific vs. subspecific variation. We used facial measurements from 179 adult crania to address the systematics of this group and included a related species that is known to include multiple subspecies. Our sample included three taxa from the S. oedipus group of the bare-face tamarins (S. oedipus, S. geoffroyi, and S. leucopus) and six subspecies from the related hairy-face tamarin species S. fuscicollis. Comparisons to S. leucopus provided a relative measure of species-level differences. Analyses that included S. fuscicollis provided a measure of subspecific variation. There was no evidence of facial sexual dimorphism in any of these taxa. A variety of multivariate statistical analyses including discriminant function and cluster analysis suggest that S. oedipus and S. geoffroyi differ morphologically at a level consistent with species-level distinctions. The extent of differences between these taxa is large. The differences in their facial morphology was on the order of differences between S. oedipus or S. geoffroyi and S. leucopus rather than the extent of variation among S. fuscicollis subspecies. Furthermore, a comparison of collecting localities revealed that the variation we observed among S. oedipus and S. geoffroyi was not clinal but presented a large morphological discontinuity at the boundary between taxa. Our analyses also suggested that S. leucopus is more similar to S. oedipus than is either to S. geoffroyi. Finally, it may be that there are some distinct species within the S. fuscicollis group. However, this hypothesis, along with other phylogenetic relationships suggested by this study, will require more data and further study. © 1992 Wiley-Liss, Inc.     

C4 gene polymorphism in primates: evolution,generation, and Chido and Rodgers antigenicity

Eleven new C4d genomic primate sequences of the fourth complement factor (C4) have been obtained. Seven of them belong to five species not yet explored for this gene: Pan paniscus (pygmy chimpanzee), Cercopithecus aethiops (green monkey), Macaca mulatta (rhesus monkey), Macaca fascicularis (cynomolgus), and Saguinus oedipus (cotton top tamarin). The New World monkeys (tamarins, four individuals) sequenced for C4 have a single C4d sequence only, which shows a B isotypic specificity and a Rodgers 3 (Rg3), Chido 1 (Ch1) antigenicity. Rg3 and Ch1 could thus be the oldest Rg/Ch specificity (at least 50 million years old) and Rg1, Rg2, Ch3, and Ch6 could be more recent human-specific antigens. Mechanisms of C4d polymorphism generation were analyzed by compiling all the presently available sequences. Examples of both point mutations and crossing-over events among C4d primate sequences could be detected. The problem of a possible trans-species inheritance of C4d polymorphism was addressed and two apparently contradicting dendrograms were obtained. One of them, constructed by using both exon and intron sequences, does not support trans-species evolution, but supports the proposed theory of extensive homogenization of the C4 genes occurring within each species, because alleles from each primate species cluster together. Another completely different dendrogram, obtained by using exon sequences only, suggests the existence of trans-species evolution for C4d polymorphism, because alleles belonging to different species cluster together in a way similar to that found for HLA class I or II alleles. However, orangutan sequences group together in both kinds of C4d sequence dendrograms and seem to have arisen from an ancestor different from that of chimpanzee, gorilla and man C4d sequences. Finally, further data have been obtained that support trans-species conservation of A-ness and B-ness and the existence of trans-specifically conserved allelic motifs, both in intronic and exonic sequences.     

Molecular evolution of IgG subclass among nonhuman primates: implication of differences in antigenic determinants among Apes

The cross-reactivity of five different rabbit polyclonal antibodies to human IgG and IgG subclass (IgG1, IgG2, IgG3, and IgG4) was determined by competitive ELISA with nine nonhuman primate species including five apes, three Old World monkeys, and one New World monkey. As similar to those previously reported, the reactivity of anti-human IgG antibody with plasma from different primate species was closely related with phylogenic distance from human. Every anti-human IgG subclass antibody showed low cross-reactivity with plasma from Old World and New World monkeys. The plasma from all apes except for gibbons (Hylobates spp.) showed 60 to 100% of cross-reactivity with anti-human IgG2 and IgG3 antibodies. On the other hand, chimpanzee (Pan troglodytes andPan paniscus) and orangutan (Pongo pygmaeus) plasma showed 100% cross-reactivity with anti-human IgG1 antibody, but gorilla (Gorilla gorilla) and gibbon plasma showed no cross-reactivity. The chimpanzee and gorilla plasma cross-reacted with anti-human IgG4 antibody at different reactivity, 100% in chimpanzee and 50% in gorilla, but no cross-reactivity was observed in orangutan and gibbon plasma. These results suggest the possibilities that the divergence of “human-type” IgG subclasses might occur at the time of divergence ofHomo sapience fromHylobatidae, and that the molecular evolution of IgG1 as well as IgG4 is different from that of IgG2 and IgG3 in great apes, this is probably caused by different in development of immune function in apes during the course of evolution.     

Distinctive Patterns of Evolution of the δ-Globin Gene (HBD) in Primates

In most vertebrates, hemoglobin (Hb) is a heterotetramer composed of two dissimilar globin chains, which change during development according to the patterns of expression of α- and β-globin family members. In placental mammals, the β-globin cluster includes three early-expressed genes, ε(HBE)-γ(HBG)-ψβ(HBBP1), and the late expressed genes, δ (HBD) and β (HBB). While HBB encodes the major adult β-globin chain, HBD is weakly expressed or totally silent. Paradoxically, in human populations HBD shows high levels of conservation typical of genes under strong evolutionary constraints, possibly due to a regulatory role in the fetal-to-adult switch unique of Anthropoid primates. In this study, we have performed a comprehensive phylogenetic and comparative analysis of the two adult β-like globin genes in a set of diverse mammalian taxa, focusing on the evolution and functional divergence of HBD in primates. Our analysis revealed that anthropoids are an exception to a general pattern of concerted evolution in placental mammals, showing a high level of sequence conservation at HBD, less frequent and shorter gene conversion events. Moreover, this lineage is unique in the retention of a functional GATA-1 motif, known to be involved in the control of the developmental expression of the β-like globin genes. We further show that not only the mode but also the rate of evolution of the δ-globin gene in higher primates are strictly associated with the fetal/adult β-cluster developmental switch. To gain further insight into the possible functional constraints that have been shaping the evolutionary history of HBD in primates, we calculated dN/dS (ω) ratios under alternative models of gene evolution. Although our results indicate that HBD might have experienced different selective pressures throughout primate evolution, as shown by different ω values between apes and Old World Monkeys + New World Monkeys (0.06 versus 0.43, respectively), these estimates corroborated a constrained evolution for HBD in Anthropoid lineages, which is unlikely to be related to protein function. Collectively, these findings suggest that sequence change at the δ-globin gene has been under strong selective constraints over 65 Myr of primate evolution, likely due to a regulatory role in ontogenic switches of gene expression.     

Ubiquitin genes as a paradigm of concerted evolution of tandem repeats

Summary Ubiquitin is remarkable for its ubiquitous distribution and its extreme protein sequence conservation. Ubiquitin genes comprise direct repeats of the ubiquitin coding unit with no spacers. The nucleotide sequences of several ubiquitin repeats from each of humans, chicken,Xenopus, Drosophila, barley, and yeast have recently been determined. By analysis of these data we show that ubiquitin is evolving more slowly than any other known protein, and that this (together with its gene organization) contributes to an ideal situation for the occurrence of concerted evolution of tandem repeats. By contrast, there is little evidence of between-cluster concerted evolution. We deduce that in ubiquitin genes, concerted evolution involves both unequal crossover and gene conversion, and that the average time since two repeated units within the polyubiquitin locus most recently shared a common ancestor is approximately 38 million years (Myr) in mammals, but perhaps only 11 Myr inDrosophila. The extreme conservatism of ubiquitin evolution also allows the inference that certain synonymous serine codons differing at the first two positions were probably mutated at single steps.     

Natural selection in the <Emphasis Type="Italic">TLR-</Emphasis>related genes in the course of primate evolution

The innate immune system constitutes the front line of host defense against pathogens. Toll-like receptors (TLRs) recognize molecules derived from pathogens and play crucial roles in the innate immune system. Here, we provide evidence that the TLR-related genes have come under natural selection pressure in the course of primate evolution. We compared the nucleotide sequences of 16 TLR-related genes, including TLRs (TLR1–10), MYD88, TILAP, TICAM1, TICAM2, MD2, and CD14, among seven primate species. Analysis of the non-synonymous/synonymous substitution ratio revealed the presence of both strictly conserved and rapidly evolving regions in the TLR-related genes. The genomic segments encoding the intracellular Toll/interleukin 1 receptor domains, which exhibited lower rates of non-synonymous substitution, have undergone purifying selection. In contrast, TLR4, which carried a high proportion of non-synonymous substitutions in the part of extracellular domain spanning 200 amino acids, was found to have been the suggestive target of positive Darwinian selection in primate evolution. However, sequence analyses from 25 primate species, including eight hominoids, six Old World monkeys, eight New World monkeys, and three prosimians, showed no evidence that the pressure of positive Darwinian selection has shaped the pattern of sequence variations in TLR4 among New World monkeys and prosimians. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

MHC class I genes in a New World primate, the cotton-top tamarin (Saguinus oedipus), have evolved by an active process of loci turnover

 Lymphocytes of a New World primate, the cotton-top tamarin (Saguinus oedipus), express classical G–related major histocompatibility complex (MHC) class I molecules with unusually limited polymorphism and variability. Three G-related loci, an F locus, an E locus, and two pseudogenes (So-N1 and So-N3) have been identified by cDNA library screening and extensive PCR analysis of both cDNA and genomic DNA from the cotton-top tamarin. Furthermore, each genus of the subfamily Callitrichinae (tamarins and marmosets) appears to express its own unique set of MHC class I genes, likely due to a rapid turnover of loci. The rapid emergence of unique MHC class I genes in the Callitrichinae genera, resulting from an active process of duplication and inactivation of loci, may account for the limited diversity of the MHC class I genes in the cotton-top tamarin. To determine the nature of the entire complement of MHC class I genes in the cotton-top tamarin, we synthesized a genomic DNA library and screened it with MHC class I-specific probes. We isolated nine new MHC class I pseudogenes from this library. These newly isolated tamarin G–related MHC class I pseudogenes are not closely related to any of their functional counterparts in the tamarin, suggesting that they do not share a recent common ancestral gene with the tamarin's currently expressed MHC class I loci. In addition, these tamarin sequences display a high rate of nonsynonymous substitutions in their putative peptide binding region. This indicates that the genes from which they have derived were likely subject to positive selection and, therefore, were once functional. Our data support the notion that an extremely high rate of loci turnover is largely responsible for the limited diversity of the MHC class I genes in the cotton-top tamarin. Received: 15 September 1997 / Revised: 2 July 1998     

Unusually limited nucleotide sequence variation of the expressed major histocompatibility complex class I genes of a New World primate species (Saguinus oedipus)

Although major histocompatibility complex (MHC) class I molecules are, as a rule, highly polymorphic in mammalian species, those of the New World primate Saguinus oedipus (cotton-top tamarin) exhibit limited polymorphism. We have cloned and sequenced twelve MHC class I cDNAs from this species. Since cloned cotton-top tamarin cell lines express three to six MHC class I molecules, this species must have at least three functional MHC class I loci. There was, however, no evidence of locus-specific substitutions in the tamarin cDNAs. Unlike all other species studied, tamarin MHC class I cDNAs displayed limited nucleotide sequence variation. The sequence similarity between the two most divergent tamarin cDNAs was 95%. To ensure that the polymerase chain reaction (PCR) primers employed in these studies had amplified all of the tamarins' expressed MHC class I genes, we used another set of primers to amplify only exons 2 and 3 from RNA and DNA. PCR of genomic DNA resulted in the amplification of six distinct clones, of which only three were well expressed. Two of these nonexpressed genes were pseudogenes and the other was a nonclassical gene. Southern blot analysis demonstrated that the tamarin has 8–11 MHC class I genes, suggesting we had indeed cloned the majority of these genes. Cotton-top tamarins are, therefore, unique among mammalian species studied to date in that they express MHC class I molecules with limited nucleotide sequence variation.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers M38403-15.     

Polymorphic microsatellite loci for the mustached tamarin (Saguinus mystax) and their cross‐species amplification in other New World monkeys

We isolated and characterized 10 polymorphic microsatellite loci for the mustached and cotton‐top tamarin (Saguinus mystax and S. oedipus) and also tested their applicability in other New World monkey genera. Six loci proved to be variable in Cebus apella and Saimiri sciureus/S. boliviensis, respectively. We also show that faeces‐derived DNA is suitable for polymerase chain reaction (PCR) analyses of a primate taxon that gives birth to twins with individuals exhibiting a haematopoetic chimaerism. The microsatellite fingerprint patterns do not display a complex mixture of the PCR‐products derived from the alleles of both dizygotic twins, but exhibit the constitutive alleles of one sibling only.     

Differential Evolutionary Fate of an Ancestral Primate Endogenous Retrovirus Envelope Gene,the EnvV Syncytin,Captured for a Function in Placentation

Syncytins are envelope genes of retroviral origin that have been co-opted for a role in placentation. They promote cell–cell fusion and are involved in the formation of a syncytium layer—the syncytiotrophoblast—at the materno-fetal interface. They were captured independently in eutherian mammals, and knockout mice demonstrated that they are absolutely required for placenta formation and embryo survival. Here we provide evidence that these “necessary” genes acquired “by chance” have a definite lifetime with diverse fates depending on the animal lineage, being both gained and lost in the course of evolution. Analysis of a retroviral envelope gene, the envV gene, present in primate genomes and belonging to the endogenous retrovirus type V (ERV-V) provirus, shows that this captured gene, which entered the primate lineage >45 million years ago, behaves as a syncytin in Old World monkeys, but lost its canonical fusogenic activity in other primate lineages, including humans. In the Old World monkeys, we show—by in situ analyses and ex vivo assays—that envV is both specifically expressed at the level of the placental syncytiotrophoblast and fusogenic, and that it further displays signs of purifying selection based on analysis of non-synonymous to synonymous substitution rates. We further show that purifying selection still operates in the primate lineages where the gene is no longer fusogenic, indicating that degeneracy of this ancestral syncytin is a slow, lineage-dependent, and multi-step process, in which the fusogenic activity would be the first canonical property of this retroviral envelope gene to be lost.     

Molecular Evolution of Prolactin in Primates

Pituitary prolactin, like growth hormone (GH) and several other protein hormones, shows an episodic pattern of molecular evolution in which sustained bursts of rapid change contrast with long periods of slow evolution. A period of rapid change occurred in the evolution of prolactin in primates, leading to marked sequence differences between human prolactin and that of nonprimate mammals. We have defined this burst more precisely by sequencing the coding regions of prolactin genes for a prosimian, the slow loris (Nycticebus pygmaeus), and a New World monkey, the marmoset (Callithrix jacchus). Slow loris prolactin is very similar in sequence to pig prolactin, so the episode of rapid change occurred during primate evolution, after the separation of lines leading to prosimians and higher primates. Marmoset prolactin is similar in sequence to human prolactin, so the accelerated evolution occurred before divergence of New World monkeys and Old World monkeys/apes. The burst of change was confined largely to coding sequence (nonsynonymous sites) for mature prolactin and is not marked in other components of the gene sequence. This and the observations that (1) there was no apparent loss of function during the episode of rapid evolution, (2) the rate of evolution slowed toward the basal rate after this burst, and (3) the distribution of substitutions in the prolactin molecule is very uneven support the idea that this episode of rapid change was due to positive adaptive selection. In the slow loris and marmoset there is no evidence for duplication of the prolactin gene, and evidence from another New World monkey (Cebus albifrons) and from the chimpanzee and human genome sequences, suggests that this is the general position in primates, contrasting with the situation for GH genes. The chimpanzee prolactin sequence differs from that of human at two residues and comparison of human and chimpanzee prolactin gene sequences suggests that noncoding regions associated with regulating expression may be evolving differently from other noncoding regions.     

Different subcellular localisations of TRIM22 suggest species-specific function

The B30.2/SPRY domain is present in many proteins, including various members of the tripartite motif (TRIM) protein family such as TRIM5α, which mediates innate intracellular resistance to retroviruses in several primate species. This resistance is dependent on the integrity of the B30.2 domain that evolves rapidly in primates and exhibits species-specific anti-viral activity. TRIM22 is another positively selected TRIM gene. Particularly, the B30.2 domain shows rapid evolution in the primate lineage and recently published data indicate an anti-viral function of TRIM22. We show here that human and rhesus TRIM22 localise to different subcellular compartments and that this difference can be assigned to the positively selected B30.2 domain. Moreover, we could demonstrate that amino acid changes in two variable loops (VL1 and VL3) are responsible for the different subcellular localisations.     

Mhc-DRB Genes and the Origin of New World Monkeys

The major histocompatibility complex (Mhc) is a family of loci characterized by its relatively rapid evolutionary turnover, large genetic distances between genes, and long persistence of allelic lineages effected by balancing selection. These features render the Mhc highly suitable for answering questions concerning speciation and adaptive radiation. The aim of the present study was to use Mhc-DRB genes to make inferences about the founding population of the Platyrrhini. Three segments, each approximately 300 base pairs in length, of the platyrrhine DRB genes were amplified by the polymerase chain reaction and sequenced. The segments were derived from intron 2, exon 3, and exon 6 of DRB genes from different species of New World monkeys. The results of the study have revealed that on a phylogenetic tree, all of the tested platyrrhine genes appear to form a single cluster, while all catarrhine DRB genes form a distinct cluster, although the bootstrap values fail to provide statistically significant support for the separation of these two clades. This observation suggests that the multiple platyrrhine genes originated from a single ancestral gene after the divergence of the Platyrrhini and Catarrhini and thus contradicts the results of an earlier study in which some exon 2 DRB sequences appeared to predate the split of the two primate groups. The inconsistency in the DRB gene phylogeny can be explained by postulating convergent evolution for the peptide-binding region of the DRB exon 2 sequences. The phylogeny of the platyrrhine DRB genes (except for exon 2) is relatively "shallow"; the distances between genes are relatively short (in comparison to the catarrhine DRB genes), and there is a tendency for sequences of individual species to cluster together. The phylogeny of the platyrrhine DRB genes is consistent with the postulate that a small population founded the group and that there is an ongoing adaptive radiation from small, relatively isolated founding populations.     

True handedness in the cotton-top tamarin (Saguinus oedipus)?

Hand use for 8 activities was studied in 20 captive cotton-top tamarins (Saguinus oedipus). The purpose of the study was to discover if hand preferences existed across tasks and across subjects, and, if so, whether these followed the patterns of preference inMacNeilage et al.'s (1987) “postural origins” theory of the evolution of hand preference and hemispheric specialization in primates. This theory suggests that, for haplorines, the right hand is used preferentially for manipulative acts while the left hand is used preferentially for visually guided acts. The study showed statistically significant right hand preferences for six of the seven actions which produced sufficient data to be tested. The strength and consistency of this preference makes this study the first to suggest true handedness in a non-human primate species. Deviations from the pattern were seen in scratching and other actions when performed in a vertical posture, but these deviations took the form of a weakening of the right hand preference rather than a reversal. The findings of the study may be seen as supportingMacNeilage et al.'s (1987) theoryonly if the actions studied are all considered to be manipulative, which can be argued to be the case, though some of the actions were also visually guided.     

Evolutionary rate of immunoglobulin alpha noncoding region is greater in hominoids than in Old World monkeys.

Recent studies on the molecular evolution of primates show that the evolutionary rate among hominoids is considerably slower than that among nonhominoid primates. However, this observation at the nucleotide-sequence level is restricted to the beta-globin family region. In this study, we sequenced orthologous immunoglobulin alpha (C alpha) genes of chimpanzee, gorilla, orangutan, and crab-eating macaque (an Old World monkey) and compared them with that of the human by using noncoding regions for analysis. Since significant differences in rates among hominoids were not found by using the relative rate test, we evaluated the ratio (R) of the evolutionary distance between Old World monkey and human to the distance between orangutan and human. The R value (1.12) for the C alpha gene was much smaller than the expected value (1.38-2.33), showing that the nucleotide substitution rate (= mutation rate per year under selective neutrality) of the C alpha gene is greater in the human lineage than in the Old World monkey lineage. We also did a similar analysis for the gamma 1-, gamma 2-, psi eta-, and delta-globin genes and found a considerable heterogeneity (1.12-2.37) among the R values, including that for the C alpha gene. This indicates that the hominoid slowdown of the evolutionary rate is not a universal phenomenon in primate evolution.     

Phylogenetic Analysis of the Cytochrome P450 3 (CYP3) Gene Family

Cytochrome P450 genes (CYP) constitute a superfamily with members known from the Bacteria, Archaea, and Eukarya. The CYP3 gene family includes the CYP3A and CYP3B subfamilies. Members of the CYP3A subfamily represent the dominant CYP forms expressed in the digestive and respiratory tracts of vertebrates. The CYP3A enzymes metabolize a wide variety of chemically diverse lipophilic organic compounds. To understand vertebrate CYP3 diversity better, we determined the killifish (Fundulus heteroclitus) CYP3A30 and CYP3A56 and the ball python (Python regius) CYP3A42 sequences. We performed phylogenetic analyses of 45 vertebrate CYP3 amino acid sequences using a Bayesian approach. Our analyses indicate that teleost, diapsid, and mammalian CYP3A genes have undergone independent diversification and that the ancestral vertebrate genome contained a single CYP3A gene. Most CYP3A diversity is the product of recent gene duplication events. There is strong support for placement of the guinea pig CYP3A genes within the rodent CYP3A diversification. The rat, mouse, and hamster CYP3A genes are mixed among several rodent CYP3A subclades, indicative of a complex history involving speciation and gene duplication. Phylogenetic analyses suggest two CYP3A gene duplication events early in rodent history, with the rat CYP3A9 and mouse Cyp3a13 clade having a sister relationship to all other rodent CYP3A genes. In primate history, the human CYP3A43 gene appears to have a sister relationship to all other known primate CYP3A genes. Other, more recent gene duplications are hypothesized to have occurred independently within the human, pig, rat, mouse, guinea pig, and fish genomes. Functional analyses suggest that gene duplication is strongly tied to acquisition of new function and that convergent evolution of CYP3A function may be frequent among independent gene copies. Current address (Rachel L. Cox): Laboratory of Aquatic Biomedicine, Marine Biology Laboratory, Woods Hole, MA 02543, USA