Decomposition of $\rm R\times C$ contingency table chi -square: Application to binned DNA fragment size data and population structure analysis

The purpose of this research is to develop a new statistical method to determine the minimum set of rows (R) in a R x C contingency table of discrete data that explains the dependence of observations. The statistical power of the method will be empirically determined by computer simulation to judge its efficiency over the presently existing methods. The method will be applied to data on DNA fragment length variation at six VNTR loci in over 72 populations from five major racial groups of human (total sample size is over 15,000 individuals; each sample having at least 50 individuals). DNA fragment lengths grouped in bins will form the basis of studying inter-population DNA variation within the racial groups are significant, will provide a rigorous re-binning procedure for forensic computation of DNA profile frequencies that takes into account intra-racial DNA variation among populations. ^

Metabolism and action of 9-$\beta$-D-arabinofuranosyl-2-fluoroadenine

9-$\beta$-D-arabinofuranosyl-2-fluoroadenine (F-ara-A) is an analogue of adenosine and 2$\sp\prime$-deoxyadenosine with potent antitumor activity both in vitro and in vivo. The mechanism of action of F-ara-A was evaluated both in whole cells and in experimental systems with purified enzymes. F-ara-A was converted to its 5$\sp\prime$-triphosphate F-ara-ATP in cells and then incorporated into DNA in a self-limiting manner. About 98% of the incorporated F-ara-AMP residues were located at the 3$\sp\prime$-termini of DNA strands, suggesting a chain termination property of this compound. DNA synthesis in CEM cells was inhibited by F-ara-A treatment with an IC$\sb{50}$ value of 1 $\mu$M. Cells were not able to restore the normal level of DNA synthesis even after being cultured in drug-free medium for 40 h. A DNA primer extension assay with M13mp18(+) single-stranded DNA template using purified human DNA polymerases $\alpha$ and further revealed that F-ara-ATP competed with dATP for incorporation into the A sites of the elongating DNA strands. The incorporation of F-ara-AMP into DNA resulted in a termination of DNA synthesis at the incorporated A sites. Pol $\alpha$ and $\delta$ were not able to efficiently extend the DNA primer with F-ara-AMP at its 3$\sp\prime$-end. Furthermore, the presence of F-ara-AMP at the 3$\sp\prime$-end of an oligodeoxyribonucleotide impaired its ligation with an adjacent DNA fragment by human and T4 ligases. Human DNA polymerase $\alpha$ incorporated more F-ara-AMP into DNA than polymerase $\delta$ and was more sensitive to the inhibition by F-ara-ATP, suggesting that polymerase $\alpha$ may be a preferred target for this analogue. On the other hand, DNA-dependent nucleotide turnover experiments and sequencing gel analysis demonstrated that DNA polymerase $\delta$ was able to remove the incorporated F-ara-AMP residue from the 3$\sp\prime$-end of the DNA strand with its 3$\sp\prime$-5$\sp\prime$ exonuclease activity in vitro, subsequently permitting further elongation of the DNA strand.^ The incorporation of F-ara-AMP into DNA was linearly correlated both with the inhibition of DNA synthesis and with the loss of clonogenicity. Termination of DNA synthesis and deletion of genetic material resulted from F-ara-AMP incorporation may be the mechanism responsible for cytotoxicity of F-ara-A. (Abstract shortened with permission of author.) ^

Characterization of the human tissue transglutaminase gene promoter

Transglutaminases are a family of calcium-dependent enzymes, that catalyze the covalent cross-linking of proteins by forming $\varepsilon(\gamma$-glutamyl)lysine isopeptide bonds. In order to investigate the molecular mechanisms regulating the expression of the tissue transglutaminase gene and to determine its biological functions, the goal of this research has been to clone and characterize the human tissue transglutaminase promoter. Thirteen clones of the tissue transglutaminase gene were obtained from the screening of a human placental genomic DNA library. A 1.74 Kb fragment derived from DNA located immediately upstream of the translation start site was subcloned and sequenced. Sequence analysis of this DNA fragment revealed that it contains a TATA box (TATAA), a CAAT box (GGACAAT), and a series of potential transcription factor binding sites and hormone response elements. Four regions of significant homology, a GC-rich region, a TG-rich region, an AG-rich region, and HR1, were identified by aligning 1.8 Kb of DNA flanking the human, mouse, and guinea pig tissue transglutaminase genes.^ To measure promoter activity, we subcloned the 1.74 Kb fragment of the tissue transglutaminase gene into a luciferase reporter vector to generate transglutaminase promoter/luciferase reporter constructs. Transfection experiments showed that this DNA segment includes a functional promoter with high constitutive activity. Deletion analysis revealed that the SP1 sites or corresponding sequences contribute to this activity. We investigated the role of DNA methylation in regulating the activity of the promoter and found that in vitro methylation of tissue transglutaminase promoter/luciferase reporter constructs suppressed their basal activity. Methylation of the promoter is inversely correlated with the expression of the tissue transglutaminase gene in vivo. These results suggest that DNA methylation may be one of the mechanisms regulating the expression of the gene. The tumor suppressor gene product p53 was also shown to inhibit the activity of the promoter, suggesting that induction of the tissue transglutaminase gene is not involved in the p53-dependent programmed cell death pathway. Although retinoids regulate the expression of the tissue transglutaminase gene in vivo, retinoid-inducible activity can not be identified in 3.7 Kb of DNA 5$\sp\prime$ to the tissue transglutaminase gene.^ The structure of the 5$\sp\prime$ end of the tissue transglutaminase gene was mapped. Alignment analysis of the human tissue transglutaminase gene with other human transglutaminases showed that tissue transglutaminase is the simplest member of transglutaminase superfamily. Transglutaminase genes show a conserved core of exons and introns but diverse N-terminuses and promoters. These observations suggest that key regulatory sequences and promoter elements have been appended upstream of the core transglutaminase gene to generate the diversity of regulated expression and regulated activity characteristic of the transglutaminase gene family. ^

Gene transfer by viral vectors

To initiate our clinical trial for chemotherapy protection, I established the retroviral vector system for human MDR1 cDNA gene transfer. The human MDR1 cDNA continued to be expressed in the transduced bone marrow cells after four cohorts of serial transplants, 17 months after the initial transduction and transplant. In addition, we used this retroviral vector pVMDR1 to transduce human bone marrow and peripheral blood CD34$\sp+$ cells on stromal monolayer in the presence of hematopoietic growth factors. These data suggest that the retroviral vector pVMDR1 could modify hematopoietic precursor cells with a capacity for long-term self renewal. Thus, it may be possible to use the MDR1 retroviruses to confer chemotherapeutic protection on human normal hematopoietic precursor cells of ovarian and breast cancer patients in whom high doses of MDR drugs may be required to control the diseases.^ Another promising vector system is recombinant adeno-associated virus (rAAV) vector. An impediment to use rAAV vectors is that production of rAAV vectors for clinical use is extremely cumbersome and labor intensive. First I set up the rAAV vector system in our laboratory and then, I focused on studies related to the production of rAAV vectors for clinical use. By using a self-inactivating retroviral vector carrying a selection marker under the control of the CMV immediate early promoter and an AAV genome with the deletion of both ITRs, I have developed either a transient or a stable method to produce rAAV vectors. These methods involve infection only and can generate high-titer rAAV vectors (up to 2 x 10$\sp5$ cfu/ml of CVL) with much less work.^ Although recombinant adenoviral vectors hardly infect early hematopoietic precursor cells lacking $\alpha\sb v\beta\sb5$ or $\alpha\sb v\beta\sb3$ integrin on their surface, but efficiently infect other cells, we can use these properties of adenoviral vectors for bone marrow purging as well as for development of new viral vectors such as pseudotyped retroviral vectors and rAAV vectors. Replacement of self-inactivating retroviral vectors by recombinant adenoviral vectors will facilitate the above strategies for production of new viral vectors. In order to accomplish these goals, I developed a new method which is much more efficient than the current methods to construct adenoviral vectors. This method involves a cosmid vector system which is utilized to construct the full-length recombinant adenoviral vectors in vitro.^ First, I developed an efficient and flexible method for in vitro construction of the full-length recombinant adenoviral vectors in the cosmid vector system by use of a three-DNA fragment ligation. Then, this system was improved by use of a two-DNA fragment ligation. The cloning capacity of recombinant adenoviral vectors constructed by this method to develop recombinant adenoviral vectors depends on the efficiency of transfection only. No homologous recombination is required for development of infectious adenoviral vectors. Thus, the efficiency of generating the recombinant adenoviral vectors by the cosmid method reported here was much higher than that by the in vitro direct ligation method or the in vivo homologous recombination method reported before. This method of the in vitro construction of recombinant adenoviral vectors in the cosmid vector system may facilitate the development of adenoviral vector for human gene therapy. (Abstract shortened by UMI.) ^

PHENOTYPIC ASSOCIATION OF GENE AMPLIFICATION WITH MULTIPLE DRUG RESISTANCE IN VINCRISTINE RESISTANT CHINESE HAMSTER OVARY CELLS

Resistance of tumors to pharmacologic agents poses a significant problem in the treatment of human malignancies. This study overviews the scope of clinical resistance and focuses upon current research attempts toward investigation of the phenomenon of multidrug resistance (MDR).^ The objective of this investigation was to determine whether gene amplification had a role in the development of the MDR phenotype in Chinese hamster ovary cells (CHO) primarily selected for resistance to vincristine (VCR). A DNA fragment, previously shown to be amplified in two independently derived Chinese hamster cell lines exhibiting the MDR phenotype, was also amplified in VCR hamster lines. Sequences flanking this fragment were shown to contain coding information for a 4.3 kb transcript overproduced in VCR cells. These sequences were not enriched in double minute DNA preparations isolated from VCR cells. There was an approximately forty-fold increase in both the level of gene amplification and transcript overproduction in the VCR cell lines, independent of the level of primary resistance. This DNA amplification and overproduction of the 4.3 kb transcript was also demonstrated in CHO cells independently selected for resistance to Adriamycin and vinblastine.^ All the DNA sequences of two hamster cDNA clones containing 785 and 932 base pair inserts showed direct homology to the published mouse mdr sequences (about 90%). This sequence conservation held for only portions of the gene when the human mdr1 sequences were compared with those from either the mouse or hamster.^ Somatic cell hybrids, constructed between VCR CHO cells and sensitive murine cells, were used to determine whether there was a functional relationship between the chromosome bearing the amplified sequences and the MDR phenotype. Concordant segregation between vincristine resistance, the MDR phenotype, the presence of MDR-associated amplified sequences, overexpression of the mRNA encoded by these sequences, overexpression of the mRNA encoded by these sequences, and CHO chromosome Z1 was consistent with the hypothesis that there is an amplified gene on chromosome Z1 of the VCR CHO cells which is responsible for MDR in these cells. ^

{\it cis\/} -acting elements and {\it trans\/} -acting factors that control serum amyloid A gene expression during inflammation

To understand how the serum amyloid A (SAA) genes are regulated, the cis-acting elements and trans-acting factors involved in the regulation of mouse SAA3 and rat SAA1 genes expression during inflammation were analyzed.^ To identify DNA sequences involved in the liver-specific expression of the mouse SAA3 gene, the 5$\sp\prime$ flanking region of this gene was analyzed by transient transfection studies. Results suggest that C/EBP, a liver-enriched transcription factor, plays an important role for the enhanced expression of the mouse SAA3 gene in hepatocytes.^ Transfection studies of the regulation of the expression of rat SAA1 gene indicated that a 322 bp fragment ($-$304 to +18) of the gene contains sufficient information for cytokine-induced expression of the reporter gene in a liver cell-specific manner. Further functional analysis of the 5$\sp\prime$ flanking region of the rat SAA1 gene demonstrated that a 65 bp DNA fragment ($-$138/$-$73) can confer cytokine-inducibility onto a heterologous promoter both in liver and nonliver cells. DNase I footprint and gel retardation assays identified five putative cis-regulatory elements within the 5$\sp\prime$ flanking region of the gene: one inducible element, a NF$\kappa$B binding site and four constitutive elements. Two constitutive elements, footprint regions I and III, were identified as C/EBP binding sites with region III having over a 10-fold higher affinity for C/EBP binding than region I. Functional analysis of the cis-elements indicated that C/EBP(I) and C/EBP(III) confer liver cell-specific activation onto a heterologous promoter, while sequences corresponding to the NF$\kappa$B element and C/EBP(I) impart cytokine responsiveness onto the heterologous promoter. These results suggest that C/EBP(I) possesses two functions: liver-specific activation and cytokine responsiveness. The identification of two cytokine responsive elements (NF$\kappa$B and C/EBP(I)), and two liver-specific elements (C/EBP(I) and C/EBP(III)) implies that multiple cis-acting elements are involved in the regulation of the expression of the rat SAA1 gene. The tissue-specific and cytokine-induced expression of rat SAA1 gene is likely the result of the interactions of these cis-acting elements with their cognate trans-acting factors as well as the interplay between the different cis-acting elements and their binding factors. (Abstract shortened with permission of author.) ^

Repression of {\it neu} oncogene by E1A: Mechanisms and biological functions

The potential effects of the E1A gene products on the promoter activities of neu were investigated. Transcription of the neu oncogene was found to be strongly repressed by the E1A gene products and this requires that conserved region 2 of the E1A proteins. The target for E1A repression was localized within a 140 base pair (bp) DNA fragment in the upstream region of the neu promoter. To further study if this transcriptional repression of neu by E1A can inhibit the transforming ability of the neu transformed cells, the E1A gene was introduced into the neu oncogene transformed B104-1-1 cells and developed B-E1A cell lines that express E1A proteins. These B-E1A stable transfectants have reduced transforming activity compared to the parental B104-1-1 cell line and we conclude that E1A can suppress the transformed phenotypes of the neu oncogene transformed cells via transcriptional repression of neu.^ To study the effects of E1A on metastasis, we first introduced the mutation-activated rat neu oncogene into 3T3 cells and showed that both the neu oncogene transformed NIH3T3 cells and Swiss Webster 3T3 cells exhibited metastatic properties in vitro and in vivo, while their parental 3T3 cells did not. Additionally, the neu-specific monoclonal antibody 7.16.4, which can down regulate neu-encoded p185 protein, effectively reduced the metastatic properties induced by neu. To investigate if E1A can reduce the metastatic potential of neu-transformed cells, we also compared the metastatic properties of B-E1A cell lines and B104-1-1 cell. B-E1A cell lines showed reduced invasiveness and lung colonization than the parental neu transformed B104-1-1 cells. We conclude that E1A gene products also have inhibitory effect on the metastatic phenotypes of the neu oncogene transformed cells.^ The product of human retinoblastoma (RB) susceptibility gene has been shown to complex with E1A gene products and is speculated to regulate gene expression. We therefore investigated in E1A-RB interaction might be involved in the regulation of neu oncogene expression. We found that the RB gene product can decrease the E1A-mediated repression of neu oncogene and the E1A binding region of the RB protein is required for the derepression function. ^

Functional analysis of the sea urchin {\it orthodenticle\/} -related gene, {\it SPOTX\/}, in aboral ectoderm-specific gene activation and aboral ectoderm cell fate specification

An important question in developmental biology is how embryonic cell types are derived from a fertilized egg. To address this question, this thesis investigates the mechanisms by which the aboral ectoderm-specific Spec2a gene is spatially and temporally regulated during sea urchin embryogenesis. The Spec2a gene of the sea urchin Strongylocentratus purpuratus has served as a valuable maker to understand the basis of lineage-specific gene activation and the role of transcription factors in cell fate specification. The hypothesis is that transcription factors responsible for cell type-specific gene activation are key components in the initial cell specification step. The Spec2a gene, which encodes a small cytosolic calcium-binding protein, is expressed exclusively in aboral ectoderm cell lineages. The 1516-bp control region of the Spec2a gene contains a 188-bp enhancer element required for temporal activation and aboral ectoderm/mesenchyme cell expression, while an unidentified element upstream of the enhancer represses expression in mesenchyme cells. Using an enhancer activation assay, combined with site-directed mutagenesis, I showed that three TAATCC/T sites within the enhancer are responsible for enhancer activity. Mutagenizing these sites and a fourth one just upstream abolished all activity from the Spec2a control region. A 77-bp DNA fragment from the Spec2a enhancer containing two of the TAATCC/T sites is sufficient for aboral ectoderm/mesenchyme cell expression. A cDNA encoding SpOtx, an orthodenticle-related protein, was cloned from S. purpuratus and shown to bind with high affinity to the TAATCC/T sequences within the Spec2a control region. SpOtx transcripts were found initially in all cells of the cleaving embryo, but they gradually became restricted to oral ectoderm and endoderm cells, suggesting that SpOtx might play a role in the initial temporal activation of the Spec2a gene and most likely has additional functions in the developing embryo. To reveal the broader biological functions of SpOtx, I injected SpOtx mRNA into living sea urchin eggs to determine what effects overexpressing the SpOtx protein might have on embryo development. SpOtx mRNA-injected embryos displayed dramatic alterations in development. Instead of developing into pluteus larvae with 15 different cell types, uniform epithelia balls were formed. These balls consisted of a thin layer of squamous cells with short cilia highly reminiscent of aboral ectoderm. Immunohistochemical staining and RT-PCR demonstrated that the SpOtx-injected embryoids expressed aboral ectoderm markers uniformly, but showed very weak or no expression of markers for non-aboral ectoderm cell types. These data strongly suggested that overexpression of SpOtx redirected the normal fate of non-aboral ectoderm cells to that of aboral ectoderm. These results show that SpOtx is involved in aboral ectoderm differentiation by activating aboral ectoderm-specific genes and that modulating its expression can lead to changes in cell fate. ^

Axl -Gas6 signaling counteracts adenovirus type 5 E1A-mediated cell growth suppression and proapoptotic activity

The adenovirus type 5 E1A (abbreviated E1A) has previously been known as an immortalization oncogene because E1A is required for transforming oncogenes, such as ras and E1B, to transform cells in primary cultures. However, E1A has also been shown to downregulate the overexpression of the Her-2/neu oncogene, resulting in suppression of transformation and tumorigenesis induced by that oncogene. In addition, E1A is able to promote apoptosis induced by anticancer drugs, irradiation, and serum deprivation. Many tyrosine kinases, such as the EGF receptor, Her-2/Neu, Src, and Axl are known to play a role in oncogenic signals in transformed cells. To study the mechanism underlying the E1A-mediated tumor-suppressing function, we exploited a modified tyrosine kinase profile assay (Proc. Natl. Acad. Sci, 93, 5958–5962, 1996) to identify potential tyrosine kinases regulated by E1A. RT-PCR products were synthesized with two degenerate primers derived from the conserved motifs of various tyrosine kinases. A tyrosine kinase downregulated by E1A was identified as Axl by analyzing the Alu I-digested RT-PCR products. We isolated the DNA fragment of interest, and found that E1A negatively regulated the expression of the transforming receptor tyrosine kinase Axl at the transcriptional level. To study whether downregulation of the Axl receptor is involved in E1A-mediated growth suppression, we transfected axl cDNA into E1A-expressing cells (ip1-E1A) to establish cells that overexpressed Axl (ip1-E1A-Axl). The Axl ligand Gas6 triggered a greater mitogenic effect in these ip1-E1A-Axl cells than in the control cells ip1-E1A and protected the Axl-expressing cells from serum deprivation-induced apoptosis. Further study showed that Akt is required for Axl-Gas6 signaling to prevent ip1-E1A-Axl cells from serum deprivation-induced apoptosis. These results indicate that downregulation of the Axl receptor by E1A is involved in E1A-mediated growth suppression and E1A-induced apoptosis, and thereby contributes to E1A's anti-tumor activities. ^

MATHEMATICAL STUDIES ON THE EVOLUTIONARY CHANGE OF DNA SEQUENCES

With the aim of understanding the mechanism of molecular evolution, mathematical problems on the evolutionary change of DNA sequences are studied. The problems studied and the results obtained are as follows: (1) Estimation of evolutionary distance between nucleotide sequences. Studying the pattern of nucleotide substitution for the case of unequal substitution rates, a new mathematical formula for estimating the average number of nucleotide substitutions per site between two homologous DNA sequences is developed. It is shown that this formula has a wider applicability than currently available formulae. A statistical method for estimating the number of nucleotide changes due to deletion and insertion is also developed. (2) Biases of the estimates of nucleotide substitutions obtained by the restriction enzyme method. The deviation of the estimate of nucleotide substitutions obtained by the restriction enzyme method from the true value is investigated theoretically. It is shown that the amount of the deviation depends on the nucleotides in the recognition sequence of the restriction enzyme used, unequal rates of substitution among different nucleotides, and nucleotide frequences, but the primary factor is the unequal rates of nucleotide substitution. When many different kinds of enzymes are used, however, the amount of average deviation is generally small. (3) Distribution of restriction fragment lengths. To see the effect of undetectable restriction fragments and fragment differences on the estimate of nucleotide differences, the theoretical distribution of fragment lengths is studied. This distribution depends on the type of restriction enzymes used as well as on the relative frequencies of four nucleotides. It is shown that undetectability of small fragments or fragment differences gives a serious underestimate of nucleotide substitutions when the length-difference method of estimation is used, but the extent of underestimation is small when the site-difference method is used. (4) Evolutionary relationships of DNA sequences in finite populations. A mathematical theory on the expected evolutionary relationships among DNA sequences (nucleons) randomly chosen from the same or different populations is developed under the assumption that the evolutionary change of nucleons is determined solely by mutation and random genetic drift. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author). UMI ^

The effects of DNA cytosine methylation on H -ras promoter structure and function

The effect of DNA cytosine methylation on H-ras promoter activity was assessed using a transient expression system employing the plasmid H-rasCAT (NaeI H-ras promoter linked to the chloramphenicol acetyltransferase (CAT) gene). This 551 bp promoter is 80% GC rich, enriched with 168 CpG dinucleotides, and contains six functional GC box elements which represent major DNA methylation target sites. Prokaryotic methyltransferases HhaI (CGm$\sp5$CG) and HpaII (Cm$\sp5$CGG) alone or in combination with a human placental methyltransferase (HP MTase) were used to introduce methyl groups at different CpG sites within the promoter. To test for functional promoter activity, the methylated plasmids were introduced into CV-1 cells and CAT activity assessed 48 h post-transfection. Methylation at specific HhaI and HpaII sites reduced CAT expression by 70%, whereas more extensive methylation at generalized CpG sites with HP MTase inactivated the promoter $>$95%. The inhibition of H-ras promoter activity was not attributable to methylation-induced differences in DNA uptake or stability in the cell, topological form of the plasmid, or methylation effects in nonpromoter regions. We also observed that DNA cytosine methylation of a 360 bp promoter fragment by HP MTase induced a local change in DNA conformation. Using three independent methodologies (nitrocellulose filter binding assays, gel mobility shifts, and Southwestern blots), we determined that this change in promoter conformation affected the interaction of nuclear proteins with cis-regulatory sequences residing in the promoter region. The results provide evidence to suggest that DNA methylation may regulate gene expression by inducing changes in local promoter conformation which in turn alters the interactions between DNA and protein factors required for transcription. The results provide supportive evidence for the hypothesis of Cedar and Riggs, who postulated that DNA methylation may regulate gene expression by altering the binding affinities of proteins for DNA. ^