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.

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.

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.

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.

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 missense mutation in DUSP6 is associated with Class III malocclusion.

Class III malocclusion is a common dentofacial phenotype with a variable prevalence according to ethnic background. The etiology of Class III malocclusion has been attributed mainly to interactions between susceptibility genes and environmental factors during the morphogenesis of the mandible and maxilla. Class III malocclusion shows familial recurrence, and family-based studies support a predominance of an autosomal-dominant mode of inheritance. We performed whole-exome sequencing on five siblings from an Estonian family affected by Class III malocclusion. We identified a rare heterozygous missense mutation, c.545C>T (p.Ser182Phe), in the DUSP6 gene, a likely causal variant. This variant co-segregated with the disease following an autosomal-dominant mode of inheritance with incomplete penetrance. Transcriptional activation of DUSP6 has been presumed to be regulated by FGF/FGFR and MAPK/ERK signaling during fundamental processes at early stages of skeletal development. Several candidate genes within a linkage region on chromosome 12q22-q23--harboring DUSP6--are implicated in the regulation of maxillary or mandibular growth. The current study reinforces that the 12q22-q23 region is biologically relevant to craniofacial development and may be genetically linked to the Class III malocclusion.

Non-syndromic tooth agenesis associated with a nonsense mutation in ectodysplasin-A (EDA).

Mutations in the ectodysplasin-A (EDA) gene have been generally associated with X-linked hypohidrotic ectodermal dysplasia (XLHED). Recently, missense mutations in EDA have been reported to cause familial non-syndromic tooth agenesis. In this study, we report a novel EDA mutation in an Estonian family segregating non-syndromic tooth agenesis with variable expressivity. Affected individuals had no associated defects in other ectodermal organs. Using whole-exome sequencing, we identified a heterozygous nonsense mutation c.874G>T (p.Glu292X) in the TNF homology domain of EDA in all affected female patients. This protein-altering variant arose de novo, and the potentially causative allele was transmitted to affected offspring from the affected mother. We suggest that the dental phenotype variability described in heterozygous female carriers of EDA mutation may occur because of the differential pattern of X-chromosome inactivation, which retains reduced levels of EDA-receptor signaling in tissues involved in tooth morphogenesis. This results in selective tooth agenesis rather than XLHED phenotype. The present study broadens the mutation spectrum for this locus and demonstrates that EDA mutations may result in non-syndromic tooth agenesis in heterozygous females.

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.

Replication and meta-analysis of TMEM132D gene variants in panic disorder.

A recent genome-wide association study in patients with panic disorder (PD) identified a risk haplotype consisting of two single-nucleotide polymorphisms (SNPs) (rs7309727 and rs11060369) located in intron 3 of TMEM132D to be associated with PD in three independent samples. Now we report a subsequent confirmation study using five additional PD case-control samples (n = 1670 cases and n = 2266 controls) assembled as part of the Panic Disorder International Consortium (PanIC) study for a total of 2678 cases and 3262 controls in the analysis. In the new independent samples of European ancestry (EA), the association of rs7309727 and the risk haplotype rs7309727-rs11060369 was, indeed, replicated, with the strongest signal coming from patients with primary PD, that is, patients without major psychiatric comorbidities (n = 1038 cases and n = 2411 controls). This finding was paralleled by the results of the meta-analysis across all samples, in which the risk haplotype and rs7309727 reached P-levels of P = 1.4e-8 and P = 1.1e-8, respectively, when restricting the samples to individuals of EA with primary PD. In the Japanese sample no associations with PD could be found. The present results support the initial finding that TMEM132D gene contributes to genetic susceptibility for PD in individuals of EA. Our results also indicate that patient ascertainment and genetic background could be important sources of heterogeneity modifying this association signal in different populations.

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.

Meta-analysis of genome-wide association studies identifies common variants in CTNNA2 associated with excitement-seeking.

The tendency to seek stimulating activities and intense sensations define excitement-seeking, a personality trait akin to some aspects of sensation-seeking. This trait is a central feature of extraversion and is a component of the multifaceted impulsivity construct. Those who score high on measures of excitement-seeking are more likely to smoke, use other drugs, gamble, drive recklessly, have unsafe/unprotected sex and engage in other risky behaviors of clinical and social relevance. To identify common genetic variants associated with the Excitement-Seeking scale of the Revised NEO Personality Inventory, we performed genome-wide association studies in six samples of European ancestry (N=7860), and combined the results in a meta-analysis. We identified a genome-wide significant association between the Excitement-Seeking scale and rs7600563 (P=2 × 10(-8)). This single-nucleotide polymorphism maps within the catenin cadherin-associated protein, alpha 2 (CTNNA2) gene, which encodes for a brain-expressed α-catenin critical for synaptic contact. The effect of rs7600563 was in the same direction in all six samples, but did not replicate in additional samples (N=5105). The results provide insight into the genetics of excitement-seeking and risk-taking, and are relevant to hyperactivity, substance use, antisocial and bipolar disorders.

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.

[The Estonian Genome Project in the context of European genome research]. / Das estnische Genomprojekt im Kontext der europäischen Genomforschung.

It is the aim of the Estonian Genome Project to establish a database which compiles phenotype and genotype data of a large part of the Estonian population. The Gene Bank will only be used for scientific and public health research. Researchers hope that it will help identify disease genes and prepare the ground for the personalized medicine of the future. Additionally, the project will improve Estonian's international competitiveness in high technology and have a strong educational effect on the population. The legal framework, the Human Genes Research Act, was passed by the Estonian parliament in December 2000. It had been drafted by a group of experts who took into account all available international guidance documents on genetic research. With the pilot project involving three selected regions successfully finished, the main project has now started.

Sequence variations in the NDUFA1 gene encoding a subunit of complex I of the respiratory chain.

NDUFA1 is one of the 36 nuclear genes encoding subunits of the mitochondrial complex I involved in the respiratory chain. The human NDUFA1 has been cloned, completely sequenced and mapped to Xq24. In the present study, we searched for sequence variations in NDUFA1 as causative defects in complex I deficiency using genomic DNA of 152 patients with various clinical phenotypes. The patient sample consisted of 54 patients (46 male and 8 female) with Leber heriditary optic neuropathy (LHON) from 48 unrelated families from Germany and 98 patients (72 male and 26 female) with biochemically proven complex I deficiency including Leigh syndrome. Patient DNA was used to amplify all three exons, including the exon/intron boundaries and the promoter region of NDUFA1 for heteroduplex analysis and direct sequencing. In the 152 patients tested, no mutation was found that could be related to any of the disease phenotypes included. However, three single-nucleotide polymorphisms (SNPs) located in the promoter region (SNP G/C at nt -71 and SNP T/C at nt -189) and in intron 1 (SNP T/G nt 1454) were discovered. Allele frequencies of the SNPs were estimated in a German and Estonian control population and compared to complex I-deficient patients. There was no significant difference between the control population, the LHON patients, or the severely affected patients with complex I deficiency, excluding an association of the polymorphisms with the diseases. Our results suggest that mutations in NDUFA1 do not cause the gender difference observed in clinically severe and complex phenotypes with complex I deficiency.