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.

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.

sumSTAAR: A flexible framework for gene-based association studies using GWAS summary statistics.

Gene-based association analysis is an effective gene-mapping tool. Many gene-based methods have been proposed recently. However, their power depends on the underlying genetic architecture, which is rarely known in complex traits, and so it is likely that a combination of such methods could serve as a universal approach. Several frameworks combining different gene-based methods have been developed. However, they all imply a fixed set of methods, weights and functional annotations. Moreover, most of them use individual phenotypes and genotypes as input data. Here, we introduce sumSTAAR, a framework for gene-based association analysis using summary statistics obtained from genome-wide association studies (GWAS). It is an extended and modified version of STAAR framework proposed by Li and colleagues in 2020. The sumSTAAR framework offers a wider range of gene-based methods to combine. It allows the user to arbitrarily define a set of these methods, weighting functions and probabilities of genetic variants being causal. The methods used in the framework were adapted to analyse genes with large number of SNPs to decrease the running time. The framework includes the polygene pruning procedure to guard against the influence of the strong GWAS signals outside the gene. We also present new improved matrices of correlations between the genotypes of variants within genes. These matrices estimated on a sample of 265,000 individuals are a state-of-the-art replacement of widely used matrices based on the 1000 Genomes Project data.

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.

Varicose veins of lower extremities: Insights from the first large-scale genetic study.

Varicose veins of lower extremities (VVs) are a common multifactorial vascular disease. Genetic factors underlying VVs development remain largely unknown. Here we report the first large-scale study of VVs performed on a freely available genetic data of 408,455 European-ancestry individuals. We identified the 12 reliably associated loci that explain 13% of the SNP-based heritability, and prioritized the most likely causal genes CASZ1, PIEZO1, PPP3R1, EBF1, STIM2, HFE, GATA2, NFATC2, and SOX9. VVs-associated variants within these loci exhibited pleiotropic effects on several phenotypes including blood pressure/hypertension and blood cell traits. Gene set enrichment analysis revealed gene categories related to abnormal vasculogenesis. Genetic correlation analysis confirmed known epidemiological associations between VVs and deep venous thrombosis, weight, rough labor, and standing job, and found a genetic overlap with multiple traits that have not been previously suspected to share common genetic background with VVs. These traits included educational attainment, fluid intelligence and prospective memory scores, walking pace (negative correlation with VVs), smoking, height, number of operations, pain, and gonarthrosis (positive correlation with VVs). Finally, Mendelian randomization analysis provided evidence for causal effects of plasma levels of MICB and CD209 proteins, and anthropometric traits such as waist and hip circumference, height, weight, and both fat and fat-free mass. Our results provide novel insight into both VVs genetics and etiology. The revealed genes and proteins can be considered as good candidates for follow-up functional studies and might be of interest as potential drug targets.

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.

Defining the genetic control of human blood plasma N-glycome using genome-wide association study.

Glycosylation is a common post-translational modification of proteins. Glycosylation is associated with a number of human diseases. Defining genetic factors altering glycosylation may provide a basis for novel approaches to diagnostic and pharmaceutical applications. Here we report a genome-wide association study of the human blood plasma N-glycome composition in up to 3811 people measured by Ultra Performance Liquid Chromatography (UPLC) technology. Starting with the 36 original traits measured by UPLC, we computed an additional 77 derived traits leading to a total of 113 glycan traits. We studied associations between these traits and genetic polymorphisms located on human autosomes. We discovered and replicated 12 loci. This allowed us to demonstrate an overlap in genetic control between total plasma protein and IgG glycosylation. The majority of revealed loci contained genes that encode enzymes directly involved in glycosylation (FUT3/FUT6, FUT8, B3GAT1, ST6GAL1, B4GALT1, ST3GAL4, MGAT3 and MGAT5) and a known regulator of plasma protein fucosylation (HNF1A). However, we also found loci that could possibly reflect other more complex aspects of glycosylation process. Functional genomic annotation suggested the role of several genes including DERL3, CHCHD10, TMEM121, IGH and IKZF1. The hypotheses we generated may serve as a starting point for further functional studies in this research area.

Genome-wide association study identifies RNF123 locus as associated with chronic widespread musculoskeletal pain.

BACKGROUND AND OBJECTIVES: Chronic widespread musculoskeletal pain (CWP) is a symptom of fibromyalgia and a complex trait with poorly understood pathogenesis. CWP is heritable (48%-54%), but its genetic architecture is unknown and candidate gene studies have produced inconsistent results. We conducted a genome-wide association study to get insight into the genetic background of CWP. METHODS: Northern Europeans from UK Biobank comprising 6914 cases reporting pain all over the body lasting >3 months and 242 929 controls were studied. Replication of three independent genome-wide significant single nucleotide polymorphisms was attempted in six independent European cohorts (n=43 080; cases=14 177). Genetic correlations with risk factors, tissue specificity and colocalisation were examined. RESULTS: Three genome-wide significant loci were identified (rs1491985, rs10490825, rs165599) residing within the genes Ring Finger Protein 123 (RNF123), ATPase secretory pathway Ca2+transporting 1 (ATP2C1) and catechol-O-methyltransferase (COMT). The RNF123 locus was replicated (meta-analysis p=0.0002), the ATP2C1 locus showed suggestive association (p=0.0227) and the COMT locus was not replicated. Partial genetic correlation between CWP and depressive symptoms, body mass index, age of first birth and years of schooling were identified. Tissue specificity and colocalisation analysis highlight the relevance of skeletal muscle in CWP. CONCLUSIONS: We report a novel association of RNF123 locus and a suggestive association of ATP2C1 locus with CWP. Both loci are consistent with a role of calcium regulation in CWP. The association with COMT, one of the most studied genes in chronic pain field, was not confirmed in the replication analysis.

Genome-wide meta-analysis of 158,000 individuals of European ancestry identifies three loci associated with chronic back pain.

Back pain is the #1 cause of years lived with disability worldwide, yet surprisingly little is known regarding the biology underlying this symptom. We conducted a genome-wide association study (GWAS) meta-analysis of chronic back pain (CBP). Adults of European ancestry were included from 15 cohorts in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and from the UK Biobank interim data release. CBP cases were defined as those reporting back pain present for ≥3-6 months; non-cases were included as comparisons ("controls"). Each cohort conducted genotyping using commercially available arrays followed by imputation. GWAS used logistic regression models with additive genetic effects, adjusting for age, sex, study-specific covariates, and population substructure. The threshold for genome-wide significance in the fixed-effect inverse-variance weighted meta-analysis was p<5×10(-8). Suggestive (p<5×10(-7)) and genome-wide significant (p<5×10(-8)) variants were carried forward for replication or further investigation in the remaining UK Biobank participants not included in the discovery sample. The discovery sample comprised 158,025 individuals, including 29,531 CBP cases. A genome-wide significant association was found for the intronic variant rs12310519 in SOX5 (OR 1.08, p = 7.2×10(-10)). This was subsequently replicated in 283,752 UK Biobank participants not included in the discovery sample, including 50,915 cases (OR 1.06, p = 5.3×10(-11)), and exceeded genome-wide significance in joint meta-analysis (OR 1.07, p = 4.5×10(-19)). We found suggestive associations at three other loci in the discovery sample, two of which exceeded genome-wide significance in joint meta-analysis: an intergenic variant, rs7833174, located between CCDC26 and GSDMC (OR 1.05, p = 4.4×10(-13)), and an intronic variant, rs4384683, in DCC (OR 0.97, p = 2.4×10(-10)). In this first reported meta-analysis of GWAS for CBP, we identified and replicated a genetic locus associated with CBP (SOX5). We also identified 2 other loci that reached genome-wide significance in a 2-stage joint meta-analysis (CCDC26/GSDMC and DCC).

Candidate genes for age at menarche are associated with endometriosis.

RESEARCH QUESTION: Are the candidate genes for age at menarche associated with a risk of endometriosis? DESIGN: Fifty-two candidate single nucleotide polymorphisms (SNP) for age at menarche, their gene-gene and gene-environment interactions were analysed for possible association with endometriosis in a sample of 395 patients and 981 controls. Association of the polymorphisms was analysed using logistic regression according to three main genetic models (additive, recessive and dominant). The gene-gene and gene-environment interactions were analysed for the second-, third- and fourth-order models with adjustment for covariates and multiple comparisons with subsequent cross-validation. RESULTS: Sixteen SNP for age at menarche out of the 52 studied were associated with endometriosis. Polymorphism rs6589964 BSX was associated with endometriosis according to the additive and recessive models (OR 1.27-1.47, Pperm ≤ 0.006). Fourteen SNP were associated with the disease within 12 most significant models of gene-gene interactions (Pperm ≤ 0.008). Twelve SNP involved in 10 most significant models of SNP-induced abortion interactions are associated with endometriosis. Fourteen of the 16 polymorphisms associated with endometriosis demonstrated pleiotropic effects: they were also associated with either age at menarche (7 SNP) or height and/or body mass index (10 SNP) in the studied sample. The 16 SNP associated with endometriosis and 316 SNP linked to them have regulatory and expression quantitative trait locus significance for 28 genes contributing to the G alpha signal pathway (fold enrichment 31.09, PFDR = 0.001) and responses to endogenous stimuli (fold enrichment 16.01, PFDR = 0.027). CONCLUSIONS: Sixteen SNP for age at menarche out of the 52 studied were associated with endometriosis.

Male Meiotic Recombination in the Steppe Agama, Trapelus sanguinolentus (Agamidae, Iguania, Reptilia).

Meiotic recombination rates and patterns of crossover distributions along the chromosomes vary considerably even between closely related species. The adaptive significance of these differences is still unclear due to the paucity of empirical data. Most data on recombination come from mammalian species, while other vertebrate clades are poorly explored. Using immunolocalization of the protein of the lateral element of the synaptonemal complex (SYCP3) and the mismatch-repair protein MLH1, which marks mature recombination nodules, we analyzed recombination rates and crossover distribution in meiotic prophase chromosomes of the steppe agama (Trapelus sanguinolentus, Agamidae, Acrodonta, Iguania) and compared them with data obtained for the genus Anolis (Dactyloidae, Pleurodonta, Iguania). We found that, despite a smaller genome size, the total SC length and the MLH1 focus number per cell are much higher in the agama than in the anoles. The distributions of the MLH1 foci in the agama are multimodal in larger chromosomes and bimodal in smaller chromosomes without a significant centromere effect, resembling the patterns known for birds. A possible relationship between karyotype remodeling and the evolution of recombination in Iguania is discussed.

The association of lumbar intervertebral disc degeneration with low back pain is modified by underlying genetic propensity to pain.

BACKGROUND CONTEXT: Associations between magnetic resonance imaging (MRI)-detected lumbar intervertebral disc degeneration (LDD) and LBP are often of modest magnitude. This association may be larger in specific patient subgroups. PURPOSE: To examine whether the association between LDD and LBP is modified by underlying genetic predispositions to pain. STUDY DESIGN: Cross-sectional study in UK Biobank (UKB) and TwinsUK. PATIENT SAMPLES: A genome-wide association study (GWAS) of the number of anatomical chronic pain locations was conducted in 347,538 UKB participants. The GWAS was used to develop a genome-wide polygenic risk score (PRS) in a holdout sample of 30,000 UKB participants. The PRS model was then used in analyses of 645 TwinsUK participants with standardized LDD MRI assessments. OUTCOME MEASURES: Ever having had LBP associated with disability lasting ≥1 month (LBP1). METHODS: Using the PRS as a proxy for "genetically-predicted propensity to pain", we stratified TwinsUK participants into PRS quartiles. A "basic" model examined the association between an LDD summary score (LSUM) and LBP1, adjusting for covariates. A "fully-adjusted" model also adjusted for PRS quartile and LSUM x PRS quartile interaction terms. RESULTS: In the basic model, the odds ratio (OR) of LBP1 was 1.8 per standard deviation of LSUM (95% confidence interval [CI] 1.4 -2.3). In the fully-adjusted model, there was a statistically significant LSUM-LBP1 association in quartile 4, the highest PRS quartile (OR = 2.5 [95% CI 1.7-3.7], p=2.6×10-6), and in quartile 3 (OR=2.0, [95% CI 1.3-3.0]; p=0.002), with small-magnitude and/or non-significant associations in the lowest two PRS quartiles. PRS quartile was a significant effect modifier of the LSUM-LBP1 association (interaction p≤0.05). CONCLUSIONS: Genetically-predicted propensity to pain modifies the LDD-LBP association, with the strongest association present in people with the highest genetic propensity to pain. Lumbar MRI findings may have stronger connections to LBP in specific subgroups of people.

Noncoding rare variants in PANX3 are associated with chronic back pain.

ABSTRACT: Back pain is the leading cause of years lived with disability worldwide, yet surprisingly, little is known regarding the biology underlying this condition. The impact of genetics is known for chronic back pain: its heritability is estimated to be at least 40%. Large genome-wide association studies have shown that common variation may account for up to 35% of chronic back pain heritability; rare variants may explain a portion of the heritability not explained by common variants. In this study, we performed the first gene-based association analysis of chronic back pain using UK Biobank imputed data including rare variants with moderate imputation quality. We discovered 2 genes, SOX5 and PANX3 , influencing chronic back pain. The SOX5 gene is a well-known back pain gene. The PANX3 gene has not previously been described as having a role in chronic back pain. We showed that the association of PANX3 with chronic back pain is driven by rare noncoding intronic polymorphisms. This result was replicated in an independent sample from UK Biobank and validated using a similar phenotype, dorsalgia, from FinnGen Biobank. We also found that the PANX3 gene is associated with intervertebral disk disorders. We can speculate that a possible mechanism of action of PANX3 on back pain is due to its effect on the intervertebral disks.

Repurposing Antihypertensive and Statin Medications for Spinal Pain: A Mendelian Randomization Study.

STUDY DESIGN: Mendelian randomization (MR) study. OBJECTIVE: To examine whether antihypertensive medications (beta-blockers, calcium channel blockers, and angiotensin-converting enzyme inhibitors) and statins can be repurposed to prevent or treat spinal pain (back or neck pain). SUMMARY OF BACKGROUND DATA: Observational studies and a recent MR study have found associations between elevated blood pressure and a greater risk of back pain. Observational studies have found associations between hyperlipidemia and statin use and greater risk of back pain. No prior MR studies have examined the effects of antihypertensives or statins on spinal pain. MATERIALS AND METHODS: This was a two-sample MR study using publicly available summary statistics from large-scale genome-wide association studies (GWAS). Sample sizes in exposure GWASs were n=757,601 (systolic blood pressure) and n=173,082 (low-density lipoprotein cholesterol), and n=1,028,947 for the outcome GWAS of spinal pain defined as health care seeking for any spinal pain-related diagnosis. Genes and cis-acting variants were identified as proxies for the drug targets of interest. MR analyses used inverse-variance weighted meta-analysis. The threshold for statistical significance after correction for multiple testing was P <0.0125. RESULTS: No statistically significant associations of these medications with spinal pain were found. However, findings were suggestive of a protective effect of beta-blockers on spinal pain risk (odds ratio [OR] 0.84, 95% confidence interval [CI] 0.72-0.98; P =0.03), and calcium channel blockers on greater spinal pain risk (OR 1.12, 95% CI 1.02-1.24; P =0.02). CONCLUSIONS: A protective effect of beta-blockers on spinal pain was suggested in the current study, consistent with findings from observational studies of various other pain phenotypes. The detrimental effect of calcium channel blockers on spinal pain suggested in the current study must be interpreted in the context of conflicting directions of effect on nonspinal pain phenotypes in other observational studies.This Mendelian randomization study examined whether antihypertensive medications (beta-blockers, calcium channel blockers, and angiotensin-converting enzyme inhibitors) and statins can be repurposed to prevent or treat spinal.This was a two-sample MR study using publicly available summary statistics from large-scale genome-wide association studies ranging size from 173,082 to 1,028,947 adults.While no statistically significant associations were found, a protective effect of beta-blockers on spinal pain was suggested (odds ratio [OR] 0.84, 95% confidence interval [CI] 0.72 to 0.98; p= 0.03), as was a detrimental effect of calcium channel blockers on spinal pain (OR 1.12, 95% CI 1.02 to 1.24; p= 0.02).

Loci Associated with Negative Heterosis for Viability and Meat Productivity in Interspecific Sheep Hybrids.

Negative heterosis can occur on different economically important traits, but the exact biological mechanisms of this phenomenon are still unknown. The present study focuses on determining the genetic factors associated with negative heterosis in interspecific hybrids between domestic sheep (Ovis aries) and argali (Ovis ammon). One locus (rs417431015) associated with viability and two loci (rs413302370, rs402808951) associated with meat productivity were identified. One gene (ARAP2) was prioritized for viability and three for meat productivity (PDE2A, ARAP1, and PCDH15). The loci associated with meat productivity were demonstrated to fit the overdominant inheritance model and could potentially be involved int negative heterosis mechanisms.

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.

Genotype-by-environment interactions in chronic back pain.

BACKGROUND CONTEXT: Chronic back pain (CBP) is a common debilitating condition with substantial societal impact. While understanding genotype-by-environment (GxE) interactions may be crucial to achieving the goals of personalized medicine, there are few large-scale studies investigating this topic for CBP. None of them systematically explore multiple CBP risk factors. PURPOSE: To estimate the extent to which genetic effects on CBP are modified by known demographic and clinical risk factors. RESEARCH DESIGN: Case-control study, genome-wide GxE interaction study. PATIENT SAMPLE: Data on up to 331,610 unrelated participants (57,881 CBP cases and 273,729 controls) from the UK Biobank cohort were used. UK Biobank is a prospective cohort with collected deep genetic and phenotypic data on approximately 500,000 individuals across the UK. OUTCOME MEASURES: Self-reported chronic back pain. METHODS: We applied a whole-genome approach to estimate the proportion of phenotypic variance explained by interactions between genotype and 12 known risk factors. We also analyzed if effects of common single-nucleotide polymorphisms on CBP are changed in presence of known risk factors. RESULTS: The results indicate a modest, if any, modification of genetic effects by examined risk factors in CBP. Our estimates suggest that detecting such weak effects would require a sample size of millions of individuals. CONCLUSIONS: The GxE interactions with examined common risk factors for CBP are either weak or absent. Interactions of such magnitude are unlikely to have the potential to inform and influence treatment strategies. Risk estimation models may use common genetic variation and the considered risk factors as independent predictors, without accounting for GxE.

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.

Population Structure and Genetic Diversity of the 175 Soybean Breeding Lines and Varieties Cultivated in West Siberia and Other Regions of Russia.

Soybean is a leguminous plant cultivated in many countries and is considered important in the food industry due to the high levels of oil and protein content in the beans. The high demand for soybeans and its products in the industry requires the expansion of cultivation areas. Despite climatic restrictions, West Siberia is gradually expanding its area of soybean cultivation. In this study, we present the first analysis of the population structure and genetic diversity of the 175 soybean Glycine max breeding lines and varieties cultivated in West Siberia (103 accessions) and other regions of Russia (72 accessions), and we compare them with the cultivated soybean varieties from other geographical locations. Principal component analysis revealed several genetic clusters with different levels of genetic heterogeneity. Studied accessions are genetically similar to varieties from China, Japan, and the USA and are genetically distant to varieties from South Korea. Admixture analysis revealed four ancestry groups based on genetic ancestry and geographical origin, which are consistent with the regions of cultivation and origin of accessions and correspond to the principal component analysis result. Population statistics, including nucleotide diversity, Tajima's D, and linkage disequilibrium, are comparatively similar to those observed for studied accessions of a different origin. This study provides essential population and genetic information about the unique collection of breeding lines and varieties cultivated in West Siberia and other Russian regions to foster further evolutionary, genome-wide associations and functional breeding studies.

Multivariate Genome-Wide Association Study of Concentrations of Seven Elements in Seeds Reveals Four New Loci in Russian Wheat Lines.

Wheat is a cereal grain that plays an important role in the world's food industry. The identification of the loci that change the concentration of elements in wheat seeds is an important challenge nowadays especially for genomic selection and breeding of novel varieties. In this study, we performed a multivariate genome-wide association study (GWAS) of the seven traits-concentrations of Zn, Mg, Mn, Ca, Cu, Fe, and K in grain-of the Russian collection of common wheat Triticum aestivum (N = 149 measured in two years in two different fields). We replicated one known locus associated with the concentration of Zn (IAAV1375). We identified four novel loci-BS00022069_51 (associated with concentrations of Ca and K), RFL_Contig6053_3082 (associated with concentrations of Fe and Mn), Kukri_rep_c70864_329 (associated with concentrations of all elements), and IAAV8416 (associated with concentrations of Fe and Mn)-three of them were located near the genes TraesCS6A02G375400, TraesCS7A02G094800, and TraesCS5B02G325400. Our result adds novel information on the loci involved in wheat grain element contents and may be further used in genomic selection.