Genotyping for genetic association studies: Methods and applications

In this thesis, two separate single nucleotide polymorphism (SNP) genotyping techniques were set up at the Finnish Genome Center, pooled genotyping was evaluated as a screening method for large-scale association studies, and finally, the former approaches were used to identify genetic factors predisposing to two distinct complex diseases by utilizing large epidemiological cohorts and also taking environmental factors into account. The first genotyping platform was based on traditional but improved restriction-fragment-length-polymorphism (RFLP) utilizing 384-microtiter well plates, multiplexing, small reaction volumes (5 µl), and automated genotype calling. We participated in the development of the second genotyping method, based on single nucleotide primer extension (SNuPeTM by Amersham Biosciences), by carrying out the alpha- and beta tests for the chemistry and the allele-calling software. Both techniques proved to be accurate, reliable, and suitable for projects with thousands of samples and tens of markers. Pooled genotyping (genotyping of pooled instead of individual DNA samples) was evaluated with Sequenom s MassArray MALDI-TOF, in addition to SNuPeTM and PCR-RFLP techniques. We used MassArray mainly as a point of comparison, because it is known to be well suited for pooled genotyping. All three methods were shown to be accurate, the standard deviations between measurements being 0.017 for the MassArray, 0.022 for the PCR-RFLP, and 0.026 for the SNuPeTM. The largest source of error in the process of pooled genotyping was shown to be the volumetric error, i.e., the preparation of pools. We also demonstrated that it would have been possible to narrow down the genetic locus underlying congenital chloride diarrhea (CLD), an autosomal recessive disorder, by using the pooling technique instead of genotyping individual samples. Although the approach seems to be well suited for traditional case-control studies, it is difficult to apply if any kind of stratification based on environmental factors is needed. Therefore we chose to continue with individual genotyping in the following association studies. Samples in the two separate large epidemiological cohorts were genotyped with the PCR-RFLP and SNuPeTM techniques. The first of these association studies concerned various pregnancy complications among 100,000 consecutive pregnancies in Finland, of which we genotyped 2292 patients and controls, in addition to a population sample of 644 blood donors, with 7 polymorphisms in the potentially thrombotic genes. In this thesis, the analysis of a sub-study of pregnancy-related venous thromboses was included. We showed that the impact of factor V Leiden polymorphism on pregnancy-related venous thrombosis, but not the other tested polymorphisms, was fairly large (odds ratio 11.6; 95% CI 3.6-33.6), and increased multiplicatively when combined with other risk factors such as obesity or advanced age. Owing to our study design, we were also able to estimate the risks at the population level. The second epidemiological cohort was the Helsinki Birth Cohort of men and women who were born during 1924-1933 in Helsinki. The aim was to identify genetic factors that might modify the well known link between small birth size and adult metabolic diseases, such as type 2 diabetes and impaired glucose tolerance. Among ~500 individuals with detailed birth measurements and current metabolic profile, we found that an insertion/deletion polymorphism of the angiotensin converting enzyme (ACE) gene was associated with the duration of gestation, and weight and length at birth. Interestingly, the ACE insertion allele was also associated with higher indices of insulin secretion (p=0.0004) in adult life, but only among individuals who were born small (those among the lowest third of birth weight). Likewise, low birth weight was associated with higher indices of insulin secretion (p=0.003), but only among carriers of the ACE insertion allele. The association with birth measurements was also found with a common haplotype of the glucocorticoid receptor (GR) gene. Furthermore, the association between short length at birth and adult impaired glucose tolerance was confined to carriers of this haplotype (p=0.007). These associations exemplify the interaction between environmental factors and genotype, which, possibly due to altered gene expression, predisposes to complex metabolic diseases. Indeed, we showed that the common GR gene haplotype associated with reduced mRNA expression in thymus of three individuals (p=0.0002).

Activator Protein-1 and Epidermal Growth Factor Receptor Interplay : In Vivo & In Vitro Studies

Critical cellular decisions such as should the cell proliferate, migrate or differentiate, are regulated by stimulatory signals from the extracellular environment, like growth factors. These signals are transformed to cellular responses through their binding to specific receptors present at the surface of the recipient cell. The epidermal growth factor receptor (EGF-R/ErbB) pathway plays key roles in governing these signals to intracellular events and cell-to-cell communication. The EGF-R forms a signaling network that participates in the specification of cell fate and coordinates cell proliferation. Ligand binding triggers receptor dimerization leading to the recruitment of kinases and adaptor proteins. This step simultaneously initiates multiple signal transduction pathways, which result in activation of transcription factors and other target proteins, leading to cellular alterations. It is known that mutations of EGF-R or in the components of these pathways, such as Ras and Raf, are commonly involved in human cancer. The four best characterized signaling pathways induced by EGF-R are the mitogen-activated protein kinase cascades (MAPKs), the lipid kinase phosphatidylinositol 3 kinase (PI3K), a group of transcription factors called Signal Transducers and Activator of Transcription (STAT), and the phospholipase Cγ; (PLCγ) pathways. The activation of each cascade culminates in kinase translocation to the nucleus to stimulate various transcription factors including activator protein 1 (AP-1). AP-1 family proteins are basic leucine zipper (bZIP) transcription factors that are implicated in the regulation of a variety of cellular processes (proliferation and survival, growth, differentiation, apoptosis, cell migration, transformation). Therefore, the regulation of AP-1 activity is critical for the decision of cell fate and their deregulated expression is widely associated with many types of cancers, such as breast and prostate cancers. The aims of this study were to characterize the roles of EGF-R signaling during normal development and malignant growth in vitro and in vivo using different cell lines and tissue samples. We show here that EGF-R regulates cell proliferation but is also required for regulation of AP-1 target gene expression in fibroblasts in a MAP-kinase mediated manner. Furthermore, EGF-R signaling is essential for enterocyte proliferation and migration during intestinal maturation. EGF-R signaling network, especially PI3-K-Akt pathway mediated AP-1 activity is involved in cellular survival in response to ionizing radiation. Taken together, these results elucidate the connection of EGF-R and AP-1 in various cellular contexts and show their importance in the regulation of cellular behaviour presenting new treatment cues for intestinal perforations and cancer therapy.

Search For Susceptibility Genes In Schizophrenia

Schizophrenia is a severe mental disorder affecting 0.4-1% of the population worldwide. It is characterized by impairments in the perception of reality and by significant social or occupational dysfunction. The disorder is one of the major contributors to the global burden of diseases. Studies of twins, families, and adopted children point to strong genetic components for schizophrenia, but environmental factors also play a role in the pathogenesis of disease. Molecular genetic studies have identified several potential positional candidate genes. The strongest evidence for putative schizophrenia susceptibility loci relates to the genes encoding dysbindin (DTNBP1) and neuregulin (NRG1), but studies lack impressive consistency in the precise genetic regions and alleles implicated. We have studied the role of three potential candidate genes by genotyping 28 single nucleotide polymorphisms in the DNTBP1, NRG1, and AKT1 genes in a large schizophrenia family sample consisting of 441 families with 865 affected individuals from Finland. Our results do not support a major role for these genes in the pathogenesis of schizophrenia in Finland. We have previously identified a region on chromosome 5q21-34 as a susceptibility locus for schizophrenia in a Finnish family sample. Recently, two studies reported association between the γ-aminobutyric acid type A receptor cluster of genes in this region and one study showed suggestive evidence for association with another regional gene encoding clathrin interactor 1 (CLINT1, also called Epsin 4 and ENTH). To further address the significance of these genes under the linkage peak in the Finnish families, we genotyped SNPs of these genes, and observed statistically significant association of variants between GABRG2 and schizophrenia. Furthermore, these variants also seem to affect the functioning of the working memory. Fetal events and obstetric complications are associated with schizophrenia. Rh incompatibility has been implicated as a risk factor for schizophrenia in several epidemiological studies. We conducted a family-based candidate-gene study that assessed the role of maternal-fetal genotype incompatibility at the RhD locus in schizophrenia. There was significant evidence for an RhD maternal-fetal genotype incompatibility, and the risk ratio was estimated at 2.3. This is the first candidate-gene study to explicitly test for and provide evidence of a maternal-fetal genotype incompatibility mechanism in schizophrenia. In conclusion, in this thesis we found evidence that one GABA receptor subunit, GABRG2, is significantly associated with schizophrenia. Furthermore, it also seems to affect to the functioning of the working memory. In addition, an RhD maternal-fetal genotype incompatibility increases the risk of schizophrenia by two-fold.

Genetic Heterogeneity in Autism Spectrum Disorders in a Population Isolate

Positional cloning has enabled hypothesis-free, genome-wide scans for genetic factors contributing to disorders or traits. Traditionally linkage analysis has been used to identify regions of interest, followed by meticulous fine mapping and candidate gene screening using association methods and finally sequencing of regions of interest. More recently, genome-wide association analysis has enabled a more direct approach to identify specific genetic variants explaining a part of the variance of the phenotype of interest. Autism spectrum disorders (ASDs) are a group of childhood onset neuropsychiatric disorders with shared core symptoms but varying severity. Although a strong genetic component has been established in ASDs, genetic susceptibility factors have largely eluded characterization. Here, we have utilized modern molecular genetic methods combined with the advantages provided by the special population structure in Finland to identify genetic risk factors for ASDs. The results of this study show that numerous genetic risk factors exist for ASDs even within a population isolate. Stratification based on clinical phenotype resulted in encouraging results, as previously identified linkage to 3p14-p24 was replicated in an independent family set of families with Asperger syndrome, but no other ASDs. Fine-mapping of the previously identified linkage peak for ASDs at 3q25-q27 revealed association between autism and a subunit of the 5-hydroxytryptamine receptor 3C (HTR3C). We also used dense, genome-wide single nucleotide polymorphism (SNP) data to characterize the population structure of Finns. We observed significant population substructure which correlates with the known history of multiple consecutive bottle-necks experienced by the Finnish population. We used this information to ascertain a genetically homogenous subset of autism families to identify possible rare, enriched risk variants using genome-wide SNP data. No rare enriched genetic risk factors were identified in this dataset, although a subset of families could be genealogically linked to form two extended pedigrees. The lack of founder mutations in this isolated population suggests that the majority of genetic risk factors are rare, de novo mutations unique to individual nuclear families. The results of this study are consistent with others in the field. The underlying genetic architecture for this group of disorders appears highly heterogeneous, with common variants accounting for only a subset of genetic risk. The majority of identified risk factors have turned out to be exceedingly rare, and only explain a subset of the genetic risk in the general population in spite of their high penetrance within individual families. The results of this study, together with other results obtained in this field, indicate that family specific linkage, homozygosity mapping and resequencing efforts are needed to identify these rare genetic risk factors.

Associations of gene polymorphisms and nutrition with calcium homeostasis and bone mineral density : Studies on skeletal nutrigenetics

Heredity explains a major part of the variation in calcium homeostasis and bone strength, and the susceptibility to osteoporosis is polygenetically regulated. Bone phenotype results from the interplay between lifestyle and genes, and several nutritional factors modulate bone health throughout life. Thus, nutrigenetics examining the genetic variation in nutrient intake and homeostatic control is an important research area in the etiology of osteoporosis. Despite continuing progress in the search for candidate genes for osteoporosis, the results thus far have been inconclusive. The main objective of this thesis was to investigate the associations of lactase, vitamin D receptor (VDR), calcium sensing receptor (CaSR) and parathyroid hormone (PTH) gene polymorphisms and lifestyle factors and their interactions with bone health in Finns at varying stages of the skeletal life span. Markers of calcium homeostasis and bone remodelling were measured from blood and urine samples. Bone strength was measured at peripheral and central bone sites. Lifestyle factors were assessed with questionnaires and interviews. Genetic lactase non-persistence (the C/C-13910 genotype) was associated with lower consumption of milk from childhood, predisposing females in particular to inadequate calcium intake. Consumption of low-lactose milk and milk products was shown to decrease the risk for inadequate calcium intake. In young adulthood, bone loss was more common in males than in females. Males with the lactase C/C-13910 genotype may be more susceptible to bone loss than males with the other lactase genotypes, although calcium intake predicts changes in bone mass more than the lactase genotype. The BsmI and FokI polymorphisms of the VDR gene were associated with bone mass in growing adolescents, but the associations weakened with age. In young adults, the A986S polymorphism of the calcium sensing receptor gene was associated with serum ionized calcium concentrations, and the BstBI polymorphism of the parathyroid gene was related to bone strength. The FokI polymorphism and sodium intake showed an interaction effect on urinary calcium excretion. A novel gene-gene interaction between the VDR FokI and PTH BstBI gene polymorphisms was found in the regulation of PTH secretion and urinary calcium excretion. Further research should be carried out with more number of Finns at varying stages of the skeletal life span and more detailed measurements of bone strength. Research should concern mechanisms by which genetic variants affect calcium homeostasis and bone strength, and the role of diet-gene and gene-gene interactions in the pathogenesis of osteoporosis.

Genetics of T cell co-stimulatory receptors -CD28, CTLA4, ICOS and PDCD1 in immunity and transplantation

Co-stimulatory signals are essential for the activation of naïve T cells and productive immune response. Naïve T cells receive first, antigen-specific signal through T cell receptor. Co-stimulatory receptors provide the second signal which can be either activating or inhibitory. The balance between signals determines the outcome of an immune response. CD28 is crucial for T cell activation; whereas cytotoxic T lymphocyte associated antigen 4 (CTLA4) mediates critical inhibitory signal. Inducible co-stimulator (ICOS) augments cytokine expression and plays role in immunoglobulin class switching. Programmed cell death 1 (PDCD1) acts as negative regulator of T cell proliferation and cytokine responses. The co-stimulatory receptor pathways are potentially involved in self-tolerance and thus, they provide a promising therapeutic strategy for autoimmune diseases and transplantation. The genes encoding CD28, CTLA4 and ICOS are located adjacently in the chromosome region 2q33. The PDCD1 gene maps further, to the region 2q37. CTLA4 and PDCD1 are associated with the risk of a few autoimmune diseases. There is strong linkage disequilibrium (LD) on the 2q33 region; the whole gene of CD28 exists in its own LD block but CTLA4 and the 5' part of ICOS are within a same LD block. The 3' part of ICOS and PDCD1 are in their own separate LD blocks. Extended haplotypes covering the 2q33 region can be identified. This study focuses on immune related conditions like coeliac disease (CD) which is a chronic inflammatory disease with autoimmune features. Immunoglobulin A deficiency (IgAD) belongs to the group of primary antibody deficiencies characterised by reduced levels of immunoglobulins. IgAD co-occurs often with coeliac disease. Renal transplantation is needed in the end stage kidney diseases. Transplantation causes strong immune response which is tried to suppress with drugs. All these conditions are multifactorial with complex genetic background and multiple environmental factors affecting the outcome. We have screened ICOS for polymorphisms by sequencing the exon regions. We detected 11 new variants and determined their frequencies in Finnish population. We have measured linkage disequilibrium on the 2q33 region in Finnish as well as other European populations and observed conserved haplotypes. We analysed genetic association and linkage of the co-stimulatory receptor gene region aiming to study if it is a common risk locus for immune diseases. The 2q33 region was replicated to be linked to coeliac disease in Finnish population and CTLA4-ICOS haplotypes were found to be associated with CD and IgAD being the first non-HLA risk locus common for CD and immunodeficiencies. We also showed association between ICOS and the outcome of kidney transplantation. Our results suggest new evidence for CTLA4-ICOS gene region to be involved in susceptibility of coeliac disease. The earlier published contradictory association results can be explained by involvement of both CTLA4 and ICOS in disease susceptibility. The pattern of variants acting together rather than a single polymorphism may confer the disease risk. These genes may predispose also to immunodeficiencies as well as decreased graft survival and delayed graft function. Consequently, the present study indicates that like the well established HLA locus, the co-stimulatory receptor genes predispose to variety of immune disorders.

GABAa Receptor Mediated Signalling in the Brain: Inhibition, Shunting and Excitation

Within central nervous system, the simple division of chemical synaptic transmission to depolarizing excitation mediated by glutamate and hyperpolarizing inhibition mediated by γ-amino butyric acid (GABA), is evidently an oversimplification. The GABAa receptor (GABAaR) mediated responses can be of opposite sign within a single resting cell, due to the compartmentalized distribution of cation chloride cotransporters (CCCs). The K+/Cl- cotransporter 2 (KCC2), member of the CCC family, promotes K+ fuelled Cl- extrusion and sets the reversal potential of GABA evoked anion currents typically slightly below the resting membrane potential. The interesting ionic plasticity property of GABAergic signalling emerges from the short-term and long-term alterations in the intraneuronal concentrations of GABAaR permeable anions (Cl- and HCO3-). The short-term effects arise rapidly (in the time scale of hundreds of milliseconds) and are due to the GABAaR activation dependent shifts in anion gradients, whereas the changes in expression, distribution and kinetic regulation of CCCs are underlying the long-term effects, which may take minutes or even hours to develop. In this Thesis, the differences in the reversal potential of GABAaR mediated responses between dopaminergic and GABAergic cell types, located in the substantia nigra, were shown to be attributable to the differences in the chloride extrusion mechanisms. The stronger inhibitory effect of GABA on GABAergic neurons was due to the cell type specific expression of KCC2 whereas the KCC2 was absent from dopaminergic neurons, leading to a less prominent inhibition brought by GABAaR activation. The levels of KCC2 protein exhibited activity dependent alterations in hippocampal pyramidal neurons. Intense neuronal activity, leading to a massive release of brain derived neurotrophic factor (BDNF) in vivo, or applications of tyrosine receptor kinase B (TrkB) agonists BDNF or neurotrophin-4 in vitro, were shown to down-regulate KCC2 protein levels which led to a reduction in the efficacy of Cl- extrusion. The GABAergic transmission is interestingly involved in an increase of extracellular K+ concentration. A substantial increase in interstitial K+ tends to depolarize the cell membrane. The effects that varying ion gradients had on the generation of biphasic GABAaR mediated responses were addressed, with particular emphasis on the novel idea that the K+/Cl- extrusion via KCC2 is accelerated in response to a rapid accumulation of intracellular Cl-. The KCC2 inhibitor furosemide produced a large reduction in the GABAaR dependent extracellular K+ transients. Thus, paradoxically, both the inefficient KCC2 activity (via increased intracellular Cl-) and efficient KCC2 activity (via increased extracellular K+) may promote excitation.

ORP2 is a sterol receptor that regulates cellular lipid metabolism

ORP2 is a member of mammalian oxysterol binding protein (OSBP)-related protein/gene family (ORPs), which is found in almost every eukaryotic organism. ORPs have been suggested to participate in the regulation of cellular lipid metabolism, vesicle trafficking and cellular signaling. ORP2 is a cytosolic protein that is ubiquitously expressed and most abundant in the brain. In previous studies employing stable cell lines with constitutive ORP2 overexpression ORP2 was shown to affect cellular cholesterol metabolism. The aim of this study was to characterize the properties and function of ORP2 further. ORP2 ligands were searched for among sterols and phosphoinositides using purified ORP2 and in vitro binding assays. As expected, ORP2 bound several oxysterols and cholesterol, the highest affinity ligand being 22(R)hydroxycholesterol. In addition, affinity for anionic membrane phospholipids, phosphoinositides was observed, which may assist in the membrane targeting of ORP2. Intracellular localization of ORP2 was also investigated. ORP2 was observed on the surface of cytoplasmic lipid droplets, which are storage organelles for neutral lipids. Lipid droplet targeting of ORP2 was inhibited when 22(R)hydroxycholesterol was added to the cells or when the N-terminal FFAT-motif of ORP2 was mutated, suggesting that oxysterols and the N-terminus of ORP2 regulate the localization and the function of ORP2. The role of ORP2 in cellular lipid metabolism was studied using HeLa cell lines that can be induced to overexpress ORP2. Overexpression of ORP2 was shown to enhance cholesterol efflux from the cells resulting in a decreased amount of cellular free cholesterol. ORP2 overexpressing cells responded to the loss of cholesterol by upregulating cholesterol synthesis and uptake. Intriguingly, also cholesterol esterification was increased in ORP2 overexpressing cells. These results may be explained by the ability of ORP2 to bind and thus transport cholesterol, which most likely leads to changes in cholesterol metabolism when ORP2 is overexpressed. ORP2 function was further investigated by silencing the endogenous ORP2 expression with short interfering RNAs (siRNA) in A431 cells. Silencing of ORP2 led to a delayed break-down of triglycerides under lipolytic conditions and an increased amount of cholesteryl esters in the presence of excess triglycerides. Together these results suggest that ORP2 is a sterol-regulated protein that functions on the surface of cytoplasmic lipid droplets to regulate the metabolism of triglycerides and cholesteryl esters. Although the exact mode of ORP2 action still remains unclear, this study serves as a good basis to investigate the molecular mechanisms and possible cell type specific functions of ORP2.