Transcriptional regulation and functional analysis of the Wilms' tumor gene WT1

The Wilms' tumor gene, WT1, encodes a zinc finger transcription factor which functions as a tumor suppressor. Defects in the WT1 gene can result in the development of nephroblastoma. WT1 is expressed during development, primarily in the metanephric kidney, the mesothelial lining of the abdomen and thorax, and the developing gonads. WT1 expression is tightly regulated and is essential for renal development. The WT1 gene encodes a protein with a proline-rich N-terminus which functions as a transcriptional repressor and C-terminus contains 4 zinc fingers that mediate DNA binding. WT1 represses transcription from a number of growth factors and growth factor receptors. WT1 mRNA undergoes alternative splicing at two sites, resulting in 4 mRNA species and polypeptide products. Exon 5, encoding 17 amino acids is alternatively spliced, and is located between the transcriptional repression domain and the DNA binding domain. The second alternative splice is the terminal 9 nucleotides of zinc finger 3, encoding the tripeptide Lys-Thr-Ser (KTS). The presence or absence of KTS within the zinc fingers of WT1 alters DNA binding.^ I have investigated transcriptional regulation of WT1, characterizing two means of repressing WT1 transcription. I have cloned a transcriptional silencer of the WT1 promoter which is located in the third intron of the WT1 gene. The silencer is 460 bp in length and contains an Alu repeat. The silencer functions in cells of non-renal origin.^ I have found that WT1 protein can autoregulate the WT1 promoter. Using the autoregulation of the WT1 promoter as a functional assay, I have defined differential consensus DNA binding motifs of WT1 isoforms lacking and containing the KTS tripeptide insertion. With these refined consensus DNA binding motifs, I have identified two additional targets of WT1 transcriptional repression, the proto-oncogenes bcl-2 and c-myc.^ I have investigated the ability of the alternatively spliced exon 5 to influence cell growth. In cell proliferation assays, isoforms of WT1 lacking exon 5 repress cell growth. WT1 isoforms containing exon 5 fail to repress cell growth to the same extent, but alter the morphology of the cells. These experiments demonstrate that the alternative splice isoforms of WT1 have differential effects on the function of WT1. These findings suggest a role for the alternative splicing of WT1 in metanephric development. ^

Transcriptional regulation and functional analysis of the human Wilms' tumor 1 gene

The Wilms' tumor 1 gene (WT1) encodes a zinc-finger transcription factor and is expressed in urogenital, hematopoietic and other tissues. It is expressed in a temporal and spatial manner in both embryonic and adult stages. To obtain a better understanding of the biological function of WT1, we studied two aspects of WT1 regulation: one is the identification of tissue-specific cis-regulatory elements that regulate its expression, the other is the downstream genes which are modulated by WT1.^ My studies indicate that in addition to the promoter, other regulatory elements are required for the tissue specific expression of this gene. A 259-bp hematopoietic specific enhancer in intron 3 of the WT1 gene increased the transcriptional activity of the WT1 promoter by 8- to 10-fold in K562 and HL60 cells. Sequence analysis revealed both GATA and c-Myb motifs in the enhancer fragment. Mutation of the GATA motif decreased the enhancer activity by 60% in K562 cells. Electrophoretic mobility shift assays showed that both GATA-1 and GATA-2 proteins in K562 nuclear extracts bind to this motif. Cotransfection of the enhancer containing reporter construct with a GATA-1 or GATA-2 expression vector showed that both GATA-1 and GATA-2 transactivated this enhancer, increasing the CAT reporter activity 10-15 fold and 5-fold respectively. Similar analysis of the c-Myb motif by cotransfection with the enhancer CAT reporter construct and a c-Myb expression vector showed that c-Myb transactivated the enhancer by 5-fold. A DNase I-hypersensitive site has been identified in the 258 bp enhancer region. These data suggest that GATA-1 and c-Myb are responsible for the activity of this enhancer in hematopoietic cells and may bind to the enhancer in vivo. In the process of searching for cis-regulatory elements in transgenic mice, we have identified a 1.0 kb fragment that is 50 kb downstream from the promoter and is required for the central nervous system expression of WT1.^ In the search for downstream target genes of WT1, we noted that the proto-oncogene N-myc is coexpressed with the tumor suppressor gene WT1 in the developing kidney and is overexpressed in many Wilms' tumors. Sequence analysis revealed eleven consensus WT1 binding sites located in the 1 kb mouse N-myc promoter. We further showed that the N-myc promoter was down-regulated by WT1 in transient transfection assays. Electrophoretic mobility shift assays showed that oligonucleotides containing the WT1 motifs could bind WT1 protein. Furthermore, a Denys-Drash syndrome mutant of WT1, R394W, that has a mutation in the DNA binding domain, failed to repress the N-myc promoter. This suggests that the repression of the N-myc promoter is mediated by DNA binding of WT1. This finding helps to elucidate the relationship of WT1 and N-myc in tumorigenesis and renal development. ^

Long range physical map of the Wilms' tumor - aniridia region on human chromosome 11

Nephroblastoma or Wilms' tumor is a pediatric renal malignancy that is the most frequently occurring childhood solid tumor. Approximately 1-2% of children with Wilms' tumor also present with aniridia, a congenital absence of all or part of the iris of the eye. These children also have high rates of genitourinary anomalies and mental retardation resulting in what is called the WAGR (Wilms' tumor, aniridia, genitourinary anomaly, mental retardation) syndrome. Cytogenetic analysis of metaphase chromosomes from these patients revealed a consistent deletion of band P13 on chromosome 11. These observations suggest close physical linkage between the disease-related loci, and further imply that development of each phenotype results from the loss of normal gene function.^ The objective of this work is to understand the molecular events at chromosome band 11p13 that are essential to the development of sporadic Wilms' tumor and sporadic aniridia. Two human/hamster somatic cell hybrids have been used to identify sixteen independent DNA probes that map to this segment of the human genome. These newly identified DNA probes and four previously reported probes (CAT, FSHB, D11S16, and HBVIS) have been used to subdivide 11p13 into five intervals defined by overlapping constitutional deletions from several WAGR patients. A long-range physical map of 11p13 has been constructed using each of these probes in Southern blot analysis of genomic DNA after digestion with infrequently cutting restriction enzymes and pulse-field gel electrophoresis. This map, established primarily with MluI and NotI, spans approximately 13 $\times$ 10$\sp{6}$ bp and encompasses deletion and translocation breakpoints associated with genitourinary anomalies, aniridia, and sporadic Wilms' tumor. This complete physical map of human chromosome band 11p13 enables us to localize the genes for sporadic Wilms' tumor and sporadic aniridia to a small number of specific NotI fragments. ^

GENETIC EVENTS PREDISPOSING TO WILMS' TUMOR

Wilms' tumor (WT) is a childhood embryonic tumor of the kidney. In some cases, WT has been associated with a chromosome deletion in the region 11p13. The majority of WT cases, however, have normal karyotypes with no discernable deletions or rearrangements of chromosome 11.^ To study the genetic events predisposing to the development of WT, I have used a number of gene markers specific for chromosome 11. Gene probes for human catalase and apolipoprotein A1 were localized to chromosome 11 by in situ hybridization. A number of other probes previously mapped to chromosome 11 were also used. Nine WT patients who were heterozygous for at least one 11p marker were shown to lose heterozygosity in their tumor DNA. Gene dosage experiments demonstrated that two chromosomes 11 were present although loss of heterozygosity had occurred in all but two cases. By using gene probes from the short and long arms of chromosome 11, I discerned that loss of heterozygosity was due to somatic recombination in four cases, chromosome deletion in two cases, and chromosome loss and reduplication or somatic recombination in these cases. Examination of DNAs from the parents of six of these patients indicated that the alleles that were lost in tumor tissues were alleles inherited from the mother. In sporadic WT cases one would expect the loss of alleles to be random. These data suggest that the loss of alleles resulting in the development of WT is not a random event, however, the significance of this is not known. ^

DYSREGULATION OF MEOX2 FOLLOWING WT1 MUTATION IN KIDNEY DEVELOPMENT AND WILMS TUMORIGENESIS

Wilms tumor (WT) is a childhood tumor of the kidney and a productive model for understanding the role of genetic alteration and interactions in tumorigenesis. The Wilms tumor gene 1 (WT1) is a transcriptional factor and one of the few genes known to have genetic alterations in WT and has been shown be inactivated in 20% of WTs. However, the mechanisms of how WT1 mutations lead to Wilms tumorigenesis and its influence on downstream genes are unknown. Since it has been established that WT1 is a transcriptional regulator, it has been hypothesized that the loss of WT1 leads to the dysregulation of downstream genes, in turn result in the formation of WTs. To identify the dysregulated downstream genes following WT1 mutations, an Affymetrix GeneChip Human Genome Array was previously conducted to assess the differentially expressed genes in the WT1-wildtype human and WT1-mutant human WTs. Approximately 700 genes were identified as being significantly dysregulated. These genes were further prioritized based on their statistical significance, fold change, chromosomal region, spatial pattern of gene expression and known or putative cellular functions. Mesenchyme homeobox 2 (MEOX2) was one of the most significantly upregulated genes in WT1-mutant WT. MEOX2 is known to play a role in cell proliferation, apoptosis, and differentiation. In addition to its biological roles, it is expressed during early kidney development in the condensed mesenchyme similar to WT1. Furthermore, the use of the Match® web-based tool from the BIOBASE Biological Data base identified a significant predicted WT1 binding site within the first intron of MEOX2. The similarity in spatial gene expression in the developing kidney and the significant predicted WT1 binding site found in the first intron of MEOX2 lead to the development of my hypothesis that MEOX2 is upregulated via a WT1-dependent manner. Here as a part of my master’s work, I have validated the Affymetrix GeneChip Human Genome Array data using an independent set of Wilms tumors. MEOX2 remained upregulated in the mutant WT1 Wilms tumor by 41-fold. Wt1 and Meox2 gene expression were assessed in murine newborn kidney; both Wt1 and Meox2 were expressed in the condensed, undifferentiated metanephric mesenchyme. I have shown that the in vivo ablation of Wt1 during embryonic development at embryonic day (E) 13.5 resulted in the slight increase of Meox2 gene expression by two fold. In order to functionally demonstrate the effect of the loss of Wt1 on Meox2 gene expression in undifferentiated metanephric mesenchyme, I have generated a kidney mesenchymal cell line to genetically ablate Wt1 in vitro by adenoviral infection. The ablation of Wt1 in the kidney mesenchymal cell line resulted in the upregulation of Meox2 by 61-fold. Moreover, the upregulation of Meox2 resulted in the significant induction of p21 and Itgb5. In addition to the dysregulation of these genes the ablation of Wt1 in the kidney mesenchymal cells resulted in decrease in cell growth and loss of cellular adherence. However, it is uncertain whether the upregulation of Meox2 caused this particular cellular phenotype. Overall, I have demonstrated that the upregulation of Meox2 is Wt1-dependent during early kidney development.

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

Analysis of primary Wilms tumors for WT1 RNA processing alterations

Wilms tumor (WT) is an embryonal renal tumor with a heterogeneous genetic etiology that serves as a valuable model for studying tumorigenesis. Biallelic inactivation of the tumor suppressor gene WT1, a zinc-finger transcriptional regulator located at 11p13, is critical for the development of some Wilms tumors. Interestingly, WT1 genomic analysis has demonstrated mutations in less than 20% of WT cases. This suggests either other genes play a more major role in Wilms tumorigenesis or WT1 is functionally altered by mechanisms other than DNA mutation. Previous observations in rat and in WT xenograft cell lines have suggested that abnormal WT1 RNA processing (exon 6 RNA editing and aberrant exon 2 splicing, respectively) is a potential mechanism of altering WT1 function in the absence of a WT1 DNA mutation. However, the role of this abnormal RNA processing has not previously been assessed in primary Wilms tumors. ^ To test the hypothesis that abnormal WT1 RNA processing is a mechanism of WT1alteration during tumor development, WT1 RNA from 85 primary tumors was analyzed using reverse transcription and polymerase chain reaction amplification (RT-PCR). Although no evidence for WT1 RNA editing was observed, variable levels (5% to 50%) of aberrant WT1 exon 2 splicing were detected for 11 tumors in the absence of a detectable WT1 DNA mutation. Also, alteration of normal WT1 alternative splicing, observed as RNA isoform loss, was detected in five tumors with no apparent WT1 genomic alteration, although no consistent pattern of RNA isoform loss was detected. This abnormal WT1 splicing, detected by either loss of exon 2 from some of the transcripts or loss of RNA isoforms, is statistically correlated with relapse (p = 0.005). These studies demonstrate that abnormal WT1 RNA processing is not a common mechanism of abrogating normal WT1 function in primary tumors. However, in those cases in which abnormal WTI splicing is present, these data indicate that it may serve as a useful prognostic marker for relapse in WT patients. ^

Regulation of Bcl-2 expression by the product of the prostate apoptosis response gene -4, Par-4

Elevated expression levels of the bcl-2 proto-oncogene have been correlated with the appearance of androgen independence in prostate cancer. Although bcl-2 was first cloned as the t (14:18) translocation breakpoint from human follicular B cell lymphoma, the mechanism of overexpression of bcl-2 is largely undefined for advanced prostate cancer, there being no gross alterations in the gene structure. We investigated the role of the product of the prostate apoptosis response gene-4 (Par-4) and the product of the Wilms' tumor 1 gene (WT1) in the regulation of Bcl-2 expression in prostate cancer cell lines. We observed growth arrest and apoptosis, upon decreasing Bcl-2 protein and transcript in the high Bcl-2 expressing, androgen-independent prostate cancer cell lines, by all trans-retinoic acid treatment but this did not occur in the androgen-dependent cell lines expressing low levels of Bcl-2. Changes in localization of Par-4, and an induction in the expression of WT1 protein accompanied the decrease in the Bcl-2 protein and transcript following all trans-retinoic acid treatment, in the androgen-independent prostate cancer cell line. In stable clones expressing ectopic Par-4 we observed decreased Bcl-2 protein and transcript. This was accompanied by an induction in WT1 expression. Finally, we detected Par-4 and WT1 proteins binding to a previously identified WT1 binding site on the bcl-2 promoter both in vitro and in vivo leading to a decrease in transcription from the bcl-2 promoter. We conclude that Par-4 regulates Bcl-2 through a WT1 binding site on the bcl-2 promoter. ^

Application of stellar photometry to the analysis of microarray images

Improvements in the analysis of microarray images are critical for accurately quantifying gene expression levels. The acquisition of accurate spot intensities directly influences the results and interpretation of statistical analyses. This dissertation discusses the implementation of a novel approach to the analysis of cDNA microarray images. We use a stellar photometric model, the Moffat function, to quantify microarray spots from nylon microarray images. The inherent flexibility of the Moffat shape model makes it ideal for quantifying microarray spots. We apply our novel approach to a Wilms' tumor microarray study and compare our results with a fixed-circle segmentation approach for spot quantification. Our results suggest that different spot feature extraction methods can have an impact on the ability of statistical methods to identify differentially expressed genes. We also used the Moffat function to simulate a series of microarray images under various experimental conditions. These simulations were used to validate the performance of various statistical methods for identifying differentially expressed genes. Our simulation results indicate that tests taking into account the dependency between mean spot intensity and variance estimation, such as the smoothened t-test, can better identify differentially expressed genes, especially when the number of replicates and mean fold change are low. The analysis of the simulations also showed that overall, a rank sum test (Mann-Whitney) performed well at identifying differentially expressed genes. Previous work has suggested the strengths of nonparametric approaches for identifying differentially expressed genes. We also show that multivariate approaches, such as hierarchical and k-means cluster analysis along with principal components analysis, are only effective at classifying samples when replicate numbers and mean fold change are high. Finally, we show how our stellar shape model approach can be extended to the analysis of 2D-gel images by adapting the Moffat function to take into account the elliptical nature of spots in such images. Our results indicate that stellar shape models offer a previously unexplored approach for the quantification of 2D-gel spots. ^

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

K1 protein of human herpesvirus 8 suppresses lymphoma cell Fas-mediated apoptosis.

Expression of the K1 gene of human herpesvirus 8 activates nuclear factor-kappaB and induces lymph node hyperplasia and lymphomas in transgenic mice. To further delineate its role in cell survival, we determined whether K1 altered apoptosis of lymphoma cells. K1 protein is expressed in Kaposi sarcoma and primary effusion lymphoma. We retrovirally transfected BJAB lymphoma, THP-1, U937, and Kaposi sarcoma SLK cells to express K1 and a K1 mutant with the deleted immunoreceptor tyrosine-based activation motif (K1m). We challenged cells with an agonistic anti-Fas antibody, Fas ligand, irradiation, and tumor necrosis factor-related apoptosis-inducing ligand. K1 transfectants but not K1m transfectants exhibited reduced levels of apoptosis induced by the anti-Fas antibody but not apoptosis induced by the tumor necrosis factor-related apoptosis-inducing ligand or irradiation. K1 expression resulted in reduced apoptosis rates as shown in several assays. K1 induced a modest reduction in levels of Fas-associated death domain protein, and procaspase 8 recruited to the death-inducing signaling complex. Finally, K1 transfectants cleaved procaspase 8 at significantly lower rates than did K1m transfectants. K1-transfected mice, compared with vector-transfected mice, showed lower death rates after challenge with anti-Fas antibody. K1 may contribute to lymphoma development by stimulating cell survival by selectively blocking Fas-mediated apoptosis.

Activation of heat shock and antioxidant responses by the natural product celastrol: transcriptional signatures of a thiol-targeted molecule.

Stress response pathways allow cells to sense and respond to environmental changes and adverse pathophysiological states. Pharmacological modulation of cellular stress pathways has implications in the treatment of human diseases, including neurodegenerative disorders, cardiovascular disease, and cancer. The quinone methide triterpene celastrol, derived from a traditional Chinese medicinal herb, has numerous pharmacological properties, and it is a potent activator of the mammalian heat shock transcription factor HSF1. However, its mode of action and spectrum of cellular targets are poorly understood. We show here that celastrol activates Hsf1 in Saccharomyces cerevisiae at a similar effective concentration seen in mammalian cells. Transcriptional profiling revealed that celastrol treatment induces a battery of oxidant defense genes in addition to heat shock genes. Celastrol activated the yeast Yap1 oxidant defense transcription factor via the carboxy-terminal redox center that responds to electrophilic compounds. Antioxidant response genes were likewise induced in mammalian cells, demonstrating that the activation of two major cell stress pathways by celastrol is conserved. We report that celastrol's biological effects, including inhibition of glucocorticoid receptor activity, can be blocked by the addition of excess free thiol, suggesting a chemical mechanism for biological activity based on modification of key reactive thiols by this natural product.

A NEW TUMOR SUPPRESSOR GENE CANDIDATE REGULATED BY THE NON-CODING RNA PCA3 IN HUMAN PROSTATE CANCER

Prostate cancer is the second leading cause of cancer-related death and the most common non-skin cancer in men in the USA. Considerable advancements in the practice of medicine have allowed a significant improvement in the diagnosis and treatment of this disease and, in recent years, both incidence and mortality rates have been slightly declining. However, it is still estimated that 1 man in 6 will be diagnosed with prostate cancer during his lifetime, and 1 man in 35 will die of the disease. In order to identify novel strategies and effective therapeutic approaches in the fight against prostate cancer, it is imperative to improve our understanding of its complex biology since many aspects of prostate cancer initiation and progression still remain elusive. The study of tumor biomarkers, due to their specific altered expression in tumor versus normal tissue, is a valid tool for elucidating key aspects of cancer biology, and may provide important insights into the molecular mechanisms underlining the tumorigenesis process of prostate cancer. PCA3, is considered the most specific prostate cancer biomarker, however its biological role, until now, remained unknown. PCA3 is a long non-coding RNA (ncRNA) expressed from chromosome 9q21 and its study led us to the discovery of a novel human gene, PC-TSGC, transcribed from the opposite strand and in an antisense orientation to PCA3. With the work presented in this thesis, we demonstrate that PCA3 exerts a negative regulatory role over PC-TSGC, and we propose PC-TSGC to be a new tumor suppressor gene that contrasts the transformation of prostate cells by inhibiting Rho-GTPases signaling pathways. Our findings provide a biological role for PCA3 in prostate cancer and suggest a new mechanism of tumor suppressor gene inactivation mediated by non-coding RNA. Also, the characterization of PCA3 and PC-TSGC led us to propose a new molecular pathway involving both genes in the transformation process of the prostate, thus providing a new piece of the jigsaw puzzle representing the complex biology of prostate cancer.

XENOESTROGEN-SPECIFIC MECHANISMS OF DEVELOPMENTAL REPROGRAMMING CORRELATE WITH GENE EXPRESSION AND TUMOR DEVELOPMENT

Environmental exposures during sensitive windows of development can reprogram normal physiological responses and alter disease susceptibility later in life in a process known as developmental reprogramming. We have shown that neonatal exposure to the xenoestrogen diethylstilbestrol (DES) can developmentally reprogram the reproductive tract in genetically susceptible Eker rats giving rise to complete penetrance of uterine leiomyoma. Based on this, we hypothesized that xenoestrogens, including genistein (GEN) and bisphenol A (BPA), reprogram estrogen-responsive gene expression in the myometrium and promote the development of uterine leiomyoma. We proposed the mechanism that is responsible for the developmental reprogramming of gene expression was through estrogen (E2)/ xenoestrogen inducedrapid ER signaling, which modifies the histone methyltransferase Enhancer of Zeste homolog 2 (EZH2) via activation of the PI3K/AKT pathway. We further hypothesized that there is a xenostrogen-specific effect on this pathway altering patterns of histone modification, DNA methylation and gene expression. In addition to our novel finding that E2/DES-induced phosphorylation of EZH2 by AKT reduces the levels of H3K27me3 in vitro and in vivo, this work demonstrates in vivo that a brief neonatal exposure to GEN, in contrast to BPA, activates the PI3K/AKT pathway to regulate EZH2 and decreases H3K27me3 levels in the neonatal uterus. Given that H3K27me3 is a repressive mark that has been shown to result in DNA methylation and gene silencing we investigated the methylation of developmentally reprogrammed genes. In support of this evidence, we show that neonatal DES exposure in comparison to VEH, leads to hypomethylation of the promoter of a developmentally reprogrammed gene, Gria2, that become hyper-responsive to estrogen in the adult myometrium indicating vi that DES exposure alter gene expression via chromatin remodeling and loss of DNA methylation. In the adult uterus, GEN and BPA exposure developmentally reprogrammed expression of estrogen-responsive genes in a manner opposite of one another, correlating with our previous data. Furthermore, the ability of GEN and BPA to developmental reprogram gene expression correlated with tumor incidence and multiplicity. These data show that xenoestrogens have unique effects on the activation of non-genomic signaling in the developing uterus that promotes epigenetic and genetic alterations, which are predictive of developmental reprogramming and correlate with their ability to modulate hormone-dependent tumor development.

SURVIVAL PREDICTION FOR BRAIN TUMOR PATIENTS USING GENE EXPRESSION DATA

Brain tumor is one of the most aggressive types of cancer in humans, with an estimated median survival time of 12 months and only 4% of the patients surviving more than 5 years after disease diagnosis. Until recently, brain tumor prognosis has been based only on clinical information such as tumor grade and patient age, but there are reports indicating that molecular profiling of gliomas can reveal subgroups of patients with distinct survival rates. We hypothesize that coupling molecular profiling of brain tumors with clinical information might improve predictions of patient survival time and, consequently, better guide future treatment decisions. In order to evaluate this hypothesis, the general goal of this research is to build models for survival prediction of glioma patients using DNA molecular profiles (U133 Affymetrix gene expression microarrays) along with clinical information. First, a predictive Random Forest model is built for binary outcomes (i.e. short vs. long-term survival) and a small subset of genes whose expression values can be used to predict survival time is selected. Following, a new statistical methodology is developed for predicting time-to-death outcomes using Bayesian ensemble trees. Due to a large heterogeneity observed within prognostic classes obtained by the Random Forest model, prediction can be improved by relating time-to-death with gene expression profile directly. We propose a Bayesian ensemble model for survival prediction which is appropriate for high-dimensional data such as gene expression data. Our approach is based on the ensemble "sum-of-trees" model which is flexible to incorporate additive and interaction effects between genes. We specify a fully Bayesian hierarchical approach and illustrate our methodology for the CPH, Weibull, and AFT survival models. We overcome the lack of conjugacy using a latent variable formulation to model the covariate effects which decreases computation time for model fitting. Also, our proposed models provides a model-free way to select important predictive prognostic markers based on controlling false discovery rates. We compare the performance of our methods with baseline reference survival methods and apply our methodology to an unpublished data set of brain tumor survival times and gene expression data, selecting genes potentially related to the development of the disease under study. A closing discussion compares results obtained by Random Forest and Bayesian ensemble methods under the biological/clinical perspectives and highlights the statistical advantages and disadvantages of the new methodology in the context of DNA microarray data analysis.