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.

Theoretical studies on the conformations of higher gangliosides

The possible conformations of higher gangliosides (GD3, GT1a. GT1b, GQ1b) have been determined by computing their potential energy using semi-empirical potential functions. The favoured conformation of the disialic acid fragment in these gangliosides is independent of its position (internal or terminal). The favoured conformations of these gangliosides have also been correlated to their biological activity. The results suggest that tetanus toxin and sendai virus may have a large binding site which can accommodate at least four sugar residues.

Theoretical studies on the conformations of higher gangliosides

The possible conformations of higher gangliosides (GD3, GT1a. GT1b, GQ1b) have been determined by computing their potential energy using semi-empirical potential functions. The favoured conformation of the disialic acid fragment in these gangliosides is independent of its position (internal or terminal). The favoured conformations of these gangliosides have also been correlated to their biological activity. The results suggest that tetanus toxin and sendai virus may have a large binding site which can accommodate at least four sugar residues.

Fluorescence And Nmr-Studies Of The Binding Of Cholinergic Fluorescent-Probes To Horse Serum-Cholinesterase

The interaction of the cholinergic fluorescent probes, 1-(5-dimethyl-aminoaphthalene-1-sulfonamido) ethane-2-trimethylammonium perchlorate, 1-(5-dimethylaminonaphthalene-1-sulfonamido) pentane-5-trimethylammonium tartarate and 1-(5-dimethylaminonaphthalene-1-sulfonamido) decane-10- trimethylammonium tartarate with horse serum cholinesterase has been examined by fluorescence and n.m.r. methods. Fluorescence titrations show binding of the decane derivative to two sites on the protein whereas the lower homologs bind largely to one site. Active site inhibitors like curbamylcholine and decamethonium abolish binding of the decane derivative to the high affinity site. The inhibitors are largely without effect on the binding of the lower homologs. N.m.r. studies clearly establish immobilization of both ends of the molecule on binding in the case of the decane derivative, whereas in the lower homologs the dimethylamino group on the naphthalene ring is significantly more affected in the presence of enzyme. The probes are effective inhibitors of the enzyme with the decane derivative being two orders of magnitude more effective than its lower homologs. Based on the n.m.r., fluorescence and inhibition studies, a model for probe binding to the enzyme is advanced. It appears that the decane derivative binds with high affinity to the catalytic anionic site while the lower affinity site is assigned to a peripheral anionic site. The lower homologs probe only the peripheral site. A comparison of fluorescence, n.m.r. and inhibition studies with acetylcholinesterases from electric eel and bovine erythrocytes is presented.

Studies of phosphazenes. Part 10. Spirocyclic derivatives of cyclotriphosphazatrienes

The reactions of hexachlorocyclotriphosphazatriene, N3P3Cl6, and its geminal bis-t-butylamino- and diphenyl derivatives, with ethylenediamine and ethanolamine are reported. In each case, both chlorine atoms attached to the same phosphorus atoms are replaced, giving rise to spirocyclic derivatives. A small quantity of a bis spirocyclic derivative, N3P3(NHCH2CH2O)2Cl2, is also obtained; this compound occurs in both cis and trans forms. Attempts to prepare fully substituted tris spirocyclic derivatives have been unsuccessful and only resinous materials were obtained. The 1H and 31P n.m.r. spectra of the products are discussed.

1H and 13C NMR spectral studies of C-2 substituted isomeric exo- and endo-5-methyl-bicyclo[3.2.1]octane-6,8-diones

A detailed analysis of the 1H and 13C NMR spectra of C-2 aryl and alkyl/desalkyl substituted isomeric exo- and endo-5-methylbicyclo[3.2.1]octane-6,8-diones is presented. The chemical shift of the C-5 angular methyl, the C-2 alkyl/olefinic (C-10)/C-2 methine protons, the aromatic proton shieldings and the characteristic AMX and ABX spectral pattern of the ketomethylene and bridgehead protons were found to be sensitive to the phenyl ring orientation (anisotropy). These distinctive features could be used for configurational distinction for this class of compounds. With increasing ortho-methoxy substitution on the phenyl ring, considerable deshilelding of the bridgehead proton was observed (ca. 0.6 ppm). Absence of the C-2 alkyl group in the desalkyl isomers resulted in substantial changes in the chemical shifts of different protons. A study of the NMR spectra of the corresponding bicyclic compounds with C-2 methoxy/hydroxy substitution instead of the aryl group revealed that the anisotropy of the phenyl ring and the electronegative oxygen substituents have opposite effects. The 13C NMR spectral assignment of each carbon resonance of C-2 aryl and alkyl/desalkyl substituted isomeric exo- and endo-5-methylbicyclo[3.2.1]octane-6,8-diones and the corresponding C-2 methoxy/hydroxy/chloro and methyl bicyclic compounds are reported. Additional ortho-methoxy substitution on the phenyl ring was found to produce considerable high field shifts of the C-10 and C-1 carbon resonances. A high-field shift was observed for the C-6 and C-8 carbonyl carbons, presumably due to 1,3-dicarbonyl interactions. The chemical shifts of C-1 aromatic, C-10 alkyl and C-2 carbons, which are sensitive to exo/endo isomerism, could be utilized in differentiating a pair of isomers.

A versatile gene-based test for genome-wide association studies

We have derived a versatile gene-based test for genome-wide association studies (GWAS). Our approach, called VEGAS (versatile gene-based association study), is applicable to all GWAS designs, including family-based GWAS, meta-analyses of GWAS on the basis of summary data, and DNA-pooling-based GWAS, where existing approaches based on permutation are not possible, as well as singleton data, where they are. The test incorporates information from a full set of markers (or a defined subset) within a gene and accounts for linkage disequilibrium between markers by using simulations from the multivariate normal distribution. We show that for an association study using singletons, our approach produces results equivalent to those obtained via permutation in a fraction of the computation time. We demonstrate proof-of-principle by using the gene-based test to replicate several genes known to be associated on the basis of results from a family-based GWAS for height in 11,536 individuals and a DNA-pooling-based GWAS for melanoma in approximately 1300 cases and controls. Our method has the potential to identify novel associated genes; provide a basis for selecting SNPs for replication; and be directly used in network (pathway) approaches that require per-gene association test statistics. We have implemented the approach in both an easy-to-use web interface, which only requires the uploading of markers with their association p-values, and a separate downloadable application.

A simple and fast two-locus quality control test to detect false positives due to batch effects in genome-wide association studies

The impact of erroneous genotypes having passed standard quality control (QC) can be severe in genome-wide association studies, genotype imputation, and estimation of heritability and prediction of genetic risk based on single nucleotide polymorphisms (SNP). To detect such genotyping errors, a simple two-locus QC method, based on the difference in test statistic of association between single SNPs and pairs of SNPs, was developed and applied. The proposed approach could detect many problematic SNPs with statistical significance even when standard single SNP QC analyses fail to detect them in real data. Depending on the data set used, the number of erroneous SNPs that were not filtered out by standard single SNP QC but detected by the proposed approach varied from a few hundred to thousands. Using simulated data, it was shown that the proposed method was powerful and performed better than other tested existing methods. The power of the proposed approach to detect erroneous genotypes was approximately 80% for a 3% error rate per SNP. This novel QC approach is easy to implement and computationally efficient, and can lead to a better quality of genotypes for subsequent genotype-phenotype investigations.

Thermodynamic and kinetic studies on saccharide binding to soya-bean agglutinin.

The fluorescence of N-dansylgalactosamine [N-(5-dimethylaminonaphthalene-1-sulphonyl)galactosamine] was enhanced 11-fold with a 25 nm blue-shift in the emission maximum upon binding to soya-bean agglutinin (SBA). This change was used to determine the association constants and thermodynamic parameters for this interaction. The association constant of 1.51 X 10(6) M-1 at 20 degrees C indicated a very strong binding, which is mainly due to a relatively small entropy value, as revealed by the thermodynamic parameters: delta G = -34.7 kJ X mol-1, delta H = -37.9 kJ X mol-1 and delta S = -10.9 J X mol-1 X K-1. The specific binding of this sugar to SBA shows that the lectin can accommodate a large hydrophobic substituent on the C-2 of galactose. Binding of non-fluorescent ligands, studied by monitoring the fluorescence changes when they are added to a mixture of SBA and N-dansylgalactosamine, indicates that a hydrophobic substituent at the anomeric position increases the affinity of the interaction. The C-6 hydroxy group also stabilizes the binding considerably. Kinetics of binding of N-dansylgalactosamine to SBA studied by stopped-flow spectrofluorimetry are consistent with a single-step mechanism and yielded k+1 = 2.4 X 10(5) M-1 X s-1 and k-1 = 0.2 s-1 at 20 degrees C. The activation parameters indicate an enthalpicly controlled association process.

Chemical Modification Studies On Ricinus-Communis (Castor Bean) Agglutinin

Ricinus communis agglutinin was subjected to various chemical treatments and the effect on its hemagglutinating and saccharide-binding properties was studied. Acetylation, succinylation and citraconylation led to a complete loss in the activity of the agglutinin, whereas reductive methylation had no effect on the activity, showing that charged amino groups were involved in the hemagglutinating and saccharide-binding activity of Ricinus agglutinin. Modification of tryptophyl, arginyl and carboxyl-group-containing residues did not lead to any loss in the activity of the agglutinin. Acetylation of tyrosyl groups with N-acetylimidazole strongly reduced the hemagglutinating and saccharide-binding property of Ricinus agglutinin. The loss in activity was restored on deacetylation of the tyrosyl groups. Modification of tyrosyl residues also led to a change in the immunological properties of the agglutinin. The initial rate of modification of tyrosyl and amino groups and the concomitant loss of activity was reduced in the presence of lactose.