Joint modeling of a number of phenotypes using multivariate methods has often been neglected in genome-wide association studies and if used, replication has not been sought. Modern omics technologies ...allow characterization of functional phenomena using a large number of related phenotype measures, which can benefit from such joint analysis. Here, we report a multivariate genome-wide association studies of 23 immunoglobulin G (IgG) N-glycosylation phenotypes. In the discovery cohort, our multi-phenotype method uncovers ten genome-wide significant loci, of which five are novel (IGH, ELL2, HLA-B-C, AZI1, FUT6-FUT3). We convincingly replicate all novel loci via multivariate tests. We show that IgG N-glycosylation loci are strongly enriched for genes expressed in the immune system, in particular antibody-producing cells and B lymphocytes. We empirically demonstrate the efficacy of multivariate methods to discover novel, reproducible pleiotropic effects.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.
The biological and clinical relevance of glycosylation is becoming increasingly recognized, leading to a growing interest in large-scale clinical and population-based studies. In the past few years, ...several methods for high-throughput analysis of glycans have been developed, but thorough validation and standardization of these methods is required before significant resources are invested in large-scale studies. In this study, we compared liquid chromatography, capillary gel electrophoresis, and two MS methods for quantitative profiling of N-glycosylation of IgG in the same data set of 1201 individuals. To evaluate the accuracy of the four methods we then performed analysis of association with genetic polymorphisms and age. Chromatographic methods with either fluorescent or MS-detection yielded slightly stronger associations than MS-only and multiplexed capillary gel electrophoresis, but at the expense of lower levels of throughput. Advantages and disadvantages of each method were identified, which should inform the selection of the most appropriate method in future studies.
Poorer glycemic control in type 1 diabetes may alter N-glycosylation patterns on circulating glycoproteins, and these alterations may be linked with diabetic kidney disease (DKD). We investigated ...associations between N-glycans and glycemic control and renal function in type 1 diabetes.
Using serum samples from 818 adults who were considered to have extreme annual loss in estimated glomerular filtration rate (eGFR; i.e., slope) based on retrospective clinical records, from among 6,127 adults in the Scottish Diabetes Research Network Type 1 Bioresource Study, we measured total and IgG-specific N-glycan profiles. This yielded a relative abundance of 39 total (GP) and 24 IgG (IGP) N-glycans. Linear regression models were used to investigate associations between N-glycan structures and HbA
, albumin-to-creatinine ratio (ACR), and eGFR slope. Models were adjusted for age, sex, duration of type 1 diabetes, and total serum IgG.
Higher HbA
was associated with a lower relative abundance of simple biantennary N-glycans and a higher relative abundance of more complex structures with more branching, galactosylation, and sialylation (GP12, 26, 31, 32, and 34, and IGP19 and 23; all
< 3.79 × 10
). Similar patterns were seen for ACR and greater mean annual loss of eGFR, which were also associated with fewer of the simpler N-glycans (all
< 3.79 × 10
).
Higher HbA
in type 1 diabetes is associated with changes in the serum N-glycome that have elsewhere been shown to regulate the epidermal growth factor receptor and transforming growth factor-β pathways that are implicated in DKD. Furthermore, N-glycans are associated with ACR and eGFR slope. These data suggest that the role of altered N-glycans in DKD warrants further investigation.
Obesity remains an unmet global health burden. Detrimental anatomical distribution of body fat is a major driver of obesity-mediated mortality risk and is demonstrably heritable. However, our ...understanding of the full genetic contribution to human adiposity is incomplete, as few studies measure adiposity directly. To address this, we impute whole-body imaging adiposity phenotypes in UK Biobank from the 4,366 directly measured participants onto the rest of the cohort, greatly increasing our discovery power. Using these imputed phenotypes in 392,535 participants yielded hundreds of genome-wide significant associations, six of which replicate in independent cohorts. The leading causal gene candidate, ADAMTS14, is further investigated in a mouse knockout model. Concordant with the human association data, the Adamts14
mice exhibit reduced adiposity and weight-gain under obesogenic conditions, alongside an improved metabolic rate and health. Thus, we show that phenotypic imputation at scale offers deeper biological insights into the genetics of human adiposity that could lead to therapeutic targets.
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.
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.
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 < 0.0002 and LASSO analysis. The associations remained nominally significant after adjustment for major vascular risk factors. In prospective analyses, lactate was nominally associated with ABI measured at years 4 and 10 after adjustment for baseline ABI. The four ABI-associated metabolites were all positively associated with prevalent CVD odds ratios (ORs) ranged from 1.29 (1.13, 1.47) to 1.49 (1.29, 1.74), all p < 0.0002, and they were also positively associated with incident CVD ORs (95% CI) ranged from 1.19 (1.02, 1.39) to 1.35 (1.17, 1.56), all p < 0.05.
Serum metabolites relating to glycolysis, fluid balance and inflammation were independently associated with both a marker of subclinical atherosclerosis and with symptomatic CVD in people with type 2 diabetes. Additional investigation is warranted to determine their roles as possible etiological and/or predictive biomarkers for atherosclerotic CVD.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
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.,
, and
), other genetic loci have strong influences on the IgG glycosylation patterns. A novel locus on chromosome 1, harboring
, which encodes for a transcription factor of the runt domain-containing family, is associated with decreased galactosylation. Interestingly, members of the
family are cross-regulated, and
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
, 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.
Immunoglobulin G (IgG) is a major effector molecule of the human immune response, and aberrations in IgG glycosylation are linked to various diseases. However, the molecular mechanisms underlying ...protein glycosylation are still poorly understood. We present a data-driven approach to infer reactions in the IgG glycosylation pathway using large-scale mass-spectrometry measurements. Gaussian graphical models are used to construct association networks from four cohorts. We find that glycan pairs with high partial correlations represent enzymatic reactions in the known glycosylation pathway, and then predict new biochemical reactions using a rule-based approach. Validation is performed using data from a GWAS and results from three in vitro experiments. We show that one predicted reaction is enzymatically feasible and that one rejected reaction does not occur in vitro. Moreover, in contrast to previous knowledge, enzymes involved in our predictions colocalize in the Golgi of two cell lines, further confirming the in silico predictions.
Effector functions of immunoglobulin G (IgG) are regulated by the composition of a glycan moiety, thus affecting activity of the immune system. Aberrant glycosylation of IgG has been observed in many ...diseases, but little is understood about the underlying mechanisms. We performed a genome-wide association study of IgG N-glycosylation (
= 8090) and, using a data-driven network approach, suggested how associated loci form a functional network. We confirmed in vitro that knockdown of
decreases the expression of fucosyltransferase FUT8, resulting in increased levels of fucosylated glycans, and suggest that RUNX1 and RUNX3, together with SMARCB1, regulate expression of glycosyltransferase MGAT3. We also show that variants affecting the expression of genes involved in the regulation of glycoenzymes colocalize with variants affecting risk for inflammatory diseases. This study provides new evidence that variation in key transcription factors coupled with regulatory variation in glycogenes modifies IgG glycosylation and has influence on inflammatory diseases.
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.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
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.
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.
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.
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.
•Meta-analysis of gene-based variant scores leverages population-specific allelic heterogeneity.•We study rare variants influencing 250 circulating serum proteins in two Greek population isolates.•4 cis and 1 trans coding variant signals are study-wide significant and pass quality control.•Five rare coding variants influence protein expression of myeloperoxidase (MPO).•Different isolates are enriched for different rare coding MPO variants.