Genetic risk scores and family history as predictors of schizophrenia in Nordic registers.

BACKGROUND: Family history is a long-standing and readily obtainable risk factor for schizophrenia (SCZ). Low-cost genotyping technologies have enabled large genetic studies of SCZ, and the results suggest the utility of genetic risk scores (GRS, direct assessments of inherited common variant risk). Few studies have evaluated family history and GRS simultaneously to ask whether one can explain away the other. METHODS: We studied 5959 SCZ cases and 8717 controls from four Nordic countries. All subjects had family history data from national registers and genome-wide genotypes that were processed through the quality control procedures used by the Psychiatric Genomics Consortium. Using external training data, GRS were estimated for SCZ, bipolar disorder (BIP), major depression, autism, educational attainment, and body mass index. Multivariable modeling was used to estimate effect sizes. RESULTS: Using harmonized genomic and national register data from Denmark, Estonia, Norway, and Sweden, we confirmed that family history of SCZ and GRS for SCZ and BIP were risk factors for SCZ. In a joint model, the effects of GRS for SCZ and BIP were essentially unchanged, and the effect of family history was attenuated but remained significant. The predictive capacity of a model including GRS and family history neared the minimum for clinical utility. CONCLUSIONS: Combining national register data with measured genetic risk factors represents an important investigative approach for psychotic disorders. Our findings suggest the potential clinical utility of combining GRS and family history for early prediction and diagnostic improvements.

An epigenome-wide association study meta-analysis of educational attainment.

The epigenome is associated with biological factors, such as disease status, and environmental factors, such as smoking, alcohol consumption and body mass index. Although there is a widespread perception that environmental influences on the epigenome are pervasive and profound, there has been little evidence to date in humans with respect to environmental factors that are biologically distal. Here we provide evidence on the associations between epigenetic modifications-in our case, CpG methylation-and educational attainment (EA), a biologically distal environmental factor that is arguably among the most important life-shaping experiences for individuals. Specifically, we report the results of an epigenome-wide association study meta-analysis of EA based on data from 27 cohort studies with a total of 10 767 individuals. We find nine CpG probes significantly associated with EA. However, robustness analyses show that all nine probes have previously been found to be associated with smoking. Only two associations remain when we perform a sensitivity analysis in the subset of never-smokers, and these two probes are known to be strongly associated with maternal smoking during pregnancy, and thus their association with EA could be due to correlation between EA and maternal smoking. Moreover, the effect sizes of the associations with EA are far smaller than the known associations with the biologically proximal environmental factors alcohol consumption, body mass index, smoking and maternal smoking during pregnancy. Follow-up analyses that combine the effects of many probes also point to small methylation associations with EA that are highly correlated with the combined effects of smoking. If our findings regarding EA can be generalized to other biologically distal environmental factors, then they cast doubt on the hypothesis that such factors have large effects on the epigenome.

16p11.2 Locus modulates response to satiety before the onset of obesity.

BACKGROUND: The 600 kb BP4-BP5 copy number variants (CNVs) at the 16p11.2 locus have been associated with a range of neurodevelopmental conditions including autism spectrum disorders and schizophrenia. The number of genomic copies in this region is inversely correlated with body mass index (BMI): the deletion is associated with a highly penetrant form of obesity (present in 50% of carriers by the age of 7 years and in 70% of adults), and the duplication with being underweight. Mechanisms underlying this energy imbalance remain unknown. OBJECTIVE: This study aims to investigate eating behavior, cognitive traits and their relationships with BMI in carriers of 16p11.2 CNVs. METHODS: We assessed individuals carrying a 16p11.2 deletion or duplication and their intrafamilial controls using food-related behavior questionnaires and cognitive measures. We also compared these carriers with cohorts of individuals presenting with obesity, binge eating disorder or bulimia. RESULTS: Response to satiety is gene dosage-dependent in pediatric CNV carriers. Altered satiety response is present in young deletion carriers before the onset of obesity. It remains altered in adolescent carriers and correlates with obesity. Adult deletion carriers exhibit eating behavior similar to that seen in a cohort of obesity without eating disorders such as bulimia or binge eating. None of the cognitive measures are associated with eating behavior or BMI. CONCLUSIONS: These findings suggest that abnormal satiety response is a strong contributor to the energy imbalance in 16p11.2 CNV carriers, and, akin to other genetic forms of obesity, altered satiety responsiveness in children precedes the increase in BMI observed later in adolescence.

The association between lower educational attainment and depression owing to shared genetic effects? Results in ~25,000 subjects.

An association between lower educational attainment (EA) and an increased risk for depression has been confirmed in various western countries. This study examines whether pleiotropic genetic effects contribute to this association. Therefore, data were analyzed from a total of 9662 major depressive disorder (MDD) cases and 14,949 controls (with no lifetime MDD diagnosis) from the Psychiatric Genomics Consortium with additional Dutch and Estonian data. The association of EA and MDD was assessed with logistic regression in 15,138 individuals indicating a significantly negative association in our sample with an odds ratio for MDD 0.78 (0.75-0.82) per standard deviation increase in EA. With data of 884,105 autosomal common single-nucleotide polymorphisms (SNPs), three methods were applied to test for pleiotropy between MDD and EA: (i) genetic profile risk scores (GPRS) derived from training data for EA (independent meta-analysis on ~120,000 subjects) and MDD (using a 10-fold leave-one-out procedure in the current sample), (ii) bivariate genomic-relationship-matrix restricted maximum likelihood (GREML) and (iii) SNP effect concordance analysis (SECA). With these methods, we found (i) that the EA-GPRS did not predict MDD status, and MDD-GPRS did not predict EA, (ii) a weak negative genetic correlation with bivariate GREML analyses, but this correlation was not consistently significant, (iii) no evidence for concordance of MDD and EA SNP effects with SECA analysis. To conclude, our study confirms an association of lower EA and MDD risk, but this association was not because of measurable pleiotropic genetic effects, which suggests that environmental factors could be involved, for example, socioeconomic status.

An epidemiological perspective of personalized medicine: the Estonian experience.

The Estonian Biobank and several other biobanks established over a decade ago are now starting to yield valuable longitudinal follow-up data for large numbers of individuals. These samples have been used in hundreds of different genome-wide association studies, resulting in the identification of reliable disease-associated variants. The focus of genomic research has started to shift from identifying genetic and nongenetic risk factors associated with common complex diseases to understanding the underlying mechanisms of the diseases and suggesting novel targets for therapy. However, translation of findings from genomic research into medical practice is still lagging, mainly due to insufficient evidence of clinical validity and utility. In this review, we examine the different elements required for the implementation of personalized medicine based on genomic information. First, biobanks and genome centres are required and have been established for the high-throughput genomic screening of large numbers of samples. Secondly, the combination of susceptibility alleles into polygenic risk scores has improved risk prediction of cardiovascular disease, breast cancer and several other diseases. Finally, national health information systems are being developed internationally, to combine data from electronic medical records from different sources, and also to gradually incorporate genomic information. We focus on the experience in Estonia, one of several countries with national goals towards more personalized health care based on genomic information, where the unique combination of elements required to accomplish this goal are already in place.

A genome-wide association study of anorexia nervosa.

Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome-wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2907 cases with AN from 14 countries (15 sites) and 14 860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery data sets. Seventy-six (72 independent) single nucleotide polymorphisms were taken forward for in silico (two data sets) or de novo (13 data sets) replication genotyping in 2677 independent AN cases and 8629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication data sets comprised 5551 AN cases and 21 080 controls. AN subtype analyses (1606 AN restricting; 1445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01 × 10(-7)) in SOX2OT and rs17030795 (P=5.84 × 10(-6)) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76 × 10(-)(6)) between CUL3 and FAM124B and rs1886797 (P=8.05 × 10(-)(6)) near SPATA13. Comparing discovery with replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P=4 × 10(-6)), strongly suggesting that true findings exist but our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field.

A K(ATP) channel gene effect on sleep duration: from genome-wide association studies to function in Drosophila.

Humans sleep approximately a third of their lifetime. The observation that individuals with either long or short sleep duration show associations with metabolic syndrome and psychiatric disorders suggests that the length of sleep is adaptive. Although sleep duration can be influenced by photoperiod (season) and phase of entrainment (chronotype), human familial sleep disorders indicate that there is a strong genetic modulation of sleep. Therefore, we conducted high-density genome-wide association studies for sleep duration in seven European populations (N=4251). We identified an intronic variant (rs11046205; P=3.99 × 10(-8)) in the ABCC9 gene that explains ≈5% of the variation in sleep duration. An influence of season and chronotype on sleep duration was solely observed in the replication sample (N=5949). Meta-analysis of the associations found in a subgroup of the replication sample, chosen for season of entry and chronotype, together with the discovery results showed genome-wide significance. RNA interference knockdown experiments of the conserved ABCC9 homologue in Drosophila neurons renders flies sleepless during the first 3 h of the night. ABCC9 encodes an ATP-sensitive potassium channel subunit (SUR2), serving as a sensor of intracellular energy metabolism.

Meta-analysis of genome-wide association studies for personality.

Personality can be thought of as a set of characteristics that influence people's thoughts, feelings and behavior across a variety of settings. Variation in personality is predictive of many outcomes in life, including mental health. Here we report on a meta-analysis of genome-wide association (GWA) data for personality in 10 discovery samples (17,375 adults) and five in silico replication samples (3294 adults). All participants were of European ancestry. Personality scores for Neuroticism, Extraversion, Openness to Experience, Agreeableness and Conscientiousness were based on the NEO Five-Factor Inventory. Genotype data of ≈ 2.4M single-nucleotide polymorphisms (SNPs; directly typed and imputed using HapMap data) were available. In the discovery samples, classical association analyses were performed under an additive model followed by meta-analysis using the weighted inverse variance method. Results showed genome-wide significance for Openness to Experience near the RASA1 gene on 5q14.3 (rs1477268 and rs2032794, P=2.8 × 10(-8) and 3.1 × 10(-8)) and for Conscientiousness in the brain-expressed KATNAL2 gene on 18q21.1 (rs2576037, P=4.9 × 10(-8)). We further conducted a gene-based test that confirmed the association of KATNAL2 to Conscientiousness. In silico replication did not, however, show significant associations of the top SNPs with Openness and Conscientiousness, although the direction of effect of the KATNAL2 SNP on Conscientiousness was consistent in all replication samples. Larger scale GWA studies and alternative approaches are required for confirmation of KATNAL2 as a novel gene affecting Conscientiousness.

A 5S rRNA/L5 complex is a precursor to ribosome assembly in mammalian cells.

A novel 5S RNA-protein (RNP) complex in human and mouse cells has been analyzed using patient autoantibodies. The RNP is small (approximately 7S) and contains most of the nonribosome-associated 5S RNA molecules in HeLa cells. The 5S RNA in the particle is matured at its 3' end, consistent with the results of in vivo pulse-chase experiments which indicate that this RNP represents a later step in 5S biogenesis than a previously described 5S*/La protein complex. The protein moiety of the 5S RNP has been identified as ribosomal protein L5, which is known to be released from ribosomes in a complex with 5S after various treatments of the 60S subunit. Indirect immunofluorescence indicates that the L5/5S complex is concentrated in the nucleolus. L5 may therefore play a role in delivering 5S rRNA to the nucleolus for assembly into ribosomes.

Androgen receptor gene haplotype is associated with male infertility.

The purpose of the current study was to evaluate the importance of androgen receptor (AR) gene haplotypes and polymorphic CAG/GGN microsatellites in the aetiology of male infertility. We genotyped six haplotype-tagging single nucleotide polymorphisms and CAG/GGN microsatellites of the AR gene in 112 infertile and 212 control Estonian men. A total of 13 AR haplotypes (HAP1-13) were identified, among which HAP4 was found to confer increased risk for male infertility (OR = 5.15, 95% CI = 1.75-15.15, p = 0.003). However, infertile patients and controls had similar lengths and distributions of both AR CAG (mean +/- SD number of repeats 21.1 +/- 2.5 vs. 21.2 +/- 2.3, respectively) and GGN (mean +/- SD number of repeats 22.5 +/- 1.5 vs. 22.4 +/- 1.9, respectively) repeats. In addition, HAP2 was associated with more CAG repeats (r = 1.17, p = 0.033) and HAP3 with fewer CAG repeats (r = -2.93, p < 0.001) than the major haplotype HAP1. HAP3 and HAP4 were associated with more GGN repeats (r = 1.35, p = 0.001 and r = 1.36, p = 0.002, respectively) than HAP1. In conclusion, our results implicated the AR-HAP4 gene haplotype in increased risk for male infertility, while no association was found between AR CAG/GGN microsatellites and impaired spermatogenesis.

Genome-wide analysis of nuclear magnetic resonance metabolites revealed parent-of-origin effect on triglycerides in medium very low-density lipoprotein in PTPRD gene.

AIM: The aim of the study was to explore the parent-of-origin effects (POEs) on a range of human nuclear magnetic resonance metabolites. MATERIALS & METHODS: We search for POEs in 14,815 unrelated individuals from Estonian and Finnish cohorts using POE method for the genotype data imputed with 1000 G reference panel and 82 nuclear magnetic resonance metabolites. RESULTS: Meta-analysis revealed the evidence of POE for the variant rs1412727 in PTPRD gene for the metabolite: triglycerides in medium very low-density lipoprotein. No POEs were detected for genetic variants that were previously known to have main effect on circulating metabolites. CONCLUSION: We demonstrated possibility to detect POEs for human metabolites, but the POEs are weak, and therefore it is hard to detect those using currently available sample sizes.

Whole-exome sequencing identifies a potential TTN mutation in a multiplex family with inguinal hernia.

PURPOSE: Inguinal hernia repair is one of the most common procedures in general surgery. Males are seven times more likely than females to develop a hernia and have a 27 % lifetime 'risk' of inguinal hernia repair. Several studies have demonstrated that a positive family history is an important risk factor for the development of primary inguinal hernia, which indicates that genetic factors may play important roles in the etiology of the disease. So far, the contribution of genetic factors and underlying mechanisms for inguinal hernia remain largely unknown. The aim of this study was to investigate a multiplex Estonian family with inguinal hernia across four generations. METHODS: The whole-exome sequencing was carried out in three affected family members and subsequent mutation screening using Sanger sequencing was performed in ten family members (six affected and four unaffected). RESULTS: Whole-exome sequencing in three affected family members revealed a heterozygous missense mutation c.88880A>C (p.Lys29627Thr; RefSeq NM_001256850.1) in the highly conserved myosin-binding A-band of the TTN gene. Sanger sequencing demonstrated that this mutation cosegregated with the disease in this family and was not present in ethnically matched control subjects. CONCLUSION: We report that missense variant in the A-band of TTN is the strongest candidate mutation for autosomal-dominant inguinal hernia with incomplete penetrance.

Polymorphic variation in TPMT is the principal determinant of TPMT phenotype: A meta-analysis of three genome-wide association studies.

Thiopurine-related hematotoxicity in pediatric acute lymphoblastic leukemia (ALL) and inflammatory bowel diseases has been linked to genetically defined variability in thiopurine S-methyltransferase (TPMT) activity. While gene testing of TPMT is being clinically implemented, it is unclear if additional genetic variation influences TPMT activity with consequences for thiopurine-related toxicity. To examine this possibility, we performed a genome-wide association study (GWAS) of red blood cell TPMT activity in 844 Estonian individuals and 245 pediatric ALL cases. Additionally, we correlated genome-wide genotypes to human hepatic TPMT activity in 123 samples. Only genetic variants mapping to chromosome 6, including the TPMT gene region, were significantly associated with TPMT activity (P < 5.0 × 10-8 ) in each of the three GWAS and a joint meta-analysis of 1,212 cases (top hit P = 1.2 × 10-72 ). This finding is consistent with TPMT genotype being the primary determinant of TPMT activity, reinforcing the rationale for genetic testing of TPMT alleles in routine clinical practice to individualize mercaptopurine dosage.

U82, a novel snoRNA identified from the fifth intron of human and mouse nucleolin gene.

A novel snoRNA, designated as U82, was found from the sequence analysis of the 5th intron of human and mouse nucleolin gene. The snoRNA U82 has characteristic boxes C, D and D' and 11 nucleotides (nt) antisense complementarity to the 18S rRNA. Presumably U82 functions as a guide in the methylation of residue A1678 in human 18S rRNA. Northern blot analysis with various oligodeoxynucleotide probes showed that human and mouse U82 is expressed as RNA variants with length of 70 (+/- 1) and 67 (+/- 1) nt in HeLa and mouse C127 cells. Most probably, the 70 nt variant of U82 is encoded by nucleolin gene 5th intron. The 67 nt variant of U82 could be a transcript of another gene, the genomic locus of which remains unknown.

A human cDNA encoding the homologue of NADH: ubiquinone oxidoreductase subunit B13.

A cDNA encoding the human homologue of bovine NADH:ubiquinone oxidoreductase (complex I of mitochondrial respiratory chain) subunit B13 has been isolated. The clone contains an open reading frame of 348 bp, 23 bp of 5'-untranslated sequence (UTR) and a long 3'UTR of 1088 bp. The deduced amino-acid sequence is 87% identical to bovine B13. Human B13 mRNA expression was observed in all tissues examined with highest levels in heart, skeletal muscle and brain. Southern analysis of human genomic DNA revealed the presence of multigene family.

Mutation 985A>G in the MCAD gene shows low incidence in Estonian population.

Medium-chain acyl-CoA dehydrogenase (MCAD) is an enzyme responsible for large part of mitochondrial beta-oxidation of fatty acids and therefore stays on key position of cellular energy supply. In case of its deficiency, starvation, rapid growth periods or infections may cause fatal lack of energy, especially in the first years of life. MCAD deficiency is inherited in an autosomal recessive manner and it has been shown to be rather common in some European countries (Great Britain 1 in 6,000, Switzerland 1 in 10,000). In Caucasoid populations one mutation, the 985A>G transition, causing the amino acid substitution K329E, accounts for about 90% of all mutant MCAD alleles. Here we present data about screening the Estonian population for this mutation. We analyzed the DNA from 1,098 persons from all regions of Estonia (all newborns born in one month) and found 5 heterozygotes for 985A>G, that makes the carrier frequency 1 in 220 and the frequency of possibly affected homozygotes 1 out of 193,000. No mutant alleles were found among the samples of the children, who had unclear diagnosis for death during the years 1994 and 1995.

Phenylalanine hydroxylase gene mutation R408W is present on 84% of Estonian phenylketonuria chromosomes.

Phenylalanine hydroxylase (PAH) is the enzyme which converts phenylalanine into tyrosine. In case of its deficiency, hyperphenylalaninemia is observed, which leads to phenylketonuria (PKU), a disease causing mental retardation, unless treated with a low-phenylalanine diet since early childhood. In Estonia, PKU is among the most common inherited metabolic diseases. The data from retrospective study and newborn screening show an approximate incidence of 1 in 6,000 newborns. Molecular analysis of 34 Estonian patients has revealed high genotypic homogeneity in this group, as 84% of the mutant alleles carry the R408W mutation. The high rate of this mutation in the Estonian population rises the speculation of Finno-Ugric contribution to the East European pool of mutant PAH alleles. Five more mutations-IVS12nt1, R261Q, R252W, R158Q, S349P-have been detected. The mutation detection rate was 92% among the studied patients.

High carrier frequency of the 35delG deafness mutation in European populations. Genetic Analysis Consortium of GJB2 35delG.

Congenital deafness accounts for about 1 in 1000 infants and approximately 80% of cases are inherited as an autosomal recessive trait. Recently, it has been demonstrated that connexin 26 (GJB2) gene is a major gene for congenital sensorineural deafness. A single mutation (named 35delG) was found in most recessive families and sporadic cases of congenital deafness, among Caucasoids, with relative frequencies ranging from 28% to 63%. We present here the analysis of the 35delG mutation in 3270 random controls from 17 European countries. We have detected a carrier frequency for 35delG of 1 in 35 in southern Europe and 1 in 79 in central and northern Europe. In addition, 35delG was detected in five out of 376 Jewish subjects of different origin, but was absent in other non-European populations. The study suggests either a single origin for 35delG somewhere in Europe or in the Middle East, and the possible presence of a carrier advantage together with a founder effect. The 35delG carrier frequency of 1 in 51 in the overall European population clearly indicates that this genetic alteration is a major mutation for autosomal recessive deafness in Caucasoids. This finding should facilitate diagnosis of congenital deafness and allow early treatment of the affected subjects.

Structural characterization of the mouse ribosomal protein S6-encoding gene.

The gene encoding mouse ribosomal protein (r-protein) S6 is 2.7 kb in length, and is composed of five exons. The intron positions of the mouse S6 (Rps6) coincide exactly to those of the homologous human S6 (RPS6), but the last intron present in the human is absent in the mouse gene. The latter displays higher G + C content than the RPS6, both in the overall sequenced region and at the 3rd codon position. The promoter area is highly conserved between mouse and human, and contains several putative cis-acting elements. Comparison of the intronic sequences of both genes revealed surprisingly a high degree of identity (63%) within 350 bp of the first intron. Besides the single-copy Rsp6 there are up to 15 S6 family members, most likely processed pseudogenes. Characterization of the Rps6 provides a basis to study the functions of the mammalian S6 by gene targeting.