Bidirectional Mendelian Randomization Study of Personality Traits Reveals a Positive Feedback Loop Between Neuroticism and Back Pain.

We conducted a bidirectional Mendelian randomization study to examine the causal effects of six personality traits (anxiety, neuroticism, extraversion, openness to experience, agreeableness, and conscientiousness) on back pain associated with health care use and the causal effect of back pain on the same risk factors. Genetic instruments for the personality traits and back pain were obtained from the largest published genome-wide association studies conducted in individuals of European ancestry. We used inverse weighted variance meta-analysis and Causal Analysis Using Summary Effect for primary analyses and sensitivity analyses to examine evidence for causal associations. We interpreted exposure-outcome associations as being consistent with a causal relationship if results of at least one primary analysis were statistically significant after accounting for multiple statistical testing (P-value < .0042), and the direction and magnitude of effect estimates were concordant between primary and sensitivity analyses. We found evidence for statistically significant bidirectional causal associations between neuroticism and back pain, with odds ratio 1.51 (95% confidence interval 1.37; 1.67) of back pain per neuroticism sum score standard deviation, P-value = 7.80e-16; and beta = .12, se = .04 of neuroticism sum score standard deviation per log odds of back pain, P-value = 2.48e-03. Other relationships did not meet our predefined criteria for causal association. PERSPECTIVE: The significant positive feedback loop between neuroticism and back pain highlights the importance of considering neuroticism in the management of patients with back pain.

Evidence of causal effects of blood pressure on back pain and back pain on type II diabetes provided by a bidirectional Mendelian randomization study.

BACKGROUND CONTEXT: Cardiovascular risk factors (hypertension, dyslipidemia, and type II diabetes) have been proposed as risk factors for back pain. However, few longitudinal studies have found significant associations between cardiovascular risk factors and back pain, and these may be explained by confounding or reverse causation. PURPOSE: To examine potential causal effects of cardiovascular risk factors on back pain, and vice versa. STUDY DESIGN: Bidirectional Mendelian randomization (MR) study. PATIENT SAMPLES: Genome-wide association studies (GWAS) with sample sizes between 173,082 and 1,028,947 participants. OUTCOME MEASURES: Outcomes included (1) back pain associated with health care use (BP-HC) in the forward MR; and (2) seven cardiovascular phenotypes in the reverse MR, including 2 measurements used for the evaluation of hypertension (diastolic blood pressure and systolic blood pressure), 4 phenotypes related to dyslipidemia (LDL cholesterol, HDL cholesterol, total cholesterol, and triglycerides), and type II diabetes. METHODS: We used summary statistics from large, publicly available GWAS for BP-HC and the 7 cardiovascular phenotypes to obtain genetic instrumental variables. We examined MR evidence for causal associations using inverse-variance weighted (IVW) analysis, Causal Analysis Using Summary Effect (CAUSE), and sensitivity analyses. RESULTS: In forward MR analyses of seven cardiovascular phenotypes, diastolic blood pressure was associated with BP-HC across all analyses (IVW estimate: OR = 1.10 per 10.5 mm Hg increase [1.04-1.17], p-value = .001), and significant associations of systolic blood pressure with BP-HC were also found (IVW estimate: OR = 1.09 per 19.3 mm Hg increase [1.04-1.15], p-value = .0006). In reverse MR analyses, only type II diabetes was associated with BP-HC across all analyses (IVW estimate: OR = 1.40 [1.13-1.73], p-value = .002). CONCLUSIONS: These findings from analyses of large, population-based samples indicate that higher blood pressure increases the risk of BP-HC, and BP-HC itself increases the risk of type II diabetes.

Causal effects of psychosocial factors on chronic back pain: a bidirectional Mendelian randomisation study.

PURPOSE: Risk factors for chronic back pain (CBP) may share underlying genetic factors, making them difficult to study using conventional methods. We conducted a bi-directional Mendelian randomisation (MR) study to examine the causal effects of risk factors (education, smoking, alcohol consumption, physical activity, sleep and depression) on CBP and the causal effect of CBP on the same risk factors. METHODS: Genetic instruments for risk factors and CBP were obtained from the largest published genome-wide association studies (GWAS) of risk factor traits conducted in individuals of European ancestry. We used inverse weighted variance meta-analysis (IVW), Causal Analysis Using Summary Effect (CAUSE) and sensitivity analyses to examine evidence for causal associations. We interpreted exposure-outcome associations as being consistent with a causal relationship if results with IVW or CAUSE were statistically significant after accounting for multiple statistical testing (p < 0.003), and the direction and magnitude of effect estimates were concordant between IVW, CAUSE, and sensitivity analyses. RESULTS: We found evidence for statistically significant causal associations between greater education (OR per 4.2 years of schooling = 0.54), ever smoking (OR = 1.27), greater alcohol consumption (OR = 1.29 per consumption category increase) and major depressive disorder (OR = 1.41) and risk of CBP. Conversely, we found evidence for significant causal associations between CBP and greater alcohol consumption (OR = 1.19) and between CBP and smoking (OR = 1.21). Other relationships did not meet our pre-defined criteria for causal association. CONCLUSION: Fewer years of schooling, smoking, greater alcohol consumption, and major depressive disorder increase the risk of CBP. CBP increases the risk of greater alcohol consumption and smoking.

Investigation of the causal relationships between human IgG N-glycosylation and 12 common diseases associated with changes in the IgG N-glycome.

Changes in the N-glycosylation of immunoglobulin G (IgG) are often observed in pathological states, such as autoimmune, inflammatory, neurodegenerative, cardiovascular diseases and some types of cancer. However, in most cases, it is not clear if the disease onset causes these changes, or if the changes in IgG N-glycosylation are among the risk factors for the diseases. The aim of this study was to investigate the casual relationships between IgG N-glycosylation traits and 12 diseases, in which the alterations of IgG N-glycome were previously reported, using two sample Mendelian randomization (MR) approach. We have performed two sample MR using publicly available summary statistics of genome-wide association studies of IgG N-glycosylation and disease risks. Our results indicate positive causal effect of systemic lupus erythematosus (SLE) on the abundance of N-glycans with bisecting N-acetylglucosamine in the total IgG N-glycome. Therefore, we suggest regarding this IgG glycosylation trait as a biomarker of SLE. We also emphasize the need for more powerful GWAS studies of IgG N-glycosylation to further elucidate the causal effect of IgG N-glycome on the diseases.

Genetic Regulation of Immunoglobulin G Glycosylation.

Defining the genetic components that control glycosylation of the human immunoglobulin G (IgG) is an ongoing effort, which has so far been addressed by means of heritability, linkage and genome-wide association studies (GWAS). Unlike the synthesis of proteins, N-glycosylation biosynthesis is not a template-driven process, but rather a complex process regulated by both genetic and environmental factors. Current heritability studies have shown that while up to 75% of the variation in levels of some IgG glycan traits can be explained by genetics, some glycan traits are completely defined by environmental influences. Advances in both high-throughput genotyping and glycan quantification methods have enabled genome-wide association studies that are increasingly used to estimate associations of millions of single-nucleotide polymorphisms and glycosylation traits. Using this method, 18 genomic regions have so far been robustly associated with IgG N-glycosylation, discovering associations with genes encoding glycosyltransferases, but also transcription factors, co-factors, membrane transporters and other genes with no apparent role in IgG glycosylation. Further computational analyses have shown that IgG glycosylation is likely to be regulated through the expression of glycosyltransferases, but have also for the first time suggested which transcription factors are involved in the process. Moreover, it was also shown that IgG glycosylation and inflammatory diseases share common underlying causal genetic variants, suggesting that studying genetic regulation of IgG glycosylation helps not only to better understand this complex process but can also contribute to understanding why glycans are changed in disease. However, further studies are needed to unravel whether changes in IgG glycosylation are causing these diseases or the changes in the glycome are caused by the disease.

Multivariate genome-wide analysis of immunoglobulin G N-glycosylation identifies new loci pleiotropic with immune function.

The N-glycosylation of immunoglobulin G (IgG) affects its structure and function. It has been demonstrated that IgG N-glycosylation patterns are inherited as complex quantitative traits. Genome-wide association studies identified loci harboring genes encoding enzymes directly involved in protein glycosylation as well as loci likely to be involved in regulation of glycosylation biochemical pathways. Many of these loci could be linked to immune functions and risk of inflammatory and autoimmune diseases. The aim of the present study was to discover and replicate new loci associated with IgG N-glycosylation and to investigate possible pleiotropic effects of these loci onto immune function and the risk of inflammatory and autoimmune diseases. We conducted a multivariate genome-wide association analysis of 23 IgG N-glycosylation traits measured in 8090 individuals of European ancestry. The discovery stage was followed up by replication in 3147 people and in silico functional analysis. Our study increased the total number of replicated loci from 22 to 29. For the discovered loci, we suggest a number of genes potentially involved in the control of IgG N-glycosylation. Among the new loci, two (near RNF168 and TNFRSF13B) were previously implicated in rare immune deficiencies and were associated with levels of circulating immunoglobulins. For one new locus (near AP5B1/OVOL1), we demonstrated a potential pleiotropic effect on the risk of asthma. Our findings underline an important link between IgG N-glycosylation and immune function and provide new clues to understanding their interplay.

Replication of 15 loci involved in human plasma protein N-glycosylation in 4802 samples from four cohorts.

Human protein glycosylation is a complex process, and its in vivo regulation is poorly understood. Changes in glycosylation patterns are associated with many human diseases and conditions. Understanding the biological determinants of protein glycome provides a basis for future diagnostic and therapeutic applications. Genome-wide association studies (GWAS) allow to study biology via a hypothesis-free search of loci and genetic variants associated with a trait of interest. Sixteen loci were identified by three previous GWAS of human plasma proteome N-glycosylation. However, the possibility that some of these loci are false positives needs to be eliminated by replication studies, which have been limited so far. Here, we use the largest set of samples so far (4802 individuals) to replicate the previously identified loci. For all but one locus, the expected replication power exceeded 95%. Of the 16 loci reported previously, 15 were replicated in our study. For the remaining locus (near the KREMEN1 gene), the replication power was low, and hence, replication results were inconclusive. The very high replication rate highlights the general robustness of the GWAS findings as well as the high standards adopted by the community that studies genetic regulation of protein glycosylation. The 15 replicated loci present a good target for further functional studies. Among these, eight loci contain genes encoding glycosyltransferases: MGAT5, B3GAT1, FUT8, FUT6, ST6GAL1, B4GALT1, ST3GAL4 and MGAT3. The remaining seven loci offer starting points for further functional follow-up investigation into molecules and mechanisms that regulate human protein N-glycosylation in vivo.

Fc-Linked IgG N-Glycosylation in FcγR Knock-Out Mice.

Immunoglobulin G (IgG) is the most abundant immunoglobulin isotype in the blood and is involved in the pathogenesis and progression of various diseases. Glycosylation of the IgG fragment crystallizable (Fc) region is shown to vary in different physiological and pathological states. Fc N-glycan composition can alter the effector functions of IgG by modulating its affinity for ligands, such as Fcγ receptors (FcγRs). However, it is not known whether IgG glycosylation is affected by the available repertoire of FcγRs, and if the Fc-linked N-glycome can compensate for modulation of the IgG-FcγR interaction. To explore this, we examined the subclass-specific Fc IgG glycoprofiles of healthy male and female FcγR knock-out mice on C57BL/6 and BALB/c backgrounds. We observed slight changes in IgG Fc N-glycan profiles in different knock-outs; however, it seems that the strain background and sex have a stronger effect on N-glycosylation of IgG Fc regions than the FcγR repertoire.

Heritability of Human Plasma N-Glycome.

The N-glycosylation profile of total human plasma proteins could be a useful biomarker for various pathological states. Reliable high-throughput methods for such profiling have been developed. However, studies of relative importance of genetic and environmental factors in regulating plasma N-glycome are scarce. The aim of our study was to determine the role of genetic factors in phenotypic variation of plasma N-glycan profile through the estimates of its heritability. Thirty-nine total plasma N-glycome traits were analyzed in 2816 individuals from the TwinsUK data set. For the majority of the traits, high heritability estimates (>50%) were obtained pointing at a significant contribution of genetic factors in plasma N-glycome variation, especially for glycans mostly attached to immunoglobulins. We have also found several structures with higher environmental contribution to their variation.

MIgGGly (mouse IgG glycosylation analysis) - a high-throughput method for studying Fc-linked IgG N-glycosylation in mice with nanoUPLC-ESI-MS.

Immunoglobulin G (IgG) N-glycosylation is crucial for its effector functions. It is a complex trait, and large sample sets are needed to discover multiple genetic factors that underlie it. While in humans such high-throughput studies of IgG N-glycans became usual, only one has been carried out in mice. Here we describe and validate a method for the relative quantification of IgG Fc-linked N-glycans in a subclass-specific manner using nano-reverse phase liquid chromatography coupled with mass-spectrometry (nanoRP-LC-MS) applied to murine IgG. High-throughput data processing is ensured by the LaCyTools software. We have shown that IgG isolation procedure is the main source of technical variation in the current protocol. The major glycoforms were quantified reliably with coefficients of variation below 6% for all the analytes with relative abundances above 5%. We have applied our method to a sample set of 3 inbred strains: BALB/c, C57BL/6 and C3H and observed differences in subclass-specific and strain-specific N-glycosylation of IgG, suggesting a significant genetic component in the regulation of Fc-linked IgG N-glycosylation.

Profiling and genetic control of the murine immunoglobulin G glycome.

Immunoglobulin G (IgG) glycosylation is essential for function of the immune system, but the genetic and environmental factors that underlie its inter-individual variability are not well defined. The Collaborative Cross (CC) genetic resource harnesses over 90% of the common genetic variation of the mouse. By analyzing the IgG glycome composition of 95 CC strains, we made several important observations: (i) glycome variation between mouse strains was higher than between individual humans, despite all mice having the same environmental influences; (ii) five genetic loci were found to be associated with murine IgG glycosylation; (iii) variants outside traditional glycosylation site motifs affected glycome variation; (iv) bisecting N-acetylglucosamine (GlcNAc) was produced by several strains although most previous studies have reported the absence of glycans containing the bisecting GlcNAc on murine IgGs; and (v) common laboratory mouse strains are not optimal animal models for studying effects of glycosylation on IgG function.

Nuclear-cytoplasmic conflict in pea (Pisum sativum L.) is associated with nuclear and plastidic candidate genes encoding acetyl-CoA carboxylase subunits.

In crosses of wild and cultivated peas (Pisum sativum L.), nuclear-cytoplasmic incompatibility frequently occurs manifested as decreased pollen fertility, male gametophyte lethality, sporophyte lethality. High-throughput sequencing of plastid genomes of one cultivated and four wild pea accessions differing in cross-compatibility was performed. Candidate genes for involvement in the nuclear-plastid conflict were searched in the reconstructed plastid genomes. In the annotated Medicago truncatula genome, nuclear candidate genes were searched in the portion syntenic to the pea chromosome region known to harbor a locus involved in the conflict. In the plastid genomes, a substantial variability of the accD locus represented by nucleotide substitutions and indels was found to correspond to the pattern of cross-compatibility among the accessions analyzed. Amino acid substitutions in the polypeptides encoded by the alleles of a nuclear locus, designated as Bccp3, with a complementary function to accD, fitted the compatibility pattern. The accD locus in the plastid genome encoding beta subunit of the carboxyltransferase of acetyl-coA carboxylase and the nuclear locus Bccp3 encoding biotin carboxyl carrier protein of the same multi-subunit enzyme were nominated as candidate genes for main contribution to nuclear-cytoplasmic incompatibility in peas. Existence of another nuclear locus involved in the accD-mediated conflict is hypothesized.

Divergence and population traits in evolution of the genus Pisum L. as reconstructed using genes of two histone H1 subtypes showing different phylogenetic resolution.

Two histone H1 subtype genes, His7 and His5, were sequenced in a set of 56 pea accessions. Phylogenetic reconstruction based on concatenated His5 and His7 sequences had three main clades. First clade corresponded to Pisum fulvum, the next divergence separated a clade inside Pisum sativum in the broad sense that did not correspond strictly to any proposed taxonomical subdivisions. According to our estimations, the earliest divergence separating P. fulvum occurred 1.7±0.4MYA. The other divergence with high bootstrap support that separated two P. sativum groups took place approximately 1.3±0.3MYA. Thus, the main divergences in the genus took place either in late Pliocene or in early Pleistocene, the time of onset of the profound climate cooling in the northern hemisphere. The ω=K(a)/K(s) ratio was 2.5 times higher for His5 sequences than for His7. Thus, His7 gene, coding for a unique subtype specific for actively growing tissues, might have evolved under stricter evolutionary constraints than His5, that codes for a minor H1 subtype with less specific expression pattern. For this reason phylogenetic reconstructions separately obtained from His5 sequences resolved tree topology much better than those obtained from His7 sequences. Computational estimation of population dynamic parameters in the genus Pisum L. from His5-His7 sequences using IMa2 software revealed a decrease of effective population size on the early stage of Pisum evolution.

Polymorphism in a histone H1 subtype with a short N-terminal domain in three legume species (Fabaceae, Fabaeae).

A number of alleles of an orthologous gene His6 encoding histone H1 subtype f (H1-6 in pea) accumulated in chromatin of old tissues were sequenced in three legume species: seven alleles in Pisum sativum, four in Vicia unijuga and eight in Lathyrus gmelinii. In the total of 19 alleles sequenced in the three species, 29 non-synonymous substitutions and six indels were found in the coding region; most of amino acid substitutions (26 of 29) and all indels occurred in the C-terminal hydrophilic domain of the encoded protein. All species were polymorphic for some non-synonymous substitutions, V. unijuga was also polymorphic for one and P. sativum for two indels. Three near-isogenic lines of P. sativum bearing different alleles showed differences in many quantitative traits; that in the growth dynamic could be tentatively attributed to the allelic substitution of subtype H1-6. The frequencies of four electromorphs in a sampled locality of V. unijuga were found to be close to those observed 25 years ago, although their rapid change in the past was supposed in the previous study.

Phylogenetic reconstruction at the species and intraspecies levels in the genus Pisum (L.) (peas) using a histone H1 gene.

A phylogenetic analysis of the genus Pisum (peas), embracing diverse wild and cultivated forms, which evoke problems with species delimitation, was carried out based on a gene coding for histone H1, a protein that has a long and variable functional C-terminal domain. Phylogenetic trees were reconstructed on the basis of the coding sequence of the gene His5 of H1 subtype 5 in 65 pea accessions. Early separation of a clear-cut wild species Pisum fulvum is well supported, while cultivated species Pisum abyssinicum appears as a small branch within Pisum sativum. Another robust branch within P. sativum includes some wild and almost all cultivated representatives of P. sativum. Other wild representatives form diverse but rather subtle branches. In a subset of accessions, PsbA-trnH chloroplast intergenic spacer was also analysed and found less informative than His5. A number of accessions of cultivated peas from remote regions have a His5 allele of identical sequence, encoding an electrophoretically slow protein product, which earlier attracted attention as likely positively selected in harsh climate conditions. In PsbA-trnH, a 8bp deletion was found, which marks cultivated representatives of P. sativum.