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

cDNA cloning and expression of a novel cell surface-expressed \b heparan sulfate/heparin \b interacting \b protein (HIP) from human cell lines and functional studies of a synthetic HIP peptide

Heparan sulfate proteoglycans and their corresponding binding sites have been suggested to play an important role during the initial attachment of blastocysts to uterine epithelium and human trophoblastic cell lines to uterine epithelial cell lines. Previous studies on RL95 cells, a human uterine epithelial cell line, characterized a single class of cell surface heparin/heparan sulfate (HP/HS)-binding sites. Three major HP/HS-binding peptide fragments were isolated from RL95 cell surfaces by tryptic digestion and partial amino-terminal amino acid sequence from each peptide fragment was obtained. In the current study, using the approaches of reverse transcription-polymerase chain reaction and cDNA library screening, a novel cell surface $\rm\underline{H}$P/HS $\rm\underline{i}$nteracting $\rm\underline{p}$rotein (HIP) has been isolated from RL95 cells. The full-length cDNA of HIP encodes a protein of 259 amino acids with a calculated molecular weight of 17,754 Da and pI of 11.75. Transfection of HIP cDNA into NIH-3T3 cells demonstrated cell surface expression and a size similar to that of HIP expressed by human cells. Predicted amino acid sequence indicates that HIP lacks a membrane spanning region and has no consensus sites for glycosylation. Northern blot analysis detected a single transcript of 1.3 kb in both total RNA and poly(A$\sp+$) RNA. Examination of human cell lines and normal tissues using both Northern blot and Western blot analysis revealed that HIP is differentially expressed in a variety of human cell lines and normal tissues, but absent in some cell lines examined. HIP has about 80% homology, at the level of both mRNA and protein, to a rodent protein, designated as ribosomal protein L29. Thus, members of the L29 family may be displayed on cell surfaces where they participate in HP/HS binding events. Studies on a synthetic peptide derived from HIP demonstrate that HIP peptide binds HS/HP with high selectivity and has high affinity (Kd = 10 nM) for a subset of polysaccharides found in commercial HIP preparations. Moreover, HIP peptide also binds certain forms of cell surface, but not secreted or intracellular. HS expressed by RL95 and JAR cells. This peptide supports the attachment of several human trophoblastic cell lines and a variety of mammalian adherent cell lines in a HS-dependent fashion. Furthermore, studies on the subset of HP specifically recognized by HIP peptide indicate that this high-affinity HP (HA-HP) has a larger median MW and a greater negative charge density than bulk HP. The minimum size of oligosaccharide required to bind to HIP peptide with high affinity is a septa- or octasaccharide. HA-HP also quantitatively binds to antithrombin-III (AT-III) with high affinity, indicating that HIP peptide and AT-III may recognize the same or similar oligosaccharide structure(s). Furthermore, HIP peptide antagonizes HP action and promotes blood coagulation in both factor Xa- and thrombin-dependent assays. Finally, HA-HP recognized by HP peptide is highly enriched with anticoagulant activity relative to bulk HP. Collectively, these results demonstrate that HIP may play a role in the HP/HS-involved cell-cell and cell-matrix interactions and recognizes a motif in HP similar or identical to that recognized by AT-III and therefore, may modulate blood coagulation. ^

Elevated albumin in retinas of monkeys with experimental glaucoma.

PURPOSE: To establish the identity of a prominent protein, approximately 70 kDa, that is markedly increased in the retina of monkeys with experimental glaucoma compared with the fellow control retina, the relationship to glaucoma severity, and its localization in the retina. METHODS: Retinal extracts were subjected to 2-D gel electrophoresis to identify differentially expressed proteins. Purified peptides from the abundant 70 kDa protein were analyzed and identified by liquid chromatography/mass spectrometry/mass spectrometry (LC/MS/MS) separation, and collision-induced dissociation sequencing. Protein identity was performed on MASCOT (Matrix Science, Boston, MA) and confirmed by Western blot. The relationship between the increase in this protein and glaucoma severity was investigated by regression analyses. Protein localization in retina was evaluated by immunohistochemistry with confocal imaging. RESULTS: The abundant protein was identified as Macaca mulatta serum albumin precursor (67 kDa) from eight non-overlapping proteolytic fragments, and the identity was confirmed by Western blot. The average increase in retinal albumin content was 2.3 fold (P = 0.015). In glaucoma eyes, albumin was localized to some neurons of the inner nuclear layer, in the inner plexiform layer, and along the vitreal surface, but it was only found in blood vessels in control retinas. CONCLUSIONS: Albumin is the abundant protein found in the glaucomatous monkey retinas. The increased albumin is primarily localized to the inner retina where oxidative damage associated with experimental glaucoma is known to be prominent. Since albumin is a major antioxidant, the increase of albumin in the retinas of eyes with experimental glaucoma may serve to protect the retina against oxidative damage.

LDOC1, A NOVEL BIOMARKER OF PROGNOSIS IN CHRONIC LYMPHOCYTIC LEUKEMIA

In chronic lymphocytic leukemia (CLL), one of the best predictors of outcome is the somatic mutation status of the immunoglobulin heavy chain variable region (IGHV) genes. Patients whose CLL cells have unmutated IGHV genes have a median survival of 8 years; those with mutated IGHV genes have a median survival of 25 years. To identify new prognostic biomarkers and molecular targets for therapy in untreated CLL patients, we reanalyzed the raw data from four published gene expression profiling microarray studies. Of 88 candidate biomarkers associated with IGHV somatic mutation status, we identified LDOC1 (Leucine Zipper, Down-regulated in Cancer 1), as one of the most significantly differentially expressed genes that distinguished mutated from unmutated CLL cases. LDOC1 is a putative transcription factor of unknown function in B-cell development and CLL pathophysiology. Using a highly sensitive quantitative RT-PCR (QRT-PCR) assay, we confirmed that LDOC1 mRNA was dramatically down-regulated in mutated compared to unmutated CLL cases. Expression of LDOC1 mRNA was also vii strongly associated with other markers of poor prognosis, including ZAP70 protein and cytogenetic abnormalities of poor prognosis (deletions of chromosomes 6q21, 11q23, and 17p13.1, and trisomy 12). CLL cases positive for LDOC1 mRNA had significantly shorter overall survival than negative cases. Moreover, in a multivariate model, LDOC1 mRNA expression predicted overall survival better than IGHV mutation status or ZAP70 protein, among the best markers of prognosis in CLL. We also discovered LDOC1S, a new LDOC1 splice variant. Using isoform-specific QRT-PCR assays that we developed, we found that both isoforms were expressed in normal B cells (naïve > memory), unmutated CLL cells, and in B-cell non-Hodgkin lymphomas with unmutated IGHV genes. To investigate pathways in which LDOC1 is involved, we knocked down LDOC1 in HeLa cells and performed global gene expression profiling. GFI1 (Growth Factor-Independent 1) emerged as a significantly up-regulated gene in both HeLa cells and CLL cells that expressed high levels of LDOC1. GFI1 oncoprotein is implicated in hematopoietic stem cell maintenance, lymphocyte development, and lymphomagenesis. Our findings indicate that LDOC1 mRNA is an excellent biomarker of overall survival in CLL, and may contribute to B-cell differentiation and malignant transformation.

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.

Differential gene expression in multistep tumorigenesis using an {\it in vitro\/} human cell model

Using a human terato-carcinoma cell line, PA-1, the functional role of the oncogenes and tumor suppressor gene involved in the multistep process of carcinogenesis have been analyzed. The expression of AP-2 was strongly correlated with the susceptibility to ras transformation. The differential responsiveness to growth factors between stage 1 ras resistant cells and stage 2 ras susceptible cells was observed, indicating that the ability of stage 2 cells to respond to the mutated ras oncogenes in transformation correlated with the ability to be stimulated by certain growth factors. Using differential screening of cDNA libraries, a number of differentially expressed cDNA clones was isolated. One of those, clone 12, is overexpressed in ras transformed stage 3 cells. The amino acid sequence of clone 12 is almost identical to a mouse LLrep3 gene that was growth-regulated, and 78% similar to a yeast ribosomal protein S4. These results suggest that the S4 gene may be involved in regulation of growth. Clone 9 is expressed in stage 1 ras resistant cells (3.5-kb and 3.0-kb transcripts) but the expression of this clone in stage 2 ras susceptible cells and stage 3 ras-transformed cells is greatly diminished. The expression of this cDNA clone was increased to at least five fold in ras resistant cells and nontumorigenic hybrids treated with retinoic acid but not increased in retinoic acid treated ras susceptible cells, ras transformed cells and the tumorigenic segregants. Partial sequence of this clone showed no homology to the sequences in Genbank. These findings suggest that clone 9 could be a suppressor gene or the genes that are involved in the biochemical pathway of tumor suppression or neurogenic differentiation. The apparent pleiotropic effect of the loss of this suppressor gene function support Harris' proposal that tumor suppressor genes regulate differentiation. The tumor suppressor gene may act as negative regulator of tumor growth by controlling gene expression in differentiation. ^

The molecular cloning and characterization of human heparanase cDNA and the immunochemical localization of heparanase in metastatic melanomas

Heparanase, an endo-$\beta$-D-glucuronidase, has been associated with melanoma metastasis. Polyclonal antibodies directed against the murine N-terminal heparanase peptide detected a M$\sb{\rm r}\sim 97,000$ protein upon SDS-polyacrylamide gel electrophoresis of mouse melanoma and human melanoma cell lysates. In an indirect immunocytochemical study, metastatic human A375-SM and mouse B16-BL6 melanoma cells were stained with the anti-heparanase antibodies. Heparanase antigen was localized in the cytoplasm of permeabilized melanoma cells as well as at the cell surface of unpermeabilized cells. Immunohistochemical staining of frozen sections from syngeneic mouse organs containing micrometastases of B16-BL6 melanoma demonstrated heparanase localized in metastatic melanoma cells, but not in adjacent normal tissues. Similar studies using frozen sections of malignant melanomas resected from patients indicated that heparanase is localized in invading melanoma cells, but not in adjacent connective tissues.^ Monoclonal antibodies directed against murine heparanase were developed and characterized. Monoclonal antibody 10E5, an IgM, precipitated and inhibitated the enzymatic activity of heparanase. A 2.6 kb cDNA was isolated from a human melanoma $\lambda$gt11 cDNA library using the monoclonal antibody 10E5. Heparan sulfate cleavage activity was detected in the lysogen lysates from E. Coli Y1089 infected with the $\lambda$gt11 cDNA and this activity was inhibited in the presence of 10-fold excess of heparin, a potent inhibitor of heparanase. The nucleotide sequence of the cDNA was determined and insignificant homology was found with the gene sequences currently known. The cDNA hybridized to a 3.2-3.4 kb mRNA in human A375 melanoma, WI-38 fibroblast, and THP-1 leukemia cells using Northern blots.^ Heparanase expression was examined using Western and Northern blots. In comparison to human A375-P melanoma cells, the quantity of 97,000 protein recognized by the polyclonal anti-heparanase antibodies doubled in the metastatic variant A375-SM cells and the quantity of 3.2-3.4 kb mRNA doubled in A375MetMix, a metastatic variant similar to A375-SM cells. In B16 murine melanoma cell, the intensity of the 97,000 protein increased more than 2 times comparing with B16-F1 cells. The extent in the increase of the protein and the mRNA levels is comparable to the change of heparanase activity observed in those cells.^ In summary, the studies suggest that (a) the N-terminus of the heparanase molecule in mouse and human is antigenically related; (b) heparanase antigens are localized at the cell surface and in the cytoplasm of metastatic human and mouse melanoma cells; (c) heparanase antigens are localized in invasive and metastatic murine and human melanomas in vivo, but not in adjacent normal tissues; (d) heparanase molecule appeared to be differentially expressed at the transcriptional as well as at the translational level; and (e) the size of human heparanase mRNA is 3.2-3.4 kilobase. ^

Studies on the adhesion properties of the liver-metastatic murine RAW117 large-cell lymphoma cells: Roles of $\beta\sb3$ integrin

Metastasis is the major cause of death in cancer patients. Since many cancers show organ-preference of metastasis, elucidation of the underlying mechanisms of metastasis will benefit diagnosis or treatment of metastatic diseases. Adhesion mechanisms are thought to be involved in organ-preference of metastasis, because metastatic cells show organ preference in adhering to organ-derived microvascular endothelial cells. The adhesion molecules in this process remain largely unidentified. I have examined a series of murine RAW117 large-cell lymphoma cells variants selected in vivo for liver-colonizing properties ($\rm{H10{>>}L17>P}$). The highly liver-metastatic H10 cells were found to differentially express much higher levels of integrin $\alpha\rm\sb{v}\beta\sb3$ than L17 or P cells. H10 cells also adhered at higher rates to vitronectin and fibronectin than to fibrinogen, fibrin, laminin and type I collagen, and adhered at significantly higher rates to (GRGDS)$\sb4$ than to monomeric RGD-peptides. In contrast, P and L17 cells did not adhere well to the above substrates. H10 cells also spread well on vitronectin and migrated toward vitronectin concentration gradients. Pretreament of H10 cells with anti-$\beta\sb3$ monoclonal antibodies resulted in significant decreases in adhesion of H10 cells to vitronectin and immobilized (GRGDS)$\sb4$, and reduced the formation of experimental liver metastases in syngeneic Balb/c mice.^ Adhesion of RAW117 cells under hydrodynamic shear stresses was also studied because tumor cell adhesion occurs under fluid shear stresses in target organ microvessels. Similar to their properties found with static adhesion assays, H10 cells stabilized their hydrodynamic adhesion to vitronectin, fibronectin and (GRGDS)$\sb4$ much more quickly than P or L17 cells. Unlike their static adhesion properties, RAW117 cells showed differential adhesion stabilization to liver-sinusoidal endothelial cell-derived extracellular matrix ($\rm{H10{>>}L17>P}$). Although not supporting static adhesion of RAW117 cells, monomeric RGD-peptides mediated adhesion stabilization of H10 cells but not L17 or P cells. Integrin $\rm\alpha\sb{v}\beta\sb3$ was found to be involved in stabilizing H10 cell adhesion to vitronectin, (GRGDS)$\sb4$, monomeric RGD-peptide R1, and liver sinusoidal endothelial cell-derived extracellular matrix.^ This study is the first to provide evidence that integrin $\rm\alpha\sb{v}\beta\sb3$ is differentially expressed in liver-metastatic lymphoma cells and involved in differential adhesion of these cells. The results indicate that strong static adhesion and especially the unique hydrodynamic adhesion of RAW117 cells to the RGD-containing substrates correlate with liver-metastatic potentials. Thus, integrin $\rm\alpha\sb{v}\beta\sb3$ may play an important role in liver-preferential metastasis of RAW117 large-cell lymphoma cells. ^

A study in the molecular and cellular functions of GSK3beta in prostate adenocarcinoma

Kinases are part of a complex network of signaling pathways that enable a cell to respond to changes in environmental conditions in a regulated and coordinated way. For example, Glycogen Synthase Kinase 3 beta (GSK3β) modulates conformational changes, protein-protein interaction, protein degradation, and activation of unique domains in proteins that transduce signals from the extracellular milieu to the nucleus. ^ In this project, I investigated the expression and function that GSK3β exhibits in prostate cells. The capacity of GSK3β to regulate two transcription factors (JUN and CREB), which are known to be inversely utilized in prostate tumor cells, was measured. JUN/AP1 is constitutively activated in PC-3 cells; whereas, CREB/CRE activity is ∼20 fold less than the former. GSK3β overexpression obliterates JUN/AP1 activity. With respect to CREB GSK3β increases CREB/CRE activity. Cellular levels of active GSK3β can determine whether JUN or CREB is preferentially active in the PC-3s. Theoretically, in response to a particular cellular context or stimulus, a cell may coordinate JUN and CREB function by regulating GSK3β.^ A comparison of various prostate cell lines showed that active GSK3β is less expressed in normal prostate epithelial cells than in tumor cells. Differentially expressed active (GSK3β) may correlate with progression of prostate carcinoma. If a known marker associated with carcinoma of the prostate could be shown to be regulated by GSK3β then, further study of GSK3β may lead to a better understanding of both possible prevention of the disease and improved therapy for advanced stages. ^ The androgen receptor (AR) is an intriguing phosphoprotein whose regulation is potentially determined by a variety of kinases. One of these is (GSK3β) I found that (GSK3β) is a regulator of the androgen receptor in both the unliganded and liganded states. It can inhibit AR function as measured by reporter assays. Also, GSK3β associates with the AR at the DNA binding domain because deletion constructs expressing either the n-terminus or the c-terminus (both having the DBD in common) immunoprecipitated with GSK3β. Increased understanding of how GSK3β functions in prostate cancer would provide clues into how (1) certain signal pathways are coordinated and (2) the androgen receptor may be regulated. ^

Towards the identification of a tumor suppressor gene on chromosome 3p12: Physical mapping and candidate gene screening

Renal cell carcinoma (RCC) is the most common malignant tumor of the kidney. Characterization of RCC tumors indicates that the most frequent genetic event associated with the initiation of tumor formation involves a loss of heterozygosity or cytogenetic aberration on the short arm of human chromosome 3. A tumor suppressor locus Nonpapillary Renal Carcinoma-1 (NRC-1, OMIM ID 604442) has been previously mapped to a 5–7 cM region on chromosome 3p12 and shown to induce rapid tumor cell death in vivo, as demonstrated by functional complementation experiments. ^ To identify the gene that accounts for the tumor suppressor activities of NRC-1, fine-scale physical mapping was conducted with a novel real-time quantitative PCR based method developed in this study. As a result, NRC-1 was mapped within a 4.6-Mb region defined by two unique sequences within UniGene clusters Hs.41407 and Hs.371835 (78,545Kb–83,172Kb in the NCBI build 31 physical map). The involvement of a putative tumor suppressor gene Robo1/Dutt1 was excluded as a candidate for NRC-1. Furthermore, a transcript map containing eleven candidate genes was established for the 4.6-Mb region. Analyses of gene expression patterns with real-time quantitative RT-PCR assays showed that one of the eleven candidate genes in the interval (TSGc28) is down-regulated in 15 out of 20 tumor samples compared with matched normal samples. Three exons of this gene have been identified by RACE experiments, although additional exon(s) seem to exist. Further gene characterization and functional studies are required to confirm the gene as a true tumor suppressor gene. ^ To study the cellular functions of NRC-1, gene expression profiles of three tumor suppressive microcell hybrids, each containing a functional copy of NRC-1, were compared with those of the corresponding parental tumor cell lines using 16K oligonucleotide microarrays. Differentially expressed genes were identified. Analyses based on the Gene Ontology showed that introduction of NRC-1 into tumor cell lines activates genes in multiple cellular pathways, including cell cycle, signal transduction, cytokines and stress response. NRC-1 is likely to induce cell growth arrest indirectly through WEE1. ^

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

Involvement of histone H3 methylation in Tup1-Ssn6 repression

The Tup1-Ssn6 complex regulates the expression of diverse classes of genes in Saccharomyces cerevisiae including those regulated by mating type, DNA damage, glucose, and anaerobic stress. The complex is recruited to target genes by sequence-specific repressor proteins. Once recruited to particular promoters, it is not completely clear how it functions to block transcription. Repression probably occurs through interactions with both the basal transcriptional machinery and components of chromatin. Tup1 interactions with chromatin are strongly influenced by acetylation of histories H3 and H4. Tup1 binds to underacetylated histone tails and requires multiple histone deacetylases (HDACs) for its repressive functions. Like acetylation, histone methylation is involved in regulation of gene expression. The possible role of histone methylation in Tup1 repression is not known. Here we examined possible roles of histone methyltransferases in Tup1-Ssn6 functions. We found that like other genes, Tup1-Ssn6 target genes exhibit increases in the levels of histone H3 lysine 4 methylation upon activation. However, deletion of individual or multiple histone methyltransferases (HMTs) and other SET-domain containing genes has no apparent effect on Tup1-Ssn6 mediated repression of a number of well-defined targets. Interestingly, we discovered that Ssn6 interacts with Set2. Since deletion of SET2 does not affect Tup1-Ssn6 repression, Ssn6 may utilize Set2 in other contexts to regulate gene repression. In order examine if the two components of the Tup1-Ssn6 complex have independent functions in the cell, we identified genes differentially expressed in tup1Δ and ssn6Δ mutants using DNA microarrays. Our data indicate that ∼4% of genes in the cell are regulated by Ssn6 independently of Tup1. In addition, expression of genes regulated by Tup1-Ssn6 seems to be differently affected by deletion of Ssn6 and deletion of Tup1, which indicates that these components might have separate functions. Our data shed new light on the classical view of Tup1-Ssn6 functions, and indicate that Ssn6 might have repressive functions as well. ^

Metastatic lung cancer in a new mouse model inheriting p53 and latent K-ras mutations

Lung cancer is a devastating disease with very poor prognosis. The design of better treatments for patients would be greatly aided by mouse models that closely resemble the human disease. The most common type of human lung cancer is adenocarcinoma with frequent metastasis. Unfortunately, current models for this tumor are inadequate due to the absence of metastasis. Based on the molecular findings in human lung cancer and metastatic potential of osteosarcomas in mutant p53 mouse models, I hypothesized that mice with both K-ras and p53 missense mutations might develop metastatic lung adenocarcinomas. Therefore, I incorporated both K-rasLA1 and p53RI72HΔg alleles into mouse lung cells to establish a more faithful model for human lung adenocarcinoma and for translational and mechanistic studies. Mice with both mutations ( K-rasLA1/+ p53R172HΔg/+) developed advanced lung adenocarcinomas with similar histopathology to human tumors. These lung adenocarcinomas were highly aggressive and metastasized to multiple intrathoracic and extrathoracic sites in a pattern similar to that seen in lung cancer patients. This mouse model also showed gender differences in cancer related death and developed pleural mesotheliomas in 23.2% of them. In a preclinical study, the new drug Erlotinib (Tarceva) decreased the number and size of lung lesions in this model. These data demonstrate that this mouse model most closely mimics human metastatic lung adenocarcinoma and provides an invaluable system for translational studies. ^ To screen for important genes for metastasis, gene expression profiles of primary lung adenocarcinomas and metastases were analyzed. Microarray data showed that these two groups were segregated in gene expression and had 79 highly differentially expressed genes (more than 2.5 fold changes and p<0.001). Microarray data of Bub1b, Vimentin and CCAM1 were validated in tumors by quantitative real-time PCR (QPCR). Bub1b , a mitotic checkpoint gene, was overexpressed in metastases and this correlated with more chromosomal abnormalities in metastatic cells. Vimentin, a marker of epithelial-mesenchymal transition (EMT), was also highly expressed in metastases. Interestingly, Twist, a key EMT inducer, was also highly upregulated in metastases by QPCR, and this significantly correlated with the overexpression of Vimentin in the same tumors. These data suggest EMT occurs in lung adenocarcinomas and is a key mechanism for the development of metastasis in K-ras LA1/+ p53R172HΔg/+ mice. Thus, this mouse model provides a unique system to further probe the molecular basis of metastatic lung cancer.^

Proteomic identification and functional characterization of prohibitin 1 and 2 in T cell activation, expansion, survival and disease

Several immune pathologies are the result of aberrant regulation of T lymphocytes. Pronounced T cell proliferation can result in autoimmunity or hematologic malignancy, whereas loss of T cell activity can manifest as immunodeficiency. Thus, there is a critical need to characterize the signal transduction pathways that mediate T cell activation so that novel and rational strategies to detect and effectively control T cell mediated disease can be achieved. ^ The first objective of this dissertation was to identify and characterize novel T cell regulatory proteins that are differentially expressed upon antigen induced activation. Using a functional proteomics approach, two members of the prohibitin (Phb) family of proteins, Phb1 and Phb2, were determined to be upregulated upon activation of primary human T cells. Furthermore, their regulated expression was dependent upon CD3 and CD28 signaling pathways which synergistically increased their expression. In contrast to previous reports of Phb nuclear localization, both proteins were determined to localize to the mitochondrial inner membrane of human T cells. Additionally, novel Phb phosphorylation sites were identified and characterized using mass spectrometry, phosphospecific antibodies and site directed mutagenesis. ^ Prohibitins have been proposed to play important roles in cancer development however the mechanism of action has not been elucidated. The second objective of this dissertation was to define the functional role of Phbs in T cell activity, survival and disease. Compared to levels in normal human T cells, Phb expression was higher in the human tumor T cell line Kit225 and subcellularly localized to the mitochondrion. Ablation of Phb expression by siRNA treatment of Kit225 cells resulted in disruption of mitochondrial membrane potential and significantly enhanced their sensitivity to cell death, suggesting they serve a protective function in T cells. Furthermore, Q-RT-PCR analysis of human oncology cDNA expression libraries indicated the Phbs may represent hematological cancer biomarkers. Indeed, Phb1 and Phb2 protein levels were 6-10 fold higher in peripheral blood mononuclear cells isolated from malignant lymphoma and multiple myeloma patients compared to healthy individuals. ^ Taken together, Phb1 and Phb2 are novel phosphoproteins upregulated during T cell activation and transformation to function in the maintenance of mitochondrial integrity and perhaps energy metabolism, thus representing previously unrecognized intracellular biomarkers and therapeutic targets for regulating T cell activation and hematologic malignancies. ^

Glioma-specific alternative RNA splicing: A role for polypyrimidine tract binding protein-1

Alternative RNA splicing plays an integral role in cell fate determination and function, especially in the cells of the brain. Errors in RNA processing contribute to diseases such as cancer, where it leads to the production of oncogenic proteins or the loss of tumor suppressors. In silica mining suggests that hundreds of splice isoforms are misexpressed in the glial cell-derived glioma. However, there is little experimental evidence of the prevalence and contribution of these changes and whether they contribute to the formation and progression of this devastating malignancy. To determine the frequency of these aberrant events, global profiling of alternative RNA splice patterns in glioma and nontumor brain was conducted using an exon array. Most splicing changes were less than 5-fold in magnitude and 14 cassette exon events were validated, including 7 previously published events. To determine the possible causes of missplicing, the differential expression levels of splicing factors in these two tissues were also analyzed. Six RNA splicing factors had greater than 2-fold changes in expression. The highest differentially expressed factor was polypyrimidine tract binding protein-1 (PTB). Evaluation by immunohistochemistry determined that this factor was elevated in both early and late stages of glioma. Glial cell-specific PTB expression in the adult brain led me to examine the role of PTB in gliomagenesis. Downregulation of PTB slowed glioma cell proliferation and migration and enhanced cell adhesion to fibronectin and vitronectin. To determine whether PTB was affecting these processes through splicing, genome-wide exon expression levels were correlated with PTB levels. Surprisingly, previously reported PTB target transcripts were insensitive to changes in PTB levels in both patient samples and PTB-depleted glioma cells. Only one validated glioma-specific splice target, RTN4/Nogo, had a significant PTB-mediated splicing change. Downregulation of PTB enhanced inclusion of its alternative exon 3, which encodes an auxiliary domain within a neurite inhibitor protein. Overexpression of this splice isoform in glioma cells slowed proliferation in a manner similar to that observed in PTB knockdown cells. In summary, aberrant expression of splicing factors such as PTB in glioma may elicit changes in splicing patterns that enhance tumorigenesis. ^