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.

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

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.

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

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

An extracytoplasmic function sigma factor operon regulates Myxococcus xanthus developmental gene expression

Myxococcus xanthus is a Gram-negative soil bacterium that undergoes multicellular development when high-density cells are starved on a solid surface. Expression of the 4445 gene, predicted to encode a periplasmic protein, commences 1.5 h after the initiation of development and requires starvation and high density conditions. Addition of crude or boiled supernatant from starving high-density cells restored 4445 expression to starving low-density cells. Addition of L-threonine or L-isoleucine to starving low-density cells also restored 4445 expression, indicating that the high-density signaling activity present in the supernatant might be composed of extracellular amino acids or small peptides. To investigate the circuitry integrating these starvation and high-density signals, the cis- and trans-acting elements controlling 4445 expression were identified. The 4445 transcription start site was determined by primer extension analysis to be 58 by upstream of the predicted translation start site. The promoter region contained a consensus sequence characteristic of e&barbelow;xtrac&barbelow;ytoplasmic f&barbelow;unction (ECF) sigma factor-dependent promoters, suggesting that 4445 expression might be regulated by an ECF sigma factor-dependent pathway, which are known to respond to envelope stresses. The small size of the minimum regulatory region, identified by 5′-end deletion analysis as being only 66 by upstream of the transcription start site, suggests that RNA polymerase could be the sole direct regulator of 4445 expression. To identify trans-acting negative regulators of 4445 expression, a strain containing a 4445-lacZ was mutagenized using the Himar1-tet transposon. The four transposon insertions characterized mapped to an operon encoding a putative ECF sigma factor, ecfA; an anti-sigma factor, reaA; and a negative regulator, reaB. The reaA and the reaB mutants expressed 4445 during growth and development at levels almost 100-fold higher than wild type, indicating that these genes encode negative regulators. The ecfA mutant expressed 4445-lacZ at basal levels, indicating that ecfA is a positive regulator. High Mg2+ concentrations over-stimulated this ecfA pathway possibly due to the depletion of exopolysaccharides and assembled type IV pili. These data indicate that the ecfA operon encodes a new regulatory stress pathway that integrates and transduces starvation and cell density cues during early development and is also responsive to cell-surface alterations.^

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

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 {\it neu\/} gene expression by the simian virus 40 large T antigen and tumor suppressors Rb and p53

The neu gene (also c-erbB-2 or HER2) encodes a 185 kilodalton protein that is frequently overexpressed in breast, ovarian and non-small cell lung cancers. Study of the regulation of neu indicates that neu gene expression can be modulated by c-myc or by the adenovirus 5 E1a gene product. This study demonstrates that the transforming protein, large T antigen, of the simian virus 40 represses neu promoter activity. Repression of neu by large T antigen is mediated through the region $-$172 to $-$79 (relative to first ATG) of the neu promoter--unlike through $-$312 to $-$172 for c-myc or E1a. This suggests a different pathway for repression of neu by large T antigen. The 10 amino acid region of large T required for binding the tumor suppressor, retinoblastoma gene product, Rb, is not necessary for repression of neu. Moreover, the tumor suppressors, Rb and p53 can independently inhibit neu promoter activity. Rb inhibits neu through a 10 base pair G-rich enhancer (GTG element) ($-$243 to $-$234) and also through regions close to transcription initiation sites ($-$172 to $-$79). Mutant Rb unable to complex large T is able to repress the region close to transcription initiation but not the GTG enhancer. Thus, Rb inhibits the two regulatory domains of the neu gene by different mechanisms. Both Rb and p53 can repress the transforming activity of activated neu in focus forming assays. These data provide evidence that tumor suppressors regulate expression of growth stimulatory genes such as neu. Therefore, one reason for the overexpression of neu that is frequently seen in breast cancer cells may be due to functional inactivation of Rb and p53 which is also a common occurrence in breast cancer cells. ^

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.

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.