BIOCHEMICAL CHARACTERIZATION OF BINDING PARTNERS OF TWO HSP70 CO-CHAPERONES IN SACCHAROMYCES CEREVISIAE

Cells are exposed to a variety of environmental and physiological changes including temperature, pH and nutrient availability. These changes cause stress to cells, which results in protein misfolding and altered cellular protein homeostasis. How proteins fold into their three-dimensional functional structure is a fundamental biological process with important relevance to human health. Misfolded and aggregated proteins are linked to multiple neurodegenerative diseases, cardiovascular disease and cystic fibrosis. To combat proteotoxic stress, cells deploy an array of molecular chaperones that assist in the repair or removal of misfolded proteins. Hsp70, an evolutionarily conserved molecular chaperone, promotes protein folding and helps maintain them in a functional state. Requisite co-chaperones, including nucleotide exchange factors (NEFs) strictly regulate and serve to recruit Hsp70 to distinct cellular processes or locations. In yeast and human cells, three structurally non-related cytosolic NEFs are present: Sse1 (Hsp110), Fes1 (HspBP1) and Snl1 (Bag-1). Snl1 is unique among the cytosolic NEFs as it is localized at the ER membrane with its Hsp70 binding (BAG) domain exposed to the cytosol. I discovered that Snl1 distinctly interacts with assembled ribosomes and several lines of evidence indicate that this interaction is both independent of and concurrent with binding to Hsp70 and is not dependent on membrane localization. The ribosome-binding site is identified as a short lysine-rich motif within the amino terminus of the Snl1 BAG domain distinct from the Hsp70 interaction region. In addition, I demonstrate ribosome association with the Snl1 homolog in the pathogenic fungus, Candida albicans and localize this putative NEF to a perinuclear/ER membrane, suggesting functional conservation in fungal BAG domain-containing proteins. As a first step in determining specific domain architecture in fungal BAG proteins, I present the preliminary steps of protein purification and analysis of the minimal Hsp70 binding region in in both S.cerevisiae and C. albicans Snl1. Contrary to previous in vitro evidence which showed the Fes1 NEF to interact with both cytosolic Hsp70s, Ssa and Ssb, Fes1 is shown to interact specifically with Ssa when expressed under normal cellular conditions in S. cerevisiae. This is the first reported evidence of Hsp70 binding selectivity for a cytosolic NEF, and suggests a possible mechanism to achieve specificity in Hsp70-dependent functions. Taken together, the work presented in this dissertation highlights the striking divergence among Hsp70 co-chaperones in selecting binding partners, which may correlate with their specific roles in the cell.

The Role of the Arched Helicases in Exosome-Mediated Function

RNA processing and degradation are two important functions that control gene expression and promote RNA fidelity in the cell. A major ribonuclease complex, called the exosome, is involved in both of these processes. The exosome is composed of ten essential proteins with only one catalytically active subunit, called Rrp44. While the same ten essential subunits make up both the nuclear and cytoplasmic exosome, there are nuclear and cytoplasmic exosome cofactors that promote specific exosome functions in each of the cell compartments. To date, it is unclear how the exosome distinguishes between RNA substrates. We hypothesize that compartment specific cofactors may promote the substrate specificity of the exosome. In this work, I characterize several cofactors of the exosome, both nuclear and cytoplasmic. First, I describe the arch domain, which is a unique domain in a nuclear and a cytoplasmic cofactor of the exosome. Specifically, I show that the arch domain of the nuclear exosome cofactor, Mtr4, is required for specific exosome-mediated activities and overlaps functionally with the exosome-associated exonuclease, Rrp6. Further, I show that the arch domain of Ski2 is required for the degradation of normal and aberrant mRNAs. Additionally, this work describes in detail the Mtr4 domains involved in the physical association with other RNA processing proteins. Further, I characterize the minimal Mtr4-binding region in a third exosome cofactor, Trf5. Understanding how exosome cofactors synergistically promote exosome function will provide us a better understanding of how the exosome complex precisely regulates its catalytic activities. As described here, cofactors play a major role in determining the substrate specificity of the nuclear and cytoplasmic exosome. Moreover, specific accessory domains, which are not involved in the catalytic function of the cofactor, are required for substrate targeting of the eukaryotic RNA exosome.

Functional Analysis of Drosophila Integrator Complex in snRNA 3' End Processing

Uridine-rich small nuclear RNAs (U snRNAs) play essential roles in eukaryotic gene expression by facilitating the removal of introns from mRNA precursors and the processing of the replication-dependent histone pre-mRNAs. Formation of the 3’ end of these snRNAs is carried out by a poorly characterized, twelve-membered protein complex named Integrator Complex. In the effort to understand Integrator Complex function in the formation of the snRNA 3’ end, we performed a functional RNAi screen in Drosophila S2 cells to identify protein factors required for snRNA 3’ end formation. This screen was conducted by using a fluorescence-based reporter that elicits GFP expression in response to a deficiency in snRNA processing. Besides scoring the known Integrator subunits, we identified Asunder and CG4785 as additional core members of the Integrator Complex. Additionally, we also found a conserved requirement for Cyclin C and Cdk8 in both fly and human snRNA 3’ end processing. We have further demonstrated that the kinase activity of Cdk8 is critical for snRNA 3’ end processing and is likely to function independent of its well-documented function within the Mediator Cdk8 module. Taken together, this work functionally defines the Drosophila Integrator Complex and demonstrates a novel function for Cyclin C/Cdk8 in snRNA 3’ end formation. This thesis work has also characterized an important functional interaction mediated by a microdomain within Integrator subunit 12 (IntS12) and IntS1 that is required for the activity of the Integrator Complex in processing the snRNA 3’ end. Through the development of a reporter-based functional RNAi-rescue assay in Drosophila S2 cells, we analyzed domains within IntS12 required for snRNA 3’ end formation. This analysis unexpectedly revealed that an N-terminal 30 amino acid region and not the highly conserved central PHD finger domain, is required for snRNA 3’ end cleavage. The IntS12 microdomain (1-45) functions autonomously, and is sufficient to interact and stabilize the putative scaffold protein IntS1. Our findings provide more details of the Integrator Complex for understanding the molecular mechanism of snRNA 3’ end processing. Moreover, these results lay the foundation for future studies of the complex through the identification of a novel functional domain within one subunit and the identification of additional subunits.

Modeling Techniques for the High-Resolution Interpretation of Cryo-Electron Microscopy Reconstructions

Essential biological processes are governed by organized, dynamic interactions between multiple biomolecular systems. Complexes are thus formed to enable the biological function and get dissembled as the process is completed. Examples of such processes include the translation of the messenger RNA into protein by the ribosome, the folding of proteins by chaperonins or the entry of viruses in host cells. Understanding these fundamental processes by characterizing the molecular mechanisms that enable then, would allow the (better) design of therapies and drugs. Such molecular mechanisms may be revealed trough the structural elucidation of the biomolecular assemblies at the core of these processes. Various experimental techniques may be applied to investigate the molecular architecture of biomolecular assemblies. High-resolution techniques, such as X-ray crystallography, may solve the atomic structure of the system, but are typically constrained to biomolecules of reduced flexibility and dimensions. In particular, X-ray crystallography requires the sample to form a three dimensional (3D) crystal lattice which is technically di‑cult, if not impossible, to obtain, especially for large, dynamic systems. Often these techniques solve the structure of the different constituent components within the assembly, but encounter difficulties when investigating the entire system. On the other hand, imaging techniques, such as cryo-electron microscopy (cryo-EM), are able to depict large systems in near-native environment, without requiring the formation of crystals. The structures solved by cryo-EM cover a wide range of resolutions, from very low level of detail where only the overall shape of the system is visible, to high-resolution that approach, but not yet reach, atomic level of detail. In this dissertation, several modeling methods are introduced to either integrate cryo-EM datasets with structural data from X-ray crystallography, or to directly interpret the cryo-EM reconstruction. Such computational techniques were developed with the goal of creating an atomic model for the cryo-EM data. The low-resolution reconstructions lack the level of detail to permit a direct atomic interpretation, i.e. one cannot reliably locate the atoms or amino-acid residues within the structure obtained by cryo-EM. Thereby one needs to consider additional information, for example, structural data from other sources such as X-ray crystallography, in order to enable such a high-resolution interpretation. Modeling techniques are thus developed to integrate the structural data from the different biophysical sources, examples including the work described in the manuscript I and II of this dissertation. At intermediate and high-resolution, cryo-EM reconstructions depict consistent 3D folds such as tubular features which in general correspond to alpha-helices. Such features can be annotated and later on used to build the atomic model of the system, see manuscript III as alternative. Three manuscripts are presented as part of the PhD dissertation, each introducing a computational technique that facilitates the interpretation of cryo-EM reconstructions. The first manuscript is an application paper that describes a heuristics to generate the atomic model for the protein envelope of the Rift Valley fever virus. The second manuscript introduces the evolutionary tabu search strategies to enable the integration of multiple component atomic structures with the cryo-EM map of their assembly. Finally, the third manuscript develops further the latter technique and apply it to annotate consistent 3D patterns in intermediate-resolution cryo-EM reconstructions. The first manuscript, titled An assembly model for Rift Valley fever virus, was submitted for publication in the Journal of Molecular Biology. The cryo-EM structure of the Rift Valley fever virus was previously solved at 27Å-resolution by Dr. Freiberg and collaborators. Such reconstruction shows the overall shape of the virus envelope, yet the reduced level of detail prevents the direct atomic interpretation. High-resolution structures are not yet available for the entire virus nor for the two different component glycoproteins that form its envelope. However, homology models may be generated for these glycoproteins based on similar structures that are available at atomic resolutions. The manuscript presents the steps required to identify an atomic model of the entire virus envelope, based on the low-resolution cryo-EM map of the envelope and the homology models of the two glycoproteins. Starting with the results of the exhaustive search to place the two glycoproteins, the model is built iterative by running multiple multi-body refinements to hierarchically generate models for the different regions of the envelope. The generated atomic model is supported by prior knowledge regarding virus biology and contains valuable information about the molecular architecture of the system. It provides the basis for further investigations seeking to reveal different processes in which the virus is involved such as assembly or fusion. The second manuscript was recently published in the of Journal of Structural Biology (doi:10.1016/j.jsb.2009.12.028) under the title Evolutionary tabu search strategies for the simultaneous registration of multiple atomic structures in cryo-EM reconstructions. This manuscript introduces the evolutionary tabu search strategies applied to enable a multi-body registration. This technique is a hybrid approach that combines a genetic algorithm with a tabu search strategy to promote the proper exploration of the high-dimensional search space. Similar to the Rift Valley fever virus, it is common that the structure of a large multi-component assembly is available at low-resolution from cryo-EM, while high-resolution structures are solved for the different components but lack for the entire system. Evolutionary tabu search strategies enable the building of an atomic model for the entire system by considering simultaneously the different components. Such registration indirectly introduces spatial constrains as all components need to be placed within the assembly, enabling the proper docked in the low-resolution map of the entire assembly. Along with the method description, the manuscript covers the validation, presenting the benefit of the technique in both synthetic and experimental test cases. Such approach successfully docked multiple components up to resolutions of 40Å. The third manuscript is entitled Evolutionary Bidirectional Expansion for the Annotation of Alpha Helices in Electron Cryo-Microscopy Reconstructions and was submitted for publication in the Journal of Structural Biology. The modeling approach described in this manuscript applies the evolutionary tabu search strategies in combination with the bidirectional expansion to annotate secondary structure elements in intermediate resolution cryo-EM reconstructions. In particular, secondary structure elements such as alpha helices show consistent patterns in cryo-EM data, and are visible as rod-like patterns of high density. The evolutionary tabu search strategy is applied to identify the placement of the different alpha helices, while the bidirectional expansion characterizes their length and curvature. The manuscript presents the validation of the approach at resolutions ranging between 6 and 14Å, a level of detail where alpha helices are visible. Up to resolution of 12 Å, the method measures sensitivities between 70-100% as estimated in experimental test cases, i.e. 70-100% of the alpha-helices were correctly predicted in an automatic manner in the experimental data. The three manuscripts presented in this PhD dissertation cover different computation methods for the integration and interpretation of cryo-EM reconstructions. The methods were developed in the molecular modeling software Sculptor (http://sculptor.biomachina.org) and are available for the scientific community interested in the multi-resolution modeling of cryo-EM data. The work spans a wide range of resolution covering multi-body refinement and registration at low-resolution along with annotation of consistent patterns at high-resolution. Such methods are essential for the modeling of cryo-EM data, and may be applied in other fields where similar spatial problems are encountered, such as medical imaging.

Interview with Julie Knobil

An oral interview with Juile (Hotchkiss) Knobil, research professor of physiology and then integrative biology at the Medical School, where she lectured on mammalian physiology and perinatal endocrinology.

Interview with Dianna Milewicz

An oral interview with Dianna Milewicz, M.D., Ph.D., Professor and Vice Chairman of the Department of Internal Medicine at The University of Texas Medical School at Houston and Director of the M.D./Ph.D. Program and co-Director of the Biomedical Engineering Center. Her research interests include the genetic basis of cardiovascular diseases, and understanding the effect of identified mutations on protein function. She has recently established a genetic core laboratory to provide molecular biology and genetic expertise to clinicians who want to initiate genetic studies on their patient populations.

GENOME-WIDE PROFILING UNVEILS CRITICIAL FUNCTIONS OF p53 IN HUMAN EMBRYONIC STEM CELLS

Embryonic stem cells (ESCs) possess two unique characteristics: infinite self-renewal and the potential to differentiate into almost every cell type (pluripotency). Recently, global expression analyses of metastatic breast and lung cancers revealed an ESC-like expression program or signature, specifically for cancers that are mutant for p53 function. Surprisingly, although p53 is widely recognized as the guardian of the genome, due to its roles in cell cycle checkpoints, programmed cell death or senescence, relatively little is known about p53 functions in normal cells, especially in ESCs. My hypothesis is that p53 has specific transcription regulatory functions in human ESCs (hESCs) that a) oppose pluripotency and b) protect the stem cell genome in response to DNA damage and stress signaling. In mouse ESCs, these roles are believed to coincide, as p53 promotes differentiation in response to DNA damage, but this is unexplored in hESCs. To determine the biological roles of p53, specifically in hESCs, we mapped genome-wide chromatin interactions of p53 by chromatin immunoprecipitation and massively parallel tag sequencing (ChIP-Seq), and did so under three VIdifferent conditions of hESC status: pluripotency, differentiation-initiated and DNA-damage-induced. ChIP-Seq showed that p53 is enriched at distinct, induction-specific gene loci during each of these different conditions. Microarray gene expression analysis and functional annotation of the distinct p53-target genes revealed that p53 regulates specific genes encoding developmental regulators, which are expressed in differentiation-initiated but not DNA- damaged hESCs. We further discovered that, in response to differentiation signaling, p53 binds regions of chromatin that are repressed but also poised for rapid activation by core pluripotency factors OCT4 and NANOG in pluripotent hESCs. In response to DNA damage, genes associated with migration and motility are targeted by p53; whereas, the prime targets of p53 in control of cell death are conserved for p53 regulation in both differentiation and DNA damage. Our genome-wide profiling and bioinformatics analyses show that p53 occupies a special set of developmental regulatory genes during early differentiation of hESCs and functions in an induction-specific manner. In conclusion, our research unveiled previously unknown functions of p53 in ESC biology, which augments our understanding of one of the most deregulated proteins in human cancers.

INTERACTION BETWEEN BRK AND HER2 IN BREAST CANCER

INTERACTION BETWEEN BRK AND HER2 IN BREAST CANCER Midan Ai, Ph.D. Supervisory Professor: Zhen Fan, M.D. Breast tumor kinase (Brk) is a nonreceptor protein-tyrosine kinase that is highly expressed in approximately two thirds of breast cancers but is not detectable or is expressed at very low levels in normal mammary epithelium. Brk plays important roles in promoting proliferation, survival, invasion, and metastasis of breast cancer cells, but the mechanism(s) of which remain largely unknown. Recent studies showed that Brk is frequently co-overexpressed with human epidermal growth factor receptor-2 (HER2) and is physically associated with HER2 in breast cancer. The mechanism needs to be determined. In my studies, I found that high expression of HER2 is correlated with high expression of Brk in breast cancer cell lines. Silencing HER2 expression via RNA interference in HER2 over-expressed breast cancer cells resulted in Brk protein decrease and overexpression of HER2 in HER2 low-expressed breast cancer cells up-regulated Brk expression. The mechanism study indicated that overexpression of HER2 increased Brk protein stability. Brk was degraded through a Ca2+-dependent protease pathway involving calpain and HER2 stimulated Brk expression via inhibiting calpain activity. Calpastatin is a calpain endogenous inhibitor and the calpain-calpastatin system has been implicated in a number of cell physiological functions. HER2 restrained calpain activation via up-regulating calpastatin expression and HER2 downstream signaling, MAPK pathway, was involved in the regulation. Furthermore, silencing of Brk expression by RNA interference in HER2-overexpressing breast cancer cells decreased HER2-mediated cell proliferation, survival, invasion/metastasis potential and increased cell sensitivity to HER2 kinase inhibitor, lapatinib, treatment, indicating that Brk plays important roles in regulating and mediating the oncogenic functions of HER2. The Stat3 pathway played important roles in Brk mediated cell survival and invasion/metastasis in the context of HER2-overexpressing breast cancer cells. However, transgenic mice with inducible expression of constitutively active Brk (CA) in the mammary epithelium failed to develop malignant change in the mammary glands after Brk induction for 15 months which indicated that expression of Brk protein alone was not sufficiently to induce spontaneous breast tumor. Bitransgenic mice with co-expression of HER2/neu and inducible expression of Brk in the mammary epithelium developed multifocal mammary tumors, but there were no significant difference in the tumor occurring time, tumor size, tumor weight and tumor multiplicity between the mouse group with co-expression of Brk and HER2/neu and the mouse group with HER2/neu expression only.

CHARACTERIZATION OF DIFFERENTIATION AND PROGNOSTIC BIOMARKERS ON CD8+ TUMOR-INFILTRATING LYMPHOCYTES IN METASTATIC MELANOMA

CD8+ cytotoxic T lymphocytes (CTL) frequently infiltrate tumors, yet most melanoma patients fail to undergo tumor regression. We studied the differentiation of the CD8+ tumor-infiltrating lymphocytes (TIL) from 44 metastatic melanoma patients using known T-cell differentiation markers. We also compared CD8+ TIL against the T cells from matched melanoma patients’ peripheral blood. We discovered a novel subset of CD8+ TIL co-expressing early-differentiation markers, CD27, CD28, and a late/senescent CTL differentiation marker, CD57. This CD8+CD57+ TIL expressed a cytolytic enzyme, granzyme B (GB), yet did not express another cytolytic pore-forming molecule, perforin (Perf). In contrast, the CD8+CD57+ T cells in the periphery were CD27-CD28-, and GBHi and PerfHi. We found this TIL subset was not senescent and could be induced to proliferate and differentiate into CD27-CD57+, perforinHi, mature CTL. This further differentiation was arrested by TGF-β1, an immunosuppressive cytokine known to be produced by many different kinds of tumors. Therefore, we have identified a novel subset of incompletely differentiated CD8+ TIL that resembled those found in patients with uncontrolled chronic viral infections. In a related study, we explored prognostic biomarkers in metastatic melanoma patients treated in a Phase II Adoptive Cell Therapy (ACT) trial, in which autologous TIL were expanded ex vivo with IL-2 and infused into lymphodepleted patients. We unexpectedly found a significant positive clinical association with the infused CD8+ TIL expressing B- and T- lymphocyte attentuator (BTLA), an inhibitory T-cell receptor. We found that CD8+BTLA+ TIL had a superior proliferative response to IL-2, and were more capable of autocrine IL-2 production in response to TCR stimulation compared to the CD8+BTLA- TIL. The CD8+BTLA+ TIL were less differentiated and resembled the incompletely differentiated CD8+ TIL described above. In contrast, CD8+BTLA- TIL were poorly proliferative, expressed CD45RA and killer-cell immunoglobulin-like receptors (KIRs), and exhibited a gene expression signature of T cell deletion. Surprisingly, ligation of BTLA by its cognate receptor, HVEM, enhanced the survival of CD8+BTLA+ TIL by activating Akt/PKB. Our studies provide a comprehensive characterization of CD8+ TIL differentiation in melanoma, and revealed BTLA as a novel T-cell differentiation marker along with its role in promoting T cell survival.

ROLE OF MIR-19A RELEASED BY INFLAMMATORY BREAST CANCER CELLS IN THE REGULATION OF DENDRITIC CELL FUNCTIONS: IN VITRO MODEL OF CROSSTALK IN THE TUMOR MICROENVIRONMENT OF INFLAMMATORY BREAST CANCER

Inflammatory breast cancer (IBC) is a rare but very aggressive form of locally advanced breast cancer (1-6% of total breast cancer patients in United States), with a 5-year overall survival rate of only 40.5%, compared with 85% of the non-IBC patients. So far, a unique molecular signature for IBC able to explain the dramatic differences in the tumor biology between IBC and non-IBC has not been identified. As immune cells in the tumor microenvironment plays an important role in regulating tumor progression, we hypothesized that tumor-associated dendritic cells (TADC) may be responsible for regulating the development of the aggressive characteristics of IBC. MiRNAs can be released into the extracellular space and mediate the intercellular communication by regulating target gene expression beyond their cells of origin. We hypothesized that miRNAs released by IBC cells can induce an increased activation status, secretion of pro-inflammatory cytokines and migration ability of TADC. In an in vitro model of IBC tumor microenvironment, we found that the co-cultured of the IBC cell line SUM-149 with immature dendritic cells (iDCSUM-149) induced a higher degree of activation and maturation of iDCSUM-149 upon stimulation with lipopolysaccharide (LPS) compared with iDCs co-cultured with the non-IBC cell line SUM-159 (iDCSUM-159), resulting in: increased expression of the costimulatory and activation markers; higher production of pro-inflammatory cytokines (TNF-a, IL-6); and 3) higher migratory ability. These differences were due to the exosome-mediated transfer of miR-19a and miR-146a from SUM-149 and SUM-159, respectively, to iDCs, causing the downregulation of the miR-19a target genes PTEN, SOCS-1 and the miR-146a target genes IRAK1, TRAF6. PTEN, SOCS-1 and IRAK1, TRAF6 are important negative and positive regulator of cytokine- and TLR-mediated activation/maturation signaling pathway in DCs. Increased levels of IL-6 induced the upregulation of miR-19a synthesis in SUM-149 cells that was associated with the induction of CD44+CD24-ALDH1+ cancer stem cells (CSCs) with epithelial-to-mesenchymal transition (EMT) characteristics. In conclusion, in IBC tumor microenvironment IL-6/miR-19a axis can represent a self-sustaining loop able to maintain a pro-inflammatory status of DCs, leading to the development of tumor cells with high metastatic potential (EMT CSCs) responsible of the poor prognosis in IBC patients.

PRKCa: Identification of a Novel Downstream Target of WT1

Wilms tumor is a childhood tumor of the kidney arising from the undifferentiated metanephric mesenchyme. Tumorigenesis is attributed to a number of genetic and epigenetic alterations. In 20% of Wilms tumors, Wilms tumor gene 1 (WT1) undergoes inactivating homozygous mutations causing loss of function of the zinc finger transcription factor it encodes. It is hypothesized that mutations in WT1 result in dysregulation of downstream target genes, leading to aberrant kidney development and/or Wilms tumor. These downstream target genes are largely unknown, and identification is important for further understanding Wilms tumor development. Heatmap data of human Wilms tumor protein expression, generated by reverse phase protein assay analysis (RPPA), show significant correlation between WT1 mutation status and low PRKCα expression (p= 0.00013); additionally, p-PRKCα (S657) also shows decreased expression in these samples (p= 0.00373). These data suggest that the WT1 transcription factor regulates PRKCα expression, and that PRKCα plays a potential role in Wilms tumor tumorigenesis. We hypothesize that the WT1 transcription factor directly/indirectly regulates PRKCα and mutations occurring in WT1 lead to decreased expression of PRKCα. Prkcα and Wt1 have been shown to co-localize in E14.5 mesenchymal cells of the developing kidney. siRNA knockdown, in-vivo ablation, and tet-inducible expression of Wt1 each independently confirm regulation of Prkcα expression by Wt1 at both RNA and protein levels, and investigation into possible WT1 binding sites in PRKCα regulatory regions has identified multiple sites to be confirmed by luciferase reporter constructs. With the goal of identifying WT1 and PRKCα downstream targets, RPPA analysis of protein expression in mesenchymal cell culture, following lentiviral delivered shRNA knockdown of Wt1 and shRNA knockdown of Prkcα, will be carried out. Apart from Wilms tumor, WT1 also plays an important role in Acute Myeloid Leukemia (AML). WT1 mutation status has been implicated, controversially, as an independent poor-prognosis factor in leukemia, leading to decreased probability of overall survival, complete remission, and disease free survival. RPPA analysis of AML patient samples showed significant decreases in PRKCα/p-PRKCα protein expression in a subset of patients (Kornblau, personal communication); therefore, the possible role of WT1 and PRKCα in leukemia disease progression is an additional focus of this study. WT1 mutation analysis of diploid leukemia patient samples revealed two patients with mutations predicted to affect WT1 activity; of these two samples, only one corresponded to the low PRKCα expression cohort. Further characterization of the role of WT1 in AML, and further understanding of WT1 regulated PRKCα expression, will be gained following RPPA analysis of protein expression in HL60 leukemia cell lines with lentiviral delivered shRNA knockdown of WT1 and shRNA knockdown of PRKCα.

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.

Investigation of I -compounds (newly discovered indigenous DNA modifications) by different approaches

I-compounds are newly discovered covalent DNA modifications detected by the $\sp{32}$P-postlabeling assay. They are age-dependent, tissue-specific and sex-different. The origin(s), chemistry and function(s) of I-compounds are unknown. The total level of I-compounds in 8-10 month old rat liver is 1 adduct in 10$\sp7$ nucleotides, which is not neglectable. It is proposed that I-compounds may play a role in spontaneous tumorigenesis and aging.^ In the present project, I-compounds were investigated by several different approaches. (1) Dietary modulation of I-compounds. (2) Comparison of I-compounds with persistent carcinogen DNA adducts and 5-methylcytosine. (3) Strain differences of I-compounds in relation to organ site spontaneous tumorigenesis. (4) Effects of nongenotoxic hepatocarcinogenes on I-compounds.^ It was demonstrated that the formation of I-compounds is diet-related. Rats fed natural ingredient diet exhibited more complex I-spot patterns and much higher levels than rats fed purified diet. Variation of major nutrients (carbohydrate, protein and fat) in the diet, produced quantitative differences in I-compounds of rat liver and kidney DNAs. Physiological level of vitamin E in the diet reduced intensity of one I-spot compared with vitamin E deficient diet. However, extremely high level of vitamin E in the diet gave extra spot and enhanced the intensities of some I-spots.^ In regenerating rat liver, I-compounds levels were reduced, as carcinogen DNA adducts, but not 5-methylcytosine, i.e. a normal DNA modification.^ Animals with higher incidences of spontaneous tumor or degenerative diseases tended to have a lower level of I-compounds.^ Choline devoid diet induced a drastic reduction of I-compound level in rat liver compared with choline supplemented diet. I-compound levels were reduced after multi-doses of carbon tetrachloride (CCl$\sb4$) exposure in rats and single dose exposure in mice. An inverse relationship was observed between I-compound level and DNA replication rate. CCl$\sb4$-related DNA adduct was detected in mice liver and intensities of some I-spots were enhanced 24 h after a single dose exposure.^ The mechanisms and explanations of these observations will be discussed. I-compounds are potentially useful indicators in carcinogenesis studies. ^

Estrogen induced desensitization of c -{\it fos\/} messenger RNA expression {\it in vivo\/}

In this thesis, we investigated the regulation of the nuclear proto-oncogene, c-fos by estrogen in vivo. In the uterus, estrogen causes a rapid, dramatic and transient induction of c-fos mRNA and this occurs by transcriptional activation. We have discovered a previously unrecognized regulatory mechanism by which fos becomes desensitized to estrogen following the transient induction. We investigated three aspects of this desensitization: (1) the kinetics and general characteristics of the phenomenon; (2) the molecular mechanism of the desensitization; and (3) the relationship of desensitization to estrogen stimulated DNA synthesis. The desensitization occurs between 3-24 hours after initial hormonal stimulation and is reversible within 72 hours. The desensitization is not species specific, in that it occurs in both the rat and mouse. The desensitization also occurs in at least two estrogen responsive tissues, the uterus and vagina. The desensitization is not unique to c-fos, since both c-myc and c-jun show similar patterns of desensitization. However, the desensitization is not observed with creatine kinase B (CKB), indicating that not all estrogen inducible genes become desensitized. In the second general area, we determined the desensitization is at the transcriptional level. The desensitization is homologous, but not heterologous, since estrogen induction does not desensitize c-fos to other agents. Other studies show that the desensitization is not due to the lack of functional estrogen receptors. Taken together, these findings suggest that the desensitization occurs at the level of the estrogen responsive element. In the third major area, we demonstrated that the desensitization appears to be related to estrogen induced DNA synthesis. Support for this suggestion comes from the observation that short acting estrogens which induce fos, but not DNA synthesis, do not produce desensitization. ^

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