Genotypic spectrum and genotype-phenotype correlation of trimethylaminuria

Trimethylaminuria (TMAU) or Fish odor syndrome is an autosomal recessive disease that is characterized by pungent body odor with subsequent psychosocial complications. There are limited studies of the sequence variants causing TMAU in the literature with most studies describing only one or two patients and lacking genotype-phenotype correlations. Also to date, there is no laboratory in the US or Europe that offers TMA genetic testing on a clinical basis. We have recently validated genetic testing in the University of Colorado DNA Diagnostic Laboratory. We have a database of a few dozen patients with a biochemical diagnosis of TMA at the University of Colorado at Denver Health Sciences Center (UCDHSC) which includes a few patients with the classical form of the disease. We have used the newly established clinical test in our institution to attempt to characterize the genotype (sequence variants including mutations and polymorphisms) of classical TMAU patients and to establish a genotype-phenotype (biochemical and clinical) association. The questionnaire results confirmed most of the previously reported epidemiological findings of TMAU and also indicated that TMAU patients use multiple intervention measures in attempt to control their symptoms with dietary control being most effective. Despite the complexity of intervention, most patients did not have any medical follow up and there was underutilization of specialist care. In a set of our patients, two deleterious mutations were identified in 4/12 patients including a novel T237P sequence variant, while the majority of our patients (8/12) did not reveal any mutations. Some of the latter were double heterozygous for the E158K and E308G polymorphisms which could explain a mild phenotype while others had only the E158K variant which raised the question of undetected mutations. These results indicate that further experiments are needed to further delineate the full mutational spectrum of the FMO3 gene. ^

Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

Extending in vitro conditioning in Aplysia to analyze operant and classical processes in the same preparation.

Operant and classical conditioning are major processes shaping behavioral responses in all animals. Although the understanding of the mechanisms of classical conditioning has expanded significantly, the understanding of the mechanisms of operant conditioning is more limited. Recent developments in Aplysia are helping to narrow the gap in the level of understanding between operant and classical conditioning, and have raised the possibility of studying the neuronal processes underlying the interaction of operant and classical components in a relatively complex learning task. In the present study, we describe a first step toward realizing this goal, by developing a single in vitro preparation in which both operant and classical conditioning can be studied concurrently. The new paradigm reproduced previously published results, even under more conservative and homogenous selection criteria and tonic stimulation regime. Moreover, the observed learning was resistant to delay, shortening, and signaling of reinforcement.

In vitro analog of classical conditioning of feeding behavior in aplysia.

The feeding behavior of Aplysia californica can be classically conditioned using tactile stimulation of the lips as a conditioned stimulus (CS) and food as an unconditioned stimulus (US). Moreover, several neural correlates of classical conditioning have been identified. The present study extended previous work by developing an in vitro analog of classical conditioning and by investigating pairing-specific changes in neuronal and synaptic properties. The preparation consisted of the isolated cerebral and buccal ganglia. Electrical stimulation of a lip nerve (AT4) and a branch of the esophageal nerve (En2) served as the CS and US, respectively. Three protocols were used: paired, unpaired, and US alone. Only the paired protocol produced a significant increase in CS-evoked fictive feeding. At the cellular level, classical conditioning enhanced the magnitude of the CS-evoked synaptic input to pattern-initiating neuron B31/32. In addition, paired training enhanced both the magnitude of the CS-evoked synaptic input and the CS-evoked spike activity in command-like neuron CBI-2. The in vitro analog of classical conditioning reproduced all of the cellular changes that previously were identified following behavioral conditioning and has led to the identification of several new learning-related neural changes. In addition, the pairing-specific enhancement of the CS response in CBI-2 indicates that some aspects of associative plasticity may occur at the level of the cerebral sensory neurons.

Predictive oncology: a review of multidisciplinary, multiscale in silico modeling linking phenotype, morphology and growth.

Empirical evidence and theoretical studies suggest that the phenotype, i.e., cellular- and molecular-scale dynamics, including proliferation rate and adhesiveness due to microenvironmental factors and gene expression that govern tumor growth and invasiveness, also determine gross tumor-scale morphology. It has been difficult to quantify the relative effect of these links on disease progression and prognosis using conventional clinical and experimental methods and observables. As a result, successful individualized treatment of highly malignant and invasive cancers, such as glioblastoma, via surgical resection and chemotherapy cannot be offered and outcomes are generally poor. What is needed is a deterministic, quantifiable method to enable understanding of the connections between phenotype and tumor morphology. Here, we critically assess advantages and disadvantages of recent computational modeling efforts (e.g., continuum, discrete, and cellular automata models) that have pursued this understanding. Based on this assessment, we review a multiscale, i.e., from the molecular to the gross tumor scale, mathematical and computational "first-principle" approach based on mass conservation and other physical laws, such as employed in reaction-diffusion systems. Model variables describe known characteristics of tumor behavior, and parameters and functional relationships across scales are informed from in vitro, in vivo and ex vivo biology. We review the feasibility of this methodology that, once coupled to tumor imaging and tumor biopsy or cell culture data, should enable prediction of tumor growth and therapy outcome through quantification of the relation between the underlying dynamics and morphological characteristics. In particular, morphologic stability analysis of this mathematical model reveals that tumor cell patterning at the tumor-host interface is regulated by cell proliferation, adhesion and other phenotypic characteristics: histopathology information of tumor boundary can be inputted to the mathematical model and used as a phenotype-diagnostic tool to predict collective and individual tumor cell invasion of surrounding tissue. This approach further provides a means to deterministically test effects of novel and hypothetical therapy strategies on tumor behavior.

Conversion of myoblasts to physiologically active neuronal phenotype.

Repressor element 1 (RE1)-silencing transcription factor (REST)/neuron-restrictive silencer factor (NRSF) can repress several terminal neuronal differentiation genes by binding to a specific DNA sequence (RE1/neuron-restrictive silencer element [NRSE]) present in their regulatory regions. REST-VP16 binds to the same RE1/NRSE, but activates these REST/NRSF target genes. However, it is unclear whether REST-VP16 expression is sufficient to cause formation of functional neurons either from neural stem cells or from heterologous stem cells. Here we show that the expression of REST-VP16 in myoblasts grown under muscle differentiation conditions blocked entry into the muscle differentiation pathway, countered endogenous REST/NRSF-dependent repression, activated the REST/NRSF target genes, and, surprisingly, activated other neuronal differentiation genes and converted the myoblasts to a physiologically active neuronal phenotype. Furthermore, in vitro differentiated neurons produced by REST-VP16-expressing myoblasts, when injected into mouse brain, survived, incorporated into the normal brain, and did not form tumors. This is the first instance in which myoblasts were converted to a neuronal phenotype. Our results suggest that direct activation of REST/NRSF target genes with a single transgene, REST-VP16, is sufficient to activate other terminal neuronal differentiation genes and to override the muscle differentiation pathways, and they suggest that this approach provides an efficient way of triggering neuronal differentiation in myoblasts and possibly other stem cells.

FoxM1B regulates NEDD4-1 expression, leading to cellular transformation and full malignant phenotype in immortalized human astrocytes.

Our recent studies have shown that the FoxM1B transcription factor is overexpressed in human glioma tissues and that the level of its expression correlates directly with glioma grade. However, whether FoxM1B plays a role in the early development of glioma (i.e., in transformation) is unknown. In this study, we found that the FoxM1B molecule causes cellular transformation and tumor formation in normal human astrocytes (NHA) immortalized by p53 and pRB inhibition. Moreover, brain tumors that arose from intracranial injection of FoxM1B-expressing immortalized NHAs displayed glioblastoma multiforme (GBM) phenotypes, suggesting that FoxM1B overexpression in immortalized NHAs not only transforms the cells but also leads to GBM formation. Mechanistically, our results showed that overexpression of FoxM1B upregulated NEDD4-1, an E3 ligase that mediates the degradation and downregulation of phosphatase and tensin homologue (PTEN) in multiple cell lines. Decreased PTEN in turn resulted in the hyperactivation of Akt, which led to phosphorylation and cytoplasmic retention of FoxO3a. Blocking Akt activation with phosphoinositide 3-kinase/Akt inhibitors inhibited the FoxM1B-induced transformation of immortalized NHAs. Furthermore, overexpression of FoxM1B in immortalized NHAs increased the expression of survivin, cyclin D1, and cyclin E, which are important molecules for tumor growth. Collectively, these results indicate that overexpression of FoxM1B, in cooperation with p53 and pRB inhibition in NHA cells, promotes astrocyte transformation and GBM formation through multiple mechanisms.

Contribution of Ectodomain Mutations in Epidermal Growth Factor Receptor to Signaling in Glioblastoma Multiforme

CONTRIBUTION OF ECTODOMAIN MUTATIONS IN EPIDERMAL GROWTH FACTOR RECEPTOR TO SIGNALING IN GLIOBLASTOMA MULTIFORME Publication No._________ Marta Rojas, M.S. Supervisory Professor: Oliver Bögler, Ph.D. The Cancer Genome Atlas (TCGA) has conducted a comprehensive analysis of a large tumor cohort and has cataloged genetic alterations involving primary sequence variations and copy number aberrations of genes involved in key signaling pathways in glioblastoma (GBM). This dataset revealed missense ectodomain point mutations in epidermal growth factor receptor (EGFR), but the biological and clinical significance of these mutations is not well defined in the context of gliomas. In our study, we focused on understanding and defining the molecular mechanisms underlying the functions of EGFR ectodomain mutants. Using proteomic approaches to broadly analyze cell signaling, including antibody array and mass spectrometry-based methods, we found a differential spectrum of tyrosine phosphorylation across the EGFR ectodomain mutations that enabled us to stratify them into three main groups that correlate with either wild type EGFR (EGFR) or the long-studied mutant, EGFRvIII. Interestingly, one mutant shared characteristics of both groups suggesting a continuum of behaviors along which different mutants fall. Surprisingly, no substantial differences were seen in activation of classical downstream signaling pathways such as Akt and S6 pathways between these classes of mutants. Importantly, we demonstrated that ectodomain mutations lead to differential tumor growth capabilities in both in vitro (anchorage independent colony formation) and in vivo conditions (xenografts). Our data from the biological characterization allowed us to categorize the mutants into three main groups: the first group typified by EGFRvIII are mutations with a more aggressive phenotype including R108K and A289T; a second group characterized by a less aggressive phenotype exemplified by EGFR and the T263P mutation; and a third group which shared characteristics from both groups and is exemplified by the mutation A289D. In addition, we treated cells overexpressing the mutants with various agents employed in the clinic including temozolomide, cisplatin and tarceva. We found that cells overexpressing the mutants in general displayed resistance to the treatments. Our findings yield insights that help with the molecular characterization of these mutants. In addition, our results from the drug studies might be valuable in explaining differential responses to specific treatments in GBM patients.

Cerebellar cortex involvement in classical eyelid conditioning

A model for cerebellar involvement in motor learning was tested using classical eyelid conditioning in the rabbit. Briefly, we assume that modifications of the strength of granule cell synapses at Purkinje cells in the cerebellar cortex and mossy fiber (MF) synapses at cerebellar interpositus nuclei are responsible for the acquisition, adaptively-timed expression, and extinction of conditioned eyelid responses (CRs). A corollary of these assumptions is that the cerebellar cortex is necessary for acquisition and extinction. This model also suggests a mechanism whereby the cerebellar cortex can discriminate different times during a conditioned stimulus (CS) and thus mediate the learned timing of CRs. Therefore, experiments were done to determine the role of the cerebellar cortex in the timing, extinction, and acquisition of CRs. Lesions of the cerebellar cortex that included the anterior lobe disrupted the learned timing of CRs such that they occurred at extremely short latencies. Stimulation of MFs in the middle cerebellar peduncle as the CS could support differently timed CRs in the same animal. These data indicate that synaptic plasticity in the cerebellar cortex mediates the learned timing of CRs. These short-latency CRs which resulted from anterior lobe damage did not extinguish, while CRs in animals receiving lesions which did not include the anterior lobe extinguished normally. Preliminary data suggests that lesions of the anterior lobe which produce short-latency responses prevent the acquisition of CRs to a novel CS. These findings indicate that the anterior lobe of cerebellar cortex is necessary for eyelid conditioning. The involvement of the anterior lobe in eyelid conditioning has not been previously reported, however, the anterior lobe has generally been spared in lesion studies examining cerebellar cortex involvement in eyelid conditioning due to its relatively inaccessible location. The observation that the anterior lobe of the cerebellar cortex is not always required for the basic expression of CRs, but is necessary for response timing, extinction, and acquisition, is consistent with the hypothesis that eyelid conditioning can involve plasticity in both the cerebellar cortex and interpositus nucleus and that plasticity in the nucleus is controlled by Purkinje cell activity. ^

Involvement of human chromosome 17 in the control of the metastatic phenotype in melanoma

Fusion of nonmetastatic murine melanoma K1735 C19H cells with metastatic human melanoma A375 C15N cells resulted in a hybrid (termed H7) which was highly metastatic in a nude mouse model. The H7 hybrid retained chromosome 17 as the sole intact human chromosome in the cell. A lung bioassay showed that the K1735 C19H cells were present in the lungs of nude mice with s.c. tumors, yet at 6-weeks after tumor resection, no cells remained in the lung and therefore did not form lung metastases. Examination of various phenotypic properties such as in vivo and in vitro growth demonstrated that phenotypically the H7 hybrid was most like the K1735 C19H cell line except for its metastatic ability. In contrast the H7 hybrid cells containing single or multiple copies of human chromosome 17 with a point mutation at codon 249 (arg-gly) of the p53 gene, readily formed lung metastases. A plasmid containing the human p53 from the H7 hybrid and four other contructs with mutations at codon 143 (val-arg), 175 (arg-his), 249 (arg-ser) and 273 (arg-his) were transfected into K1735 C19H cells. K1735 C19H cells expressing human p53 genes with mutations at codons 249, both the arg-ser mutation and the mutation from the H7 hybrid and 273 produced significantly more lung metastases.^ In vitro assays demonstrated that responses to various cytotoxic and DNA damaging agents varied with the presence of mutant p53 and with the type of agent used. When cultured in mouse lung-conditioned medium, the K1735 C19H cell line was growth-inhibited, while cells containing a mutant human p53 (either on the whole chromosome 17, as in the H7 hybrid cells or from a stably transfected construct) were growth stimulated. Western blot analysis of lung-conditioned media derived from either 6-month or 15-month old mice has detected high levels of soluble Fas ligand in the medium from older animals. Comparison of the levels of Fas receptor on the K1735 C19H cell line and the H7 hybrid were almost identical, but 50% of the K1735 C19H cells were killed in the presence of anti-Fas antibody as opposed to 7% of the H7 hybrid cells. The growth-inhibitory effects of the lung-conditioned medium on the K1735 C19H cells were abrogated by coculture with Fas-Fc, which competes with the Fas ligand for receptor binding. Growth-inhibition of the K1735 C19H was 54% when cultured in 60 $\mu$g/0.2 ml lung-conditioned medium and a control Fc, with only 9% inhibition in 60 $\mu$g/0.2 ml lung-conditioned medium and Fas-Fc. Growth of the H7 cells and K1735 C19H cells transfected with various mutant human p53 genes were unchanged by the presence of either the control Fc or the Fas-Fc. This indicates that the presence of human chromosome 17, and mutant p53 in part protects the cells from Fas:Fas ligand induced apoptosis, and allows the growth of lung metastases. ^

Effects of butyrate on the colon cancer cell phenotype-chemosensitization and upregulation of proteins implicated in cancer progression

Colon cancer is the second leading cause of cancer mortality in the U.S. Surgery is the only truly effective human colon cancer (HCC) therapy due to marked intrinsic drug resistance. The inefficacy of therapies developed for metastatic HCC suggests that advances in colon cancer chemoprevention and chemotherapy will be needed to reduce HCC mortality. The dietary fiber metabolite butyrate (NaB) is a candidate cancer chemopreventive agent that inhibits growth, promotes differentiation and stimulates apoptosis of HCC cells. Epidemiological and experimental studies suggest that dietary fiber protects against the development of HCC, however, recent large prospective trials have not found significant protection. ^ The first central hypothesis of this dissertation project is that the diversity of phenotypic changes induced by NaB in HCC cells includes molecular alterations that oppose its chemopreventive action and thereby limit its efficacy. We investigated the effect of NaB on the expression/activity of epidermal growth factor receptor (EGFR) and cyclooxygenase-2 (COX-2) in HCC HT29 cells. NaB treatment induced a 13-fold increase in EGFR expression in concert with its chemopreventive action in vitro, i.e., induction of growth suppression and G1 arrest, apoptosis and a differentiated phenotype. NaB-induced EGFR was active based on multiple lines of evidence. The EGFR was: (1) heavily phosphorylated at Tyrosine (P-Tyr); (2) associated with the cytoskeleton; (3) localized at the cell surface, and activated in response to EGF; and (4) NaB treatment of the cells induced activation of the EGFR effector Erk1/2. NaB treatment also induced a 7-fold increase in COX-2 expression. The NaB-induced COX-2 was active based on significantly increased PGE2 production. ^ The second central hypothesis is that NaB treatment would render HCC cells more chemosensitive to chemotherapy agents based on the increased apoptotic index induced by NaB. NaB treatment chemosensitized HT29 cells to 5-FU and doxorubicin, despite increases in the expression of P-glycoprotein and a related drug resistance protein (MRP). ^ These results raise the intriguing possibility that the chemopreventive effects of fiber may require concomitant treatment with EGFR and/or COX-2 inhibitors. Similarly, NaB may be a rational drug to combine with existing chemotherapeutic agents for the management of advanced HCC patients. ^

TUMOR PROGRESSION: ANALYSIS OF THE INSTABILITY OF THE METASTATIC PHENOTYPE, SENSITIVITY TO RADIATION AND CHEMOTHERAPY (PHENOTYPIC DRIFT, BREAST CANCER)

The major complications for tumor therapy are (i) tumor spread (metastasis); (ii) the mixed nature of tumors (heterogeneity); and (iii) the capacity of tumors to evolve (progress). To study these tumor characteristics, the rat 13762NF mammary adenocarcinoma was cloned and studied for metastatic properties and sensitivities to therapy (chemotherapy, radiation and hyperthermia). The cell clones were heterogeneous and no correlation between metastatic potential and therapeutic sensitivities was observed. Further, these phenotypes were unstable during passage in vitro; yet, the changes were clone dependent and reproducible using different cryoprotected cell stocks. To understand the phenotypic instability, subclones were isolated from low and high passage cell clones. Each subclone possessed a unique composite phenotype. Again, no apparent correlation was seen between metastatic potential and sensitivity to therapy. The results demonstrated that (1) tumor cells are heterogeneous for multiple phenotypes; (2) tumor cells are unstable for multiple phenotypes; (3) the magnitude, direction and time of occurrence of phenotypic drift is clone dependent; (4) the sensitivity of cell clones to ionizing radiation (gamma or heat) and chemotherapy agents is independent of their metastatic potential; (5) shifts in metastatic potential and sensitivity to therapy may occur simultaneously but are not linked; and (6) tumor cells independently diverge to form several subpopulations with unique phenotypic profiles. ^

Behavioral and cellular analysis of classical conditioning of feeding behavior in Aplysia

In classical conditioning, an associative form of learning, animals learn to associate two stimuli. Cellular and molecular mechanisms for the induction and consolidation of associative learning and memory at the level of single cells and synaptic connections have been studied in both vertebrate and invertebrate animals. The majority of studies, however, relied on aversive stimuli to induce learning. This bias may limit the extent to which identified mechanisms generalize to other forms of associative learning and memory, such as appetitive forms. The goal of the present study was to develop a classical conditioning procedure for the marine mollusk Aplysia californica using appetitive reinforcement, and to analyze associative learning using behavioral and electrophysiological techniques. ^ Using tactile stimulation of the lips as the conditional stimulus (CS) and food as the unconditional stimulus (US) a training protocol was developed that reliably induced classical conditioning of feeding behavior. Memory persisted for at least 24 hours. The gross organization of reinforcement-mediating pathways was analyzed in additional behavioral experiments. Moreover, neurophysiological correlates of classical conditioning were identified and characterized in an in vitro preparation containing the circuitry for feeding behavior. In vitro stimulation of a nerve (AT4) that may mediate the CS during training, resulted in a greater number of buccal motor patterns (BMPs) in brains from conditioned animals, as compared to control animals. The majority of these BMPs were ingestion-like, consistent with the increased number of bites in response to the CS after classical conditioning. Moreover, classical conditioning correlated with increased excitatory synaptic input to BMP-initiating neuron B31/32, in response to stimulation of AT 4, as compared to controls. The expression of the correlates of classical conditioning identified in this study was specific to stimulation of AT 4, which is consistent the stimulus specificity that is characteristic for classical conditioning. ^ The identification of cellular correlates of classical conditioning documented here provides the basis for future, more detailed analyses of an appetitive form of associative learning and memory, that may extend the working knowledge of the cellular and molecular mechanisms for associative plasticity in general. ^

Loss of activator protein-2alpha results in overexpression of the thrombin receptor and correlates with the malignant phenotype of human melanoma

Increasing evidence demonstrates that the thrombin receptor (protease activated receptor-1, PAR-1) plays a major role in tumor invasion and contributes to the metastatic phenotype of human melanoma. We demonstrate that the metastatic potential of human melanoma cells correlates with overexpression of PAR-1. The promoter of the PAR-1 gene contains multiple putative AP-2 and Sp1 consensus elements. We provide evidence that an inverse correlation exists between the expression of AP-2 and the expression of PAR-1 in human melanoma cells. Re-expression of AP-2 in WM266-4 melanoma cells (AP-2 negative) resulted in decreased mRNA and protein expression of PAR-1 and significantly reduced the tumor potential in nude mice. ChIP analysis of the PAR-1 promoter regions bp −365 to −329 (complex 1) and bp −206 to −180 (complex 2) demonstrates that in metastatic cells Sp1 is predominantly binding to the PAR-1 promoter, while in nonmetastatic cells AP-2 is bound. In vitro analysis of complex 1 demonstrates that AP-2 and Sp1 bind to this region in a mutually exclusive manner. Transfection experiments with full-length and progressive deletions of the PAR-1 promoter luciferase constructs demonstrated that metastatic cells had increased promoter activity compared to low and nonmetastatic melanoma cells. Our data shows that exogenous AP-2 expression decreased promoter activity, while transient expression of Sp1 further activated expression of the reporter gene. Mutational analysis of complex 1 within PAR-1 luciferase constructs further demonstrates that the regulation of PAR-1 is mediated through interactions with AP-2 and Sp1. Moreover, loss of AP-2 in metastatic cells alters the AP-2 to Sp1 ratio and DNA-binding activity resulting in overexpression of PAR-1. In addition, we evaluated the expression of AP-2 and PAR-1 utilizing a tissue microarray of 93 melanocytic lesions spanning from benign nevi to melanoma metastasis. We report loss of AP-2 expression in malignant tumors compared to benign tissue while PAR-1 was expressed more often in metastatic melanoma cells than in benign melanocytes. We propose that loss of AP-2 results in increased expression of PAR-1, which in turn results in upregulation of gene products that contribute to the metastatic phenotype of melanoma. ^

Single nucleotide polymorphism (SNP) selection for genotype-phenotype association studies

With hundreds of single nucleotide polymorphisms (SNPs) in a candidate gene and millions of SNPs across the genome, selecting an informative subset of SNPs to maximize the ability to detect genotype-phenotype association is of great interest and importance. In addition, with a large number of SNPs, analytic methods are needed that allow investigators to control the false positive rate resulting from large numbers of SNP genotype-phenotype analyses. This dissertation uses simulated data to explore methods for selecting SNPs for genotype-phenotype association studies. I examined the pattern of linkage disequilibrium (LD) across a candidate gene region and used this pattern to aid in localizing a disease-influencing mutation. The results indicate that the r2 measure of linkage disequilibrium is preferred over the common D′ measure for use in genotype-phenotype association studies. Using step-wise linear regression, the best predictor of the quantitative trait was not usually the single functional mutation. Rather it was a SNP that was in high linkage disequilibrium with the functional mutation. Next, I compared three strategies for selecting SNPs for application to phenotype association studies: based on measures of linkage disequilibrium, based on a measure of haplotype diversity, and random selection. The results demonstrate that SNPs selected based on maximum haplotype diversity are more informative and yield higher power than randomly selected SNPs or SNPs selected based on low pair-wise LD. The data also indicate that for genes with small contribution to the phenotype, it is more prudent for investigators to increase their sample size than to continuously increase the number of SNPs in order to improve statistical power. When typing large numbers of SNPs, researchers are faced with the challenge of utilizing an appropriate statistical method that controls the type I error rate while maintaining adequate power. We show that an empirical genotype based multi-locus global test that uses permutation testing to investigate the null distribution of the maximum test statistic maintains a desired overall type I error rate while not overly sacrificing statistical power. The results also show that when the penetrance model is simple the multi-locus global test does as well or better than the haplotype analysis. However, for more complex models, haplotype analyses offer advantages. The results of this dissertation will be of utility to human geneticists designing large-scale multi-locus genotype-phenotype association studies. ^