Different approaches to identify markers for meat quality and production traits in pigs: analysis of the transcriptome in post mortem skeletal muscle and investigation of candidate genes

Due to the growing attention of consumers towards their food, improvement of quality of animal products has become one of the main focus of research. To this aim, the application of modern molecular genetics approaches has been proved extremely useful and effective. This innovative drive includes all livestock species productions, including pork. The Italian pig breeding industry is unique because needs heavy pigs slaughtered at about 160 kg for the production of high quality processed products. For this reason, it requires precise meat quality and carcass characteristics. Two aspects have been considered in this thesis: the application of the transcriptome analysis in post mortem pig muscles as a possible method to evaluate meat quality parameters related to the pre mortem status of the animals, including health, nutrition, welfare, and with potential applications for product traceability (chapters 3 and 4); the study of candidate genes for obesity related traits in order to identify markers associated with fatness in pigs that could be applied to improve carcass quality (chapters 5, 6, and 7). Chapter three addresses the first issue from a methodological point of view. When we considered this issue, it was not obvious that post mortem skeletal muscle could be useful for transcriptomic analysis. Therefore we demonstrated that the quality of RNA extracted from skeletal muscle of pigs sampled at different post mortem intervals (20 minutes, 2 hours, 6 hours, and 24 hours) is good for downstream applications. Degradation occurred starting from 48 h post mortem even if at this time it is still possible to use some RNA products. In the fourth chapter, in order to demonstrate the potential use of RNA obtained up to 24 hours post mortem, we present the results of RNA analysis with the Affymetrix microarray platform that made it possible to assess the level of expression of more of 24000 mRNAs. We did not identify any significant differences between the different post mortem times suggesting that this technique could be applied to retrieve information coming from the transcriptome of skeletal muscle samples not collected just after slaughtering. This study represents the first contribution of this kind applied to pork. In the fifth chapter, we investigated as candidate for fat deposition the TBC1D1 [TBC1 (tre-2/USP6, BUB2, cdc16) gene. This gene is involved in mechanisms regulating energy homeostasis in skeletal muscle and is associated with predisposition to obesity in humans. By resequencing a fragment of the TBC1D1 gene we identified three synonymous mutations localized in exon 2 (g.40A>G, g.151C>T, and g.172T>C) and 2 polymorphisms localized in intron 2 (g.219G>A and g.252G>A). One of these polymorphisms (g.219G>A) was genotyped by high resolution melting (HRM) analysis and PCR-RFLP. Moreover, this gene sequence was mapped by radiation hybrid analysis on porcine chromosome 8. The association study was conducted in 756 performance tested pigs of Italian Large White and Italian Duroc breeds. Significant results were obtained for lean meat content, back fat thickness, visible intermuscular fat and ham weight. In chapter six, a second candidate gene (tribbles homolog 3, TRIB3) is analyzed in a study of association with carcass and meat quality traits. The TRIB3 gene is involved in energy metabolism of skeletal muscle and plays a role as suppressor of adipocyte differentiation. We identified two polymorphisms in the first coding exon of the porcine TRIB3 gene, one is a synonymous SNP (c.132T> C), a second is a missense mutation (c.146C> T, p.P49L). The two polymorphisms appear to be in complete linkage disequilibrium between and within breeds. The in silico analysis of the p.P49L substitution suggests that it might have a functional effect. The association study in about 650 pigs indicates that this marker is associated with back fat thickness in Italian Large White and Italian Duroc breeds in two different experimental designs. This polymorphisms is also associated with lactate content of muscle semimembranosus in Italian Large White pigs. Expression analysis indicated that this gene is transcribed in skeletal muscle and adipose tissue as well as in other tissues. In the seventh chapter, we reported the genotyping results for of 677 SNPs in extreme divergent groups of pigs chosen according to the extreme estimated breeding values for back fat thickness. SNPs were identified by resequencing, literature mining and in silico database mining. analysis, data reported in the literature of 60 candidates genes for obesity. Genotyping was carried out using the GoldenGate (Illumina) platform. Of the analyzed SNPs more that 300 were polymorphic in the genotyped population and had minor allele frequency (MAF) >0.05. Of these SNPs, 65 were associated (P<0.10) with back fat thickness. One of the most significant gene marker was the same TBC1D1 SNPs reported in chapter 5, confirming the role of this gene in fat deposition in pig. These results could be important to better define the pig as a model for human obesity other than for marker assisted selection to improve carcass characteristics.

A novel specific genetic translocation in epithelioid hemangioendotelioma, showing a fusion of the WWTR1 and CAMTA1 genes, supports the monoclonality of multifocal epithelioid hemangioendotelioma.

Like other vascular tumors, epithelioid hemangioendothelioma (EHE) is multifocal in approximately 50% of cases, and it is unclear whether the separate lesions represent multifocal disease or metastases. We hypothesized that the identification of an identical WWTR1-CAMTA1 rearrangement in different EHEs from the same patient supports the monoclonal origin of EHE. To test our hypothesis, we undertook a molecular analysis of two multicentric EHEs of the liver, including separate tumor samples from each patient. Matherial and Methods: We retrieved two cases of EHE with available tissue for molecular analysis. In both cases, fluorescence in situ hybridization (FISH) was performed to identify the presence of the WWTR1-CAMTA1 rearrangement to confirm the histologic diagnosis of EHE, as previously described. The reverse transcription-polymerase chain reaction (RT-PCR) products were analyzed by electrophoresis and the RT-PCR–amplified products were sequenced using the Sanger method. Results: FISH analysis revealed signal abnormalities in both WWTR1 and CAMTA1. Combined results confirmed the presence of the t(1;3)(1p36.23;3q25.1) translocation in both cases of EHE. Using RT-PCR analysis, we found that the size of the rearranged bands was identical in the different tumors from each patient. The sequence of the fusion gene confirmed a different WWTR1-CAMTA1 rearrangement in each patient, but an identical WWTR1-CAMTA1 rearrangement in the different lesions from each patient. Discussion: Because of its generally indolent clinical course, EHE is commonly classified as a multifocal, rather than metastatic, disease. In this study, we examined two cases of multifocal liver EHE and found an identical WWTR1-CAMTA1 rearrangement in each lesion from the same patient, but not between the two patients. These findings suggest that multifocal EHE arises from metastasis of the same neoplastic clone rather than from the simultaneous formation of multiple neoplastic clones, which supports the monoclonal origin of multifocal EHE.

IGF system, CD99 and tumor-specific chromosomal translocations: evaluation of their cross-talk and definition of their role in Ewing sarcoma and prostate cancer

Aberrant expression of ETS transcription factors, including FLI1 and ERG, due to chromosomal translocations has been described as a driver event in initiation and progression of different tumors. In this study, the impact of prostate cancer (PCa) fusion gene TMPRSS2-ERG was evaluated on components of the insulin-like growth factor (IGF) system and the CD99 molecule, two well documented targets of EWS-FLI1, the hallmark of Ewing sarcoma (ES). The aim of this study was to identify common or distinctive ETS-related mechanisms which could be exploited at biological and clinical level. The results demonstrate that IGF-1R represents a common target of ETS rearrangements as ERG and FLI1 bind IGF-1R gene promoter and their modulation causes alteration in IGF-1R protein levels. At clinical level, this mechanism provides basis for a more rationale use of anti-IGF-1R inhibitors as PCa cells expressing the fusion gene better respond to anti-IGF-1R agents. EWS-FLI1/IGF-1R axis provides rationale for combination of anti-IGF-1R agents with trabectedin, an alkylator agent causing enhanced EWS-FLI1 occupancy on the IGF-1R promoter. TMPRSS2-ERG also influences prognosis relevance of IGF system as high IGF-1R correlates with a better biochemical progression free survival (BPFS) in PCa patients negative for the fusion gene while marginal or no association was found in the total cases or TMPRSS2-ERG-positive cases, respectively. This study indicates CD99 is differentially regulated between ETS-related tumors as CD99 is not a target of ERG. In PCa, CD99 did not show differential expression between TMPRSS2-ERG-positive and –negative cells. A direct correlation was anyway found between ERG and CD99 proteins both in vitro and in patients putatively suggesting that ERG target genes comprehend regulators of CD99. Despite a little trend suggesting a correlation between CD99 expression and a better BPFS, no clinical relevance for CD99 was found in the field of prognostic biomarkers.

Die Interaktion von VHL und PTEN bei der Tumorinitiation und Progression

Das VHL-Syndrom umfasst Erkrankungen, die mit einem Funktionsverlust von VHL einhergehen. Das Tumorspektrum umfasst retinale und zerebrale Hämangioblastome, Nierenzysten und klarzellige Nierenkarzinome, Zysten und Tumore des Pankreas, Phäochromocytome, Adenome der Hoden und Tumore des Mittelohrs. Obwohl aufgrund klinischer Studien bekannt ist, welche VHL-Mutation mit welchen Neoplasien assoziiert werden können, konnte bisher kein VHL-Mausmodell das Krankheitsbild des VHL-Syndroms widerspiegeln. Daher ist vermutlich eine zusätzliche Fehlregulation weiterer Gene nötig ist, um die Tumorgenese in den verschiedenen Geweben zu induzieren. In mehreren klarzelligen Nierenkarzinomen konnte bereits eine PTEN-Defizienz nachgewiesen werden, der Verlust von PTEN wird außerdem auch mit der Tumorgenese von Phäochromocytomen assoziiert. Möglicherweise wirken VHL und PTEN also in der Tumorsuppression in der Niere und der Nebenniere zusammen.rnIm Rahmen dieser Arbeit konnte erstmals eine VHL-vermittelte Stabilisierung der PTEN-Konzentration sowohl in embryonalen als auch in Tumor-Zellen der Niere nachgewiesen werden. Die Analyse des Regulationsmechanismus ergab erstens eine Hypoxie-abhängige Abnahme der Transkription von PTEN. Des Weiteren konnte eine VHL-vermittelte Ubiquitinylierung von NEDD4-1, welches als E3-Ligase von PTEN dessen Degradation und Kerntransport reguliert, ermittelt werden. rnIn Nierenkarzinom-Zellen wurde weiterhin eine VHL- bzw. PTEN-Restitution induziert, um die Auswirkungen der beiden Tumorsuppressoren auf das Zellverhalten in vitro und in vivo zu untersuchen. Sowohl VHL als auch PTEN hatten dieselben Effekte lediglich in unterschiedlicher Intensität auf das Verhalten der Zellen. So konnte VHL- und PTEN-abhängig eine Verstärkung der Adhäsion, eine Inhibierung der Migration und eine Verminderung der Überlebens- und Metastasierungsfähigkeit nachgewiesen werden. Des Weiteren wurden Mausmodelle mit einem ubiquitären, heterozygoten Pten-Verlust generiert, die teilweise eine zusätzliche Haploinsuffizienz von Vhl bzw. eine heterozygote VHL Typ II-Mutation (V2B oder V2C) trugen. Sporadisch entwickelten diese Mäuse Vhl-abhängig Lebertumore und Pten-abhängig Lymphome und Ovarialkarzinome. Einige Mäuse mit einer kombinierten Vhl- und Pten-Defizienz bildeten zusätzlich Nierenzysten aus, die teilweise das gesamte Volumen der Niere einnahmen. Besonders häufig entstanden in Pten-haploinsuffizienten Mäusen Phäochromocytome, die durch eine zusätzliche V2B- oder V2C-Mutation in gleichaltrigen Mäusen deutlich weiterentwickelt waren. Demnach induziert erst der gemeinsame Verlust von Vhl und Pten die Bildung von Nierenzysten und Phäochromocytomen, welche dem Krankheitsbild des VHL-Syndroms zugeordnet werden.rnDie Untersuchungen innerhalb dieser Arbeit zeigen erstmalig die Interaktion und Kooperation von VHL und PTEN in der Tumorsuppression. Die Resultate bieten außerdem die Grundlage für weitere Analysen der Auswirkung der VHL-vermittelten PTEN-Stabilisierung und für detailliertere Untersuchungen der durch die kombinierte Vhl- und Pten-Defizienz induzierten Neoplasien der Niere und der Nebennieren-Tumore in in vivo Mausmodellen.rn

Inactivation of the hypermethylated in cancer 1 tumour suppressor--not just a question of promoter hypermethylation?

The chromosomal region 17p13.3 is frequently deleted or epigenetically silenced in a variety of human cancers. It includes the hypermethylated in cancer 1 (HIC1) gene placed telomerically to the p53 tumour suppressor gene. HIC1 encodes a transcriptional repressor, and its targets identified to date are genes involved in proliferation, tumour growth and angiogenesis. In addition, HIC1 functionally cooperates with p53 to suppress cancer development. Frequent allelic loss at position 17p13.1 in human cancers often points to mutations of the tumour suppressor p53. However, in a variety of cancer types, allelic loss of the short arm of chromosome 17 may hit regions distal to p53 and, interestingly, without leading to p53 mutations. Furthermore, the neighbouring region 17p13.3 often shows loss of heterozygosity or DNA hypermethylation in various types of solid tumours and leukaemias. In line with this concept, Wales et al. described a new potential tumour suppressor in this region and named it hypermethylated in cancer 1 (HIC1). Further, it was shown that in the majority of cases hypermethylation of this chromosomal region leads to epigenetic inactivation of HIC1. A role for HIC1 in tumour development is further supported by a mouse model, since various spontaneous, age- and gender-specific malignant tumours occur in heterozygous Hic1+/- knockout mice. Furthermore, exogenously delivered HIC1 leads to a significant decrease in clonogenic survival in cancer cell lines. This review highlights the role of HIC1 inactivation in solid tumours and particularly in leukaemia development.

Disabling c-Myc in childhood medulloblastoma and atypical teratoid/rhabdoid tumor cells by the potent G-quadruplex interactive agent S2T1-6OTD

We investigated here the effects of S2T1-6OTD, a novel telomestatin derivative that is synthesized to target G-quadruplex-forming DNA sequences, on a representative panel of human medulloblastoma (MB) and atypical teratoid/rhabdoid (AT/RT) childhood brain cancer cell lines. S2T1-6OTD proved to be a potent c-Myc inhibitor through its high-affinity physical interaction with the G-quadruplex structure in the c-Myc promoter. Treatment with S2T1-6OTD reduced the mRNA and protein expressions of c-Myc and hTERT, which is transcriptionally regulated by c-Myc, and decreased the activities of both genes. In remarkable contrast to control cells, short-term (72-hour) treatment with S2T1-6OTD resulted in a dose- and time-dependent antiproliferative effect in all MB and AT/RT brain tumor cell lines tested (IC(50), 0.25-0.39 micromol/L). Under conditions where inhibition of both proliferation and c-Myc activity was observed, S2T1-6OTD treatment decreased the protein expression of the cell cycle activator cyclin-dependent kinase 2 and induced cell cycle arrest. Long-term treatment (5 weeks) with nontoxic concentrations of S2T1-6OTD resulted in a time-dependent (mainly c-Myc-dependent) telomere shortening. This was accompanied by cell growth arrest starting on day 28 followed by cell senescence and induction of apoptosis on day 35 in all of the five cell lines investigated. On in vivo animal testing, S2T1-6OTD may well represent a novel therapeutic strategy for childhood brain tumors.

Genetic evidence supporting the association of protease and protease inhibitor genes with inflammatory bowel disease: a systematic review

As part of the European research consortium IBDase, we addressed the role of proteases and protease inhibitors (P/PIs) in inflammatory bowel disease (IBD), characterized by chronic mucosal inflammation of the gastrointestinal tract, which affects 2.2 million people in Europe and 1.4 million people in North America. We systematically reviewed all published genetic studies on populations of European ancestry (67 studies on Crohn's disease [CD] and 37 studies on ulcerative colitis [UC]) to identify critical genomic regions associated with IBD. We developed a computer algorithm to map the 807 P/PI genes with exact genomic locations listed in the MEROPS database of peptidases onto these critical regions and to rank P/PI genes according to the accumulated evidence for their association with CD and UC. 82 P/PI genes (75 coding for proteases and 7 coding for protease inhibitors) were retained for CD based on the accumulated evidence. The cylindromatosis/turban tumor syndrome gene (CYLD) on chromosome 16 ranked highest, followed by acylaminoacyl-peptidase (APEH), dystroglycan (DAG1), macrophage-stimulating protein (MST1) and ubiquitin-specific peptidase 4 (USP4), all located on chromosome 3. For UC, 18 P/PI genes were retained (14 proteases and 4 protease inhibitors), with a considerably lower amount of accumulated evidence. The ranking of P/PI genes as established in this systematic review is currently used to guide validation studies of candidate P/PI genes, and their functional characterization in interdisciplinary mechanistic studies in vitro and in vivo as part of IBDase. The approach used here overcomes some of the problems encountered when subjectively selecting genes for further evaluation and could be applied to any complex disease and gene family.

Tenascin-C downregulates wnt inhibitor dickkopf-1, promoting tumorigenesis in a neuroendocrine tumor model

The extracellular matrix molecule tenascin-C (TNC) is a major component of the cancer-specific matrix, and high TNC expression is linked to poor prognosis in several cancers. To provide a comprehensive understanding of TNC's functions in cancer, we established an immune-competent transgenic mouse model of pancreatic β-cell carcinogenesis with varying levels of TNC expression and compared stochastic neuroendocrine tumor formation in abundance or absence of TNC. We show that TNC promotes tumor cell survival, the angiogenic switch, more and leaky vessels, carcinoma progression, and lung micrometastasis. TNC downregulates Dickkopf-1 (DKK1) promoter activity through the blocking of actin stress fiber formation, activates Wnt signaling, and induces Wnt target genes in tumor and endothelial cells. Our results implicate DKK1 downregulation as an important mechanism underlying TNC-enhanced tumor progression through the provision of a proangiogenic tumor microenvironment.

Loss of Raf-1 kinase inhibitor protein (RKIP) is strongly associated with high-grade tumor budding and correlates with an aggressive phenotype in pancreatic ductal adenocarcinoma (PDAC)

BACKGROUND Raf-1 kinase inhibitor protein (RKIP) has emerged as a significant metastatic suppressor in a variety of human cancers and is known to inhibit Ras/Raf/MEK/ERK signaling. By suppressing the activation of the NFkB/SNAIL circuit, RKIP can regulate the induction of epithelial-mesenchymal transition (EMT). The aim of this study was to evaluate RKIP expression and to determine its association with clinicopathological features, including EMT in form of tumor budding in pancreatic ductal adenocarcinoma (PDAC). METHODS Staining for RKIP was performed on a multipunch Tissue Microarray (TMA) of 114 well-characterized PDACs with clinico-pathological, follow-up and adjuvant therapy information. RKIP-expression was assessed separately in the main tumor body and in the tumor buds. Another 3 TMAs containing normal pancreatic tissue, precursor lesions (Pancreatic Intraepithelial Neoplasia, PanINs) and matched lymph node metastases were stained in parallel. Cut-off values were calculated by receiver operating characteristic (ROC) curve analysis. RESULTS We found a significant progressive loss of RKIP expression between normal pancreatic ductal epithelia (average: 74%), precursor lesions (PanINs; average: 37%), PDAC (average 20%) and lymph node metastases (average 8%, p<0.0001). RKIP expression was significantly lower in tumor buds (average: 6%) compared to the main tumor body (average 20%; p<0.005). RKIP loss in the tumor body was marginally associated with advanced T-stage (p=0.0599) as well as high-grade peritumoral (p=0.0048) and intratumoral budding (p=0.0373). RKIP loss in the buds showed a clear association with advanced T stage (p=0.0089). CONCLUSIONS The progressive loss of RKIP seems to play a major role in the neoplastic transformation of pancreas, correlates with aggressive features in PDAC and is associated with the presence of EMT in form of tumor budding.

Early gene expression analysis in 9L orthotopic tumor-bearing rats identifies immune modulation in molecular response to synchrotron microbeam radiation therapy

Synchrotron Microbeam Radiation Therapy (MRT) relies on the spatial fractionation of the synchrotron photon beam into parallel micro-beams applying several hundred of grays in their paths. Several works have reported the therapeutic interest of the radiotherapy modality at preclinical level, but biological mechanisms responsible for the described efficacy are not fully understood to date. The aim of this study was to identify the early transcriptomic responses of normal brain and glioma tissue in rats after MRT irradiation (400Gy). The transcriptomic analysis of similarly irradiated normal brain and tumor tissues was performed 6 hours after irradiation of 9 L orthotopically tumor-bearing rats. Pangenomic analysis revealed 1012 overexpressed and 497 repressed genes in the irradiated contralateral normal tissue and 344 induced and 210 repressed genes in tumor tissue. These genes were grouped in a total of 135 canonical pathways. More than half were common to both tissues with a predominance for immunity or inflammation (64 and 67% of genes for normal and tumor tissues, respectively). Several pathways involving HMGB1, toll-like receptors, C-type lectins and CD36 may serve as a link between biochemical changes triggered by irradiation and inflammation and immunological challenge. Most immune cell populations were involved: macrophages, dendritic cells, natural killer, T and B lymphocytes. Among them, our results highlighted the involvement of Th17 cell population, recently described in tumor. The immune response was regulated by a large network of mediators comprising growth factors, cytokines, lymphokines. In conclusion, early response to MRT is mainly based on inflammation and immunity which appear therefore as major contributors to MRT efficacy.

Characterization of the 9L gliosarcoma implanted in the Fischer rat: an orthotopic model for a grade IV brain tumor.

Among rodent models for brain tumors, the 9L gliosarcoma is one of the most widely used. Our 9L-European Synchrotron Radiation Facility (ESRF) model was developed from cells acquired at the Brookhaven National Laboratory (NY, USA) in 1997 and implanted in the right caudate nucleus of syngeneic Fisher rats. It has been largely used by the user community of the ESRF during the last decade, for imaging, radiotherapy, and chemotherapy, including innovative treatments based on particular irradiation techniques and/or use of new drugs. This work presents a detailed study of its characteristics, assessed by magnetic resonance imaging (MRI), histology, immunohistochemistry, and cytogenetic analysis. The data used for this work were from rats sampled in six experiments carried out over a 3-year period in our lab (total number of rats = 142). The 9L-ESRF tumors were induced by a stereotactic inoculation of 10(4) 9L cells in the right caudate nucleus of the brain. The assessment of vascular parameters was performed by MRI (blood volume fraction and vascular size index) and by immunostaining of vessels (rat endothelial cell antigen-1 and type IV collagen). Immunohistochemistry and regular histology were used to describe features such as tumor cell infiltration, necrosis area, nuclear pleomorphism, cellularity, mitotic characteristics, leukocytic infiltration, proliferation, and inflammation. Moreover, for each of the six experiments, the survival of the animals was assessed and related to the tumor growth observed by MRI or histology. Additionally, the cytogenetic status of the 9L cells used at ESRF lab was investigated by comparative genomics hybridization analysis. Finally, the response of the 9L-ESRF tumor to radiotherapy was estimated by plotting the survival curves after irradiation. The median survival time of 9L-ESRF tumor-bearing rats was highly reproducible (19-20 days). The 9L-ESRF tumors presented a quasi-exponential growth, were highly vascularized with a high cellular density and a high proliferative index, accompanied by signs of inflammatory responses. We also report an infiltrative pattern which is poorly observed on conventional 9 L tumor. The 9L-ESRF cells presented some cytogenetic specificities such as altered regions including CDK4, CDKN2A, CDKN2B, and MDM2 genes. Finally, the lifespan of 9L-ESRF tumor-bearing rats was enhanced up to 28, 35, and 45 days for single doses of 10, 20, and 2 × 20 Gy, respectively. First, this report describes an animal model that is used worldwide. Second, we describe few features typical of our model if compared to other 9L models worldwide. Altogether, the 9L-ESRF tumor model presents characteristics close to the human high-grade gliomas such as high proliferative capability, high vascularization and a high infiltrative pattern. Its response to radiotherapy demonstrates its potential as a tool for innovative radiotherapy protocols.

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.

Generation and characterization of antigenic variants induced by exposure of tumor cells to UV radiation in vitro

Antigenic changes present in nonantigenic tumor cells exposed to UV radiation (UV) in vitro were investigated by addressing the following questions: (1) Are antigenic variants (AV) produced that are rejected in normal but not immunosuppressed mice? (2) Does generation of AV depend upon intrinsic properties of the cells exposed or result from the action of UV? (3) Is antigenic modification induced by UV due to increased histocompatibility antigen expression? (4) Do AV crossreact immunologically with parental tumor or with other AV? and (5) Is the UV-associated common antigen expressed on UV-induced tumors present on UV-irradiated tumor cells? AV were generated at different frequencies following in vitro UV irradiation of a spontaneous murine fibrosarcoma (51% of cell lines tested), a murine melanoma (56%), and two melanoma clones (100% and 11%). This indicated that the percentage of AV produced is an intrinsic property of the cell line exposed. The increased antigenicity did not correlate with an increased expression of class I histocompatibility antigens. Immunological experiments demonstrated that the AV and parental cells shared a determinant that was susceptible to immune recognition, but incapable of inducing immunity. In contrast, the AV were noncrossreactive, suggesting that variant-specific antigens were also expressed. Finally, the AV were recognized by UV-induced suppressor cells, indicating that the UV-associated common antigen expressed by UV-induced tumors was also present. This investigation provides new information on the susceptibility of tumors to antigenic modification by UV and on the relationship between tumor antigens and neoplastic transformation. Furthermore, it suggests an immunological approach for cancer therapy. ^

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