Molecular and biochemical analysis of the {\it Drosophila\/} segmentation gene {\it runt\/}

Genetic evidence has indicated that the segmentation gene runt plays a key role in regulating gene expression of the pair-rule genes hairy, even-skipped, and fushi tarazu. In contrast to other pair-rule genes, sequence data of the runt open reading frame did not reveal homologies to DNA-binding motifs of known transcriptional regulatory proteins. This thesis project examined several properties of the runt gene based on the sequence of the transcription unit, including the subcellular localization of the protein in vivo, its ability to bind DNA, and the functionality of a putative nucleotide binding domain.^ A runt-specific antibody was generated and used to demonstrate that runt is localized in the nucleus. Since the precise overlap of the pair-rule stripes is thought to be critical for the determination of cellular identity along the anterior-posterior axis, phasing of early runt expression in the blastoderm was examined with regard to the segmentation genes hairy, even-skipped, and fushi tarazu. runt was also expressed at later stages of embryogenesis, including expression in neuroblasts, and ganglion mother cells of the developing nervous system. Expression at this stage was required for the subsequent formation of specific neurons and runt was extensively expressed in the central and peripheral nervous systems.^ Several experiments were done to address the biochemical function of the runt protein. A direct interaction of runt with DNA was first examined. Although bacterial expressed runt was found to bind dsDNA-cellulose, subsequent experiments failed to detect sequence-specific interactions with DNA. Inter-species conservation of the putative nucleotide binding domain suggested that this region was functionally important, and runt protein bound a labeled ATP analog with high affinity in vitro. Finally, the effect of substitution of a critical residue of the nucleotide binding domain on runt activity was examined in vivo. Ectopic expression of the mutant protein indicated that this conserved substitution altered, but did not eliminate, runt activity as evaluated by segmentation phenotype and viability. ^

The cloning and expression of myogenin, a gene regulating myogenesis

Expression of the differentiated skeletal muscle phenotype is a process that appears to occur in at least two stages. First, pluripotent stem cells become committed to the myogenic lineage. Although undifferentiated and capable of continued proliferation, determined myoblasts are restricted to a single developmental fate. Upon receiving the appropriate environmental signals, these determined myoblasts withdraw from the cell cycle, fuse to form multi-nucleated myotubes, and begin to express a battery of muscle-specific gene products that make up the functional and contractile apparatus of the muscle. This project is aimed at the identification and characterization of factors that control the determination and differentiation of myogenic cells. We have cloned a cDNA, called myogenin, that plays an important role in these processes. Myogenin is expressed exclusively in skeletal muscle in vivo and myogenic cell lines in vitro. Its expression is sharply upregulated during differentiation. When constitutively expressed in fibroblasts, myogenin converts these cells to the myogenic lineage. Transfected cells behave as myogenic tissue culture cells with respect to the genes they express, the way they respond to environmental cues, and are capable of fusing to form multinucleated myotubes. Sequence analysis showed that this cDNA has homology to a family of transcription factors in a region of 72 amino acids known as the basic helix-loop-helix motif. This domain appears to mediate binding to a DNA sequence element known as an E-box (CANNTG) essential for the activity of the enhancers of many muscle-specific genes.^ Analysis of myogenin in tissue culture cells showed that its expression is responsive to many of the environmental cues, such as the presence of growth factors and oncogenes, that modulate myogenesis. In an attempt to identify the cis- and trans-elements that control myogenin expression and thereby understand what factors are responsible for the establishment of the myogenic lineage, we have cloned the myogenin gene. After analysis of the gene structure, we constructed a series of reporter constructs from the 5$\prime$ upstream sequence of the myogenin gene to determine which cis-acting sequences might be important in myogenin regulation. We found that 184 nucleotides of the 5$\prime$ sequence was sufficient to direct high-level muscle-specific expression of the reporter gene. Two sequence elements present in the 184 fragment, an E-box and a MEF-2 site, have been shown previously to be important in muscle-specific transcription. Mutagenesis of these sites revealed that both sites are necessary for full activity of the myogenin promoter, and suggests that a complex hierarchy of transcription factors control myogenic differentiation. ^

A cladistic analysis of apolipoprotein B gene variation and plasma lipid and apolipoprotein B levels, using familial data

Any functionally important mutation is embedded in an evolutionary matrix of other mutations. Cladistic analysis, based on this, is a method of investigating gene effects using a haplotype phylogeny to define a set of tests which localize causal mutations to branches of the phylogeny. Previous implementations of cladistic analysis have not addressed the issue of analyzing data from related individuals, though in human studies, family data are usually needed to obtain unambiguous haplotypes. In this study, a method of cladistic analysis is described in which haplotype effects are parameterized in a linear model which accounts for familial correlations. The method was used to study the effect of apolipoprotein (Apo) B gene variation on total-, LDL-, and HDL-cholesterol, triglyceride, and Apo B levels in 121 French families. Five polymorphisms defined Apo B haplotypes: the signal peptide Insertion/deletion, Bsp 1286I, XbaI, MspI, and EcoRI. Eleven haplotypes were found, and a haplotype phylogeny was constructed and used to define a set of tests of haplotype effects on lipid and apo B levels.^ This new method of cladistic analysis, the parametric method, found significant effects for single haplotypes for all variables. For HDL-cholesterol, 3 clusters of evolutionarily-related haplotypes affecting levels were found. Haplotype effects accounted for about 10% of the genetic variance of triglyceride and HDL-cholesterol levels. The results of the parametric method were compared to those of a method of cladistic analysis based on permutational testing. The permutational method detected fewer haplotype effects, even when modified to account for correlations within families. Simulation studies exploring these differences found evidence of systematic errors in the permutational method due to the process by which haplotype groups were selected for testing.^ The applicability of cladistic analysis to human data was shown. The parametric method is suggested as an improvement over the permutational method. This study has identified candidate haplotypes for sequence comparisons in order to locate the functional mutations in the Apo B gene which may influence plasma lipid levels. ^

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. ^

Identification of two vitamin D response elements in the rat osteopontin gene

Studies to elucidate the function of vitamin D have demonstrated an important role in regulating bone-related cells, including osteoblasts and osteoclasts. A seemingly paradoxical observation is that 1,25(OH)$\sb2$D$\sb3$, the active metabolite of vitamin D, stimulates bone resorption, yet regulates transcription of genes expressed by osteoblasts. One mechanism that could explain these actions is the upregulation of transcription of osteoblast-specific genes. These gene products could then act as effectors to influence osteoclastic activity. We hypothesized that molecular signals could be deposited directly into the mineralized matrix in the form of noncollagenous proteins, such as osteopontin (OPN). The structure, biosynthesis and localization of OPN suggest that it could function to mediate the molecular "cross talk" between osteoblasts and osteoclasts in response to 1,25(OH)$\sb2$D$\sb3$. To begin to address this hypothesis, elucidation of the molecular mechanisms of action involved in the transactivation of OPN by 1,25(OH)$\sb2$D$\sb3$ is essential.^ In the present study, the rat opn gene was isolated and characterized. Functional analysis by transient transfection of the 5$\sp\prime$ flanking sequences of the rat opn gene fused to the luciferase gene demonstrated that OPN is transcriptionally upregulated by 1,25(OH)$\sb2$D$\sb3$, mediated through two vitamin D response elements (VDRE). Both proximal and distal VDREs are structurally similar (two imperfect direct repeats separated by a 3 nucleotide spacer) and bind protein complexes that include the VDR and retinoid-X receptor (RXR). Isolated VDRE expression constructs produce functional activity of equivalent magnitude of responsiveness to 1,25(OH)$\sb2$D$\sb3$. However, expression constructs containing either VDRE and at least 200 bp of 5$\sp\prime$ and 3$\sp\prime$ flanking sequence demonstrated that the distal VDRE produces an amplitude of response significantly higher than the proximal VDRE. We conclude that the transcriptional upregulation of the opn gene by 1,25(OH)$\sb2$D$\sb3$ involves the transactivation of two VDREs, while maximal responsiveness requires interaction of the VDREs with additional cis-elements contained in the 5$\sp\prime$ sequence. ^

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. ^

Models of DNA sequence evolution

Models of DNA sequence evolution and methods for estimating evolutionary distances are needed for studying the rate and pattern of molecular evolution and for inferring the evolutionary relationships of organisms or genes. In this dissertation, several new models and methods are developed.^ The rate variation among nucleotide sites: To obtain unbiased estimates of evolutionary distances, the rate heterogeneity among nucleotide sites of a gene should be considered. Commonly, it is assumed that the substitution rate varies among sites according to a gamma distribution (gamma model) or, more generally, an invariant+gamma model which includes some invariable sites. A maximum likelihood (ML) approach was developed for estimating the shape parameter of the gamma distribution $(\alpha)$ and/or the proportion of invariable sites $(\theta).$ Computer simulation showed that (1) under the gamma model, $\alpha$ can be well estimated from 3 or 4 sequences if the sequence length is long; and (2) the distance estimate is unbiased and robust against violations of the assumptions of the invariant+gamma model.^ However, this ML method requires a huge amount of computational time and is useful only for less than 6 sequences. Therefore, I developed a fast method for estimating $\alpha,$ which is easy to implement and requires no knowledge of tree. A computer program was developed for estimating $\alpha$ and evolutionary distances, which can handle the number of sequences as large as 30.^ Evolutionary distances under the stationary, time-reversible (SR) model: The SR model is a general model of nucleotide substitution, which assumes (i) stationary nucleotide frequencies and (ii) time-reversibility. It can be extended to SRV model which allows rate variation among sites. I developed a method for estimating the distance under the SR or SRV model, as well as the variance-covariance matrix of distances. Computer simulation showed that the SR method is better than a simpler method when the sequence length $L>1,000$ bp and is robust against deviations from time-reversibility. As expected, when the rate varies among sites, the SRV method is much better than the SR method.^ The evolutionary distances under nonstationary nucleotide frequencies: The statistical properties of the paralinear and LogDet distances under nonstationary nucleotide frequencies were studied. First, I developed formulas for correcting the estimation biases of the paralinear and LogDet distances. The performances of these formulas and the formulas for sampling variances were examined by computer simulation. Second, I developed a method for estimating the variance-covariance matrix of the paralinear distance, so that statistical tests of phylogenies can be conducted when the nucleotide frequencies are nonstationary. Third, a new method for testing the molecular clock hypothesis was developed in the nonstationary case. ^

Transcriptional regulation of the mouse $\alpha$2(I) collagen gene

The mouse $\alpha$2(I) collagen gene is specifically expressed in a limited number of cell types in the body including fibroblasts and osteoblasts. We had previously shown that a promoter containing the sequences between $-$350 and +54 bp was expressed at low levels in a cell- and tissue-specific fashion in transgenic mice. Further studies suggested that the sequence between $-$315 and $-$284 bp could mediate cell- and tissue-specific expression of reporter genes in cell culture and in transgenic mice. We report here characterization of the proteins binding to this segment and propose a model for the cell-specific expression conferred by this sequence. In this study we also identified a strong enhancer for the mouse $\alpha$2(I) collagen gene located approximately 13.5 to 19.5 kb upstream of the transcriptional start site. This enhancer segment is characterized by the presence of three cell-specific hypersensitive sites and can drive high levels of cell-specific expression of a heterologous 220-bp mouse $\alpha$1(I) collagen promoter. In the course of this study, we identified a novel zinc finger transcription factor (designated murine epithelial zinc finger, mEZF) which was transiently expressed in the mesenchymal cells which give rise to the skeletal primordia and the metanephric kidney during the early stages of embryogenesis. In newborn mice, the mEZF gene is expressed at high levels in differentiated epithelial cells of the skin, oral mucosa, tongue, esophagus, stomach and colon. Chromosomal mapping suggested that the mEZF gene mapped to mouse Chromosome 4 and that the human homolog of mEZF would likely map to human Chromosome 9q31. This region of the human genome contains tumor suppressor genes for basal cell carcinomas of the skin as well as for squamous cell carcinomas of various organs. We cloned and characterized the human homolog of mEZF and mapped its chromosomal position as a first step in determining whether or not this gene plays a role in the development of these tumors. ^

The story of RP10: A positional cloning approach to the identification of an autosomal dominant retinitis pigmentosa gene

Retinitis pigmentosa (RP) is an inherited retinal degenerative disease that is the leading cause of inherited blindness worldwide. Characteristic features of the disease include night blindness, progressive loss of visual fields, and deposition of pigment on the retina in a bone spicule-like pattern. RP is marked by extreme genetic heterogeneity with at least 19 autosomal dominant, autosomal recessive and X-linked loci identified. RP10, which maps to chromosome 7q, was the fifth autosomal dominant RP locus identified, and accounts for the early-onset disease in two independent families. Extensive linkage and haplotype analyses have been performed in these two families which have allowed the assignment of the disease locus to a 5-cM region on chromosome 7q31.3. In collaboration with Dr. Eric Green (National Center for Human Genome Research, National Institutes of Health), a well-characterized physical map of the region was constructed which includes YAC, BAC and cosmid coverage. The entire RP10 critical region resides within a 9-Mb well-characterized YAC contig. These physical maps not only provided the resources to undertake the CAIGES (cDNA amplification for identification of genomic expressed sequences) procedure for identification of retinal candidate genes within the critical region, but also identified a number of candidate genes, including transducin-$\gamma$ and blue cone pigment genes. All candidate genes examined were excluded. In addition, a number of ESTs were mapped within the critical region. EST20241, which was isolated from an eye library, corresponded to the 3$\sp\prime$ region of the ADP-ribosylation factor (ARF) 5 gene. ARF5, with its role in vesicle transport and possible participation in the regulation of the visual transduction pathway, became an extremely interesting candidate gene. Using a primer walking approach, the entire 3.2 kb genomic sequence of the ARF5 gene was generated and developed intronic primers to screen for coding region mutations in affected family members. No mutations were found in the ARF5 gene, however, a number of additional ESTs have been mapped to the critical region, and, as the large-scale sequencing projects get underway, megabases of raw sequence data from the RP10 region are becoming available. These resources will hasten the isolation and characterization of the RP10 gene. ^

Accessing novel developmental mechanisms in the mouse by gene trapping

Genetic analysis is a powerful method for analyzing the function of specific genes in development. I sought to identify novel genes in the mouse using a genetic analysis relying on the expression pattern and phenotype of mutated genes. To this end, I have conducted a gene trap screen using the vector $\rm SA\beta geo,$ a promoterless DNA construct that encodes a fusion protein with lacZ and neomycin resistance activities. Productive integration and expression of the $\beta$geo protein in embryonic stem (ES) cells requires integration into an active transcription unit. The endogenous regulatory elements direct reporter gene expression which reflects the expression of the endogenous gene. Of eight mouse lines generated from gene trap ES cell clones, four showed differential regulation of $\beta$geo activity during embryogenesis. These four were analyzed in more detail.^ Three of the lines RNA 1, RNA2 and RNA 3 had similar expression patterns, within subsets of cells in sites of embryonic hematopoiesis. Cloning of the trapped genes revealed that all three integrations had occurred within 45S rRNA precursor transcription units. These results imply that there exists in these cells some mechanism responsible for the efficient production of the $\beta$geo protein from an RNA polymerase I transcript that is not present in most of the cells in the embryo.^ The fourth line, GT-2, showed widespread, dynamic expression. Many of the sites of expression were important classic embryonic induction systems. Cloning of the sequences fused to the $5\sp\prime$ end of the $\beta$geo sequence revealed that the trapped gene contained significant sequence homology with a previously identified human sequence HumORF5. An open reading frame of this sequence is homologous to a group of eukaryotic proteins that are members of the RNA helicase superfamily I.^ Analysis of the gene trap lines suggests that potentially novel developmental mechanisms have been uncovered. In the case of RNA 1, 2 and 3, the differential production of ribosomal RNAs may be required for differentiation or function of the $\beta$geo positive hematopoietic cells. In the GT-2 line, a previously unsuspected temporal and spatial regulation of a putative RNA helicase implies a role for this activity during specific aspects of mouse development. ^

Fibrillins: Sequence, phylogeny and carboxy terminal domain proteolytic processing

Fibrillin-1 and -2 are large secreted glycoproteins that are known to be components of extracellular matrix microfibrils located in the vasculature, basement membrane and various connective tissues. These microfibrils are often associated with a superstructure known as the elastic fiber. During the development of elastic tissues, fibrillin microfibrils precede the appearance of elastin and may provide a scaffolding for the deposition and crosslinking of elastin. Using RT/PCR, we cloned and sequenced 3.85Kbp of the FBN2 gene. Five differences were found between our contig sequence and that published by Zhang et al. (1995). Like many extracellular matrix proteins, the fibrillins are modular proteins. We compared analogous domains of the two fibrillins and also members of the latent TGF-$\beta$ binding protein (LTBP) family to determine their phylogenetic relationship. We found that the two families are homologous. LTBP-2 is the most similar to the fibrillin family while FBN-1 is the most similar to the LTBP family. The fibrillin-1 carboxy terminal domain is proteolytically processed. Two eukaryotic protein expression systems, baculoviral and CHO-K1, were developed to examine the proteolytic processing of the carboxy terminal domain of the fibrillin-1 protein. Both expression systems successfully processed the domain and both processed a mutant less efficiently. In the CHO-K1 cells, processing occurred intracellularly. ^

Molecular mechanisms of chondrocyte-specific expression on the mouse pro$\alpha$1(II) collagen gene

Type II collagen is a major chondrocyte-specific component of the cartilage extracellular matrix and it represents a typical differentiation marker of mature chondrocytes. In order to delineate cis-acting elements of the mouse pro$\alpha1$(II) collagen gene that control chondrocyte-specific expression in intact mouse embryos, we generated transgenic mice harboring chimeric constructions in which varying lengths of the promoter and intron 1 sequences were linked to a $\beta$-galactosidase reporter gene. A construction containing a 3000-bp promoter and a 3020-bp intron 1 fragment directed high levels of $\beta$-galactosidase expression specifically to chondrocytes. Successive deletions of intron 1 delineated a 48-bp fragment which targeted $\beta$-galactosidase expression to chondrocytes with the same specificity as the larger intron 1 fragment. When the Col2a1 promoter was replaced with a minimal $\beta$-globin promoter, the 48-bp intron 1 sequence was still able to target expression of the transgene to chondrocytes, specifically. Therefore a 48-bp intron 1 DNA segment of the mouse Col2a1 gene contains the necessary information to confer high-level, temporally correct, chondrocyte expression to a reporter gene in intact mouse embryos and that Col2a1 promoter sequences are dispensable for chondrocyte expression. Nuclear proteins present selectively in mouse primary chondrocytes and rat chondrosarcoma cells bind to the three putative HMG (High-Mobility-Group) domain protein binding sites in this 48-bp sequence and the chondrocyte-specific proteins likely bind the DNA through minor groove. Together, my results indicate that a 48-bp sequence in Col2a1 intron 1 controls chondrocyte-specific expression in vivo and suggest that chondrocytes contain specific nuclear proteins involved in enhancer activity. ^

Transcriptional regulation and functional analysis of the human Wilms' tumor 1 gene

The Wilms' tumor 1 gene (WT1) encodes a zinc-finger transcription factor and is expressed in urogenital, hematopoietic and other tissues. It is expressed in a temporal and spatial manner in both embryonic and adult stages. To obtain a better understanding of the biological function of WT1, we studied two aspects of WT1 regulation: one is the identification of tissue-specific cis-regulatory elements that regulate its expression, the other is the downstream genes which are modulated by WT1.^ My studies indicate that in addition to the promoter, other regulatory elements are required for the tissue specific expression of this gene. A 259-bp hematopoietic specific enhancer in intron 3 of the WT1 gene increased the transcriptional activity of the WT1 promoter by 8- to 10-fold in K562 and HL60 cells. Sequence analysis revealed both GATA and c-Myb motifs in the enhancer fragment. Mutation of the GATA motif decreased the enhancer activity by 60% in K562 cells. Electrophoretic mobility shift assays showed that both GATA-1 and GATA-2 proteins in K562 nuclear extracts bind to this motif. Cotransfection of the enhancer containing reporter construct with a GATA-1 or GATA-2 expression vector showed that both GATA-1 and GATA-2 transactivated this enhancer, increasing the CAT reporter activity 10-15 fold and 5-fold respectively. Similar analysis of the c-Myb motif by cotransfection with the enhancer CAT reporter construct and a c-Myb expression vector showed that c-Myb transactivated the enhancer by 5-fold. A DNase I-hypersensitive site has been identified in the 258 bp enhancer region. These data suggest that GATA-1 and c-Myb are responsible for the activity of this enhancer in hematopoietic cells and may bind to the enhancer in vivo. In the process of searching for cis-regulatory elements in transgenic mice, we have identified a 1.0 kb fragment that is 50 kb downstream from the promoter and is required for the central nervous system expression of WT1.^ In the search for downstream target genes of WT1, we noted that the proto-oncogene N-myc is coexpressed with the tumor suppressor gene WT1 in the developing kidney and is overexpressed in many Wilms' tumors. Sequence analysis revealed eleven consensus WT1 binding sites located in the 1 kb mouse N-myc promoter. We further showed that the N-myc promoter was down-regulated by WT1 in transient transfection assays. Electrophoretic mobility shift assays showed that oligonucleotides containing the WT1 motifs could bind WT1 protein. Furthermore, a Denys-Drash syndrome mutant of WT1, R394W, that has a mutation in the DNA binding domain, failed to repress the N-myc promoter. This suggests that the repression of the N-myc promoter is mediated by DNA binding of WT1. This finding helps to elucidate the relationship of WT1 and N-myc in tumorigenesis and renal development. ^

Regulatory factors and control of anthrax toxin gene expression

Bacillus anthracis plasmid pXO1 carries genes for three anthrax toxin proteins, pag (protective antigen), cya (edema factor), and lef (lethal factor). Expression of the toxin genes is enhanced by two signals: CO$\sb2$/bicarbonate and temperature. The CO$\sb2$/bicarbonate effect requires the presence of pXO1. I hypothesized that pXO1 harbors a trans-acting regulatory gene(s) required for CO$\sb2$/bicarbonate-enhanced expression of the toxin genes. Characterization of such a gene(s) will lead to increased understanding of the mechanisms by which B. anthracis senses and responds to host environments.^ A regulatory gene (atxA) on pXO1 was identified. Transcription of all three toxin genes is decreased in an atxA-null mutant. There are two transcriptional start sites for pag. Transcription from the major site, P1, is enhanced in elevated CO$\sb2$. Only P1 transcripts are significantly decreased in the atxA mutant. Deletion analysis of the pag upstream region indicates that the 111-bp region upstream of the P1 site is sufficient for atxA-mediated increase of this transcript. The cya and lef genes each have one apparent transcriptional start site. The cya and lef transcripts are significantly decreased in the atxA mutant. The atxA mutant is avirulent in mice. The antibody response to all three toxin proteins is significantly decreased in atxA mutant-infected mice. These data suggest that the atxA gene product activates expression of the toxin genes and is essential for virulence.^ Since expression of the toxin genes is dependent on atxA, whether increased toxin gene expression in response to CO$\sb2$/bicarbonate and temperature is associated with increased atxA expression was investigated. I monitored steady state levels of atxA mRNA and AtxA protein in different growth conditions. The results indicate that expression of atxA is not influenced by CO$\sb2$/bicarbonate. Steady state levels of atxA mRNA and AtxA protein are higher at 37$\sp\circ$C than 28$\sp\circ$C. However, increased pag expression at high temperature can not be attributed directly to increased atxA expression.^ There is evidence that an additional factor(s) may be involved in regulation of pag. Expression of pag in strains overproducing AtxA is significantly decreased compared to the wildtype strain. A specific interaction of tagged-AtxA with the pag upstream DNA has not been demonstrated. Furthermore, four proteins in B. anthracis extract can be co-immunoprecipitated with tagged-AtxA. Amino-terminal sequence of one protein has been determined and found highly homologous to chaperonins of GroEL family. Studies are under way to determine if this GroEL-like protein interactions with AtxA and plays any role in atxA-mediated activation of toxin genes. ^

Transcriptional regulation of the human {\it PAX6\/} gene

PAX6, a member of the paired-type homeobox gene family, is expressed in a partially and temporally restricted pattern in the developing central nervous system, and its mutation is responsible for human aniridia (AN) and mouse small eye (Sey). The objective of this study was to characterize the PAX6 gene regulation at the transcriptional level, and thereby gain a better understanding of the molecular basis of the dynamic expression pattern and the diversified function of the human PAX6 gene.^ Initially, we examined the transcriptional regulation of the PAX6 gene by transient transfection assays and identified multiple cis-regulatory elements that function differently in different cell lines. The transcriptional initiation site was identified by RNase protection and primer extension assays. Examination of the genomic DNA sequence indicated that the PAX6 promoter has a TATA like-box (ATATTTT) at $-$26 bp, and two CCAAT-boxes are located at positions $-$70 and $-$100 bp. A 38 bp ply (CA) sequence was located 992 bp upstream from the initiation site. Transient transfection assays in glioblastoma cells and leukemia cells indicate that a 92 bp region was required for basal level PAX6 promoter activity. Gel retardation assays showed that this 92 bp sequence can form four DNA-protein complexes which can be specifically competed by a 31-mer oligonucleotide containing a PAX6 TATA-like sequence or an adenovirus TATA box. The activation of the promoter is positively correlated with the expression of PAX6 transcripts in cells tested.^ Based on the results obtained from the in vitro transfection assays, we did further dissection assay and functional analysis in both cell-culture and transgenic mice. We found that a 5 kb upstream promoter sequence is required for the tissue specific expression in the forebrain region which is consistent with that of the endogenous PAX6 gene. A 267 bp cell-type specific repressor located within the 5 kb fragment was identified and shown to direct forebrain specific expression. The cell-type specific repressor element has been narrowed to a 30 bp region which contains a consensus E-box by in vitro transfection assays. The third regulatory element identified was contained in a 162 bp sequence (+167 to +328) which functions as a midbrain repressor, and it appeared to be required for establishing the normal expression pattern of the PAX6 gene. Finally, a highly conserved 216 bp sequence identified in intron 4 exhibited as a spinal cord specific enhancer. And this 216 bp cis-regulatory element can be used as a marker to trace the differentiation and migration of progenitor cells in the developing spinal cord. These studies show that the concerted action of multiple cis-acting regulatory elements located upstream and downstream of the transcription initiation site determines the tissue specific expression of PAX6 gene. ^