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.

Regulation of survivin gene expression in the human endometrium and endometrial cancer

In the United States, endometrial cancer is the leading cancer of the female reproductive tract. There are 40,100 new cases and 7,470 deaths from endometrial cancer estimated for 2008 (47). The average five year survival rate for endometrial cancer is 84% however, this figure is substantially lower in patients diagnosed with late stage, advanced disease and much higher for patients diagnosed in early stage disease (47). Endometrial cancer (EC) has been associated with several risk factors including obesity, diabetes, hypertension, previously documented occurrence of hereditary non-polyposis colorectal cancer (HNPCC), and heightened exposure to estrogen (25). As of yet, there has not been a dependable molecular predictor of endometrial cancer occurrence in women with these predisposing factors. The goal of our lab is to identify genes that are aberrantly expressed in EC and may serve as molecular biomarkers of EC progression. One candidate protein that we are exploring as a biomarker of EC progression is the cell survival protein survivin.

XENOESTROGEN-SPECIFIC MECHANISMS OF DEVELOPMENTAL REPROGRAMMING CORRELATE WITH GENE EXPRESSION AND TUMOR DEVELOPMENT

Environmental exposures during sensitive windows of development can reprogram normal physiological responses and alter disease susceptibility later in life in a process known as developmental reprogramming. We have shown that neonatal exposure to the xenoestrogen diethylstilbestrol (DES) can developmentally reprogram the reproductive tract in genetically susceptible Eker rats giving rise to complete penetrance of uterine leiomyoma. Based on this, we hypothesized that xenoestrogens, including genistein (GEN) and bisphenol A (BPA), reprogram estrogen-responsive gene expression in the myometrium and promote the development of uterine leiomyoma. We proposed the mechanism that is responsible for the developmental reprogramming of gene expression was through estrogen (E2)/ xenoestrogen inducedrapid ER signaling, which modifies the histone methyltransferase Enhancer of Zeste homolog 2 (EZH2) via activation of the PI3K/AKT pathway. We further hypothesized that there is a xenostrogen-specific effect on this pathway altering patterns of histone modification, DNA methylation and gene expression. In addition to our novel finding that E2/DES-induced phosphorylation of EZH2 by AKT reduces the levels of H3K27me3 in vitro and in vivo, this work demonstrates in vivo that a brief neonatal exposure to GEN, in contrast to BPA, activates the PI3K/AKT pathway to regulate EZH2 and decreases H3K27me3 levels in the neonatal uterus. Given that H3K27me3 is a repressive mark that has been shown to result in DNA methylation and gene silencing we investigated the methylation of developmentally reprogrammed genes. In support of this evidence, we show that neonatal DES exposure in comparison to VEH, leads to hypomethylation of the promoter of a developmentally reprogrammed gene, Gria2, that become hyper-responsive to estrogen in the adult myometrium indicating vi that DES exposure alter gene expression via chromatin remodeling and loss of DNA methylation. In the adult uterus, GEN and BPA exposure developmentally reprogrammed expression of estrogen-responsive genes in a manner opposite of one another, correlating with our previous data. Furthermore, the ability of GEN and BPA to developmental reprogram gene expression correlated with tumor incidence and multiplicity. These data show that xenoestrogens have unique effects on the activation of non-genomic signaling in the developing uterus that promotes epigenetic and genetic alterations, which are predictive of developmental reprogramming and correlate with their ability to modulate hormone-dependent tumor development.

SURVIVAL PREDICTION FOR BRAIN TUMOR PATIENTS USING GENE EXPRESSION DATA

Brain tumor is one of the most aggressive types of cancer in humans, with an estimated median survival time of 12 months and only 4% of the patients surviving more than 5 years after disease diagnosis. Until recently, brain tumor prognosis has been based only on clinical information such as tumor grade and patient age, but there are reports indicating that molecular profiling of gliomas can reveal subgroups of patients with distinct survival rates. We hypothesize that coupling molecular profiling of brain tumors with clinical information might improve predictions of patient survival time and, consequently, better guide future treatment decisions. In order to evaluate this hypothesis, the general goal of this research is to build models for survival prediction of glioma patients using DNA molecular profiles (U133 Affymetrix gene expression microarrays) along with clinical information. First, a predictive Random Forest model is built for binary outcomes (i.e. short vs. long-term survival) and a small subset of genes whose expression values can be used to predict survival time is selected. Following, a new statistical methodology is developed for predicting time-to-death outcomes using Bayesian ensemble trees. Due to a large heterogeneity observed within prognostic classes obtained by the Random Forest model, prediction can be improved by relating time-to-death with gene expression profile directly. We propose a Bayesian ensemble model for survival prediction which is appropriate for high-dimensional data such as gene expression data. Our approach is based on the ensemble "sum-of-trees" model which is flexible to incorporate additive and interaction effects between genes. We specify a fully Bayesian hierarchical approach and illustrate our methodology for the CPH, Weibull, and AFT survival models. We overcome the lack of conjugacy using a latent variable formulation to model the covariate effects which decreases computation time for model fitting. Also, our proposed models provides a model-free way to select important predictive prognostic markers based on controlling false discovery rates. We compare the performance of our methods with baseline reference survival methods and apply our methodology to an unpublished data set of brain tumor survival times and gene expression data, selecting genes potentially related to the development of the disease under study. A closing discussion compares results obtained by Random Forest and Bayesian ensemble methods under the biological/clinical perspectives and highlights the statistical advantages and disadvantages of the new methodology in the context of DNA microarray data analysis.

The recombinant adeno-associated virus vector (rAAV2)-mediated apolipoprotein B mRNA-specific hammerhead ribozyme: a self-complementary AAV2 vector improves the gene expression.

BACKGROUND: In humans, overproduction of apolipoprotein B (apoB) is positively associated with premature coronary artery diseases. To reduce the levels of apoB mRNA, we have designed an apoB mRNA-specific hammerhead ribozyme targeted at nucleotide sequences GUA6679 (RB15) mediated by adenovirus, which efficiently cleaves and decreases apoB mRNA by 80% in mouse liver and attenuates the hyperlipidemic condition. In the current study, we used an adeno-associated virus vector, serotype 2 (AAV2) and a self-complementary AAV2 vector (scAAV2) to demonstrate the effect of long-term tissue-specific gene expression of RB15 on the regulation apoB mRNA in vivo. METHODS: We constructed a hammerhead ribozyme RB15 driven by a liver-specific transthyretin (TTR) promoter using an AAV2 vector (rAAV2-TTR-RB15). HepG2 cells and hyperlipidemic mice deficient in both the low density lipoprotein receptor and the apoB mRNA editing enzyme genes (LDLR-/-Apobec1-/-; LDb) were transduced with rAAV2-TTR-RB15 and a control vector rAAV-TTR-RB15-mutant (inactive ribozyme). The effects of ribozyme RB15 on apoB metabolism and atherosclerosis development were determined in LDb mice at 5-month after transduction. A self-complementary AAV2 vector expressing ribozyme RB15 (scAAV2-TTR-RB15) was also engineered and used to transduce HepG2 cells. Studies were designed to compare the gene expression efficiency between rAAV2-TTR-RB15 and scAAV2-TTR-RB15. RESULTS: The effect of ribozyme RB15 RNA on reducing apoB mRNA levels in HepG2 cells was observed only on day-7 after rAAV2-TTR-RB15 transduction. And, at 5-month after rAAV2-TTR-RB15 treatment, the apoB mRNA levels in LDb mice were significantly decreased by 43%, compared to LDb mice treated with control vector rAAV2-TTR-RB15-mutant. Moreover, both the rAAV2-TTR-RB15 viral DNA and ribozyme RB15 RNA were still detectable in mice livers at 5-month after treatment. However, this rAAV2-TTR-RB15 vector mediated a prolonged but low level of ribozyme RB15 gene expression in the mice livers, which did not produce the therapeutic effects on alteration the lipid levels or the inhibition of atherosclerosis development. In contrast, the ribozyme RB15 RNA mediated by scAAV2-TTR-RB15 vector was expressed immediately at day-1 after transduction in HepG2 cells. The apoB mRNA levels were decreased 47% (p = 0.001), compared to the control vector scAAV2-TTR-RB15-mutant. CONCLUSION: This study provided evidence that the rAAV2 single-strand vector mediated a prolonged but not efficient transduction in mouse liver. However, the scAAV2 double-strand vector mediated a rapid and efficient gene expression in liver cells. This strategy using scAAV2 vectors represents a better approach to express small molecules such as ribozyme.

Regulation of mucin gene expression by CREB via a nonclassical retinoic acid signaling pathway.

Vitamin A and its metabolite retinoic acid (RA) are essential elements for normal lung development and the differentiation of lung epithelial cells. We previously showed that RA rapidly activated cyclic AMP response element-binding protein (CREB) in a nonclassical manner in normal human tracheobronchial epithelial (NHTBE) cells. In the present study, we further demonstrated that this nonclassical signaling of RA on the activation of CREB plays a critical role in regulating the expression of airway epithelial cell differentiation markers, the MUC2, MUC5AC, and MUC5B genes. We found that RA rapidly activates the protein kinase Calpha isozyme and transmits the activation signal to CREB via the Raf/MEK/extracellular signal-regulated kinase/p90 ribosomal S6 kinase (RSK) pathway. Activated RSK translocated from the cytoplasm to the nucleus, where it phosphorylates CREB. Activated CREB then binds to a cis-acting replication element motif on the promoter (at nucleotides [nt] -878 to -871) of the MUC5AC gene. The depletion of CREB using small interfering RNA abolished not only the RA-induced MUC5AC but also RA-induced MUC2 and MUC5B. Taken together, our findings demonstrate that CREB activation via this nonclassical RA signaling pathway may play an important role in regulating the expression of mucin genes and mediating the early biological effects of RA during normal mucous differentiation in NHTBE cells.

Enterococcus faecalis rnjB is required for pilin gene expression and biofilm formation.

Pili in Gram-positive bacteria play a major role in the colonization of host tissue and in the development of biofilms. They are promising candidates for vaccines or drug targets since they are highly immunogenic and share common structural and functional features among various Gram-positive pathogens. Numerous publications have helped build a detailed understanding of pilus surface assembly, yet regulation of pilin gene expression has not been well defined. Utilizing a monoclonal antibody developed against the Enterococcus faecalis major pilus protein EbpC, we identified mutants from a transposon (Tn) insertion library which lack surface-exposed Ebp pili. In addition to insertions in the ebp regulon, an insertion in ef1184 (dapA) significantly reduced levels of EbpC. Analysis of in-frame dapA deletion mutants and mutants with the downstream gene rnjB deleted further demonstrated that rnjB was responsible for the deficiency of EbpC. Sequence analysis revealed that rnjB encodes a putative RNase J2. Subsequent quantitative real-time PCR (qRT-PCR) and Northern blotting demonstrated that the ebpABC mRNA transcript level was significantly decreased in the rnjB deletion mutant. In addition, using a reporter gene assay, we confirmed that rnjB affects the expression of the ebpABC operon. Functionally, the rnjB deletion mutant was attenuated in its ability to produce biofilm, similar to that of an ebpABC deletion mutant which lacks Ebp pili. Together, these results demonstrate the involvement of rnjB in E. faecalis pilin gene expression and provide insight into a novel mechanism of regulation of pilus production in Gram-positive pathogens.

The molecular mechanism of retinoic acid action: Regulation of tissue transglutaminase gene expression

Retinoic acid is a small lipophilic molecule that exerts profound effects on the growth and differentiation of both normal and transformed cells. It is also a natural morphogen that is critical in the development of embryonic structures. The molecular effects of retinoic acid involve alterations in the expression of several proteins and these changes are presumably mediated in part by alterations in gene expression. For instance, retinoic acid causes a rapid induction of tissue transglutaminase, an enzyme involved in protein cross-linking. The molecular mechanisms responsible for the effects of retinoic acid on gene expression have not been characterized. To approach this question, I have isolated and characterized tissue transglutaminase of cDNA clones. The deduced amino acid sequences of tissue transglutaminase and of factor XIIIa showed a relatively high degree of homology in their putative calcium binding domains.^ To explore the mechanism of induction of this enzyme, both primary (macrophages) and cultured cells (Swiss 3T3-C2 and CHO fibroblasts) were used. I found that retinoic acid is a general inducer of tissue transglutaminase mRNA in these cells. In murine peritoneal macrophages retinoic acid causes a rapid accumulation of this mRNA and this effect is independent of concurrent protein synthesis. The retinoic acid effect is not mediated by a post-transcriptional increase in the stability of the tissue transglutaminase mRNA, but appears to involve an increase in the transcription rate of the tissue transglutaminase gene. This provides the first example of regulation by retinoic acid of a specific gene, supporting the hypothesis that these molecules act by directly regulating the transcriptional activity of specific genes. A molecular model for the effects of retinoic acid on the expression of genes linked to cellular proliferation and differentiation is proposed. ^

Structure-function relation in retinoid regulated gene expression

Retinoids are known to inhibit proliferation of and induce terminal differentiation of many normal and transformed cells. It has been postulated that retinoids exert their effect by altering gene expression. HL-60 cells and macrophages both respond to retinoic acid action by the rapid induction of the enzyme tissue transglutaminase. The induction has been shown to be due to increased transcription of the transglutaminase gene. The first part of the dissertation studied the structure-function relationship of retinoid-regulated transglutaminase induction, differentiation and proliferation in HL-60 cells using retinoid analogs. The results indicated strict structural constraints and a strong structure-function correlation between transglutaminase induction and differentiation; those retinoids that induced transglutaminase also induced differentiation, those analogs that did not induce transglutaminase could not induce differentiation. The ability of the retinoids to induce transglutaminase in HL-60 cells was paralleled in macrophages. However, the antiproliferative effect of the retinoids displayed less stringent structural constraints than their differentiation- and transglutaminase-inducing properties. Specifically all the retinoids were able to inhibit proliferation to varying extents. It is concluded that the induction of transglutaminase and of differentiation by retinoids is mediated by receptors. While receptor mediation cannot be entirely ruled out, with the current data no definitive statement can be made about the antiproliferative activity of retinoids. Also, the concordance in the ability of the retinoids to induce transglutaminase and the ability to induce differentiation of HL-60 cells suggests that the former is an early response of the cells to retinoids and differentiation a later consequence on the same pathway. Using the induction of transglutaminase as an index of the direct, or primary, effect of retinoids on gene expression, the second part of the dissertation investigates, by 2D gel electrophoresis, the alteration in the rates of synthesis of other proteins in macrophages and HL-60 cells in response to short incubations with retinoic acid. Any changes in parallel with transglutaminase were taken to indicate proteins directly under the control of retinoic acid. It is concluded that retinoic acid regulates the expression of a circumscribed set of genes in a cell-specific manner. The results support the hypothesis that retinoids exert their multiple effects on myeloid cells, in part, by receptor-mediated alternations in gene expression. ^

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

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

H2B histone gene expression during {\it Xenopus laevis\/} development

Histone gene expression is replication-independent during oogenesis and early embryogenesis in amphibians; however, it becomes replication-dependent during later embryogenesis and remains replication-dependent through adulthood. In order to understand the mechanism for this switch in transcriptional regulation of histone gene expression during amphibian development, linker-scanning mutations were made in a Xenopus laevis H2B histone gene promoter by oligonucleotide site-directed mutagenesis and assayed by microinjection into oocytes and embryos. The Xenopus H2B gene has a relatively simple promoter containing several transcriptional regulatory elements, including TFIID, CCAAT, and ATF motifs, required for maximal transcription in both oocytes and embryos. Factors binding to the CCAAT and ATF motifs are present in oocytes and embryos and increase slightly in abundance during early development. A sequence (CTTTACAT) in the frog H2B promoter resembling the conserved octamer motif (ATTTGCAT), the target for cell-cycle regulation of a human H2B gene, is additionally required for maximal H2B transcription in frog embryos. Oocytes and embryos contain multiple octamer-binding proteins that are expressed in a sequential manner during early development. Sequences encoding three novel octamer-binding proteins were isolated from Xenopus cDNA libraries by virtue of their similarity with the DNA binding (POU) domain of the ubiquitously expressed transcription factor Oct-1. The protein encoded by one of these genes, termed Oct-60, was localized mainly in the cytoplasm of oocytes and was also present in early embryos until the gastrula stage of development. Proteins encoded by the other two genes, Oct-25 and Oct-91, were present in embryos after the mid-blastula stage of development and decreased by early neurula stage. The activity of the Xenopus H2B octamer motif in embryos is not specifically associated with increased binding by Oct-1 or the appearance of novel octamer-binding proteins but requires the presence of an intact CCAAT motif. We found that synergistic interactions among promoter elements are important for full H2B promoter activity. The results suggest that transcription of the Xenopus H2B gene is replication-dependent when it is activated at the mid-blastula stage of development and that replication-dependent H2B transcription is mediated by Oct-1. ^

Gene expression in sea urchin development: Study of {\it LpS1\/} gene regulation and cloning and characterization of the {\it SpKrox}-{\it 1\/} gene

Cell differentiation are associated with activation of cell lineage-specific genes. The $LpS{\it 1}\beta$ gene of Lytechinus pictus is activated at the late cleavage stage. $LpS{\it 1}\beta$ transcripts accumulate exclusively in aboral ectoderm lineages. Previous studies demonstrated two G-string DNA-elements, proximal and distal G-strings, which bind to an ectoderm-enriched nuclear factor. In order to define the cis-elements which control positive expression of the $LpS{\it 1}\beta$ gene, the regulatory region from $-$108 to +17 bp of the $LpS{\it 1}\beta$ gene promoter was characterized. The ectoderm G-string factor binds to a G/C-rich region larger than the G-string itself and the binding of the G-string factor requires sequences immediately downstream from the G-string. These downstream sequences are essential for full promoter activity. In addition, only 108 bp of $LpS{\it 1}\beta\ 5\sp\prime$ flanking DNA drives $LpS{\it 1}\beta$ gene expression in aboral ectoderm/mesenchyme cells. Therefore, for positive control of $LpS{\it 1}\beta$ gene expression, two regions of 5$\sp\prime$ flanking DNA are required: region I from base pairs $-$762 to $-$511, and region II, which includes the G/C-rich element, from base pairs $-$108 to $-$61. A mesenchyme cell repressor element is located within region I.^ DNA-binding proteins play key roles in determination of cell differentiation. The zinc finger domain is a DNA-binding domain present in many transcription factors. Based on homologies in zinc fingers, a zinc finger-encoding gene, SpKrox-1, was cloned from S. purpuratus. The putative SpKrox-1 protein has all structural characteristics of a transcription factor: four zinc fingers for DNA binding; acidic domain for transactivation; basic domain for nuclear targeting; and leucine zipper for dimerization. SpKrox-1 RNA transcripts showed a transient expression pattern which correlates largely with early embryonic development. The spatial expression of SpKrox-1 mRNA was distributed throughout the gastrula and larva ectodermal wall. However, SpKrox-1 was not expressed in pigment cells. The SpKrox-1 gene is thus a marker of a subset of SMCs or ectoderm cells. The structural features, and the transient temporal and restricted spatial expression patterns suggest that SpKrox-1 plays a role in a specific developmental event. ^

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

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

A urokinase receptor promoter element (-152/-135) bound with an AP-2-alpha-related protein and SP1 is required for the induction of gene expression by PMA and SRC

The urokinase-type plasminogen activator receptor (u-PAR) promotes extracellular matrix degradation, invasion and metastasis. A first objective of this dissertation was to identify cis-elements and trans-acting factors activating u-PAR gene expression through a previously footprinted (–148/–124) promoter region. Mobility shifting experiments on nuclear extracts of a high u-PAR-expressing colon cancer cell line (RKO) indicated Sp1, Sp3 and a factor similar to, but distinct from, AP-2α bound to an oligonucleotide spanning –152/–135. Mutations preventing the binding of the AP-2α-related factor reduced u-PAR promoter activity. In RKO, the expression of a dominant negative AP-2 (AP-2αB) diminished u-PAR promoter activity, protein and u-PAR mediated laminin degradation. Conversely, u-PAR promoter activity in low u-PAR-expressing GEO cells was increased by AP-2αA expression. PMA treatment, which induces u-PAR expression, caused an increased amount of the AP-2α-related factor-containing complex in GEO, and mutations preventing AP-2α-like and Sp1/Sp3 binding reduced the u-PAR promoter stimulation by PMA. In resected colon cancers, u-PAR protein amounts were related to the amount of the AP-2α-related factor-containing complex. In conclusion, constitutive and PMA- inducible u-PAR gene expression and -proteolysis are mediated partly through transactivation via a promoter sequence (–152/435) bound with an AP-2α-related factor and Sp1/Sp3. ^ A second interest of this dissertation was to determine if a constitutively active Src regulates the transcription of the u-PAR gene, since c-src expression increases invasion in colon cancer. Increased u-PAR protein and laminin degradation paralleling elevated Src activity was evident in SW480 colon cancer cells stably expressing a constitutively active Src (Y- c-src527F). Nuclear run-on experiments indicated that this was due largely to transcriptional activation. While transient transfection of SW480 cells with Y-c-src527F induced a u-PAR-CAT-reporter, mutations preventing Sp1-binding to promoter region –152/435 abolished this induction. Mobility shift assays revealed increased Sp1 binding to region –152/135 with nuclear extracts of Src-transfected SW480 cells. Finally, the amounts of endogenous u-PAR in resected colon cancers significantly correlated with Src-activity. These data suggest that u-PAR gene expression and proteolysis are regulated by Src, this requiring the promoter region (–152/–135) bound with Sp1, thus, demonstrating for the first time that transcription factor Sp1 is a downstream effector of Src. ^