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.

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 < 10-9), but they were not replicated when accounting for multiple testing. In total, 27 associations were nominally replicated with a P < 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Uncovering the complex genetics of human character.

Human personality is 30-60% heritable according to twin and adoption studies. Hundreds of genetic variants are expected to influence its complex development, but few have been identified. We used a machine learning method for genome-wide association studies (GWAS) to uncover complex genotypic-phenotypic networks and environmental interactions. The Temperament and Character Inventory (TCI) measured the self-regulatory components of personality critical for health (i.e., the character traits of self-directedness, cooperativeness, and self-transcendence). 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 five clusters of people with distinct profiles of character traits regardless of genotype. Third, we found 42 SNP sets that identified 727 gene loci and were significantly associated with one or more of the character profiles. Each character profile was related to different SNP sets with distinct molecular processes and neuronal functions. Environmental influences measured in childhood and adulthood had small but significant effects. We confirmed the replicability of 95% of the 42 SNP sets in healthy Korean and German samples, as well as their associations with character. The identified SNPs explained nearly all the heritability expected for character in each sample (50 to 58%). We conclude that self-regulatory personality traits are strongly influenced by organized interactions among more than 700 genes despite variable cultures and environments. These gene sets modulate specific molecular processes in brain for intentional goal-setting, self-reflection, empathy, and episodic learning and 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.

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.

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.

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).