Rare Variants in PLXNA4 and Parkinson’s Disease

Approximately 20% of individuals with Parkinson’s disease (PD) report a positive family history. Yet, a large portion of causal and disease-modifying variants is still unknown. We used exome sequencing in two affected individuals from a family with late-onset familial PD followed by frequency assessment in 975 PD cases and 1014 ethnically-matched controls and linkage analysis to identify potentially causal variants. Based on the predicted penetrance and the frequencies, a variant in PLXNA4 proved to be the best candidate and PLXNA4 was screened for additional variants in 862 PD cases and 940 controls, revealing an excess of rare non-synonymous coding variants in PLXNA4 in individuals with PD. Although we cannot conclude that the variant in PLXNA4 is indeed the causative variant, these findings are interesting in the light of a surfacing role of axonal guidance mechanisms in neurodegenerative disorders but, at the same time, highlight the difficulties encountered in the study of rare variants identified by next-generation sequencing in diseases with autosomal dominant or complex patterns of inheritance.     

Nonadditive Effects of Genes in Human Metabolomics

Genome-wide association studies (GWAS) are widely applied to analyze the genetic effects on phenotypes. With the availability of high-throughput technologies for metabolite measurements, GWAS successfully identified loci that affect metabolite concentrations and underlying pathways. In most GWAS, the effect of each SNP on the phenotype is assumed to be additive. Other genetic models such as recessive, dominant, or overdominant were considered only by very few studies. In contrast to this, there are theories that emphasize the relevance of nonadditive effects as a consequence of physiologic mechanisms. This might be especially important for metabolites because these intermediate phenotypes are closer to the underlying pathways than other traits or diseases. In this study we analyzed systematically nonadditive effects on a large panel of serum metabolites and all possible ratios (22,801 total) in a population-based study [Cooperative Health Research in the Region of Augsburg (KORA) F4, N = 1,785]. We applied four different 1-degree-of-freedom (1-df) tests corresponding to an additive, dominant, recessive, and overdominant trait model as well as a genotypic model with two degree-of-freedom (2-df) that allows a more general consideration of genetic effects. Twenty-three loci were found to be genome-wide significantly associated (Bonferroni corrected P ≤ 2.19 × 10⁻¹²) with at least one metabolite or ratio. For five of them, we show the evidence of nonadditive effects. We replicated 17 loci, including 3 loci with nonadditive effects, in an independent study (TwinsUK, N = 846). In conclusion, we found that most genetic effects on metabolite concentrations and ratios were indeed additive, which verifies the practice of using the additive model for analyzing SNP effects on metabolites.     

Multiple loci are associated with white blood cell phenotypes

White blood cell (WBC) count is a common clinical measure from complete blood count assays, and it varies widely among healthy individuals. Total WBC count and its constituent subtypes have been shown to be moderately heritable, with the heritability estimates varying across cell types. We studied 19,509 subjects from seven cohorts in a discovery analysis, and 11,823 subjects from ten cohorts for replication analyses, to determine genetic factors influencing variability within the normal hematological range for total WBC count and five WBC subtype measures. Cohort specific data was supplied by the CHARGE, HeamGen, and INGI consortia, as well as independent collaborative studies. We identified and replicated ten associations with total WBC count and five WBC subtypes at seven different genomic loci (total WBC count-6p21 in the HLA region, 17q21 near ORMDL3, and CSF3; neutrophil count-17q21; basophil count- 3p21 near RPN1 and C3orf27; lymphocyte count-6p21, 19p13 at EPS15L1; monocyte count-2q31 at ITGA4, 3q21, 8q24 an intergenic region, 9q31 near EDG2), including three previously reported associations and seven novel associations. To investigate functional relationships among variants contributing to variability in the six WBC traits, we utilized gene expression- and pathways-based analyses. We implemented gene-clustering algorithms to evaluate functional connectivity among implicated loci and showed functional relationships across cell types. Gene expression data from whole blood was utilized to show that significant biological consequences can be extracted from our genome-wide analyses, with effect estimates for significant loci from the meta-analyses being highly corellated with the proximal gene expression. In addition, collaborative efforts between the groups contributing to this study and related studies conducted by the COGENT and RIKEN groups allowed for the examination of effect homogeneity for genome-wide significant associations across populations of diverse ancestral backgrounds.     

Automatic Scoring and Quality Assessment Using Accuracy Bounds for FP-TDI SNP Genotyping Data

BACKGROUND: Human diversity, namely single nucleotide polymorphisms (SNPs), is becoming a focus of biomedical research. Despite the binary nature of SNP determination, the majority of genotyping assay data need a critical evaluation for genotype calling. We applied statistical models to improve the automated analysis of 2-dimensional SNP data. METHODS: We derived several quantities in the framework of Gaussian mixture models that provide figures of merit to objectively measure the data quality. The accuracy of individual observations is scored as the probability of belonging to a certain genotype cluster, while the assay quality is measured by the overlap between the genotype clusters. RESULTS: The approach was extensively tested with a dataset of 438 nonredundant SNP assays comprising >150,000 datapoints. The performance of our automatic scoring method was compared with manual assignments. The agreement for the overall assay quality is remarkably good, and individual observations were scored differently by man and machine in 2.6% of cases, when applying stringent probability threshold values. CONCLUSION: Our definition of bounds for the accuracy for complete assays in terms of misclassification probabilities goes beyond other proposed analysis methods. We expect the scoring method to minimise human intervention and provide a more objective error estimate in genotype calling.     

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.

     

Epigenetic Signatures at AQP3 and SOCS3 Engage in Low-Grade Inflammation across Different Tissues

BackgroundElevated levels of C-reactive protein (CRP, determined by a high-sensitivity assay) indicate low-grade inflammation which is implicated in many age-related disorders. Epigenetic studies on CRP might discover molecular mechanisms underlying CRP regulation. We aimed to identify DNA methylation sites related to CRP concentrations in cells and tissues regulating low-grade inflammation.ResultsGenome-wide DNA methylation was measured in peripheral blood in 1,741 participants of the KORA F4 study using Illumina HumanMethylation450 BeadChip arrays. Four CpG sites (located at BCL3, AQP3, SOCS3, and cg19821297 intergenic at chromosome 19p13.2, P ≤ 1.01E-07) were significantly hypomethylated at high CRP concentrations independent of various confounders including age, sex, BMI, smoking, and white blood cell composition. Findings were not sex-specific. CRP-related top genes were enriched in JAK/STAT pathways (Benjamini-Hochberg corrected P < 0.05). Results were followed-up in three studies using DNA from peripheral blood (EPICOR, n = 503) and adipose tissue (TwinsUK, n = 368) measured as described above and from liver tissue (LMU liver cohort, n = 286) measured by MALDI-TOF mass spectrometry using EpiTYPER. CpG sites at the AQP3 locus (significant p-values in peripheral blood = 1.72E-03 and liver tissue = 1.51E-03) and the SOCS3 locus (p-values in liver < 2.82E-05) were associated with CRP in the validation panels.ConclusionsEpigenetic modifications seem to engage in low-grade inflammation, possibly via JAK/STAT mediated pathways. Results suggest a shared relevance across different tissues at the AQP3 locus and highlight a role of DNA methylation for CRP regulation at the SOCS3 locus.