Identification of a novel tumor suppressor gene on human chromosome 3p14-p12 involved in renal cell carcinoma

Nonpapillary renal cell carcinoma (RCC) is an adult cancer of the kidney which occurs both in familial and sporadic forms. The familial form of RCC is associated with translocations involving chromosome 3 with a breakpoint at 3p14-p13. Studies focused on sporadic RCC have shown two commonly deleted regions at 3p14.3-p13 and 3p21.3. In addition, a more distal region mapping to 3p26-p25 has been linked to the Von Hippel Lindau (VHL) disease gene. A large proportion of VHL patients develop RCC. The short arm of human chromosome 3 can, therefore, be dissected into three distinct regions which could encode tumor suppressor genes for RCC. Loss or inactivation of one or more of these loci may be an important step in the genesis of RCC.^ I have used the technique of microcell-mediated chromosome transfer to introduce an intact, normal human chromosome 3 and defined fragments of 3p, dominantly marked with pSV2neo, into the highly malignant RCC cell line SN12C.19. The introduction of chromosome 3 and of a centric fragment of 3p, encompassing 3p14-q11, into SN12C.19 resulted in dramatic suppression of tumor growth in nude mice. Another defined deletion hybrid contained the region 3p12-q24 of the introduced human chromosome and failed to suppress tumorigenicity. These data define the region 3p14-p12, the most proximal region of high frequency allele loss in sporadic RCC as well as the region containing the translocation breakpoint in familial RCC, to contain a novel tumor suppressor locus involved in RCC. We have designated this locus nonpapillary renal cell carcinoma-1 (NRC-1). Furthermore, we have functional evidence that NRC-1 controls the growth of RCC cells by inducing rapid cell death in vivo. ^

Network-based genome-wide association study of Crohn's disease

Genome-wide association studies (GWAS) have rapidly become a standard method for disease gene discovery. Many recent GWAS indicate that for most disorders, only a few common variants are implicated and the associated SNPs explain only a small fraction of the genetic risk. The current study incorporated gene network information into gene-based analysis of GWAS data for Crohn's disease (CD). The purpose was to develop statistical models to boost the power of identifying disease-associated genes and gene subnetworks by maximizing the use of existing biological knowledge from multiple sources. The results revealed that Markov random field (MRF) based mixture model incorporating direct neighborhood information from a single gene network is not efficient in identifying CD-related genes based on the GWAS data. The incorporation of solely direct neighborhood information might lead to the low efficiency of these models. Alternative MRF models looking beyond direct neighboring information are necessary to be developed in the future for the purpose of this study.^

Theoretical and experimental studies of linkage disequilibrium in human populations

Linkage disequilibrium (LD) is defined as the nonrandom association of alleles at two or more loci in a population and may be a useful tool in a diverse array of applications including disease gene mapping, elucidating the demographic history of populations, and testing hypotheses of human evolution. However, the successful application of LD-based approaches to pertinent genetic questions is hampered by a lack of understanding about the forces that mediate the genome-wide distribution of LD within and between human populations. Delineating the genomic patterns of LD is a complex task that will require interdisciplinary research that transcends traditional scientific boundaries. The research presented in this dissertation is predicated upon the need for interdisciplinary studies and both theoretical and experimental projects were pursued. In the theoretical studies, I have investigated the effect of genotyping errors and SNP identification strategies on estimates of LD. The primary importance of these two chapters is that they provide important insights and guidance for the design of future empirical LD studies. Furthermore, I analyzed the allele frequency distribution of 26,530 single nucleotide polymorphisms (SNPs) in three populations and generated the first-generation natural selection map of the human genome, which will be an important resource for explaining and understanding genomic patterns of LD. Finally, in the experimental study, I describe a novel and simple, low-cost, and high-throughput SNP genotyping method. The theoretical analyses and experimental tools developed in this dissertation will facilitate a more complete understanding of patterns of LD in human populations. ^

A full pedigree based method for the statistical assessment of genetic anticipation

Genetic anticipation is defined as a decrease in age of onset or increase in severity as the disorder is transmitted through subsequent generations. Anticipation has been noted in the literature for over a century. Recently, anticipation in several diseases including Huntington's Disease, Myotonic Dystrophy and Fragile X Syndrome were shown to be caused by expansion of triplet repeats. Anticipation effects have also been observed in numerous mental disorders (e.g. Schizophrenia, Bipolar Disorder), cancers (Li-Fraumeni Syndrome, Leukemia) and other complex diseases. ^ Several statistical methods have been applied to determine whether anticipation is a true phenomenon in a particular disorder, including standard statistical tests and newly developed affected parent/affected child pair methods. These methods have been shown to be inappropriate for assessing anticipation for a variety of reasons, including familial correlation and low power. Therefore, we have developed family-based likelihood modeling approaches to model the underlying transmission of the disease gene and penetrance function and hence detect anticipation. These methods can be applied in extended families, thus improving the power to detect anticipation compared with existing methods based only upon parents and children. The first method we have proposed is based on the regressive logistic hazard model. This approach models anticipation by a generational covariate. The second method allows alleles to mutate as they are transmitted from parents to offspring and is appropriate for modeling the known triplet repeat diseases in which the disease alleles can become more deleterious as they are transmitted across generations. ^ To evaluate the new methods, we performed extensive simulation studies for data simulated under different conditions to evaluate the effectiveness of the algorithms to detect genetic anticipation. Results from analysis by the first method yielded empirical power greater than 87% based on the 5% type I error critical value identified in each simulation depending on the method of data generation and current age criteria. Analysis by the second method was not possible due to the current formulation of the software. The application of this method to Huntington's Disease and Li-Fraumeni Syndrome data sets revealed evidence for a generation effect in both cases. ^

Investigations of stress responses in Saccharomyces cerevisiae: Transcriptional control and heat shock protein function

The baker's yeast, Saccharomyces cerevisiae responds to the cytotoxic effects of elevated temperature (37-42°C) by activating transcription of ∼150 genes, termed heat shock genes, collectively required to compensate for the abundance of misfolded and aggregated proteins and various physiological modifications necessary for the cell to survive and grow at heat shock temperatures. An intriguing facet of the yeast heat shock response is the remarkable similarity it shares with the global remodeling that occurs in mammalian cells in response to numerous pathophysiological conditions including cancer and cardiovascular disease and thus provides an ideal model system. I have therefore investigated several novel features of stress signaling, transcriptional regulation, and physiology. Initial work focused on the characterization of SYM1, a novel heat shock gene in yeast which was demonstrated to be required for growth on the nonfermentable carbon source ethanol at elevated temperature, and to be the functional ortholog of the mammalian kidney disease gene, Mpv17. Additional work addressed the role of two proteins, the Akt-related kinase, Sch9, and Sse1, the yeast Hsp110 protein chaperone homolog, in signaling by protein kinase A, establishing Sse1 as a critical negative regulator of this pathway. Furthermore, I have demonstrated a role for Sse1 in biogenesis and stability of the stress-response transcription factor, Msn2; a finding that has been extended to include a select subset of additional high molecular weight proteins, suggesting a more global role for this chaperone in stabilizing the cellular proteome. The final emphasis of my doctoral work has included the finding that celastrol, a compound isolated from the plant family Celasfraceae, a component of traditional Chinese herbal medicine, can activate heat shock transcription factor (Hsf1) in yeast and mammalian cells through an oxidative stress mechanism. Celastrol treatment simultaneously activates both heat shock and oxidative stress response pathways, resulting in increased cytoprotection. ^

Identification of germ cell tumor modifier genes from the 129.MOLF-Chr19 consomic mouse strain

Naturally occurring genetic variants confer susceptibility to disease in the human population, including in testicular germ cell tumor development. Disease susceptibility loci for testicular germ cell tumors have been identified by genetic mapping in humans and mice. However, the identity of many of the susceptibility genes remains unclear. My study utilized a chromosome substitution strain, the 129.MOLF-Chr 19 (or M19 strain), to identify candidate testicular germ cell tumor susceptibility genes. Males of this strain have a high incidence of germ cell tumors in the testes. By forward genetic approaches, five susceptibility loci were fine-mapped and the genetic interactions were dissected. In addition, I identified three protein-coding genes and one micro-RNA as testicular tumor susceptibility genes by genomic screening. Using reverse genetic approaches, I verified one of the candidates, Splicing factor 1, as a modifier of testicular tumor. Deficiency of SF1 significantly reduces the incidence of testicular tumors in mice. This study highlights the advantage of the 129.MOLF-Chr 19 consomic strain in disease gene identification and validation. It also sets the stage to elucidate the molecular mechanisms of tumorigenesis in the testis. ^

From gene to protein: Investigating the disease etiology of retinitis pigmentosa

Retinitis pigmentosa (RP) is a name given to a group of inherited retinal dystrophies that lead to progressive photoreceptor degeneration, and thus, visual impairment. It is evident at both the clinical and the molecular level that these are heterogeneous disorders, with wide variation in severity, mode of inheritance, and phenotype. The genetics of RP are not simple; the disease can be inherited in dominant, recessive, X-linked, and digenic modes. Autosomal dominant RP (adRP) results from mutations in at least ten mapped loci, but there may be dozens of genetic loci where mutations can cause RP. To date, there are over a hundred genes known to cause retinal degenerative diseases, and less than half of these have been cloned (RetNet). Among the dozens of retinitis pigmentosa loci known to exist, only a few have been identified and the remainders are inferred from linkage studies. Today, the genes for seven of the twelve-adRP loci have been identified, and these are rhodopsin, peripherin/RDS, NRL, ROM1, CRX, RP13 and RP1. My research projects involved a combination of the continued search for genes involved in retinal dystrophies, as well the investigation into the role of peripherin/RDS and RP1 in the disease etiology of autosomal dominant RP. ^ Most of the mutations leading to inherited retinal disorders have been identified in predominately retina expressed genes like rhodopsin, peripherin/RDS, and RP1. Expressed sequence tags (ESTs) that were retina-specific were culled from sequence databases and, together with laboratory analysis, were analyzed as potential candidate genes for retinal dystrophies. Thirteen of the fifty-five identified retina-specific ESTs mapped to within candidate regions for inherited retinopathies. One of these is RP1L1, a homologue of RP1 and a potential cause of adRP. ^ Once a disease-associated gene has been identified, elucidating the role of that gene in the visual process is essential for understanding what happens when the process is defective as it is in adRP. My next projects involved investigating the role of a novel 5′ donor +3 splice site mutation on the mRNA of peripherin/RDS in adRP affected individuals, and comparative sequencing in RP1 to define conserved regions of the protein. Comparative sequencing is a powerful way to delineate critical regions of a sequence because different regions of a gene have different functions, and each region is subject to different levels of functional or structural constraints. Establishing a framework of conserved domains is beneficial not only for structural or functional studies, but can also aid in determining the potential effects of mutations. With the completion of sequencing of human genome, and other organisms such as Saccharomyces cerevisiae, Caenorhabditis elegans , and Drosophila, the facility of comparative sequencing will only increase in the future. Comparative sequencing has already become an established procedure for pinpointing conserved regions of a protein, and is an efficient way to target regions of a protein for experimental and/or evolutionary analysis. ^

Adiposity and the perilipin gene: The role of single locus and multilocus variation, obesity-related quantitative traits, and disease-related phenotypes in the Atherosclerosis Risk in Communities (ARIC) study

Obesity is a complex multifactorial disease and is a public health priority. Perilipin coats the surface of lipid droplets in adipocytes and is believed to stabilize these lipid bodies by protecting triglyceride from early lipolysis. This research project evaluated the association between genetic variation within the human perilipin (PLIN) gene and obesity-related quantitative traits and disease-related phenotypes in Non-Hispanic White (NHW) and African American (AA) participants from the Atherosclerosis Risk in Communities (ARIC) Study. ^ Multivariate linear regression, multivariate logistic regression, and Cox proportional hazards models evaluated the association between single gene variants (rs2304794, rs894160, rs8179071, and rs2304795) and multilocus variation (rs894160 and rs2304795) within the PLIN gene and both obesity-related quantitative traits (body weight, body mass index [BMI], waist girth, waist-to-hip ratio [WHR], estimated percent body fat, and plasma total triglycerides) and disease-related phenotypes (prevalent obesity, metabolic syndrome [MetS], prevalent coronary heart disease [CHD], and incident CHD). Single variant analyses were stratified by race and gender within race while multilocus analyses were stratified by race. ^ Single variant analyses revealed that rs2304794 and rs894160 were significantly related to plasma triglyceride levels in all NHWs and NHW women. Among AA women, variant rs8179071 was associated with triglyceride levels and rs2304794 was associated with risk-raising waist circumference (>0.8 in women). The multilocus effects of variants rs894160 and rs2304795 were significantly associated with body weight, waist girth, WHR, estimated percent body fat, class II obesity (BMI ≥ 35 kg/m2), class III obesity (BMI ≥ 35 kg/m2), and risk-raising WHR (>0.9 in men and >0.8 in women) in AAs. Variant rs2304795 was significantly related to prevalent MetS among AA males and prevalent CHD in NHW women; multilocus effects of the PLIN gene were associated with prevalent CHD among NHWs. Rs2304794 was associated with incident CHD in the absence of the MetS among AAs. These findings support the hypothesis that variation within the PLIN gene influences obesity-related traits and disease-related phenotypes. ^ Understanding these effects of the PLIN genotype on the development of obesity can potentially lead to tailored health promotion interventions that are more effective. ^

Genome-wide gene-gene interaction analysis for cardiovascular disease

Numerous studies have been carried out to try to better understand the genetic predisposition for cardiovascular disease. Although it is widely believed that multifactorial diseases such as cardiovascular disease is the result from effects of many genes which working alone or interact with other genes, most genetic studies have been focused on identifying of cardiovascular disease susceptibility genes and usually ignore the effects of gene-gene interactions in the analysis. The current study applies a novel linkage disequilibrium based statistic for testing interactions between two linked loci using data from a genome-wide study of cardiovascular disease. A total of 53,394 single nucleotide polymorphisms (SNPs) are tested for pair-wise interactions, and 8,644 interactions are found to be significant with p-values less than 3.5×10-11. Results indicate that known cardiovascular disease susceptibility genes tend not to have many significantly interactions. One SNP in the CACNG1 (calcium channel, voltage-dependent, gamma subunit 1) gene and one SNP in the IL3RA (interleukin 3 receptor, alpha) gene are found to have the most significant pair-wise interactions. Findings from the current study should be replicated in other independent cohort to eliminate potential false positive results.^

Genetic epidemiology of gallbladder disease in Mexican Americans and cholesterol 7$\alpha$-hydroxylase gene sequence variation

Among Mexican Americans, the second largest minority group in the United States, the prevalence of gallbladder disease is markedly elevated. Previous data from both genetic admixture and family studies indicate that there is a genetic component to the occurrence of gallbladder disease in Mexican Americans. However, prior to this thesis no formal genetic analysis of gallbladder disease had been carried out nor had any contributing genes been identified.^ The results of complex segregation analysis in a sample of 232 Mexican American pedigrees documented the existence of a major gene having two alleles with age- and gender-specific effects influencing the occurrence of gallbladder disease. The estimated frequency of the allele increasing susceptibility was 0.39. The lifetime probabilities that an individual will be affected by gallbladder disease were 1.0, 0.54, and 0.00 for females of genotypes "AA", "Aa", and "aa", respectively, and 0.68, 0.30, and 0.00 for males, respectively. This analysis provided the first conclusive evidence for the existence of a common single gene having a large effect on the occurrence of gallbladder disease.^ Human cholesterol 7$\alpha$-hydroxylase is the rate-limiting enzyme in bile acid synthesis. The results of an association study in both a random sample and a matched case/control sample showed that there is a significant association between cholesterol 7$\alpha$-hydroxylase gene variation and the occurrence of gallbladder disease in Mexican Americans males but not in females. These data have implicated a specific gene, 7$\alpha$-hydroxylase, in the etiology of gallbladder disease in this population.^ Finally, I asked whether the inferred major gene from complex segregation analysis is genetically linked to the cholesterol 7$\alpha$-hydroxylase gene. Three pedigrees predicted to be informative for linkage analysis by virtue of supporting the major gene hypothesis and having parents with informative genotypes and multiple offspring were selected for this linkage analysis. In each of these pedigrees, the recombination fractions maximized at 0 with a positive, albeit low, LOD score. The results of this linkage analysis provide preliminary and suggestive evidence that the cholesterol 7$\alpha$-hydroxylase gene and the inferred gallbladder disease susceptibility gene are genetically linked. ^

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.

Mutual Information for Testing Gene-Environment Interaction

Despite current enthusiasm for investigation of gene-gene interactions and gene-environment interactions, the essential issue of how to define and detect gene-environment interactions remains unresolved. In this report, we define gene-environment interactions as a stochastic dependence in the context of the effects of the genetic and environmental risk factors on the cause of phenotypic variation among individuals. We use mutual information that is widely used in communication and complex system analysis to measure gene-environment interactions. We investigate how gene-environment interactions generate the large difference in the information measure of gene-environment interactions between the general population and a diseased population, which motives us to develop mutual information-based statistics for testing gene-environment interactions. We validated the null distribution and calculated the type 1 error rates for the mutual information-based statistics to test gene-environment interactions using extensive simulation studies. We found that the new test statistics were more powerful than the traditional logistic regression under several disease models. Finally, in order to further evaluate the performance of our new method, we applied the mutual information-based statistics to three real examples. Our results showed that P-values for the mutual information-based statistics were much smaller than that obtained by other approaches including logistic regression models.

Regulation of survivin gene expression in the human endometrium and endometrial cancer

In the United States, endometrial cancer is the leading cancer of the female reproductive tract. There are 40,100 new cases and 7,470 deaths from endometrial cancer estimated for 2008 (47). The average five year survival rate for endometrial cancer is 84% however, this figure is substantially lower in patients diagnosed with late stage, advanced disease and much higher for patients diagnosed in early stage disease (47). Endometrial cancer (EC) has been associated with several risk factors including obesity, diabetes, hypertension, previously documented occurrence of hereditary non-polyposis colorectal cancer (HNPCC), and heightened exposure to estrogen (25). As of yet, there has not been a dependable molecular predictor of endometrial cancer occurrence in women with these predisposing factors. The goal of our lab is to identify genes that are aberrantly expressed in EC and may serve as molecular biomarkers of EC progression. One candidate protein that we are exploring as a biomarker of EC progression is the cell survival protein survivin.

Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

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.