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.     

The cytochrome P450 gene superfamily in Drosophila melanogaster: annotation, intron-exon organization and phylogeny

The cytochrome P450 gene superfamily is represented by 90 sequences in the Drosophila melanogaster genome. Of these 90 P450 sequences, 83 code for apparently functional genes whereas seven are apparent pseudogenes. More than half of the genes belong to only two families, CYP4 and CYP6. The CYP6 family is insect specific whereas the CYP4 family includes sequences from vertebrates. There are eight genes coding for mitochondrial P450s as deduced from their homology to CYP12A1 from the house fly. The genetic map of the distribution of D. melanogaster P450 genes shows (a) the absence of P450 genes on the chromosome 4 and Y, (b) more than half of the P450 genes are found on chromosome 2, and (c) the largest cluster contains nine genes. Sequence alignments were used to draw phylogenetic trees and to analyze the intron-exon organization of each functional P450 gene. Only five P450 genes are intronless. We found 57 unique intron positions, of which 23 were phase zero, 19 were phase one and 15 were phase two. There was a relatively good correlation between intron conservation and phylogenetic relationship between members of the P450 subfamilies. Although the function of many P450 proteins from vertebrates, fungi, plants and bacteria is known, only a single P450 from D. melanogaster, CYP6A2, has been functionally characterized. Gene organization appears to be a useful tool in the study of the regulation, the physiological role and the function of these P450s.     

A Revised Evolutionary History of the CYP1A Subfamily: Gene Duplication, Gene Conversion, and Positive Selection

Members of cytochrome P450 subfamily 1A (CYP1As) are involved in detoxification and bioactivation of common environmental pollutants. Understanding the functional evolution of these genes is essential to predicting and interpreting species differences in sensitivity to toxicity caused by such chemicals. The CYP1A gene subfamily comprises a single ancestral representative in most fish species and two paralogs in higher vertebrates, including birds and mammals. Phylogenetic analysis of complete coding sequences suggests that mammalian and bird paralog pairs (CYP1A1/2 and CYP1A4/5, respectively) are the result of independent gene duplication events. However, comparison of vertebrate genome sequences revealed that CYP1A genes lie within an extended region of conserved fine-scale synteny, suggesting that avian and mammalian CYP1A paralogs share a common genomic history. Algorithms designed to detect recombination between nucleotide sequences indicate that gene conversion has homogenized most of the length of the chicken CYP1A genes, as well as the 5′ end of mammalian CYP1As. Together, these data indicate that avian and mammalian CYP1A paralog pairs resulted from a single gene duplication event and that extensive gene conversion is responsible for the exceptionally high degree of sequence similarity between CYP1A4 and CYP1A5. Elevated nonsynonymous/synonymous substitution ratios within a putatively unconverted stretch of ∼250 bp suggests that positive selection may have reduced the effective rate of gene conversion in this region, which contains two substrate recognition sites. This work significantly alters our understanding of functional evolution in the CYP1A subfamily, suggesting that gene conversion and positive selection have been the dominant processes of sequence evolution. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Yves Van de Peer]     

A cytochrome P4501B gene from a fish, Pleuronectes platessa

Tetrapod cytochrome P4501 family (CYP1A1, CYP1A2 and CYP1B1) enzymes are most active in hydroxylating a variety of environmental contaminants including polyaromatic hydrocarbons (PAH), planar polychlorinated biphenyls and arylamines and thus play a pivotal role in the toxicology of these compounds. Mammalian CYP1A1 and CYP1A2 genes appear to have diverged after the evolutionary emergence of mammals, whereas fish species apparently possess only one CYP1A family gene, and fish CYP1A enzymes exhibit properties of both of the mammalian isoforms. We have isolated a further CYP1 family gene from a marine flatfish (plaice; Pleuronectes platessa), which, on the basis of exon organisation and sequence similarity, can be assigned as a piscine CYP1B. Its deduced amino acid sequence shows the closest (54%) identity to mammalian CYP1B1 proteins and, on the basis of molecular modeling studies, shows a high degree of positional and structural conservation of the substrate contacting amino acid residues in its putative active site when compared to other CYP1 enzymes. Phylogenetic analysis of fish and mammalian CYP1 family sequences indicates that the plaice CYP1B and mammalian CYP1B1 genes share a common ancestry. Plaice CYP1B has a more restricted tissue expression profile than the previously isolated plaice CYP1A, only being detectable, by Northern blotting, in gill tissue. In contrast to CYP1A, which shows extensive PAH-dependent induction in a variety of tissues, plaice CYP1B appears unresponsive to treatment with a prototypical PAH-type inducer, beta-naphthoflavone (BNF).     

哺乳动物抗凝血酶基因的分子特征、微共线性与适应性进化

抗凝血酶(AT)是哺乳动物体内重要的天然抗凝血因子之一,它隶属于丝氨酸蛋白酶抑制物家族,主要参与并调节复杂的血凝过程。基于比较基因组学手段,共挖掘出17个哺乳动物抗凝血酶基因(AT),并剖析了它们的基因结构、微共线性、保守基序、功能结构域、以及系统进化关系。基因结构与微共线分析表明,哺乳动物AT基因具有5-12个外显子,大多数是7个外显子;不同物种AT基因所处区段之间具有较好的共线性。哺乳动物AT蛋白特征分析显示,SERPIN功能结构域与保守基序1、2、3、4、5、8和9存在相互重叠的部分。AT基因进化树揭示基因进化和通常认为的物种进化几乎一致。同时,利用PAML中Codeml的位点特异模型在AT基因中发掘了1个正选择位点328D。     

High-mobility-group family genes from Lampetra japonica reveal their early origin and molecular evolution in the vertebrate lineage

High-mobility group family (HMG) genes are ubiquitous in vertebrates, including mammals, birds, amphibians and fishes. To elucidate the molecular phylogeny of the HMG genes in the primitive vertebrate, we have cloned three homologues of HMG-box genes, called Lj-HMGB1, Lj-HMGB2 and Lj-HMGBX, from a cDNA library generated from lymphocyte-like cells of the Japanese lamprey (Lampetra japonica), an Agnathan that occupies a critical phylogenetic position between invertebrates and vertebrates. The open reading frames of Lj-HMGB1, Lj-HMGB2 and Lj-HMGBX contained 627 bp, 585 bp and 678 bp, respectively. The analysis of the deduced amino acid sequences indicated that these three putative Lj-HMGB proteins contain four domains: HMG-box A, HMG-box B, an acidic carboxyl-terminal tail and a linker. A phylogenetic analysis revealed that the Lj-HMGB proteins fall outside the vertebrate clade; Lj-HMGBX is descended from a gene ancestral to the mammalian HMGB1/2/3. This discovery implies that there was a gene duplication event in the HMGB1/2/3 gene family that occurred after the divergence of the vertebrates (Cyclostomata) from the Cephalochordata and Urochordata at least 450 million years ago (MYA). The Lj-HMGB1, Lj-HMGB2 and Lj-HMGBX genes were detected in most tissues of the lamprey by RT-PCR. Our findings provide insight into the phylogeny of the HMGB genes in vertebrates.     

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.     

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     

Phylogenetic and Functional Analysis of the Vertebrate Cytochrome P450 2 Family

Cytochrome P450 (CYP) proteins compose a highly diverse superfamily found in all domains of life. These proteins are enzymes involved in metabolism of endogenous and exogenous compounds. In vertebrates, the CYP2 family is one of the largest, most diverse and plays an important role in mammalian drug metabolism. However, there are more than 20 vertebrate CYP2 subfamilies with uncertain evolution and fairly discrete subfamily composition within vertebrate classes, hindering extrapolation of knowledge across subfamilies. To better understand CYP2 diversity, a phylogenetic analysis of 196 CYP2 protein sequences from 16 species was performed using a maximum likelihood approach and Bayesian inference. The analyses included the CYP2 compliment from human, fugu, zebrafish, stickleback, medaka, cow, and dog genomes. Additional sequences were included from rabbit, marsupial, platypus, chicken, frog, and salmonid species. Three CYP2 sequences from the tunicate Ciona intestinalis were utilized as the outgroup. Results indicate a single ancestral vertebrate CYP2 gene and monophyly of all CYP2 subfamilies. Two subfamilies (CYP2R and CYP2U) pre-date vertebrate diversification, allowing direct comparison across vertebrate classes, while all other subfamilies originated during vertebrate diversification, often within specific vertebrate lineages. Analysis of site-specific evolution indicates that some substrate recognition sites (SRS) previously proposed for CYP genes do not have elevated rates of evolution, suggesting that these regions of the protein are not necessarily important in recognition of CYP2 substrates. Type II functional divergence analysis identified multiple residues in the active site of CYP2F, CYP2A, and CYP2B proteins that have undergone radical biochemical changes and may be functionally important.     

The mammalian alphaD-globin gene lineage and a new model for the molecular evolution of alpha-globin gene clusters at the stem of the mammalian radiation

We have explored the evolution of the alpha-globin gene family by comparative sequence and phylogenetic analyses of mammalian alpha-globin genes. Our analyses reveal the existence of a new alpha-globin gene lineage in mammals that is related to the alpha(D)-globin genes of birds, squamates and turtles. The gene is located in the middle of the alpha-globin gene cluster of a marsupial, Sminthopsis macroura and of humans. It exists in a wide variety of additional mammals, including pigs, cows, cats, and dogs, but is a pseudogene in American marsupials. Evolutionary analyses suggest that the gene has generally evolved under purifying selection, indicative of a functional gene. The presence of mRNA products in humans, pigs, and cows also suggest that the gene is expressed and likely to be functional. The analyses support the hypothesis that the alpha(D)-globin gene lineage has an ancient evolutionary origin that predates the divergence of amniotes. The structural similarity of alpha-globin gene clusters of marsupials and humans suggest that an eight gene cluster (5'-zeta2-zeta1-alpha(D)-alpha3-alpha2-alpha1-theta-omega-3'), including seven alpha-like genes and one beta-like globin gene (omega-globin) existed in the common ancestor of all marsupial and eutherian mammals. This basic structure has remained relatively stable in marsupials and in the lineage leading to humans, although omega-globin has been lost from the alpha-globin gene cluster of humans.     

Molecular evolution of the polyamine oxidase gene family in Metazoa

ABSTRACT: BACKGROUND: Polyamine oxidase enzymes catalyze the oxidation of polyamines and acetylpolyamines. Since polyamines are basic regulators of cell growth and proliferation, their homeostasis is crucial for cell life. Members of the polyamine oxidase gene family have been identified in a wide variety of animals, including vertebrates, arthropodes, nematodes, placozoa, as well as in plants and fungi. Polyamine oxidases (PAOs) from yeast can oxidize spermine, N1-acetylspermine, and N1-acetylspermidine, however, in vertebrates two different enzymes, namely spermine oxidase (SMO) and acetylpolyamine oxidase (APAO), specifically catalyze the oxidation of spermine, and N1-acetylspermine/N1-acetylspermidine, respectively. Little is known about the molecular evolutionary history of these enzymes. However, since the yeast PAO is able to catalyze the oxidation of both acetylated and non acetylated polyamines, and in vertebrates these functions are addressed by two specialized polyamine oxidase subfamilies (APAO and SMO), it can be hypothesized an ancestral reference for the former enzyme from which the latter would have been derived. RESULTS: We analysed 36 SMO, 26 APAO and 14 PAO homologue protein sequences from 54 taxa including various vertebrates and invertebrates. The analysis of the full-length sequences and the principal domains of vertebrate and invertebrate PAOs yielded consensus primary protein sequences for vertebrate SMOs and APAOs, and invertebrate PAOs. This analysis, coupled to molecular modeling techniques, also unveiled sequence regions that confer specific structural and functional properties, including substrate specificity, by the different PAO subfamilies. Molecular phylogenetic trees revealed a basal position of all the invertebrates PAO enzymes relative to vertebrate SMOs and APAOs. PAOs from insects constitute a monophyletic clade. Two PAO variants sampled in the amphioxus are basal to the dichotomy between two well supported monophyletic clades including, respectively, all the SMOs and APAOs from vertebrates. The two vertebrate monophyletic clades clustered strictly mirroring the organismal phylogeny of fishes, amphibians, reptiles, birds, and mammals. Evidences from comparative genomic analysis, structural evolution and functional divergence in a phylogenetic framework across Metazoa suggested an evolutionary scenario where the ancestor PAO coding sequence, present in invertebrates as an orthologous gene, has been duplicated in the vertebrate branch to originate the paralogous SMO and APAO genes. A further genome evolution event concerns the SMO gene of placental, but not marsupial and monotremate, mammals which increased its functional variation following an alternative splicing (AS) mechanism. CONCLUSIONS: In this study the explicit integration in a phylogenomic framework of phylogenetic tree construction, structure prediction, and biochemical function data/prediction, allowed inferring the molecular evolutionary history of the PAO gene family and to disambiguate paralogous genes related by duplication event (SMO and APAO) and orthologous genes related by speciation events (PAOs, SMOs/APAOs). Further, while in vertebrates experimental data corroborate SMO and APAO molecular function predictions, in invertebrates the finding of a supported phylogenetic clusters of insect PAOs and the co-occurrence of two PAO variants in the amphioxus urgently claim the need for future structure-function studies.     

Relaxin gene family in teleosts: phylogeny, syntenic mapping, selective constraint, andexpression analysis

Background  

In recent years, the relaxin family of signaling molecules has been shown to play diverse roles in mammalian physiology, but little is known about its diversity or physiology in teleosts, an infraclass of the bony fishes comprising ~ 50% of all extant vertebrates. In this paper, 32 relaxin family sequences were obtained by searching genomic and cDNA databases from eight teleost species; phylogenetic, molecular evolutionary, and syntenic data analyses were conducted to understand the relationship and differential patterns of evolution of relaxin family genes in teleosts compared with mammals. Additionally, real-time quantitative PCR was used to confirm and assess the tissues of expression of five relaxin family genes in Danio rerio and in situ hybridization used to assess the site-specific expression of the insulin 3-like gene in D. rerio testis.     

Phylogenetic analysis of vertebrate kininogen genes

Kininogens, the precursors of bradykinins, vary extremely in both structure and function among different taxa of animals, in particular between mammals and amphibians. This includes even the most conserved bradykinin domain in terms of biosynthesis mode and structure. To elucidate the evolutionary dynamics of kininogen genes, we have identified 19 novel amino acid sequences from EST and genomic databases (for mammals, birds, and fishes) and explored their phylogenetic relationships using combined amino acid sequence and gene structure as markers. Our results show that there were initially two paralogous kininogen genes in vertebrates. During their evolution, the original gene was saved with frequent multiplication in amphibians, but lost in fishes, birds, and mammals, while the novel gene was saved with multiple functions in fishes, birds, and mammals, but became a pseudogene in amphibians. We also propose that the defense mechanism against specific predators in amphibian skin secretions has been bradykinin receptor dependent. Our findings may provide a foundation for identification and structural, functional, and evolutionary analyses of more kininogen genes and other gene families.     

Evolutionary origins of the transforming growth factor-beta gene family.

A molecular phylogeny for the transforming growth factor-beta (TGF-beta) gene family based on a comparison of nucleotide sequences is proposed. A phylogenetic tree constructed from these sequences shows that the family evolved from a common ancestral gene that came into existence at about the time of arthropod and chordate divergence. This model suggests that the present day TGF-beta gene family consists of four members: TGF-beta 1 (= TGF-beta 4), TGF-beta 2, TGF-beta 3, and TGF-beta 5. The molecular phylogeny and Southern hybridization data also suggest that the proteins for mammalian TGF-beta 1 and chicken TGF-beta 4 are the products of homologous rather than duplicated genes. If the gene duplication event that produced the ancestral gene for TGF-beta 1 occurred before the divergence of birds and mammals, then sufficient time would have elapsed to generate these quite distinct avian and mammalian TGF-beta 1 proteins. Therefore, the TGF-beta family contains four distinct proteins, TGF-beta 1, 2, 3, and 5.     

Identification and characterisation of a cytochrome P450 gene and processed pseudogene from an arachnid: the cattle tick, Boophilus microplus.

We isolated and sequenced the first known cytochrome P450 gene and pseudogene from an arachnid, the cattle tick, Boophilus microplus. Both the gene and pseudogene belong to the family CYP4, but a new subfamily, CYP4W, had to be created for these genes because they are substantially different to other CYP4 genes. The gene, CYP4W1, has greatest homology with CYP4C1 from a cockroach, Blaberus discoidalis. The predicted molecular weight of the protein encoded by CYP4W1 (63 KDa) is greater than that of the other CYP4 genes. The pseudogene, CYP4W1P, is probably a processed pseudogene derived from the functional gene CYP4W1. This is only the third CYP processed pseudogene to be identified. The pseudogene is 98% identical to the functional gene, CYP4W1, therefore we hypothesise that this pseudogene evolved recently from the functional gene. The CYP4 genes from arthropods have diverged from each other more than those of mammals; consequently the phylogeny of the arthropod genes could not be resolved.     

Identification of the mammalian Not gene via a phylogenomic approach

Despite the great morphological diversity of early embryos, the underlying mechanisms of gastrulation are known to be broadly conserved in vertebrates. However, a number of genes characterized as fulfilling an essential function in this process in several model organisms display no clear ortholog in mammalian genomes. We have devised an in silico phylogenomic approach, based on exhaustive similarity searches in vertebrate genomes and subsequent bayesian phylogenetic analyses, to identify such missing genes, presumed to be highly divergent. This approach has been used to identify mammalian orthologs of Not, an homeodomain containing gene previously characterized in Xenopus, chick and zebrafish as playing a critical role in the formation of the notochord. This attempt led to the identification of a highly divergent mammalian Not-related gene in the mouse, human and rat. The results from phylogenetic reconstructions, synteny analyses, expression pattern analyses in wild-type and mutant mouse embryos, and overexpression experiments in Xenopus embryos converge to confirm these genes as representatives of the Not family in mammals. The identification of the mammalian Not gene delivers an important component for the understanding of the genetics underlying notochord formation in mammals and its evolution among vertebrates. The phylogenomic method used to retrieve this gene thus provides a tool, which can complement or validate genome annotations in situations when they are weakly supported.     

Origin and evolution of the regulatory gene male-specific lethal-3

Dosage compensation in Drosophila is mediated by genes known as "male-specific lethals" (msls). Several msls, including male-specific lethal-3 (msl-3), encode proteins of unknown function. We cloned the Drosophila virilis msl-3 gene. Using the information provided by the sequences of the Drosophila melanogaster and D. virilis genes, we found that sequences of other species can be aligned along their entire lengths with msl-3. Among them, there are genes in yeasts (the Schizosaccharomyces pombe Alp13 gene, as well as a putative Alp13 homolog, found in Saccharomyces cerevisae) and in mammals (MRG15 and MSL3L1 and their relatives) plus uncharacterized sequences of the nematode Caenorhabditis elegans and the plants Arabidopsis thaliana, Lycopersicon esculentum, and Zea mays. A second Drosophila gene of this family has also been found. It is thus likely that msl-3-like genes are present in all eukaryotes. Phylogenetic analyses suggest that msl-3 is orthologous to the mammalian MSL3L1 genes, while the second Drosophila melanogaster gene (which we have called Dm MRG15) is orthologous to mammalian MRG15. These analyses also suggest that the msl-3/MRG15 duplication occurred after the fungus/animal split, while an independent duplication occurred in plants. The proteins encoded by these genes have similar structures, including a putative chromodomain close to their N-terminal end and a putative leucine zipper at their C-terminus. The possible functional roles of these proteins are discussed.