Estimation of effect size distribution from genome-wide association studies and implications for future discoveries. Nat Genet. 2010;42:570-575. [PMID: 20562874 DOI: 10.1038/ng.610]

Association between kidney stones and urological cancers: results from the NHANES 2007-2020 and Mendelian randomization study

BACKGROUND

Kidney stones is a common urological disease with a rising incidence in global. The association between kidney stones and urological cancers remains controversial. This study utilized the data from the 2007-2020 National Health and Nutrition Examination Survey (NHANES) and Mendelian randomization (MR) analysis to evaluate the association and potential causal relationship between kidney stones and renal cell carcinoma, bladder cancer, and prostate cancer.

METHODS

Multivariate logistic regression was used to examine the association between kidney stones history and urological cancers, followed by stratified analyses. Subsequently, causal relationships were explored via the inverse variance weighted (IVW), weighted median, and MR-Egger methods. Sensitivity analyses were performed to ensure the robustness of the findings.

RESULTS

Data from 13,013 individuals (5,138 males) were analyzed. Kidney stones was significantly associated with an increased risk of renal cell carcinoma (OR = 1.92, 95% CI 1.90-1.95, P < 0.001), bladder cancer (OR = 2.749, 95% CI 2.71-2.78, P < 0.001), and prostate cancer (OR = 2.03, 95% CI 2.02-2.04, P < 0.001). However, MR analysis did not provide evidence for a genetic causal relationship between kidney stones and these cancers. Sensitivity analyses confirmed the stability and reliability of the MR results.

CONCLUSION

Kidney stones increased the risk of renal cell carcinoma, bladder cancer, and prostate cancer in the US population. MR analysis did not establish a genetic causal relationship between kidney stones and renal cell carcinoma, bladder cancer, and prostate cancer in the European population.

Bidirectional Two-Sample Mendelian Randomization Study Reveals Causal Associations Between Aging and Endometriosis

Background

Previous studies have suggested that aging may influence reproductive functions of female. Nonetheless, the causal relationship between aging and endometriosis has yet to be completely understood.

Objective

This study aims to determine whether aging had a causal association with the incidence of endometriosis.

Methods

We conducted bidirectional MR analyses to evaluate the causal relationship between aging biomarkers, particularly leukocyte telomere length (LTL), and endometriosis risk. Instrumental variables for LTL were derived from the UK Biobank GWAS, while endometriosis-associated variants were obtained from the FinnGen GWAS dataset. Subgroup analyses were performed to investigate the association between LTL and endometriosis subtypes. Additionally, validation was performed using independent GWAS meta-analysis datasets.

Results

Inverse variance-weighted (IVW) analysis revealed a significant association between longer LTL and an increased risk of endometriosis (OR-IVW = 1.276, 95% CI: 1.143 to 1.424, FDR-adjusted P = 7.00E-5), with consistent findings across multiple MR methods. Sensitivity analysis using an independent GWAS meta-analysis dataset did not confirm the LTL-endometriosis association (OR-IVW = 1.128, 95% CI: 0.140 to 9.115, P = 0.910). Bidirectional MR analysis found no causal relationship between endometriosis and LTL. Subgroup analyses indicated that longer LTL was significantly associated with endometriosis of the ovary (OR-IVW = 1.343, 95% CI: 1.143 to 1.577, P = 3.00E-4) and endometriosis of the rectovaginal septum and vagina (OR-IVW = 1.336, 95% CI: 1.064 to 1.676, P = 0.013), while no significant association was found with endometriosis of the pelvic peritoneum.

Conclusion

Our findings suggest that longer LTL may contribute to an increased risk of endometriosis, particularly in ovarian and rectovaginal subtypes. However, no causal effect of endometriosis on aging was observed. The lack of replication in independent datasets highlights the potential influence of population heterogeneity and dataset-specific factors, warranting further validation in diverse cohorts.

The Pathogenesis of Inflammatory Bowel Diseases

Inflammatory bowel diseases (IBDs) are relapsing, remitting inflammatory diseases of the intestinal tract. Familial aggregation and genome-wide association studies revealed susceptibility variants that point toward a combination of innate immune and adaptive immune dysregulation that in concert with environmental factors, such as our microbiome, can initiate and perpetuate inflammation. Innate immune perturbations include functional abnormalities in the intestinal barrier, endoplasmic reticulum stress, and abnormal recognition of microbes. Adaptive immune changes include dysregulation of cytokines, regulatory T cells, and leukocyte migration. IBD is linked with an abnormal wound-healing response leading to fibrosis. This article summarizes key pathogenic mechanisms in the pathogenesis of IBDs.

Uncovering the genetic architecture and evolutionary roots of androgenetic alopecia in African men

Androgenetic alopecia is a highly heritable trait. However, much of our understanding about the genetics of male-pattern baldness comes from individuals of European descent. Here, we examined a dataset comprising 2,136 men from Ghana, Nigeria, Senegal, and South Africa that were genotyped using the Men of African Descent and Carcinoma of the Prostate Array. We first tested how genetic predictions of baldness generalize from Europe to Africa and found that polygenic scores from European genome-wide association studies (GWASs) yielded area under the curve statistics that ranged from 0.513 to 0.546, indicating that genetic predictions of baldness generalized poorly from European to African populations. Subsequently, we conducted an African GWAS of androgenetic alopecia, focusing on self-reported baldness patterns at age 45. After correcting for age at recruitment, population structure, and study site, we identified 266 moderately significant associations, 51 of which were independent (p < 10-5, r2 < 0.2). Most baldness associations were autosomal, and the X chromosome does not seem to have a large impact on baldness in African men. Although Neanderthal alleles have previously been associated with skin and hair phenotypes, within the limits of statistical power, we did not find evidence that continental differences in the genetic architecture of baldness are due to Neanderthal introgression. While most loci that are associated with androgenetic alopecia do not have large integrative haplotype scores or fixation index statistics, multiple baldness-associated SNPs near the EDA2R and AR genes have large allele frequency differences between continents. Collectively, our findings illustrate how population genetic differences contribute to the limited portability of polygenic predictions across ancestries.

Association Between Hypoxia-Inducible Factor-1α and Neurological Diseases: A Bidirectional Two-Sample Mendelian Randomization Analysis

BACKGROUND

Previous studies have suggested that hypoxia-inducible factor 1-α (HIF-1α) exerted multiple effects on different central nervous system disorders. However, it is still uncertain whether plasma HIF-1α can be a causal indicator for the relevant diseases. This study aimed to test the causality relationship between plasma HIF-1α and neurological diseases, including cerebrovascular diseases, migraines, and neurodegenerative diseases with a Mendelian randomization (MR) method.

METHODS

Single-nucleotide polymorphisms (SNPs) genetically representing plasma HIF-1α were screened as instrumental variables (IVs). Summary-level data for neurological disorder from genome-wide association studies (GWAS) were identified as outcomes. The causal effects between the IVs and outcomes were determined via the major analysis of inverse-variance-weighted (IVW) method. The reverse causal direction was also performed to investigate the possibility of reverse causation.

RESULTS

The findings revealed that plasma HIF-1α was identified to be genetically associated with cardioembolic stroke (CES) (OR = 0.885; 95% confidence interval [CI] = 0.796-0.985, p = 0.026), migraine (OR = 0.941, 95% CI = 0.888-0.998, p = 0.041), and drug-induced migraine without aura (MOA) (OR = 0.586, 95% CI = 0.375-0.916, p = 0.019). There was no association identified in plasma HIF-1α with subarachnoid hemorrhage (SAH), other stroke and migraine subtype, and neurodegenerative disorders. The reverse-MR analysis revealed that the above-stated neurological diseases did not have a causal effect on plasma HIF-1α levels. Sensitivity and validation analyses support that the above results are stable.

CONCLUSIONS

Our research indicated that plasma HIF-1α may have a causal effect on the risk of CES, migraine and drug-induced MOA, providing new insights for those disease prevention and therapeutic approaches.

Surfactant protein levels and genetic variants as biomarkers for COVID-19 severity in children

Since its outbreak, the novel coronavirus (COVID-19) has significantly impacted the pediatric population. Pulmonary surfactant dysfunction has been linked to other respiratory diseases in children and COVID-19 in adults, but its role in COVID-19 severity remains unclear. We hypothesized that elevated surfactant protein (SP) levels and single nucleotide polymorphisms (SNPs) of SP genes are associated with severe COVID-19 in children. We enrolled 325 COVID-19 positive children and categorized them as having mild or severe disease. Plasma SP-A, SP-B, and SP-D levels were measured. DNA was extracted and genotyped for SNPs in five SP genes, SFTPA1, SFTPA2, SFTPB, SFTPC, and SFTPD. Quantile regression was used to compare SP levels between groups, and receiver operating curve analysis determined an optimal cutoff value of SP level for predicting severe COVID-19. Logistic regression evaluated the odds ratio (OR) for severe disease and associations between SNPs and COVID-19 severity. We found that increased plasma SP-A levels, but not SP-B or SP-D, were significantly associated with severe COVID-19. No significant correlation was observed between age and SP levels. A plasma SP-A level of 10 ng/mL was identified as the optimal cutoff for predicting severe COVID-19, with an OR of 5.9, indicating that children with SP-A levels above this threshold are nearly six times more likely to develop severe COVID-19 disease. In addition, the rs8192340 of SFTPC was associated with decreased risk of severe COVID-19 before, but not after, Bonferroni correction. These findings suggest that plasma SP-A may serve as a potential biomarker for severe COVID-19 in children.NEW & NOTEWORTHY Surfactant dysfunction is linked to other pulmonary diseases, but its role in pediatric coronavirus (COVID-19) is unclear. We found elevated plasma surfactant protein (SP)-A levels, but not SP-B or SP-D, significantly associated with severe COVID-19. A plasma SP-A threshold of 10 ng/mL predicted severe COVID-19. The rs8192340 of SFTPC was associated with decreased risk of severe COVID-19 before, but not after, Bonferroni correction. These findings suggest plasma SP-A may serve as a potential biomarker for pediatric COVID-19 severity.

Genetic architecture in Greenland is shaped by demography, structure and selection

Greenlandic Inuit and other indigenous populations are underrepresented in genetic research1,2, leading to inequity in healthcare opportunities. To address this, we performed analyses of sequenced or imputed genomes of 5,996 Greenlanders with extensive phenotypes. We quantified their historical population bottleneck and how it has shaped their genetic architecture to have fewer, but more common, variable sites. Consequently, we find twice as many high-impact genome-wide associations to metabolic traits in Greenland compared with Europe. We infer that the high-impact variants arose after the population split from Native Americans and thus are Arctic-specific, and show that some of them are common due to not only genetic drift but also selection. We also find that European-derived polygenic scores for metabolic traits are only half as accurate in Greenlanders as in Europeans, and that adding Arctic-specific variants improves the overall accuracy to the same level as in Europeans. Similarly, lack of representation in public genetic databases makes genetic clinical screening harder in Greenlandic Inuit, but inclusion of Greenlandic data remedies this by reducing the number of non-causal candidate variants by sixfold. Finally, we identify pronounced genetic fine structure that explains differences in prevalence of monogenic diseases in Greenland and, together with recent changes in mobility, leads to a predicted future reduction in risk for certain recessive diseases. These results illustrate how including data from Greenlanders can greatly reduce inequity in genomic-based healthcare.

High-dimensional, outcome-dependent missing data problems: Models for the human KIR loci

Missing data problems are common in biological, high-dimensional data, where data can be partially or completely missing. Algorithms have been developed to reconstruct the missing values by means of imputation or expectation-maximization algorithms. For missing data problems, it has been suggested that the regression model of interest should be incorporated into the imputation procedure to reduce bias of the regression coefficients. We here consider a challenging missing data problem, where diplotypes of the KIR loci are to be reconstructed. These loci are difficult to genotype, resulting in ambiguous genotype calls. We extend a previously proposed expectation-maximization algorithm by incorporating a potentially high-dimensional regression model to model the outcome. Three strategies are evaluated: (1) only allelic predictors, (2) allelic predictors and forward-backward selection on haplotype predictors, and (3) penalized regression on a saturated model. In a simulation study, we compared these strategies with a baseline expectation-maximization algorithm without outcome model. For extreme choices of effect sizes and missingness levels, the outcome-based expectation-maximization algorithms outperformed the no-outcome expectation-maximization algorithm. However, in all other cases, the no-outcome expectation-maximization algorithm performed either superior or comparable to the three strategies, suggesting the outcome model can have a harmful effect. In a data analysis concerning death after allogeneic hematopoietic stem cell transplantation as a function of donor KIR genes, expectation-maximization algorithms with and without outcome showed very similar results. In conclusion, outcome based missing data models in the high-dimensional setting have to be used with care and are likely to lead to biased results.

A Review of Mendelian Randomization: Assumptions, Methods, and Application to Obesity-Related Diseases

Mendelian randomization (MR) is a statistical method that uses genetic variants as instrumental variables to estimate the causal effect of exposure on an outcome in the presence of unmeasured confounding. In this review, we argue that it is crucial to acknowledge the instrumental variable assumptions in MR analysis. We describe widely used MR methods, using an example from obesity-related metabolic disorders. We describe situations in which instrumental variable assumptions are violated and explain how to evaluate these violations and employ robust methods for accommodating such violations.

The relationship between smoking and recurrent aphthous stomatitis: A Mendelian randomization study

INTRODUCTION

Existing research suggests an association between smoking and the incidence of recurrent aphthous stomatitis (RAS); however, the causal relationship remains ambiguous. We employed Mendelian randomization (MR) to clarify the potential causal association between smoking and the risk of developing RAS.

METHODS

We utilized genome-wide association study (GWAS) sequencing data related to smoking from the Finnish database as instrumental variables (IVs) and GWAS data for RAS from the UK Biobank (UKB) as the outcome to perform a two-sample MR analysis. The selection of IVs was rigorously controlled according to the three principal assumptions of relevance, independence, and exclusivity. The primary analytical methods utilized were inverse variance weighting (IVW) and weighted median (WM), supplemented by MR-Egger, simple mode, and weighted mode techniques to infer causality between smoking and RAS. Sensitivity analyses were conducted using MR-PRESSO, Cochran's Q, and the MR-Egger intercept to ensure the robustness of the findings.

RESULTS

The findings from the IVW and WM analyses suggest a causal association between smoking and an elevated risk of RAS (IVW: OR=1.003; 95% CI: 1.0002-1.005, p=0.033; WM: OR=1.003; 95% CI: 1.00006-1.007, p=0.044). Compared to non-smokers, smokers have a 0.3% increase in the risk of RAS. Furthermore, the sensitivity analysis did not reveal any inconsistencies that would contradict the MR results.

CONCLUSIONS

Our findings provide preliminary evidence of a potential causal relationship between smoking and the risk of RAS, which may contribute to a deeper understanding of the underlying mechanisms. Further research is needed to confirm these results and explore their implications for clinical practice.

The genomic architecture of circulating cytokine levels points to drug targets for immune-related diseases

Circulating cytokines orchestrate immune reactions and are promising drug targets for immune-mediated and inflammatory diseases. Exploring the genetic architecture of circulating cytokine levels could yield key insights into causal mediators of human disease. Here, we performed genome-wide association studies (GWAS) for 40 circulating cytokines in meta-analyses of 74,783 individuals. We detected 359 significant associations between cytokine levels and variants in 169 independent loci, including 150 trans- and 19 cis-acting loci. Integration with transcriptomic data point to key regulatory mechanisms, such as the buffering function of the Atypical Chemokine Receptor 1 (ACKR1) acting as scavenger for multiple chemokines and the role of tumor necrosis factor receptor-associated factor 1 (TRAFD1) in modulating the cytokine storm triggered by TNF signaling. Applying Mendelian randomization (MR), we detected a network of complex cytokine interconnections with TNF-b, VEGF, and IL-1ra exhibiting pleiotropic downstream effects on multiple cytokines. Drug target cis-MR using 2 independent proteomics datasets paired with colocalization revealed G-CSF/CSF-3 and CXCL9/MIG as potential causal mediators of asthma and Crohn's disease, respectively, but also a potentially protective role of TNF-b in multiple sclerosis. Our results provide an overview of the genetic architecture of circulating cytokines and could guide the development of targeted immunotherapies.

Association between bone mineral density and scoliosis: a two-sample mendelian randomization study in european populations

BACKGROUND

Previous studies have shown that bone mineral density (BMD) has a certain impact on scoliosis. However, up to now, there is no clear evidence that there is a causal association between the two. The aim of this study is to investigate whether there is a causal association between BMD at different body positions and scoliosis by two-sample Mendelian randomization (MR).

METHODS

Genetic variants (SNPS) strongly associated with BMD (total body BMD (TB-BMD), lumbar spine BMD (LS-BMD), femoral neck BMD (FN-BMD), heel BMD (HE-BMD), and forearm BMD (FA-BMD)) were extracted from GEFOS and genome-wide association analysis (GWAS) databases SNPs) were used as instrumental variables (IVs). Scoliosis was also selected from the Finnish database as the outcome. Inverse variance weighting (IVW) method was used as the main analysis method, and multiple sensitivity analysis was performed by combining weighted median, MR-Egger, MR Multi-effect residuals and outliers.

RESULTS

IVW results showed that TB-BMD (OR = 0.83, 95%CI: 0.66-1.55 P = 0.13), LS-BMD (OR = 0.72, 95%CI: 0.52-0.99, P = 0.04), FN-BMD (OR = 0.74, 95%CI: 0.50-1.09, P = 0.13), FA-BMD (OR = 0.95,95%CI: 0.70-1.28, P = 0.75), HE-BMD (OR = 0.91, 95%CI: 0.77-1.08, P = 0.29). Sensitivity analyses showed no evidence of pleiotropy or heterogeneity (p > 0.05) (MR-PRESSO and Cochrane). The results were further validated by leave-one-out test and MR-Egger intercept, which confirmed the robustness of the study results.

CONCLUSION

In conclusion, the present study demonstrates that the causal role of genetic prediction of scoliosis increases with decreasing lumbar BMD. There was no evidence that BMD at the remaining sites has a significant causal effect on scoliosis. Our results suggest that the lumbar spine BMD should be routinely measured in the population at high risk of scoliosis. If osteoporosis occurs, appropriate treatment should be given to reduce the incidence of scoliosis.

CLINICAL TRIAL NUMBER

Not applicable.

The emergence of nonlinear evolutionary trade-offs and the maintenance of genetic polymorphisms

Evolutionary models of quantitative traits often assume trade-offs between beneficial and detrimental traits, requiring modellers to specify a function linking trait values. The choice of trade-off function can be consequential; functions that assume diminishing returns (accelerating costs) typically lead to single equilibrium genotypes, while decelerating costs often lead to genetic polymorphisms. Despite their importance, our current theory has little to say on which trade-off functions are the most biologically plausible. To address this gap, we explored how the genetic determination of quantitative traits can lead to different trade-off functions, using resistance to infectious diseases as an example trait. We developed a model where alleles at separate loci pleiotropically increase resistance while decreasing fecundity. We then used this model to generate genotype landscapes and investigate how epistasis effects the trade-off function. Regardless of the strength of epistasis, our model consistently led to accelerating costs. We then incorporated our genotype model into an eco-evolutionary model of disease resistance. Unlike other models with accelerating costs, our approach often led to genetic polymorphisms. Our results suggest that accelerating costs are a strong null model for evolutionary trade-offs and that the eco-evolutionary conditions required for polymorphism may be more nuanced than previously thought.

Causal association between metabolites and upper gastrointestinal tumors: A Mendelian randomization study

Upper gastrointestinal (UGI) tumors, notably gastric cancer (GC) and esophageal cancer (EC), are significant global health concerns due to their high morbidity and mortality rates. However, only a limited number of metabolites have been identified as biomarkers for these cancers. To explore the association between metabolites and UGI tumors, the present study conducted a comprehensive two‑sample Mendelian randomization (MR) analysis using publicly available genetic data. In the present study, the causal relationships were examined between 1,400 metabolites and UGI cancer using methods such as inverse variance weighting and weighted medians, along with sensitivity analyses for heterogeneity and pleiotropy. Functional experiments were conducted to validate the MR results. The analysis identified 57 metabolites associated with EC and 58 with GC. Key metabolites included fructosyllysine [EC: Odds ratio (OR)=1.450, 95% confidence interval (CI)=1.087‑1.934, P=0.011; GC: OR=1.728, 95% CI=1.202‑2.483, P=0.003], 2'‑deoxyuridine to cytidine ratio (EC: OR=1.464, 95% CI=1.111‑1.929, P=0.007; GC: OR=1.464, 95% CI=1.094‑1.957, P=0.010) and carnitine to protonylcarnitine (C3) ratio (EC: OR=0.655, 95% CI=0.499‑0.861, P=0.002; GC: OR=0.664, 95% CI=0.486‑0.906, P=0.010). Notably, fructosyllysine levels and the 2'‑deoxyuridine to cytidine ratio were identified as risk factors for both EC and GC, while the C3 ratio served as a protective factor. Functional experiments demonstrated that fructosyllysine and the 2'‑deoxyuridine to cytidine ratio promoted the proliferation of EC and GC cells, whereas carnitine inhibited their proliferation. In conclusion, the present findings provide insights into the causal factors and biomarkers associated with UGI tumors, which may be instrumental in guiding targeted dietary and pharmacological interventions, thereby contributing to the prevention and treatment of UGI cancer.

Towards functional maps of non-coding variants in cancer

Large scale cancer genomic studies in patients have unveiled millions of non-coding variants. While a handful have been shown to drive cancer development, the vast majority have unknown function. This review describes the challenges of functionally annotating non-coding cancer variants and understanding how they contribute to cancer. We summarize recently developed high-throughput technologies to address these challenges. Finally, we outline future prospects for non-coding cancer genetics to help catalyze personalized cancer therapy.

Unifying approaches from statistical genetics and phylogenetics for mapping phenotypes in structured populations

In both statistical genetics and phylogenetics, a major goal is to identify correlations between genetic loci or other aspects of the phenotype or environment and a focal trait. In these 2 fields, there are sophisticated but disparate statistical traditions aimed at these tasks. The disconnect between their respective approaches is becoming untenable as questions in medicine, conservation biology, and evolutionary biology increasingly rely on integrating data from within and among species, and once-clear conceptual divisions are becoming increasingly blurred. To help bridge this divide, we lay out a general model describing the covariance between the genetic contributions to the quantitative phenotypes of different individuals. Taking this approach shows that standard models in both statistical genetics (e.g., genome-wide association studies; GWAS) and phylogenetic comparative biology (e.g., phylogenetic regression) can be interpreted as special cases of this more general quantitative-genetic model. The fact that these models share the same core architecture means that we can build a unified understanding of the strengths and limitations of different methods for controlling for genetic structure when testing for associations. We develop intuition for why and when spurious correlations may occur analytically and conduct population-genetic and phylogenetic simulations of quantitative traits. The structural similarity of problems in statistical genetics and phylogenetics enables us to take methodological advances from one field and apply them in the other. We demonstrate by showing how a standard GWAS technique-including both the genetic relatedness matrix (GRM) as well as its leading eigenvectors, corresponding to the principal components of the genotype matrix, in a regression model-can mitigate spurious correlations in phylogenetic analyses. As a case study, we re-examine an analysis testing for coevolution of expression levels between genes across a fungal phylogeny and show that including eigenvectors of the covariance matrix as covariates decreases the false positive rate while simultaneously increasing the true positive rate. More generally, this work provides a foundation for more integrative approaches for understanding the genetic architecture of phenotypes and how evolutionary processes shape it.

Heterogeneous genetic architectures of prostate cancer susceptibility in sub-Saharan Africa

Men of African descent have the highest prostate cancer incidence and mortality rates, yet the genetic basis of prostate cancer in African men has been understudied. We used genomic data from 3,963 cases and 3,509 controls from Ghana, Nigeria, Senegal, South Africa and Uganda to infer ancestry-specific genetic architectures and fine-map disease associations. Fifteen independent associations at 8q24.21, 6q22.1 and 11q13.3 reached genome-wide significance, including four new associations. Intriguingly, multiple lead associations are private alleles, a pattern arising from recent mutations and the out-of-Africa bottleneck. These African-specific alleles contribute to haplotypes with odds ratios above 2.4. We found that the genetic architecture of prostate cancer differs across Africa, with effect size differences contributing more to this heterogeneity than allele frequency differences. Population genetic analyses reveal that African prostate cancer associations are largely governed by neutral evolution. Collectively, our findings emphasize the utility of conducting genetic studies that use diverse populations.

First Application of Whole Genome Sequencing in Myelinated Retinal Nerve Fibers (MRNF)

Genetic features are currently unknown in myelinated retinal nerve fibers (MRNF). For a 20-year-old asymptomatic female with unilateral MRNF, we performed whole genome sequencing (WGS) by standard workflow protocol to produce contiguous long-read sequences with Illumina DNA PCR-Free Prep. After tagmentation, libraries were sequenced on separate runs via NovaSeq 6000 platform at 2 x 150bp read length. Gene variants included rs2248799, rs2672589, rs7555070, rs247616_T and rs2043085_C all associated with an increased macular degeneration risk, and seven novel variants of uncertain significance. For optic disc enlargement, variants rs9988687_A, rs11079419_T, rs6787363 and rs10862708_A suggested an increased risk for this condition. In contrast, modeling revealed retinal detachment risk was reduced by variants identified at rs9651980_T, rs4373767_T, and rs7940691_T which were among five other previously unreported variants. WGS data placed proband at the 66th and 64th percentiles for disc anomaly and retinal detachment risk, respectively. Additionally, risk determined from 16 loci associated with age-related macular degeneration found the patient to be at the 18th percentile for this diagnosis (i.e., below average genetic predisposition). Fundoscopic findings showed mean RNFL thickness was lower with MRNF (77 OS vs. 96?m OD) and RNFL symmetry was impaired (43 %) but stable between 2020 and 2023. Rim area and cup volume were also substantially different (2.33 OS vs. 1.34mm2 OD, and 0.001 OS vs. 0.151mm3 OD, respectively). As the first known evaluation of MRNF via WGS, these data reveal a mixed picture with variants associated with different risks for potentially related ocular pathologies. In addition, we identify multiple new variants of unknown significance. Factors affecting gene expression in MRNF require further study. Key words: Whole genome sequencing, Retina, Myelination, Anatomy, Gene variants.

Coeliac disease and postpartum depression: are they linked? A two-sample Mendelian randomization study

Background

To explore the potential causal associations between coeliac disease(CD) and postpartum depression(PPD) by using two-sample Mendelian randomization(MR) analysis.

Methods

The IEU OPEN GWAS project was utilized to identify genetic loci strongly associated with CD as instrumental variables (IVs), and MR analysis was performed using inverse variance weighting(IVW), weighted median, weighted model, and MR-Egger. MR analyses were used to examine whether there was a link between CD and PPD, with an OR and 95% CI. Meanwhile, the relationship between CD and depression(DP) was analyzed using MR. The sensitivity analysis was conducted using MR-Egger intercept analysis, Cochran's Q test, and leave-one-out analysis.

Results

From the GWAS online database, 13 single-nucleotide polymorphisms (SNPs) were chosen as IVs. The IVW results showed a relationship between PPD and a genetically predicted risk of developing CD (OR = 1.022, 95% CI: 1.001-1.044, P = 0.043). However, the presence of DP was not linked with CD (OR=0.991, 95% CI: 0.978-1.003, P=0.151). Potential horizontal pleiotropy was not discovered using MR-Egger intercept analysis (PPD: P=0.725; DP: P=0.785), and Cochran's Q test for heterogeneity revealed no significant heterogeneity (PPD: P=0.486; DP: P=0.909). A leave-one-out analysis found that individual SNPs had minimal effect on overall causal estimations.

Conclusion

MR research discovered a link between CD and PPD.

An externally validated resting-state brain connectivity signature of pain-related learning

Pain can be conceptualized as a precision signal for reinforcement learning in the brain and alterations in these processes are a hallmark of chronic pain conditions. Investigating individual differences in pain-related learning therefore holds important clinical and translational relevance. Here, we developed and externally validated a novel resting-state brain connectivity-based predictive model of pain-related learning. The pre-registered external validation indicates that the proposed model explains 8-12% of the inter-individual variance in pain-related learning. Model predictions are driven by connections of the amygdala, posterior insula, sensorimotor, frontoparietal, and cerebellar regions, outlining a network commonly described in aversive learning and pain. We propose the resulting model as a robust and highly accessible biomarker candidate for clinical and translational pain research, with promising implications for personalized treatment approaches and with a high potential to advance our understanding of the neural mechanisms of pain-related learning.

Comparative analysis of genetic risk scores for predicting biochemical recurrence in prostate cancer patients after radical prostatectomy

BACKGROUND

In recent years, Genome-Wide Association Studies (GWAS) has identified risk variants related to complex diseases, but most genetic variants have less impact on phenotypes. To solve the above problems, methods that can use variants with low genetic effects, such as genetic risk score (GRS), have been developed to predict disease risk.

METHODS

As the GRS model with the most incredible prediction power for complex diseases has not been determined, our study used simulation data and prostate cancer data to explore the disease prediction power of three GRS models, including the simple count genetic risk score (SC-GRS), the direct logistic regression genetic risk score (DL-GRS), and the explained variance weighted GRS based on directed logistic regression (EVDL-GRS).

RESULTS AND CONCLUSIONS

We used 26 SNPs to establish GRS models to predict the risk of biochemical recurrence (BCR) after radical prostatectomy. Combining clinical variables such as age at diagnosis, body mass index, prostate-specific antigen, Gleason score, pathologic T stage, and surgical margin and GRS models has better predictive power for BCR. The results of simulation data (statistical power = 0.707) and prostate cancer data (area under curve = 0.8462) show that DL-GRS has the best prediction performance. The rs455192 was the most relevant locus for BCR (p = 2.496 × 10-6) in our study.

The Emergence of Non-Linear Evolutionary Trade-offs and the Maintenance of Genetic Polymorphisms

Evolutionary models of quantitative traits often assume trade-offs between beneficial and detrimental traits, requiring modelers to specify a function linking costs to benefits. The choice of trade-off function is often consequential; functions that assume diminishing returns (accelerating costs) typically lead to single equilibrium genotypes, while decelerating costs often lead to evolutionary branching. Despite their importance, we still lack a strong theoretical foundation to base the choice of trade-off function. To address this gap, we explore how trade-off functions can emerge from the genetic architecture of a quantitative trait. We developed a multi-locus model of disease resistance, assuming each locus had random antagonistic pleiotropic effects on resistance and fecundity. We used this model to generate genotype landscapes and explored how additive versus epistatic genetic architectures influenced the shape of the trade-off function. Regardless of epistasis, our model consistently led to accelerating costs. We then used our genotype landscapes to build an evolutionary model of disease resistance. Unlike other models with accelerating costs, our approach often led to genetic polymorphisms at equilibrium. Our results suggest that accelerating costs are a strong null model for evolutionary trade-offs and that the eco-evolutionary conditions required for polymorphism may be more nuanced than previously believed.

Association between n-3 PUFA and lung function: results from the NHANES 2007-2012 and Mendelian randomisation study

This study aimed to investigate the association between n-3 PUFA and lung function. First, a cross-sectional study was conducted based on the National Health and Nutrition Examination Survey (NHANES) 2007-2012 data. n-3 PUFA intake was obtained from 24-h dietary recalls. A multivariable linear regression model was used to assess the observational associations of n-3 PUFA intake with lung function. Subsequently, a two-sample Mendelian randomisation (MR) was performed to estimate the potential causal effect of n-3 PUFA on lung function. Genetic instrumental variables were extracted from published genome-wide association studies. Summary statistics about n-3 PUFA was from UK Biobank. Inverse variance weighted was the primary analysis approach. The observational study did not demonstrate a significant association between n-3 PUFA intake and most lung function measures; however, a notable exception was observed with significant findings in the highest quartile for forced vital capacity (FVC) and % predicted FVC. The MR results also showed no causal effect of circulating n-3 PUFA concentration on lung function (forced expiratory volume in one second (FEV1), β = 0·01301, se = 0·01932, P = 0·5006; FVC, β = -0·001894, se = 0·01704, P = 0·9115; FEV1:FVC, β = 0·03118, se = 0·01743, P = 0·07359). These findings indicate the need for further investigation into the impact of higher n-3 PUFA consumption on lung health.

Metabolome-wide Mendelian randomization for age at menarche and age at natural menopause

BACKGROUND

The role of metabolism in the variation of age at menarche (AAM) and age at natural menopause (ANM) in the female population is not entirely known. We aimed to investigate the causal role of circulating metabolites in AAM and ANM using Mendelian randomization (MR).

METHODS

We combined MR with genetic colocalization to investigate potential causal associations between 658 metabolites and AAM and between 684 metabolites and ANM. We extracted genetic instruments for our exposures from four genome-wide association studies (GWAS) on circulating metabolites and queried the effects of these variants on the outcomes in two large GWAS from the ReproGen consortium. Additionally, we assessed the mediating role of the body mass index (BMI) in these associations, identified metabolic pathways implicated in AAM and ANM, and sought validation for selected metabolites in the Avon Longitudinal Study of Parents and Children (ALSPAC).

RESULTS

Our analysis identified 10 candidate metabolites for AAM, but none of them colocalized with AAM. For ANM, 76 metabolites were prioritized (FDR-adjusted MR P-value ≤ 0.05), with 17 colocalizing, primarily in the glycerophosphocholines class, including the omega-3 fatty acid and phosphatidylcholine (PC) categories. Pathway analyses and validation in ALSPAC mothers also highlighted the role of omega and polyunsaturated fatty acids levels in delaying age at menopause.

CONCLUSIONS

Our study suggests that metabolites from the glycerophosphocholine and fatty acid families play a causal role in the timing of both menarche and menopause. This underscores the significance of specific metabolic pathways in the biology of female reproductive longevity.

The causal role of multiple psycho-emotional disorders in gastroesophageal reflux disease: A two-sample Mendelian randomized study

BACKGROUND

Observational studies have previously shown a potential link between psycho-emotional disorders, such as mood swings, highly strung, anxious feelings, and gastroesophageal reflux disease (GERD). However, the credibility of these associations could be influenced by various confounding factors. Consequently, our study sought to employ a Mendelian randomization (MR) approach to elucidate a potential causal relationship between psycho-emotional disorders and GERD.

METHOD

Information on independent genetic variants linked to mood swings, highly strung, and anxious feelings was gathered from European populations participating in the IEU Open GWAS research. The FinnGen Consortium provided the genome-wide association study (GWAS) summary statistics for GERD. Our analysis employed the inverse variance weighted (IVW) method under the random effects model as the main analytical method. To further bolster our findings, we employed the weighted median and MR Egger methods. In addition, we conducted a series of sensitivity analyses.

RESULTS

Our study supports the existence of a causal relationship between psycho-emotional disorders and GERD. Mood swings, highly strung, and anxious feelings adversely affected GERD risk (mood swings: OR 2.21, 95% CI 1.19-5.59, p = 3.09 × 10-2; highly strung: OR 5.63, 95% CI 1.77-17.94, p = 3.42 × 10-3; anxious feelings: OR 2.48, 95% CI 1.08-4.33, p = 2.89 × 10-2).

CONCLUSION

This Mendelian randomization study provides robust support for the notion that mood swings, highly strung and anxious feelings, are associated with an increased risk of developing GERD.

Association between migraine and venous thromboembolism: a Mendelian randomization and genetic correlation study

Objective

Previous observational studies have reported an increased risk of venous thromboembolism (VTE) among individuals with migraine. This study aimed to investigate the causal effect of migraine on the development of VTE, as well as explore the genetic correlation between them.

Methods

We conducted a two-sample Mendelian randomization (MR) study using publicly available summary statistics from large-scale genome-wide association studies for migraine and VTE. Linkage disequilibrium score regression analysis was performed to estimate the genetic correlation between migraine and VTE.

Results

There were several shared risk variants (p-value < 5 × 10-8) between migraine and VTE. Linkage disequilibrium score regression analysis found a significant positive genetic correlation between migraine and VTE. The genetic correlations based on two migraine datasets were 0.208 (se = 0.031, p-value = 2.91 × 10-11) and 0.264 (se = 0.040, p-value = 4.82 × 10-11), respectively. Although main MR analysis showed that migraine was associated with an increased risk of VTE (odds ratio = 1.069, 95% confidence interval = 1.022-1.118, p-value = 0.004), the association attenuated to non-significance when using several other MR methods and using another set of genetic instruments. In addition, evidence of heterogeneity was found. Reverse MR analysis showed VTE was associated with increased risk of migraine with aura (odds ratio = 1.137, 95% confidence interval = 1.062-1.218, p-value = 2.47 × 10-4) with no evidence of pleiotropy and heterogeneity.

Conclusion

We showed suggestive evidence indicating an association between migraine and increased risk of VTE. Additionally, we found robust evidence suggesting that VTE is associated with an increased risk of migraine. The positive genetic correlation indicates that migraine and VTE has shared genetic basis. Further investigations will be necessary to address potential sex-specific effects in the analysis.

Circulating cytokines and vascular dementia: A bi-directional Mendelian randomization study

Inflammatory responses are associated with the development of vascular dementia (VaD). Circulating cytokines modulate the inflammatory response and are important for the immune system. To further elucidate the role of the immune system in VaD, we used Mendelian randomization (MR) to comprehensively and bi-directionally assess the role of circulating cytokines in VaD. Using state-of-the-art genome-wide association studies, we primarily assessed whether different genetic levels of 41 circulating cytokines affect the risk of developing VaD and, in turn, whether the genetic risk of VaD affects these circulating cytokines. We used inverse variance weighting (IVW) and several other MR methods to assess the bidirectional causality between circulating cytokines and VaD, and performed sensitivity analyses. Tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) was inversely associated with VaD risk [odds ratio (OR): 0.74, 95 % confidence interval (CI): 0.60-0.92, P = 0.007, 0.007]. VaD was associated with seven circulating cytokines: macrophage inflammatory protein 1b (MIP-1 beta) [OR: 1.05, 95 % CI: 1.01-1.08, P = 0.009], Interleukin-12p70 (IL-12) [OR: 1.04, 95 % CI: 1.00-1.08, P = 0.047], Interleukin-17 (IL-17) [OR: 1.04, 95 % CI: 1.00-1.07, P = 0.038], Interleukin-7 (IL-7) [OR: 1.07, 95 % CI: 1.02-1.12, P = 0.009], Interferon gamma (IFN-γ) [OR: 1.03, 95 % CI: 1.00-1.07, P = 0.046], Granulocyte-colony stimulating factor (GCSF) [OR: 1.06, 95 % CI: 1.02-1.09, P = 0.001], Fibroblast growth factor (FGF) [P = 0.001], and Fibroblast growth factor (FGF) [P = 0.001]. Fibroblast growth factor basic (FGF-Basic) [OR: 1.04, 95 % CI: 1.01-1.08, P = 0.02] were positively correlated. Circulating cytokines are associated with VaD, and further studies are needed to determine whether they are effective targets for intervention to prevent or treat VaD.

Serum electrolyte concentrations and risk of atrial fibrillation: an observational and mendelian randomization study

BACKGROUND

Atrial fibrillation (AF) is a prevalent arrhythmic condition resulting in increased stroke risk and is associated with high mortality. Electrolyte imbalance can increase the risk of AF, where the relationship between AF and serum electrolytes remains unclear.

METHODS

A total of 15,792 individuals were included in the observational study, with incident AF ascertainment in the Atherosclerosis Risk in Communities (ARIC) study. The Cox regression models were applied to calculate the hazard ratio (HR) and 95% confidence interval (CI) for AF based on different serum electrolyte levels. Mendelian randomization (MR) analyses were performed to examine the causal association.

RESULTS

In observational study, after a median 19.7 years of follow-up, a total of 2551 developed AF. After full adjustment, participants with serum potassium below the 5th percentile had a higher risk of AF relative to participants in the middle quintile. Serum magnesium was also inversely associated with the risk of AF. An increased incidence of AF was identified in individuals with higher serum phosphate percentiles. Serum calcium levels were not related to AF risk. Moreover, MR analysis indicated that genetically predicted serum electrolyte levels were not causally associated with AF risk. The odds ratio for AF were 0.999 for potassium, 1.044 for magnesium, 0.728 for phosphate, and 0.979 for calcium, respectively.

CONCLUSIONS

Serum electrolyte disorders such as hypokalemia, hypomagnesemia and hyperphosphatemia were associated with an increased risk of AF and may also serve to be prognostic factors. However, the present study did not support serum electrolytes as causal mediators for AF development.

Causal relationship between telomere length and risk of intracranial aneurysm: a bidirectional Mendelian randomization study

Background

Telomere length is closely linked to the aging phenotype, where cellular aging results in the production of a cascade of cell factors and the senescence-associated secretory phenotype (SASP), leading to an inflammatory response. The presence of inflammation plays a crucial role in the formation of intracranial aneurysms. Nevertheless, the relationship between telomere length and intracranial aneurysms remains unclear. This study aims to explore the causal connection between telomere length and intracranial aneurysms through the utilization of Mendelian randomization (MR) analysis.

Methods

Data on telomere length were obtained from the genome-wide association studies conducted on the UK Biobank, comprising a total of 472,174 participants. Data on intracranial aneurysms were obtained from the summary dataset of the Global Genome-wide Association Study (GWAS) conducted by the International Stroke Genetics Consortium. The dataset consisted of 7,495 cases and 71,934 controls, all of European descent. Initially, the linkage disequilibrium score was used to investigate the connection between telomere length and intracranial aneurysms. Subsequently, a bidirectional MR was conducted using two-sample analysis to assess whether there is a causal connection between telomere length and intracranial aneurysm risk. The results were analyzed utilizing five MR methods, with the inverse variance weighted method serving as the main methodology. In addition, we did various analyses to evaluate the presence of heterogeneity, pleiotropy, and sensitivity in the study results. A reverse MR analysis was conducted to investigate potential reverse causal links.

Results

In the forward MR analysis, it was observed that both the inverse variance-weighted and weighted median analyses implied a potential causal relationship between longer telomere length and a decreased incidence of intracranial aneurysms (IVW: OR = 0.66, 95% CI: 0.47-0.92, p = 1.49 × 10-2). There was no heterogeneity or horizontal pleiotropy. The findings were verified to be robust through the utilization of leave-one-out analysis. The use of reverse MR analysis did not establish a potential causal link between the occurrence of intracranial aneurysms and telomere length.

Conclusion

There may exist a potential correlation between longer telomere length and a decreased likelihood of intracranial aneurysms within the European population. The present study offers novel insights into the correlation between telomere length and intracranial aneurysms. Additional research is required to clarify the underlying mechanisms and validate our discoveries in diverse populations.

Unifying approaches from statistical genetics and phylogenetics for mapping phenotypes in structured populations

In both statistical genetics and phylogenetics, a major goal is to identify correlations between genetic loci or other aspects of the phenotype or environment and a focal trait. In these two fields, there are sophisticated but disparate statistical traditions aimed at these tasks. The disconnect between their respective approaches is becoming untenable as questions in medicine, conservation biology, and evolutionary biology increasingly rely on integrating data from within and among species, and once-clear conceptual divisions are becoming increasingly blurred. To help bridge this divide, we derive a general model describing the covariance between the genetic contributions to the quantitative phenotypes of different individuals. Taking this approach shows that standard models in both statistical genetics (e.g., Genome-Wide Association Studies; GWAS) and phylogenetic comparative biology (e.g., phylogenetic regression) can be interpreted as special cases of this more general quantitative-genetic model. The fact that these models share the same core architecture means that we can build a unified understanding of the strengths and limitations of different methods for controlling for genetic structure when testing for associations. We develop intuition for why and when spurious correlations may occur using analytical theory and conduct population-genetic and phylogenetic simulations of quantitative traits. The structural similarity of problems in statistical genetics and phylogenetics enables us to take methodological advances from one field and apply them in the other. We demonstrate this by showing how a standard GWAS technique-including both the genetic relatedness matrix (GRM) as well as its leading eigenvectors, corresponding to the principal components of the genotype matrix, in a regression model-can mitigate spurious correlations in phylogenetic analyses. As a case study of this, we re-examine an analysis testing for co-evolution of expression levels between genes across a fungal phylogeny, and show that including covariance matrix eigenvectors as covariates decreases the false positive rate while simultaneously increasing the true positive rate. More generally, this work provides a foundation for more integrative approaches for understanding the genetic architecture of phenotypes and how evolutionary processes shape it.

Human genetic determinants of COVID-19 in Brazil: challenges and future plans

COVID-19 pandemic represented a worldwide major challenge in different areas, and efforts undertaken by the scientific community led to the understanding of some of the genetic determinants that influence the different COVID-19 outcomes. In this paper, we review the studies about the role of human genetics in COVID-19 severity and how Brazilian studies also contributed to those findings. Rare variants in genes related to Inborn Errors of Immunity (IEI) in the type I interferons pathway, and its phenocopies, have been described as being causative of severe outcomes. IEI and its phenocopies are present in Brazil, not only in COVID-19 patients, but also in autoimmune conditions and severe reactions to yellow fever vaccine. In addition, studies focusing on common variants and GWAS studies encompassing worldwide patients have found several loci associated with COVID-19 severity. A GWAS study including only Brazilian COVID-19 patients identified a new locus 1q32.1 associated with COVID-19 severity. Thus, more comprehensive studies considering the Brazilian genomic diversity should be performed, since they can help to reveal not only what are the genetic determinants that contribute to the different outcomes for COVID-19 in the Brazilian population, but in the understanding of human genetics in different health conditions.

Investigating Causal Effects of Hematologic Traits on Lung Cancer: A Mendelian Randomization Study

BACKGROUND

Observational studies have suggested blood cell counts may act as predictors of cancer. It is not known whether these hematologic traits are causally associated with lung cancer.

METHODS

Two-sample bidirectional univariable Mendelian randomization (MR) and multivariable MR (MVMR) were performed to investigate the causal association between hematologic traits and the overall risk of lung cancer and three histologic subtypes [lung adenocarcinoma, squamous cell lung cancer, and small cell lung cancer (SCLC)]. The instrumental variables of 23 hematologic traits were strictly selected from large-scale genome-wide association studies. Inverse-variance weighted method and five extra methods were used to obtain robust causal estimates.

RESULTS

We found evidence that genetically influenced higher hematocrit [OR, 0.845; 95% confidence interval (CI), 0.783-0.913; P = 1.68 × 10-5] and hemoglobin concentration (OR, 0.868; 95% CI, 0.804-0.938; P = 3.20 × 10-4) and reticulocyte count (OR, 0.923; 95% CI, 0.872-0.976; P = 5.19 × 10-3) decreased lung carcinoma risk, especially in ever smokers. MVMR further identified hematocrit independently of smoking as an independent predictor. Subgroup analysis showed that a higher plateletcrit level increased the risk of small cell lung carcinoma (OR, 1.288; 95% CI, 1.126-1.474; P = 2.25 × 10-4).

CONCLUSIONS

Genetically driven higher levels of reticulocyte count and hematocrit decreased lung cancer risk. Higher plateletcrit had an adverse effect on SCLC. Hematologic traits may act as low-cost factors for lung cancer risk stratification.

IMPACT

Further studies are required to elucidate the potential mechanisms underlying the dysregulation of homeostasis related to hematologic traits, such as subclinical inflammation.

Mendelian randomization identifies circulating proteins as biomarkers for age at menarche and age at natural menopause

Age at menarche (AAM) and age at natural menopause (ANM) are highly heritable traits and have been linked to various health outcomes. We aimed to identify circulating proteins associated with altered ANM and AAM using an unbiased two-sample Mendelian randomization (MR) and colocalization approach. By testing causal effects of 1,271 proteins on AAM, we identified 22 proteins causally associated with AAM in MR, among which 13 proteins (GCKR, FOXO3, SEMA3G, PATE4, AZGP1, NEGR1, LHB, DLK1, ANXA2, YWHAB, DNAJB12, RMDN1 and HPGDS) colocalized. Among 1,349 proteins tested for causal association with ANM using MR, we identified 19 causal proteins among which 7 proteins (CPNE1, TYMP, DNER, ADAMTS13, LCT, ARL and PLXNA1) colocalized. Follow-up pathway and gene enrichment analyses demonstrated links between AAM-related proteins and obesity and diabetes, and between AAM and ANM-related proteins and various types of cancer. In conclusion, we identified proteomic signatures of reproductive ageing in women, highlighting biological processes at both ends of the reproductive lifespan.

PharmGWAS: a GWAS-based knowledgebase for drug repurposing

Leveraging genetics insights to promote drug repurposing has become a promising and active strategy in pharmacology. Indeed, among the 50 drugs approved by FDA in 2021, two-thirds have genetically supported evidence. In this regard, the increasing amount of widely available genome-wide association studies (GWAS) datasets have provided substantial opportunities for drug repurposing based on genetics discoveries. Here, we developed PharmGWAS, a comprehensive knowledgebase designed to identify candidate drugs through the integration of GWAS data. PharmGWAS focuses on novel connections between diseases and small-molecule compounds derived using a reverse relationship between the genetically-regulated expression signature and the drug-induced signature. Specifically, we collected and processed 1929 GWAS datasets across a diverse spectrum of diseases and 724 485 perturbation signatures pertaining to a substantial 33609 molecular compounds. To obtain reliable and robust predictions for the reverse connections, we implemented six distinct connectivity methods. In the current version, PharmGWAS deposits a total of 740 227 genetically-informed disease-drug pairs derived from drug-perturbation signatures, presenting a valuable and comprehensive catalog. Further equipped with its user-friendly web design, PharmGWAS is expected to greatly aid the discovery of novel drugs, the exploration of drug combination therapies and the identification of drug resistance or side effects. PharmGWAS is available at https://ngdc.cncb.ac.cn/pharmgwas.

Associations of the circulating levels of cytokines with risk of ankylosing spondylitis: a Mendelian randomization study

Background

Observational studies have shown that changes in circulating cytokine/growth factor levels occur throughout the initiation and progression of ankylosing spondylitis (AS), yet whether they are etiologic or downstream effects remains unclear. In this study, we performed a summarized-level bidirectional Mendelian randomization (MR) analysis to shed light on the causal relationship between the two.

Methods

Genetic instrumental-variables (IVs) associated with circulating cytokine/growth factor levels were derived from a genome-wide association study (GWAS) of 8,293 European individuals, whereas summary data for the AS were obtained from a FinnGen GWAS of 166,144 participants. We used the inverse-variance-weighted (IVW) method as the main analysis for causal inference. Furthermore, several sensitivity analyses (MR-Egger, weighted median, MR-PRESSO and Cochran's Q test) were utilized to examine the robustness of the results. Finally, reverse MR analysis was performed to assess reverse causality between AS and circulating cytokine/growth factor levels.

Results

After Bonferroni correction, circulating levels of Cutaneous T-cell attracting (CTACK) and Monocyte specific chemokine 3 (MCP-3) were positively associated with a higher risk of AS (odds ratio [OR]: 1.224, 95% confidence interval [95% Cl]: 1.022 ~ 1.468, P = 0.028; OR: 1.250, 95% Cl: 1.016 ~ 1.539, P = 0.035). In addition, elevated circulating levels of Basic fibroblast growth factor (FGF-basic), Granulocyte colony-stimulating factor (G-CSF) and MCP-3 was considered a consequence of AS disease (β = 0.023, P = 0.017; β = 0.017, P = 0.025; β = 0.053, P = 0.025). The results of the sensitivity analysis were generally consistent.

Conclusion

The present study supplies genetic evidence for the relationship between circulating cytokine levels and AS. Targeted interventions of specific cytokines may help to reduce the risk of AS initiation and progression.

Interrogating Causal Effects of Body Composition and Puberty-Related Risk Factors on Adolescent Idiopathic Scoliosis: A Two-Sample Mendelian Randomization Study

Adolescent idiopathic scoliosis (AIS) is the most common form of pediatric musculoskeletal disorder. Observational studies have pointed to several risk factors for AIS, but almost no evidence exists to support their causal association with AIS. Here, we applied Mendelian randomization (MR), known to limit bias from confounding and reverse causation, to investigate causal associations between body composition and puberty-related exposures and AIS risk in Europeans and Asians. For our two-sample MR studies, we used single nucleotide polymorphisms (SNPs) associated with body mass index (BMI), waist-hip ratio, lean mass, childhood obesity, bone mineral density (BMD), 25-hydroxyvitamin D (25OHD), age at menarche, and pubertal growth in large European genome-wide association studies (GWAS), and with adult osteoporosis risk and age of menarche in Biobank Japan. We extracted estimates of the aforementioned SNPs on AIS risk from the European or Asian subsets of the largest multiancestry AIS GWAS (N = 7956 cases/88,459 controls). The results of our inverse variance-weighted (IVW) MR estimates suggest no causal association between the aforementioned risk factors and risk of AIS. Pleiotropy-sensitive MR methods yielded similar results. However, restricting our analysis to European females with AIS, we observed a causal association between estimated BMD and the risk of AIS (IVW odds ratio for AIS = 0.1, 95% confidence interval 0.01 to 0.7, p = 0.02 per SD increase in estimated BMD), but this association was no longer significant after adjusting for BMI, body fat mass, and 25OHD and remained significant after adjusting for age at menarche in multivariable MR. In conclusion, we demonstrated a protective causal effect of BMD on AIS risk in females of European ancestry, but this effect was modified by BMI, body fat mass, and 25OHD levels. Future MR studies using larger AIS GWAS are needed to investigate small effects of the aforementioned exposures on AIS. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Causal associations between circulation β-carotene and cardiovascular disease: A Mendelian randomization study

The causal association between circulating β-carotene concentrations and cardiovascular disease (CVD) remains controversial. We conducted a Mendelian randomization study to explore the effects of β-carotene on various cardiovascular diseases, including myocardial infarction, atrial fibrillation, heart failure, and stroke. Three single nucleotide polymorphisms (SNPs) associated with the β-carotene levels were obtained by searching published data and used as instrumental variables. Genetic association estimates for 4 CVDs (including myocardial infarction, atrial fibrillation, heart failure, and stroke) in the primary analysis, blood pressure and serum lipids (high-density lipoprotein [HDL] cholesterol, LDL cholesterol, and triglycerides) in the secondary analysis were obtained from large-scale genome-wide association studies (GWASs). We applied inverse variance-weighted as the primary analysis method, and 3 others were used to verify as sensitivity analysis. Genetically predicted circulating β-carotene levels (natural log-transformed, µg/L) were positively associated with myocardial infarction (odds ratio [OR] 1.10, 95% confidence interval [CI] 1.02-1.18, P = .011) after Bonferroni correction. No evidence supported the causal effect of β-carotene on atrial fibrillation (OR 1.02, 95% CI 0.96-1.09, P = .464), heart failure (OR 1.07, 95% CI 0.97-1.19, P = .187), stroke (OR 1.03, 95% CI 0.93-1.15, P = .540), blood pressure (P > .372) and serum lipids (P > .239). Sensitivity analysis produced consistent results. This study provides evidence for a causal relationship between circulating β-carotene and myocardial infarction. These findings have important implications for understanding the role of β-carotene in CVD and may inform dietary recommendations and intervention strategies for preventing myocardial infarction.

Genetic Atrial Cardiomyopathies: Common Features, Specific Differences, and Broader Relevance to Understanding Atrial Cardiomyopathy

Atrial cardiomyopathy is a condition that causes electrical and contractile dysfunction of the atria, often along with structural and functional changes. Atrial cardiomyopathy most commonly occurs in conjunction with ventricular dysfunction, in which case it is difficult to discern the atrial features that are secondary to ventricular dysfunction from those that arise as a result of primary atrial abnormalities. Isolated atrial cardiomyopathy (atrial-selective cardiomyopathy [ASCM], with minimal or no ventricular function disturbance) is relatively uncommon and has most frequently been reported in association with deleterious rare genetic variants. The genes involved can affect proteins responsible for various biological functions, not necessarily limited to the heart but also involving extracardiac tissues. Atrial enlargement and atrial fibrillation are common complications of ASCM and are often the predominant clinical features. Despite progress in identifying disease-causing rare variants, an overarching understanding and approach to the molecular pathogenesis, phenotypic spectrum, and treatment of genetic ASCM is still lacking. In this review, we aim to analyze the literature relevant to genetic ASCM to understand the key features of this rather rare condition, as well as to identify distinct characteristics of ASCM and its arrhythmic complications that are related to specific genotypes. We outline the insights that have been gained using basic research models of genetic ASCM in vitro and in vivo and correlate these with patient outcomes. Finally, we provide suggestions for the future investigation of patients with genetic ASCM and improvements to basic scientific models and systems. Overall, a better understanding of the genetic underpinnings of ASCM will not only provide a better understanding of this condition but also promises to clarify our appreciation of the more commonly occurring forms of atrial cardiomyopathy associated with ventricular dysfunction.

A novel missense compound heterozygous variant in TLR1 gene is associated with susceptibility to rheumatoid arthritis - structural perspective and functional annotations

INTRODUCTION

Besides human leukocyte antigen (HLA-DRB1) locus, more than 100 loci across the genome have been identified and linked with the onset, expression and/or progression of rheumatoid arthritis (RA). However, there are still grey areas in our understanding of the key genetic contributors of the disease, particularly in familial cases.

METHODS

In the present study, we have performed the whole exome sequencing (WES) of RA patients from two consanguineous families of Pakistan in a quest to identify novel, high-impact, RA-susceptibility genetic variants.

RESULTS

Through stepwise filtering, around 17,000 variants (common in the affected members) were recognized, out of which 2651 were predicted to be deleterious. Of these, 196 had direct relevance to RA. When selected for homozygous recessive mode of inheritance, two novel pathogenic variants (c.1324T>C, p.Cys442→Arg442; c.2036T>C, p.Ile679→Thr679) in the TLR1 gene displayed the role of compound heterozygosity in modulating the phenotypic expression and penetrance of RA. The structural and functional consequences of the TLR1 missense single nucleotide mutations (Cys442→Arg442; Ile679→Thr679) were evaluated through molecular dynamic simulation (MDS) studies. Analysis showed domain's rigidification, conferring stability to mutant TLR1-TIR/TIRAP-TIR complex with concomitant increase in molecular interactions with pro-inflammatory cytokines, compared to the wild-type conformation. Gene co-expression network analysis highlighted interlinked partnering genes along with interleukin-6 production of TLR1 (corrected p-value 2.98e-4) and acetylcholine receptor activity of CHRNG (corrected p-value 6.12e-2) as highly enriched associated functions.

CONCLUSION

The results, validated through case-control study subjects, suggested that the variants identified through WES and confirmed through Sanger sequencing and MDS are the novel disease variants and are likely to confer RA-susceptibility, independently and/or in a family-specific context. Key Points • Exploration of population based/ethno-specific big data is imperative to identify novel causal variants of RA. • Two new deleterious missense mutations in mutational hotspot exon 4 of TLR1 gene have been identified in Pakistani RA patients. • MD simulation data provides evidence for domain's rigidification, conferring stability to mutant TLR1-TIR/TIRAP-TIR complex, with concomitant increase in production of pro-inflammatory cytokines, thus adding to the onset/erosive outcome of RA.

Circulating Copper and Liver Cancer

The association between circulating copper and the risk of liver cancer has been investigated by previous studies, while the findings were inconsistent. Thus, we aimed to evaluate the association between circulating copper and liver cancer by using meta-analysis and Mendelian randomization (MR). For meta-analysis, PubMed and Web of Science were searched to identify eligible studies published before April 4, 2022. Standardized mean difference (SMD) with 95% confidence interval (CI) in circulating copper level between liver cancer patients and controls were pooled. Furthermore, we selected genetic instruments for circulating copper from a genome-wide association study (GWAS) to conduct MR analysis. The summary statistics related to liver cancer were obtained from two large independent cohorts, UKBB and FinnGen, respectively. MR analysis was performed mainly by inverse-variance weighted (IVW) approach, followed by maximum-likelihood method as sensitivity analysis. In meta-analysis of eight studies, circulating copper was found to be higher in liver cancer patients (SMD: 1.65; 95% CI: 0.65 to 2.65) with high heterogeneity (I2 = 96.40%, P = 0.001). However, inconsistent findings were observed among subgroups with high evidence. In MR analysis, genetically predicted circulating copper was not significantly associated with the risk of liver cancer by IVW in UKBB (OR: 1.38; 95% CI: 0.72 to 2.65) and FinnGen (OR: 1.10; 95% CI: 0.69 to 1.73) separately, and the pooled results produced similar results (OR: 1.18, 95% CI: 0.81 to 1.72). Moreover, non-significant finding was confirmed by using maximum-likelihood method. There is no sufficient evidence to demonstrate that high levels of circulating copper increase the risks of liver cancer.

Germline predisposition to clonal hematopoiesis

We now recognize that with aging, hematopoietic stem and progenitor cells (HSPCs) acquire mutations that confer a fitness advantage and clonally expand in a process now termed clonal hematopoiesis (CH). Because CH predisposes to a variety of health problems, including cancers, cardiovascular diseases, and inflammatory conditions, there is intense interest in the inherited alleles associated with the development of CH. DNA variants near TERT, SMC4, KPNA4, IL12A, CD164, and ATM confer the strongest associations. In this review, we discuss our current state of knowledge regarding germline predisposition to CH.

A novel penalized inverse-variance weighted estimator for Mendelian randomization with applications to COVID-19 outcomes

Mendelian randomization utilizes genetic variants as instrumental variables (IVs) to estimate the causal effect of an exposure variable on an outcome of interest even in the presence of unmeasured confounders. However, the popular inverse-variance weighted (IVW) estimator could be biased in the presence of weak IVs, a common challenge in MR studies. In this article, we develop a novel penalized inverse-variance weighted (pIVW) estimator, which adjusts the original IVW estimator to account for the weak IV issue by using a penalization approach to prevent the denominator of the pIVW estimator from being close to zero. Moreover, we adjust the variance estimation of the pIVW estimator to account for the presence of balanced horizontal pleiotropy. We show that the recently proposed debiased IVW (dIVW) estimator is a special case of our proposed pIVW estimator. We further prove that the pIVW estimator has smaller bias and variance than the dIVW estimator under some regularity conditions. We also conduct extensive simulation studies to demonstrate the performance of the proposed pIVW estimator. Furthermore, we apply the pIVW estimator to estimate the causal effects of five obesity-related exposures on three coronavirus disease 2019 (COVID-19) outcomes. Notably, we find that hypertensive disease is associated with an increased risk of hospitalized COVID-19; and peripheral vascular disease and higher body mass index are associated with increased risks of COVID-19 infection, hospitalized COVID-19, and critically ill COVID-19.

The causal relationship between genetically determined telomere length and meningiomas risk

Background

Studies have shown that longer leukocyte telomere length (LTL) is significantly associated with increased risk of meningioma. However, there is limited evidence concerning the causal association of LTL with benign and malignant meningiomas or with the location of benign tumors.

Methods

We used three LTL datasets from different sources, designated by name and sample size as LTL-78592, LTL-9190, and LTL-472174. The linkage disequilibrium score (LDSC) was used to explore the association between LTL and meningioma. We utilized two-sample bidirectional Mendelian randomization (TSMR) to evaluate whether LTL is causally related to meningioma risk. We adjusted for confounders by conducting multivariable Mendelian randomization (MVMR).

Results

In the LTL-78592, longer LTL was significantly associated with increased risk of malignant [odds ratio (OR) = 5.14, p = 1.04 × 10-5], benign (OR = 4.81, p < 0.05), benign cerebral (OR = 5.36, p < 0.05), and benign unspecified meningioma (OR = 8.26, p < 0.05). The same results were obtained for the LTL-9190. In the LTL-472174, longer LTL was significantly associated with increased risk of malignant (OR = 4.94, p < 0.05), benign (OR = 3.14, p < 0.05), and benign cerebral meningioma (OR = 3.59, p < 0.05). Similar results were obtained in the MVMR. In contrast, only benign cerebral meningioma displayed a possible association with longer LTL (OR = 1.01, p < 0.05). No heterogeneity or horizontal pleiotropy was detected.

Conclusion

In brief, genetically predicted longer LTL may increase the risk of benign, malignant, and benign cerebral meningiomas, regardless of the LTL measure, in European populations.

Increasing serum iron levels and their role in the risk of infectious diseases: a Mendelian randomization approach

OBJECTIVES

Increased iron stores have been associated with elevated risks of different infectious diseases, suggesting that iron supplementation may increase the risk of infections. However, these associations may be biased by confounding or reverse causation. This is important, since up to 19% of the population takes iron supplementation. We used Mendelian randomization (MR) to bypass these biases and estimate the causal effect of iron on infections.

METHODS

As instrumental variables, we used genetic variants associated with iron biomarkers in two genome-wide association studies (GWASs) of European ancestry participants. For outcomes, we used GWAS results from the UK Biobank, FinnGen, the COVID-19 Host Genetics Initiative or 23andMe, for seven infection phenotypes: 'any infections', combined, COVID-19 hospitalization, candidiasis, pneumonia, sepsis, skin and soft tissue infection (SSTI) and urinary tract infection (UTI).

RESULTS

Most of our analyses showed increasing iron (measured by its biomarkers) was associated with only modest changes in the odds of infectious outcomes, with all 95% odds ratios confidence intervals within the 0.88 to 1.26 range. However, for the three predominantly bacterial infections (sepsis, SSTI, UTI), at least one analysis showed a nominally elevated risk with increased iron stores (P <0.05).

CONCLUSION

Using MR, we did not observe an increase in risk of most infectious diseases with increases in iron stores. However for bacterial infections, higher iron stores may increase odds of infections. Hence, using genetic variation in iron pathways as a proxy for iron supplementation, iron supplements are likely safe on a population level, but we should continue the current practice of conservative iron supplementation during bacterial infections or in those at high risk of developing them.

Causal association between serum total bilirubin and cholelithiasis: a bidirectional two-sample Mendelian randomization study

Background

Observational studies about the association between serum total bilirubin and cholelithiasis are inconsistent. Hence, it is essential to reevaluate the association between serum total bilirubin and cholelithiasis and to verify whether such association is causal or not.

Methods

We selected single-nucleotide polymorphisms (SNPs) that are strongly associated with exposure as instrumental variable and conducted a bidirectional two-sample Mendelian randomization (MR) study to explore the causal association between serum total bilirubin and cholelithiasis. We implemented the inverse-variance weighted approach as a primary analysis to combine the Wald ratio estimates. Four additional analyses, namely, MR-Egger regression, weighted median, weighted mode, and MR-pleiotropy residual sum and outlier (PRESSO), were utilized to investigate the causal association and the influence of potential pleiotropy.

Results

A total of 116 SNPs were selected as valid instrumental variables to estimate the causal association of serum total bilirubin on cholelithiasis, and causal association between genetically determined serum total bilirubin and cholelithiasis was demonstrated [beta = 0.10; 95% confident interval (CI), 0.07 to 0.14; p < 0.001]. Likewise, the other methods, namely, the weighted median (beta = 0.12; 95% CI, 0.08 to 0.15; p < 0.001), MR-Egger (beta = 0.11; 95% CI, 0.08 to 0.15; p < 0.001), weighted mode (beta = 0.11; 95% CI, 0.08 to 0.15; p < 0.001), and MR-PRESSO approaches, further confirmed that this result (p = 0.054) indicates similar results. In addition, seven SNPs were selected as instrumental variable to estimate causal association of cholelithiasis on serum total bilirubin, and the result supported the causal effect of cholelithiasis to serum total bilirubin (beta = 0.12; 95% CI, 0.09 to 0.15; p < 0.001). At the same time, the other methods, namely, the weighted median (beta = 0.10; 95% CI, 0.06 to 0.13; p < 0.001), MR-Egger (beta = 0.12; 95% CI, 0.07 to 0.18; p = 0.007), weighted mode (beta = 0.09; 95% CI, 0.03 to 0.14, p = 0.019), and MR-PRESSO methods, further confirmed this result (p < 0.001).

Conclusion

Our MR study revealed that the serum total bilirubin was causally associated with the risk of cholelithiasis, and the genetic predisposition to cholelithiasis was causally associated with the increased serum total bilirubin levels.

A Mendelian randomization study on causal effects of 25(OH) vitamin D levels on diabetic nephropathy

BACKGROUND

Vitamin D supplementation is associated with a lower incidence of diabetic nephropathy (DN); however, whether this association is causative is uncertain.

METHODS

We used two-sample Mendelian randomization to examine the causal influence of vitamin D on diabetic nephropathy in 7,751 individuals with type I diabetes-related nephropathy (T1DN) and 9,933 individuals with type II diabetes-related nephropathy (T2DN). Meanwhile, we repeated some previous studies on the influence of KIM-1 (kidney injury molecule 1) and body mass index (BMI) on DN. Additionally, to test the validity of the instruments variable for vitamin D, we conducted two negative controls Mendelian randomization (MR) on breast and prostate cancer, and a positive control MR on multiple sclerosis.

RESULTS

Results of the MR analysis showed that there was no causal association between 25(OH)D with the early/later stage of T1DN (early: OR = 0.903, 95%CI: 0.229 to 3.555; later: OR = 1.213, 95%CI: 0.367 to 4.010) and T2DN (early: OR = 0.588, 95%CI: 0.182 to 1.904; later: OR = 0.904, 95%CI: 0.376 to 2.173), nor with the kidney function of patients with diabetes mellitus: eGFRcyea (creatinine-based estimated GFR) (Beta = 0.007, 95%CI: -0.355 to 0.369)) or UACR (urinary albumin creatinine ratio) (Beta = 0.186, 95%CI: -0.961 to 1.333)).

CONCLUSIONS

We found no evidence that Vitamin D was causally associated with DN or kidney function in diabetic patients.

Evidence of causality of low body mass index on risk of adolescent idiopathic scoliosis: a Mendelian randomization study

Introduction

Adolescent idiopathic scoliosis (AIS) is a disorder with a three-dimensional spinal deformity and is a common disease affecting 1-5% of adolescents. AIS is also known as a complex disease involved in environmental and genetic factors. A relation between AIS and body mass index (BMI) has been epidemiologically and genetically suggested. However, the causal relationship between AIS and BMI remains to be elucidated.

Material and methods

Mendelian randomization (MR) analysis was performed using summary statistics from genome-wide association studies (GWASs) of AIS (Japanese cohort, 5,327 cases, 73,884 controls; US cohort: 1,468 cases, 20,158 controls) and BMI (Biobank Japan: 173430 individual; meta-analysis of genetic investigation of anthropometric traits and UK Biobank: 806334 individuals; European Children cohort: 39620 individuals; Population Architecture using Genomics and Epidemiology: 49335 individuals). In MR analyses evaluating the effect of BMI on AIS, the association between BMI and AIS summary statistics was evaluated using the inverse-variance weighted (IVW) method, weighted median method, and Egger regression (MR-Egger) methods in Japanese.

Results

Significant causality of genetically decreased BMI on risk of AIS was estimated: IVW method (Estimate (beta) [SE] = -0.56 [0.16], p = 1.8 × 10-3), weighted median method (beta = -0.56 [0.18], p = 8.5 × 10-3) and MR-Egger method (beta = -1.50 [0.43], p = 4.7 × 10-3), respectively. Consistent results were also observed when using the US AIS summary statistic in three MR methods; however, no significant causality was observed when evaluating the effect of AIS on BMI.

Conclusions

Our Mendelian randomization analysis using large studies of AIS and GWAS for BMI summary statistics revealed that genetic variants contributing to low BMI have a causal effect on the onset of AIS. This result was consistent with those of epidemiological studies and would contribute to the early detection of AIS.

The genetic basis of major depressive disorder

The genetic dissection of major depressive disorder (MDD) ranks as one of the success stories of psychiatric genetics, with genome-wide association studies (GWAS) identifying 178 genetic risk loci and proposing more than 200 candidate genes. However, the GWAS results derive from the analysis of cohorts in which most cases are diagnosed by minimal phenotyping, a method that has low specificity. I review data indicating that there is a large genetic component unique to MDD that remains inaccessible to minimal phenotyping strategies and that the majority of genetic risk loci identified with minimal phenotyping approaches are unlikely to be MDD risk loci. I show that inventive uses of biobank data, novel imputation methods, combined with more interviewer diagnosed cases, can identify loci that contribute to the episodic severe shifts of mood, and neurovegetative and cognitive changes that are central to MDD. Furthermore, new theories about the nature and causes of MDD, drawing upon advances in neuroscience and psychology, can provide handles on how best to interpret and exploit genetic mapping results.

Missing Causality and Heritability of Autoimmune Hepatitis

BACKGROUND

Autoimmune hepatitis has an unknown cause and genetic associations that are not disease-specific or always present. Clarification of its missing causality and heritability could improve prevention and management strategies.

AIMS

Describe the key epigenetic and genetic mechanisms that could account for missing causality and heritability in autoimmune hepatitis; indicate the prospects of these mechanisms as pivotal factors; and encourage investigations of their pathogenic role and therapeutic potential.

METHODS

English abstracts were identified in PubMed using multiple key search phases. Several hundred abstracts and 210 full-length articles were reviewed.

RESULTS

Environmental induction of epigenetic changes is the prime candidate for explaining the missing causality of autoimmune hepatitis. Environmental factors (diet, toxic exposures) can alter chromatin structure and the production of micro-ribonucleic acids that affect gene expression. Epistatic interaction between unsuspected genes is the prime candidate for explaining the missing heritability. The non-additive, interactive effects of multiple genes could enhance their impact on the propensity and phenotype of autoimmune hepatitis. Transgenerational inheritance of acquired epigenetic marks constitutes another mechanism of transmitting parental adaptations that could affect susceptibility. Management strategies could range from lifestyle adjustments and nutritional supplements to precision editing of the epigenetic landscape.

CONCLUSIONS

Autoimmune hepatitis has a missing causality that might be explained by epigenetic changes induced by environmental factors and a missing heritability that might reflect epistatic gene interactions or transgenerational transmission of acquired epigenetic marks. These unassessed or under-evaluated areas warrant investigation.

Associations of the circulating levels of cytokines with risk of amyotrophic lateral sclerosis: a Mendelian randomization study

BACKGROUND

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disorder that is accompanied by muscle weakness and muscle atrophy, typically resulting in death within 3-5 years from the disease occurrence. Though the cause of ALS remains unclear, increasing evidence has suggested that inflammation is involved in the pathogenesis of ALS. Thus, we performed two-sample Mendelian randomization (MR) analyses to estimate the associations of circulating levels of cytokines and growth factors with the risk of ALS.

METHODS

Genetic instrumental variables for circulating cytokines and growth factors were identified from a genome-wide association study (GWAS) of 8293 European participants. Summary statistics of ALS were obtained from a GWAS including 20,806 ALS cases and 59,804 controls of European ancestry. We used the inverse-variance weighted (IVW) method as the primary analysis. To test the robustness of our results, we further performed the simple-median method, weighted-median method, MR-Egger regression, and MR pleiotropy residual sum and outlier test. Finally, a reverse MR analysis was performed to assess the possibility of reverse causation between ALS and the cytokines that we identified.

RESULTS

After Bonferroni correction, genetically predicted circulating level of basic fibroblast growth factor (FGF-basic) was suggestively associated with a lower risk of ALS [odds ratio (OR): 0.74, 95% confidence interval (95% CI): 0.60-0.92, P = 0.007]. We also observed suggestive evidence that interferon gamma-induced protein 10 (IP-10) was associated with a 10% higher risk of ALS (OR: 1.10, 95% CI: 1.03-1.17, P = 0.005) in the primary study. The results of sensitivity analyses were consistent.

CONCLUSIONS

Our systematic MR analyses provided suggestive evidence to support causal associations of circulating FGF-basic and IP-10 with the risk of ALS. More studies are warranted to explore how these cytokines may affect the development of ALS.