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

Systems Biology Approaches to Probe Gene Regulation in Bacteria

Mechanisms that allow pathogens to colonize the host are not the product of isolated genes, but instead emerge from the concerted operation of regulatory networks. Therefore, identifying components and the systemic behavior of networks is necessary to a better understanding of gene regulation and pathogenesis. To this end, I have developed systems biology approaches to study transcriptional and post-transcriptional gene regulation in bacteria, with an emphasis in the human pathogen Mycobacterium tuberculosis (Mtb). First, I developed a network response method to identify parts of the Mtb global transcriptional regulatory network utilized by the pathogen to counteract phagosomal stresses and survive within resting macrophages. As a result, the method unveiled transcriptional regulators and associated regulons utilized by Mtb to establish a successful infection of macrophages throughout the first 14 days of infection. Additionally, this network-based analysis identified the production of Fe-S proteins coupled to lipid metabolism through the alkane hydroxylase complex as a possible strategy employed by Mtb to survive in the host. Second, I developed a network inference method to infer the small non-coding RNA (sRNA) regulatory network in Mtb. The method identifies sRNA-mRNA interactions by integrating a priori knowledge of possible binding sites with structure-driven identification of binding sites. The reconstructed network was useful to predict functional roles for the multitude of sRNAs recently discovered in the pathogen, being that several sRNAs were postulated to be involved in virulence-related processes. Finally, I applied a combined experimental and computational approach to study post-transcriptional repression mediated by small non-coding RNAs in bacteria. Specifically, a probabilistic ranking methodology termed rank-conciliation was developed to infer sRNA-mRNA interactions based on multiple types of data. The method was shown to improve target prediction in Escherichia coli, and therefore is useful to prioritize candidate targets for experimental validation.

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.

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

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

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

Gene transfer by viral vectors

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

Cloning, characterization and regulation of the soluble guanylyl cyclase alpha1 and beta1 genes

Cell signaling by nitric oxide (NO) through soluble guanylyl cyclase (sGC) and cGMP production regulates physiological responses such as smooth muscle relaxation, neurotransmission, and cell growth and differentiation. Although the NO receptor, sGC, has been studied extensively at the protein level, information on regulation of the sGC genes remains elusive. In order to understand the molecular mechanisms involved at the level of gene expression, cDNA and genomic fragments of the murine sGCα1 subunit gene were obtained through library screenings. Using the acquired clones, the sGCα 1 gene structure was determined following primer extension, 3 ′RACE and intron/exon boundary analyses. The basal activity of several 5′-flanking regions (putative promoter regions) for both the α1 and β1 sGC subunits were determined following their transfection into mouse N1E-115 neuroblastoma and rat RENE1Δ14 uterine epithelial cells using a luciferase reporter plasmid. Using the sGC sequences, real-time RT-PCR assays were designed to measure mRNA levels of the sGC α1 and β1 genes in rat, mouse and human. Subsequent studies found that uterine sGC mRNA and protein levels decreased rapidly in response to 17β-estradiol (estrogen) in an in vivo rat model. As early as 1 hour following treatment, mRNA levels of both sGC mRNAs decreased, and reached their lowest level of expression after 3 hours. This in vivo response was completely blocked by the pure estrogen receptor antagonist, ICI 182,780, was not seen in several other tissues examined, did not occur in response to other steroid hormones, and was due to a post-transcriptional mechanism. Additional studies ex vivo and in various cell culture models suggested that the estrogen-mediated decreased sGC mRNA expression did not require signals from other tissues, but may require cell communication or paracrine factors between different cell types within the uterus. Using chemical inhibitors and molecular targeting in other related studies, it was revealed that c-Jun-N-terminal kinase (JNK) signaling was responsible for decreased sGC mRNA expression in rat PC12 and RFL-6 cells, two models previously determined to exhibit rapid decreased sGC mRNA expression in response to different stimuli. To further investigate the post-transcriptional gene regulation, the full length sGCα1 3′-untranslated region (3′UTR) was cloned from rat uterine tissue and ligated downstream of the rabbit β-globin gene and expressed as a chimeric mRNA in the rat PC12 and RFL-6 cell models. Expression studies with the chimeric mRNA showed that the sGCα 1 3′UTR was not sufficient to mediate the post-transcriptional regulation of its mRNA by JNK or cAMP signaling in PC12 and RFL-6 cells. This study has provided numerous valuable tools for future studies involving the molecular regulation of the sGC genes. Importantly, the present results identified a novel paradigm and a previously unknown signaling pathway for sGC mRNA regulation that could potentially be exploited to treat diseases such as uterine cancers, neuronal disorders, hypertension or various inflammatory conditions. ^

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

FUNCTIONS OF DEADENYLATION FACTORS IN MRNA DECAY AND MRNA PROCESSING BODY FORMATION

Most newly synthesized messenger RNAs possess a 5’ cap and a 3’ poly(A) tail. The process of poly(A) tail shortening, also termed deadenylation, is important for post-transcriptional gene regulation, because deadenylation not only leads to mRNA translational inhibition but also is the first step of major mRNA degradation. Translationally inhibited mRNAs can be stored and/or degraded in dynamic cytoplasmic foci termed mRNA processing bodies, or P bodies, which are conserved in eukaryotes. To shed new light on the mechanisms of P body formation and P body functions, I focused on the link between deadenylation factors and P bodies. I found that the two major deadenylation complexes, Pan3-Pan2 and Ccr4-Caf1, can both be enriched in P bodies. The deadenylase activity of the Ccr4-Caf1 complex is prerequisite for P body formation. Pan3, but not the deadenylase Pan2, is essential for P body formation. While the C-terminal domain of Pan3 is important for interaction with Pan2, Pan3 N-terminal domain is important for Pan3 to form cytoplasmic foci colocalizing with P bodies and to promote mRNA decay. Interestingly, Pan3 N-terminal domain may be phosphorylated to regulate Pan3 localization and functions. Aside from the functions of the two deadenylation complexes in P bodies, I also studied all reported human P body proteins as a whole using bioinformatics. This effort not only has generated a comprehensive picture of the functions of and interactions among human P body proteins, but also has predicted proteins that may regulate P body formation and/or functions. In summary, my study has established a direct link between mRNA deadenylation and P body formation and has also led to new hypotheses to guide future research on how P body dynamics are controlled.

Mammalian miRNA RISC recruits CAF1 and PABP to affect PABP-dependent deadenylation.

MicroRNAs (miRNAs) inhibit mRNA expression in general by base pairing to the 3'UTR of target mRNAs and consequently inhibiting translation and/or initiating poly(A) tail deadenylation and mRNA destabilization. Here we examine the mechanism and kinetics of miRNA-mediated deadenylation in mouse Krebs-2 ascites extract. We demonstrate that miRNA-mediated mRNA deadenylation occurs subsequent to initial translational inhibition, indicating a two-step mechanism of miRNA action, which serves to consolidate repression. We show that a let-7 miRNA-loaded RNA-induced silencing complex (miRISC) interacts with the poly(A)-binding protein (PABP) and the CAF1 and CCR4 deadenylases. In addition, we demonstrate that miRNA-mediated deadenylation is dependent upon CAF1 activity and PABP, which serves as a bona fide miRNA coactivator. Importantly, we present evidence that GW182, a core component of the miRISC, directly interacts with PABP via its C-terminal region and that this interaction is required for miRNA-mediated deadenylation.

TRANSCRIPTIONAL AND POST-TRANSLATIONAL MECHANISMS CONTRIBUTE TO MAINTENANCE OF REST IN NEURAL TUMORS

The RE-1 silencing transcription factor (REST) is an important regulator of normal nervous system development. It negatively regulates neuronal lineage specification in neural progenitors by binding to its consensus RE-1 element(s) located in the regulatory region of its target neuronal differentiation genes. The developmentally coordinated down-regulation of REST mRNA and protein in neural progenitors triggers terminal neurogenesis. REST is overexpressed in pediatric neural tumors such as medulloblastoma and neuroblastoma and is associated with poor neuronal differentiation. High REST protein correlate with poor prognosis for patients with medulloblastoma, however similar studies have not been done with neuroblastoma patients. Mechanism(s) underlying elevated REST levels medulloblastoma and neuroblastoma are unclear, and is the focus of this thesis project. We discovered that transcriptional and post-translational mechanisms govern REST mis-regulation in medulloblastoma and neuroblastoma. In medulloblastoma, REST transcript is aberrantly elevated in a subset of patient samples. Using loss of function and gain of function experiments, we provide evidence that the Hairy Enhancer of Split (HES1) protein represses REST transcription in medulloblastoma cell lines, modulates the expression of neuronal differentiation genes, and alters the survival potential of these cells in vitro. We also show that REST directly represses its own expression in an auto-regulatory feedback loop. Interestingly, our studies identified a novel interaction between REST and HES1. We also observed their co-occupancy at the RE-1 sites, thereby suggesting potential for co-regulation of REST expression. Our pharmacological studies in neuroblastoma using retinoic acid revealed that REST levels are controlled by transcriptional and post-transcriptional mechanisms. Post-transcriptional mechanisms are mediated by modulation of E3 ligase or REST, SCFβ-TRCP, and contribute to resistance of some cells to retinoic acid treatment.

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.

The transcriptional regulation of the {\it neu\/} gene: Examples of activation, repression, and cell-specific expression

The neu oncogene encodes a growth factor receptor-like protein, p185, with an intrinsic tyrosine kinase activity. A single point mutation, an A to T transversion resulting in an amino acid substitution from valine to glutamic acid, in the transmembrane domain of the rat neu gene was found to be responsible for the transforming and tumorigenic phenotype of the cells that carry it. In contrast, the human proto-neu oncogene is frequently amplified in tumors and cell lines derived from tumors and the human neu gene overexpression/amplification in breast and ovarian cancers is known to correlate with poor patient prognosis. Examples of the human neu gene overexpression in the absence of gene amplification have been observed, which may suggest the significant role of the transcriptional and/or post-transcriptional control of the neu gene in the oncogenic process. However, little is known about the transcriptional mechanisms which regulate the neu gene expression. In this study, three examples are presented to demonstrate the positive and negative control of the neu gene expression.^ First, by using band shift assays and methylation interference analyses, I have identified a specific protein-binding sequence, AAGATAAAACC ($-$466 to $-$456), that binds a specific trans-acting factor termed RVF (for EcoRV factor on the neu promoter). The RVF-binding site is required for maximum transcriptional activity of the rat neu promoter. This same sequence is also found in the corresponding regions of both human and mouse neu promoters. Furthermore, this sequence can enhance the CAT activity driven by a minimum promoter of the thymidine kinase gene in an orientation-independent manner, and thus it behaves as an enhancer. In addition, Southwestern (DNA-protein) blot analysis using the RVF-binding site as a probe points to a 60-kDa polypeptide as a potential candidate for RVF.^ Second, it has been reported that the E3 region of adenovirus 5 induces down-regulation of epidermal growth factor (EGF) receptor through endocytosis. I found that the human neu gene product, p185, (an EGF receptor-related protein) is also down-regulated by adenovirus 5, but via a different mechanism. I demonstrate that the adenovirus E1a gene is responsible for the repression of the human neu gene at the transcriptional level.^ Third, a differential expression of the neu gene has been found in two cell model systems: between the mouse fibroblast Swiss-Webster 3T3 (SW3T3) and its variant NR-6 cells; and between the mouse liver tumor cell line, Hep1-a, and the mouse pancreas tumor cell line, 266-6. Both NR-6 and 266-6 cell lines are not able to express the neu gene product, p185. I demonstrate that, in both cases, the transcriptional repression of the neu gene may account for the lack of the p185 expression in these two cell lines. ^

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