Mapping of the gene network that regulates glycan clock of ageing.

Glycans are an essential structural component of immunoglobulin G (IgG) that modulate its structure and function. However, regulatory mechanisms behind this complex posttranslational modification are not well known. Previous genome-wide association studies (GWAS) identified 29 genomic regions involved in regulation of IgG glycosylation, but only a few were functionally validated. One of the key functional features of IgG glycosylation is the addition of galactose (galactosylation), a trait which was shown to be associated with ageing. We performed GWAS of IgG galactosylation (N=13,705) and identified 16 significantly associated loci, indicating that IgG galactosylation is regulated by a complex network of genes that extends beyond the galactosyltransferase enzyme that adds galactose to IgG glycans. Gene prioritization identified 37 candidate genes. Using a recently developed CRISPR/dCas9 system we manipulated gene expression of candidate genes in the in vitro IgG expression system. Upregulation of three genes, EEF1A1, MANBA and TNFRSF13B, changed the IgG glycome composition, which confirmed that these three genes are involved in IgG galactosylation in this in vitro expression system.

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.

Development and Replication of a Genome-Wide Polygenic Risk Score for Chronic Back Pain.

Chronic back pain (CBP) is a complex heritable trait and a major cause of disability worldwide. We developed and validated a genome-wide polygenic risk score (PRS) for CBP using a large-scale GWAS based on UK Biobank participants of European ancestry (N = 265,000). The PRS showed poor overall predictive ability (AUC = 0.56 and OR = 1.24 per SD, 95% CI: 1.22-1.26), but individuals from the 99th percentile of PRS distribution had a nearly two-fold increased risk of CBP (OR = 1.82, 95% CI: 1.60-2.06). We validated the PRS on an independent TwinsUK sample, obtaining a similar magnitude of effect. The PRS was significantly associated with various ICD-10 and OPCS-4 diagnostic codes, including chronic ischemic heart disease (OR = 1.1, p-value = 4.8 × 10-15), obesity, metabolism-related traits, spine disorders, disc degeneration, and arthritis-related disorders. PRS and environment interaction analysis with twelve known CBP risk factors revealed no significant results, suggesting that the magnitude of G × E interactions with studied factors is small. The limited predictive ability of the PRS that we developed is likely explained by the complexity, heterogeneity, and polygenicity of CBP, for which sample sizes of a few hundred thousand are insufficient to estimate small genetic effects robustly.

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.

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.

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.

Method for the Isolation of "RNA-seq-Quality" RNA from Human Intervertebral Discs after Mortar and Pestle Homogenization.

The problem of isolating high-quality total RNA from intervertebral discs has no recognized solution yet. This is due to the extremely low content of live cells in the samples and the voluminous intercellular matrix. A variety of published protocols focused on isolating RNA from articular cartilage have recommended the use of expensive equipment, enzymatic matrix cleavage, or cell culture. In our study, we used a combination of the traditional QIAzol protocol (Qiagen, Germany) and RNEasy column purification (Qiagen, Germany) to obtain high-quality RNA from post-surgical intervertebral disc fragments. Only a mortar and a pestle were used for grinding, making our method particularly accessible. The isolated RNA with a RIN of ~7 is suitable for studying the expression profile of chondrocytes in situ. RNA-seq analysis of three samples demonstrated cell type ratios to be mostly relevant to intervertebral disc tissues, with over 70% of the chondrocytes of the three subtypes having an admixture of blood-related cells.

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.

Association Between Human Gut Microbiome and N-Glycan Composition of Total Plasma Proteome.

Being one of the most dynamic entities in the human body, glycosylation of proteins fine-tunes the activity of the organismal machinery, including the immune system, and mediates the interaction with the human microbial consortium, typically represented by the gut microbiome. Using data from 194 healthy individuals, we conducted an associational study to uncover potential relations between the gut microbiome and the blood plasma N-glycome, including N-glycome of immunoglobulin G. While lacking strong linkages on the multivariate level, we were able to identify associations between alpha and beta microbiome diversity and the blood plasma N-glycome profile. Moreover, for two bacterial genera, namely, Bilophila and Clostridium innocuum, significant associations with specific glycans were also shown. The study's results suggest a non-trivial, possibly weak link between the total plasma N-glycome and the gut microbiome, predominantly involving glycans related to the immune system proteins, including immunoglobulin G. Further studies of glycans linked to microbiome-related proteins in well-selected patient groups are required to conclusively establish specific associations.

A protocol for recruiting and analyzing the disease-oriented Russian disc degeneration study (RuDDS) biobank for functional omics studies of lumbar disc degeneration.

Lumbar intervertebral disc degeneration (DD) disease is one of the main risk factors for low back pain and a leading cause of population absenteeism and disability worldwide. Despite a variety of biological studies, lumbar DD is not yet fully understood, partially because there are only few studies that use systematic and integrative approaches. This urges the need for studies that integrate different omics (including genomics and transcriptomics) measured on samples within a single cohort. This protocol describes a disease-oriented Russian disc degeneration study (RuDDS) biobank recruitment and analyses aimed to facilitate further omics studies of lumbar DD integrating genomic, transcriptomic and glycomic data. A total of 1,100 participants aged over 18 with available lumbar MRI scans, medical histories and biological material (whole blood, plasma and intervertebral disc tissue samples from surgically treated patients) will be enrolled during the three-year period from two Russian clinical centers. Whole blood, plasma and disc tissue specimens will be used for genotyping with genome-wide SNP-arrays, glycome profiling and RNA sequencing, respectively. Omics data will be further used for a genome-wide association study of lumbar DD with in silico functional annotation, analysis of plasma glycome and lumbar DD disease interactions and transcriptomic data analysis including an investigation of differential expression patterns associated with lumbar DD disease. Statistical tests applied in each of the analyses will meet the standard criteria specific to the attributed study field. In a long term, the results of the study will expand fundamental knowledge about lumbar DD development and contribute to the elaboration of novel personalized approaches for disease prediction and therapy. Additionally to the lumbar disc degeneration study, a RuDDS cohort could be used for other genetic studies, as it will have unique omics data. Trial registration number NCT04600544.

ANANASTRA: annotation and enrichment analysis of allele-specific transcription factor binding at SNPs.

We present ANANASTRA, https://ananastra.autosome.org, a web server for the identification and annotation of regulatory single-nucleotide polymorphisms (SNPs) with allele-specific binding events. ANANASTRA accepts a list of dbSNP IDs or a VCF file and reports allele-specific binding (ASB) sites of particular transcription factors or in specific cell types, highlighting those with ASBs significantly enriched at SNPs in the query list. ANANASTRA is built on top of a systematic analysis of allelic imbalance in ChIP-Seq experiments and performs the ASB enrichment test against background sets of SNPs found in the same source experiments as ASB sites but not displaying significant allelic imbalance. We illustrate ANANASTRA usage with selected case studies and expect that ANANASTRA will help to conduct the follow-up of GWAS in terms of establishing functional hypotheses and designing experimental verification.

Genetic regulation of post-translational modification of two distinct proteins.

Post-translational modifications diversify protein functions and dynamically coordinate their signalling networks, influencing most aspects of cell physiology. Nevertheless, their genetic regulation or influence on complex traits is not fully understood. Here, we compare the genetic regulation of the same PTM of two proteins - glycosylation of transferrin and immunoglobulin G (IgG). By performing genome-wide association analysis of transferrin glycosylation, we identify 10 significantly associated loci, 9 of which were not reported previously. Comparing these with IgG glycosylation-associated genes, we note protein-specific associations with genes encoding glycosylation enzymes (transferrin - MGAT5, ST3GAL4, B3GAT1; IgG - MGAT3, ST6GAL1), as well as shared associations (FUT6, FUT8). Colocalisation analyses of the latter suggest that different causal variants in the FUT genes regulate fucosylation of the two proteins. Glycosylation of these proteins is thus genetically regulated by both shared and protein-specific mechanisms.

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.

Quantitative Genetics of Human Protein N-Glycosylation.

Although changes in protein glycosylation are observed in a wide range of diseases and pathological states, the examples of use of glycans as biomarkers and therapeutic targets are limited. This is not in small part because the understanding of human glycome regulation in vivo is incomplete and fragmented. Combination of human glycomics and genomics offers a powerful "data-driven hypotheses" approach to dissect the complex human glycobiology in vivo in an agnostic manner.In this chapter we review a decade of quantitative genetic studies of human N-glycome, including studies of its heritability and gene-mapping via genome-wide association studies (GWASs). We show that GWASs of human N-glycome start revealing regulators of the biochemical network of N-glycosylation. Some of these regulators demonstrate pleiotropic effects on human disease, especially autoimmune and inflammatory. We emphasize the use of in silico functional methods and multi-omics approaches to prioritize functional mechanisms to be further validated in laboratory experiments. This combined approach will lead to better understanding of mechanisms of regulation of human protein glycosylation and will provide a rich source of etiologic insight, therapeutic interventions, and biomarkers.

Identification of tissue-specific and common methylation quantitative trait loci in healthy individuals using MAGAR.

BACKGROUND: Understanding the influence of genetic variants on DNA methylation is fundamental for the interpretation of epigenomic data in the context of disease. There is a need for systematic approaches not only for determining methylation quantitative trait loci (methQTL), but also for discriminating general from cell type-specific effects. RESULTS: Here, we present a two-step computational framework MAGAR ( https://bioconductor.org/packages/MAGAR ), which fully supports the identification of methQTLs from matched genotyping and DNA methylation data, and additionally allows for illuminating cell type-specific methQTL effects. In a pilot analysis, we apply MAGAR on data in four tissues (ileum, rectum, T cells, B cells) from healthy individuals and demonstrate the discrimination of common from cell type-specific methQTLs. We experimentally validate both types of methQTLs in an independent data set comprising additional cell types and tissues. Finally, we validate selected methQTLs located in the PON1, ZNF155, and NRG2 genes by ultra-deep local sequencing. In line with previous reports, we find cell type-specific methQTLs to be preferentially located in enhancer elements. CONCLUSIONS: Our analysis demonstrates that a systematic analysis of methQTLs provides important new insights on the influences of genetic variants to cell type-specific epigenomic variation.

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.

The limits of normal approximation for adult height.

Adult height inspired the first biometrical and quantitative genetic studies and is a test-case trait for understanding heritability. The studies of height led to formulation of the classical polygenic model, that has a profound influence on the way we view and analyse complex traits. An essential part of the classical model is an assumption of additivity of effects and normality of the distribution of the residuals. However, it may be expected that the normal approximation will become insufficient in bigger studies. Here, we demonstrate that when the height of hundreds of thousands of individuals is analysed, the model complexity needs to be increased to include non-additive interactions between sex, environment and genes. Alternatively, the use of log-normal approximation allowed us to still use the additive effects model. These findings are important for future genetic and methodologic studies that make use of adult height as an exemplar trait.

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.

Genome-wide association studies of low back pain and lumbar spinal disorders using electronic health record data identify a locus associated with lumbar spinal stenosis.

ABSTRACT: Identifying genetic risk factors for lumbar spine disorders may lead to knowledge regarding underlying mechanisms and the development of new treatments. We conducted a genome-wide association study involving 100,811 participants with genotypes and longitudinal electronic health record data from the Electronic Medical Records and Genomics Network and Geisinger Health. Cases and controls were defined using validated algorithms and clinical diagnostic codes. Electronic health record-defined phenotypes included low back pain requiring healthcare utilization (LBP-HC), lumbosacral radicular syndrome (LSRS), and lumbar spinal stenosis (LSS). Genome-wide association study used logistic regression with additive genetic effects adjusting for age, sex, site-specific factors, and ancestry (principal components). A fixed-effect inverse-variance weighted meta-analysis was conducted. Genetic variants of genome-wide significance (P < 5 × 10-8) were carried forward for replication in an independent sample from UK Biobank. Phenotype prevalence was 48.8% for LBP-HC, 19.8% for LSRS, and 7.9% for LSS. No variants were significantly associated with LBP-HC. One locus was associated with LSRS (lead variant rs146153280:C>G, odds ratio [OR] = 1.17 for G, P = 2.1 × 10-9), but was not replicated. Another locus on chromosome 2 spanning GFPT1, NFU1, and AAK1 was associated with LSS (lead variant rs13427243:G>A, OR = 1.10 for A, P = 4.3 × 10-8) and replicated in UK Biobank (OR = 1.11, P = 5.4 × 10-5). This was the first genome-wide association study meta-analysis of lumbar spinal disorders using electronic health record data. We identified 2 novel associations with LSRS and LSS; the latter was replicated in an independent sample.

Nontrivial Replication of Loci Detected by Multi-Trait Methods.

The ever-growing genome-wide association studies (GWAS) have revealed widespread pleiotropy. To exploit this, various methods that jointly consider associations of a genetic variant with multiple traits have been developed. Most efforts have been made concerning improving GWAS discovery power. However, how to replicate these discovered pleiotropic loci has yet to be discussed thoroughly. Unlike a single-trait scenario, multi-trait replication is not trivial considering the underlying genotype-multi-phenotype map of the associations. Here, we evaluate four methods for replicating multi-trait associations, corresponding to four levels of replication strength. Weak replication cannot justify pleiotropic genetic effects, whereas strong replication using our developed correlation methods can inform consistent pleiotropic genetic effects across the discovery and replication samples. We provide a protocol for replicating multi-trait genetic associations in practice. The described methods are implemented in the free and open-source R package MultiABEL.