key: cord-1044068-e5fdh69h
authors: Zhang, Min; Bai, Ye; Wang, Yutong; Cui, Huijie; Tang, Mingshuang; Wang, Lanbing; Wang, Xin; Gu, Dongqing
title: Cumulative Evidence for Associations Between Genetic Variants in Interleukin 6 Receptor Gene and Human Diseases and Phenotypes
date: 2022-04-14
journal: Front Immunol
DOI: 10.3389/fimmu.2022.860703
sha: 3f7243f58f365bfad9704868d7057ba6c4e923ca
doc_id: 1044068
cord_uid: e5fdh69h

BACKGROUND: Genetic studies have linked polymorphisms in the interleukin 6 receptor (IL6R) gene to the risk of multiple human diseases and phenotypes, yet have reported inconsistent results. We aimed to synthesize current knowledge of variants in the IL6R gene on the risk of diseases and phenotypes. METHODS: We searched the Medline and Embase databases to identify relevant publications. Meta-analysis was performed utilizing DerSimonian and Laird random-effects model. We also graded cumulative evidence for significant associations. Furthermore, phenome-wide analyses and functional annotations were performed for variants with strong evidence. RESULTS: We included 155 studies for evaluating the associations between 80 polymorphisms in the IL6R gene and the risk of 102 human diseases and 98 phenotypes. We conducted 58 main meta-analyses, and 41 significant associations were identified. Strong evidence was assigned to 29 associations that investigated ten variants (rs2228145, rs4129267, rs7529229, rs4537545, rs7518199, rs4845625, rs4553185, rs4845618, rs4845371, and rs6667434) related to the risk of four cardiovascular diseases (coronary heart disease, coronary artery disease, atherosclerosis, and abdominal aortic aneurysms), four inflammatory diseases (rheumatoid arthritis, Crohn’s disease, dermatitis, and asthma), and concentration of four phenotypes (C-reactive protein, fibrinogen, IL-6, and sIL-6R). Furthermore, phenome-wide analysis verified that rs2228145 associated with asthma and dermatitis risk. Functional analyses indicated that these polymorphisms fall within exon, enhancer regions. CONCLUSIONS: Our study comprehensively summarizes current data on the genetic architecture of the IL6R gene and highlights the pharmacological targeting potential of IL-6R on cardiovascular and inflammatory diseases.

Interleukin (IL)-6 receptor (IL6R) gene is located in chromosome 1q21.3 (Gene ID: 3570) (1). Its gene production, IL-6R, could specifically bind to its ligand (IL-6) and subsequently recruit two molecules of glycoprotein 130 (gp130) to initiate intercellular signaling (2, 3) . IL-6/IL-6R/gp130 signaling is implicated in the development of inflammatory and autoimmune diseases, including novel coronavirus disease 2019 (COVID-19) (4) .

Genetic studies have indicated that variants in the IL6R gene are associated with a range of common diseases and phenotypes. In 2003, Kim et al. firstly identified seven single-nucleotide polymorphisms (SNPs) in the IL6R gene in the Korean population (5) . One of these variants, Asp358Ala (rs2228145, previously rs8192284), a non-synonymous located in exon 9, has been shown to be associated with the index of obesity in Indians (6) in the same year. Subsequent studies have implicated that Asp358Ala was associated with a series of inflammatory/ autoimmune diseases, including inflammatory bowel disease (7), rheumatoid arthritis (8) , cardiovascular diseases (CVDs) (9, 10), and type 2 diabetes (11) . In 2012, rs2228145 was implicated to play a protective effect in the etiology of CVD diseases, including coronary heart disease (CHD), atrial fibrillation (AF), and abdominal aortic aneurysms (AAA) (12) . This effect was verified by findings from a large phenome-wide association analysis among populations from the Million Veteran Program, which showed that rs2228145 was significantly associated with 22 diseases (13) . In addition to CVDs, rs2228145 was found to decrease the risk of ulcerative colitis (UC) or Crohn's disease (CD) by a genome-wide meta-analysis of 38, 197 patients with IBD and more than 40,000 controls (7) . A genomewide association study (GWAS) of European ancestry identified that rs2228145 also decreased risk of rheumatoid arthritis (8) , and this association has been replicated by subsequent studies (14) . Additionally, rs2228145 plays a protective effect in the etiology of type I diabetes (T1D) (12) while positively associated with dermatitis and renal function disorder (13) .

Other polymorphisms in the IL6R gene, such as rs4129267 (T allele), rs4845625 (C allele), or rs7529229 (C allele), were also reported to be associated with the risk of many diseases, such as asthma (15) and CVDs (16) . A GWAS in Australia identified that rs4129267 (T allele) significantly increased risk of asthma (OR = 1.09, combined p = 2.4×10 −8 ) (15) . Then, another GWAS in the Caucasian population indicated that this variant was also associated with an increased risk of AAA. Other large-scale genetic analyses identified that rs4845625 (C allele) and rs7529229 (C allele) were implicated in risk of CVDs, such as atrial fibrillation (17) , coronary artery disease (CAD) (16) , AAA, and CHD (10) .

In addition to diseases, variants in the IL6R gene were also suggested to be related to hematological inflammatory biomarkers and some conventional CVD risk factors. Several publications showed that the C allele of rs2228145 increased circulating soluble IL-6R (sIL-6R) and IL-6 concentrations by 34% and 15%, respectively (7, 9) . Furthermore, rs4845371, rs6667434, or rs7529229 was associated with lower levels of circulation of C-reactive protein (CRP), fibrinogen, troponin I, creatine kinase-MB (CK-MB), pro-B-type natriuretic peptide, and systolic blood pressure, which were consistent with a reduced risk of CVDs (9, 13). With the increasing number of genetic studies published since 2003, several metaanalyses have been performed to investigate the role of polymorphisms in IL6R gene on the risk of human diseases and phenotypes. However, existing meta-analyses were limited to a single variant or illness, and the results might not be adequate partly due to the small sample size (18) (19) (20) (21) . These drawbacks may limit the implementation of the evidence in the clinical practice and disease prevention strategies. Therefore, we aimed to provide a systematic synopsis on the associations between polymorphisms in the IL6R gene and human diseases and phenotypes. Firstly, we conducted a systematic review and meta-analysis to evaluate the epidemiological associations between variants in the IL6R gene and human diseases and phenotypes. Secondly, we graded levels of cumulative evidence for significant associations combining Venice criteria and falsepositive report probability (FPRP) tests. Thirdly, we performed a phenome-wide analysis for independent variants with strong evidence using data from UK Biobank to uncover new relationships. Fourthly, we performed functional annotations for the variants in the IL6R gene using data from the Encyclopedia of DNA Elements Project (ENCODE), the 1000 Genomes Project, the Genotype-Tissue Expression (GTEx) Project, and other sources to provide clues for the pathogenesis of these associations.

The methods used in the present study were based on guidelines proposed by the Human Genome Epidemiology Network (HuGENet) for a systematic review of genetic association studies and the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) (22, 23) . This study has been registered to the International Prospective Register of Systematic Reviews (PROSPERO; registration ID: CRD42020201735).

We systematically searched PubMed and Embase using keywords ("interleukin 6 receptor" OR "IL-6R" OR "IL6R" OR "interleukin-6 receptor" OR "IL-6 receptor" OR "rs2228145" OR "Asp358Ala" OR "D358A") to identify the original genetic association studies of IL6R published in English up to August 1, 2021. We also retrieved published reviews, meta-analyses, pooled analyses, and the included studies to identify additional publications. Articles were eligible to be included if they met the following inclusion criteria: (i) they were observational studies conducted in humans; (ii) they investigated associations of genetic variants in IL6R gene and human diseases or phenotypes; (iii) they directly provided risk estimates such as odds ratio (OR), relative risk (RR), hazard ratio (HR), and 95% confidence interval (CI), coefficient b and standard error (SE), or provided valid data to calculate these risk estimates under additive genetic model; (iv) or for publications investigating continuous quantitative phenotypes, they should provide mean and standard deviation (SD). When multiple studies use the same population, the most recent or the one with the largest sample size was selected. Studies without full text or sufficient information to calculate effect size were excluded. Two investigators (MZ and YB) independently assessed the eligibility of each publication, and disagreement was settled through discussion with the principal author (DG).

Three authors (MZ, YW, and XW) independently extracted the data using a predesigned collection sheet including PMID, first author, publishing year, variants, major and minor alleles, source of population, ethnicity of the study population, study design, sample size of cases and controls, genotype and allele counts, matched information, effects of associations and the corresponding 95% CIs, regression coefficients, SE or p-value (for studies using multiple adjusted models, the most fully adjusted estimates were extracted), mean concentrations of quantitative traits and SD for different genotypes, Hardy-Weinberg equilibrium (HWE) among controls, and genotyping methods.

In our study, the primary outcome was the risk of disease or categorical phenotypes, measured by OR/RR/HR and 95% CI under an additive genetic model. Due to the relatively low incidence of the events mentioned above, the effect size of OR, HR, and RR was deemed identical in this meta-analysis (24) . For continuous quantitative phenotypes, the primary outcome was the mean deviation, measured by standardized mean difference (SMD) and 95% CI after pooling mean concentration and SD under a co-dominant model (mutant-type homozygote vs. wildtype homozygote, heterozygote vs. wild-type homozygote).

Meta-analysis was conducted for variants with at least three independent studies or datasets. We estimated the effects of minor alleles (in Caucasians) on the outcome for each variant. We adopted the methods used in original studies if data were extracted from GWASs or large GWAS meta-analyses from collaborative studies. For heterogeneity assessment, we performed the Cochran's Q test, and the I 2 statistic was used to evaluate and quantify heterogeneity between studies (25, 26) . A random-effects model was performed if heterogeneity was evident (I 2 ≥ 50%), while the fixed effects model was used if I 2 < 50%. For variants that showed significant associations with relevant diseases or phenotypes, sensitivity analyses were conducted to examine whether the associations would be lost when a single study was excluded or when the first positive report or studies with deviation from HWE in controls were excluded. Begg's test (27) and Egger's test (28) were conducted to evaluate potential publication bias and small-study bias, respectively.

Statistical analysis was done using STATA, version 15 (Stata, College Station, TX). Tests of heterogeneity and bias were twotailed, and a p-value of less than 0.10 was considered significant as recommended. All the other statistical tests were two-tailed, and a p-value of < 0.05 was considered statistically significant unless otherwise stated.

Venice criteria were applied to evaluate the epidemiological credibility of significant associations identified by primary metaanalyses (see Supplementary Notes) (29) . Briefly, credibility was determined by three criteria: (i) amount of evidence, (ii) replication, and (iii) protection from bias. Based on the grade of A, B, or C in these criteria, epidemiological credibility was defined as strong when all three grades were A, moderate when all three grades were A or B, and weak when any grade was C.

We also estimated the noteworthiness of the significant associations by calculating FPRP (30) . A prior probability was set at 0.05 to estimate the FPRP value based on the OR obtained from a meta-analysis, and we adopted an FPRP cutoff value of 0.20. Briefly, FPRP value < 0.05, 0.05 < FPRP value < 0.20, and FPRP value > 0.20 were defined as strong, moderate, and weak evidence of true association, respectively. The cumulative epidemiological evidence was adjusted based on the FPRP results: the evidence may be upgraded if FPRP evidence was strong and downgraded if the FPRP result was weak.

We performed linkage disequilibrium (LD) analysis for variants significantly associated with risk of diseases and phenotypes in the present meta-analyses, and r 2 > 0.80 was considered as significant using data from the 1000 Genomes Project (31) . We then conducted a phenome-wide analysis to estimate associations of the independent variants with 778 phenotypes from UK Biobank in GeneATLAS through querying summary data (32) . We considered p-values < 6.43 × 10 -5 (0.05/778) as significant, adjusted for multiple comparisons of variants and 778 phenotypes.

We conducted analyses to investigate the potential regulatory effect of polymorphisms significantly associated with relative diseases or phenotypes using the Encyclopedia of DNA Elements (ENCODE) tool HaploReg v4.1 (33) and the University of California, Santa Cruz (UCSC) Genome Browser (34) . We examined genome-wide cis-eQTL data in multiple tissues from the Blood eQTL browser (35) , and the GTEx Project (36) . The significance threshold for these analyses was set to p-values of less than 0.005 accounting for 10 tests. In addition, we searched the published literature for the IL6R gene in relation to diseases from PubMed.

The literature search strategy is shown in Figure 1 . A total of 155 studies were included for evaluating the associations between 80 polymorphisms in the IL6R gene and the risk of human diseases and phenotypes. Of the 155 publications, 114 regarding the association between 67 variants in the IL6R gene and 102 diseases (Supplementary Table 1 ) and 67 regarding 50 variants associated with 98 phenotypes (Supplementary Tables 2, 3) were included. A total of 131 articles were cross-sectional or casecontrol studies, and 24 articles were longitudinal studies. As for the study population, 94 were conducted in populations of European ancestry, 34 were in Asians, four were in the African population, four were in the population of mixed ancestry, and 19 were on other ancestries. The mean number of cases and controls in these articles was 3,855 (range, 17 to 51,441) and 33,415 (range, 20 to 897,488), respectively. Of the 155 studies, 134 were published since 2010 (86.345%).

Gene and Human Diseases/Phenotypes by the Meta-Analyses

We conducted 27 main meta-analyses to estimate associations of six variants in IL6R gene with risk of 20 kinds of diseases under an additive model. Summary findings of meta-analyses are shown in Table 1 . Of the 27 main analyses, 21 associations between six variants and 14 diseases indicated statistical significance (p < 0.05). Specifically, the coronary heart disease risk was negatively associated with minor allele of rs7529229 (OR = 0.954, 95% CI: 0.932-0.977), rs4537545 (OR = 0.940, 95% CI: 0.903-0.978), and rs2228145 (OR = 0.957, 95% CI: 0.948-0.967). Atrial fibrillation risk was inversely associated with minor allele of rs7529229 (OR = 0.900, 95% CI: 0.850-0.960), rs4537545 (OR = 0.901, 95% CI: 0.847-0.952), and rs28638007 (OR = 0.900, 95% CI: 0.850-0.950), but positively associated with rs4845625 (OR = 1.179, 95% CI: 1.044-1.330). The abdominal aortic aneurysms risk was negatively associated with rs7529229 (OR = 0.841, 95% CI: 0.800-0.884), but it was positively associated with rs4129267 (OR = 1.141, 95% CI: 1.101-1.183). As for asthma, rs4129267 (OR = 1.073, 95% CI: 1.035-1.112) and rs2228145 (OR = 1.053, 95% CI: 1.023-1.084) were associated with increased risk of asthma. In addition, the CAD risk was positively associated with rs4845625 (OR = 1.060, 95% CI: 1.042-1.079).

Furthermore, rs2228145 (C allele) was associated with 1.026-1.071). However, we did not find significant associations between rs2228145 and risk of obesity, dengue, mental disorder, multiple myeloma, chronic obstructive pulmonary disease (COPD), and anterior cruciate ligament injury (p > 0.05).

As shown in Table 2 , we performed 31 meta-analyses for associations between ten variants in IL6R gene and 13 kinds of categorical phenotypes under an additive model. Twenty significant associations were identified, which investigated ten variants associated with five kinds of phenotypes (p < 0.05) ( We also conducted meta-analyses for rs2228145 and 12 continuous phenotypes under a co-dominant model (Supplementary Table 4 ). The results showed that rs2228145 was significantly associated with the levels of IL-6 and CRP.

As shown in Tables 1 and 2, 24 associations revealed no or little heterogeneity (I 2 < 25%), ten associations revealed moderate heterogeneity (25% > I 2 <50%) and 24 associations revealed strong heterogeneity (I 2 > 50%). The results indicated that publication bias was found in three associations (p < 0.10), and a small study bias was found in nine associations (p < 0.10) (Supplementary Table 5 ). In addition, sensitivity analysis suggested that most associations did not significantly change when a single or first positive study, small studies, or studies with deviation from HWE in controls were excluded (data not shown).

Using Venice criteria and FPRP tests, we graded cumulative evidence level for the 41 significant associations, which evaluated associations between 11 variants in the IL6R gene and 14 diseases and five phenotypes ( Tables 1, 2 and Supplementary Table 5 ). According to Venice criteria, strong, moderate, and weak evidence were assigned to 29, four, and 12 associations, respectively. The FPRP test suggested that 34 of them had an FPRP value less than 0.05, two associations had an FPRP value between 0.05 and 0.20, and five associations had an FPRP value greater than 0.20. Combining Venice Criteria and FPRP results, strong, moderate, and weak evidence were finally assigned to 29, 8, and four associations, respectively. The 29 strong associations evaluated ten polymorphisms related to nine diseases and four phenotypes.

Linkage disequilibrium analysis (Supplementary Table 6 ) showed that five variants (rs2228145, rs7518199, rs4537545, rs7529229, and rs4129267) were highly correlated with each other, while another five SNPs (rs4845625, rs4845618, rs4845371, rs6667434, and rs4553185) were also highly correlated with each other according to r 2 (r 2 > 0.8) in the European population.

Therefore, phenome-wide association analysis was just conducted for two SNPs (rs2228145 and rs4845625) using data from UK Biobank. The results are shown in Table 3 and Figure 2 . Variant rs2228145 and its highly correlated SNPs were found to be significantly associated with three binary phenotypes (asthma, eczema/dermatitis, and other dermatitis) and 11 non-binary phenotypes (monocyte percentage, monocyte count, mean corpuscular hemoglobin, mean corpuscular volume, mean platelet volume, platelet count, red blood cell distribution width, hemoglobin concentration, platelet distribution width, lymphocyte count, and mean sphered cell volume) at p < 6.43 × 10 -6 . Variant rs4845625 and its highly correlated SNPs were significantly associated with four non-binary phenotypes of mean platelet (thrombocyte) volume, platelet count, mean corpuscular hemoglobin, and monocyte count.

Bioinformatics analyses suggested that the ten variants with strong evidence and variants highly correlated with them might fall within an exon, a DNase I hypersensitivity site, a strong prompter, and an enhancer activity region ( Table 4 and Supplementary Figure 1) . Of those SNPs, rs2228145 is a nonsynonymous variant located in exon 9 of the IL6R gene. In variant rs2228145, A to C transition is observed at position 358 in the IL6R gene, resulting in the alteration of amino acid from alanine to aspartate or valine. The C allele of rs2228145 is associated with inflammatory/autoimmune diseases as well as quantitative traits (7, 9) , and the mechanisms might be involved in regulating the balance of transmembrane IL-6R (mIL-6R) and sIL-6R concentrations and responding to IL-6 (12). Additionally, rs2228145 is in LD with four intronic SNPs (rs4129267, rs4537545, rs7529229, and rs7518199) in Europeans (r 2 > 0.85 for all tests), indicating similar functional mechanisms in this ethnic group. According to LD analysis in the European population, cis-eQTL analysis was conducted for rs2228145 and rs4845625 (Supplementary Table 7 ). The findings indicate that rs2228145 and rs4845625 are associated with a decrease in IL-6R in an artery or heart tissue. In addition, rs2228145 is an eQTL for SHE, TDRD10, and UBE2Q1 genes, and rs4845625 is an eQTL for PSMD8P1, TDRD10, UBAP2L, and AQP10 genes.

Our study summarizes the published data thus far on relationships between variants in the IL6R gene and the risk of human diseases and phenotypes. The meta-analyses results show that 12 variants in the IL6R gene are significantly associated with risk of CVDs (CAD, coronary heart disease, atrial fibrillation, aneurysms, and atherosclerosis), inflammatory or autoimmune diseases (rheumatoid arthritis, asthma, Crohn's disease, ulcerative colitis, dermatitis, and type 1 diabetes), and levels of inflammatory biomarkers (CRP, fibrinogen, IL-6, IL-6R, and sIL-6R). We then assigned strong cumulative evidence to 29 significant associations, which assessed ten variants in this gene related to the risk of five cardiovascular diseases, four inflammatory diseases, and the concentration of four phenotypes. Further analysis of functional annotation indicates that the strong variants might fall in putative functional regions.

IL6R gene encodes a subunit of an IL-6R protein complex that exists in two forms, an mIL-6R and sIL-6R, resulting in different signal transduction mechanisms known as IL-6/IL-6R classic signaling and trans-signaling, respectively. These two models of signaling can mediate the anti-inflammatory (classic signaling) and pro-inflammatory (trans-signaling) response of IL-6 (37, 38) . By specifically binding to IL-6, IL-6/IL-6R signaling plays an essential role in the pathogenesis of many diseases, including COVID-2019 (4, 7, 9, 10) . Previous studies have suggested that rs2228145 has similar synergistic effects as tocilizumab on a series of cardiovascular and inflammatory biomarkers (10, 12, 19) . Consistent with this, our meta-analysis found that this variant was significantly related to a decreased serum concentration of CRP and fibrinogen. There were also consistent associations with strong cumulative evidence between rs2228145 and a reduced risk of CHD, atherosclerosis, AAA, rheumatoid arthritis, and inflammatory bowel disease identified by this meta-analysis, findings that suggest targeting the IL-6R is a plausible strategy in cardiovascular and inflammatory diseases. On the contrary, the C allele of rs2228145 is strongly related to higher levels of sIL-6R and IL-6, which seems paradoxical. This apparent paradox might be interpreted that the C allele reduces mIL-6R level either by alternative splicing or increasing the membrane-bound form shedding rather than by regulating the production of IL-6R or IL-6 (9). Such mechanisms could account for the accumulation of both sIL-6R and IL-6 in the circulation of individuals carrying the C allele of rs2228145.

The mutation of rs2228145 promotes the proteolysis of the IL-6R, which increases the level of sIL-6R and decreases the level of mIL-6R (39) . Thus, this variant exerts its functional mechanism by regulating the balance between mIL-6R and sIL-6R (12) . The minor allele of rs2228145 strongly increases the concentration of sIL-6R, whereas it reduces the surface expression of mIL-6R in individual immune cells. Reducing mIL-6R expression leads to impaired IL-6 responsiveness, as measured by decreased phosphorylation of the transcription factors STAT3 and STAT1 following stimulation with IL-6 (12, 19) . Therefore, it is speculated that rs2228145 impairs classical IL-6/IL-6R signaling and reduces the inflammatory response, thus leading to a lower risk of CVDs, inflammatory, or autoimmune diseases (19) . Our findings also link rs2228145 with an increased risk of asthma and dermatitis, but the exact mechanisms for these associations remain unclear. Experimental evidence indicates that IL-6/sIL-6R trans-signaling might play an essential role in the development of asthma (40, 41) and dermatitis (42) . Further studies are warranted to elucidate these mechanisms.

LD analysis using data from the 1000 Genomes Project reveals that four variants (rs7518199, rs4537545, rs7529229, and rs4129267) are highly correlated with rs2228145 in the European population (r 2 > 0.90 for all tests), indicating that the functional mechanisms of these SNPs associated with diseases and phenotypes might be similar. In addition, the variant rs4845625 (T allele) strongly associated with CAD is in weak or moderate LD with rs2228145 (r 2 = 0.04, 0.462, and 0.63 in African, European, and Asian populations, respectively), reflecting that functional mechanisms of IL6R variants associated with CVD risk might be distinct. Additional studies are needed to investigate rs4845625 and its related SNPs affecting the development of CVD. Clinically, targeting the IL-6/IL-6R pathway is a plausible therapeutic strategy in several cardiovascular and inflammatory diseases, including COVID-2019 (37, (43) (44) (45) . For example, clinical trials suggest that tocilizumab, a humanized monoclonal biologic drug that blocks both sIL-6R and mIL-6R, is an effective treatment for rheumatoid arthritis (44) . In addition, sgp130Fc, a protein that specifically blocks trans-signaling, is an effective therapy for various preclinical chronic and autoimmune disease models (37, 42) . Therefore, our findings have clinical implications and may provide tools for identifying patients who can benefit from the therapeutic intervention in the IL-6/IL-6R pathway.

Some limitations should be noted in this study. Firstly, although we have tried our best to perform a comprehensive meta-analysis to include all the relevant literature, a small part of articles may still be missed. Secondly, several studies were omitted without valid data for extraction. For example, some articles that just provide a histogram of the average level of phenotypes corresponding to different genotypes were excluded. Thirdly, meta-analysis was not performed for some diseases and phenotypes, such as breast cancer, colorectal cancer, stroke, and systemic lupus erythematosus because fewer than three datasets were included. Fourthly, we could not define a unified analytical standard across studies since we could not obtain raw data from the original studies. To minimize bias, we adopted age, sex, or multivariable-adjusted risk estimates. Fifthly, moderate and weak evidence of significant associations needs to be explained with caution. Future studies with a larger sample size are required to refute or confirm these associations. FIGURE 2 | Plot of phenome-wide analysis for rs2228145 using data from UB Biobank. We estimated rs2228145 associated with 778 phenotypes in UK Biobank cohort. Phenotypes are grouped along the x-axis by categorization within the PheWAS code hierarchy. We considered a p-value < 6.43 × 10 −5 (the red line) as significant. In conclusion, our study provides comprehensive evidence that variants in the IL6R gene are associated with the risk of inflammatory or autoimmune diseases and levels of inflammatory biomarkers. Also, it highlights the significant role of the IL-6/IL-6R pathway in the pathogenesis of cardiovascular and inflammatory diseases. Phenome-wide association analysis confirmed that variants in the IL6R gene are associated with asthma and dermatitis. Bioinformatics analyses showed that most of these variants fall in putative functional regions. Those findings highlight that the variations in the IL6R gene may be considered a useful genetic tool for investigating the pharmacological targeting potential of IL-6R. More efforts should be made to understand the biological pathways further and apply these lines of evidences for clinical practice and public health for risk evaluation and management.

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Ethical review and approval were not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

DG initiated and designed this study and approved the final draft manuscript. MZ and YB performed study selection, data extraction, and statistical analysis. LW, HC, and MT performed the cumulative evidence assessments and phenomewide analysis. YW and XW conducted the functional analyses. MZ also drafted the first manuscript. All authors contributed to the article and approved the submitted version.

This study was supported by the National Natural Science Foundation of China (81903393 and 81903398). The funding agencies of this study had no role in study design, data collection, data management, data analysis, data interpretation, writing of the manuscript, or the decision for submission.

The Contribution of the Functional IL6R Polymorphism Rs2228145, eQTLs and Other Genome-Wide SNPs to the Heritability of Plasma sIL-6R Levels

Molecular Characterization of Grass Carp Interleukin-6 Receptor and the Agonistic Activity of Its Soluble Form in Head Kidney Leucocytes

Receptor System and its Role in Physiological and Pathological Conditions

Correlation Between Interleukin Gene Polymorphisms and Current Prevalence and Mortality Rates Due to Novel Coronavirus Disease 2019 (COVID-2019) in 23 Countries

Identification of Novel SNPs in the Interleukin 6 Receptor Gene (IL6R)

Variants in the Interleukin 6 Receptor Gene Are Associated With Obesity in Pima Indians

Variation in Interleukin 6 Receptor Gene Associates With Risk of Crohn's Disease and Ulcerative Colitis

High-Density Genetic Mapping Identifies New Susceptibility Loci for Rheumatoid Arthritis

Interleukin-6 Receptor Pathways in Coronary Heart Disease: A Collaborative Meta-Analysis of 82 Studies

The Interleukin-6 Receptor as a Target for Prevention of Coronary Heart Disease: A Mendelian Randomisation Analysis

Variation in the Interleukin-6 Receptor Gene Associates With Type 2 Diabetes in Danish Whites

Functional IL6R 358Ala Allele Impairs Classical IL-6 Receptor Signaling and Influences Risk of Diverse Inflammatory Diseases

Association of Interleukin 6 Receptor Variant With Cardiovascular Disease Effects of Interleukin 6 Receptor Blocking Therapy: A Phenome-Wide Association Study

Interleukin 6 Receptor (IL6-R) Gene Polymorphisms Underlie Susceptibility to Rheumatoid Arthritis

Identification of IL6R and Chromosome 11q13.5 as Risk Loci for Asthma

Large-Scale Association Analysis Identifies New Risk Loci for Coronary Artery Disease

IL-6R Variants Associated With Diseases/Phenotypes

Large-Scale Candidate Gene Analysis in Whites and African Americans Identifies IL6R Polymorphism in Relation to Atrial Fibrillation: The National Heart, Lung, and Blood Institute's Candidate Gene Association Resource (CARe) Project

Association of IL-6 Promoter and Receptor Polymorphisms With Multiple Myeloma Risk: A Systematic Review and Meta-Analysis

Interleukin-6 Receptor Pathways in Abdominal Aortic Aneurysm

An Association Study Between Genetic Polymorphism in the Interleukin-6 Receptor Gene and Coronary Heart Disease

Meta-Analysis of Cytokine and Chemokine Genes in Schizophrenia

Systematic Reviews of Genetic Association Studies. Human Genome Epidemiology Network

Preferred Reporting Items for Systematic Review and Meta-Analysis Protocols (PRISMA-P) 2015: Elaboration and Explanation

Meta-Analysis of Soy Consumption and Gastrointestinal Cancer Risk

Quantitative Synthesis in Systematic Reviews

Quantifying Heterogeneity in a Meta-Analysis

Operating Characteristics of a Rank Correlation Test for Publication Bias

Bias in Meta-Analysis Detected by a Simple

Assessment of Cumulative Evidence on Genetic Associations: Interim Guidelines

Assessing the Probability That a Positive Report Is False: An Approach for Molecular Epidemiology Studies

An Integrated Map of Genetic Variation From 1,092 Human Genomes

An Atlas of Genetic Associations in UK Biobank

HaploReg: A Resource for Exploring Chromatin States, Conservation, and Regulatory Motif Alterations Within Sets of Genetically Linked Variants

Systematic Identification of Trans eQTLs as Putative Drivers of Known Disease Associations

The Genotype-Tissue Expression (GTEx) Pilot Analysis: Multitissue Gene Regulation in Humans

The Soluble Interleukin 6 Receptor: Advanced Therapeutic Options in Inflammation

Activation of the JAK/STAT3 and PI3K/AKT Pathways Are Crucial for IL-6 Trans-Signaling-Mediated Pro-Inflammatory Response in Human Vascular Endothelial Cells

The Interleukin-6 Receptor Asp358Ala Single Nucleotide Polymorphism Rs2228145 Confers Increased Proteolytic Conversion Rates by ADAM Proteases

Therapeutic Targeting of the IL-6 Trans-Signaling/Mechanistic Target of Rapamycin Complex 1 Axis in Pulmonary Emphysema

The IL6R Gene Polymorphisms Are Associated With sIL-6r, IgE and Lung Function in Chinese Patients With Asthma

Oral Administration of Lactobacillus Delbrueckii Subspecies Bulgaricus OLL1073R-1 Suppresses Inflammation by Decreasing Interleukin-6 Responses in a Murine Model of Atopic Dermatitis

Tocilizumab in the Treatment of Rheumatoid Arthritis: An Evidence-Based Review and Patient Selection

Targeting Interleukin-6 in Inflammatory Autoimmune Diseases and Cancers

Tocilizumab for the Treatment of Severe COVID-19 Pneumonia With Hyperinflammatory Syndrome and Acute Respiratory Failure: A Single Center Study of 100 Patients in Brescia, Italy

Conflict of Interest: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022. 860703/full#supplementary-material