Cytochrome P450 4F isoforms in different species: Identification, gene regulation and putative roles in inflammation

The cytochrome P450 4F subfamily comprises a group of enzymes that metabolize derivatives of arachidonic acid such as prostaglandins, lipoxins leukotrienes and hydroxyeicosatetraenoic acids, which are important mediators involved in the inflammatory response. Therefore, we speculate that CYP4Fs might be able to modulate the extent of the inflammation by controlling of the tissue levels of these inflammatory mediators, especially, leukotriene B4. One way to provide support for this hypothesis is to test whether the expression of CYP4Fs changes under inflammatory conditions, since these changes are required to adjust the levels of inflammatory mediators. ^ A lipopolysacchride (LPS) induced rat inflammation model was used to analyze the expressions of rat CYP4F4 and CYP4F5 in liver and kidney. LPS administration did not change the constitutive expression level of CYP4F4 and CYP4F5. In liver, the expressions of CYP4F4 and CYP4F5 decreased to 50–60% of the untreated level. The same effect of LPS on CYP4F4 and CYP4F5 expression can be mimicked in hepatocyte primary cultures treated with LPS, indicating a direct of effect of LPS on hepatocytes. LPS treatment also decreased the activity of liver microsomes towards chlorpromazine, however, antibody inhibition study revealed that liver CYP4Fs are not the only players in metabolizing chlorpromazine. To study further the underlying mechanism, CYP4F5 gene was isolated, characterized, and the promoter region was defined. ^ Accumulating evidence showed that peroxisome proliferator-activated receptors (PPARs) play an active role in inflammation. To investigate the possible role of PPARα in regulating CYP4F expression by inflammation or by clofibrate treatment, the expressions of two new mouse 4F isoforms were analyzed in PPARα knockout mice upon LPS or clofibrate challenge. A novel induction of CYP4F15 by LPS and clofibrate was observed in kidney, and this effect is totally dependent on the presence of PPARα. Renal CYP4F16 expression was not affected by LPS or clofibrate in both (+/+) and (−/−) mice. In contrast, hepatic expressions of CYP4F15 and CYP4F16 were reduced significantly in (+/+) mice, but much less in (−/−) mice, suggesting that PPARα is partially responsible for this down-regulation. Clofibrate treatment reduced the expression of CYP4F16 in liver, but has no effect on CYP4F15 and PPARα does not have a role in hepatic CYP4F expression regulated by clofibrate. In general, CYP4Fs are regulated in an isoform-, tissue- and species-specific manner. ^ A human CYP4F isoform, CYP4F11, was isolated. The genomic structure was also solved by using database mining and bioinformatics tools. Localization of CYP4F11 to chromosome 19, 16 kb upstream of CYP4F2, suggests that human CYP4F genes may form a cluster on chromosome 19. This novel human 4F is highly expressed in liver, as well as in kidney, heart and skeletal muscle. Further study of the activity and gene regulation on CYP4F11 will provide us more insights into the physiological functions of CYP4F subfamily. ^

Towards the identification of a tumor suppressor gene on chromosome 3p12: Physical mapping and candidate gene screening

Renal cell carcinoma (RCC) is the most common malignant tumor of the kidney. Characterization of RCC tumors indicates that the most frequent genetic event associated with the initiation of tumor formation involves a loss of heterozygosity or cytogenetic aberration on the short arm of human chromosome 3. A tumor suppressor locus Nonpapillary Renal Carcinoma-1 (NRC-1, OMIM ID 604442) has been previously mapped to a 5–7 cM region on chromosome 3p12 and shown to induce rapid tumor cell death in vivo, as demonstrated by functional complementation experiments. ^ To identify the gene that accounts for the tumor suppressor activities of NRC-1, fine-scale physical mapping was conducted with a novel real-time quantitative PCR based method developed in this study. As a result, NRC-1 was mapped within a 4.6-Mb region defined by two unique sequences within UniGene clusters Hs.41407 and Hs.371835 (78,545Kb–83,172Kb in the NCBI build 31 physical map). The involvement of a putative tumor suppressor gene Robo1/Dutt1 was excluded as a candidate for NRC-1. Furthermore, a transcript map containing eleven candidate genes was established for the 4.6-Mb region. Analyses of gene expression patterns with real-time quantitative RT-PCR assays showed that one of the eleven candidate genes in the interval (TSGc28) is down-regulated in 15 out of 20 tumor samples compared with matched normal samples. Three exons of this gene have been identified by RACE experiments, although additional exon(s) seem to exist. Further gene characterization and functional studies are required to confirm the gene as a true tumor suppressor gene. ^ To study the cellular functions of NRC-1, gene expression profiles of three tumor suppressive microcell hybrids, each containing a functional copy of NRC-1, were compared with those of the corresponding parental tumor cell lines using 16K oligonucleotide microarrays. Differentially expressed genes were identified. Analyses based on the Gene Ontology showed that introduction of NRC-1 into tumor cell lines activates genes in multiple cellular pathways, including cell cycle, signal transduction, cytokines and stress response. NRC-1 is likely to induce cell growth arrest indirectly through WEE1. ^

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.

Identification of the block in targeted retroviral-mediated gene transfer

A chimeric retroviral vector (33E67) containing a CD33-specific single-chain antibody was generated in an attempt to target cells displaying the CD33 surface antigen. The chimeric envelope protein was translated, processed, and incorporated into viral particles as efficiently as wild-type envelope protein. The viral particles carrying the 33E67 envelope protein could bind efficiently to the CD33 receptor on target cells and were internalized, but no gene transfer occurred. A unique experimental approach was used to examine the basis for this postbinding block. Our data indicate that the chimeric envelope protein itself cannot participate in the fusion process, the most reasonable explanation being that this chimeric protein cannot undergo the appropriate conformational change that is thought to be triggered by receptor binding, a suggested prerequisite to subsequent fusion and core entry. These results indicate that the block to gene transfer in this system, and probably in most of the current chimeric retroviral vectors to date, is the inability of the chimeric envelope protein to undergo this obligatory conformational change.

Vibrio cholerae O1 El Tor: Identification of a gene cluster required for the rugose colony type, exopolysaccharide production, chlorine resistance, and biofilm formation

The rugose colony variant of Vibrio cholerae O1, biotype El Tor, is shown to produce an exopolysaccharide, EPSETr, that confers chlorine resistance and biofilm-forming capacity. EPSETr production requires a chromosomal locus, vps, that contains sequences homologous to carbohydrate biosynthesis genes of other bacterial species. Mutations within this locus yield chlorine-sensitive, smooth colony variants that are biofilm deficient. The biofilm-forming properties of EPSETr may enable the survival of V. cholerae O1 within environmental aquatic habitats between outbreaks of human disease.

Identification and characterization of the ARP1 gene, a target for the human acute leukemia ALL1 gene

ALL1, the human homologue of Drosophila trithorax, is directly involved in human acute leukemias associated with abnormalities at 11q23. Using the differential display method, we isolated a gene that is down-regulated in All1 double-knockout mouse embryonic stem (ES) cells. The gene, designated ARP1 (also termed RIEG, Ptx2, or Otlx2), is a member of a family of homeotic genes containing a short motif shared with several homeobox genes. Using a bacterially synthesized All1 polypeptide encompassing the AT-hook motifs, we identified a 0.5-kb ARP1 DNA fragment that preferentially bound to the polypeptide. Within this DNA, a region of ≈100 bp was protected by the polypeptide from digestion with ExoIII and DNase I. Whole-mount in situ hybridization to early mouse embryos of 9.5–10.5 days indicated a complex pattern of Arp1 expression spatially overlapping with the expression of All1. Although the ARP1 gene is expressed strongly in bone marrow cells, no transcripts were detected in six leukemia cell lines with 11q23 translocations. These results suggest that ARP1 is up-regulated by the All1 protein, possibly through direct interaction with an upstream DNA sequence of the former. The results are also consistent with the suggestion that ALL1 chimeric proteins resulting from 11q23 abnormalities act in a dominant negative fashion.

Identification of a new mammalian centrin gene, more closely related to Saccharomyces cerevisiae CDC31 gene

Among the numerous centrin isoforms identified by two-dimensional gel electrophoresis in human cells, an acidic and slow-migrating isoform is particularly enriched in a centrosome fraction. We report here that this isoform specifically reacts with antibodies raised against Saccharomyces cerevisiae Cdc31p and is present, as other centrin isoforms, in the distal lumen of centrioles. It is encoded by a new centrin gene, which we propose to name HsCEN3 (Homo sapiens centrin gene 3). This gene is more closely related to the yeast CDC31 gene, and shares less identity with algae centrin than HsCEN1 and HsCEN2. A murine CDC31-related gene was also found that shows 98% identity and 100% similarity with HsCEN3, demonstrating a higher interspecies conservation than the murine centrin gene MmCEN1 (Mus musculus centrin gene 1) with either HsCEN1, or HsCEN2. Finally, immunological data suggest that a CDC31-related gene could exist in amphibians and echinoderms as well. All together, our data suggest the existence of two divergent protein subfamilies in the current centrin family, which might be involved in distinct centrosome-associated functions. The possible implication of this new mammalian centrin gene in centrosome duplication is discussed.

Screening for imprinted genes by allelic message display: Identification of a paternally expressed gene Impact on mouse chromosome 18

A systematic screen termed the allelic message display (AMD) was developed for the hunting of imprinted genes. In AMD, differential display PCR is adopted to image allelic expression status of multiple polymorphic transcripts in two parental mouse strains, reciprocal F1 hybrids and pooled backcross progenies. From the displayed patterns, paternally and maternally expressed transcripts can be unequivocally identified. The effectiveness of AMD screening was clearly demonstrated by the identification of a paternally expressed gene Impact on mouse chromosome 18, the predicted product of which belongs to the YCR59c/yigZ hypothetical protein family composed of yeast and bacterial proteins with currently unknown function. In contrast with previous screening methods necessitating positional cloning efforts or generation of parthenogenetic embryos, this approach requires nothing particular but appropriately crossed mice and can be readily applied to any tissues at various developmental stages. Hence, AMD would considerably accelerate the identification of imprinted genes playing pivotal roles in mammalian development and the pathogenesis of various diseases.

Identification of a gene encoding an acyl CoA:diacylglycerol acyltransferase, a key enzyme in triacylglycerol synthesis

Triacylglycerols are quantitatively the most important storage form of energy for eukaryotic cells. Acyl CoA:diacylglycerol acyltransferase (DGAT, EC 2.3.1.20) catalyzes the terminal and only committed step in triacylglycerol synthesis, by using diacylglycerol and fatty acyl CoA as substrates. DGAT plays a fundamental role in the metabolism of cellular diacylglycerol and is important in higher eukaryotes for physiologic processes involving triacylglycerol metabolism such as intestinal fat absorption, lipoprotein assembly, adipose tissue formation, and lactation. DGAT is an integral membrane protein that has never been purified to homogeneity, nor has its gene been cloned. We identified an expressed sequence tag clone that shared regions of similarity with acyl CoA:cholesterol acyltransferase, an enzyme that also uses fatty acyl CoA as a substrate. Expression of a mouse cDNA for this expressed sequence tag in insect cells resulted in high levels of DGAT activity in cell membranes. No other acyltransferase activity was detected when a variety of substrates, including cholesterol, were used as acyl acceptors. The gene was expressed in all tissues examined; during differentiation of NIH 3T3-L1 cells into adipocytes, its expression increased markedly in parallel with increases in DGAT activity. The identification of this cDNA encoding a DGAT will greatly facilitate studies of cellular glycerolipid metabolism and its regulation.

Identification of differentially expressed mRNA in prokaryotic organisms by customized amplification libraries (DECAL): The effect of isoniazid on gene expression in Mycobacterium tuberculosis

Understanding the effects of the external environment on bacterial gene expression can provide valuable insights into an array of cellular mechanisms including pathogenesis, drug resistance, and, in the case of Mycobacterium tuberculosis, latency. Because of the absence of poly(A)+ mRNA in prokaryotic organisms, studies of differential gene expression currently must be performed either with large amounts of total RNA or rely on amplification techniques that can alter the proportional representation of individual mRNA sequences. We have developed an approach to study differences in bacterial mRNA expression that enables amplification by the PCR of a complex mixture of cDNA sequences in a reproducible manner that obviates the confounding effects of selected highly expressed sequences, e.g., ribosomal RNA. Differential expression using customized amplification libraries (DECAL) uses a library of amplifiable genomic sequences to convert total cellular RNA into an amplified probe for gene expression screens. DECAL can detect 4-fold differences in the mRNA levels of rare sequences and can be performed on as little as 10 ng of total RNA. DECAL was used to investigate the in vitro effect of the antibiotic isoniazid on M. tuberculosis, and three previously uncharacterized isoniazid-induced genes, iniA, iniB, and iniC, were identified. The iniB gene has homology to cell wall proteins, and iniA contains a phosphopantetheine attachment site motif suggestive of an acyl carrier protein. The iniA gene is also induced by the antibiotic ethambutol, an agent that inhibits cell wall biosynthesis by a mechanism that is distinct from isoniazid. The DECAL method offers a powerful new tool for the study of differential gene expression.

Identification and characterization of Saccharomyces cerevisiae EXO1, a gene encoding an exonuclease that interacts with MSH2

A two-hybrid screen was used to identify Saccharomyces cerevisiae genes encoding proteins that interact with MSH2. One gene was found to encode a homologue of Schizosaccharomyces pombe EXO1, a double-stranded DNA-specific 5′–3′ exonuclease. S. cerevisiae EXO1 interacted with both S. cerevisiae and human MSH2 in two-hybrid and coimmunoprecipitation experiments. exo1 mutants showed a mutator phenotype, and epistasis analysis was consistent with EXO1 functioning in the MSH2-dependent mismatch repair pathway. exo1 mutations were lethal in combination with rad27 mutations, and overexpression of EXO1 suppressed both the temperature sensitive and mutator phenotypes of rad27 mutants.

Identification and characterization of a retinoid-induced class II tumor suppressor/growth regulatory gene

Retinoids, synthetic and natural analogs of retinoic acid, exhibit potent growth inhibitory and cell differentiation activities that account for their beneficial effects in treating hyperproliferative diseases such as psoriasis, actinic keratosis, and certain neoplasias. Tazarotene is a synthetic retinoid that is used in the clinic for the treatment of psoriasis. To better understand the mechanism of retinoid action in the treatment of hyperproliferative diseases, we used a long-range differential display–PCR to isolate retinoid-responsive genes from primary human keratinocytes. We have identified a cDNA, tazarotene-induced gene 3 (TIG3; Retinoic Acid Receptor Responder 3) showing significant homology to the class II tumor suppressor gene, H-rev 107. Tazarotene treatment increases TIG3 expression in primary human keratinocytes and in vivo in psoriatic lesions. Increased TIG3 expression is correlated with decreased proliferation. TIG3 is expressed in a number of tissues, and expression is reduced in cancer cell lines and some primary tumors. In breast cancer cell lines, retinoid-dependent TIG3 induction is observed in lines that are growth suppressed by retinoids but not in nonresponsive lines. Transient over-expression of TIG3 in T47D or Chinese hamster ovary cells inhibits colony expansion. Finally, studies in 293 cells expressing TIG3 linked to an inducible promoter demonstrated decreased proliferation with increased TIG3 levels. These studies suggest that TIG3 may be a growth regulator that mediates some of the growth suppressive effects of retinoids.

Identification of semaphorin E as a non-MDR drug resistance gene of human cancers

To improve cancer chemotherapy, a better understanding of the molecular mechanisms of drug resistance is essential. To identify the molecules responsible for drug resistance that is unrelated to MDR1 or MRP gene products, a eukaryotic expression cDNA library of cis-diamminedichloroplatinum(II) (CDDP)-resistant ovarian cancer TYKnuR cells was introduced into Cos-7 cells. After repeated CDDP selection, cDNA homologous to murine semaphorin E was isolated from surviving cells. Human semaphorin E (H-sema E) was overexpressed in CDDP-resistant cell lines and was readily induced not only by diverse chemotherapeutic drugs but also by x-ray and UV irradiation. Transfection of H-sema E conferred a drug-resistant phenotype to CDDP-sensitive cells. In addition, the aberrant expression of H-sema E protein was detected immunohistochemically in 14 of 42 (33.3%) recurrent squamous cell carcinomas removed at autopsy after extensive radiochemotherapy. Recently, another member of the semaphorin family, CD100, was shown to significantly improve the viability of B lymphocytes. These results suggest the involvement of semaphorins in diverse cell survival mechanisms.

Identification of a novel vertebrate circadian clock-regulated gene encoding the protein nocturnin

Photoreceptors of the Xenopus laevis retina are the site of a circadian clock. As part of a differential display screen for rhythmic gene products in this system, we have identified a photoreceptor-specific mRNA expressed in peak abundance at night. cDNA cloning revealed an open reading frame encoding a putative 388 amino acid protein that we have named “nocturnin” (for night-factor). This protein has strong sequence similarity to the C-terminal domain of the yeast transcription factor, CCR4, as well as a leucine zipper-like dimerization motif. Nocturnin mRNA levels exhibit a high amplitude circadian rhythm and nuclear run-on analysis indicates that it is controlled by the retinal circadian clock at the level of transcription. Our observations suggest that nocturnin may function through protein–protein interaction either as a component of the circadian clock or as a downstream effector of clock function.

Identification of genes induced by factor deprivation in hematopoietic cells undergoing apoptosis using gene-trap mutagenesis and site-specific recombination

A strategy employing gene-trap mutagenesis and site-specific recombination (Cre/loxP) has been developed to isolate genes that are transcriptionally activated during programmed cell death. Interleukin-3 (IL-3)-dependent hematopoietic precursor cells (FDCP1) expressing a reporter plasmid that codes for herpes simplex virus–thymidine kinase, neomycin phosphotransferase, and murine IL-3 were transduced with a retroviral gene-trap vector carrying coding sequences for Cre-recombinase (Cre) in the U3 region. Activation of Cre expression from integrations into active genes resulted in a permanent switching between the selectable marker genes that converted the FDCP1 cells to factor independence. Selection for autonomous growth yielded recombinants in which Cre sequences in the U3 region were expressed from upstream cellular promoters. Because the expression of the marker genes is independent of the trapped cellular promoter, genes could be identified that were transiently induced by IL-3 withdrawal.