Cis-acting elements and developmental regulators govern toxin gene expression in Bacillus anthracis

Expression of the structural genes for the anthrax toxin proteins is coordinately controlled by host-related signals such as elevated CO2 , and the trans-acting positive regulator, AtxA. No specific binding of AtxA to the toxin gene promoters has been demonstrated and no sequence-based similarities are apparent in the promoter regions of toxin genes. We hypothesized that the toxin genes possess common structural features that are required for positive regulation. To test this hypothesis, I performed an extensive characterization of the toxin gene promoters. I determined the minimal sequences required for atxA-mediated toxin gene expression and compared these sequences for structural similarities. In silico modeling and in vitro experiments indicated significant curvature within these regions. Random mutagenesis revealed that point mutations associated with reduced transcriptional activity, mostly mapped to areas of high curvature. This work enabled the identification of two potential cis-acting elements implicated in AtxA-mediated regulation of the toxin genes. In addition to the growth condition requirements and AtxA, toxin gene expression is under growth phase regulation. The transition state regulator AbrB represses atxA expression to influence toxin synthesis. Here I report that toxin gene expression also requires sigH, a gene encoding the RNA polymerase sigma factor associated with development in B. subtilis. In the well-studied B. subtilis system, σH is part of a feedback control pathway that involves AbrB and the major response regulator of sporulation initiation, Spo0A. My data indicate that in B. anthracis, regulatory relationships exist between these developmental regulators and atxA . Interestingly, during growth in toxin-inducing conditions, sigH and abrB expression deviates from that described for B. subtilis, affecting expression of the atxA gene. These findings, combined with previous observations, suggest that the steady state level of atxA expression is critical for optimal toxin gene transcription. I propose a model whereby, under toxin-inducing conditions, control of toxin gene expression is fine-tuned by the independent effects of the developmental regulators on the expression of atxA . The growth condition-dependent changes in expression of these regulators may be crucial for the correct timing and uninterrupted expression of the toxin genes during infection. ^

The molecular basis for beta-lactamase gene expression in B. anthracis, B. cereus and B. thuringiensis

The susceptibility of most Bacillus anthracis strains to β-lactam antibiotics is intriguing considering that the B. anthracis genome harbors two β-lactamase genes, bla1 and bla2, and closely-related species, Bacillus cereus and Bacillus thuringiensis, typically produce β-lactamases. This work demonstrates that B. anthracis bla expression is affected by two genes, sigP and rsp, predicted to encode an extracytoplasmic function sigma factor and an antisigma factor, respectively. Deletion of the sigP/rsp locus abolished bla expression in a penicillin-resistant clinical isolate and had no effect on bla expression in a prototypical penicillin-susceptible strain. Complementation with sigP/rsp from the penicillin-resistant strain, but not the penicillin-susceptible strain, conferred β-lactamase activity upon both mutants. These results are attributed to a nucleotide deletion near the 5' end of rsp in the penicillin-resistant strain that is predicted to result in a nonfunctional protein. B. cereus and B. thuringiensis sigP and rsp homologues are required for inducible penicillin resistance in those species. Expression of the B. cereus or B. thuringiensis sigP and rsp genes in a B. anthracis sigP/rsp-null mutant confers resistance to β-lactam antibiotics, suggesting that while B. anthracis contains the genes necessary for sensing β-lactam antibiotics, the B. anthracis sigP/rsp gene products are insufficient for bla induction. ^ Because alternative sigma factors recognize unique promoter sequence, direct targets can be elucidated by comparing transcriptional profiling results with an in silico search using the sigma factor binding sequence. Potential σP -10 and -35 promoter elements were identified upstream from bla1 bla2 and sigP. Results obtained from searching the B. anthracis genome with the conserved sequences were evaluated against transcriptional profiling results comparing B. anthracis 32 and an isogenic sigP/rsp -null strain. Results from these analyses indicate that while the absence of the sigP gene significantly affects the transcript levels of 16 genes, only bla1, bla2 and sigP are directly regulated by σP. The genomes of B. cereus and B. thuringiensis strains were also analyzed for the potential σP binding elements. The sequence was located upstream from the sigP and bla genes, and previously unidentified genes predicted to encode a penicillin-binding protein (PBP) and a D-alanyl-D-alanine carboxypeptidase, indicating that the σ P regulon in these species responds to cell-wall stress caused by β-lactam antibiotics. ^ β-lactam antibiotics prevent attachment of new peptidoglycan to the cell wall by blocking the active site of PBPs. A B. cereus and B. thuringiensis pbp-encoding gene located near bla1 contains a potential σP recognition sequence upstream from the annotated translational start. Deletion of this gene abolished β-lactam resistance in both strains. Mutations in the active site of the PBP were detrimental to β-lactam resistance in B. cereus, but not B. thuringiensis, indicating that the transpeptidase activity is only important in B. cereus. I also found that transcript levels of the PBP-encoding gene are not significantly affected by the presence of β-lactam antibiotic. Based on these data I hypothesize that the gene product acts a sensor of β-lactam antibiotic. ^

Transcriptional and translational regulators of gene expression in germ cell development

Germ cell development is a highly coordinated process driven, in part, by regulatory mechanisms that control gene expression. Not only transcription, but also translation, is under regulatory control to direct proper germ cell development. In this dissertation, I have focused on two regulators of germ cell development. One is the homeobox protein RHOX10, which has the potential to be both a transcriptional and translational regulator in mouse male germ cell development. The other is the RNA-binding protein, Hermes, which functions as a translational regulator in Xenopus laevis female germ cell development. ^ Rhox10 is a member of reproductive homeobox gene X-(linked (Rhox) gene cluster, of which expression is developmentally regulated in developing mouse testes. To identify the cell types and developmental stages in which Rhox10 might function, I characterized its temporal and spatial expression pattern in mouse embryonic, neonatal, and adult tissues. Among other things, this analysis revealed that both the level and the subcellular localization of RHOX10 are regulated during germ cell development. To understand the role of Rhox10 in germ cell development, I generated transgenic mice expressing an artificial microRNA (miRNA) targeting Rhox10. While this artificial miRNA robustly downregulated RHOX10 protein expression in vitro, it did not significantly reduce RHOX10 expression in vivo. So I next elected to knockdown RHOX10 levels in spermatogonial stem cells (SSCs), which I found highly express both Rhox10 mRNA and RHOX10 protein. Using a recently developed in vitro culture system for SSCs combined with a short-hairpin RNA (shRNA) approach, I strongly depleted RHOX10 expression in SSCs. These RHOX10-depleted cells exhibited a defect in the ability to form stem cell clusters in vitro. Expression profiling analysis revealed many genes regulated by Rhox10, including many meiotic genes, which could be downstream of Rhox10 in a molecular pathway that controls SSC differentiation. ^ RNA recognition motif (RRM) containing protein, Hermes is localized in germ plasm, where dormant mRNAs are also located, of Xenopus oocytes, which implicates its role in translational regulator. To understand the function of Hermes in oocyte meiosis, I used a morpholino oligonucleotide (MO) based knockdown approach. Microinjection of Hermes MO into fully grown oocytes, which are arrested in meiotic prophase, caused acceleration of oocytes reentry into meiosis (i.e., maturation) upon progesterone induction. Using a candidate approach, I identified at least three targets of Hermes: Ringo/Spy, Xcat2, and Mos. Ringo/Spy and Mos are known to have functions in oocyte maturation, while Ringo/Spy, Xcat2 mRNA are localized in the germ plasm of oocytes, which drives germ cell specification after fertilization. This led me to propose that Hermes functions in both oocyte maturation and germ cell development through its ability to regulate 3 crucial target mRNAs. ^

XIPHOPHORUS ENZYME GENE EXPRESSION: TISSUE SPECIFICITY, GENETIC CONTROL, AND STABILITY IN CULTURED CELLS

Five permanent cell lines were developed from Xiphophorus maculatus, X. helleri, and their hybrids using three tissue sources, including adults and embryos of different stages. To evaluate cell line gene expression for retention of either tissue-of-origin-specific or ontogenetic stage-specific characters, the activity distribution of 44 enzyme loci was determined in 11 X. maculatus tissues, and the developmental genetics of 17 enzyme loci was charted in X. helleri and in helleri x maculatus hybrids using starch gel electrophoresis. In the process, eight new loci were discovered and characterized for Xiphophorus.^ No Xiphophorus cell line showed retention of tissue-of-origin-specific or ontogenetic stage-specific enzyme gene expressional traits. Instead, gene expression was similar among the cell lines. One enzyme, adenosine deaminase (ADA) was an exception. Two adult-origin cell lines expressed ADA, whereas, three cell lines derived independently from embryos did not. ADA('-) expression of Xiphophorus embryo-derived cell lines may represent retention of an embryonic gene expressional trait. In one cell line (T(,3)) derived from 13 pooled interspecific hybrid (F(,2)) embryos, shifts with time were observed at enzyme loci polymorphic between the two species. This suggested shifts in ratios of cells of different genotypes in the population rather than unstable gene expression in one dominant cell type.^ Verification of this hypothesis was attempted by cloning the culture--seeking clones having different genetic signatures. The large number of loci electrophoretically polymorphic between the two species and whose alleles segregated independently into the 13 progeny from which this culture originated almost guaranteed the presence of different genetic signatures (lineages) in T(,3).^ Seven lineages of cells were found within T(,3), each expressing genotypes at some loci not characteristic of the expression of the culture-as-a-whole, supporting the hypothesis tested. Quantitative studies of ADA expression in the whole culture (ADA('-)) and in clones of these seven lineages suggested the predominance in T(,3) of ADA deficient cell lineages, although moderate to high ADA output clones also occurred. Thus, T(,3) has the potential to shift phenotypes from ADA('-) to ADA('+) by simply changing proportions of its constituent cell types, demonstrating that such shifts can occur in any cell culture containing cells of mixed expressional characteristics.^

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

Estrogen-induced gene expression and patient survival in high-grade serous carcinoma of the ovary

Background. Assessment of estrogen receptor (ER) expression has inconsistent utility as a prognostic marker in epithelial ovarian carcinoma. In breast and endometrial cancers, the use of estrogen-induced gene panels, rather than ER expression alone, has shown improved prognostic capability. Specifically, over-expression of estrogen-induced genes in these tumors is associated with a better prognosis and signifies estrogen sensitivity that can be exploited with hormone antagonizing agents. It was therefore hypothesized that estrogen-induced gene expression in ovarian carcinoma would successfully predict outcomes and differentiate between tumors of varying estrogen sensitivities. Methods. Two hundred nineteen (219) patients with ovarian cancer who underwent surgery at M. D. Anderson between 2004 and 2007 were identified. Of these, eighty-three (83) patients were selected for inclusion because they had advanced stage, high-grade serous carcinoma of the ovary or peritoneum, had not received neoadjuvant chemotherapy, and had readily available frozen tissue for study. All patients had also received adjuvant treatment with platinum and taxane agents. The expression of seven genes known to be induced by estrogen in the female reproductive tract (EIG121, sFRP1, sFRP4, RALDH2, PR, IGF-1, and ER) was measured using qRT-PCR. Unsupervised cluster analyses of multiple gene permutations were used to categorize patients as high or low estrogen-induced gene expressors. QPCR gene expression results were then compared to ER and PR immunohistochemical (IHC) expression. Cox proportional hazards models were used to evaluate the effects of both individual genes and selected gene clusters on patient survival. Results. Median follow-up time was 38.7 months (range 1-68 months). In a multivariate model, overall survival was predicted by sFRP1 expression (HR 1.10 [1.02-1.19], p=0.01) and EIG121 expression (HR 1.28 [1.10-1.49], p<0.01). A cluster defined by EIG121 and ER was further examined because that combination appeared to reasonably segregate tumors into distinct groups of high and low estrogen-induced gene expressors. Shorter overall survival was associated with high estrogen-induced gene expressors (HR 2.84 [1.11-7.30], p=0.03), even after adjustment for race, age, body mass index, and residual disease at debulking. No difference in IHC ER or PR expression was noted between gene clusters. Conclusion. In sharp contrast to breast and endometrial cancers, high estrogen-induced gene expression predicts shorter overall survival in patients with high-grade serous ovarian carcinoma. An estrogen-induced gene biomarker panel may have utility as prognostic indicator and may be useful to guide management with estrogen antagonists in this population.^

DETERMINING THE GENOTYPE-PHENOTYPE CONNECTION IN SYNTHETIC INDUCIBLE GENE EXPRESSION SYSTEMS

Introduction Gene expression is an important process whereby the genotype controls an individual cell’s phenotype. However, even genetically identical cells display a variety of phenotypes, which may be attributed to differences in their environment. Yet, even after controlling for these two factors, individual phenotypes still diverge due to noisy gene expression. Synthetic gene expression systems allow investigators to isolate, control, and measure the effects of noise on cell phenotypes. I used mathematical and computational methods to design, study, and predict the behavior of synthetic gene expression systems in S. cerevisiae, which were affected by noise. Methods I created probabilistic biochemical reaction models from known behaviors of the tetR and rtTA genes, gene products, and their gene architectures. I then simplified these models to account for essential behaviors of gene expression systems. Finally, I used these models to predict behaviors of modified gene expression systems, which were experimentally verified. Results Cell growth, which is often ignored when formulating chemical kinetics models, was essential for understanding gene expression behavior. Models incorporating growth effects were used to explain unexpected reductions in gene expression noise, design a set of gene expression systems with “linear” dose-responses, and quantify the speed with which cells explored their fitness landscapes due to noisy gene expression. Conclusions Models incorporating noisy gene expression and cell division were necessary to design, understand, and predict the behaviors of synthetic gene expression systems. The methods and models developed here will allow investigators to more efficiently design new gene expression systems, and infer gene expression properties of TetR based systems.

Post-transcriptional Regulation of Mammalian Gene Expression in Non-coding Region of Target RNA

Tumor Suppressor Candidate 2 (TUSC2) is a novel tumor suppressor gene located in the human chromosome 3p21.3 region. TUSC2 mRNA transcripts could be detected on Northern blots in both normal lung and some lung cancer cell lines, but no endogenous TUSC2 protein could be detected in a majority of lung cancer cell lines. Mechanisms regulating TUSC2 protein expression and its inactivation in primary lung cancer cells are largely unknown. We investigated the role of the 5’- and 3’-untranslated regions (UTRs) of the TUSC2 gene in the regulation of TUSC2 protein expression. We found that two small upstream open-reading frames (uORFs) in the 5’UTR of TUSC2 could markedly inhibit the translational initiation of TUSC2 protein by interfering with the “scanning” of the ribosome initiation complexes. Site-specific stem-loop array reverse transcription-polymerase chain reaction (SLA-RT-PCR) verified several micoRNAs (miRNAs) targeted at 3’UTR and directed TUSC2 cleavage and degradation. In addition, we used the established let-7-targeted high mobility group A2 (Hmga2) mRNA as a model system to study the mechanism of regulation of target mRNA by miRNAs in mammalian cells under physiological conditions. There have been no evidence of direct link between mRNA downregulation and mRNA cleavages mediated by miRNAs. Here we showed that the endonucleolytic cleavages on mRNAs were initiated by mammalian miRNA in seed pairing style. Let-7 directed cleavage activities among the eight predicted potential target sites have varied efficiency, which are influenced by the positional and the structural contexts in the UTR. The 5’ cleaved RNA fragments were mostly oligouridylated at their 3’-termini and accumulated for delayed 5’–3’ degradation. RNA fragment oligouridylation played important roles in marking RNA fragments for delayed bulk degradation and in converting RNA degradation mode from 3’–5’ to 5’–3’ with cooperative efforts from both endonucleolytic and non-catalytic miRNA-induced silencing complex (miRISC). Our findings point to a mammalian miRNA-mediated mechanism for the regulation of mRNA that miRNA can decrease target mRNA through target mRNA cleavage and uridine addition

MUTATIONS IN STAT3 ASSOCIATED WITH HUMAN HYPER IgE SYNDROME ENHANCE NFκB AND MAPK-MEDIATED GENE EXPRESSION

Hyper IgE syndrome (HIES) is a multisystem disorder resulting in bone and immune system abnormalities. It is associated with mutations in STAT3, which disrupt protein domains responsible for transcriptional function. Patients with HIES display osteoporosis and enhanced inflammatory cytokine production similar to hematopoietic Stat3-deficient mice. Since osteoclast and inflammatory cytokine genes are NFκB targets, these observations indicate a possible deregulation of NFκB signaling in both mice and humans with STAT3-deficiency. Here, we sought to examine the role of STAT3 in the regulation of NFκB-mediated gene expression through analysis of three HIES STAT3 point mutations in both hematopoietic and non- hematopoietic cells. We found that IL-6-induced tyrosine phosphorylation of STAT3 was partially or completely abrogated by HIES mutations in the transactivation domain (V713L) or SH2 domain (V637M), respectively, in both hematopoietic and non- hematopoietic cells. By contrast, IL-6-induced tyrosine phosphorylation of an HIES mutant in the STAT3 DNA-binding domain (R382W) was intact. The R382W and V713L mutants significantly reduced IL-6-dependent STAT3 transcriptional activity in reporter gene assays. Moreover, the R382W and V637M mutants significantly diminished IL-6-responsive expression of the endogenous STAT3 target gene, Socs3, as assessed by quantitative real-time PCR (qPCR) in the RAW macrophage cell line. These observations indicate the HIES mutants dominantly suppress the transcriptional activity of wild type STAT3, albeit to varying degrees. All three HIES mutants enhanced LPS-induced expression of the NFκB target genes IL6 (IL-6), Cxcl10 (IP- 10), and Tnf (TNFα) in RAW cells, as indicated by qPCR. Furthermore, overexpression of wild type STAT3 in Stat3-deficient murine embryonic fibroblasts significantlyreduced LPS-stimulated expression of IL6, Cxcl10, and IL12p35. In addition, in aprimary murine osteoclast differentiation assay, a STAT3-specific SH2 domain inhibitor led to significantly increased levels of osteoclast-specific gene expression. These results suggest that STAT3 serves as a negative regulator of NFκB-mediated gene expression, and furthermore imply that STAT3 mutations associated with HIES contribute to the osteopenia and inflammation observed in HIES patients.

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.

{\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.) ^

Role and regulation of c-KIT gene expression in tumor growth and metastasis of human melanoma

The purpose of this study was to investigate the role of the c-KIT receptor in the progression of human melanoma and the mechanism(s) for the regulation of c-KIT gene expression in human melanoma.^ The molecular changes associated with the transition of melanoma cells from radial growth phase (RGP) to vertical growth phase (VGP) (metastatic phenotype) are not well-defined. Expression of the tyrosine-kinase receptor c-KIT progressively decreases during local tumor growth and invasion of human melanomas. To provide direct evidence that the metastasis of human melanoma is associated with the loss of c-KIT expression, highly metastatic A375SM cells, which express very low or undetectable levels of c-KIT, were tranduced with the human c-KIT gene. We demonstrated that enforced c-KIT expression in highly metastatic human melanoma cells significantly suppressed their tumorigenicity and metastatic propensity in nude mice. In addition, we showed that the ligand for c-KIT, SCF, induces apoptosis in human melanoma cells expressing c-KIT under both in vitro and in vivo conditions. These results suggest that loss of c-KIT receptor may allow malignant melanoma cells to escape SCF/c-KIT-mediated apoptosis, thus contributing to tumor growth and eventually metastasis.^ Furthermore, we investigated the possible mechanism(s) for the down-regulation of c-KIT gene expression in malignant melanoma. Sequence analysis of the c-KIT promoter indicated that this promoter contains several consensus binding-site sequences including three putative AP2 and two Myb sites. Although Myb was shown to be associated with c-KIT expression in human hemotopoietic cells, we found no correlation between c-KIT expression and Myb expression in human melanoma cell lines. In contrast, we showed that c-KIT expression directly correlates with expression of AP2 in human melanoma cells. We found that highly metastatic cells do not express the transcription factor AP2. Expression of AP2 in A375SM cells (c-KIT-negative and AP2-negative) was enough to restore luciferase activity driven by the c-KIT promoter in a dose-dependent manner. On the other hand, co-expression of the dominant-negative form of AP2 (AP2B) in Mel-501 cells (c-KIT-positive and AP2-positive) resulted in two-fold reduction in luciferase activity. Electrophoretic mobility shift assays revealed that the c-KIT promoter contains functional AP2 binding sites which could associate with AP2 protein. Endogenous c-KIT gene expression levels were elevated in AP2 stably-transfected human melanoma A375SM cells. Expression of exogenous AP2 in A375SM cells inhibited their tumorigenicity and metastatic potential in nude mice. The c-KIT ligand, SCF, also induced apoptosis in the AP2 stably-transfected A375SM cells. The identification of AP2 as an important regulator for c-KIT expression suggests that AP2 may have tumor growth and metastasis inhibitory properties, possibly mediated through c-KIT/SCF effects on apoptosis of human melanoma cells. Since AP2 binding sites were found in the promoters of other genes involved in the progression of human melanoma, such as MMP2 (72 kDa collagenase), MCAM/MUC18 and P21/WAF-1, our findings suggest that loss of AP2 expression might be a crucial event in the development of malignant melanoma. ^

Mechanism of molecular regulation of nuclear -encoded mitochondrial gene expression by electrical stimulation in the cultured neonatal rat cardiac myocytes

To answer the question whether increased energy demand resulting from myocyte hypertrophy and enhanced $\beta$-myosin heavy chain mRNA, contractile protein synthesis and assembly leads to mitochondrial proliferation and differentiation, we set up an electrical stimulation model of cultured neonatal rat cardiac myocytes. We describe, as a result of increased contractile activity, increased mitochondrial profiles, cytochrome oxidase mRNA, and activity, as well as a switch in mitochondrial carnitine palmitoyltransferase-I (CPT-I) from the liver to muscle isoform. We investigate physiological pathways that lead to accumulation of gene transcripts for nuclear encoded mitochondrial proteins in the heart. Cardiomyocytes were stimulated for varying times up to 72 hr in serum-free culture. The mRNA contents for genes associated with transcriptional activation (c-fos, c-jun, junB, nuclear respiratory factor 1 (Nrf-1)), mitochondrial proliferation (cytochrome c (Cyt c), cytochrome oxidase), and mitochondrial differentiation (carnitine palmitonyltransferase I (CPT-I) isoforms) were measured. The results establish a temporal pattern of mRNA induction beginning with c-fos (0.25-3 hr) and followed by c-jun (0.5-3 hr), junB (0.5-6 hr), NRF-1 (1-12 hr), Cyt c (12-72 hr), cytochrome c oxidase (12-72 hr). Induction of the latter was accompanied by a marked decrease in the liver-specific CPT-I mRNA. Electrical stimulation increased c-fos, $\beta$-myosin heavy chain, and Cyt c promoter activities. These increases coincided with a rise in their respective endogenous gene transcripts. NRF-1, cAMP response element (CRE), and Sp-1 site mutations within the Cyt c promoter reduced luciferase expression in both stimulated and nonstimulated myocytes. Mutations in the Nrf-1 and CRE sites inhibited the induction by electrical stimulation or by transfection of c-jun into non-paced cardiac myocytes whereas mutation of the Sp-1 site maintained or increased the fold induction. This is consistent with the appearance of NRF-1 and fos/jun mRNAs prior to that of Cyt c. Overexpression of c-jun by transfection also activates the Nrf-1 and Cyt c mRNA sequentially. Electrical stimulation of cardiac myocytes activates the c-Jun-N-terminal kinase so that the fold-activation of the cyt c promoter is increased by pacing when either c-jun or c-fos/c-jun are cotransfected. We have identified physical association of Nrf-1 protein with the Nrf-1 enhancer element and of c-Jun with the CRE binding sites on the Cyt c promoter. This is the first demonstration that induction of Nrf-1 and c-Jun by pacing of cardiac myocytes directly mediates Cyt c gene expression and mitochondrial proliferation in response to hypertrophic stimuli in the heart.^ Subsequent to gene activation pathways that lead to mitochondrial proliferation, we observed an isoform switch in CPT-I from the liver to muscle mRNA. We have found that the half-life for the muscle CPT-I is not affected by electrical stimulation, but electrical decrease the T1/2 in the liver CPT-I by greater than 50%. This suggests that the liver CPT-I switch to muscle isoform is due to (1) a decrease in T1/2 of liver CPT-I and (2) activation of muscle CPT-Itranscripts by electrical stimulation. (Abstract shortened by UMI.) ^

ATXA -dependent gene expression in Bacillus anthracis

The Bacillus anthracis toxin genes, cya, lef , and pag, can be viewed as a regulon, in which transcription of all three genes is activated in trans by the same regulatory gene, atxA, in response to the same signal, CO2. I determined that several phenotypes are associated with the atxA gene. In addition to being toxin-deficient, an atxA -null mutant grows poorly on minimal media and sporulates early compared to the parent strain. Furthermore, an atxA-null mutant has an altered 2-D gel protein profile. I used a genetic approach to find additional atxA-regulated genes. Random transcriptional lacZ fusions were generated in B. anthracis using transposon Tn 917-LTV3. Transposon-insertion libraries were screened for mutants expressing increased β-galactosidase activity in 5% CO2. Introduction of an atxA-null mutation in these mutants revealed that 79% of the CO2-regulated fusions were also atxA-dependent. DNA sequence analysis of transposon insertion sites in mutants carrying CO 2/atxA-regulated fusions revealed ten mutants harboring transposon insertions in loci distinct from the toxin genes. The majority of the tcr (toxin co-regulated) loci mapped within the pXO1 pathogenicity island. These results indicate a clear association of atxA with CO2-enhanced gene expression in B. anthracis and provide evidence that atxA regulates genes other than the structural genes for the anthrax toxin proteins. ^ Characterization of one tcr locus revealed a new regulatory gene, pagR. The pagR gene (300 nt) is located downstream of pag. pagR is cotranscribed with pag and is responsible for autogenous control of the operon. pagR also represses expression of cya and lef. Repression of toxin gene expression by pagR may be mediated by atxA. The steady state level of atxA mRNA is increased in a pagR mutant. Recombinant PagR protein purified from Escherichia coli did not specifically bind the promoter regions of pagA or atxA. An unidentified factor in B. anthracis crude extracts, however, was able to bind the atxA promoter in the absence of PagR or AtxA. These investigations increase our knowledge of virulence regulation in B. anthracis and ultimately will lead to a better understanding of anthrax disease. ^

Regulation of tyrosine hydroxylase gene expression by mRNA stabilization

Tyrosine hydroxylase (TH) expression increases in adrenal chromaffin cells treated with the nicotinic agonist, dimethylphenylpiperazinium (DMPP; 1 μM). We are using this response as a model of the changes in TH level that occur during increased cholinergic neural activity. Here we report a 4-fold increase in TH mRNA half-life in DMPP-treated chromaffin cells that is apparent when using a pulse-chase analysis to measure TH mRNA half-life. No increase is apparent using actinomycin D to measure half-life, indicating a requirement for ongoing transcription. Characterization of protein binding to the TH 3′UTR using RNA electro-mobility shift assays show the presence of two complexes both of which are increased by DMPP-treatment. The faster migrating complex (FMC) increases 2.5-fold and the slower migrating complex (SMC) increases 1.5-fold. Separation of UV crosslinked RNA-protein complexes on SDS polyacrylamide gels shows FMC to contain a single protein whereas SMC contains two proteins. Northwesterns yielded similar results. Transfection studies reveal an increase in expression of the full-length TH transcript due to DMPP-treatment similar to that of endogenous TH mRNA. This finding suggests the increased expression is due primarily to mRNA stabilization. Transfection of luciferase reporter constructs containing regions of the TH 3′UTR reveal only the full-length 3′UTR influenced the expression level of reporter transcripts. ^