Genome-wide association and genetic functional studies identify autism susceptibility candidate 2 gene (AUTS2) in the regulation of alcohol consumption

Alcohol consumption is a moderately heritable trait, but the genetic basis in humans is largely unknown, despite its clinical and societal importance. We report a genome-wide association study meta-analysis of approximately 2.5 million directly genotyped or imputed SNPs with alcohol consumption (gram per day per kilogram body weight) among 12 population-based samples of European ancestry, comprising 26,316 individuals, with replication genotyping in an additional 21,185 individuals. SNP rs6943555 in autism susceptibility candidate 2 gene (AUTS2) was associated with alcohol consumption at genome-wide significance (P = 4 x 10(-8) to P = 4 x 10(-9)). We found a genotype-specific expression of AUTS2 in 96 human prefrontal cortex samples (P = 0.026) and significant (P < 0.017) differences in expression of AUTS2 in whole-brain extracts of mice selected for differences in voluntary alcohol consumption. Down-regulation of an AUTS2 homolog caused reduced alcohol sensitivity in Drosophila (P < 0.001). Our finding of a regulator of alcohol consumption adds knowledge to our understanding of genetic mechanisms influencing alcohol drinking behavior.

Variation in BMPR1B, TGFRB1 and BMPR2 and control of dizygotic twinning

Genes in the TGF9 signaling pathway play important roles in the regulation of ovarian follicle growth and ovulation rate. Mutations in three genes in this pathway, growth differentiation factor 9 (GDF9), bone morphogenetic protein 15 (BMP15) and the bone morphogenetic protein receptor B 1 (BMPRB1), influence dizygotic (DZ) twinning rates in sheep. To date, only variants in GDF9 and BMP15, but not their receptors transforming growth factor ss receptor 1 (TGFBR1), bone morphogenetic protein receptor 2 (BMPR2) and BMPR1B, have been investigated with respect to their roles in human DZ twinning. We screened for rare and novel variants in TGFBR1, BMPR2 and BMPR1B in mothers of dizygotic twins (MODZT) from twin-dense families, and assessed association between genotyped and imputed variants and DZ twinning in another large sample of MODZT. Three novel variants were found: a deep intronic variant in BMPR2, and one intronic and one non-synonymous exonic variant in BMPRB1 which would result in the replacement of glutamine by glutamic acid at amino acid position 294 (p.Gln294Glu). None of these variants were predicted to have major impacts on gene function. However, the p.Gln294Glu variant changes the same amino acid as a sheep BMPR1B functional variant and may have functional consequences. Six BMPR1B variants were marginally associated with DZ twinning in the larger case-control sample, but these were no longer significant once multiple testing was taken into account. Our results suggest that variation in the TGF9 signaling pathway type II receptors has limited effects on DZ twinning rates in humans.

A quantitative-trait genome-wide association study of alcoholism risk in the community: findings and implications

BACKGROUND Given moderately strong genetic contributions to variation in alcoholism and heaviness of drinking (50% to 60% heritability) with high correlation of genetic influences, we have conducted a quantitative trait genome-wide association study (GWAS) for phenotypes related to alcohol use and dependence. METHODS Diagnostic interview and blood/buccal samples were obtained from sibships ascertained through the Australian Twin Registry. Genome-wide single nucleotide polymorphism (SNP) genotyping was performed with 8754 individuals (2062 alcohol-dependent cases) selected for informativeness for alcohol use disorder and associated quantitative traits. Family-based association tests were performed for alcohol dependence, dependence factor score, and heaviness of drinking factor score, with confirmatory case-population control comparisons using an unassessed population control series of 3393 Australians with genome-wide SNP data. RESULTS No findings reached genome-wide significance (p = 8.4 x 10(-8) for this study), with lowest p value for primary phenotypes of 1.2 x 10(-7). Convergent findings for quantitative consumption and diagnostic and quantitative dependence measures suggest possible roles for a transmembrane protein gene (TMEM108) and for ANKS1A. The major finding, however, was small effect sizes estimated for individual SNPs, suggesting that hundreds of genetic variants make modest contributions (1/4% of variance or less) to alcohol dependence risk. CONCLUSIONS We conclude that: - 1) meta-analyses of consumption data may contribute usefully to gene discovery; - 2) translation of human alcoholism GWAS results to drug discovery or clinically useful prediction of risk will be challenging, and; - 3) through accumulation across studies, GWAS data may become valuable for improved genetic risk differentiation in research in biological psychiatry (e.g., prospective high-risk or resilience studies).

GWAS of butyrylcholinesterase activity identifies four novel loci, independent effects within BCHE and secondary associations with metabolic risk factors

Serum butyrylcholinesterase (BCHE) activity is associated with obesity, blood pressure and biomarkers of cardiovascular and diabetes risk. We have conducted a genome-wide association scan to discover genetic variants affecting BCHE activity, and to clarify whether the associations between BCHE activity and cardiometabolic risk factors are caused by variation in BCHE or whether BCHE variation is secondary to the metabolic abnormalities. We measured serum BCHE in adolescents and adults from three cohorts of Australian twin and family studies. The genotypes from approximately 2.4 million single-nucleotide polymorphisms (SNPs) were available in 8791 participants with BCHE measurements. We detected significant associations with BCHE activity at three independent groups of SNPs at the BCHE locus (P = 5.8 x 10(-262), 7.8 x 10(-47), 2.9 x 10(-12)) and at four other loci: RNPEP (P = 9.4 x 10(-16)), RAPH1-ABI2 (P = 4.1 x 10(-18)), UGT1A1 (P = 4.0 x 10(-8)) and an intergenic region on chromosome 8 (P = 1.4 x 10(-8)). These loci affecting BCHE activity were not associated with metabolic risk factors. On the other hand, SNPs in genes previously associated with metabolic risk had effects on BCHE activity more often than can be explained by chance. In particular, SNPs within FTO and GCKR were associated with BCHE activity, but their effects were partly mediated by body mass index and triglycerides, respectively. We conclude that variation in BCHE activity is due to multiple variants across the spectrum from uncommon/large effect to common/small effect, and partly results from (rather than causes) metabolic abnormalities.

Search for genetic variants influencing human height

Human growth and attained height are determined by a combination of genetic and environmental effects and in modern Western societies > 80% of the observed variation in height is determined by genetic factors. Height is a fundamental human trait that is associated with many socioeconomic and psychosocial factors and health measures, however little is known of the identity of the specific genes that influence height variation in the general population. This thesis work aimed to identify the genetic variants that influence height in the general population by genome-wide linkage analysis utilizing large family samples. The study focused on analysis of three separate sets of families consisting of: 1) 1,417 individuals from 277 Finnish families (FinnHeight), 2) 8,450 individuals from 3,817 families from Australia and Europe (EUHeight) and 3) 9,306 individuals from 3,302 families from the United States (USHeight). The most significant finding in this study was found in the Finnish family sample where we a locus in the chromosomal region 1p21 was linked to adult height. Several regions showed evidence for linkage in the Australian, European and US families with 8q21 and 15q25 being the most significant. The region on 1p21 was followed up with further studies and we were able to show that the collagen 11-alpha-1 gene (COL11A1) residing at this location was associated with adult height. This association was also confirmed in an independent Finnish population cohort (Health 2000) consisting of 6,542 individuals. From this population sample, we estimated that homozygous males and females for this gene variant were 1.1 and 0.6 cm taller than the respective controls. In this thesis work we identified a gene variant in the COL11A1 gene that influences human height, although this variant alone explains only 0.1% of height variation in the Finnish population. We also demonstrated in this study that special stratification strategies such as performing sex-limited analyses, focusing on dizygous twin pairs, analyzing ethnic groups within a population separately and utilizing homogenous populations such as the Finns can improve the statistical power of finding QTL significantly. Also, we concluded from the results of this study that even though genetic effects explain a great proportion of height variance, it is likely that there are tens or even hundreds of genes with small individual effects underlying the genetic architecture of height.

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.