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     

Cytochrome P450 1 genes in early deuterostomes (tunicates and sea urchins) and vertebrates (chicken and frog): origin and diversification of the CYP1 gene family

Cytochrome P450 family 1 (CYP1) proteins are important in a large number of toxicological processes. CYP1A and CYP1B genes are well known in mammals, but the evolutionary history of the CYP1 family as a whole is obscure; that history may provide insight into endogenous functions of CYP1 enzymes. Here, we identify CYP1-like genes in early deuterostomes (tunicates and echinoderms), and several new CYP1 genes in vertebrates (chicken, Gallus gallus and frog, Xenopus tropicalis). Profile hidden Markov models (HMMs) generated from vertebrate CYP1A and CYP1B protein sequences were used to identify 5 potential CYP1 homologs in the tunicate Ciona intestinalis genome. The C. intestinalis genes were cloned and sequenced, confirming the predicted sequences. Orthologs of 4 of these genes were found in the Ciona savignyi genome. Bayesian phylogenetic analyses group the tunicate genes in the CYP1 family, provisionally in 2 new subfamilies, CYP1E and CYP1F, which fall in the CYP1A and CYP1B/1C clades. Bayesian and maximum likelihood analyses predict functional divergence between the tunicate and vertebrate CYP1s, and regions within CYP substrate recognition sites were found to differ significantly in position-specific substitution rates between tunicates and vertebrates. Subsequently, 10 CYP1-like genes were found in the echinoderm Strongylocentrotus purpuratus (sea urchin) genome. Several of the tunicate and echinoderm CYP1-like genes are expressed during development. Canonical xenobiotic response elements are present in the upstream genomic sequences of most tunicate and sea urchin CYP1s, and both groups are predicted to possess an aryl hydrocarbon receptor (AHR), suggesting possible regulatory linkage of AHR and these CYPs. The CYP1 family has undergone multiple rounds of gene duplication followed by functional divergence, with at least one gene lost in mammals. This study provides new insight into the origin and evolution of CYP1 genes.     

Phylogenetic and functional analyses of the cytochrome P450 family 4

Cytochrome P450 family 4 (CYP4) proteins metabolize fatty acids, eicosanoids, and vitamin D and are important for chemical defense. The purpose of this study was to determine the evolutionary relationships between vertebrate CYP4 subfamilies and raise functional hypotheses regarding CYP4 subfamilies with little empirical data. 132 CYP4 sequences from 28 species were utilized for phylogenetic reconstructions by maximum likelihood and Bayesian inference. Monophyly was not found with the CYP4T and CYP4B subfamilies. CYP4V clustered with invertebrate subfamilies. Evolutionary rates of functional divergence were high in pairwise comparison with CYP4X yet, comparisons with mammalian CYP4F22 genes generally had no statistically significant divergence. Radical biochemical changes were detected in regions associated with substrate binding and the active site in comparisons among the CYP4A, CYP4X, and CYP4B subfamilies. Lastly, gene expression patterns, determined in silico with EST libraries from human, chicken, frog and fish, for CYP4V was markedly different between human and actinopterygian species. Further consideration should be given to the nomenclature of the CYP4T and CYP4B subfamily genes. Strong support was seen for the placement of CYP4A as a basal subfamily to CYP4X and CYP4Z. The B, B', J', K', K″ helices and a region at the end of C-terminus were suggested as conserved regions in CYP4 genes. The function of CYP4X was hypothesized to specialize in metabolism of long chain fatty acids. CYP4F22 genes may share a similar function to other CYP4F genes, although gene expression sites were different.     

Post-Messinian evolutionary relationships across the Sicilian channel: Mitochondrial and nuclear markers link a new green toad from Sicily to African relatives

Background

The major birch pollen allergen, Bet v 1, is a member of the ubiquitous PR-10 family of plant pathogenesis-related proteins. In recent years, a number of diverse plant proteins with low sequence similarity to Bet v 1 was identified. In addition, determination of the Bet v 1 structure revealed the existence of a large superfamily of structurally related proteins. In this study, we aimed to identify and classify all Bet v 1-related structures from the Protein Data Bank and all Bet v 1-related sequences from the Uniprot database.

Results

Structural comparisons of representative members of already known protein families structurally related to Bet v 1 with all entries of the Protein Data Bank yielded 47 structures with non-identical sequences. They were classified into eleven families, five of which were newly identified and not included in the Structural Classification of Proteins database release 1.71. The taxonomic distribution of these families extracted from the Pfam protein family database showed that members of the polyketide cyclase family and the activator of Hsp90 ATPase homologue 1 family were distributed among all three superkingdoms, while members of some bacterial families were confined to a small number of species. Comparison of ligand binding activities of Bet v 1-like superfamily members revealed that their functions were related to binding and metabolism of large, hydrophobic compounds such as lipids, hormones, and antibiotics. Phylogenetic relationships within the Bet v 1 family, defined as the group of proteins with significant sequence similarity to Bet v 1, were determined by aligning 264 Bet v 1-related sequences. A distance-based phylogenetic tree yielded a classification into 11 subfamilies, nine exclusively containing plant sequences and two subfamilies of bacterial proteins. Plant sequences included the pathogenesis-related proteins 10, the major latex proteins/ripening-related proteins subfamily, and polyketide cyclase-like sequences.

Conclusion

The ubiquitous distribution of Bet v 1-related proteins among all superkingdoms suggests that a Bet v 1-like protein was already present in the last universal common ancestor. During evolution, this protein diversified into numerous families with low sequence similarity but with a common fold that succeeded as a versatile scaffold for binding of bulky ligands.     

Cytochrome P450 diversification and hostplant utilization patterns in specialist and generalist moths: Birth,death and adaptation

Across insect genomes, the size of the cytochrome P450 monooxygenase (CYP) gene superfamily varies widely. CYPome size variation has been attributed to reciprocal adaptive radiations in insect detoxification genes in response to plant biosynthetic gene radiations driven by co‐evolution between herbivores and their chemically defended hostplants. Alternatively, variation in CYPome size may be due to random “birth‐and‐death” processes, whereby exponential increase via gene duplications is limited by random decay via gene death or transition via divergence. We examined CYPome diversification in the genomes of seven Lepidoptera species varying in host breadth from monophagous (Bombyx mori) to highly polyphagous (Amyelois transitella). CYPome size largely reflects the size of Clan 3, the clan associated with xenobiotic detoxification, and to some extent phylogenetic age. Consistently across genomes, families CYP6, CYP9 and CYP321 are most diverse and CYP6AB, CYP6AE, CYP6B, CYP9A and CYP9G are most diverse among subfamilies. Higher gene number in subfamilies is due to duplications occurring primarily after speciation and specialization (“P450 blooms”), and the genes are arranged in clusters, indicative of active duplicating loci. In the parsnip webworm, Depressaria pastinacella, gene expression levels in large subfamilies are high relative to smaller subfamilies. Functional and phylogenetic data suggest a correlation between highly dynamic loci (reflective of extensive gene duplication, functionalization and in some cases loss) and the ability of enzymes encoded by these genes to metabolize hostplant defences, consistent with an adaptive, nonrandom process driven by ecological interactions.     

Substrate recognition sites in cytochrome P450 family 2 (CYP2) proteins inferred from comparative analyses of amino acid and coding nucleotide sequences.

The substrate recognition regions in cytochrome P450 family 2 (CYP2) proteins were inferred by group-to-group alignment of CYP2 sequences and those of bacterial P450s, including Pseudomonas putida P450 101A (P450cam), whose substrate-binding residues have been definitely identified by x-ray crystallography of a substrate-bound form (Poulos T. L., Finzel, B. C., and Howard, A. J. (1987) J. Mol. Biol. 195, 687-700). The six putative substrate recognition sites, SRSs, thus identified are dispersively located along the primary structure and constitute about 16% of the total residues. All the reported point mutations and chimeric fragments that significantly affect the substrate specificities of the parental CYP2 enzymes fell within or overlapped some of the six SRSs. Analysis of nucleotide substitution patterns in closely related members in four subfamilies, CYP2A, 2B, 2C, and 2D, consistently indicated that the SRSs have accumulated more nonsynonymous (amino acid-changing) substitutions than the rest of the sequence. This observation supports the idea that diversification of duplicate genes of drug-metabolizing P450s occurs primarily in substrate recognition regions to cope with an increasing number of foreign compounds.     

Analysis of Species-Selectivity of Human,Mouse and Rat Cytochrome P450 1A and 2B Subfamily Enzymes using Molecular Modeling,Docking and Dynamics Simulations

Cytochrome P450 (CYP) 1A and 2B subfamily enzymes are important drug metabolizing enzymes, and are highly conserved across species in terms of sequence homology. However, there are major to minor structural and macromolecular differences which provide for species-selectivity and substrate-selectivity. Therefore, species-selectivity of CYP1A and CYP2B subfamily proteins across human, mouse and rat was analyzed using molecular modeling, docking and dynamics simulations when the chiral molecules quinine and quinidine were used as ligands. The three-dimensional structures of 17 proteins belonging to CYP1A and CYP2B subfamilies of mouse and rat were predicted by adopting homology modeling using the available structures of human CYP1A and CYP2B proteins as templates. Molecular docking and dynamics simulations of quinine and quinidine with CYP1A subfamily proteins revealed the existence of species-selectivity across the three species. On the other hand, in the case of CYP2B subfamily proteins, no role for chirality of quinine and quinidine in forming complexes with CYP2B subfamily proteins of the three species was indicated. Our findings reveal the roles of active site amino acid residues of CYP1A and CYP2B subfamily proteins and provide insights into species-selectivity of these enzymes across human, mouse, and rat.     

Evolution of the vertebrate ABC gene family: analysis of gene birth and death

Vertebrate evolution has been largely driven by the duplication of genes that allow for the acquisition of new functions. The ATP-binding cassette (ABC) proteins constitute a large and functionally diverse family of membrane transporters. The members of this multigene family are found in all cellular organisms, most often engaged in the translocation of a wide variety of substrates across lipid membranes. Because of the diverse function of these genes, their large size, and the large number of orthologs, ABC genes represent an excellent tool to study gene family evolution. We have identified ABC proteins from the sea squirt (Ciona intestinalis), zebrafish (Danio rerio), and chicken (Gallus gallus) and, using phylogenetic analysis, identified those genes with a one-to-one orthologous relationship to human ABC proteins. All ABC protein subfamilies found in Ciona and zebrafish correspond to the human subfamilies, with the exception of a single ABCH subfamily gene found only in zebrafish. Multiple gene duplication and deletion events were identified in different lineages, indicating an ongoing process of gene evolution. As many ABC genes are involved in human genetic diseases, and important drug transport phenotypes, the understanding of ABC gene evolution is important to the development of animal models and functional studies.     

cytochrome p450 superfamily in Arabidopsis thaliana: isolation of cDNAs,Differential Expression,and RFLP mapping of multiple cytochromes P450

We have isolated multiple cDNAs encoding cytochromes P450 (P450s) from Arabidopsis thaliana employing a PCR strategy. Degenerate oligonucleotide primers were designed from amino acid sequences conserved between two plant P450s, CYP71A1 and CYP73A2, including the heme-binding site and the proline-rich motif found in the N-terminal region, and 11 putative P450 fragments were amplified from first-strand cDNA from 7-day-old Arabidopsis as a template. With these PCR fragments as hybridization probes, 13 full-length and 3 partial cDNAs encoding different P450s have been isolated from an Arabidopsis cDNA library. These P450s have been assigned to either one of the established subfamilies: CYP71B, CYP73A, and CYP83A; or novel subfamilies: CYP76C, CYP83B, and CYP91A. The primary protein structures predicted from the cDNA sequences revealed that the regions around both the heme-binding site and the proline-rich motif were highly conserved among all these P450s. The N-terminal structures of the predicted P450 proteins suggested that these Arabidopsis P450s were located at the endoplasmic reticulum membrane. The loci of four P450 genes were determined by RFLP mapping. One of the clones, CYP71B2, was located at a position very close to the ga4 and gai mutations. RNA blot analysis showed expression patterns unique to each of the P450s in terms of tissue specificity and responsiveness to wounding and light/dark cycle, implicating involvement of these P450s in diverse metabolic processes.     

Evolution of EF-hand calcium-modulated proteins. III. Exon sequences confirm most dendrograms based on protein sequences: Calmodulin dendrograms show significant lack of parallelism

Summary In the first report in this series we presented dendrograms based on 152 individual proteins of the EF-hand family. In the second we used sequences from 228 proteins, containing 835 domains, and showed that eight of the 29 subfamilies are congruent and that the EF-hand domains of the remaining 21 subfamilies have diverse evolutionary histories. In this study we have computed dendrograms within and among the EF-hand subfamilies using the encoding DNA sequences. In most instances the dendrograms based on protein and on DNA sequences are very similar. Significant differences between protein and DNA trees for calmodulin remain unexplained. In our fourth report we evaluate the sequences and the distribution of introns within the EF-hand family and conclude that exon shuffling did not play a significant role in its evolution.     

Exceptional among-lineage variation in diversification rates during the radiation of Australia's most diverse vertebrate clade

The disparity in species richness among groups of organisms is one of the most pervasive features of life on earth. A number of studies have addressed this pattern across higher taxa (e.g. 'beetles'), but we know much less about the generality and causal basis of the variation in diversity within evolutionary radiations at lower taxonomic scales. Here, we address the causes of variation in species richness among major lineages of Australia's most diverse vertebrate radiation, a clade of at least 232 species of scincid lizards. We use new mitochondrial and nuclear intron DNA sequences to test the extent of diversification rate variation in this group. We present an improved likelihood-based method for estimating per-lineage diversification rates from combined phylogenetic and taxonomic (species richness) data, and use the method in a hypothesis-testing framework to localize diversification rate shifts on phylogenetic trees. We soundly reject homogeneity of diversification rates among members of this radiation, and find evidence for a dramatic rate increase in the common ancestor of the genera Ctenotus and Lerista. Our results suggest that the evolution of traits associated with climate tolerance may have had a role in shaping patterns of diversity in this group.     

Hepatic cytochrome P450 enzymes belonging to the CYP2C subfamily from an Australian marsupial, the koala (Phascolarctos cinereus)

Cytochromes P450 (CYPs) are critically important in the oxidative metabolism of a diverse array of xenobiotics and endogenous substrates. We have previously reported that the obligate Eucalyptus feeder koala (Phascolarctos cinereus) exhibits a higher hepatic CYP2C activity as compared to non-Eucalyptus feeders human or rat, with stimulation of CYP2C activity by cineole. In the present study, we examine CYP2C expression by immunohistochemistry and describe the identification and cloning of koala CYP2Cs. Utilising anti-rat CYP2C6 antibody, the expression of CYP2C was found to be uniform across the hepatic sections, being consistent with that observed in human and rat. Two 1647 and 1638 bp koala liver CYP2C complete cDNAs, designated CYP2C47 and CYP2C48 respectively, were cloned by cDNA library screening. The koala CYP2C cDNAs encode a protein of 495 amino acids. Three additional partial CYP2C sequences were also identified from the koala, indicating the multiplicity of the CYP2C subfamily in this unique marsupial species. The results of this study demonstrate the presence of koala hepatic CYP2Cs that share several common features with other published CYP2Cs; however CYP2C47 and CYP2C48 contain four extra amino acid residues at the NH2-terminal, a transmembrane anchor which was reported being a fundamentally conserved structure core of all eukaryote CYP enzymes.     

Time scale for cyclostome evolution inferred with a phylogenetic diagnosis of hagfish and lamprey cDNA sequences

The Cyclostomata consists of the two orders Myxiniformes (hagfishes) and Petromyzoniformes (lampreys), and its monophyly has been unequivocally supported by recent molecular phylogenetic studies. Under this updated vertebrate phylogeny, we performed in silico evolutionary analyses using currently available cDNA sequences of cyclostomes. We first calculated the GC-content at four-fold degenerate sites (GC(4)), which revealed that an extremely high GC-content is shared by all the lamprey species we surveyed, whereas no striking pattern in GC-content was observed in any of the hagfish species surveyed. We then estimated the timing of diversification in cyclostome evolution using nucleotide and amino acid sequences. We obtained divergence times of 470-390 million years ago (Mya) in the Ordovician-Silurian-Devonian Periods for the interordinal split between Myxiniformes and Petromyzoniformes; 90-60 Mya in the Cretaceous-Tertiary Periods for the split between the two hagfish subfamilies, Myxininae and Eptatretinae; 280-220 Mya in the Permian-Triassic Periods for the split between the two lamprey subfamilies, Geotriinae and Petromyzoninae; and 30-10 Mya in the Tertiary Period for the split between the two lamprey genera, Petromyzon and Lethenteron. This evolutionary configuration indicates that Myxiniformes and Petromyzoniformes diverged shortly after the common ancestor of cyclostomes split from the future gnathostome lineage. Our results also suggest that intra-subfamilial diversification in hagfish and lamprey lineages (especially those distributed in the northern hemisphere) occurred in the Cretaceous or Tertiary Periods.     

Toxin-resistant sodium channels: parallel adaptive evolution across a complete gene family

Approximately 75% of vertebrate proteins belong to protein families encoded by multiple evolutionarily related genes, a pattern that emerged as a result of gene and genome duplications over the course of vertebrate evolution. In families of genes with similar or related functions, adaptation to a strong selective agent should involve multiple adaptive changes across the entire gene family. However, we know of no evolutionary studies that have explicitly addressed this point. Here, we show how 4 taxonomically diverse species of pufferfishes (Tetraodontidae) each evolved resistance to the guanidinium toxins tetrodotoxin (TTX) and saxitoxin (STX) via parallel amino acid replacements across all 8 sodium channels present in teleost fish genomes. This resulted in diverse suites of coexisting sodium channel types that all confer varying degrees of toxin resistance, yet show remarkable convergence among genes and phylogenetically diverse species. Using site-directed mutagenesis and expression of a vertebrate sodium channel, we also demonstrate that resistance to TTX/STX is enhanced up to 15-fold by single, frequently observed replacements at 2 sites that have not previously been implicated in toxin binding but show similar or identical replacements in pufferfishes and in distantly related vertebrate and nonvertebrate animals. This study presents an example of natural selection acting upon a complete gene family, repeatedly arriving at a diverse but limited number of adaptive changes within the same genome. To be maximally informative, we suggest that future studies of molecular adaptation should consider all functionally similar paralogs of the affected gene family.     

Amino Acid Patterns Within Short Consensus Repeats Define Conserved Duplicons Shared by Genes of the RCA Complex

Complement control proteins (CCPs) contain repeated protein domains, short consensus repeats (SCRs), which must be relevant to diverse functions such as complement activation, coagulation, viral binding, fetal implantation, and self-nonself recognition. Although SCRs share some discontinuous and imperfect motifs, there are many variable positions and indels making classification in subfamilies extremely difficult. Using domain-by-domain phylogenetic analysis, we have found that most domains can be classified into only 11 subfamilies, designated a, b, c, d, e, f, g, h, i, j, or k and identified by critical residues. Each particular CCP is characterized by the order of representatives of the subfamilies. Human complement receptor 1 (CR1) has ajefbkd repeated four times and followed by ch. The classification crosses CCPs and indicates that a particular CCP is a function of the mix of SCRs. The aje set is a feature of several CCPs including human CR1 and DAF and murine Crry and appears to be associated with the success or failure of implantation inter alia. This approach facilitates genomic analysis of available sequences and suggests a framework for the evolution of CCPs. Units of duplication range from single SCRs, to septamers such as efbkdaj, to extensive segments such as MCP-CR1L. Imperfections of duplication with subsequent deletion have contributed to diversification.     

Isolation of Resistance Gene Candidates (RGCs) and characterization of an RGC cluster in cassava

Plant disease resistance genes (R genes) show significant similarity amongst themselves in terms of both their DNA sequences and structural motifs present in their protein products. Oligonucleotide primers designed from NBS (Nucleotide Binding Site) domains encoded by several R-genes have been used to amplify NBS sequences from the genomic DNA of various plant species, which have been called Resistance Gene Analogues (RGAs) or Resistance Gene Candidates (RGCs). Using specific primers from the NBS and TIR (Toll/Interleukin-1 Receptor) regions, we identified twelve classes of RGCs in cassava ( Manihot esculenta Crantz). Two classes were obtained from the PCR-amplification of the TIR domain. The other 10 classes correspond to the NBS sequences and were grouped into two subfamilies. Classes RCa1 to RCa5 are part of the first subfamily and were linked to a TIR domain in the N terminus. Classes RCa6 to RCa10 corresponded to non-TIR NBS-LRR encoding sequences. BAC library screening with the 12 RGC classes as probes allowed the identification of 42 BAC clones that were assembled into 10 contigs and 19 singletons. Members of the two TIR and non-TIR NBS-LRR subfamilies occurred together within individual BAC clones. The BAC screening and Southern hybridization analyses showed that all RGCs were single copy sequences except RCa6 that represented a large and diverse gene family. One BAC contained five NBS sequences and sequence analysis allowed the identification of two complete RGCs encoding two highly similar proteins. This BAC was located on linkage group J with three other RGC-containing BACs. At least one of these genes, RGC2, is expressed constitutively in cassava tissues.Communicated by M.-A. Grandbastien     

Ancestral Sequence Reconstruction of a Cytochrome P450 Family Involved in Chemical Defense Reveals the Functional Evolution of a Promiscuous,Xenobiotic-Metabolizing Enzyme in Vertebrates

The cytochrome P450 family 1 enzymes (CYP1s) are a diverse family of hemoprotein monooxygenases, which metabolize many xenobiotics including numerous environmental carcinogens. However, their historical function and evolution remain largely unstudied. Here we investigate CYP1 evolution via the reconstruction and characterization of the vertebrate CYP1 ancestors. Younger ancestors and extant forms generally demonstrated higher activity toward typical CYP1 xenobiotic and steroid substrates than older ancestors, suggesting significant diversification away from the original CYP1 function. Caffeine metabolism appears to be a recently evolved trait of the CYP1A subfamily, observed in the mammalian CYP1A lineage, and may parallel the recent evolution of caffeine synthesis in multiple separate plant species. Likewise, the aryl hydrocarbon receptor agonist, 6-formylindolo[3,2-b]carbazole (FICZ) was metabolized to a greater extent by certain younger ancestors and extant forms, suggesting that activity toward FICZ increased in specific CYP1 evolutionary branches, a process that may have occurred in parallel to the exploitation of land where UV-exposure was higher than in aquatic environments. As observed with previous reconstructions of P450 enzymes, thermostability correlated with evolutionary age; the oldest ancestor was up to 35 °C more thermostable than the extant forms, with a ¹⁰T₅₀ (temperature at which 50% of the hemoprotein remains intact after 10 min) of 71 °C. This robustness may have facilitated evolutionary diversification of the CYP1s by buffering the destabilizing effects of mutations that conferred novel functions, a phenomenon which may also be useful in exploiting the catalytic versatility of these ancestral enzymes for commercial application as biocatalysts.