The cellular and molecular involvement of glutathione in cisplatin-induced apoptosis

The goal of this study was to investigate the cellular and molecular mechanisms by which glutathione (GSH) is involved in the process of apoptosis induced by cisplatin [cis-diamminedichloroplatinum(II), cis-DDP] in the HL60 human promyelocytic leukemia cell line. The data show that during the onset or induction of apoptosis, GSH levels in cisplatin-treated cells increased 50% compared to control cells. The increase in intracellular GSH was associated with enhanced expression of γ-glutamylcysteine synthetase (γ-GCS), the enzyme that catalyzes the rate- limiting step in the biosynthesis of glutathione. After depletion of intracellular GSH with D,L-buthionine-(S,R)-sulfoximine (BSO), an inhibitor of γ-GCS, biochemical and morphological analysis revealed that the mechanism of cell death had switched from apoptosis to necrosis. In contrast, when intracellular GSH was elevated by exposure of cells to a GSH-ethyl-ester and then treatment with cisplatin, no change in the induction and kinetics of apoptosis were observed. However, when cells were exposed to cisplatin before intracellular GSH levels were increased, apoptosis was observed to occur 6 hours earlier compared to cells without GSH elevation. To further examine the molecular aspects of these effects of GSH on the apoptotic process, changes in the expression of bcl-2 and bax, were investigated in cells with depleted and elevated GSH. Using reverse transcription polymerase chain reaction, no significant change in the expression of bcl-2 gene transcripts was observed in cells in either the GSH depleted or elevated state; however, a 75% reduction in GSH resulted in a 40% decrease in the expression of bax gene transcripts. In contrast, a 6-fold increase in GSH increased the expression of bax by 3-fold relative to controls. Similar results were obtained for bax gene expression and protein synthesis by northern analysis and immunoprecipitation, respectively. These results suggest that GSH serves a dual role in the apoptotic process. The first role which is indirect, involves the protection of the cell from extensive damage following exposure to a specific toxicant so as to prevent death by necrosis, possibly by interacting with the DNA damaging agent and/or its active metabolites. The second role involves a direct involvement of GSH in the apoptotic process that includes upregulation of bax expression. ^

INTEGRATIVE BIOMARKER IDENTIFICATION AND CLASSIFICATION USING HIGH THROUGHPUT ASSAYS

It is well accepted that tumorigenesis is a multi-step procedure involving aberrant functioning of genes regulating cell proliferation, differentiation, apoptosis, genome stability, angiogenesis and motility. To obtain a full understanding of tumorigenesis, it is necessary to collect information on all aspects of cell activity. Recent advances in high throughput technologies allow biologists to generate massive amounts of data, more than might have been imagined decades ago. These advances have made it possible to launch comprehensive projects such as (TCGA) and (ICGC) which systematically characterize the molecular fingerprints of cancer cells using gene expression, methylation, copy number, microRNA and SNP microarrays as well as next generation sequencing assays interrogating somatic mutation, insertion, deletion, translocation and structural rearrangements. Given the massive amount of data, a major challenge is to integrate information from multiple sources and formulate testable hypotheses. This thesis focuses on developing methodologies for integrative analyses of genomic assays profiled on the same set of samples. We have developed several novel methods for integrative biomarker identification and cancer classification. We introduce a regression-based approach to identify biomarkers predictive to therapy response or survival by integrating multiple assays including gene expression, methylation and copy number data through penalized regression. To identify key cancer-specific genes accounting for multiple mechanisms of regulation, we have developed the integIRTy software that provides robust and reliable inferences about gene alteration by automatically adjusting for sample heterogeneity as well as technical artifacts using Item Response Theory. To cope with the increasing need for accurate cancer diagnosis and individualized therapy, we have developed a robust and powerful algorithm called SIBER to systematically identify bimodally expressed genes using next generation RNAseq data. We have shown that prediction models built from these bimodal genes have the same accuracy as models built from all genes. Further, prediction models with dichotomized gene expression measurements based on their bimodal shapes still perform well. The effectiveness of outcome prediction using discretized signals paves the road for more accurate and interpretable cancer classification by integrating signals from multiple sources.

Bifurcation and singularity analysis of a molecular network for the induction of long-term memory.

Withdrawal reflexes of the mollusk Aplysia exhibit sensitization, a simple form of long-term memory (LTM). Sensitization is due, in part, to long-term facilitation (LTF) of sensorimotor neuron synapses. LTF is induced by the modulatory actions of serotonin (5-HT). Pettigrew et al. developed a computational model of the nonlinear intracellular signaling and gene network that underlies the induction of 5-HT-induced LTF. The model simulated empirical observations that repeated applications of 5-HT induce persistent activation of protein kinase A (PKA) and that this persistent activation requires a suprathreshold exposure of 5-HT. This study extends the analysis of the Pettigrew model by applying bifurcation analysis, singularity theory, and numerical simulation. Using singularity theory, classification diagrams of parameter space were constructed, identifying regions with qualitatively different steady-state behaviors. The graphical representation of these regions illustrates the robustness of these regions to changes in model parameters. Because persistent protein kinase A (PKA) activity correlates with Aplysia LTM, the analysis focuses on a positive feedback loop in the model that tends to maintain PKA activity. In this loop, PKA phosphorylates a transcription factor (TF-1), thereby increasing the expression of an ubiquitin hydrolase (Ap-Uch). Ap-Uch then acts to increase PKA activity, closing the loop. This positive feedback loop manifests multiple, coexisting steady states, or multiplicity, which provides a mechanism for a bistable switch in PKA activity. After the removal of 5-HT, the PKA activity either returns to its basal level (reversible switch) or remains at a high level (irreversible switch). Such an irreversible switch might be a mechanism that contributes to the persistence of LTM. The classification diagrams also identify parameters and processes that might be manipulated, perhaps pharmacologically, to enhance the induction of memory. Rational drug design, to affect complex processes such as memory formation, can benefit from this type of analysis.

A model of late long-term potentiation simulates aspects of memory maintenance.

Late long-term potentiation (L-LTP) denotes long-lasting strengthening of synapses between neurons. L-LTP appears essential for the formation of long-term memory, with memories at least partly encoded by patterns of strengthened synapses. How memories are preserved for months or years, despite molecular turnover, is not well understood. Ongoing recurrent neuronal activity, during memory recall or during sleep, has been hypothesized to preferentially potentiate strong synapses, preserving memories. This hypothesis has not been evaluated in the context of a mathematical model representing ongoing activity and biochemical pathways important for L-LTP. In this study, ongoing activity was incorporated into two such models - a reduced model that represents some of the essential biochemical processes, and a more detailed published model. The reduced model represents synaptic tagging and gene induction simply and intuitively, and the detailed model adds activation of essential kinases by Ca(2+). Ongoing activity was modeled as continual brief elevations of Ca(2+). In each model, two stable states of synaptic strength/weight resulted. Positive feedback between synaptic weight and the amplitude of ongoing Ca(2+) transients underlies this bistability. A tetanic or theta-burst stimulus switches a model synapse from a low basal weight to a high weight that is stabilized by ongoing activity. Bistability was robust to parameter variations in both models. Simulations illustrated that prolonged periods of decreased activity reset synaptic strengths to low values, suggesting a plausible forgetting mechanism. However, episodic activity with shorter inactive intervals maintained strong synapses. Both models support experimental predictions. Tests of these predictions are expected to further understanding of how neuronal activity is coupled to maintenance of synaptic strength. Further investigations that examine the dynamics of activity and synaptic maintenance can be expected to help in understanding how memories are preserved for up to a lifetime in animals including humans.

In vivo definition of genetic and cellular aspects of mammalian sexual differentiation

The differentiation of the reproductive organs is an essential developmental process required for the proper transmission of the genetic material. Müllerian inhibiting substance (MIS) is produced by testes and is necessary for the regression of the Müllerian ducts: the anlagen of the uterus, fallopian tubes and cervix. In vitro and standard transgenic mouse studies indicate that the nuclear hormone receptor Steroidogenic factor 1 (SF-1) and the transcription factor SOX9 play an essential role in the regulation of Mis. To test this hypothesis, mutations in the endogenous SF-1 and SOX9 binding sites in the mouse Mis promoter were introduced by gene targeting in embryonic stem (ES) cells. In disagreement with cell culture and transgenic mouse studies, male mice homozygous for the mutant SF-1 binding site correctly initiated Mis transcription in the fetal testes, although at significantly reduced levels. Surprisingly, sufficient Mis was produced for complete elimination of the Müllerian duct system. However, when the SF-1 binding site mutation was combined with an Mis -null allele, the further decrease in Mis levels led to a partial retention of uterine tissue, but only at a distance from the testes. In contrast, males homozygous for the mutant SOX9 binding site did not initiate Mis transcription, resulting in pseudohermaphrodites with a uterus and oviducts. These studies suggest an essential role for SOX9 in the initiation of Mis transcription, whereas SF-1 appears to act as a quantitative regulator of Mis transcript levels perhaps for influencing non-Müllerian duct tissues. ^ The Mis type II receptor, a member of the TGF- b superfamily, is also required for the proper regression of the Müllerian ducts. Mis type II receptor-deficient human males and their murine counterparts develop as pseudohermaphrodites. A lacZ reporter cassette was introduced into the mouse Mis type II receptor gene, by homologous recombination in ES cells. Expression studies, based on b -galactosidase activity, show marked expression of the MIS type II receptor in the postnatal Sertoli cells of the testis as well as in the prenatal and postnatal granulosa cells of the ovary. Expression is also seen in the mesenchymal cells surrounding the Müllerian duct and in the longitudinal muscle layer of the uterus. ^

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

Identification and characterization of a molecular marker in epithelial ovarian cancer

Carcinomas that arise from the ovarian surface epithelium represent a great challenge in gynecologic oncology. Although the prognosis of ovarian cancer is influenced by many factors capable of predicting clinical outcome, including tumor stage, pathological grade, and amount of residual disease following primary surgery, the biological aspects of ovarian cancer are not completely understood, thus implying that there may be other predictive indicators that could be used independently or in conjunction with these factors to provide a clearer clinical picture. The identification of additional markers with biological relevance is desirable. To identify disease-associated peptides, a phage display random peptide library was used to screen immunoglobulins derived from a patient with ovarian cancer. One peptide was markedly enriched following three rounds of affinity selection. The presence of autoantibodies against the peptide was examined in a panel of ovarian cancer patients. Stage IV patients exhibited a high percentage of positive reactivity (59%). This was in contrast to stage III patients, who only displayed 7% positive reactivity. Antibodies against the peptide were affinity purified, and heat-shock protein 90 (Hsp90) was identified as the corresponding autoantigen. The expression profile of the identified antigen was determined. Hsp90 was expressed in all sections examined regardless of degree of anaplasia. This thesis shows that utilizing the humoral response to ovarian cancer can be used to identify a tumor antigen in ovarian cancer. The data show that certain antigens may be expressed in ovarian tumors independent of the disease stage or grade, whereas circulating antibodies against such epitopes are only found in a subset of patients. ^

Signal relay between two integral membrane proteins: The sensory rhodopsin I/HtrI transducer molecular complex

The molecular complex of sensory rhodopsin I (SRI) and its transducer HtrI mediate color-sensitive phototaxis in the archaeon Halobacterium salinarum. Orange light causes an attractant response by a one-photon reaction and white light causes a repellent response by a two-photon reaction. Three aspects of this molecular complex were explored: (i) We determined the stoichiometry of SRI and HtrI to be 2:2 by gene fusion analysis. A SRI-HtrI fusion protein was expressed in H. salinarum and shown to mediate 1-photon and 2-photon phototaxis responses comparable to wild-type complex. Disulfide crosslinking demonstrated that the fusion protein is a homodimer in the membrane. Measurement of photochemical reaction kinetics and pH titration of absorption spectra established that both SRI domains are complexed to HtrI in the fusion protein, and therefore the stoichiometry is 2:2. (ii) Cytoplasmic channel closure of SRI by HtrI, an important aspect of their interaction, was investigated by incremental HtrI truncation. We found that binding of the membrane-embedded portion of HtrI is insufficient for channel closure, whereas cytoplasmic extension of the second HtrI transmembrane helix by 13 residues blocks proton conduction through the channel as well as full-length HtrI. The closure activity is localized to 5 specific residues, each of which incrementally contributes to reduction of proton conductivity. Moreover, these same residues in the dark incrementally and proportionally increase the pKa of the Asp76 counterion to the protonated Schiff base chromophore. We conclude that this critical region of HtrI alters the dark conformation of SRI as well as light-induced channel opening. (iii) We developed a procedure for reconstituting HtrI-free SRI and the SRI/HtrI complex into liposomes, which exhibit photocycles with opened and closed cytoplasmic channels, respectively, as in the membrane. This opens the way for study of the light-induced conformational change and the interaction in vitro by fluorescence and spin-labeling. Single-cysteine mutations were introduced into helix F of SRI, labeled with a nitroxide spin probe and a fluorescence probe, reconstituted into proteoliposomes, and light-induced conformational changes detected in the complex. The probe signals can now be used as the readout of signaling to analyze mutants and the kinetics of signal relay. ^

Spatial and temporal aspects of synaptic plasticity

The notion that changes in synaptic efficacy underlie learning and memory processes is now widely accepted even if definitive proof of the synaptic plasticity and memory hypothesis is still lacking. When learning occurs, patterns of neural activity representing the occurrence of events cause changes in the strength of synaptic connections within the brain. Reactivation of these altered connections constitutes the experience of memory for these events and for other events with which they may be associated. These statements summarize a long-standing theory of memory formation that we refer to as the synaptic plasticity and memory hypothesis. Since activity-dependent synaptic plasticity is induced at appropriate synapses during memory formation, and is both necessary and sufficient for the information storage, we can speculate that a methodological study of the synapse will help us understand the mechanism of learning. Random events underlie a wide range of biological processes as diverse as genetic drift and molecular diffusion, regulation of gene expression and neural network function. Additionally spatial variability may be important especially in systems with nonlinear behavior. Since synapse is a complex biological system we expect that stochasticity as well as spatial gradients of different enzymes may be significant for induction of plasticity. ^ In that study we address the question "how important spatial and temporal aspects of synaptic plasticity may be". We developed methods to justify our basic assumptions and examined the main sources of variability of calcium dynamics. Among them, a physiological method to estimate the number of postsynaptic receptors as well as a hybrid algorithm for simulating postsynaptic calcium dynamics. Additionally we studied how synaptic geometry may enhance any possible spatial gradient of calcium dynamics and how that spatial variability affect plasticity curves. Finally, we explored the potential of structural synaptic plasticity to provide a metaplasticity mechanism specific for the synapse. ^

Systematic interpretation of molecular beacon polymerase chain reaction for identifying rpoB mutations in Mycobacterium tuberculosis isolates with mixed resistant and susceptible bacteria

Detection of multidrug-resistant tuberculosis (MDR-TB), a frequent cause of treatment failure, takes 2 or more weeks to identify by culture. RIF-resistance is a hallmark of MDR-TB, and detection of mutations in the rpoB gene of Mycobacterium tuberculosis using molecular beacon probes with real-time quantitative polymerase chain reaction (qPCR) is a novel approach that takes ≤2 days. However, qPCR identification of resistant isolates, particularly for isolates with mixed RIF-susceptible and RIF-resistant bacteria, is reader dependent and limits its clinical use. The aim of this study was to develop an objective, reader-independent method to define rpoB mutants using beacon qPCR. This would facilitate the transition from a research protocol to the clinical setting, where high-throughput methods with objective interpretation are required. For this, DNAs from 107 M. tuberculosis clinical isolates with known susceptibility to RIF by culture-based methods were obtained from 2 regions where isolates have not previously been subjected to evaluation using molecular beacon qPCR: the Texas–Mexico border and Colombia. Using coded DNA specimens, mutations within an 81-bp hot spot region of rpoB were established by qPCR with 5 beacons spanning this region. Visual and mathematical approaches were used to establish whether the qPCR cycle threshold of the experimental isolate was significantly higher (mutant) compared to a reference wild-type isolate. Visual classification of the beacon qPCR required reader training for strains with a mixture of RIF-susceptible and RIF-resistant bacteria. Only then had the visual interpretation by an experienced reader had 100% sensitivity and 94.6% specificity versus RIF-resistance by culture phenotype and 98.1% sensitivity and 100% specificity versus mutations based on DNA sequence. The mathematical approach was 98% sensitive and 94.5% specific versus culture and 96.2% sensitive and 100% specific versus DNA sequence. Our findings indicate the mathematical approach has advantages over the visual reading, in that it uses a Microsoft Excel template to eliminate reader bias or inexperience, and allows objective interpretation from high-throughput analyses even in the presence of a mixture of RIF-resistant and RIF-susceptible isolates without the need for reader training.^

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

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.