Quantitative trait linkage analysis of longitudinal change in body weight

One of the great strengths of the Framingham Heart Study data, provided for the Genetic Analysis Workshop 13, is the long-term survey of phenotypic data. We used this unique data to create new phenotypes representing the pattern of longitudinal change of the provided phenotypes, especially systolic blood pressure and body weight. We performed a linear regression of body weight and systolic blood pressure on age and took the slopes as new phenotypes for quantitative trait linkage analysis using the SOLAR package. There was no evidence for heritability of systolic blood pressure change. Heritability was estimated as 0.15 for adult life "body weight change", measured as the regression slope, and "body weight gain" (including only individuals with a positive regression slope), and as 0.22 for body weight "change up to 50" (regression slope of weight on age up to an age of 50). With multipoint analysis, two regions on the long arm of chromosome 8 showed the highest LOD scores of 1.6 at 152 cM for "body weight change" and of >1.9 around location 102 cM for "body weight gain" and "change up to 50". The latter two LOD scores almost reach the threshold for suggestive linkage. We conclude that the chromosome 8 region may harbor a gene acting on long-term body weight regulation, thereby contributing to the development of the metabolic syndrome

Linkage analysis of alcohol dependence using MOD scores

Alcohol dependence is a typical example of a complex trait that is governed by several genes and for which the mode of inheritance is unknown. We analyzed the microsatellite markers and the Affymetrix single-nucleotide polymorphisms (SNPs) for a subset of the Collaborative Study on the Genetics of Alcoholism family sample, 93 pedigrees of Caucasian ancestry comprising 919 persons, 390 of whom are affected according to DSM III-R and Feighner criteria. In particular, we performed parametric single-marker linkage analysis using MLINK of the LINKAGE package (for the microsatellite data), as well as multipoint MOD-score analysis with GENEHUNTER-MODSCORE (for the microsatellite and SNP data). By use of two liability classes, different penetrances were assigned to males and females. In order to investigate parent-of-origin effects, we calculated MOD scores under trait models with and without imprinting. In addition, for the microsatellite data, the MOD-score analysis was performed with sex-averaged as well as sex-specific maps. The highest linkage peaks were obtained on chromosomes 1, 2, 7, 10, 12, 13, 15, and 21. There was evidence for paternal imprinting at the loci on chromosomes 2, 10, 12, 13, 15, and 21. A tendency to maternal imprinting was observed at two loci on chromosome 7. Our findings underscore the fact that an adequate modeling of the genotype-phenotype relation is crucial for the genetic mapping of a complex trait

Characterization of missing values in untargeted MS-based metabolomics data and evaluation of missing data handling strategies.

BACKGROUND: Untargeted mass spectrometry (MS)-based metabolomics data often contain missing values that reduce statistical power and can introduce bias in biomedical studies. However, a systematic assessment of the various sources of missing values and strategies to handle these data has received little attention. Missing data can occur systematically, e.g. from run day-dependent effects due to limits of detection (LOD); or it can be random as, for instance, a consequence of sample preparation. METHODS: We investigated patterns of missing data in an MS-based metabolomics experiment of serum samples from the German KORA F4 cohort (n = 1750). We then evaluated 31 imputation methods in a simulation framework and biologically validated the results by applying all imputation approaches to real metabolomics data. We examined the ability of each method to reconstruct biochemical pathways from data-driven correlation networks, and the ability of the method to increase statistical power while preserving the strength of established metabolic quantitative trait loci. RESULTS: Run day-dependent LOD-based missing data accounts for most missing values in the metabolomics dataset. Although multiple imputation by chained equations performed well in many scenarios, it is computationally and statistically challenging. K-nearest neighbors (KNN) imputation on observations with variable pre-selection showed robust performance across all evaluation schemes and is computationally more tractable. CONCLUSION: Missing data in untargeted MS-based metabolomics data occur for various reasons. Based on our results, we recommend that KNN-based imputation is performed on observations with variable pre-selection since it showed robust results in all evaluation schemes.This work was supported by grants from the German Federal Ministry of Education and Research (BMBF), by BMBF Grant No. 01ZX1313C (project e:Athero-MED) and Grant No. 03IS2061B (project Gani_Med). Moreover, the research leading to these results has received funding from the European Union’s Seventh Framework Programme [FP7-Health-F5-2012] under grant agreement No. 305280 (MIMOmics) and from the European Research Council (starting grant “LatentCauses”). KS is supported by Biomedical Research Program funds at Weill Cornell Medical College in Qatar, a program funded by the Qatar Foundation. The KORA Augsburg studies were financed by the Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany and supported by grants from the German Federal Ministry of Education and Research (BMBF). Analyses in the EPIC-Norfolk study were supported by funding from the Medical Research Council (MC_PC_13048 and MC_UU_12015/1)

The permanently chaperone-active small heat shock protein Hsp17 from Caenorhabditis elegans exhibits topological separation of its N-terminal regions

Small Heat shock proteins (sHsps) are a family of molecular chaperones that bind nonnative proteins in an ATP-independent manner. Caenorhabditis elegans encodes 16 different sHsps, among them Hsp17, which is evolutionarily distinct from other sHsps in the nematode. The structure and mechanism of Hsp17 and how these may differ from other sHsps remain unclear. Here, we find that Hsp17 has a distinct expression pattern, structural organization, and chaperone function. Consistent with its presence under nonstress conditions, and in contrast to many other sHsps, we determined that Hsp17 is a mono-disperse, permanently active chaperone in vitro, which interacts with hundreds of different C. elegans proteins under physiological conditions. Additionally, our cryo-EM structure of Hsp17 reveals that in the 24-mer complex, 12 N-terminal regions are involved in its chaperone function. These flexible regions are located on the outside of the spherical oligomer, whereas the other 12 N-terminal regions are engaged in stabilizing interactions in its interior. This allows the same region in Hsp17 to perform different functions depending on the topological context. Taken together, our results reveal structural and functional features that further define the structural basis of permanently active sHsps

Circulating metabolites modulated by diet are associated with depression

Metabolome reflects the interplay of genome and exposome at molecular level and thus can provide deep insights into the pathogenesis of a complex disease like major depression. To identify metabolites associated with depression we performed a metabolome-wide association analysis in 13,596 participants from five European population-based cohorts characterized for depression, and circulating metabolites using ultra high-performance liquid chromatography/tandem accurate mass spectrometry (UHPLC/MS/MS) based Metabolon platform. We tested 806 metabolites covering a wide range of biochemical processes including those involved in lipid, amino-acid, energy, carbohydrate, xenobiotic and vitamin metabolism for their association with depression. In a conservative model adjusting for life style factors and cardiovascular and antidepressant medication use we identified 8 metabolites, including 6 novel, significantly associated with depression. In individuals with depression, increased levels of retinol (vitamin A), 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1) (lecithin) and mannitol/sorbitol and lower levels of hippurate, 4-hydroxycoumarin, 2-aminooctanoate (alpha-aminocaprylic acid), 10-undecenoate (11:1n1) (undecylenic acid), 1-linoleoyl-GPA (18:2) (lysophosphatidic acid; LPA 18:2) are observed. These metabolites are either directly food derived or are products of host and gut microbial metabolism of food-derived products. Our Mendelian randomization analysis suggests that low hippurate levels may be in the causal pathway leading towards depression. Our findings highlight putative actionable targets for depression prevention that are easily modifiable through diet interventions.</p

Cell Specific eQTL Analysis without Sorting Cells

The functional consequences of trait associated SNPs are often investigated using expression quantitative trait locus (eQTL) mapping. While trait-associated variants may operate in a cell-type specific manner, eQTL datasets for such cell-types may not always be available. We performed a genome-environment interaction (GxE) meta-analysis on data from 5,683 samples to infer the cell type specificity of whole blood cis-eQTLs. We demonstrate that this method is able to predict neutrophil and lymphocyte specific cis-eQTLs and replicate these predictions in independent cell-type specific datasets. Finally, we show that SNPs associated with Crohn’s disease preferentially affect gene expression within neutrophils, including the archetypal NOD2 locus

Genome-wide analysis identifies 12 loci influencing human reproductive behavior.

The genetic architecture of human reproductive behavior-age at first birth (AFB) and number of children ever born (NEB)-has a strong relationship with fitness, human development, infertility and risk of neuropsychiatric disorders. However, very few genetic loci have been identified, and the underlying mechanisms of AFB and NEB are poorly understood. We report a large genome-wide association study of both sexes including 251,151 individuals for AFB and 343,072 individuals for NEB. We identified 12 independent loci that are significantly associated with AFB and/or NEB in a SNP-based genome-wide association study and 4 additional loci associated in a gene-based effort. These loci harbor genes that are likely to have a role, either directly or by affecting non-local gene expression, in human reproduction and infertility, thereby increasing understanding of these complex traits

Polygenic prediction of educational attainment within and between families from genome-wide association analyses in 3 million individuals

We conduct a genome-wide association study (GWAS) of educational attainment (EA) in a sample of ~3 million individuals and identify 3,952 approximately uncorrelated genome-wide-significant single-nucleotide polymorphisms (SNPs). A genome-wide polygenic predictor, or polygenic index (PGI), explains 12-16% of EA variance and contributes to risk prediction for ten diseases. Direct effects (i.e., controlling for parental PGIs) explain roughly half the PGI's magnitude of association with EA and other phenotypes. The correlation between mate-pair PGIs is far too large to be consistent with phenotypic assortment alone, implying additional assortment on PGI-associated factors. In an additional GWAS of dominance deviations from the additive model, we identify no genome-wide-significant SNPs, and a separate X-chromosome additive GWAS identifies 57

Search for neutral B meson decays to two charged leptons

The decays $\mathrm{B_d^0,\,B_s^0 \rightarrow e^+e^-,\,\mu^+\mu^-,\, e^\pm\mu^\mp}$ are searched for in 3.5 million hadronic ${\mathrm{Z}}$ events, which constitute the full LEP I data sample collected by the L3 detector. No signals are observed, therefore upper limits at the 90\%(95\%) confidence levels are set on the following branching fractions: % \begin{center}% {\setlength{\tabcolsep}{2pt} \begin{tabular}{lccccclcccc}% % Br$({\mathrm{B_d^0 \rightarrow {\mathrm{e^+e^-}}}})$ & $<$ & $1.4(1.8)$ & $\times$ & $10^{-5}$; & \hspace*{5mm} & Br$({\mathrm{B_s^0 \rightarrow {\mathrm{e^+e^-}}}})$ & $<$ & $5.4(7.0)$ & $\times$ & $10^{-5}$; \\% Br$({\mathrm{B_d^0 \rightarrow \mu^+\mu^-}})$ & $<$ & $1.0(1.4)$ & $\times$ & $10^{-5}$; & \hspace*{5mm} & Br$({\mathrm{B_s^0 \rightarrow \mu^+\mu^-}})$ & $<$ & $3.8(5.1)$ & $\times$ & $10^{-5}$; \\% Br$({\mathrm{B_d^0 \rightarrow {\mathrm{e^\pm\mu^\mp}}}})$ & $<$ & $1.6(2.0)$ & $\times$ & $10^{-5}$; & \hspace*{5mm} & Br$({\mathrm{B_s^0 \rightarrow {\mathrm{e^\pm\mu^\mp}}}})$ & $<$ & $4.1(5.3)$ & $\times$ & $10^{-5}$. \\% % \end{tabular}% } \end{center}% % The results for ${\mathrm{B_s^0\rightarrow{\mathrm{e^+e^-}}}}$ and ${\mathrm{B_s^0 \rightarrow {\mathrm{e^\pm\mu^\mp}}}}$ are the first limits set on these decay modes