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.

In silico genome-wide gene-based association analysis reveals new genes predisposing to coronary artery disease.

Genome-wide association study (GWAS) have identified more than 300 single nucleotide polymorphisms at 163 independent loci associated with coronary artery disease (CAD). However, there is no full understanding about the causal genes for CAD and the mechanisms of their action. We aimed to perform a post GWAS analysis to identify genes whose polymorphism may influence the risk of CAD. Using the UK Biobank GWAS summary statistics, we performed a gene-based association analysis. We found 63 genes significantly associated with CAD due to their within-gene polymorphisms. Many of these genes are well known. Some known CAD genes such as FURIN and SORT1 did not show the gene-based association because their variants had low GWAS signals or gene-based association was inflated by the strong GWAS signal outside the gene. For several known CAD genes, we demonstrated that their effects could be explained not only or not at all by their own variants but by the variants within the neighboring genes controlling their expression. Using several bioinformatics techniques, we suggested potential mechanisms underlying gene-CAD associations. Three genes, CDK19, NCALD, and ARHGEF12 were not previously associated with CAD. The role of these genes should be clarified in further studies.

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.

ISSLS Prize in Clinical Science 2020. Examining causal effects of body mass index on back pain: a Mendelian randomization study.

PURPOSE: Measures of body fat accumulation are associated with back pain, but a causal association is unclear. We hypothesized that BMI would have causal effects on back pain. We conducted a two-sample Mendelian randomization (MR) study to assess the causal effect of body mass index (BMI) on the outcomes of (1) back pain and (2) chronic back pain (duration > 3 months). METHODS: We identified genetic instrumental variables for BMI (n = 60 variants) from a meta-analysis of genome-wide association studies (GWAS) conducted by the Genetic Investigation of ANthropometric Traits consortium in individuals of European ancestry (n = 322,154). We conducted GWAS of back pain and chronic back pain (n = 453,860) in a non-overlapping sample of individuals of European ancestry. We used inverse-variance weighted (IVW) meta-analysis as the primary method to estimate causal effects. RESULTS: The IVW analysis showed evidence supporting a causal association of BMI on back pain, with a 1-standard deviation (4.65 kg/m2) increase in BMI conferring 1.15 times the odds of back pain (95% confidence interval [CI]: 1.06-1.25, p = 0.001]; effects were directionally consistent in secondary analysis and sensitivity analyses. The IVW analysis supported a causal association of BMI on chronic back pain (OR 1.20 per 1 SD deviation increase in BMI [95% CI 1.09-1.32; p = 0.0002]), and effects were directionally consistent in secondary analysis and sensitivity analyses. CONCLUSION: In this first MR study of BMI and back pain, we found a significant causal effect of BMI on both back pain and chronic back pain. These slides can be retrieved under Electronic Supplementary Material.

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.

A New Method for Conditional Gene-Based Analysis Effectively Accounts for the Regional Polygenic Background.

Gene-based association analysis is a powerful tool for identifying genes that explain trait variability. An essential step of this analysis is a conditional analysis. It aims to eliminate the influence of SNPs outside the gene, which are in linkage disequilibrium with intragenic SNPs. The popular conditional analysis method, GCTA-COJO, accounts for the influence of several top independently associated SNPs outside the gene, correcting the z statistics for intragenic SNPs. We suggest a new TauCOR method for conditional gene-based analysis using summary statistics. This method accounts the influence of the full regional polygenic background, correcting the genotype correlations between intragenic SNPs. As a result, the distribution of z statistics for intragenic SNPs becomes conditionally independent of distribution for extragenic SNPs. TauCOR is compatible with any gene-based association test. TauCOR was tested on summary statistics simulated under different scenarios and on real summary statistics for a 'gold standard' gene list from the Open Targets Genetics project. TauCOR proved to be effective in all modelling scenarios and on real data. The TauCOR's strategy showed comparable sensitivity and higher specificity and accuracy than GCTA-COJO on both simulated and real data. The method can be successfully used to improve the effectiveness of gene-based association analyses.

No evidence for causal effects of C-reactive protein (CRP) on chronic pain conditions: a Mendelian randomization study.

Objective: We conducted a Mendelian randomization (MR) study to examine causal associations of C-reactive protein (CRP) with (1) spinal pain; (2) extent of multisite chronic pain; and (3) chronic widespread musculoskeletal pain. Design: Two-sample MR study. Setting/Subjects: We used summary statistics from publicly available genome-wide association studies (GWAS) conducted in multiple cohorts and biobanks. Genetic instrumental variables were taken from an exposure GWAS of CRP (n=204,402). Outcome GWASs examined spinal pain (n=1,028,947), extent of multisite chronic pain defined as the number of locations with chronic pain (n=387,649), and chronic widespread pain (n=249,843). Methods: We examined MR evidence for causal associations using inverse-variance weighted (IVW) analysis and sensitivity analyses using other methods. We calculated odds ratios (ORs), 95% confidence intervals (95% CIs), and p-values, using a Bonferroni correction (p<0.0166) to account for 3 primary comparisons. Results: Greater serum CRP (mg/L) was not significantly causally associated with spinal pain (OR=1.04, 95% CI 1.00-1.08; p=0.07) in IVW analysis. Greater serum CRP also showed no significant causal association with extent of multisite chronic pain in IVW analysis (beta coefficient= 0.014, standard error=0.011; p=0.19). CRP also showed no significant causal association with chronic widespread pain in IVW analysis (OR=1.00, 95% CI 1.00-1.00; p=0.75). All secondary and sensitivity analyses also showed no significant associations. Conclusions: This MR study found no causal association of CRP on spinal pain, the extent of chronic pain, or chronic widespread pain. Future studies examining mechanistic biomarkers for pain conditions should consider other candidates besides CRP.

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 Twins UK. 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=.002), with small-magnitude and/or nonsignificant associations in the lowest 2 PRS quartiles. PRS quartile was a significant effect modifier of the LSUM-LBP1 association (interaction p≤.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.

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).

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.

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.

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.