Differential gene expression in multistep tumorigenesis using an {\it in vitro\/} human cell model

Using a human terato-carcinoma cell line, PA-1, the functional role of the oncogenes and tumor suppressor gene involved in the multistep process of carcinogenesis have been analyzed. The expression of AP-2 was strongly correlated with the susceptibility to ras transformation. The differential responsiveness to growth factors between stage 1 ras resistant cells and stage 2 ras susceptible cells was observed, indicating that the ability of stage 2 cells to respond to the mutated ras oncogenes in transformation correlated with the ability to be stimulated by certain growth factors. Using differential screening of cDNA libraries, a number of differentially expressed cDNA clones was isolated. One of those, clone 12, is overexpressed in ras transformed stage 3 cells. The amino acid sequence of clone 12 is almost identical to a mouse LLrep3 gene that was growth-regulated, and 78% similar to a yeast ribosomal protein S4. These results suggest that the S4 gene may be involved in regulation of growth. Clone 9 is expressed in stage 1 ras resistant cells (3.5-kb and 3.0-kb transcripts) but the expression of this clone in stage 2 ras susceptible cells and stage 3 ras-transformed cells is greatly diminished. The expression of this cDNA clone was increased to at least five fold in ras resistant cells and nontumorigenic hybrids treated with retinoic acid but not increased in retinoic acid treated ras susceptible cells, ras transformed cells and the tumorigenic segregants. Partial sequence of this clone showed no homology to the sequences in Genbank. These findings suggest that clone 9 could be a suppressor gene or the genes that are involved in the biochemical pathway of tumor suppression or neurogenic differentiation. The apparent pleiotropic effect of the loss of this suppressor gene function support Harris' proposal that tumor suppressor genes regulate differentiation. The tumor suppressor gene may act as negative regulator of tumor growth by controlling gene expression in differentiation. ^

Analysis of differential gene expression in nontumorigenic glioblastoma multiforme

Loss of chromosome 10 represents the most common cytogenetic abnormality in high grade gliomas (glioblastoma multiforme). To identify genes involved in the malignant progression of human gliomas, a subtractive hybridization was performed between a tumorigenic glioblastoma cell line (LG11) and a nontumorgenic hybrid cell (LG11.3) containing an introduced chromosome 10. LG11 mRNA was subtracted from LG11.3 cDNA to produce cDNA probes enriched for sequences whose expression differs quantitatively from the parental tumorigenic cells. Both known and novel sequences were identified as a result of the subtraction. Northern blot analysis was then used to confirm differential expression of several subtracted clones. One novel clone, clone 17, identified a 2.6 kb message that showed a consistent two to four fold increase in expression in the LG11.3 nontumorigenic cells. Clone 17 (340 bp) was used successfully to screen for a near full-length version, RIG (regulated in glioma), which was 2,569 bp in size. The RIG cDNA sequence showed homology to clone 17 and to an anonymous EST (IB666), but to no previously identified genes. This screening effort also identified several independent clones representing novel sequences, most of which failed to show increased expression in the nontumorigenic GBM cells. Tissue distribution studies of RIG indicated highest levels of expression in human brain with appreciably lower levels in heart and lung. In vitro transcription and translation experiments demonstrated the ability of RIG to direct the synthesis of a 13 kD protein product. However, open reading frame analysis revealed no identify with previously described motifs or any known proteins. Using a combination of somatic cell hybrid panels and in situ hybridization, the RIG gene was mapped to chromosome 11p14-11p15. Further study of RIG and related gene products may provide insight into the negative regulation of glial oncogenesis. ^

Comparison of multiple comparison methods for identifying differential gene expression in simulated and real papillary thyroid cancer microarray data

The difficulty of detecting differential gene expression in microarray data has existed for many years. Several correction procedures try to avoid the family-wise error rate in multiple comparison process, including the Bonferroni and Sidak single-step p-value adjustments, Holm's step-down correction method, and Benjamini and Hochberg's false discovery rate (FDR) correction procedure. Each multiple comparison technique has its advantages and weaknesses. We studied each multiple comparison method through numerical studies (simulations) and applied the methods to the real exploratory DNA microarray data, which detect of molecular signatures in papillary thyroid cancer (PTC) patients. According to our results of simulation studies, Benjamini and Hochberg step-up FDR controlling procedure is the best process among these multiple comparison methods and we discovered 1277 potential biomarkers among 54675 probe sets after applying the Benjamini and Hochberg's method to PTC microarray data.^

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.

THE ROLE OF MAP KINASES ON THE FUNCTIONAL HETEROGENEITY OF HUMAN CD8+ T CELL MATURATION SUBSETS

While prior studies have focused on naïve (CD45RA+CD27+) and early stage memory (CD45RA-CD27+) CD8+ T cells, late memory CD8+ T cells (CD45RA+CD27) have received less interest because this subset of T cells is generally recognized as effectors, which produce IFNγ (but no IL-2) and perforin. However, multiple studies suggest that late memory CD8+ T cells may provide inadequate protection in infectious diseases and cancer models. To better understand the unique function of late memory CD8+ T cells, I optimized multi-color flow cytometry techniques to assess the cytokine production of each human CD8+ T cell maturation subset. I demonstrated that late memory CD8+ T cells are the predominant producer of CC chemokines (e.g. MIP-1β), but rarely produce IL-2; therefore they do not co-produce IL-2/IFNγ (polyfunctionality), which has been shown to be critical for protective immunity against chronic viral infection. These data suggest that late memory CD8+ T cells are not just cytotoxic effectors, but may have unique functional properties. Determining the molecular signature of each CD8+ T cell maturation subset will help characterize the role of late memory CD8+ T cells. Prior studies suggest that ERK1 and ERK2 play a role in cytokine production including IL-2 in T cells. Therefore, I tested whether differential expression of ERK1 and ERK2 in CD8+ T cell maturation subsets contributes to their functional signature by a novel flow cytometry technique. I found that the expression of total ERK1, but not ERK2, is significantly diminished in late memory CD8+ T cells and that ERK1 expression is strongly associated with IL-2 production and CD28 expression. I also found that IL-2 production is increased in late memory CD8+ T cells by over-expressing ERK1. Collectively, these data suggest that ERK1 is required for IL-2 production in human CD8+ T cells. In summary, this dissertation demonstrated that ERK1 is down-regulated in human late memory CD8+ T cells, leading to decreased production of IL-2. The data in this dissertation also suggested that the functional heterogeneity in human CD8+ T cell maturation subsets results from their differential ERK1 expression.

Transcription factor SEF: Purification and functional characterization

Cytokine-induced transcription of the serum amyloid A3 (SAA3) gene promoter requires a transcriptional enhancer that contains three functional elements: two C/EBP-binding sites and a third site that interacts with a constitutively expressed transcription factor, SAA3 enhancer factor (SEF). Deletion or site-specific mutations in the SEF-binding site drastically reduced SAA3 promoter activity, strongly suggesting that SEF is important in SAA3 promoter function. To further elucidate its role in the regulation of the SAA3 gene, we purified SEF from HeLa cell nuclear extracts to near homogeneity by using conventional liquid chromatography and DNA-affinity chromatography. Ultraviolet cross-linking and Southwestern experiments indicated that SEF consisted of a single polypeptide with an apparent molecular mass of 65 kDa. Protein sequencing, oligonucleotide competition and antibody supershift experiments identified SEF as transcription factor LBP-1c/CP2/LSF. Cotransfection of SEF expression plasmid with SAA3-luciferase reporter resulted in 3- to 5-fold activation of SAA3 promoter. Interestingly, when SEF-transfected cells were treated with either conditioned medium (CM) or interleukin (IL) 1, the SAA3 promoter was synergistically activated in a dose-dependent manner. Furthermore, when SEF-binding site was mutated, the response of SAA3 promoter to IL-1 or CM stimulation was abolished or drastically decreased, suggesting that SEF may functionally cooperate with an IL-1-inducible transcription factor. Indeed, our functional studies showed that NFκB is a key transcription factor that mediates the IL-1-induced expression of SAA3 gene, and that SEF can synergize with NFκBp65 to activate SAA3 promoter. By coimmunoprecipitation experiments, we found that SEF could specifically interact with NFκBp65, and that the association of these two factors was enhanced upon IL-1 and CM stimulation. This suggests that the molecular basis for the functional synergy between SEF and NFκB may be due to the ability of SEF to physically interact with NPκB. In addition to its interaction with SEF, NFκB-dependent activation also requires the weak κB site in the C element and its interaction with C/EBP. Besides its role in regulating SAA3 gene expression, we provide evidence that SEF could also bind in a sequence-specific manner to the promoters of α2-macroglobulin, Aα fibrinogen, and 6–16 genes and to an intronic enhancer of the human Wilm's tumor 1 gene, suggesting a functional role in the regulation of these genes. By coimmunoprecipitation experiments, we determined that SEF could specifically associate with both Stat3 and Stat2 upon cytokine stimulation. To examine the functional roles of such interactions, we evaluated the effects of SEF on the transcriptional regulation of two reporter genes: Aα fibrinogen and 6–16, which are IL-6- and interferon-α-responsive, respectively. Our results showed that cotransfection of SEF expression plasmid can activate the expression of Aα fibrinogen gene and 6–16 gene. Moreover, SEF can dramatically enhance the interferon-α-induced expression of 6–16 gene and IL-6-induced expression of Aα fibrinogen gene, suggesting that SEF may functionally cooperate with ISGF3 and Stat3 to mediate interferon-α and IL-6 signaling. ^ Our findings that SEF can interact with multiple cytokine-inducible transcription factors to mediate the expression of target genes open a new avenue of investigation of cooperative transcriptional regulation of gene expression, and should further our understanding of differential gene expression in response to a specific stimulus. In summary, our data provide evidence that SEF can mediate the signaling of different cytokines by interacting with various cytokine-inducible transcription factors. ^

Targeted Deletion Of Functionally Validated Enhancers Defines Their Role in Mouse Limb Development

Transcriptional enhancers are genomic DNA sequences that contain clustered transcription factor (TF) binding sites. When combinations of TFs bind to enhancer sequences they act together with basal transcriptional machinery to regulate the timing, location and quantity of gene transcription. Elucidating the genetic mechanisms responsible for differential gene expression, including the role of enhancers, during embryological and postnatal development is essential to an understanding of evolutionary processes and disease etiology. Numerous methods are in use to identify and characterize enhancers. Several high-throughput methods generate large datasets of enhancer sequences with putative roles in embryonic development. However, few enhancers have been deleted from the genome to determine their roles in the development of specific structures, such as the limb. Manipulation of enhancers at their endogenous loci, such as the deletion of such elements, leads to a better understanding of the regulatory interactions, rules and complexities that contribute to faithful and variant gene transcription – the molecular genetic substrate of evolution and disease. To understand the endogenous roles of two distinct enhancers known to be active in the mouse embryo limb bud we deleted them from the mouse genome. I hypothesized that deletion of these enhancers would lead to aberrant limb development. The enhancers were selected because of their association with p300, a protein associated with active transcription, and because the human enhancer sequences drive distinct lacZ expression patterns in limb buds of embryonic day (E) 11.5 transgenic mice. To confirm that the orthologous mouse enhancers, mouse 280 and 1442 (M280 and M1442, respectively), regulate expression in the developing limb we generated stable transgenic lines, and examined lacZ expression. In M280-lacZ mice, expression was detected in E11.5 fore- and hindlimbs in a region that corresponds to digits II-IV. M1442-lacZ mice exhibited lacZ expression in posterior and anterior margins of the fore- and hindlimbs that overlapped with digits I and V and several wrist bones. We generated mice lacking the M280 and M1442 enhancers by gene targeting. Intercrosses between M280 -/+ and M1442 -/+, respectively, generated M280 and M1442 null mice, which are born at expected Mendelian ratios and manifest no gross limb malformations. Quantitative real-time PCR of mutant E11.5 limb buds indicated that significant changes in transcriptional output of enhancer-proximal genes accompanied the deletion of both M280 and M1442. In neonatal null mice we observed that all limb bones are present in their expected positions, an observation also confirmed by histology of E18.5 distal limbs. Fine-scale measurement of E18.5 digit bone lengths found no differences between mutant and control embryos. Furthermore, when the developmental progression of cartilaginous elements was analyzed in M280 and M1442 embryos from E13.5-E15.5, transient development defects were not detected. These results demonstrate that M280 and M1442 are not required for mouse limb development. Though M280 is not required for embryonic limb development it is required for the development and/or maintenance of body size – adult M280 mice are significantly smaller than control littermates. These studies highlight the importance of experiments that manipulate enhancers in situ to understand their contribution to development.

Differential expression of urate oxidase in rat liver

An affinity-purified monospecific antibody was prepared to study the differential expression of the peroxisomal enzyme urate oxidase in rat liver during development and in various metabolic states. Monospecific antibody for urate oxidase was affinity purified from a pool of antibodies initially produced against a mixture of proteins from a Percoll density gradient fraction. Immunogold staining of samples of the gradient fraction and rat liver tissue with the affinity-purified antibody demonstrated labelling of peroxisomal core structures. Screening of liver homogenates from rats at different developmental stages using immunoblot analysis demonstrated low levels of urate oxidase prior to 20 days of age; at 20 days of age, urate oxidase levels are 2-fold greater than the 15-day old levels and approximate adult levels. Catalase expression during rat development mimicked the differential expression pattern of urate oxidase. The increase between days 15 and 20 was determined to be independent of the process of weaning. Administration of exogenous glucocorticoid hormone to 10-day old rats resulted in a precocious rise (2.5-fold) in urate oxidase levels, but adrenalectomy at 10 days of age did not cause decreased expression in the fourth week of life. In adult animals, exogenous glucocorticoid did not influence urate oxidase levels, but adrenalectomized rats had urate oxidase levels that were 40 percent of control expression 4 days post-surgery. Catalase expression was not influenced by glucocorticoid status in these studies. Glucocorticoid regulation of urate oxidase expression appears to be one part of a more complex mechanism controlling levels of the enzyme. Exogenous glucocorticoid administration influenced urate oxidase levels in an age-dependent manner; in addition, it is possible that the control mechanism for urate oxidase may include factors which can modulate expression in the absence of glucocorticoids. The effect of glucocorticoids on urate oxidase expression can not be extended to include all peroxisomal proteins, since catalase is unaffected. Glucocorticoids appear to participate in the complex regulation of urate oxidase expression; glucocorticoids influence urate oxidase specifically and do not modulate all peroxisomal proteins. ^

Transcriptional analysis of liver-specific expression and differential interleukin-6 response of rat kininogen genes

Analyses of rat T1 kininogen gene/chloramphenicol acetyltransferase (T1K/CAT) constructs revealed two regions important for tissue-specific and induced regulation of T1 kininogen.^ Although the T1 kininogen gene is inducible by inflammatory cytokines, a highly homologous K kininogen gene is minimally responsive. Moreover, the basal expression of a KK/CAT construct was 5- to 7-fold higher than that of the analogous T1K/CAT construct. To examine the molecular basis of this differential regulation, a series of promoter swapping experiments was carried out. Our transfection results showed that at least two regions in the K kininogen gene are important for its high basal expression: a distal 19-bp region (C box) constituted a binding site for CCAAT/enhancer binding protein (C/EBP) family proteins and a proximal 66-bp region contained two adjacent binding sites for hepatocyte nuclear factor-3 (HNF-3). The distal HNF-3 binding site from the K kininogen promoter demonstrated a stronger affinity than that from the T1 kininogen promoter. Since C/EBP and HNF-3 are highly enriched in the liver and known to enhance transcription of liver-specific genes, differential binding affinities of these factors accounted for the higher basal expression of the K kininogen gene.^ In contrast to the K kininogen C box, the T1 kininogen C box does not bind C/EBP presumably due to their two-nucleotide divergence. This sequence divergence, however, converts it to a consensus binding sequence for two IL-6-inducible transcription factors--IL-6 response element binding protein and acute-phase response factor. To functionally determine whether C box sequences are important for their differential acute-phase response, T1 and K kininogen C boxes were swapped and analyzed after transfection into Hep3B cells. Our results showed that the T1 kininogen C box is indeed one of the IL-6 response elements in T1 kininogen promoter. Furthermore, its function can be modulated by a 5$\sp\prime$-adjacent C/EBP-binding site (B box) whose mutation significantly reduced the overall induced activity. Moreover, this B box is the target site for binding and transactivation of another IL-6 inducible transcription factor C/EBP$\delta.$ Evolutionary divergence of a few critical nucleotides can either lead to subtle changes in the binding affinities of a given transcription factor or convert a binding sequence for a constitutive factor to a site recognized by an inducible factor. (Abstract shortened by UMI.) ^

The differential expression of the $\beta\sb2$-adrenergic receptor modifies agonist-dependent adenylylcyclase activation

There have been multiple reports which indicate that variations in $\beta$AR expression affect the V$\sb{\rm max}$ observed for the agonist-dependent activation of adenylylcyclase. This observation has been ignored by most researchers when V$\sb{\rm max}$ values obtained for wild type and mutant receptors are compared. Such an imprecise analysis may lead to erroneous conclusions concerning the ability of a receptor to activate adenylylcyclase. Equations were derived from the Cassel-Selinger model of GTPase activity and Tolkovsky and Levitzki's Collision Coupling model which predict that the EC$\sb{50}$ and V$\sb{\rm max}$ for the activation of adenylylcyclase are a function of receptor number. Experimental results for L cell clones in which either hamster or human $\beta$AR were transfected at varying levels showed that EC$\sb{50}$ decreases and V$\sb{\rm max}$ increases as receptor number increases. Comparison of these results with simulations obtained from the equations describing EC$\sb{50}$ and V$\sb{\rm max}$ showed a close correlation. This documents that the kinetic parameters of adenylylcyclase activation change with the level of receptor expression and relates this phenomenon to a theoretical framework concerning the mechanisms involved in $\beta$AR signal transduction.^ One of the terms used in the equations which expressed the EC$\sb{50}$ and V$\sb{\rm max}$ as a function of receptor number is coupling efficiency, defined as $\rm k\sb1/k\sb{-1}$. Calculation of $\rm k\sb1/k\sb{-1}$ can be accomplished for wild type receptors with the easily measured experimental values of agonist K$\sb{\rm d}$, EC$\sb{50}$ and receptor number. This was demonstrated for hamster $\beta$AR which yielded a coupling efficiency of 0.15 $\pm$ 0.003 and human $\beta$AR which yielded a coupling efficiency of 0.90 $\pm$ 0.031. $\rm k\sb1/k\sb{-1}$ replaces the traditional qualitative evaluation of the ability to activate adenylylcyclase, which utilizes V$\sb{\rm max}$ without correction for variation in receptor number, with a quantitative definition that more accurately describes the ability of $\beta$AR to couple to G$\sb{\rm s}$.^ The equations which express the EC$\sb{50}$ and V$\sb{\rm max}$ for adenylylcyclase activation as a function of receptor number and coupling efficiency were tested to determine whether they could accurately simulate the changes seen in these parameters during desensitization. Data from original desensitization experiments and data from the literature (24,25,52,54,83) were compared to simulated changes in EC$\sb{50}$ and V$\sb{\rm max}$. In a variety of systems the predictions of the equations were consistent with the changes observed in EC$\sb{50}$ and V$\sb{\rm max}$. In addition reductions in the calculated value of $\rm k\sb1/k\sb{-1}$ was shown to correlate well with $\beta$AR phosphorylation and to be minimally affected by sequestration and down-regulation. ^

A model of the roles of essential kinases in the induction and expression of late long-term potentiation.

The induction of late long-term potentiation (L-LTP) involves complex interactions among second-messenger cascades. To gain insights into these interactions, a mathematical model was developed for L-LTP induction in the CA1 region of the hippocampus. The differential equation-based model represents actions of protein kinase A (PKA), MAP kinase (MAPK), and CaM kinase II (CAMKII) in the vicinity of the synapse, and activation of transcription by CaM kinase IV (CAMKIV) and MAPK. L-LTP is represented by increases in a synaptic weight. Simulations suggest that steep, supralinear stimulus-response relationships between stimuli (e.g., elevations in [Ca(2+)]) and kinase activation are essential for translating brief stimuli into long-lasting gene activation and synaptic weight increases. Convergence of multiple kinase activities to induce L-LTP helps to generate a threshold whereby the amount of L-LTP varies steeply with the number of brief (tetanic) electrical stimuli. The model simulates tetanic, -burst, pairing-induced, and chemical L-LTP, as well as L-LTP due to synaptic tagging. The model also simulates inhibition of L-LTP by inhibition of MAPK, CAMKII, PKA, or CAMKIV. The model predicts results of experiments to delineate mechanisms underlying L-LTP induction and expression. For example, the cAMP antagonist RpcAMPs, which inhibits L-LTP induction, is predicted to inhibit ERK activation. The model also appears useful to clarify similarities and differences between hippocampal L-LTP and long-term synaptic strengthening in other systems.

EXPRESSION OF THE HUMAN PLACENTAL LACTOGEN GENES

The placenta is the site of synthesis of various peptide and steroid hormones related to pregnancy. Human placental lactogen (hPL) is the predominant peptide hormone secreted by term placenta and its synthesis is tissue-specific and coupled to placenta development. The objective of this work was to study the structure and expression of the hPL.^ Poly(A('+))RNA from human term placenta was translated in a mouse-derived cell-free system. A major band corresponding to pre-hPL and a minor band comigrating with mature hPL, represent (TURN)15% of the total radioactively labeled proteins. Analysis of the poly(A('+))RNA showed a prominent band at approximately 860 nucleotides. A corresponding band was observed in Northern blots of total RNA, hybridized with {('32)P}-labeled recombinant plasmid containing a portion of hPL cDNA. Similar analyses of nuclear RNA showed at least four additional bands at 990, 1200, 1460 and 1760 nucleotides, respectively, which are likely precursors of hPL mRNA. Poly(A('+))RNA was used to construct a cDNA library, of which approximately 5% of the clones were found to hybridize to hPL DNA sequences. Heteroduplexes constructed between a clone containing a 815 bp hPL cDNA insert and a hPL genomic DNA clone revealed four small intervening sequences which can account for the lengths observed in hnRNA molecules.^ Recombinant plasmid HCS-pBR322 containing a 550 bp insert of a cDNA transcript of human placental lactogen (hPL) mRNA was ('3)H-labeled an hybridized in situ to human chromosome preparations. These experiments allowed assignment of the hPL and growth hormone (hGH) genes, which have over 90% nucleotide homology in their coding sequences, to band q22-24 of chromosome 17. A gene copy number experiment showed that both genes are present in (TURN)3 copies per haploid genome.^ Experiments were designed to determine if all members of the hPL gene cluster, consisting of four non-allelic genes, are transcribed in term placenta. Advantage was taken of differences in restriction endonuclease sites in the coding portions of the different hPL genes, to distinguish the putative cDNAs of the transcriptionally active genes. Two genes were found to be represented in the cDNA library and their cDNA transcripts were isolated and characterized. Three independent methods showed that their corresponding mRNAs are about equally represented in the hPL mRNA population. The two cDNAs code for prehPL proteins which differ at a single amino acid position. However the secreted hPLs have identical amino acid sequences. A tetramer insertion duplication was found in a palindrome area of the 3' untranslated region of one of the hPL mRNAs. ^

INCREASED EXPRESSION OF ONE OF TWO ADENOSINE DEAMINASE ALLELES IN A HUMAN CHORIOCARCINOMA CELL LINE FOLLOWING SELECTION WITH ADENINE NUCLEOSIDES

The human choriocarcinoma cell line JEG-3 is heterozygous at the adenosine deaminase (ADA) gene locus. Both allelic genes are under strong but incomplete repression causing a very low level expression of the gene locus. Because cytotoxic adenosine analogues such as 9-(beta)-D arabinofuranosyladenine (ara-A) and 9-(beta)-D xylofuranosyladenine (xyl-A) can be specifically detoxified by the action of ADA, these analogues were used to select for JEG-3 derived cells which had increased ADA expression. When JEG-3 cells were subjected to a multi-step, successively increasing dosage of either ara-A or xyl-A, resistant cells with increased ADA expression were generated. This increased ADA expression in the resistant cells was unstable, so that when the selective pressure was removed, cellular ADA expression would decrease. Subclone analysis of xyl-A resistant cells revealed that compared to parental JEG-3 cells, individual resistant cells had either elevated ADA levels or decreased adenosine kinase (ADK) levels or both. This altered ADA and ADK expression in the resistant cells were found to be independent events. Because of high endogenous tissue conversion factor (TCF) expression in the JEG-3 cells, the allelic nature of the increased ADA expression in most of the resistant cells could not be determined. However, several resistant subcloned cells were found to have lost TCF expression. These TCF('-) cells expressed only the ADA*2 allelic gene product. Cell fusion experiments demonstrated that the ADA*1 allelic gene was intact and functional in the A3-1A7 cell line. Chromosomal analysis of the A3-1A7 cells showed that they had no double-minutes or homogeneously staining chromosomal regions, although a pair of new chromosomes were found in these cells. Segregation analysis of the hybrid cells indicated that an ADA*2 allelic gene was probably located on this new chromosome. The analysis of the A3-1A7 cell line suggested that the expression of only ADA 2 in these cells was the result of possibly a cis-deregulation of the ADA gene locus or more probably an amplification of the ADA*2 allelic gene. Two effective positive selection systems for ADA('+) cells were also developed and tested. These selection systems should eventually lead to the isolation of the ADA gene.^

Pharmacogenetics of the human glutathione S-transferase P1 gene in the metabolism and therapeutic efficacy of cis-diamminedichloroplatinum

The human GSTP1 gene has been shown, conclusively, to be polymorphic. The three main GSTP1 alleles, GSTP1*A, GSTP1*B, and GSTP1*C, encode proteins which differ in the 3-dimensional structure of their active sites and in their function in phase II metabolism of carcinogens, mutagens, and anticancer agents. Although, it is well established that GSTP1 is over expressed in many human tumors and that the levels of GSTP1 expression correlate directly with tumor resistance to chemotherapy and inversely with patient survival, the significance of the polymorphic GSTP1 gene locus on tumor response to chemotherapy remains unclear. The goal of this project was to define the role and significance of the polymorphic GSTP1 gene locus in GSTP1-based tumor drug resistance and as a determinant of patient response to chemotherapy. The hypothesis to be tested was that the polymorphic GSTP1 gene locus will confer to tumors a differential ability to metabolize cisplatin resulting in a GSTP1 genotype-based sensitivity to cisplatin. The study examined: (a) whether the different GSTP 1 alleles confer different levels of cellular protection against cisplatin-induced cytotoxicity, (b) whether the allelic GSTP1 proteins metabolize cisplatin with different efficiencies, and (c) whether the GSTP1 genotype is a determinant of tumor response to cisplatin therapy. The results demonstrate that the GSTP1 alleles differentially protect tumors against cisplatin-induced apoptosis and clonogenic cell kill in the rank order: GSTP1*C > GSTP1*B > GSTP1*A. The same rank order was observed for the kinetics of GSTP1-catalyzed cisplatin metabolism, both in cell-free and cellular systems, to the rate-limiting monoglutathionyl-platinum metabolite, which was characterized, for the first time, by mass spectral analysis. Finally, this study demonstrates that both GSTP1 genotype and the level of GSTP1 expression significantly contribute to tumor sensitivity to cisplatin treatment. Overall, the results of this project show that the polymorphic GSTP1 gene locus plays a significant role in tumor sensitivity to cisplatin treatment. Furthermore, these studies have contributed to the overall understanding of the significance of the polymorphic GSTP1 gene locus in tumor resistance to cancer chemotherapy and have provided the basis for further investigations into how this can be utilized to optimize and individualize cancer chemotherapy for cancer patients. ^

Identification of putative 1p36 tumor suppressor genes involved in oligodendroglioma development

Oligodendrogliomas are primary neoplasms of the central nervous system (CNS). One of the most common and characteristic chromosomal abnormalities observed in oligodendroglioma is allelic loss of 1p (Reifenberger et al., 1994; Bello et al., 1995). Since 1p loss has been reported for both well-differentiated and anaplastic oligodendroglioma, it is believed to occur early in tumor development (Bello et al., 1995). This allelic loss also has clinical significance, for oligodendroglioma patients with 1p loss generally respond significantly better to combination chemotherapy and have longer average survival than do oligodendroglioma patients without 1p loss (Cairncross et al., 1998). To date, no genes on 1p have been implicated as essential to the development or treatment response of oligodendroglioma. In order to localize and/or identify a gene involved in oligodendroglioma development, I tested 170 oligodendrogliomas for deletions of 1p and tested 26 tumors for differential expression of genes in the region of 1p36. Evidence obtained from these methods implicated two genes, SHREW1 and the gene encoding DNA fragmentation factor beta (DFFB). The function for the SHREW1 locus is currently not well known, but preliminary data suggests that it a novel member of adherens junctions. The DFFB gene is an enhancer for apoptosis. Thus, both SHREW1 and DFFB may be candidates for an oligodendroglioma tumor suppressor. Mutational analysis of both genes did not uncover any mutations. Future studies will evaluate other mechanisms that may be responsible for inactivation of these genes in oligodendrogliomas. ^