Search

A large-scale comparative study in two human populations analyzes human plasma transferrin and immunoglobulin G (IgG) N-glycosylation profiles and their association with age, sex, and various biochemical and physiological traits.

Post-translational modifications are known to diversify protein functions, but the effect of genetic variation on the modifications is not well known. Here, the authors find both shared and protein-specific genetic mechanisms regulating the glycosylation of two different proteins.

Abstract Background With the exponential growth in available biomedical data, there is a need for data integration methods that can extract information about relationships between the data sets. However, these data sets might have very different characteristics. For interpretable results, data-specific variation needs to be quantified. For this task, Two-way Orthogonal Partial Least Squares (O2PLS) has been proposed. To facilitate application and development of the methodology, free and open-source software is required. However, this is not the case with O2PLS. Results We introduce OmicsPLS, an open-source implementation of the O2PLS method in R. It can handle both low- and high-dimensional datasets efficiently. Generic methods for inspecting and visualizing results are implemented. Both a standard and faster alternative cross-validation methods are available to determine the number of components. A simulation study shows good performance of OmicsPLS compared to alternatives, in terms of accuracy and CPU runtime. We demonstrate OmicsPLS by integrating genetic and glycomic data. Conclusions We propose the OmicsPLS R package: a free and open-source implementation of O2PLS for statistical data integration. OmicsPLS is available at https://cran.r-project.org/package=OmicsPLS and can be installed in R via install.packages(“OmicsPLS”).

IgG glycosylation is an important factor in immune function, yet the molecular details of protein glycosylation remain poorly understood. The data-driven approach presented here uses large-scale plasma IgG mass spectrometry measurements to infer new biochemical reactions in the glycosylation pathway.

Multivariate analysis methods can uncover the relationship between phenotypic measures characterised by modern omic techniques. Here the authors conduct a multivariate GWAS on IgG N-glycosylation phenotypes and identify 5 novel loci enriched in immune system genes.

Multiple sclerosis (MS) is an inflammatory autoimmune disorder affecting the central nervous system (CNS), with unresolved aetiology. Previous studies have implicated N-glycosylation, a highly regulated enzymatic attachment of complex sugars to targeted proteins, in MS pathogenesis. We investigated individual variation in N-glycosylation of the total plasma proteome and of IgG in MS. Both plasma protein and IgG N-glycans were chromatographically profiled and quantified in 83 MS cases and 88 age- and sex-matched controls. Comparing levels of glycosylation features between MS cases and controls revealed that core fucosylation (p = 6.96 × 10−3) and abundance of high-mannose structures (p = 1.48 × 10−2) were the most prominently altered IgG glycosylation traits. Significant changes in plasma protein N-glycome composition were observed for antennary fucosylated, tri- and tetrasialylated, tri- and tetragalactosylated, high-branched N-glycans (p-value range 1.66 × 10−2–4.28 × 10−2). Classification performance of N-glycans was examined by ROC curve analysis, resulting in an AUC of 0.852 for the total plasma N-glycome and 0.798 for IgG N-glycome prediction models. Our results indicate that multiple aspects of protein glycosylation are altered in MS, showing increased proinflammatory potential. N-glycan alterations showed substantial value in classification of the disease status, nonetheless, additional studies are warranted to explore their exact role in MS development and utility as biomarkers.

Correction to: Nature Communications (2017) 8:1231. doi:10.1038/s41467-017-01525-0

Abstract Biological age captures physiological deterioration better than chronological age and is amenable to interventions. Blood-based biomarkers have been identified as suitable candidates for biological age estimation. This study aims to improve biological age estimation using machine learning models and a feature-set of 60 circulating biomarkers available from the UK Biobank (n = 306,116). We implement an Elastic-Net derived Cox model with 25 selected biomarkers to predict mortality risk (C-Index = 0.778; 95% CI [0.767–0.788]), which outperforms the well-known blood-biomarker based PhenoAge model (C-Index = 0.750; 95% CI [0.739–0.761]), providing a C-Index lift of 0.028 representing an 11% relative increase in predictive value. Importantly, we then show that using common clinical assay panels, with few biomarkers, alongside imputation and the model derived on the full set of biomarkers, does not substantially degrade predictive accuracy from the theoretical maximum achievable for the available biomarkers. Biological age is estimated as the equivalent age within the same-sex population which corresponds to an individual’s mortality risk. Values ranged between 20-years younger and 20-years older than individuals’ chronological age, exposing the magnitude of ageing signals contained in blood markers. Thus, we demonstrate a practical and cost-efficient method of estimating an improved measure of Biological Age, available to the general population.

Our understanding of the genetic contribution to human adiposity is incomplete, as few studies measure adiposity directly. Here, the authors impute whole-body imaging adiposity phenotypes in large biobanks, enhancing their power to discover genes driving human adiposity, and further investigate one such gene using a mouse model.

A large-scale GWAS provides insight on diabetes-dependent genetic effects on the glomerular filtration rate, a common metric to monitor kidney health in disease.

Abstract Background Atherosclerotic cardiovascular diseases (CVD) is the leading cause of death in diabetes, but the full range of biomarkers reflecting atherosclerotic burden and CVD risk in people with diabetes is unknown. Metabolomics may help identify novel biomarkers potentially involved in development of atherosclerosis. We investigated the serum metabolomic profile of subclinical atherosclerosis, measured using ankle brachial index (ABI), in people with type 2 diabetes, compared with the profile for symptomatic CVD in the same population. Methods The Edinburgh Type 2 Diabetes Study is a cohort of 1,066 individuals with type 2 diabetes. ABI was measured at baseline, years 4 and 10, with cardiovascular events assessed at baseline and during 10 years of follow-up. A panel of 228 metabolites was measured at baseline using nuclear magnetic resonance spectrometry, and their association with both ABI and prevalent CVD was explored using univariate regression models and least absolute shrinkage and selection operator (LASSO). Metabolites associated with baseline ABI were further explored for association with follow-up ABI and incident CVD. Results Mean (standard deviation, SD) ABI at baseline was 0.97 (0.18, N = 1025), and prevalence of CVD was 35.0%. During 10-year follow-up, mean (SD) change in ABI was + 0.006 (0.178, n = 436), and 257 CVD events occurred. Lactate, glycerol, creatinine and glycoprotein acetyls levels were associated with baseline ABI in both univariate regression [βs (95% confidence interval, CI) ranged from − 0.025 (− 0.036, − 0.015) to − 0.023 (− 0.034, − 0.013), all p

Objective: Deep sequencing offers unparalleled access to rare variants in human populations. Understanding their role in disease is a priority, yet prohibitive sequencing costs mean that many cohorts lack the sample size to discover these effects on their own. Meta-analysis of individual variant scores allows the combination of rare variants across cohorts and study of their aggregated effect at the gene level, boosting discovery power. However, the methods involved have largely not been field-tested. In this study, we aim to perform the first meta-analysis of gene-based rare variant aggregation optimal tests, applied to the human cardiometabolic proteome. Methods: Here, we carry out this analysis across MANOLIS, Pomak and ORCADES, three isolated European cohorts with whole-genome sequencing (total N = 4,422). We examine the genetic architecture of 250 proteomic traits of cardiometabolic relevance. We use a containerised pipeline to harmonise variant lists across cohorts and define four sets of qualifying variants. For every gene, we interrogate protein-damaging variants, exonic variants, exonic and regulatory variants, and regulatory only variants, using the CADD and Eigen scores to weigh variants according to their predicted functional consequence. We perform single-cohort rare variant analysis and meta-analyse variant scores using the SMMAT package. Results: We describe 5 rare variant pQTLs (RV-pQTL) which pass our stringent significance threshold (7.45 × 10−11) and quality control procedure. These were split between four cis signals for MARCO, TEK, MMP2 and MPO, and one trans association for GDF2 in the SERPINA11 gene. We show that the cis-MPO association, which was not detectable using the single-point data alone, is driven by 5 missense and frameshift variants. These include rs140636390 and rs119468010, which are specific to MANOLIS and ORCADES, respectively. We show how this kind of signal could improve the predictive accuracy of genetic factors in common complex disease such as stroke and cardiovascular disease. Conclusions: Our proof-of-concept study demonstrates the power of gene-based meta-analyses for discovering disease-relevant associations complementing common-variant signals by incorporating population-specific rare variation.

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P

Human population isolates provide a snapshot of the impact of historical demographic processes on population genetics. Such data facilitate studies of the functional impact of rare sequence variants on biomedical phenotypes, as strong genetic drift can result in higher frequencies of variants that are otherwise rare. We present the first whole genome sequencing (WGS) study of the VIKING cohort, a representative collection of samples from the isolated Shetland population in northern Scotland, and explore how its genetic characteristics compare to a mainland Scottish population. Our analyses reveal the strong contributions played by the founder effect and genetic drift in shaping genomic variation in the VIKING cohort. About one tenth of all high-quality variants discovered are unique to the VIKING cohort or are seen at frequencies at least ten fold higher than in more cosmopolitan control populations. Multiple lines of evidence also suggest relaxation of purifying selection during the evolutionary history of the Shetland isolate. We demonstrate enrichment of ultra-rare VIKING variants in exonic regions and for the first time we also show that ultra-rare variants are enriched within regulatory regions, particularly promoters, suggesting that gene expression patterns may diverge relatively rapidly in human isolates.

Immunoglobulin G (IgG), a glycoprotein secreted by plasma B-cells, plays a major role in the human adaptive immune response and are associated with a wide range of diseases. Glycosylation of the Fc binding region of IgGs, responsible for the antibody’s effector function, is essential for prompting a proper immune response. This study focuses on the general genetic impact on IgG glycosylation as well as corresponding subclass specificities. To identify genetic loci involved in IgG glycosylation, we performed a genome-wide association study (GWAS) on liquid chromatography electrospray mass spectrometry (LC–ESI-MS)—measured IgG glycopeptides of 1,823 individuals in the Cooperative Health Research in the Augsburg Region (KORA F4) study cohort. In addition, we performed GWAS on subclass-specific ratios of IgG glycans to gain power in identifying genetic factors underlying single enzymatic steps in the glycosylation pathways. We replicated our findings in 1,836 individuals from the Leiden Longevity Study (LLS). We were able to show subclass-specific genetic influences on single IgG glycan structures. The replicated results indicate that, in addition to genes encoding for glycosyltransferases (i.e., ST6GAL1, B4GALT1, FUT8, and MGAT3), other genetic loci have strong influences on the IgG glycosylation patterns. A novel locus on chromosome 1, harboring RUNX3, which encodes for a transcription factor of the runt domain-containing family, is associated with decreased galactosylation. Interestingly, members of the RUNX family are cross-regulated, and RUNX3 is involved in both IgA class switching and B-cell maturation as well as T-cell differentiation and apoptosis. Besides the involvement of glycosyltransferases in IgG glycosylation, we suggest that, due to the impact of variants within RUNX3, potentially mechanisms involved in B-cell activation and T-cell differentiation during the immune response as well as cell migration and invasion involve IgG glycosylation.