NOVEL MECHANISMS OF ANTIGEN PROCESSING THAT ENHANCE BCG VACCINE EFFICACY

Mycobacterium tuberculosis, the causative agent of tuberculosis, is the most lethal single infectious agent afflicting man today causing 2 million deaths per year. The World Health Organization recommends a vaccine as the best option to prevent this disease. The current vaccine, BCG, has a variable efficacy and does not protect adults. It is known that BCG vaccine becomes sequestered in special phagosome compartments of macrophages that do not fuse with lysosomes. Since lysosome fusion is necessary for peptide production and T cell priming leading to protective TH1 immunity, we hypothesized that vaccine efficacy is reduced and occurs perhaps due to non-lysosome dependent mechanisms. We therefore proposed an in depth analysis of phagosome environment, and its proteome to unravel mechanisms of antigen processing and presentation. We initially discovered that three mechanisms of pH regulation including vacuolar proton ATPase, phagocyte oxidase and superoxide dismutase (SOD) secretion from BCG vaccine affect antigen processing within phagosomes. These studies led to the discovery that a mutant of BCG vaccine which lacked SOD was a better vaccine. Subsequently, the proteomic analysis of vaccine phagosomes led to the discovery of novel protease (γ-secretase) enriched on BCG vaccine phagosomes. We then demonstrated that these proteases generated a peptide from the BCG vaccine which was presented through the MHC-II pathway to T cells and induced a TH1 response. The specificity of antigen production from γ-secretase was confirmed through siRNA knockdown of the components of the protease namely, nicastrin, presenilin and APH, which led to a decrease in antigen presentation. We therefore conclude that, even though BCG phagosomes are sequestered and do not fuse with lysosomes to generate peptide antigens, there are complex and novel in situ mechanisms within phagosomes that are capable of generating an immune response. We conclude that TH1 immunity to BCG vaccine arises mostly due to non-lysosome dependent immune mechanisms of macrophages and dendritic cells.

Glycogen Synthase Kinase 3 is Required for Optimal AKT Activation

The phosphatidylinositol 3-kinase (PI3K) pathway, through its major effector node AKT, is critical for the promotion of cell growth, division, motility and apoptosis evasion. This signaling axis is therefore commonly targeted in the form of mutations and amplifications in a myriad of malignancies. Glycogen synthase kinase 3 (GSK3) was first discovered as the kinase responsible for phosphorylating and inhibiting the activity of glycogen synthase, ultimately antagonizing the storage of glucose as glycogen. Its activity counteracts the effects of insulin in glucose metabolism and AKT has long been recognized as one of the key molecules capable of phosphorylating GSK3 and inhibiting its activity. However, here we demonstrate that GSK3 is required for optimal phosphorylation and activation of AKT in different malignant cell lines, and that this effect is independent of the type of growth factor stimulation and can happen even in basal states. Both GSK3 alpha and GSK3 beta isoforms are necessary for AKT to become fully active, displaying a redundant role in the setting. We also demonstrate that this effect of GSK3 on AKT phosphorylation and full activation is dependent on its kinase activity, since highly specific inhibitors targeting GSK3 catalytic activity also promote a reduction in phosphorylated AKT. Analysis of reverse phase protein array screening of MDA-MB-231 breast cancer cells treated with RNA interference targeting GSK3 unexpectedly revealed an increase in levels of phosphorylated MAPK14 (p38). Treatment with the selective p38 inhibitor SB 202190 rescued AKT activation in that cell line, corroborating the importance of unbiased proteomic analysis in exposing cross-talks between signaling networks and demonstrating a critical role for p38 in the regulation of AKT phosphorylation.

LMW-E MEDIATES MAMMARY TUMORIGENESIS BY DEREGULATING ACINAR MORPHOGENESIS & GENERATING CANCER STEM CELLS

Cyclin E is the regulatory subunit of the cyclin E/CDK2 complex that mediates the G1-S phase transition. N-terminal cleavage of cyclin E by elastase in breast cancer generates two low molecular weight (LMW) isoforms that exhibit both enhanced kinase activity and resistance to p21 and p27 inhibition compared to fulllength cyclin E. Clinically, approximately 27% of breast cancer patients overexpress LMW-E and associate with poor survival. Therefore, we hypothesize that LMW-E disrupts normal mammary acinar morphogenesis and serves as the initial route into breast tumor development. We first demonstrate that LMW-E overexpression in non-tumorigenic hMECs is sufficient to induce tumor formation in athymic mice significantly more than overexpression of full-length cyclin E and requires CDK2- associated kinase activity. Further in vivo passaging of these tumors augments LMW-E expression and tumorigenic potential. When subjected to acinar morphogenesis in vitro, LMW-E mediates significant morphological disruption by generating hyperproliferative and multi-acinar complexes. Proteomic analysis of patient tissues and tumor cells with high LMW-E expression reveals that the activation of the b-Raf-ERK1/2-mTOR pathway in concert with high LMW-E expression predicts poor patient survival. Combination treatment using roscovitine (CDK inhibitor) plus either rapamycin (mTOR inhibitor) or sorafenib (b-raf inhibitor) effectively prevented aberrant acinar formation in LMW-E-expressing cells by inducing the G1/S cell cycle arrest. In addition, the LMW-E-expressing tumor cells exhibit phenotypes characteristic of the EMT and enhanced cellular invasiveness. These tumor cells also enrich for cells with CSC phenotypes such as increased CD44hi/CD24lo population, enhanced mammosphere formation, and upregulation of ALDH expression and enzymatic activity. Furthermore, the CD44hi/CD24lo population also shows positive correlation with LMW-E expression in both the tumor cell line model and breast cancer patient samples (p<0.0001 & p=0.0435, respectively). Combination treatment using doxorubicin and salinomycin demonstrates synergistic cytotoxic effects in cells with LMW-E expression but not in those with full-length cyclin E expression. Finally, ProtoArray microarray identifies Hbo1 as a novel substrate of the cyclin E/CDK2 complex and its overexpression results in enrichment for CSCs. Collectively, these data emphasize the strong oncogenic potential of LMW-E in mammary tumorigenesis and suggest possible therapeutic strategies to treat breast cancer patients with high LMW-E expression.

Roles of IkappaB kinase alpha in centrosome duplication and skin carcinogenesis

IκB kinase α (IKKα) is one kinase subunit of the IKK complex that is responsible for NF-κB activation. Previous studies have shown that IKKα determines mouse keratinocyte terminal differentiation independent of the NF-κB pathway. Accumulating evidence suggests that IKKα functions as a tumor suppressor in skin carcinogenesis; however, the downstream pathways mediating this function are largely unknown. By using primary cultured keratinocytes, we found that Ikkα-/- cells developed aneuploidy and underwent spontaneous immortalization and transformation while wild type cells underwent terminal differentiation in the same culture condition. Using proteomic analysis we identified nucleophosmin (NPM), a centrosome duplication regulator, as an IKKα substrate. We further demonstrated that IKKα interacted with NPM and colocalized with NPM on the centrosome, suggesting that NPM is a physiological substrate of IKKα. Loss of IKKα reduced centrosome-bound NPM and promoted abnormal centrosome amplification, which contributed to aneuploidy development. Detailed analysis revealed that ablation of IKKα target site serine-125 of NPM induced destabilization of NPM hexamers, disrupted NPM association with centrosomes, and resulted in abnormal centrosome amplification. Re-introduction of IKKα rescued the defect in Ikkα-/- keratinocytes. Thus, IKKα is required for maintaining proper centrosome duplication by phosphorylating NPM. ^ UV is the major etiological agent for human skin cancer and UV-induced mouse skin carcinogenesis is one of the most relevant experimental models for human skin carcinogenesis. Thus, we further evaluated IKKα function in UV-induced skin carcinogenesis in Ikkα+/- mice. We demonstrated that IKKα is also critical in UV skin carcinogenesis, as evidenced by increased tumor multiplicity and reduced tumor latency in Ikkα+/- mice after chronic UVB treatment. Reduced expression of IKKα decreased UV-induced apoptosis and promoted accumulation of P53 mutations in the epidermis. This indicates that IKKα is critical for UV-induced apoptosis in vivo and thus prevents mutation accumulation that is important for tumor development. ^ Together, these findings uncover previously unknown in vivo functions of IKKα in centrosome duplication and apoptosis, thus providing a possible mechanism of how loss of IKKα may contribute to skin carcinogenesis. ^

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.

Identification and characterization of cellular genes differentially expressed in a subset of tumors with non -functional WT1

Wilms tumor (WT) or nephroblastoma is a genetically heterogeneous pediatric renal tumor that accounts for 6–7% of all childhood cancers in the U.S. WT1, located at 11p13, is the sole WT gene cloned to date. Additional genomic regions containing genes that play a role in the development of Wilms tumor include 11p15, 7p, 16q, 1p, 17q and 19q. This heterogeneity has made it extremely difficult to develop an understanding of the pathways involved in the development of WT, even in the 5–20% of tumors that show mutations at the WT1 locus. My research addresses this gap in our current comprehension of the development of WT. ^ I have used two complementary approaches to extend the current understanding of molecular changes involved in the development of WT. In order to minimize complexities due to genetic heterogeneity, I confined my analysis to the WT1 pathway by assessing those genetically defined tumors that carry WT1 mutations. WT1 encodes a zinc finger transcription factor, and in vitro studies have identified many genes that are potentially regulated in vivo by WT1. However, there is very little in vivo data that suggests that they are transcriptionally regulated endogenously by WT1. In one approach I assessed the role of WT1 in the in vivo regulation of PDGFA and IGF2, two genes that are strong contenders for endogenous regulation by WT1. Using primary tissue samples, I found no correlation between the level of RNA expression of WT1 with either PDGFA or IGF2, suggesting that WT1 does not play a critical role in their expression in either normal kidney or WT. ^ In a parallel strategy, using differential display analysis I compared global gene expression in a subset of tumors with known homozygous inactivating WT1 mutations (WT1-tumors) to the gene expression in a panel of appropriate control tissues (fetal kidney, normal kidney, rhabdoid tumor and pediatric renal cell carcinoma). Transcripts that are aberrantly expressed in this subset of Wilms tumors are candidates for endogenous transcriptional regulation by WT1 as well as for potentially functioning in the development of WT. By this approach I identified several differentially expressed transcripts. I further characterized two of these transcripts, identifying a candidate WT gene in the process. I then performed a detailed analysis of this WT candidate gene, which maps to 7p. Future studies will shed more light on the role of these differentially expressed genes in WT. ^

Transcriptional analysis of liver-specific expression and differential interleukin-6 response of rat kininogen genes

Analyses of rat T1 kininogen gene/chloramphenicol acetyltransferase (T1K/CAT) constructs revealed two regions important for tissue-specific and induced regulation of T1 kininogen.^ Although the T1 kininogen gene is inducible by inflammatory cytokines, a highly homologous K kininogen gene is minimally responsive. Moreover, the basal expression of a KK/CAT construct was 5- to 7-fold higher than that of the analogous T1K/CAT construct. To examine the molecular basis of this differential regulation, a series of promoter swapping experiments was carried out. Our transfection results showed that at least two regions in the K kininogen gene are important for its high basal expression: a distal 19-bp region (C box) constituted a binding site for CCAAT/enhancer binding protein (C/EBP) family proteins and a proximal 66-bp region contained two adjacent binding sites for hepatocyte nuclear factor-3 (HNF-3). The distal HNF-3 binding site from the K kininogen promoter demonstrated a stronger affinity than that from the T1 kininogen promoter. Since C/EBP and HNF-3 are highly enriched in the liver and known to enhance transcription of liver-specific genes, differential binding affinities of these factors accounted for the higher basal expression of the K kininogen gene.^ In contrast to the K kininogen C box, the T1 kininogen C box does not bind C/EBP presumably due to their two-nucleotide divergence. This sequence divergence, however, converts it to a consensus binding sequence for two IL-6-inducible transcription factors--IL-6 response element binding protein and acute-phase response factor. To functionally determine whether C box sequences are important for their differential acute-phase response, T1 and K kininogen C boxes were swapped and analyzed after transfection into Hep3B cells. Our results showed that the T1 kininogen C box is indeed one of the IL-6 response elements in T1 kininogen promoter. Furthermore, its function can be modulated by a 5$\sp\prime$-adjacent C/EBP-binding site (B box) whose mutation significantly reduced the overall induced activity. Moreover, this B box is the target site for binding and transactivation of another IL-6 inducible transcription factor C/EBP$\delta.$ Evolutionary divergence of a few critical nucleotides can either lead to subtle changes in the binding affinities of a given transcription factor or convert a binding sequence for a constitutive factor to a site recognized by an inducible factor. (Abstract shortened by UMI.) ^

Analysis of differential gene expression in nontumorigenic glioblastoma multiforme

Loss of chromosome 10 represents the most common cytogenetic abnormality in high grade gliomas (glioblastoma multiforme). To identify genes involved in the malignant progression of human gliomas, a subtractive hybridization was performed between a tumorigenic glioblastoma cell line (LG11) and a nontumorgenic hybrid cell (LG11.3) containing an introduced chromosome 10. LG11 mRNA was subtracted from LG11.3 cDNA to produce cDNA probes enriched for sequences whose expression differs quantitatively from the parental tumorigenic cells. Both known and novel sequences were identified as a result of the subtraction. Northern blot analysis was then used to confirm differential expression of several subtracted clones. One novel clone, clone 17, identified a 2.6 kb message that showed a consistent two to four fold increase in expression in the LG11.3 nontumorigenic cells. Clone 17 (340 bp) was used successfully to screen for a near full-length version, RIG (regulated in glioma), which was 2,569 bp in size. The RIG cDNA sequence showed homology to clone 17 and to an anonymous EST (IB666), but to no previously identified genes. This screening effort also identified several independent clones representing novel sequences, most of which failed to show increased expression in the nontumorigenic GBM cells. Tissue distribution studies of RIG indicated highest levels of expression in human brain with appreciably lower levels in heart and lung. In vitro transcription and translation experiments demonstrated the ability of RIG to direct the synthesis of a 13 kD protein product. However, open reading frame analysis revealed no identify with previously described motifs or any known proteins. Using a combination of somatic cell hybrid panels and in situ hybridization, the RIG gene was mapped to chromosome 11p14-11p15. Further study of RIG and related gene products may provide insight into the negative regulation of glial oncogenesis. ^