The Human CYP2C Locus: A Prototype for Intergenic and Exon Repetition Splicing Events

In human there are four known CYP2C genes that have been mapped to chromosome 10q24 with the order Cen–2C18–2C19–2C9–2C8–Tel. Previously we have shown that splicing events joining exons from the neighboring 2C18 and 2C19 genes occur in human liver and epidermis. Here evidence is presented that the terminal genes of this cluster, 2C18 and 2C8, are also involved in intergenic splicing. Most interestingly, several of these 2C18/2C8 RNAs were composed of all nine exons, thus conceivably having the potential for coding functional proteins. Moreover, chimeric RNA species consisting of exons originating not only from the CYP2C8 and CYP2C18 genes, but also from the CYP2C19 gene were detected. In all cases the exons from the different CYP2C genes were joined at the correct canonical splice sites. However, the closely linked RBP4 gene is not participating in intergenic splicing with the CYP2C genes. In addition, CYP2C8 gene expression was found to generate a variety of scrambled RNA molecules including species that contained repetitions of certain exons.     

Intergenic transcripts containing a novel human cytochrome P450 2C exon 1 spliced to sequences from the CYP2C9 gene

The cytochrome P450 2C (CYP2C) gene locus was found to includea novel exon 1 sequence with high similarity to the canonicalexon 1 of CYP2C18. Rapid amplification of cDNA ends (RACE) andPCR amplifications of human liver cDNA revealed the presenceof several intergenic species containing the CYP2C18 exon 1–likesequence spliced to different combinations of exonic and intronicsequences from the CYP2C9 gene. One splice variant was foundto have an open reading frame starting at the canonical translationinitiation codon of the CYP2C18 exon 1–like sequence.Another variant consisted of the nine typical CYP2C9 exons splicedafter the CYP2C18 exon 1–like sequence through a segmentof CYP2C9 5' flanking sequences. Moreover, analysis of bacterialartificial chromosome (BAC) clones revealed that the CYP2C18exon 1–like sequence was located in the intergenic regionbetween the CYP2C19 and CYP2C9 genes. The finding that a solitaryexon is spliced with sequences from a neighboring gene may beinterpreted as representing a general evolutionary mechanismaimed at using the full expression potential of a cell's genomicinformational content.     

Cloning and characterization of a novel gene (C8orf2), a human representative of a novel gene family with homology to C. elegans C42.C1.9.

Large-scale sequencing of selected genomic regions, coupled with in silico gene trapping, is a robust approach to identifying previously unknown genes. In this way we have found a gene (C8orf2) that is highly homologous to C. elegans C42C1.9. C8orf2 was situated on 8p11. 2 between STS markers NIB1979 (proximal) and AFMA295ZD5 (distal), oriented toward the centromere. C8orf2 consisted of 16 exons spanning more than 16.5 kb of genomic DNA, and was expressed ubiquitously in human tissues. The gene encoded 339-and 152-amino acid polypeptides by alternative splicing; the larger variant contained a region extremely rich in charged amino acids, in particular lysine and glutamic acid. C8orf2 also bore sequence homology to the human KE04p gene. Its conservation among highly divergent species suggests that C8orf2 belongs to a novel gene family.     

Mapping and characterization of the mouse and human SS18 genes, two human SS18-like genes and a mouse Ss18 pseudogene.

We have previously isolated and characterized a mouse cDNA orthologous to the human synovial sarcoma associated SS18 (formerly named SSXT and SYT) cDNA. Here, we report the characterization of the genomic structure of the mouse Ss18 gene. Through in silico methods with sequence information contained in the public databases, we did the same for the human SS18 gene and two human SS18 homologous genes, SS18L1 and SS18L2. In addition, we identified a mouse Ss18 processed pseudogene and mapped it to chromosome 1, band A2-3. The mouse Ss18 gene, which is subject to extensive alternative splicing, is made up of 11 exons, spread out over approximately 45 kb of genomic sequence. The human SS18 gene is also composed of 11 exons with similar intron-exon boundaries, spreading out over about 70 kb of genomic sequence. One alternatively spliced exon, which is not included in the published SS18 cDNA, corresponds to a stretch of sequence which we previously identified in the mouse Ss18 cDNA. The human SS18L1 gene, which is also made up of 11 exons with similar intron-exon boundaries, was mapped to chromosome 20 band q13.3. The smaller SS18L2 gene, which is composed of three exons with similar boundaries as the first three exons of the other three genes, was mapped to chromosome 3 band p21. Through sequence and mutation analyses this gene could be excluded as a candidate gene for 3p21-associated renal cell cancer. In addition, we created a detailed BAC map around the human SS18 gene, placing it unequivocally between the CA-repeat marker AFMc014wf9 and the dihydrofolate reductase pseudogene DHFRP1. The next gene in this map, located distal to SS18, was found to be the TBP associated factor TAFII-105 (TAF2C2). Further analogies between the mouse Ss18 gene, the human SS18 gene and its two homologous genes were found in the putative promoter fragments. All four promoters resemble the promoters of housekeeping genes in that they are TATA-less and embedded in canonical CpG islands, thus explaining the high and widespread expression of the SS18 genes.     

Divergent structures of Caenorhabditis elegans cytochrome P450 genes suggest the frequent loss and gain of introns during the evolution of nematodes

The Caenorhabditis elegans genome contains more than 60 cytochrome P450 (CYP) genes. The exon-intron organizations of all of the available and potentially active C. elegans CYP genes were inferred by a newly developed program for predicting protein-coding exons based on the alignment of a genomic DNA sequence and a protein profile. From the predicted amino acid sequences, all of the C. elegans CYP genes except one were classified into three groups, which were closely related to the mammalian drug-metabolizing P450 gene families CYP2, CYP3, and CYP4. The gene structures were strikingly divergent within each group; 20, 10, and 5 unique gene organizations were identified among 40, 18, and 5 genes in the CYP2-, CYP3-, and CYP4-related groups, respectively. The degrees of divergence in gene organization were strongly correlated with those in the amino acid sequences of encoding proteins, and the minimum rate of change in an intron insertion site was estimated to be about 90 times less frequent than amino acid substitutions. Parsimonious analyses suggested that frequent loss and gain of introns has occurred during the evolution of CYP genes in each group after the divergence of nematodes, arthropods, and deuterostomia. Few, if any, incidents of intron sliding were evident, and a model that did not allow intron insertions was highly inconsistent with the observations. All of these findings are explained better by the intron-late view than by the intron-early view.     

Comparison of the substrate specificities of human liver cytochrome P450s 2C9 and 2C18: application to the design of a specific substrate of CYP 2C18.

A series of 2-aroylthiophenes derived from tienilic acid by replacement of its OCH2COOH substituent with groups bearing various functions have been synthesized and studied as possible substrates of recombinant human liver cytochrome P450s 2C9 and 2C18 expressed in yeast. Whereas only compounds bearing a negative charge acted as substrates of CYP 2C9 and were hydroxylated at position 5 of their thiophene ring at a significant rate, many neutral 2-aroylthiophenes were 5-hydroxylated by CYP 2C18 with kcat values of >2 min-1. Among the various compounds that were studied, those bearing an alcohol function were the best CYP 2C18 substrates. One of them, compound 3, which bears a terminal O(CH2)3OH function, appeared to be a particularly good substrate of CYP 2C18. It was regioselectively hydroxylated by CYP 2C18 at position 5 of its thiophene ring with a KM value of 9 +/- 1 microM and a kcat value of 125 +/- 25 min-1, which are the highest described so far for a CYP 2C. A comparison of the oxidations of 3, by yeast-expressed CYP 1A1, 1A2, 2C8, 2C9, 2C18, 2C19, 2D6, 2E1, 3A4, and 3A5, showed that only CYP 2C8, 2C18, and 2C19 were able to catalyze the 5-hydroxylation of 3. However, the catalytic efficiency of CYP 2C18 for that reaction was considerably higher (kcat/KM value being 3-4 orders of magnitude larger than those found for CYP 2C8 and 2C19). Several human P450s exhibited small activities for the oxidative O-dealkylation of 3. The four recombinant CYP 2Cs were the best catalysts for that reaction (kcat between 1 and 5 min-1) when compared to all the P450s that were tested, even though it is a minor reaction in the case of CYP 2C18. All these results show that compound 3 is a new, selective, and highly efficient substrate for CYP 2C18 that should be useful for the study of this P450 in various organs and tissues. They also suggest some key differences between the active sites of CYP 2C9 and CYP 2C18 for substrate recognition.     

Comparison of P450s from human and fugu: 420 million years of vertebrate P450 evolution

The fugu (pufferfish) genome has been sequenced, and a second genome assembly was released 17 May 2002. Exhaustive searches were made to identify all P450 genes and pseudogenes from the earlier release of 26 October 2001. P450 genes assembled as completely as possible from these data were used to do additional searches of the newer assembly and all P450 genes and pseudogenes in the available fugu sequence data have been identified, compared to human P450s, and assigned names. There are 54 P450 genes in fugu and 1 nearly intact pseudogene (CYP3A50P). CYP1A is missing much of its N-terminal half; however, 45 P450 genes are completely assembled. Eight others are lacking only one or two exons or less. CYP2X4 is known only from an EST. This may be a 55th P450 gene if it represents an accurate sequence. In addition to 2X4, there are 16 other pseudogene fragments or small pieces of P450 genes. At the P450 family level, 17 of 18 mammalian families are found in fugu. CYP39 is the only CYP family missing and it is not seen in any other fish sequence data either. The CYP2 family shows the largest degree of divergence. In the CYP2 family, only CYP2R1 and CYP2U1 are conserved as recognizable subfamilies across species. Intron-exon boundaries are largely preserved across 420 million years of evolution.     

Circular RNAs are abundant,conserved, and associated with ALU repeats

Circular RNAs composed of exonic sequence have been described in a small number of genes. Thought to result from splicing errors, circular RNA species possess no known function. To delineate the universe of endogenous circular RNAs, we performed high-throughput sequencing (RNA-seq) of libraries prepared from ribosome-depleted RNA with or without digestion with the RNA exonuclease, RNase R. We identified >25,000 distinct RNA species in human fibroblasts that contained non-colinear exons (a “backsplice”) and were reproducibly enriched by exonuclease degradation of linear RNA. These RNAs were validated as circular RNA (ecircRNA), rather than linear RNA, and were more stable than associated linear mRNAs in vivo. In some cases, the abundance of circular molecules exceeded that of associated linear mRNA by >10-fold. By conservative estimate, we identified ecircRNAs from 14.4% of actively transcribed genes in human fibroblasts. Application of this method to murine testis RNA identified 69 ecircRNAs in precisely orthologous locations to human circular RNAs. Of note, paralogous kinases HIPK2 and HIPK3 produce abundant ecircRNA from their second exon in both humans and mice. Though HIPK3 circular RNAs contain an AUG translation start, it and other ecircRNAs were not bound to ribosomes. Circular RNAs could be degraded by siRNAs and, therefore, may act as competing endogenous RNAs. Bioinformatic analysis revealed shared features of circularized exons, including long bordering introns that contained complementary ALU repeats. These data show that ecircRNAs are abundant, stable, conserved and nonrandom products of RNA splicing that could be involved in control of gene expression.