Pancreatic Secretory Trypsin Inhibitor (SPINK1) Gene Mutation in Patients with Acute Alcohol Pancreatitis (AAP) Compared to Healthy Controls and Heavy Alcohol Users without Pancreatitis

Only 3-5% of heavy alcohol users develop acute alcohol pancreatitis (AAP). This suggests that additional triggers are required to initiate the inflammatory process. Genetic susceptibility contributes to the development of AAP, and SPINK1 mutation is a documented risk factor. We investigated the prevalence of the SPINK1(N34S) mutation in patients with AAP compared to heavy alcohol users who had never suffered an episode of pancreatitis. Blood samples for the mutational analysis from patients with first episode (n = 60) and recurrent AAP (n = 43) and from heavy alcohol users without a history of AAP (n = 98) as well as from a control population (n = 1914) were obtained. SPINK1 mutation was found in 8.7% of the patients with AAP. The prevalence was significantly lower in healthy controls (3.4%, OR 2.72; 1.32-5.64) and very low in alcoholics without pancreatitis (1.0%, OR 9.29; 1.15-74.74). In a comparison adjusted for potential cofounders between AAP patients and alcoholics, SPINK1 was found to be an independent marker for AAP. The prevalence of the SPINK1 mutation is overrepresented in AAP patients and very low in alcoholics without pancreatitis. This finding may play a role in understanding the variable susceptibility to AAP found in heavy alcohol users.

A genome-wide association study of total child psychiatric problems scores

Substantial genetic correlations have been reported across psychiatric disorders and numerous cross-disorder genetic variants have been detected. To identify the genetic variants underlying general psychopathology in childhood, we performed a genome-wide association study using a total psychiatric problem score. We analyzed 6,844,199 common SNPs in 38,418 school-aged children from 20 population-based cohorts participating in the EAGLE consortium. The SNP heritability of total psychiatric problems was 5.4% (SE = 0.01) and two loci reached genome-wide significance: rs10767094 and rs202005905. We also observed an association of SBF2, a gene associated with neuroticism in previous GWAS, with total psychiatric problems. The genetic effects underlying the total score were shared with common psychiatric disorders only (attention-deficit/hyperactivity disorder, anxiety, depression, insomnia) (rG > 0.49), but not with autism or the less common adult disorders (schizophrenia, bipolar disorder, or eating disorders) (rG < 0.01). Importantly, the total psychiatric problem score also showed at least a moderate genetic correlation with intelligence, educational attainment, wellbeing, smoking, and body fat (rG > 0.29). The results suggest that many common genetic variants are associated with childhood psychiatric symptoms and related phenotypes in general instead of with specific symptoms. Further research is needed to establish causality and pleiotropic mechanisms between related traits.

Afamin predicts the prevalence and incidence of nonalcoholic fatty liver disease

OBJECTIVES

In the general population, increased afamin concentrations are associated with the prevalence and incidence of metabolic syndrome as well as type 2 diabetes. Although metabolic syndrome is commonly associated with nonalcoholic fatty liver disease (NAFLD), there exist no information on afamin and NAFLD.

METHODS

Afamin concentrations were cross-sectionally measured in 146 Austrian patients with NAFLD, in 45 patients without NAFLD, and in 292 age- and sex-matched healthy controls. Furthermore, the feasibility of afamin to predict incident NAFLD was evaluated in 1,434 adult participants in the population-based Cardiovascular Risk in Young Finns Study during a 10-year follow-up.

RESULTS

Median afamin concentrations were significantly higher in NAFLD patients (83.6 mg/L) than in patients without NAFLD (61.6 mg/L, p<0.0001) or in healthy controls (63.9 mg/L, p<0.0001). In age- and sex-adjusted logistic regression analyses a 10 mg/L increase of afamin was associated with a 1.5-fold increase of having NAFLD as compared with patients without NAFLD and the risk was even two-fold when compared with healthy controls. In the population-based cohort, afamin concentrations at baseline were significantly lower in participants without NAFLD (n=1,195) than in 239 participants who developed NAFLD (56.5 vs. 66.9 mg/L, p<0.0001) during the 10-year follow up, with highest afamin values observed in individuals developing severe forms of NAFLD. After adjustment for several potentially confounding parameters, afamin remained an independent predictor for the development of NAFLD (OR=1.37 [95% CI 1.23-1.54] per 10 mg/L increase, p<0.0001).

CONCLUSIONS

Afamin concentrations are increased in patients with NAFLD and independently predict the development of NAFLD in a population-based cohort.

Mitochondrial genome-wide analysis of nuclear DNA methylation quantitative trait loci

Mitochondria have a complex communication network with the surrounding cell and can alter nuclear DNA methylation (DNAm). Variation in the mitochondrial DNA (mtDNA) has also been linked to differential DNAm. Genome-wide association studies have identified numerous DNAm quantitative trait loci, but these studies have not examined the mitochondrial genome. Herein, we quantified nuclear DNAm from blood and conducted a mitochondrial genome-wide association study of DNAm, with an additional emphasis on sex- and prediabetes-specific heterogeneity. We used the Young Finns Study (n = 926) with sequenced mtDNA genotypes as a discovery sample and sought replication in the Ludwigshafen Risk and Cardiovascular Health study (n = 2317). We identified numerous significant associations in the discovery phase (P &lt; 10−9), but they were not replicated when accounting for multiple testing. In total, 27 associations were nominally replicated with a P &lt; 0.05. The replication analysis presented no evidence of sex- or prediabetes-specific heterogeneity. The 27 associations were included in a joint meta-analysis of the two cohorts, and 19 DNAm sites associated with mtDNA variants, while four other sites showed haplogroup associations. An expression quantitative trait methylation analysis was performed for the identified DNAm sites, pinpointing two statistically significant associations. This study provides evidence of a mitochondrial genetic control of nuclear DNAm with little evidence found for sex- and prediabetes-specific effects. The lack of a comparable mtDNA data set for replication is a limitation in our study and further studies are needed to validate our results.

Large-scale cis- and trans-eQTL analyses identify thousands of genetic loci and polygenic scores that regulate blood gene expression

Trait-associated genetic variants affect complex phenotypes primarily via regulatory mechanisms on the transcriptome. To investigate the genetics of gene expression, we performed cis- and trans-expression quantitative trait locus (eQTL) analyses using blood-derived expression from 31,684 individuals through the eQTLGen Consortium. We detected cis-eQTL for 88% of genes, and these were replicable in numerous tissues. Distal trans-eQTL (detected for 37% of 10,317 trait-associated variants tested) showed lower replication rates, partially due to low replication power and confounding by cell type composition. However, replication analyses in single-cell RNA-seq data prioritized intracellular trans-eQTL. Trans-eQTL exerted their effects via several mechanisms, primarily through regulation by transcription factors. Expression of 13% of the genes correlated with polygenic scores for 1,263 phenotypes, pinpointing potential drivers for those traits. In summary, this work represents a large eQTL resource, and its results serve as a starting point for in-depth interpretation of complex phenotypes.

Three genetic-environmental networks for human personality

Phylogenetic, developmental, and brain-imaging studies suggest that human personality is the integrated expression of three major systems of learning and memory that regulate (1) associative conditioning, (2) intentionality, and (3) self-awareness. We have uncovered largely disjoint sets of genes regulating these dissociable learning processes in different clusters of people with (1) unregulated temperament profiles (i.e., associatively conditioned habits and emotional reactivity), (2) organized character profiles (i.e., intentional self-control of emotional conflicts and goals), and (3) creative character profiles (i.e., self-aware appraisal of values and theories), respectively. However, little is known about how these temperament and character components of personality are jointly organized and develop in an integrated manner. In three large independent genome-wide association studies from Finland, Germany, and Korea, we used a data-driven machine learning method to uncover joint phenotypic networks of temperament and character and also the genetic networks with which they are associated. We found three clusters of similar numbers of people with distinct combinations of temperament and character profiles. Their associated genetic and environmental networks were largely disjoint, and differentially related to distinct forms of learning and memory. Of the 972 genes that mapped to the three phenotypic networks, 72% were unique to a single network. The findings in the Finnish discovery sample were blindly and independently replicated in samples of Germans and Koreans. We conclude that temperament and character are integrated within three disjoint networks that regulate healthy longevity and dissociable systems of learning and memory by nearly disjoint sets of genetic and environmental influences.

Genetic association study of childhood aggression across raters, instruments, and age

Childhood aggressive behavior (AGG) has a substantial heritability of around 50%. Here we present a genome-wide association meta-analysis (GWAMA) of childhood AGG, in which all phenotype measures across childhood ages from multiple assessors were included. We analyzed phenotype assessments for a total of 328 935 observations from 87 485 children aged between 1.5 and 18 years, while accounting for sample overlap. We also meta-analyzed within subsets of the data, i.e., within rater, instrument and age. SNP-heritability for the overall meta-analysis (AGGoverall) was 3.31% (SE = 0.0038). We found no genome-wide significant SNPs for AGGoverall. The gene-based analysis returned three significant genes: ST3GAL3 (P = 1.6E-06), PCDH7 (P = 2.0E-06), and IPO13 (P = 2.5E-06). All three genes have previously been associated with educational traits. Polygenic scores based on our GWAMA significantly predicted aggression in a holdout sample of children (variance explained = 0.44%) and in retrospectively assessed childhood aggression (variance explained = 0.20%). Genetic correlations (rg) among rater-specific assessment of AGG ranged from rg = 0.46 between self- and teacher-assessment to rg = 0.81 between mother- and teacher-assessment. We obtained moderate-to-strong rgs with selected phenotypes from multiple domains, but hardly with any of the classical biomarkers thought to be associated with AGG. Significant genetic correlations were observed with most psychiatric and psychological traits (range [Formula: see text]: 0.19-1.00), except for obsessive-compulsive disorder. Aggression had a negative genetic correlation (rg = ~-0.5) with cognitive traits and age at first birth. Aggression was strongly genetically correlated with smoking phenotypes (range [Formula: see text]: 0.46-0.60). The genetic correlations between aggression and psychiatric disorders were weaker for teacher-reported AGG than for mother- and self-reported AGG. The current GWAMA of childhood aggression provides a powerful tool to interrogate the rater-specific genetic etiology of AGG.

Adulthood blood levels of hsa-miR-29b-3p associate with preterm birth and adult metabolic and cognitive health

Preterm birth (PTB) is associated with increased risk of type 2 diabetes and neurocognitive impairment later in life. We analyzed for the first time the associations of PTB with blood miRNA levels in adulthood. We also investigated the relationship of PTB associated miRNAs and adulthood phenotypes previously linked with premature birth. Blood MicroRNA profiling, genome-wide gene expression analysis, computer-based cognitive testing battery (CANTAB) and serum NMR metabolomics were performed for Young Finns Study subjects (aged 34-49 years, full-term n = 682, preterm n = 84). Preterm birth (vs. full-term) was associated with adulthood levels of hsa-miR-29b-3p in a fully adjusted regression model (p = 1.90 × 10^-4, FDR = 0.046). The levels of hsa-miR-29b-3p were down-regulated in subjects with PTB with appropriate birthweight for gestational age (p = 0.002, fold change [FC] = - 1.20) and specifically in PTB subjects with small birthweight for gestational age (p = 0.095, FC = - 1.39) in comparison to individuals born full term. Hsa-miR-29b-3p levels correlated with the expressions of its target-mRNAs BCL11A and CS and the gene set analysis results indicated a target-mRNA driven association between hsa-miR-29b-3p levels and Alzheimer's disease, Parkinson's disease, Insulin signaling and Regulation of Actin Cytoskeleton pathway expression. The level of hsa-miR-29b-3p was directly associated with visual processing and sustained attention in CANTAB test and inversely associated with serum levels of VLDL subclass component and triglyceride levels. In conlcusion, adult blood levels of hsa-miR-29b-3p were lower in subjects born preterm. Hsa-miR-29b-3p associated with cognitive function and may be linked with adulthood morbidities in subjects born preterm, possibly through regulation of gene sets related to neurodegenerative diseases and insulin signaling as well as VLDL and triglyceride metabolism.

Genome-wide association study of circulating interleukin 6 levels identifies novel loci

Interleukin 6 (IL-6) is a multifunctional cytokine with both pro- and anti-inflammatory properties with a heritability estimate of up to 61%. The circulating levels of IL-6 in blood have been associated with an increased risk of complex disease pathogenesis. We conducted a two-staged, discovery and replication meta genome-wide association study (GWAS) of circulating serum IL-6 levels comprising up to 67 428 (ndiscovery = 52 654 and nreplication = 14 774) individuals of European ancestry. The inverse variance fixed effects based discovery meta-analysis, followed by replication led to the identification of two independent loci, IL1F10/IL1RN rs6734238 on chromosome (Chr) 2q14, (Pcombined = 1.8 × 10-11), HLA-DRB1/DRB5 rs660895 on Chr6p21 (Pcombined = 1.5 × 10-10) in the combined meta-analyses of all samples. We also replicated the IL6R rs4537545 locus on Chr1q21 (Pcombined = 1.2 × 10-122). Our study identifies novel loci for circulating IL-6 levels uncovering new immunological and inflammatory pathways that may influence IL-6 pathobiology.

Publisher Correction: A meta-analysis of genome-wide association studies identifies multiple longevity genes

A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-22613-2

The Role of Inflammatory Cytokines as Intermediates in the Pathway from Increased Adiposity to Disease

OBJECTIVE

This study aimed to investigate the role of cytokines as intermediates in the pathway from increased adiposity to disease.

METHODS

BMI and circulating levels of up to 41 cytokines were measured in individuals from three Finnish cohort studies (n = 8,293). Mendelian randomization (MR) was used to assess the impact of BMI on circulating cytokines and the impact of BMI-driven cytokines on risk of obesity-related diseases.

RESULTS

Observationally, BMI was associated with 19 cytokines. For every SD increase in BMI, causal effect estimates were strongest for hepatocyte growth factor, monocyte chemotactic protein-1 (MCP-1), and tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) and were as ratios of geometric means 1.13 (95% CI: 1.08-1.19), 1.08 (95% CI: 1.04-1.14), and 1.13 (95% CI: 1.04-1.21), respectively. TRAIL was associated with a small increase in the odds of coronary artery disease (odds ratio: 1.03; 95% CI: 1.00-1.06). There was inconsistent evidence for a protective role of MCP-1 against inflammatory bowel diseases.

CONCLUSIONS

Observational and MR estimates of the effect of BMI on cytokine levels were generally concordant. There was little evidence for an effect of raised levels of BMI-driven cytokines on disease. These findings illustrate the challenges of MR when applied in the context of molecular mediation.

Examining the effect of mitochondrial DNA variants on blood pressure in two Finnish cohorts

High blood pressure (BP) is a major risk factor for many noncommunicable diseases. The effect of mitochondrial DNA single-nucleotide polymorphisms (mtSNPs) on BP is less known than that of nuclear SNPs. We investigated the mitochondrial genetic determinants of systolic, diastolic, and mean arterial BP. MtSNPs were determined from peripheral blood by sequencing or with genome-wide association study SNP arrays in two independent Finnish cohorts, the Young Finns Study and the Finnish Cardiovascular Study, respectively. In total, over 4200 individuals were included. The effects of individual common mtSNPs, with an additional focus on sex-specificity, and aggregates of rare mtSNPs grouped by mitochondrial genes were evaluated by meta-analysis of linear regression and a sequence kernel association test, respectively. We accounted for the predicted pathogenicity of the rare variants within protein-encoding and the tRNA regions. In the meta-analysis of 87 common mtSNPs, we did not observe significant associations with any of the BP traits. Sex-specific and rare-variant analyses did not pinpoint any significant associations either. Our results are in agreement with several previous studies suggesting that mtDNA variation does not have a significant role in the regulation of BP. Future studies might need to reconsider the mechanisms thought to link mtDNA with hypertension.

DNA methylation and lipid metabolism: an EWAS of 226 metabolic measures

BACKGROUND

The discovery of robust and trans-ethnically replicated DNA methylation markers of metabolic phenotypes, has hinted at a potential role of epigenetic mechanisms in lipid metabolism. However, DNA methylation and the lipid compositions and lipid concentrations of lipoprotein sizes have been scarcely studied. Here, we present an epigenome-wide association study (EWAS) (N = 5414 total) of mostly lipid-related metabolic measures, including a fine profiling of lipoproteins. As lipoproteins are the main players in the different stages of lipid metabolism, examination of epigenetic markers of detailed lipoprotein features might improve the diagnosis, prognosis, and treatment of metabolic disturbances.

RESULTS

We conducted an EWAS of leukocyte DNA methylation and 226 metabolic measurements determined by nuclear magnetic resonance spectroscopy in the population-based KORA F4 study (N = 1662) and replicated the results in the LOLIPOP, NFBC1966, and YFS cohorts (N = 3752). Follow-up analyses in the discovery cohort included investigations into gene transcripts, metabolic-measure ratios for pathway analysis, and disease endpoints. We identified 161 associations (p value < 4.7 × 10-10), covering 16 CpG sites at 11 loci and 57 metabolic measures. Identified metabolic measures were primarily medium and small lipoproteins, and fatty acids. For apolipoprotein B-containing lipoproteins, the associations mainly involved triglyceride composition and concentrations of cholesterol esters, triglycerides, free cholesterol, and phospholipids. All associations for HDL lipoproteins involved triglyceride measures only. Associated metabolic measure ratios, proxies of enzymatic activity, highlight amino acid, glucose, and lipid pathways as being potentially epigenetically implicated. Five CpG sites in four genes were associated with differential expression of transcripts in blood or adipose tissue. CpG sites in ABCG1 and PHGDH showed associations with metabolic measures, gene transcription, and metabolic measure ratios and were additionally linked to obesity or previous myocardial infarction, extending previously reported observations.

CONCLUSION

Our study provides evidence of a link between DNA methylation and the lipid compositions and lipid concentrations of different lipoprotein size subclasses, thus offering in-depth insights into well-known associations of DNA methylation with total serum lipids. The results support detailed profiling of lipid metabolism to improve the molecular understanding of dyslipidemia and related disease mechanisms.

Uncovering the complex genetics of human temperament

Experimental studies of learning suggest that human temperament may depend on the molecular mechanisms for associative conditioning, which are highly conserved in animals. The main genetic pathways for associative conditioning are known in experimental animals, but have not been identified in prior genome-wide association studies (GWAS) of human temperament. We used a data-driven machine learning method for GWAS to uncover the complex genotypic-phenotypic networks and environmental interactions related to human temperament. In a discovery sample of 2149 healthy Finns, we identified sets of single-nucleotide polymorphisms (SNPs) that cluster within particular individuals (i.e., SNP sets) regardless of phenotype. Second, we identified 3 clusters of people with distinct temperament profiles measured by the Temperament and Character Inventory regardless of genotype. Third, we found 51 SNP sets that identified 736 gene loci and were significantly associated with temperament. The identified genes were enriched in pathways activated by associative conditioning in animals, including the ERK, PI3K, and PKC pathways. 74% of the identified genes were unique to a specific temperament profile. Environmental influences measured in childhood and adulthood had small but significant effects. We confirmed the replicability of the 51 Finnish SNP sets in healthy Korean (90%) and German samples (89%), as well as their associations with temperament. The identified SNPs explained nearly all the heritability expected in each sample (37-53%) despite variable cultures and environments. We conclude that human temperament is strongly influenced by more than 700 genes that modulate associative conditioning by molecular processes for synaptic plasticity and long-term memory.

Multivariate Genome-wide Association Analysis of a Cytokine Network Reveals Variants with Widespread Immune, Haematological, and Cardiometabolic Pleiotropy

Cytokines are essential regulatory components of the immune system, and their aberrant levels have been linked to many disease states. Despite increasing evidence that cytokines operate in concert, many of the physiological interactions between cytokines, and the shared genetic architecture that underlies them, remain unknown. Here, we aimed to identify and characterize genetic variants with pleiotropic effects on cytokines. Using three population-based cohorts (n = 9,263), we performed multivariate genome-wide association studies (GWAS) for a correlation network of 11 circulating cytokines, then combined our results in meta-analysis. We identified a total of eight loci significantly associated with the cytokine network, of which two (PDGFRB and ABO) had not been detected previously. In addition, conditional analyses revealed a further four secondary signals at three known cytokine loci. Integration, through the use of Bayesian colocalization analysis, of publicly available GWAS summary statistics with the cytokine network associations revealed shared causal variants between the eight cytokine loci and other traits; in particular, cytokine network variants at the ABO, SERPINE2, and ZFPM2 loci showed pleiotropic effects on the production of immune-related proteins, on metabolic traits such as lipoprotein and lipid levels, on blood-cell-related traits such as platelet count, and on disease traits such as coronary artery disease and type 2 diabetes.

A meta-analysis of genome-wide association studies identifies multiple longevity genes

Human longevity is heritable, but genome-wide association (GWA) studies have had limited success. Here, we perform two meta-analyses of GWA studies of a rigorous longevity phenotype definition including 11,262/3484 cases surviving at or beyond the age corresponding to the 90th/99th survival percentile, respectively, and 25,483 controls whose age at death or at last contact was at or below the age corresponding to the 60th survival percentile. Consistent with previous reports, rs429358 (apolipoprotein E (ApoE) ε4) is associated with lower odds of surviving to the 90th and 99th percentile age, while rs7412 (ApoE ε2) shows the opposite. Moreover, rs7676745, located near GPR78, associates with lower odds of surviving to the 90th percentile age. Gene-level association analysis reveals a role for tissue-specific expression of multiple genes in longevity. Finally, genetic correlation of the longevity GWA results with that of several disease-related phenotypes points to a shared genetic architecture between health and longevity.

The role of oxytocinergic genes in the intergenerational transmission of parent-child relationship qualities

Parenting qualities are known to transmit across generations, but less is known about genetic processes that may modify how strongly parenting quality carries across generations. We examined in prospective data whether oxytocinergic genes of offspring moderate the intergenerational transmission of warm and accepting parent-child relationship qualities. The sample comprised 1167 Finnish parents (G2, 62% female) and their mothers (G1). At the study baseline, G1 mothers (Mage = 38) reported parent-child relationship qualities towards G2 children (age range 3-18). After 28-34 years, G2 offspring reported parent-child relationship qualities towards their own children using the same questionnaire. A cumulative genetic score was computed for G2 by summing up previously identified four alleles associated with non-optimal parenting or social impairments across OXTR (rs1042778, rs2254298, rs53576) and CD38 (rs3796863) genes. Results indicated no interaction effects of G2 cumulative genetic score on the transmission of parent-child relationship qualities. Among single polymorphisms in OXTR, the interaction effects of rs53576 and rs1042778 were found. G1 maternal emotional warmth was associated with higher G2 emotional warmth among G2 participants with the OXTR rs53576 AA/AG genotype, but not among those with the GG genotype. G1 maternal acceptance was associated with higher G2 acceptance among those G2 participants with the OXTR rs1042778 GG/GT genotype, but not among those with the TT genotype. Oxytocinergic genes may influence sensitivity to quality of parent-child relationship, although this needs replication in future studies.

Whole blood microRNA levels associate with glycemic status and correlate with target mRNAs in pathways important to type 2 diabetes

We analyzed the associations between whole blood microRNA profiles and the indices of glucose metabolism and impaired fasting glucose and examined whether the discovered microRNAs correlate with the expression of their mRNA targets. MicroRNA and gene expression profiling were performed for the Young Finns Study participants (n = 871). Glucose, insulin, and glycated hemoglobin (HbA1c) levels were measured, the insulin resistance index (HOMA2-IR) was calculated, and the glycemic status (normoglycemic [n = 534]/impaired fasting glucose [IFG] [n = 252]/type 2 diabetes [T2D] [n = 24]) determined. Levels of hsa-miR-144-5p, -122-5p, -148a-3p, -589-5p, and hsa-let-7a-5p associated with glycemic status. hsa-miR-144-5p and -148a-3p associated with glucose levels, while hsa-miR-144-5p, -122-5p, -184, and -339-3p associated with insulin levels and HOMA2-IR, and hsa-miR-148a-3p, -15b-3p, -93-3p, -146b-5p, -221-3p, -18a-3p, -642a-5p, and -181-2-3p associated with HbA1c levels. The targets of hsa-miR-146b-5p that correlated with its levels were enriched in inflammatory pathways, and the targets of hsa-miR-221-3p were enriched in insulin signaling and T2D pathways. These pathways showed indications of co-regulation by HbA1c-associated miRNAs. There were significant differences in the microRNA profiles associated with glucose, insulin, or HOMA-IR compared to those associated with HbA1c. The HbA1c-associated miRNAs also correlated with the expression of target mRNAs in pathways important to the development of T2D.

Discovery of mitochondrial DNA variants associated with genome-wide blood cell gene expression: a population-based mtDNA sequencing study

The effect of mitochondrial DNA (mtDNA) variation on peripheral blood transcriptomics in health and disease is not fully known. Sex-specific mitochondrially controlled gene expression patterns have been shown in Drosophila melanogaster but in humans, evidence is lacking. Functional variation in mtDNA may also have a role in the development of type 2 diabetes and its precursor state, i.e. prediabetes. We examined the associations between mitochondrial single-nucleotide polymorphisms (mtSNPs) and peripheral blood transcriptomics with a focus on sex- and prediabetes-specific effects. The genome-wide blood cell expression data of 19 637 probes, 199 deep-sequenced mtSNPs and nine haplogroups of 955 individuals from a population-based Young Finns Study cohort were used. Significant associations were identified with linear regression and analysis of covariance. The effects of sex and prediabetes on the associations between gene expression and mtSNPs were studied using random-effect meta-analysis. Our analysis identified 53 significant expression probe-mtSNP associations after Bonferroni correction, involving 7 genes and 31 mtSNPs. Eight probe-mtSNP signals remained independent after conditional analysis. In addition, five genes showed differential expression between haplogroups. The meta-analysis did not show any significant differences in linear model effect sizes between males and females but identified the association between the OASL gene and mtSNP C16294T to show prediabetes-specific effects. This study pinpoints new independent mtSNPs associated with peripheral blood transcriptomics and replicates six previously reported associations, providing further evidence of the mitochondrial genetic control of blood cell gene expression. In addition, we present evidence that prediabetes might lead to perturbations in mitochondrial control.

Genetic and environmental perturbations lead to regulatory decoherence

Correlation among traits is a fundamental feature of biological systems that remains difficult to study. To address this problem, we developed a flexible approach that allows us to identify factors associated with inter-individual variation in correlation. We use data from three human cohorts to study the effects of genetic and environmental variation on correlations among mRNA transcripts and among NMR metabolites. We first show that environmental exposures (infection and disease) lead to a systematic loss of correlation, which we define as 'decoherence'. Using longitudinal data, we show that decoherent metabolites are better predictors of whether someone will develop metabolic syndrome than metabolites commonly used as biomarkers of this disease. Finally, we demonstrate that correlation itself is under genetic control by mapping hundreds of 'correlation quantitative trait loci (QTLs)'. Together, this work furthers our understanding of how and why coordinated biological processes break down, and points to a potential role for decoherence in disease.

Editorial note

This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Gene-environment correlations in parental emotional warmth and intolerance: genome-wide analysis over two generations of the Young Finns Study

BACKGROUND

Genomic analysis of the child might offer new potential to illuminate human parenting. We examined whether offspring (G2) genome-wide genotype variation (SNPs) is associated with their mother's (G1) emotional warmth and intolerance, indicating a gene-environment correlation. If this association is stronger than between G2's genes and their emotional warmth and intolerance toward their own children, then this would indicate the presence of an evocative gene-environment correlation. To further understand how G1 mother's parenting has been evoked by genetically influenced characteristics of the child (G2), we examined whether child (G2) temperament partially accounted for the association between offspring genes and parental responses.

METHODS

Participants were from the Young Finns Study. G1 mothers (N = 2,349; mean age 39 years) self-reported the emotional warmth and intolerance toward G2 in 1980 when the participants were from 3 to 18 years old. G2 participants answered the same parenting scales in 2007/2012 (N = 1,378; mean age = 38 years in 2007; 59% female) when their children were on average 11 years old. Offspring temperament traits were self-reported in 1992 (G2 age range 15-30 years). Estimation of the phenotypic variance explained by the SNPs of G2 was done by genome-wide complex trait analysis with restricted maximum likelihood (GCTA-GREML).

RESULTS

Results showed that the SNPs of a child (G2) explained 22.6% of the phenotypic variance of maternal intolerance (G1; p-value = .039). G2 temperament trait negative emotionality explained only 2.4% points of this association. G2 genes did not explain G1 emotional warmth or G2's own emotional warmth and intolerance. However, further analyses of a combined measure of both G1 parenting scales found genetic effects. Parent or child gender did not moderate the observed associations.

CONCLUSIONS

Presented genome-wide evidence is pointing to the important role a child plays in affecting and shaping his/her family environment, though the underlying mechanisms remain unclear.

The effect of weight on labor market outcomes: An application of genetic instrumental variables

This paper contributes to the literature on the labor market consequences of obesity by using a novel instrument: genetic risk score, which reflects the predisposition to higher body mass index (BMI) across many genetic loci. We estimate instrumental variable models of the effect of BMI on labor market outcomes using Finnish data that have many strengths, for example, BMI that is measured rather than self-reported, and data on earnings and social income transfers that are from administrative tax records and are thus free of the problems associated with nonresponse, reporting error or top coding. The empirical results are sensitive to whether we use a narrower or broader genetic risk score, and to model specification. For example, models using the narrower genetic risk score as an instrument imply that a one-unit increase in BMI is associated with 6.9% lower wages, 1.8% fewer years employed, and a 3 percentage point higher probability of receiving any social income transfers. However, when we use a newer, broader genetic risk score, we cannot reject the null hypothesis of no effect. Future research using genetic risk scores should examine the sensitivity of their results to the risk score used.

Genome Analyses of >200,000 Individuals Identify 58 Loci for Chronic Inflammation and Highlight Pathways that Link Inflammation and Complex Disorders

C-reactive protein (CRP) is a sensitive biomarker of chronic low-grade inflammation and is associated with multiple complex diseases. The genetic determinants of chronic inflammation remain largely unknown, and the causal role of CRP in several clinical outcomes is debated. We performed two genome-wide association studies (GWASs), on HapMap and 1000 Genomes imputed data, of circulating amounts of CRP by using data from 88 studies comprising 204,402 European individuals. Additionally, we performed in silico functional analyses and Mendelian randomization analyses with several clinical outcomes. The GWAS meta-analyses of CRP revealed 58 distinct genetic loci (p < 5 × 10-8). After adjustment for body mass index in the regression analysis, the associations at all except three loci remained. The lead variants at the distinct loci explained up to 7.0% of the variance in circulating amounts of CRP. We identified 66 gene sets that were organized in two substantially correlated clusters, one mainly composed of immune pathways and the other characterized by metabolic pathways in the liver. Mendelian randomization analyses revealed a causal protective effect of CRP on schizophrenia and a risk-increasing effect on bipolar disorder. Our findings provide further insights into the biology of inflammation and could lead to interventions for treating inflammation and its clinical consequences.

Consortium-based genome-wide meta-analysis for childhood dental caries traits

Prior studies suggest dental caries traits in children and adolescents are partially heritable, but there has been no large-scale consortium genome-wide association study (GWAS) to date. We therefore performed GWAS for caries in participants aged 2.5-18.0 years from nine contributing centres. Phenotype definitions were created for the presence or absence of treated or untreated caries, stratified by primary and permanent dentition. All studies tested for association between caries and genotype dosage and the results were combined using fixed-effects meta-analysis. Analysis included up to 19 003 individuals (7530 affected) for primary teeth and 13 353 individuals (5875 affected) for permanent teeth. Evidence for association with caries status was observed at rs1594318-C for primary teeth [intronic within ALLC, odds ratio (OR) 0.85, effect allele frequency (EAF) 0.60, P 4.13e-8] and rs7738851-A (intronic within NEDD9, OR 1.28, EAF 0.85, P 1.63e-8) for permanent teeth. Consortium-wide estimated heritability of caries was low [h2 of 1% (95% CI: 0%: 7%) and 6% (95% CI 0%: 13%) for primary and permanent dentitions, respectively] compared with corresponding within-study estimates [h2 of 28% (95% CI: 9%: 48%) and 17% (95% CI: 2%: 31%)] or previously published estimates. This study was designed to identify common genetic variants with modest effects which are consistent across different populations. We found few single variants associated with caries status under these assumptions. Phenotypic heterogeneity between cohorts and limited statistical power will have contributed; these findings could also reflect complexity not captured by our study design, such as genetic effects which are conditional on environmental exposure.

Oxytocin receptor gene (OXTR) variant rs1042778 moderates the influence of family environment on changes in perceived social support over time

BACKGROUND

Lack of social support is an established risk factor across health outcomes, making it important to examine its family environmental and genetic determinants.

METHODS

In a 27-year follow-up of the Young Finns Study (N = 2341), we examined with a latent growth curve model whether genes involved in the oxytocin signaling pathway-namely, oxytocin receptor gene (OXTR) variants rs1042778, rs2254298, and rs53576-moderate the effect of early-life social experiences on perceived social support across the life span. Mothers reported the emotional warmth and acceptance towards their children at baseline when the participants were from 3 to 18 years old (1980). Perceived family support and support from friends and peripheral sources were assessed in five follow-ups 18 years apart (1989-2007).

RESULTS

Maternal emotional warmth and acceptance predicted the initial level of perceived social support across subscales, while the rate of change in family support was affected by the family environment only if participants carried the T-allele of OXTR rs1042778. This gene-environment interaction was not found for the rate of change in support from friends and peripheral sources and we also did not find associations between latent growth in perceived social support and OXTR variants rs53576 and rs2254298.

LIMITATIONS

Selective attrition in perceived social support, maternal emotional warmth and acceptance, gender, and SES. Family environment was assessed by a non-standardized measure.

CONCLUSIONS

OXTR rs1042778 polymorphism seems to contribute to changes in perceived family support in that way that some individuals (T-allele carriers) 'recover', to some extent, from the effects of early-life social experiences, whereas others (G/G genotype carriers) do not.

Multi-ethnic genome-wide association study for atrial fibrillation

Atrial fibrillation (AF) affects more than 33 million individuals worldwide1 and has a complex heritability2. We conducted the largest meta-analysis of genome-wide association studies (GWAS) for AF to date, consisting of more than half a million individuals, including 65,446 with AF. In total, we identified 97 loci significantly associated with AF, including 67 that were novel in a combined-ancestry analysis, and 3 that were novel in a European-specific analysis. We sought to identify AF-associated genes at the GWAS loci by performing RNA-sequencing and expression quantitative trait locus analyses in 101 left atrial samples, the most relevant tissue for AF. We also performed transcriptome-wide analyses that identified 57 AF-associated genes, 42 of which overlap with GWAS loci. The identified loci implicate genes enriched within cardiac developmental, electrophysiological, contractile and structural pathways. These results extend our understanding of the biological pathways underlying AF and may facilitate the development of therapeutics for AF.

Identification of a Genetic Variation in ERAP1 Aminopeptidase that Prevents Human Cytomegalovirus miR-UL112-5p-Mediated Immunoevasion

Herein, we demonstrate that HCMV miR-UL112-5p targets ERAP1, thereby inhibiting the processing and presentation of the HCMV pp65_495-503 peptide to specific CTLs. In addition, we show that the rs17481334 G variant, naturally occurring in the ERAP1 3' UTR, preserves ERAP1 from miR-UL112-5p-mediated degradation. Specifically, HCMV miR-UL112-5p binds the 3' UTR of ERAP1 A variant, but not the 3' UTR of ERAP1 G variant, and, accordingly, ERAP1 expression is reduced both at RNA and protein levels only in human fibroblasts homozygous for the A variant. Consistently, HCMV-infected GG fibroblasts were more efficient in trimming viral antigens and being lysed by HCMV-peptide-specific CTLs. Notably, a significantly decreased HCMV seropositivity was detected among GG individuals suffering from multiple sclerosis, a disease model in which HCMV is negatively associated with adult-onset disorder. Overall, our results identify a resistance mechanism to HCMV miR-UL112-5p-based immune evasion strategy with potential implications for individual susceptibility to infection and other diseases.

The Significance of Serum Anti-Borrelia Antibodies in the Diagnostic Work-Up of Uveitis

Purpose.

To assess the utility of testing uveitis patients for anti-Borrelia antibodies in an area endemic for Lyme borreliosis.

Methods.

We examined 161 uveitis patients for serum antibodies to Borrelia burgdorferi by Lyme ELISA. Antibodies were determined in patients with uveitis of unknown etiology and non-selectively from patients with an established diagnosis.

Results.

Concentrations of antibodies to B. burgdorferi were elevated in 26 uveitis patients (16.1%), with elevated IgG in 11 of them (6.8%). In four of these patients Lyme borreliosis was a highly suggestive cause of uveitis because of a history of tick bites, systemic symptoms, response to antibiotic therapy, and/or a positive polymerase chain reaction result. Other causes of uveitis were ruled out. All these patients had vitritis.

Conclusions.

Non-selective testing of uveitis patients for Lyme antibodies is not reasonable even in endemic areas. We recommend using the Borrelia antibody test only in cases of uveitis of unknown cause, especially in patients with vitritis or other symptoms of Lyme borreliosis.

Genome-Wide Interactions with Dairy Intake for Body Mass Index in Adults of European Descent

SCOPE

Body weight responds variably to the intake of dairy foods. Genetic variation may contribute to inter-individual variability in associations between body weight and dairy consumption.

METHODS AND RESULTS

A genome-wide interaction study to discover genetic variants that account for variation in BMI in the context of low-fat, high-fat and total dairy intake in cross-sectional analysis was conducted. Data from nine discovery studies (up to 25 513 European descent individuals) were meta-analyzed. Twenty-six genetic variants reached the selected significance threshold (p-interaction <10-7) , and six independent variants (LINC01512-rs7751666, PALM2/AKAP2-rs914359, ACTA2-rs1388, PPP1R12A-rs7961195, LINC00333-rs9635058, AC098847.1-rs1791355) were evaluated meta-analytically for replication of interaction in up to 17 675 individuals. Variant rs9635058 (128 kb 3' of LINC00333) was replicated (p-interaction = 0.004). In the discovery cohorts, rs9635058 interacted with dairy (p-interaction = 7.36 × 10-8) such that each serving of low-fat dairy was associated with 0.225 kg m-2 lower BMI per each additional copy of the effect allele (A). A second genetic variant (ACTA2-rs1388) approached interaction replication significance for low-fat dairy exposure.

CONCLUSION

Body weight responses to dairy intake may be modified by genotype, in that greater dairy intake may protect a genetic subgroup from higher body weight.

Sugar-sweetened beverage intake associations with fasting glucose and insulin concentrations are not modified by selected genetic variants in a ChREBP-FGF21 pathway: a meta-analysis

AIMS/HYPOTHESIS

Sugar-sweetened beverages (SSBs) are a major dietary contributor to fructose intake. A molecular pathway involving the carbohydrate responsive element-binding protein (ChREBP) and the metabolic hormone fibroblast growth factor 21 (FGF21) may influence sugar metabolism and, thereby, contribute to fructose-induced metabolic disease. We hypothesise that common variants in 11 genes involved in fructose metabolism and the ChREBP-FGF21 pathway may interact with SSB intake to exacerbate positive associations between higher SSB intake and glycaemic traits.

METHODS

Data from 11 cohorts (six discovery and five replication) in the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) Consortium provided association and interaction results from 34,748 adults of European descent. SSB intake (soft drinks, fruit punches, lemonades or other fruit drinks) was derived from food-frequency questionnaires and food diaries. In fixed-effects meta-analyses, we quantified: (1) the associations between SSBs and glycaemic traits (fasting glucose and fasting insulin); and (2) the interactions between SSBs and 18 independent SNPs related to the ChREBP-FGF21 pathway.

RESULTS

In our combined meta-analyses of discovery and replication cohorts, after adjustment for age, sex, energy intake, BMI and other dietary covariates, each additional serving of SSB intake was associated with higher fasting glucose (β ± SE 0.014 ± 0.004 [mmol/l], p = 1.5 × 10^-3) and higher fasting insulin (0.030 ± 0.005 [log _e pmol/l], p = 2.0 × 10^-10). No significant interactions on glycaemic traits were observed between SSB intake and selected SNPs. While a suggestive interaction was observed in the discovery cohorts with a SNP (rs1542423) in the β-Klotho (KLB) locus on fasting insulin (0.030 ± 0.011 log _e pmol/l, uncorrected p = 0.006), results in the replication cohorts and combined meta-analyses were non-significant.

CONCLUSIONS/INTERPRETATION

In this large meta-analysis, we observed that SSB intake was associated with higher fasting glucose and insulin. Although a suggestive interaction with a genetic variant in the ChREBP-FGF21 pathway was observed in the discovery cohorts, this observation was not confirmed in the replication analysis.

TRIAL REGISTRATION

Trials related to this study were registered at clinicaltrials.gov as NCT00005131 (Atherosclerosis Risk in Communities), NCT00005133 (Cardiovascular Health Study), NCT00005121 (Framingham Offspring Study), NCT00005487 (Multi-Ethnic Study of Atherosclerosis) and NCT00005152 (Nurses' Health Study).

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.

Extended Serum Lipid Profile Predicting Long-Term Survival in Patients Treated for Abdominal Aortic Aneurysms

BACKGROUND

Individuals treated for abdominal aortic aneurysms (AAAs) are high-risk patients in whom better risk prediction could improve survival. Contemporary serum lipid parameters, such as apolipoproteins and lipoprotein subfractions, may improve or complement the prognostic value of traditional serum lipids. The aim of this study was to ascertain the extended serum lipid profiles, long-term prognosis and their association in AAA patients.

METHODS

Altogether 498 patients treated for AAAs and with available serum lipid values were retrospectively analysed. Contemporary lipid parameters were estimated using a neural network model, the extended Friedewald formula.

RESULTS

Younger age, smoking and urgent or emergency surgery were associated with an unfavourable, and coronary disease and previous stroke with a favourable lipid profile. In multivariable analysis-in addition to advanced age, aneurysm rupture, smoking, pulmonary disease and diabetes-high triglycerides and traditional LDL cholesterol were significant independent risk factors for mortality, HR 1.84 (95% CI 1.20-2.81) and 1.79 (95% CI 1.18-2.73), respectively, while higher EFW-IDL cholesterol was associated with better survival, HR 0.31 (95% CI 0.19-0.65). Including serum lipid parameters improved the prediction of 5-year survival (NRI = 17.7%, p = 0.016).

CONCLUSIONS

Extended serum lipid parameters complement risk prediction of patients treated for AAAs. An unfavourable lipid profile is associated with treatment of AAA earlier in life and with inferior long-term survival.

Genome-wide meta-analysis associates HLA-DQA1/DRB1 and LPA and lifestyle factors with human longevity

Genomic analysis of longevity offers the potential to illuminate the biology of human aging. Here, using genome-wide association meta-analysis of 606,059 parents’ survival, we discover two regions associated with longevity (HLA-DQA1/DRB1 and LPA). We also validate previous suggestions that APOE, CHRNA3/5, CDKN2A/B, SH2B3 and FOXO3A influence longevity. Next we show that giving up smoking, educational attainment, openness to new experience and high-density lipoprotein (HDL) cholesterol levels are most positively genetically correlated with lifespan while susceptibility to coronary artery disease (CAD), cigarettes smoked per day, lung cancer, insulin resistance and body fat are most negatively correlated. We suggest that the effect of education on lifespan is principally mediated through smoking while the effect of obesity appears to act via CAD. Using instrumental variables, we suggest that an increase of one body mass index unit reduces lifespan by 7 months while 1 year of education adds 11 months to expected lifespan.

Variability in human longevity is genetically influenced. Using genetic data of parental lifespan, the authors identify associations at HLA-DQA/DRB1 and LPA and find that genetic variants that increase educational attainment have a positive effect on lifespan whereas increasing BMI negatively affects lifespan.

Blood pathway analyses reveal differences between prediabetic subjects with or without dyslipidaemia. The Cardiovascular Risk in Young Finns Study

BACKGROUND

Prediabetes often occurs together with dyslipidaemia, which is paradoxically treated with statins predisposing to type 2 diabetes mellitus. We examined peripheral blood pathway profiles in prediabetic subjects with (PR_D ) and without dyslipidaemia (PR₀ ) and compared these to nonprediabetic controls without dyslipidaemia (C₀ ).

METHODS

The participants were from the Cardiovascular Risk in Young Finns Study, including 1240 subjects aged 34 to 49 years. Genome-wide expression data of peripheral blood and gene set enrichment analysis were used to investigate the differentially expressed genes and enriched pathways between different subtypes of prediabetes.

RESULTS

Pathways for cholesterol synthesis, interleukin-12-mediated signalling events, and downstream signalling in naïve CD8+ T-cells were upregulated in the PR₀ group in comparison with controls (C₀ ). The upregulation of these pathways was independent of waist circumference, blood pressure, smoking status, and insulin. Adjustment for CRP left the CD8+ T-cell signalling and interleukin-12-mediated signalling event pathway upregulated. The cholesterol synthesis pathway was also upregulated when all prediabetic subjects (PR₀ and PR_D ) were compared with the nonprediabetic control group. No pathways were upregulated or downregulated when the PR_D group was compared with the C₀ group. Five genes in the PR₀ group and 1 in the PR_D group were significantly differentially expressed in comparison with the C₀ group.

CONCLUSIONS

Blood cell gene expression profiles differ significantly between prediabetic subjects with and without dyslipidaemia. Whether this classification may be used in detection of prediabetic individuals at a high risk of cardiovascular complications remains to be examined.

Plasma Concentrations of Afamin Are Associated With Prevalent and Incident Type 2 Diabetes: A Pooled Analysis in More Than 20,000 Individuals

OBJECTIVE

The human vitamin E-binding glycoprotein afamin is primarily expressed in the liver and has been associated with prevalent and incident metabolic syndrome. These data were in line with observations in transgenic mice. We thus investigated whether afamin concentrations are associated with prediabetes, type 2 diabetes, and insulin resistance (IR).

RESEARCH DESIGN AND METHODS

Individual-level baseline (n = 20,136) and follow-up data (n = 14,017) of eight prospective cohort studies were investigated. Study-level data were combined using random-effects meta-analyses. Main outcomes were prevalent and incident type 2 diabetes, prediabetes, and IR. Discrimination and reclassification of participants was analyzed for incident type 2 diabetes.

RESULTS

Mean afamin concentrations between studies ranged from 61 to 73 mg/L. The eight studies included 1,398 prevalent and 585 incident cases of type 2 diabetes. Each increase of afamin by 10 mg/L was associated with prevalent type 2 diabetes (odds ratio [OR] 1.19 [95% CI 1.12-1.26], P = 5.96 × 10^-8). Afamin was positively associated with IR assessed by HOMA-IR (β 0.110 [95% CI 0.089-0.132], P = 1.37 × 10^-23). Most importantly, afamin measured at baseline was an independent predictor for 585 incident cases of type 2 diabetes (OR 1.30 [95% CI 1.23-1.38], P = 3.53 × 10^-19) and showed a significant and valuable gain in risk classification accuracy when added to this extended adjustment model.

CONCLUSIONS

This pooled analysis in >20,000 individuals showed that afamin is strongly associated with IR, prevalence, and incidence of type 2 diabetes independent of major metabolic risk factors or parameters. Afamin might be a promising novel marker for the identification of individuals at high risk for the development of type 2 diabetes.

A genome-wide association meta-analysis on lipoprotein (a) concentrations adjusted for apolipoprotein (a) isoforms

High lipoprotein (a) [Lp(a)] concentrations are an independent risk factor for cardiovascular outcomes. Concentrations are strongly influenced by apo(a) kringle IV repeat isoforms. We aimed to identify genetic loci associated with Lp(a) concentrations using data from five genome-wide association studies (n = 13,781). We identified 48 independent SNPs in the LPA and 1 SNP in the APOE gene region to be significantly associated with Lp(a) concentrations. We also adjusted for apo(a) isoforms to identify loci affecting Lp(a) levels independently from them, which resulted in 31 SNPs (30 in the LPA, 1 in the APOE gene region). Seven SNPs showed a genome-wide significant association with coronary artery disease (CAD) risk. A rare SNP (rs186696265; MAF ∼1%) showed the highest effect on Lp(a) and was also associated with increased risk of CAD (odds ratio = 1.73, P = 3.35 × 10⁻³⁰). Median Lp(a) values increased from 2.1 to 91.1 mg/dl with increasing number of Lp(a)-increasing alleles. We found the APOE2-determining allele of rs7412 to be significantly associated with Lp(a) concentrations (P = 3.47 × 10⁻¹⁰). Each APOE2 allele decreased Lp(a) by 3.34 mg/dl corresponding to ∼15% of the population’s mean values. Performing a gene-based test of association, including suspected Lp(a) receptors and regulators, resulted in one significant association of the TLR2 gene with Lp(a) (P = 3.4 × 10⁻⁴). In summary, we identified a large number of independent SNPs in the LPA gene region, as well as the APOE2 allele, to be significantly associated with Lp(a) concentrations.

Genome-wide Association Study Identifies 27 Loci Influencing Concentrations of Circulating Cytokines and Growth Factors

Circulating cytokines and growth factors are regulators of inflammation and have been implicated in autoimmune and metabolic diseases. In this genome-wide association study (GWAS) of up to 8,293 Finns we identified 27 genome-widely significant loci (p < 1.2 × 10^-9) for one or more cytokines. Fifteen of the associated variants had expression quantitative trait loci in whole blood. We provide genetic instruments to clarify the causal roles of cytokine signaling and upstream inflammation in immune-related and other chronic diseases. We further link inflammatory markers with variants previously associated with autoimmune diseases such as Crohn disease, multiple sclerosis, and ulcerative colitis and hereby elucidate the molecular mechanisms underpinning these diseases and suggest potential drug targets.

Blood hsa-miR-122-5p and hsa-miR-885-5p levels associate with fatty liver and related lipoprotein metabolism-The Young Finns Study

MicroRNAs are involved in disease development and may be utilized as biomarkers. We investigated the association of blood miRNA levels and a) fatty liver (FL), b) lipoprotein and lipid pathways involved in liver lipid accumulation and c) levels of predicted mRNA targets in general population based cohort. Blood microRNA profiling (TaqMan OpenArray), genome-wide gene expression arrays and nuclear magnetic resonance metabolomics were performed for Young Finns Study participants aged 34–49 years (n = 871). Liver fat status was assessed ultrasonographically. Levels of hsa-miR-122-5p and -885-5p were up-regulated in individuals with FL (fold change (FC) = 1.55, p = 1.36 * 10⁻¹⁴ and FC = 1.25, p = 4.86 * 10⁻⁴, respectively). In regression model adjusted with age, sex and BMI, hsa-miR-122-5p and -885-5p predicted FL (OR = 2.07, p = 1.29 * 10⁻⁸ and OR = 1.41, p = 0.002, respectively). Together hsa-miR-122-5p and -885-5p slightly improved the detection of FL beyond established risk factors. These miRNAs may be associated with FL formation through the regulation of lipoprotein metabolism as hsa-miR-122-5p levels associated with small VLDL, IDL, and large LDL lipoprotein subclass components, while hsa-miR-885-5p levels associated inversely with XL HDL cholesterol levels. Hsa-miR-885-5p levels correlated inversely with oxysterol-binding protein 2 (OSBPL2) expression (r = −0.143, p = 1.00 * 10⁻⁴) and suppressing the expression of this lipid receptor and sterol transporter could link hsa-miR-885-5p with HDL cholesterol levels.

USF1 deficiency activates brown adipose tissue and improves cardiometabolic health

USF1 (upstream stimulatory factor 1) is a transcription factor associated with familial combined hyperlipidemia and coronary artery disease in humans. However, whether USF1 is beneficial or detrimental to cardiometabolic health has not been addressed. By inactivating USF1 in mice, we demonstrate protection against diet-induced dyslipidemia, obesity, insulin resistance, hepatic steatosis, and atherosclerosis. The favorable plasma lipid profile, including increased high-density lipoprotein cholesterol and decreased triglycerides, was coupled with increased energy expenditure due to activation of brown adipose tissue (BAT). Usf1 inactivation directs triglycerides from the circulation to BAT for combustion via a lipoprotein lipase-dependent mechanism, thus enhancing plasma triglyceride clearance. Mice lacking Usf1 displayed increased BAT-facilitated, diet-induced thermogenesis with up-regulation of mitochondrial respiratory chain complexes, as well as increased BAT activity even at thermoneutrality and after BAT sympathectomy. A direct effect of USF1 on BAT activation was demonstrated by an amplified adrenergic response in brown adipocytes after Usf1 silencing, and by augmented norepinephrine-induced thermogenesis in mice lacking Usf1. In humans, individuals carrying SNP (single-nucleotide polymorphism) alleles that reduced USF1 mRNA expression also displayed a beneficial cardiometabolic profile, featuring improved insulin sensitivity, a favorable lipid profile, and reduced atherosclerosis. Our findings identify a new molecular link between lipid metabolism and energy expenditure, and point to the potential of USF1 as a therapeutic target for cardiometabolic disease.

Genome-Wide Meta-Analysis of Sciatica in Finnish Population

Sciatica or the sciatic syndrome is a common and often disabling low back disorder in the working-age population. It has a relatively high heritability but poorly understood molecular mechanisms. The Finnish population is a genetic isolate where small founder population and bottleneck events have led to enrichment of certain rare and low frequency variants. We performed here the first genome-wide association (GWAS) and meta-analysis of sciatica. The meta-analysis was conducted across two GWAS covering 291 Finnish sciatica cases and 3671 controls genotyped and imputed at 7.7 million autosomal variants. The most promising loci (p<1x10^-6) were replicated in 776 Finnish sciatica patients and 18,489 controls. We identified five intragenic variants, with relatively low frequencies, at two novel loci associated with sciatica at genome-wide significance. These included chr9:14344410:I (rs71321981) at 9p22.3 (NFIB gene; p = 1.30x10^-8, MAF = 0.08) and four variants at 15q21.2: rs145901849, rs80035109, rs190200374 and rs117458827 (MYO5A; p = 1.34x10^-8, MAF = 0.06; p = 2.32x10^-8, MAF = 0.07; p = 3.85x10^-8, MAF = 0.06; p = 4.78x10^-8, MAF = 0.07, respectively). The most significant association in the meta-analysis, a single base insertion rs71321981 within the regulatory region of the transcription factor NFIB, replicated in an independent Finnish population sample (p = 0.04). Despite identifying 15q21.2 as a promising locus, we were not able to replicate it. It was differentiated; the lead variants within 15q21.2 were more frequent in Finland (6–7%) than in other European populations (1–2%). Imputation accuracies of the three significantly associated variants (chr9:14344410:I, rs190200374, and rs80035109) were validated by genotyping. In summary, our results suggest a novel locus, 9p22.3 (NFIB), which may be involved in susceptibility to sciatica. In addition, another locus, 15q21.2, emerged as a promising one, but failed to replicate.

Talin and vinculin are downregulated in atherosclerotic plaque; Tampere Vascular Study

BACKGROUND AND AIMS

Focal adhesions (FA) play an important role in the tissue remodeling and in the maintenance of tissue integrity and homeostasis. Talin and vinculin proteins are among the major constituents of FAs contributing to cellular well-being and intercellular communication.

METHODS

Microarray analysis (MA) and qRT-PCR low-density array were implemented to analyze talin-1, talin-2, meta-vinculin and vinculin gene expression in circulating blood and arterial plaque.

RESULTS

All analyzed genes were significantly and consistently downregulated in plaques (carotid, abdominal aortic and femoral regions) compared to left internal thoracic artery (LITA) control. The use of LITA samples as controls for arterial plaque samples was validated using immunohistochemistry by comparing LITA samples with healthy arterial samples from a cadaver. Even though the differences in expression levels between stable and unstable plaques were not statistically significant, we observed further negative tendency in the expression in unstable atherosclerotic plaques. The confocal tissue imaging revealed gradient of talin-1 expression in plaque with reduction close to the vessel lumen. Similar gradient was observed for talin-2 expression in LITA controls but was not detected in plaques. This suggests that impaired tissue mechanostability affects the tissue remodeling and healing capabilities leading to development of unstable plaques.

CONCLUSIONS

The central role of talin and vinculin in cell adhesions suggests that the disintegration of the tissue in atherosclerosis could be partially driven by downregulation of these genes, leading to loosening of cell-ECM interactions and remodeling of the tissue.

A genome-wide association meta-analysis on apolipoprotein A-IV concentrations

Apolipoprotein A-IV (apoA-IV) is a major component of HDL and chylomicron particles and is involved in reverse cholesterol transport. It is an early marker of impaired renal function. We aimed to identify genetic loci associated with apoA-IV concentrations and to investigate relationships with known susceptibility loci for kidney function and lipids. A genome-wide association meta-analysis on apoA-IV concentrations was conducted in five population-based cohorts (n = 13,813) followed by two additional replication studies (n = 2,267) including approximately 10 M SNPs. Three independent SNPs from two genomic regions were significantly associated with apoA-IV concentrations: rs1729407 near APOA4 (P = 6.77 × 10 ^{-  44}), rs5104 in APOA4 (P = 1.79 × 10^-24) and rs4241819 in KLKB1 (P = 5.6 × 10^-14). Additionally, a look-up of the replicated SNPs in downloadable GWAS meta-analysis results was performed on kidney function (defined by eGFR), HDL-cholesterol and triglycerides. From these three SNPs mentioned above, only rs1729407 showed an association with HDL-cholesterol (P = 7.1 × 10 ^{-  07}). Moreover, weighted SNP-scores were built involving known susceptibility loci for the aforementioned traits (53, 70 and 38 SNPs, respectively) and were associated with apoA-IV concentrations. This analysis revealed a significant and an inverse association for kidney function with apoA-IV concentrations (P = 5.5 × 10^-05). Furthermore, an increase of triglyceride-increasing alleles was found to decrease apoA-IV concentrations (P = 0.0078). In summary, we identified two independent SNPs located in or next the APOA4 gene and one SNP in KLKB1 The association of KLKB1 with apoA-IV suggests an involvement of apoA-IV in renal metabolism and/or an interaction within HDL particles. Analyses of SNP-scores indicate potential causal effects of kidney function and by lesser extent triglycerides on apoA-IV concentrations.

A genome-wide approach to children's aggressive behavior: The EAGLE consortium

Individual differences in aggressive behavior emerge in early childhood and predict persisting behavioral problems and disorders. Studies of antisocial and severe aggression in adulthood indicate substantial underlying biology. However, little attention has been given to genome-wide approaches of aggressive behavior in children. We analyzed data from nine population-based studies and assessed aggressive behavior using well-validated parent-reported questionnaires. This is the largest sample exploring children's aggressive behavior to date (N = 18,988), with measures in two developmental stages (N = 15,668 early childhood and N = 16,311 middle childhood/early adolescence). First, we estimated the additive genetic variance of children's aggressive behavior based on genome-wide SNP information, using genome-wide complex trait analysis (GCTA). Second, genetic associations within each study were assessed using a quasi-Poisson regression approach, capturing the highly right-skewed distribution of aggressive behavior. Third, we performed meta-analyses of genome-wide associations for both the total age-mixed sample and the two developmental stages. Finally, we performed a gene-based test using the summary statistics of the total sample. GCTA quantified variance tagged by common SNPs (10-54%). The meta-analysis of the total sample identified one region in chromosome 2 (2p12) at near genome-wide significance (top SNP rs11126630, P = 5.30 × 10(-8) ). The separate meta-analyses of the two developmental stages revealed suggestive evidence of association at the same locus. The gene-based analysis indicated association of variation within AVPR1A with aggressive behavior. We conclude that common variants at 2p12 show suggestive evidence for association with childhood aggression. Replication of these initial findings is needed, and further studies should clarify its biological meaning. © 2015 Wiley Periodicals, Inc.

Associations of functional alanine-glyoxylate aminotransferase 2 gene variants with atrial fibrillation and ischemic stroke

Asymmetric and symmetric dimethylarginines (ADMA and SDMA) impair nitric oxide bioavailability and have been implicated in the pathogenesis of atrial fibrillation (AF). Alanine–glyoxylate aminotransferase 2 (AGXT2) is the only enzyme capable of metabolizing both of the dimethylarginines. We hypothesized that two functional AGXT2 missense variants (rs37369, V140I; rs16899974, V498L) are associated with AF and its cardioembolic complications. Association analyses were conducted using 1,834 individulas with AF and 7,159 unaffected individuals from two coronary angiography cohorts and a cohort comprising patients undergoing clinical exercise testing. In coronary angiography patients without structural heart disease, the minor A allele of rs16899974 was associated with any AF (OR = 2.07, 95% CI 1.59-2.68), and with paroxysmal AF (OR = 1.98, 95% CI 1.44–2.74) and chronic AF (OR = 2.03, 95% CI 1.35–3.06) separately. We could not replicate the association with AF in the other two cohorts. However, the A allele of rs16899974 was nominally associated with ischemic stroke risk in the meta-analysis of WTCCC2 ischemic stroke cohorts (3,548 cases, 5,972 controls) and with earlier onset of first-ever ischemic stroke (360 cases) in the cohort of clinical exercise test patients. In conclusion, AGXT2 variations may be involved in the pathogenesis of AF and its age-related thromboembolic complications.

Meta-analysis of Genome-Wide Association Studies for Extraversion: Findings from the Genetics of Personality Consortium

Extraversion is a relatively stable and heritable personality trait associated with numerous psychosocial, lifestyle and health outcomes. Despite its substantial heritability, no genetic variants have been detected in previous genome-wide association (GWA) studies, which may be due to relatively small sample sizes of those studies. Here, we report on a large meta-analysis of GWA studies for extraversion in 63,030 subjects in 29 cohorts. Extraversion item data from multiple personality inventories were harmonized across inventories and cohorts. No genome-wide significant associations were found at the single nucleotide polymorphism (SNP) level but there was one significant hit at the gene level for a long non-coding RNA site (LOC101928162). Genome-wide complex trait analysis in two large cohorts showed that the additive variance explained by common SNPs was not significantly different from zero, but polygenic risk scores, weighted using linkage information, significantly predicted extraversion scores in an independent cohort. These results show that extraversion is a highly polygenic personality trait, with an architecture possibly different from other complex human traits, including other personality traits. Future studies are required to further determine which genetic variants, by what modes of gene action, constitute the heritable nature of extraversion.

Integrative approaches for large-scale transcriptome-wide association studies

Many genetic variants influence complex traits by modulating gene expression, thus altering the abundance of one or multiple proteins. Here we introduce a powerful strategy that integrates gene expression measurements with summary association statistics from large-scale genome-wide association studies (GWAS) to identify genes whose cis-regulated expression is associated with complex traits. We leverage expression imputation from genetic data to perform a transcriptome-wide association study (TWAS) to identify significant expression-trait associations. We applied our approaches to expression data from blood and adipose tissue measured in ∼ 3,000 individuals overall. We imputed gene expression into GWAS data from over 900,000 phenotype measurements to identify 69 new genes significantly associated with obesity-related traits (BMI, lipids and height). Many of these genes are associated with relevant phenotypes in the Hybrid Mouse Diversity Panel. Our results showcase the power of integrating genotype, gene expression and phenotype to gain insights into the genetic basis of complex traits.