CYP3 phylogenomics: evidence for positive selection of CYP3A4 and CYP3A7

Signatures of Co-evolution and Co-regulation in the CYP3A and CYP4F Genes in Humans

Cytochromes P450 (CYP450) are hemoproteins generally involved in the detoxification of the body of xenobiotic molecules. They participate in the metabolism of many drugs and genetic polymorphisms in humans have been found to impact drug responses and metabolic functions. In this study, we investigate the genetic diversity of CYP450 genes. We found that two clusters, CYP3A and CYP4F, are notably differentiated across human populations with evidence for selective pressures acting on both clusters: we found signals of recent positive selection in CYP3A and CYP4F genes and signals of balancing selection in CYP4F genes. Furthermore, an extensive amount of unusual linkage disequilibrium is detected in this latter cluster, indicating co-evolution signatures among CYP4F genes. Several of the selective signals uncovered co-localize with expression quantitative trait loci (eQTL), which could suggest epistasis acting on co-regulation in these gene families. In particular, we detected a potential co-regulation event between CYP3A5 and CYP3A43, a gene whose function remains poorly characterized. We further identified a causal relationship between CYP3A5 expression and reticulocyte count through Mendelian randomization analyses, potentially involving a regulatory region displaying a selective signal specific to African populations. Our findings linking natural selection and gene expression in CYP3A and CYP4F subfamilies are of importance in understanding population differences in metabolism of nutrients and drugs.

The protective role of hydrophilic tetrahydroxylated bile acids (THBA)

Bile acids are key components of bile required for human health. In humans and mice, conditions of reduced bile flow, cholestasis, induce bile acid detoxification by producing tetrahydroxylated bile acids (THBA), more hydrophilic and less cytotoxic than the usual bile acids, which are typically di- or tri-hydroxylated. Mice deficient in the Bile Salt Export Pump (Bsep, or Abcb11), the primary bile acid transporter in liver cells, produce high levels of THBA, and avoid the severe liver damage typically seen in humans with BSEP deficiencies. THBA can suppress bile acid-induced liver damage in Mdr2-deficient mice, caused by their lack of phospholipids in bile exposing their biliary tracts to unbound bile acids. Here we review THBA-related works in both animals and humans, and discuss their potential relevance and applications as a class of functional bile acids.

Diversifying selection detected in only a minority of xenobiotic-metabolizing CYP1-3 genes among primate species

1. Primates exhibit a high degree of among-species dietary diversity, which likely exposes them to varying levels of xenobiotic compounds. Here, we examined the evolution of primate CYP1-3 gene families, and we classified the 15 CYP1-3 gene subfamilies as either xenobiotic-metabolizing (XM) or endogenous-metabolizing (EM) based on sources in the P450 literature. 2. We predicted that XM P450s would show (1) greater variability in gene-copy number and (2) more evidence of diversifying selection and, especially on codons that encode the substrate-recognition sites (SRSs) for the final enzymes. 3. Counter to our first prediction, EM and XM P450s showed similar levels of variation in gene-copy number. We did find, however, that four XM P450 subfamilies (CYP2C, CYP2D, CYP2E, and CYP3A) showed evidence of diversifying selection while no EM subfamilies demonstrated any consistent signal of diversifying selection. Of these four, CYP2C, CYP2D, and CYP3A showed significant links between SRSs and diversifying selection. 4. These results reveal an amount of evolutionary dynamism that would not be expected when viewing P450 subfamilies along a simple binary EM/XM spectrum. We recommend that comparative studies of cytochrome P450 evolution should focus on the CYP2C, CYP2D, CYP2E, and CYP3A subfamilies.

The marmoset cytochrome P450 superfamily: Sequence/phylogenetic analyses, genomic structure, and catalytic function

The common marmoset (Callithrix jacchus) is a New World monkey that has attracted much attention as a potentially useful primate model for preclinical testing. A total of 36 marmoset cytochrome P450 (P450) isoforms in the P450 1-51 subfamilies have been identified and characterized by the application of genome analysis and molecular functional characterization. In this mini-review, we provide an overview of the genomic structures, sequence identities, and substrate selectivities of marmoset P450s compared with those of human P450s. Based on the sequence identity, phylogeny, and genomic organization of marmoset P450s, orthologous relationships were established between human and marmoset P450s. Twenty-four members of the marmoset P450 1A, 2A, 2B, 2C, 2D, 2E, 3A, 4A, and 4F subfamilies shared high degrees of homology in terms of cDNA (>89%) and amino acid sequences (>85%) with the corresponding human P450s; P450 2C76 was among the exceptions. Phylogenetic analysis using amino acid sequences revealed that marmoset P450s in the P450 1-51 families were located in the same clades as their human and macaque P450 homologs. This finding underlines the evolutionary closeness of marmoset P450s to their human and macaque homologs. Most marmoset P450 1-4 enzymes catalyzed the typical drug-metabolizing reactions of the corresponding human P450 homologs, except for some differences of P450 2A6 and 2B6. Consequently, it appears that the substrate specificities of enzymes in the P450 1-4 families are generally similar in marmosets and humans. The information presented here supports a better understanding of the functional characteristics of marmoset P450s and their similarities and differences with human P450s. It is hoped that this mini-review will facilitate the successful use of marmosets as primate models in drug metabolism and pharmacokinetic studies.

Transcriptome-Based Identification and Molecular Evolution of the Cytochrome P450 Genes and Expression Profiling under Dimethoate Treatment in Amur Stickleback (Pungitius sinensis)

Cytochrome P450s (CYPs) are a family of membrane-bound mono-oxygenase proteins, which are involved in cell metabolism and detoxification of various xenobiotic substances. In this study, we identified 58 putative CYP genes in Amur stickleback (Pungitius sinensis) based on the transcriptome sequencing. Conserved motif distribution suggested their functional relevance within each group. Some present recombination events have accelerated the evolution of this gene family. Moreover, a few positive selection sites were identified, which may have accelerated the functional divergence of this family of proteins. Expression patterns of these CYP genes were investigated and indicated that most were affected by dimethoate treatment, suggesting that CYPs were involved in the detoxication of dimethoate. This study will provide a foundation for the further functional investigation of CYP genes in fishes.

Is there adaptation in the human genome for taste perception and phase I biotransformation?

Background

During the modern human expansion, new environmental pressures may have driven adaptation, especially in genes related to the perception of ingested substances and their detoxification. Consequently, positive (adaptive) selection may have occurred in genes related to taste, and in those related to the CYP450 system due to its role in biotransformation of potentially toxic compounds. A total of 91 genes (taste receptors and CYP450 superfamily) have been studied using Hierarchical Boosting, a powerful combination of different selection tests, to detect signatures of recent positive selection in three continental human populations: Northern Europeans (CEU), East Asians (CHB) and Africans (YRI). Analyses have been refined with selection analyses of the 26 populations of 1000 Genomes Project Phase 3.

Results

Genes related to taste perception have not been positively selected in the three continental human populations. This finding suggests that, contrary to results of previous studies, different allele frequencies among populations in genes such as TAS2R38 and TAS2R16 are not due to positive selection but to genetic drift. CYP1 and CYP2 genes, also previously considered to be under positive selection, did not show signatures of selective sweeps. However, three genes belonging to the CYP450 system have been identified by the Hierarchical Boosting as positively selected: CYP3A4 and CYP3A43 in CEU, and CYP27A1 in CHB.

Conclusions

No main adaptive differences are found in known taste receptor genes among the three continental human populations studied. However, there are important genetic adaptations in the cytochrome P450 system related to the Out of Africa expansion of modern humans. We confirmed that CYP3A4 and CYP3A43 are under selection in CEU, and we report for the first time CYP27A1 to be under positive selection in CHB.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1366-7) contains supplementary material, which is available to authorized users.

Hepatic expression of cytochrome P450 enzymes in non-human primate species

Cytochromes P450 (P450) largely remain to be characterized in great apes. Comparative immunochemical detection of drug metabolizing forms of P450s 1A, 2A, 2B, 2C, 2D, 2E, 2J, 3A, 4A, and 4F in liver microsomes from chimpanzees, gorillas, orangutans, gibbons, cynomolgus and rhesus macaques, and common marmosets were carried out.

Marmoset Cytochrome P450 3A4 Ortholog Expressed in Liver and Small-Intestine Tissues Efficiently Metabolizes Midazolam, Alprazolam, Nifedipine, and Testosterone

Common marmosets (Callithrix jacchus), small New World primates, are increasingly attracting attention as potentially useful animal models for drug development. However, characterization of cytochrome P450 (P450) 3A enzymes involved in the metabolism of a wide variety of drugs has not investigated in marmosets. In this study, sequence homology, tissue distribution, and enzymatic properties of marmoset P450 3A4 ortholog, 3A5 ortholog, and 3A90 were investigated. Marmoset P450 3A forms exhibited high amino acid sequence identities (88-90%) to the human and cynomolgus monkey P450 3A orthologs and evolutionary closeness to human and cynomolgus monkey P450 3A orthologs compared with other P450 3A enzymes. Among the five marmoset tissues examined, P450 3A4 ortholog mRNA was abundant in livers and small intestines where P450 3A4 ortholog proteins were immunologically detected. Three marmoset P450 3A proteins heterologously expressed in Escherichia coli membranes catalyzed midazolam 1'- and 4-hydroxylation, alprazolam 4-hydroxylation, nifedipine oxidation, and testosterone 6β-hydroxylation, similar to cynomolgus monkey and human P450 3A enzymes. Among the marmoset P450 3A enzymes, P450 3A4 ortholog effectively catalyzed midazolam 1'-hydroxylation, comparable to microsomes from marmoset livers and small intestines. Correlation analyses with 23 individual marmoset liver microsomes suggested contributions of P450 3A enzymes to 1'-hydroxylation of both midazolam (human P450 3A probe) and bufuralol (human P450 2D6 probe), similar to cynomolgus monkey P450 3A enzymes. These results indicated that marmoset P450 3A forms had functional characteristics roughly similar to cynomolgus monkeys and humans in terms of tissue expression patterns and catalytic activities, suggesting marmosets as suitable animal models for P450 3A-dependent drug metabolism.

Early Neolithic genomes from the eastern Fertile Crescent

We sequenced Early Neolithic genomes from the Zagros region of Iran (eastern Fertile Crescent), where some of the earliest evidence for farming is found, and identify a previously uncharacterized population that is neither ancestral to the first European farmers nor has contributed substantially to the ancestry of modern Europeans. These people are estimated to have separated from Early Neolithic farmers in Anatolia some 46,000 to 77,000 years ago and show affinities to modern-day Pakistani and Afghan populations, but particularly to Iranian Zoroastrians. We conclude that multiple, genetically differentiated hunter-gatherer populations adopted farming in southwestern Asia, that components of pre-Neolithic population structure were preserved as farming spread into neighboring regions, and that the Zagros region was the cradle of eastward expansion.

Qualitative de novo analysis of full length cDNA and quantitative analysis of gene expression for common marmoset (Callithrix jacchus) transcriptomes using parallel long-read technology and short-read sequencing

The common marmoset (Callithrix jacchus) is a non-human primate that could prove useful as human pharmacokinetic and biomedical research models. The cytochromes P450 (P450s) are a superfamily of enzymes that have critical roles in drug metabolism and disposition via monooxygenation of a broad range of xenobiotics; however, information on some marmoset P450s is currently limited. Therefore, identification and quantitative analysis of tissue-specific mRNA transcripts, including those of P450s and flavin-containing monooxygenases (FMO, another monooxygenase family), need to be carried out in detail before the marmoset can be used as an animal model in drug development. De novo assembly and expression analysis of marmoset transcripts were conducted with pooled liver, intestine, kidney, and brain samples from three male and three female marmosets. After unique sequences were automatically aligned by assembling software, the mean contig length was 718 bp (with a standard deviation of 457 bp) among a total of 47,883 transcripts. Approximately 30% of the total transcripts were matched to known marmoset sequences. Gene expression in 18 marmoset P450- and 4 FMO-like genes displayed some tissue-specific patterns. Of these, the three most highly expressed in marmoset liver were P450 2D-, 2E-, and 3A-like genes. In extrahepatic tissues, including brain, gene expressions of these monooxygenases were lower than those in liver, although P450 3A4 (previously P450 3A21) in intestine and P450 4A11- and FMO1-like genes in kidney were relatively highly expressed. By means of massive parallel long-read sequencing and short-read technology applied to marmoset liver, intestine, kidney, and brain, the combined next-generation sequencing analyses reported here were able to identify novel marmoset drug-metabolizing P450 transcripts that have until now been little reported. These results provide a foundation for mechanistic studies and pave the way for the use of marmosets as model animals for drug development in the future.

De novo sequence analysis of cytochrome P450 1-3 genes expressed in ostrich liver with highest expression of CYP2G19

The cytochrome P450 (CYP) 1-3 families are involved in xenobiotic metabolism, and are expressed primarily in the liver. Ostriches (Struthio camelus) are members of Palaeognathae with the earliest divergence from other bird lineages. An understanding of genes coding for ostrich xenobiotic metabolizing enzyme contributes to knowledge regarding the xenobiotic metabolisms of other Palaeognathae birds. We investigated CYP1-3 genes expressed in female ostrich liver using a next-generation sequencer. We detected 10 CYP genes: CYP1A5, CYP2C23, CYP2C45, CYP2D49, CYP2G19, CYP2W2, CYP2AC1, CYP2AC2, CYP2AF1, and CYP3A37. We compared the gene expression levels of CYP1A5, CYP2C23, CYP2C45, CYP2D49, CYP2G19, CYP2AF1, and CYP3A37 in ostrich liver and determined that CYP2G19 exhibited the highest expression level. The mRNA expression level of CYP2G19 was approximately 2-10 times higher than those of other CYP genes. The other CYP genes displayed similar expression levels. Our results suggest that CYP2G19, which has not been a focus of previous bird studies, has an important role in ostrich xenobiotic metabolism.

Toxic effects of dietary exposure to T-2 toxin on intestinal and hepatic biotransformation enzymes and drug transporter systems in broiler chickens

The effects of the mycotoxin T-2 on hepatic and intestinal drug-metabolizing enzymes (cytochrome P450) and drug transporter systems (MDR1 and MRP2) in poultry were investigated during this study. Broiler chickens received either uncontaminated feed, feed contaminated with 68μg/kg or 752μg/kg T-2 toxin. After 3weeks, the animals were euthanized and MDR1, MRP2, CYP1A4, CYP1A5 and CYP3A37 mRNA expression were analyzed using qRT-PCR. Along the entire length of the small intestine no significant differences were observed. In the liver, genes coding for CYP1A4, CYP1A5 and CYP3A37 were significantly down-regulated in the group exposed to 752μg/kg T-2. For CYP1A4, even a contamination level of 68μg/kg T-2 caused a significant decrease in mRNA expression. Expression of MDR1 was not significantly decreased in the liver. In contrast, hepatic MRP2 expression was significantly down-regulated after exposure to 752μg/kg T-2. Hepatic and intestinal microsomes were prepared to test the enzymatic activity of CYP3A. In the ileum and liver CYP3A activity was significantly increased in the group receiving 752μg/kg T-2 compared to the control group. The results of this study show that drug metabolizing enzymes and drug transporter mechanisms can be influenced due to prolonged exposure to relevant doses of T-2.

Rapid birth-and-death evolution of the xenobiotic metabolizing NAT gene family in vertebrates with evidence of adaptive selection

Background

The arylamine N-acetyltransferases (NATs) are a unique family of enzymes widely distributed in nature that play a crucial role in the detoxification of aromatic amine xenobiotics. Considering the temporal changes in the levels and toxicity of environmentally available chemicals, the metabolic function of NATs is likely to be under adaptive evolution to broaden or change substrate specificity over time, making NATs a promising subject for evolutionary analyses. In this study, we trace the molecular evolutionary history of the NAT gene family during the last ~450 million years of vertebrate evolution and define the likely role of gene duplication, gene conversion and positive selection in the evolutionary dynamics of this family.

Results

A phylogenetic analysis of 77 NAT sequences from 38 vertebrate species retrieved from public genomic databases shows that NATs are phylogenetically unstable genes, characterized by frequent gene duplications and losses even among closely related species, and that concerted evolution only played a minor role in the patterns of sequence divergence. Local signals of positive selection are detected in several lineages, probably reflecting response to changes in xenobiotic exposure. We then put a special emphasis on the study of the last ~85 million years of primate NAT evolution by determining the NAT homologous sequences in 13 additional primate species. Our phylogenetic analysis supports the view that the three human NAT genes emerged from a first duplication event in the common ancestor of Simiiformes, yielding NAT1 and an ancestral NAT gene which in turn, duplicated in the common ancestor of Catarrhini, giving rise to NAT2 and the NATP pseudogene. Our analysis suggests a main role of purifying selection in NAT1 protein evolution, whereas NAT2 was predicted to mostly evolve under positive selection to change its amino acid sequence over time. These findings are consistent with a differential role of the two human isoenzymes and support the involvement of NAT1 in endogenous metabolic pathways.

Conclusions

This study provides unequivocal evidence that the NAT gene family has evolved under a dynamic process of birth-and-death evolution in vertebrates, consistent with previous observations made in fungi.

The cytochrome P450 genesis locus: the origin and evolution of animal cytochrome P450s

The neighbourhoods of cytochrome P450 (CYP) genes in deuterostome genomes, as well as those of the cnidarians Nematostella vectensis and Acropora digitifera and the placozoan Trichoplax adhaerens were examined to find clues concerning the evolution of CYP genes in animals. CYP genes created by the 2R whole genome duplications in chordates have been identified. Both microsynteny and macrosynteny were used to identify genes that coexisted near CYP genes in the animal ancestor. We show that all 11 CYP clans began in a common gene environment. The evidence implies the existence of a single locus, which we term the 'cytochrome P450 genesis locus', where one progenitor CYP gene duplicated to create a tandem set of genes that were precursors of the 11 animal CYP clans: CYP Clans 2, 3, 4, 7, 19, 20, 26, 46, 51, 74 and mitochondrial. These early CYP genes existed side by side before the origin of cnidarians, possibly with a few additional genes interspersed. The Hox gene cluster, WNT genes, an NK gene cluster and at least one ARF gene were close neighbours to this original CYP locus. According to this evolutionary scenario, the CYP74 clan originated from animals and not from land plants nor from a common ancestor of plants and animals. The CYP7 and CYP19 families that are chordate-specific belong to CYP clans that seem to have originated in the CYP genesis locus as well, even though this requires many gene losses to explain their current distribution. The approach to uncovering the CYP genesis locus overcomes confounding effects because of gene conversion, sequence divergence, gene birth and death, and opens the way to understanding the biodiversity of CYP genes, families and subfamilies, which in animals has been obscured by more than 600 Myr of evolution.

Inhibition of cytochrome P450 enzymes involved in ketamine metabolism by use of liver microsomes and specific cytochrome P450 enzymes from horses, dogs, and humans

OBJECTIVE

To identify and characterize cytochrome P450 enzymes (CYPs) responsible for the metabolism of racemic ketamine in 3 mammalian species in vitro by use of chemical inhibitors and antibodies.

SAMPLE

Human, canine, and equine liver microsomes and human single CYP3A4 and CYP2C9 and their canine orthologs.

PROCEDURES

Chemical inhibitors selective for human CYP enzymes and anti-CYP antibodies were incubated with racemic ketamine and liver microsomes or specific CYPs. Ketamine N-demethylation to norketamine was determined via enantioselective capillary electrophoresis.

RESULTS

The general CYP inhibitor 1-aminobenzotriazole almost completely blocked ketamine metabolism in human and canine liver microsomes but not in equine microsomes. Chemical inhibition of norketamine formation was dependent on inhibitor concentration in most circumstances. For all 3 species, inhibitors of CYP3A4, CYP2A6, CYP2C19, CYP2B6, and CYP2C9 diminished N-demethylation of ketamine. Anti-CYP3A4, anti-CYP2C9, and anti-CYP2B6 antibodies also inhibited ketamine N-demethylation. Chemical inhibition was strongest with inhibitors of CYP2A6 and CYP2C19 in canine and equine microsomes and with the CYP3A4 inhibitor in human microsomes. No significant contribution of CYP2D6 to ketamine biotransformation was observed. Although the human CYP2C9 inhibitor blocked ketamine N-demethylation completely in the canine ortholog CYP2C21, a strong inhibition was also obtained by the chemical inhibitors of CYP2C19 and CYP2B6. Ketamine N-demethylation was stereoselective in single human CYP3A4 and canine CYP2C21 enzymes.

CONCLUSIONS AND CLINICAL RELEVANCE

Human-specific inhibitors of CYP2A6, CYP2C19, CYP3A4, CYP2B6, and CYP2C9 diminished ketamine N-demethylation in dogs and horses. To address drug-drug interactions in these animal species, investigations with single CYPs are needed.

Pregnane X receptor and yin yang 1 contribute to the differential tissue expression and induction of CYP3A5 and CYP3A4

The hepato-intestinal induction of the detoxifying enzymes CYP3A4 and CYP3A5 by the xenosensing pregnane X receptor (PXR) constitutes a key adaptive response to oral drugs and dietary xenobiotics. In contrast to CYP3A4, CYP3A5 is additionally expressed in several, mostly steroidogenic organs, which creates potential for induction-driven disturbances of the steroid homeostasis. Using cell lines and mice transgenic for a CYP3A5 promoter we demonstrate that the CYP3A5 expression in these organs is non-inducible and independent from PXR. Instead, it is enabled by the loss of a suppressing yin yang 1 (YY1)-binding site from the CYP3A5 promoter which occurred in haplorrhine primates. This YY1 site is conserved in CYP3A4, but its inhibitory effect can be offset by PXR acting on response elements such as XREM. Taken together, the loss of YY1 binding site from promoters of the CYP3A5 gene lineage during primate evolution may have enabled the utilization of CYP3A5 both in the adaptive hepato-intestinal response to xenobiotics and as a constitutively expressed gene in other organs. Our results thus constitute a first description of uncoupling induction from constitutive expression for a major detoxifying enzyme. They also suggest an explanation for the considerable tissue expression differences between CYP3A5 and CYP3A4.

Effect of estradiol on gene expression profile in cynomolgus macaque liver: implications for drug-metabolizing enzymes

Estrogen regulation of gene expression is essential for physiological function of estrogen-responsive tissues, such as mammary glands, ovaries, and the uterus. In the liver, estrogen is responsible for sex-dependent gene expression of drug-metabolizing enzymes in rodents. However, the influence of estrogen on hepatic gene expression has not been fully investigated in primates, including human. Macaque, including cynomolgus macaque, is an important species for comparative studies aimed at understanding human physiology due to its evolutionary closeness to human. To identify estrogen-responsive genes in primate liver, therefore, hepatic gene expression was compared, by microarray analysis, in ovariectomized cynomolgus macaques treated with estradiol or solvent (control). The analysis identified 98 estradiol-responsive genes; 47 and 51 were up- and down-regulated by estradiol, respectively (≥2.0-fold, P < 0.05). Expression of drug-metabolizing enzyme genes was also influenced by estradiol treatment; estradiol enhanced expression of GSTM5 (3.8-fold, P < 0.05) and CYP3A8(4) (2.7-fold, P < 0.01), but lowered expression of CYP4F12 (2.2-fold, P < 0.01), as verified by quantitative polymerase chain reaction. In particular, CYP3A8(4), orthologous to human CYP3A4, is an essential drug-metabolizing enzyme in cynomolgus macaque liver. These results suggest that expression of hepatic genes, including drug-metabolizing enzyme genes, is at least partly regulated by estradiol in cynomolgus macaque.

Species extrapolation for the 21st century

Safety factors are used in ecological risk assessments to extrapolate from the toxic responses of laboratory test species to all species representing that group in the environment. More accurate extrapolation of species responses is important. Advances in understanding the mechanistic basis for toxicological responses and identifying molecular response pathways can provide a basis for extrapolation across species and, in part, an explanation for the variability in whole organism responses to toxicants. We highlight potential short- and medium-term development goals to meet our long-term aspiration of truly predictive in silico extrapolation across wildlife species' response to toxicants. A conceptual approach for considering cross-species extrapolation is presented. Critical information is required to establish evidence-based species extrapolation, including identification of critical molecular pathways and regulatory networks that are linked to the biological mode of action and species' homologies. A case study is presented that examines steroidogenesis inhibition in fish after exposure to fadrozole or prochloraz. Similar effects for each compound among fathead minnow, medaka, and zebrafish were attributed to similar inhibitor pharmacokinetic/pharmacodynamic distributions and sequences of cytochrome P45019A1/2 (CYP19A1/2). Rapid advances in homology modeling allow the prediction of interactions of chemicals with enzymes, for example, CYP19 aromatase, which would eventually allow a prediction of potential aromatase toxicity of new compounds across a range of species. Eventually, predictive models will be developed to extrapolate across species, although substantial research is still required. Knowledge gaps requiring research include defining differences in life histories (e.g., reproductive strategies), understanding tissue-specific gene expression, and defining the role of metabolism on toxic responses and how these collectively affect the power of interspecies extrapolation methods.

Molecular evolution and functional divergence of the cytochrome P450 3 (CYP3) Family in Actinopterygii (ray-finned fish)

Background

The cytochrome P450 (CYP) superfamily is a multifunctional hemethiolate enzyme that is widely distributed from Bacteria to Eukarya. The CYP3 family contains mainly the four subfamilies CYP3A, CYP3B, CYP3C and CYP3D in vertebrates; however, only the Actinopterygii (ray-finned fish) have all four subfamilies and detailed understanding of the evolutionary relationship of Actinopterygii CYP3 family members would be valuable.

Methods and Findings

Phylogenetic relationships were constructed to trace the evolutionary history of the Actinopterygii CYP3 family genes. Selection analysis, relative rate tests and functional divergence analysis were combined to interpret the relationship of the site-specific evolution and functional divergence in the Actinopterygii CYP3 family. The results showed that the four CYP3 subfamilies in Actinopterygii might be formed by gene duplication. The first gene duplication event was responsible for divergence of the CYP3B/C clusters from ancient CYP3 before the origin of the Actinopterygii, which corresponded to the fish-specific whole genome duplication (WGD). Tandem repeat duplication in each of the homologue clusters produced stable CYP3B, CYP3C, CYP3A and CYP3D subfamilies. Acceleration of asymmetric evolutionary rates and purifying selection together were the main force for the production of new subfamilies and functional divergence in the new subset after gene duplication, whereas positive selection was detected only in the retained CYP3A subfamily. Furthermore, nearly half of the functional divergence sites appear to be related to substrate recognition, which suggests that site-specific evolution is closely related with functional divergence in the Actinopterygii CYP3 family.

Conclusions

The split of fish-specific CYP3 subfamilies was related to the fish-specific WGD, and site-specific acceleration of asymmetric evolutionary rates and purifying selection was the main force for the origin of the new subfamilies and functional divergence in the new subset after gene duplication. Site-specific evolution in substrate recognition was related to functional divergence in the Actinopterygii CYP3 family.

Identification and developmental expression of the full complement of Cytochrome P450 genes in Zebrafish

Background

Increasing use of zebrafish in drug discovery and mechanistic toxicology demands knowledge of cytochrome P450 (CYP) gene regulation and function. CYP enzymes catalyze oxidative transformation leading to activation or inactivation of many endogenous and exogenous chemicals, with consequences for normal physiology and disease processes. Many CYPs potentially have roles in developmental specification, and many chemicals that cause developmental abnormalities are substrates for CYPs. Here we identify and annotate the full suite of CYP genes in zebrafish, compare these to the human CYP gene complement, and determine the expression of CYP genes during normal development.

Results

Zebrafish have a total of 94 CYP genes, distributed among 18 gene families found also in mammals. There are 32 genes in CYP families 5 to 51, most of which are direct orthologs of human CYPs that are involved in endogenous functions including synthesis or inactivation of regulatory molecules. The high degree of sequence similarity suggests conservation of enzyme activities for these CYPs, confirmed in reports for some steroidogenic enzymes (e.g. CYP19, aromatase; CYP11A, P450scc; CYP17, steroid 17a-hydroxylase), and the CYP26 retinoic acid hydroxylases. Complexity is much greater in gene families 1, 2, and 3, which include CYPs prominent in metabolism of drugs and pollutants, as well as of endogenous substrates. There are orthologous relationships for some CYP1 s and some CYP3 s between zebrafish and human. In contrast, zebrafish have 47 CYP2 genes, compared to 16 in human, with only two (CYP2R1 and CYP2U1) recognized as orthologous based on sequence. Analysis of shared synteny identified CYP2 gene clusters evolutionarily related to mammalian CYP2 s, as well as unique clusters.

Conclusions

Transcript profiling by microarray and quantitative PCR revealed that the majority of zebrafish CYP genes are expressed in embryos, with waves of expression of different sets of genes over the course of development. Transcripts of some CYP occur also in oocytes. The results provide a foundation for the use of zebrafish as a model in toxicological, pharmacological and chemical disease research.

Species differences in the pharmacology and toxicology of PEGylated helper-dependent adenovirus

Clinically relevant doses of helper-dependent adenoviruses (HDAds) provoke the host response against capsid proteins in primates and rodents. To determine if PEGylation truly affects this, baboons and mice were given either HDAd or PEG-HDAd expressing beta-galactosidase at 5 × 10¹¹ or 3 × 10¹² virus particles per kilogram (vp/kg) by iv infusion. Serum cytokines and blood chemistries were assessed for 96 h. PEG-HDAd reduced IL-6 6-fold in mice and 3-fold in the primate. This vector reduced IL-12 by 50% in both animal models. PEGylation reduced serum transaminases by approximately 50% at each dose in the primate and the mouse. PEGylation did not alter hepatic transduction efficiency in the mouse but did reduce transduction efficiency in the liver and the spleen of primates. Unmodified and PEGylated virus suppressed hepatic CYP3A activity in both animal models. PEGylation doubled the half-life (t(½)) of the virus in the mouse and cut plasma clearance (CL) in half without affecting the half-life in primates. These results suggest that there are notable species-specific differences in the biodistribution of and response to PEG-modified vectors which may be linked to differences in binding properties to coagulation factors, receptor density and tissue architecture in the liver.

Cytochrome P450-mediated 17beta-estradiol metabolism in zebrafish (Danio rerio)

Cytochrome P4501 (CYP1) and CYP3A proteins are primarily responsible for the metabolism of 17beta-estradiol (E(2)) in mammals. We have cloned and heterologously expressed CYP1A, CYP1B1, CYP1C1, CYP1C2, CYP1D1, and CYP3A65 from zebrafish (Danio rerio) to determine the CYP-mediated metabolism of E(2) in a non-mammalian species. Constructs of each CYP cDNA were created using a leader sequence from the bacterial ompA gene to allow appropriate expression in Escherichia coli without 5' modification of the gene. Membrane vesicles were purified, and functional CYP protein was verified using carbon monoxide difference spectra and fluorescent catalytic assays with the substrates 7-ethoxyresorufin and 7-benzyloxy-4-(trifluoromethyl)-coumarin. Rates of in vitro E(2) metabolism into 4-hydroxyE(2) (4-OHE(2)), 2-hydroxyE(2) (2-OHE(2)), and 16alpha-hydroxyE(1) (16alpha-OHE(1)) metabolites were determined by gas chromatography/mass spectrometry. The 2-OHE(2) metabolite was produced by all CYPs tested, while 4-OHE(2) was only detected following incubation with CYP1A, CYP1B1, CYP1C1, and CYP1C2. The 16alpha-OHE(1) metabolite was only produced by CYP1A. The highest rates of E(2) metabolism were from CYP1A and CYP1C1, followed by CYP1C2. CYP1B1, CYP1D1, and CYP3A65 had low rates of E(2) metabolism. E(2) metabolism by zebrafish CYP1A, CYP1C1, and CYP1C2 produced similar ratios of 4-OHE(2) to 2-OHE(2) as previous studies with mammalian CYP1As. CYP1B1 formed the highest ratio of 4-OHE(2) to 2-OHE(2) metabolites. Contrary to mammals, these results suggest that fish CYP1A and CYP1C proteins are primarily responsible for E(2) metabolism, with only minor contributions from CYP3A65 and CYP1B1. Similar to mammals, 2-OHE(2) is the predominant metabolite from CYP-mediated E(2) metabolism in fish, suggesting that all vertebrate species produce the same major E(2) metabolite.

Molecular population genetics of human CYP3A locus: signatures of positive selection and implications for evolutionary environmental medicine

BACKGROUND

The human CYP3A gene cluster codes for cytochrome P450 (CYP) subfamily enzymes that catalyze the metabolism of various exogenous and endogenous chemicals and is an obvious candidate for evolutionary and environmental genomic study. Functional variants in the CYP3A locus may have undergone a selective sweep in response to various environmental conditions.

OBJECTIVE

The goal of this study was to profile the allelic structure across the human CYP3A locus and investigate natural selection on that locus.

METHODS

From the CYP3A locus spanning 231 kb, we resequenced 54 genomic DNA fragments (a total of 43,675 bases) spanning four genes (CYP3A4, CYP3A5, CYP3A7, and CYP3A43) and two pseudogenes (CYP3AP1 and CYP3AP2), and randomly selected intergenic regions at the CYP3A locus in Africans (24 individuals), Caucasians (24 individuals), and Chinese (29 individuals). We comprehensively investigated the nucleotide diversity and haplotype structure and examined the possible role of natural selection in shaping the sequence variation throughout the gene cluster.

RESULTS

Neutrality tests with Tajima's D, Fu and Li's D* and F*, and Fay and Wu's H indicated possible roles of positive selection on the entire CYP3A locus in non-Africans. Sliding-window analyses of nucleotide diversity and frequency spectrum, as well as haplotype diversity and phylogenetically inferred haplotype structure, revealed that CYP3A4 and CYP3A7 had recently undergone or were undergoing a selective sweep in all three populations, whereas CYP3A43 and CYP3A5 were undergoing a selective sweep in non-Africans and Caucasians, respectively.

CONCLUSION

The refined allelic architecture and selection spectrum for the human CYP3A locus highlight that evolutionary dynamics of molecular adaptation may underlie the phenotypic variation of the xenobiotic disposition system and varied predisposition to complex disorders in which xenobiotics play a role.

Ligand diversity of human and chimpanzee CYP3A4: activation of human CYP3A4 by lithocholic acid results from positive selection

For currently unknown reasons, the evolution of CYP3A4 underwent acceleration in the human lineage after the split from chimpanzee. We investigated the significance of this event by comparing Escherichia coli-expressed CYP3A4 from humans, chimpanzee, and their most recent common ancestor. The expression level of chimpanzee CYP3A4 was approximately 50% of the human CYP3A4, whereas ancestral CYP3A4 did not express in E. coli. Steady-state kinetic analysis with 7-benzyloxyquinoline, 7-benzyloxy-4-(trifluoromethyl)coumarin (7-BFC), and testosterone showed no significant differences between human and chimpanzee CYP3A4. Upon addition of alpha-naphthoflavone (25 microM), human CYP3A4 showed a slightly decreased substrate concentration at which 50% of the maximal rate V(max) is reached for 7-BFC, whereas chimpanzee CYP3A4 showed a >2-fold increase. No significant differences in inhibition/activation were found for a panel of 43 drugs and endogenous compounds, suggesting that the wide substrate spectrum of human CYP3A4 precedes the human-chimpanzee split. A striking exception was the hepatotoxic secondary bile acid lithocholic acid, which at saturation caused a 5-fold increase in 7-BFC debenzylation by human CYP3A4 but not by chimpanzee CYP3A4. Mutagenesis of human CYP3A4 revealed that at least four of the six amino acids positively selected in the human lineage contribute to the activating effect of lithocholic acid. In summary, the wide functional conservation between chimpanzee and human CYP3A4 raises the prospect that phylogenetically more distant primate species such as rhesus and squirrel monkey represent suitable models of the human counterpart. Positive selection on the human CYP3A4 may have been triggered by an increased load of dietary steroids, which led to a novel defense mechanism against cholestasis.