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

GENE EXPRESSION IN CHRONIC MYELOGENOUS LEUKEMIA

The goal of the present work was to identify and characterize gene sequences that are preferentially expressed in CML in an effort to better understand the molecular basis of the disease. As high abundance mRNAs generally encode proteins that are phenotypically characteristic of cells, positive-negative screening of a CML cDNA library was used to identify cDNA clones containing sequences preferentially transcribed in CML. One cDNA sequence that fulfilled this criterion, C-A3, has been characterized in some detail. It represents a small mRNA ((TURN)496 nucleotides) that is highly abundant ((TURN)2% of the poly(A('+))RNA) in cells from the chronic phase of CML. In situ hybridization to whole cells indicates the principal leukocytes that express C-A3 sequences are eosinophils, basophils and immature myelocytes. Surprisingly, CML patients with high numbers of myeloblasts do not have an abundance of C-A3 transcripts, although transcript levels remain elevated in patients with lymphoblasts. In AML, high transcript levels are only found sporadically and occasionally different sized transcripts can be detected. Sequences from the 3' end of the C-A3 message are present in 2-5 copies per haploid genome. The 3' end of C-A3 localizes to bands 8q21.1 and 8q23 by in situ chromosomal hybridization. This is a region that is often involved in hematopoietic malignancies. Restriction digests of human genomic DNA show a correlation between the presence of a 2.3 kb Hind III fragment and certain types of leukemia. All of the leukemic DNAs tested had this fragment. In comparison, only one of five normal DNAs had a band this size. Analysis of the nucleotide sequence indicates that C-A3 probably encodes a small, hydrophobic peptide which may be part of a larger protein. ^

THE MECHANISM OF RETINOIC ACID REGULATION OF TISSUE TRANSGLUTAMINASE GENE EXPRESSION

Retinoic acid regulates cellular growth and differentiation by altering the expression of specific sets of genes, but the molecular mechanism by which this is achieved is unknown. We have used the rapid induction of a specific enzyme, tissue transglutaminase in mouse macrophages, human leukemia cells and a variety of other cell types to study the regulation of gene expression by retinoic acid. Soluble retinoic acid binding proteins, such as cellular Retinoic Acid Binding Protein (cRABP), have been proposed as specific mediators of retinoic acid regulation of gene expression. This thesis demonstrates the lack of cRABP in a number of cell lines which are sensitive to retinoic acid regulation of tissue transglutaminase expression. These cells are also devoid of other soluble retinoic acid binding activity. The level of retinoic acid binding activity that could have been detected (6 fmol) is far below that of most cells and tissues which are sensitive to the effects of retinoic acid on growth and differentiation. A mouse melanoma cell line, S91-C2, was found to contain an unusual retinoic acid binding protein which has a lower affinity for retinoic acid than mouse tissue cRABP and also behaves differently on gel filtration HPLC chromatography.^ The induction of tissue transglutaminase by retinoic acid in macrophages is specifically inhibited by pertussis toxin. Pertussis toxin ADP-riblosylates membrane GTP-binding proteins such as N(,i) and interferes with signalling from plasma membrane receptors to regulatory enzymes. Pertussis toxin inhibition of transglutaminase induction is due to inhibition of tissue transglutaminase mRNA accumulation and is paralleled by the ADP-ribosylation of a 41,000 dalton macrophage membrane protein. It is concluded that soluble retinoic acid binding proteins are not essential for retinoic acid induction of tissue transglutaminase and that a membrane GTP-binding protein is closely linked to the sensitive response of macrophages to retinoic acid. ^

New aspects of estrogen action in the endometrium: Reciprocal interplay with retinoic acid pathway and a novel estrogen-regulated gene

The uterine endometrium is a major target for the estrogen. However, the molecular basis of estrogen action in the endometrium is largely unknown. I have used two approaches to study the effects of estrogen on the endometrium. One approach involved the study of the interaction between estrogen and retinoic acid (RA) pathways in the endometrium. I have demonstrated that estrogen administration to rodents and estrogen replacement therapy (ERT) in postmenopausal women selectively induced the endometrial expression of retinaldehyde dehydrogenase II (RALDH2), a critical enzyme of RA biosynthesis. RALDH2 was expressed exclusively in the stromal cells, especially in the stroma adjacent to the luminal and glandular epithelia. The induction of RALDH2 by estrogen required estrogen receptor and occurred via a direct increase in RALDH2 transcription. Among the three RA receptors, estrogen selectively induced the expression of RARα. In parallel, estrogen also increased the utilization of all-trans retinol (the substrate for RA biosynthesis) and the expression of two RA-regulated marker genes, cellular retinoic acid binding protein II (CRABP2) and tissue transglutaminase (tTG) in the endometrium. Thus estrogen coordinately upregulated both the production and signaling of RA in both the rodent and human endometrium. This coordinate upregulation of RA system appeared to play a role in counterbalancing the stimulatory effects of estrogen on the endometrium, since the depletion of endogenous RA in mice led to an increase in estrogen-stimulated stromal proliferation and endometrial Akt phosphorylation. In addition, I have also used a systematic approach (DNA microarray) to categorize genes and pathways affected by the ERT in the endometrium of postmenopausal women and identified a novel estrogen-regulated gene EIG121. EIG121 was exclusively expressed in the glandular epithelial cells of the endometrium and induced by estrogen in vivo and in cultured cell lines. Compared with the normal endometrium, EIG121 was highly overexpressed in type 1 endometrial cancer, but profoundly suppressed in type 2 endometrial tumors. Taken together, these studies suggested that estrogen regulates the expression of many genes of both the pro-proliferative and anti-proliferative pathways and the abnormality of these pathways may increase the risks for estrogen-dependent endometrial hyperplasia and endometrial cancer. ^

Conditional expression of Sry, and Wnt4 by gene-targeting: Studies in sexual and skeletal development

Sry and Wnt4 cDNAs were individually introduced into the ubiquitously-expressed Rosa26 ( R26) locus by gene targeting in embryonic stem (ES) cells to create a conditional gene expression system in mice. In the targeted alleles, expression of these cDNAs should be blocked by a neomycin resistance selection cassette that is flanked by loxP sites. Transgene expression should be activated after the blocking cassette is deleted by Cre recombinase. ^ To test this conditional expression system, I have bred R26-stop- Sry and R26-stop-Wnt4 heterozygotes with a MisRII-Cre mouse line that expresses Cre in the gonads of both sexes. Analysis of these two types of bigenic heterozygotes indicated that their gonads developed normally like those of wild types. However, one XX R26-Sry/R26-Sry; MisR2-Cre/+ showed epididymis-like structures resembling those of males. In contrast, only normal phenotypes were observed in XY R26-Wnt4/R26-Wnt4; MisR2-Cre /+ mice. To interpret these results, I have tested for Cre recombinase activity by Southern blot and transcription of the Sry and Wnt4 transgenes by RT-PCR. Results showed that bigenic mutants had insufficient activation of the transgenes in their gonads at E12.5 and E13.5. Therefore, the failure to observe mutant phenotypes may have resulted from low activity of MisR2-Cre recombination at the appropriate time. ^ Col2a1-Cre transgenic mice express Cre in differentiating chondrocytes. R26-Wnt4; Col2a1-Cre bigenic heterozygous mice were found to exhibit a dramatic alteration in growth presumably caused by Wnt4 overexpression during chondrogenesis. R26-Wnt4; Col2a1-Cre mice exhibited dwarfism beginning approximately 10 days after birth. In addition, they also had craniofacial abnormalities, and had delayed ossification of the lumbar vertebrate and pelvic bones. Histological analysis of the growth plates of R26-Wnt4; Col2a1-Cre mice revealed less structural organization and a delay in onset of the primary and secondary ossification centers. Molecular studies confirmed that overexpression of Wnt4 causes decreased proliferation and early maturation of chondrocytes. In addition, R26-Wnt4; Col2a1-Cre mice had decreased expression of vascular endothelial growth factor (VEGF), suggesting that defects in vascularization may contribute to the dwarf phenotype. Finally, 9-month-old R26-Wnt4; Col2a1-Cre mice had significantly more fat cells in the marrow cavities of their metaphysis long bones, implying that long-term overexpression of Wnt4may cause bone marrow pathologies. In conclusion, Wnt4 was activated by Col2a1-Cre recombinase and was overexpressed in the growth plate, resulting in aberrant proliferation and differentiation of chondrocytes, and ultimately leads to dwarfism in mice. ^

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

Downstream targets of the Rhox5 homeobox gene

The X-linked mouse Rhox gene cluster contains over 30 homeobox genes that are candidates to regulate multiple steps in male and female gametogenesis. The founding member of the Rhox gene cluster, Rhox5, is an androgen-dependent gene expressed in Sertoli cells that promotes the survival and differentiation of the adjacent male germ cells. To decipher downstream signaling pathways of Rhox5, I used in vivo and in vitro microarray profiling to identify and characterize downstream targets of Rhox5 in the testis. This led to the identification of many Rhox5 -regulated genes, two of which I focused on in more detail. One of them, Unc5c, encodes a pro-apoptotic receptor with tumor suppressor activity that I found is negatively regulated by Rhox5 through a Rhox5-response element in the Unc5c 5' untranslated region (5' UTR). Examination of other mouse Rhox family members revealed that Rhox2 and Rhox3 also have the ability to downregulate Unc5c expression. The human RHOX protein RHOXF2 also had this ability, indicating that Unc5c repression is a conserved Rhox-dependent response. The repression of Unc5c expression by Rhox5 may, in part, mediate Rhox5's pro-survival function in the testis, as I found that Unc5c mutant mice have decreased germ cell apoptosis in the testis. This along with my other data leads me to propose a model in which Rhox5 is a negative regulator upstream of Unc5c in a Sertoli-cell pathway that promotes germ-cell survival. The other Rhox5-regulated gene that I studied in detail is insulin II (Ins2). Several lines of evidence, including electrophoretic mobility shift anaylsis, promoter mutagenesis, and chromatin immuoprecipitation analysis indicated that Ins2 is a direct target of Rhox5. Structure-function analysis identified homeodomain residues and the RHOX5 amino-terminal domain crucial for conferring Ins2 inducibility. Rhox5 regulates not only the Ins2 gene but also genes encoding other secreted proteins regulating metabolism (adiponectin and resistin), the rate-liming enzyme for monosaturated fatty acid biosynthesis (SCD-1), and transcription factors crucial for regulating metabolism (the nuclear hormone receptor PPARγ). I propose that the regulation of some or all of these molecules in Sertoli cells is responsible for the Rhox5-dependent survival of the adjacent germ cells. ^

Transcriptional and translational regulators of gene expression in germ cell development

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

S100A4 is a molecular mediator of endometrial carcinoma invasion

The molecular mechanisms of endometrail cancer invasion are poorly understood. S100A4, a member of the S100 Ca2+-binding protein family, was identified by oligonucleotide microarray qRT-PCR, and IHC, to be highly overexpressed in invasive endometrial carcinomas compared to non-invasive tumors. HEC-1A endometrial cancer cells transfected with S100A4 siRNA had undetectable S100A4 protein, decreased migration and invasion. The mechanism of S100A4 upregulation in endometrial cancer remains unclear. Methylation of the S100A4 gene was detected in benign endometrial glands and grade 1 tumors with no S100A4 expression. In contrast, grade 3 endometrioid tumors with high S100A4 expression showed no methylation of the gene. 5-Aza-2'-deoxycytidine, an inhibitor of DNA methyltransferase, induced the expression of S100A4 in the less invasive EC cell line, KLE, in which the S100A4 gene is hypermethylated and minimally expressed. S100A4 was induced during TGF-β1-triggered cell scattering in HEC-1A cells, in which S100A4 was demethylated. Transfection of HEC-1A cells with S100A4 siRNA significantly reduced the effect of TGF-β1 on basal migration and invasion. Our preliminary data suggested that this upregulation was mediated by the transcription factor Snail. One Snail binding consensus site was found in the region where DNA methylation was closely correlated with S100A4 gene expression. Chromatin immunoprecipitation assay confirmed the binding of Snail to this consensus site in HEC-1A cells. In SPEC2 endometrial cancer cells, loss of Snail leads to repressed S100A4 gene expression. Similar to S100A4, Snail was overexpressed in aggressive endometrial tumors. Our study suggested that the S100A4 gene was demethylated and further upregulated by the TGF-β1 and Snail pathway in invasive endometrial cancer. S100A4 could potentially serve as a good molecular marker for invasiveness and a target for therapeutic intervention for advanced endometrial cancer. ^

Statistical Methods for Differential Expressions of Genes Detected in Multiple-Condition Experiment of Microarray

Most studies of differential gene-expressions have been conducted between two given conditions. The two-condition experimental (TCE) approach is simple in that all genes detected display a common differential expression pattern responsive to a common two-condition difference. Therefore, the genes that are differentially expressed under the other conditions other than the given two conditions are undetectable with the TCE approach. In order to address the problem, we propose a new approach called multiple-condition experiment (MCE) without replication and develop corresponding statistical methods including inference of pairs of conditions for genes, new t-statistics, and a generalized multiple-testing method for any multiple-testing procedure via a control parameter C. We applied these statistical methods to analyze our real MCE data from breast cancer cell lines and found that 85 percent of gene-expression variations were caused by genotypic effects and genotype-ANAX1 overexpression interactions, which agrees well with our expected results. We also applied our methods to the adenoma dataset of Notterman et al. and identified 93 differentially expressed genes that could not be found in TCE. The MCE approach is a conceptual breakthrough in many aspects: (a) many conditions of interests can be conducted simultaneously; (b) study of association between differential expressions of genes and conditions becomes easy; (c) it can provide more precise information for molecular classification and diagnosis of tumors; (d) it can save lot of experimental resources and time for investigators.^

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

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.

{\it cis\/} -acting elements and {\it trans\/} -acting factors that control serum amyloid A gene expression during inflammation

To understand how the serum amyloid A (SAA) genes are regulated, the cis-acting elements and trans-acting factors involved in the regulation of mouse SAA3 and rat SAA1 genes expression during inflammation were analyzed.^ To identify DNA sequences involved in the liver-specific expression of the mouse SAA3 gene, the 5$\sp\prime$ flanking region of this gene was analyzed by transient transfection studies. Results suggest that C/EBP, a liver-enriched transcription factor, plays an important role for the enhanced expression of the mouse SAA3 gene in hepatocytes.^ Transfection studies of the regulation of the expression of rat SAA1 gene indicated that a 322 bp fragment ($-$304 to +18) of the gene contains sufficient information for cytokine-induced expression of the reporter gene in a liver cell-specific manner. Further functional analysis of the 5$\sp\prime$ flanking region of the rat SAA1 gene demonstrated that a 65 bp DNA fragment ($-$138/$-$73) can confer cytokine-inducibility onto a heterologous promoter both in liver and nonliver cells. DNase I footprint and gel retardation assays identified five putative cis-regulatory elements within the 5$\sp\prime$ flanking region of the gene: one inducible element, a NF$\kappa$B binding site and four constitutive elements. Two constitutive elements, footprint regions I and III, were identified as C/EBP binding sites with region III having over a 10-fold higher affinity for C/EBP binding than region I. Functional analysis of the cis-elements indicated that C/EBP(I) and C/EBP(III) confer liver cell-specific activation onto a heterologous promoter, while sequences corresponding to the NF$\kappa$B element and C/EBP(I) impart cytokine responsiveness onto the heterologous promoter. These results suggest that C/EBP(I) possesses two functions: liver-specific activation and cytokine responsiveness. The identification of two cytokine responsive elements (NF$\kappa$B and C/EBP(I)), and two liver-specific elements (C/EBP(I) and C/EBP(III)) implies that multiple cis-acting elements are involved in the regulation of the expression of the rat SAA1 gene. The tissue-specific and cytokine-induced expression of rat SAA1 gene is likely the result of the interactions of these cis-acting elements with their cognate trans-acting factors as well as the interplay between the different cis-acting elements and their binding factors. (Abstract shortened with permission of author.) ^

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