key: cord-0022837-6vhvxbv8
authors: Wang, Chih-Chien; Huang, Chih-Yun; Lee, Meng-Chang; Tsai, Dung-Jang; Wu, Chia-Chun; Su, Sui-Lung
title: Genetic association between TNF-α G-308A and osteoarthritis in Asians: A case–control study and meta-analysis
date: 2021-11-04
journal: PLoS One
DOI: 10.1371/journal.pone.0259561
sha: 3e81bdca24a38728872c0fc6b57759bc27fbaa8e
doc_id: 22837
cord_uid: 6vhvxbv8

BACKGROUND: Osteoarthritis (OA) is an important health issue in elderly people. Many studies have suggested that genetic factors are important risk factors for OA, of which tumor necrosis factor-α (TNF-α) is one of the most examined genes. Moreover, several studies have investigated the relationship between TNF-α G-308A polymorphisms and OA risk, but consistent results have not been obtained. OBJECTIVE: This study examines the association between TNF-α G-308A polymorphisms and knee OA. Moreover, meta-analysis and trial sequential analysis (TSA) was used to determine whether this is a susceptibility gene for knee OA. METHODS: Between 2015 and 2019, 591 knee OA cases and 536 healthy controls were recruited. The Kellgren–Lawrence grading system was used to identify the knee OA cases. A meta-analysis was conducted including related studies published until 2020 from PubMed, Embase, and previous meta-analysis to improve the evidence level of the current study. The results were expressed as odds ratios (ORs) with corresponding 95% confidence intervals (CI) to evaluate the effect of this polymorphism on knee OA risk. The TSA was used to estimate the sample sizes required in this issue. RESULTS: A nonsignificant association was found between the AA genotype and knee OA [adjusted OR, 0.84; 95% CI, 0.62–1.15) in the recessive model] in the present case–control study, and analysis of other genetic models showed a similar trend. After adding the critical case–control samples for Asians, the TNF-α G-308A, AA genotype exhibited 2.57 times more risk of developing arthritis when compared with the GG + GA genotype (95% CI, 1.56–4.23), and the cumulative samples for TSA (n = 2182) were sufficient to obtain a definite conclusion. CONCLUSIONS: The results of this meta-analysis revealed that the TNF-α G-308A, AA genotype is a susceptible genotype for OA in the Asian population. This study integrated all current evidence to arrive at this conclusion, suggesting that future studies on Asians are not required.

Osteoarthritis (OA) is the most common joint disease that is the main cause of disability among elderly people globally [1] . Genetic factors are particularly important, and the heritability of knee OA is around 45% [2] . Moreover, a previous study suggested that OA is primarily influenced by genetic risk factors due to common population polymorphisms in multiple genes [3] . Thus, identifying more candidate genes and evaluating their effects are desirable.

Tumor necrosis factor-α (TNF-α) is located in the major histocompatibility complex III region (6p21.3) on chromosome 6 [4] . TNF-α participates in inflammatory responses and induces similar responses in other cells [5] . Its signal is mainly transmitted through TNF receptor-1 [6, 7] , thereby regulating other cytokines and growth factors to elicit their effects on the body [8] . The increased inflammation or decreased anabolism results in cartilage structure destruction. Therefore, TNF-α is a susceptibility OA marker [9] . Commonly discussed TNF-α single nucleotide polymorphisms (SNPs) include G-308A and G-238A. Previous studies found that the expression level of TNF-α G-308A gene polymorphism has the highest correlation with OA [10] and that an increased OA risk prevails in the recessive genetic model, which only exists among the Asian population [11] . Therefore, this study aimed to examine the relationship between TNF-α G-308A gene polymorphism and OA.

Twelve studies have currently investigated the association between G-308A polymorphism and knee OA [10, [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] . However, no satisfactory consensus has been reached. Four studies considered that the minor allele (A allele) carriers had a higher risk of knee OA [10, 14, 21, 22] . One study considered that the minor allele (A allele) carriers had a lower risk of knee OA [17] . Seven studies also showed a null association between TNF-α G-308A polymorphism and knee OA [12, 13, 15, 16, [18] [19] [20] . The trial sequential analysis (TSA) provided an opportunity to evaluate whether the most recent conclusions are supported by the current cumulative samples [23] . The sample size provided by the currently available studies is only 1,533, 1,647, and 600 in Asians, Caucasians, and Egyptians [10, [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] , which is not enough for obtaining a decisive conclusion. Thus, this study aims to conduct a case-control study to validate the association between TNF-α G-308A polymorphism and knee OA in Taiwan. A meta-analysis was also performed to improve the evidence level and evaluate whether the latest conclusions are supported by the current cumulative samples using TSA.

This study enrolled 1,127 participants (591 cases and 536 controls) comprising individuals �65 years old. All received Taipei City senior medical check-ups between  January 2015 and December 2019 at the TSGH, a medical teaching hospital at the National  Defense Medical Center, Taipei, Taiwan. The check-up is a government-driven welfare program for individuals �65 years old who are registered Taipei City residents for >1 year. Patient information was examined while the participants underwent check-ups. All participants who received the study information, understood the process, and provided written consent were enrolled. Participants were excluded if sufficient blood samples were drawn. The exclusion criteria were patients who had no knee X-ray or genotype data.

2.1.3 Genomic DNA extraction and genotyping. Approximately 10 mL of peripheral blood was intravenously extracted from participants by a physician or nurse. Genomic DNA was isolated using standard procedures with proteinase K (Invitrogen, Carlsbad, CA, USA) digestion and phenol/chloroform methods. The rs1800629 SNP was genotyped using the iPLEX Gold SNP genotyping method. Genotyping was performed under blind conditions. At least 10% of the samples were randomly selected for repeated genotyping to validate the results, and the concordance rate was 99%.

Continuous variables of the general demographic data were expressed as mean and standard deviation using Student's t-test. Differences in genotype and allele frequencies between knee OA patients and healthy controls were tested using the χ 2 test. Odds ratios (ORs) and 95% confidence intervals (CIs) for the risk of knee OA were calculated using logistic regression. The primary analysis was the association test based on the recessive model assumption [11] . The other assumptions of allele type, genotype, and dominant models were also presented. A p value of <0.05 was considered significant. R 3.4.4 was used for statistical analyses.

The PRISMA checklist and Meta-analysis on Genetic Association Studies Checklist are described in S1 Table [24] . The samples examined in this study comprised the general population, and the correlation between TNF-α G-308A and OA risk was examined. This study used synonyms of "TNF α G-308A" and "osteoarthritis" to search PubMed, EMBASE, and Cochrane databases for papers until 17 November 2020 (see S2 Table for details). The language of the articles was limited to English. In addition, the papers included in the meta-analysis studies were manually examined to avoid the omission of important papers. The papers conforming to the following conditions were included in the current study: (1) case-control studies or cross-sectional studies; (2) clear diagnostic standard for the case group, for e.g., case group (KL � 2) and control group (KL < 2); and (3) samples were from adults �18 years old.

Two reviewers (Chih-Yun Huang and Sui-Lung Su) were responsible for the independent extraction of literature data in this study. The data extracted included the last name of the first author, publication year, country, ethnicity of the study population, and gene distribution in the case and control groups. All extracted papers were assessed using the Newcastle-Ottawa Scale, and the details are shown in S1 Table. 2.2.3 Statistical analysis. All included papers were described using appropriate proportions or mean values. The current meta-analysis used ORs with 95% CIs to examine the correlation between TNF-α G-308A and OA. The I 2 test was used to assess heterogeneity. I 2 >50% denoted moderate to high heterogeneity. Egger's regression and a funnel plot were used to examine the symmetry after combination. The primary analysis was the association test based on the recessive model assumption [11] . The other assumptions of allele type and dominant models were also presented. The random-effects model was used to combine the results. The significance level of this study was set as 0.05. The "metafor" [25] and "meta" [26] packages of R software version 3.3.1 were used. In addition, sensitivity analysis was used to exclude papers that did not conform to the Hardy-Weinberg equilibrium or had irrational data. Moreover, the current study employed TSA to validate whether the meta-analysis result was a definite conclusion [27] . TSA was used for stratification analysis based on race (Caucasian, Asian, and Egyptian). Type 1 error, power, and heterogeneity were set at 0.05, 0.8, and 0%, 65%, and 0% in Caucasian, Asian, and Egyptian populations, respectively. A review of the previous literature showed that the OR of correlation between TNF-α G-308A and OA was around 1.5. In this study, OR was set as 1.5 because the A allele is a possible risk factor. The Taiwan Biobank database and 1000 Genome database were used as references for minor allele frequency, which was 0.11, 0.069, and 0.12 for Asians, Caucasians, and Egyptians, respectively. Table 1 shows the distribution of general demographic variables of the case-control population. This enrolled 1,127 subjects, including 536 subjects in the control group with a mean age of 71.78 ± 6.85 years old (264 men and 272 women) and 591 subjects in the case group with a mean age of 73.57 ± 7.36 years old (214 men and 377 women). The proportion of males was lower in the case group than in the control group (p < 0.001), and the age and body mass index of the case group were higher than those of the control group (p < 0.001 and p = 0.019). Table 2 shows the correlation between TNF-α G-308A and knee OA in the case and control groups. The distribution of the A allele did not demonstrate any significant differences between the control and case groups (p = 0.091) and did not exhibit a significant correlation with OA (OR, 2.74; 95% CI, 0.85-8.81). The dominant and recessive models were further used to validate the results and similarly did not find any significant differences after correcting for covariates. Therefore, the current case-control study found that TNF-α G-308A does not have a significant correlation with OA. The samples from the case-control study were included in 

the meta-analysis and further employed TSA validation to see if definite conclusions could be obtained for Asian samples to increase the power of evidence for the meta-analysis. PubMed. Another 141 papers from Embase were then obtained, and manual searches from other meta-analyses were done in 151 papers after excluding repeated papers. Screening based on the title and abstract was conducted, resulting in the exclusion of 34 retrospective or metaanalysis papers, 75 papers unrelated to this study, eight conference abstracts, and seven animal studies. Based on the screening of the entire manuscripts, four non-knee OA papers, 10 papers containing gene polymorphisms that were not defined in this study, and one paper with repeated samples were noted. Finally, 12 papers were included in the analysis. S3-S5 Tables show the general description of the papers included in the meta-analysis, ethnicities studied in the papers, and the paper quality evaluation. 

The sample size for Caucasians was 1,533 patients after TSA estimation without significant correlation between TNF-α G-308A and OA, and the cumulative sample size did not reach the 

target sample size (2,528). Therefore, definite results could not be obtained in the meta-analysis of Caucasians (S1 Fig). The cumulative sample size for Asians was 1,647 patients before adding samples from this study, and a significant correlation was noted between TNF-α G-308A and OA. The cumulative sample size for Asians was 2,774 after the inclusion of this study, but the cumulative sample size did not reach the target sample size (4,831). Therefore, definite results could not be obtained in the meta-analysis of Asians (Fig 3) . The cumulative sample size for Egyptians was 600 patients, and TNF-α G-308A showed a significant correlation with OA. However, the cumulative sample size did not reach the target sample size (1,617). Therefore, definite results could not be obtained in the meta-analysis of Egyptians (S2 Fig) . 

TNF-α G-308A and OA

The study by Han et al. showed an OR ratio that was too high [14] . However, the studies by Moos TSA estimation was performed after exclusion. Consequently, 981 cumulative Caucasian samples were noted without significant correlation between TNF-α G-308A and OA, and the sample size was insufficient to obtain a definite conclusion (S3 Fig). The cumulative sample size for Asians was 1,055 patients before adding samples from this study, and a significant correlation was noted between TNF-α G-308A and OA. The cumulative sample size for Asians was 2,182 after the inclusion of this study, and the cumulative sample size reached the target 

TNF-α G-308A and OA sample size (1,691; Fig 5) . This result signified that TNF-α G-308A and OA were also significantly correlated in Asians, and a decisive conclusion could be confirmed. Thus, the case-control samples in this study provided critical evidence for establishing a decisive conclusion.

The current case-control study identified no significant correlation between TNF-α G-308A and OA risk. This result is consistent with some of the previous studies [12, 13, 15, 16, [18] [19] [20] . However, other studies have indicated that TNF-α G-308A is correlated with OA risk [10, 14, 17, 21, 22] . G-308A was changed from guanine (G) to adenine (A) [28] , and the performance of TNF-α is affected when the G allele mutates into the A allele [14, 29] . TNF-α enhances articular chondrocyte stiffness and impaired contraction function and is associated with OA [30] . Increased TNF-α mRNA expression in cartilage increases inflammation, thereby stimulating matrix metalloproteinase in chondrocytes, degrading extracellular matrix in cartilage, and disrupting cartilage structure integrity [31, 32] . Recent studies have shown that TNF-α is the 

most important mediator that alters the cartilage matrix degradation and balance in patients with arthritis and ultimately leads to cartilage degradation [33] . Thus, the number and sample size of the selected studies are suggested to be small and insufficient to obtain a conclusive result. Additionally, N4-acetylcytidine (ac4C) modification of mRNA enhances mRNA stability and translation efficiency which correlate with the occurrence, development, and prognosis of diseases. Escalating ac4C levels in the urine may also be associated with inflammatory responses [34] . More investigations on TNF-α and ac4C would help clarify the role of TNF-α in the mechanisms of OA.

The current meta-analysis found an increased risk by using the recessive model analysis. However, a subgroup analysis found that the relationship between TNF-α-308 polymorphism and OA risk only existed among Asians and not among Caucasians, which is consistent with the results of previous meta-analyses [11, 18, 19] . Another study also found different results in different ethnic populations. The study group of North and South Indians separately explored the relationship between TNF-α G-308A and diseases. The results showed an association between southern Indians' genotypes and diseases, but no association was found in northern Indians, confirming that TNF-α G-308A gene variants are affected by race [35] . Heterogeneity was significantly decreased and the TNF-α G-308A polymorphism was associated with OA 

TNF-α G-308A and OA risk in the Asian population when the analysis stratified by ethnicity was performed. This finding may be attributed to the AA genotype that is more important for knee OA susceptibility in Asians than in Caucasians. The TSA result showed that the current cumulative samples were sufficient to reach a decisive conclusion in Asians but not in Caucasians. This is believed to be the first meta-analysis study to explore the correlation between TNF-α G-308A and knee OA using TSA.

In the precision medicine era, genomic and environmental factors that are significantly related to disease phenotypes may be used as biomarkers for risk prediction and early disease diagnosis [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] . Well-trained deep learning models based on these biomarkers have been widely used as tools to identify diseases [47] . The next goal of the current study is to use deep learning approaches to predict OA or its progression upon genomic variants.

The current study has three major strengths: (1) none of the previous meta-analyses estimated whether the cumulative sample size is sufficient to obtain a definite conclusion. This study employed TSA estimation and found that definite conclusions can be obtained based on the Asian sample size. (2) This study used the random-effects model to combine results and avoid serious errors caused by model selection based on heterogeneity [48] . (3) In this study, two evaluators used the Newcastle-Ottawa Scale to evaluate the paper quality and excluded those that did not conform to Hardy-Weinberg equilibrium or had irrational data, thereby enhancing the reliability of the results. However, this study has some limitations. (1) The current study only included papers published in English, and those in other languages were not included in the meta-analysis. This omission may result in bias in the combined results. (2) The high heterogeneity could not be explained, which may imply potential gene-gene and gene-environment interactions. Further research is required to shed light on complete population characteristics for future meta-analysis. Moreover, the issues (e.g., publication bias from smaller-scale research or heterogeneity) can be assessed or processed by utilizing quality assessment scores [49] [50] [51] [52] [53] or stratification [54] [55] [56] .

(3) This study shows limited causality to explore the genetic variant effect on OA development. In the future, linkage disequilibrium score and Mendelian randomization should be employed to see if the TNF-α G-308A genetic polymorphism or other SNPs in this gene may causally trigger the development of OA through mediating the expression of this gene in cartilage [57] [58] [59] [60] [61] .

In conclusion, the definite null relationship between TNF-α G-308A and knee OA was validated in this case-control study. However, this meta-analysis demonstrated that the AA genotype of the TNF-α G-308A polymorphism increased the risk of OA in Asians. A decisive conclusion could be obtained, for which the current case-control samples provided critical evidence. In contrast, TNF-α G-308A polymorphism was not associated with knee OA risk in Caucasians. Further studies with a large sample size are necessary to validate whether TNF-α gene polymorphism contributes to knee OA susceptibility in Caucasians.

Supporting information S1 Table. TNF-α G-308A paper quality evaluation. (DOCX) S1 Fig. Estimation of Caucasian sample size for correlation between TNF-α G-308A and osteoarthritis. A trial sequential analysis was performed using a dominant model assumption. Detailed settings: significance level = 0.05, power = 0.8, least extreme odds ratio to be detected = 1.5, minor allele frequency = 0.069, and I 2 (heterogeneity) = 0%. 

Prevalence and development of hip and knee osteoarthritis according to American College of Rheumatology criteria in the CHECK cohort

Genetic influences on osteoarthritis in women: a twin study

The role of structural genes in the pathogenesis of osteoarthritic disorders

Action of bacterial polysaccharide on tumors. Iii. Repeated response of Sarcoma 37, in tolerant mice, to Serratia marcescens endotoxin

Mammalian MAP kinase signalling cascades

Tumor necrosis factor alpha-mediated inhibition of melanogenesis is dependent on nuclear factor kappa B activation

Treatment of TNF mediated diseases by selective inhibition of soluble TNF or TNFR1

Signal transduction by tumor necrosis factor and its relatives

Role of proinflammatory cytokines in the pathophysiology of osteoarthritis

Association analysis of two candidate polymorphisms in the tumour necrosis factor-alpha gene with osteoarthritis in a Chinese population

Meta-analysis of tumor necrosis factor alpha −308 polymorphism and knee osteoarthritis risk

Association of genotypes affecting the expression of interleukin-1beta or interleukin-1 receptor antagonist with osteoarthritis

Tumour necrosis factor alpha −308G/A gene polymorphism: lack of association with knee osteoarthritis in a Turkish population

TNF-alpha and TNF-beta polymorphisms are associated with susceptibility to osteoarthritis in a Korean population

High expression of TNF alpha is associated with −308 and −238 TNF alpha polymorphisms in knee osteoarthritis

Proinflammatory cytokine gene polymorphisms and the plasma and synovial fluid levels of cytokines in patients with post-traumatic knee osteoarthritis

The single-nucleotide polymorphism (SNP) of tumor necrosis factor alpha −-308G/A gene is associated with early-onset primary knee osteoarthritis in an Egyptian female population

Association between tumor necrosis factor alpha rs1800629 polymorphism and risk of osteoarthritis in a Chinese population

Association of TNF-alpha-308G>A and −238G>A polymorphisms with knee osteoarthritis risk: A case-control study and meta-analysis

Association of cytokine gene polymorphisms with osteoarthritis susceptibility

Functional status and severity of osteoarthritis in elderly is associated to the polymorphism of TNFA gene

Tumor necrosis factor-α: molecular assessment of gene expression, genetic variants and serum level in Egyptian patients with knee osteoarthritis

Trial sequential analysis reveals insufficient information size and potentially false positive results in many meta-analyses

Preferred Reporting Items for Systematic Reviews and Meta-Analyses: the PRISMA statement

Conducting meta-analyses in R with the metafor package

meta: Meta-Analysis with R

Copenhagen Trial Unit, Centre for Clinical Intervention Research

An allelic polymorphism within the human tumor necrosis factor alpha promoter region is strongly associated with HLA A1, B8, and DR3 alleles

Autocrine regulation of collagenase 3 (matrix metalloproteinase 13) during osteoarthritis

Interleukin-1beta and tumor necrosis factor-alpha increase stiffness and impair contractile function of articular chondrocytes

Transcriptional regulation of collagenase (MMP-1, MMP-13) genes in arthritis: integration of complex signaling pathways for the recruitment of gene-specific transcription factors

Gene expression (collagenase, tissue inhibitor of metalloproteinases, complement, and HLA-DR) in rheumatoid arthritis and osteoarthritis synovium. Quantitative analysis and effect of intraarticular corticosteroids

The role of synovitis in pathophysiology and clinical symptoms of osteoarthritis

The processing, gene regulation, biological functions, and clinical relevance of N4-acetylcytidine on RNA: A systematic review

Role of ethnic variations in TNF-α and TNF-β polymorphisms and risk of breast cancer in India

LEPR hypomethylation is significantly associated with gastric cancer in males

Immunodeficiency promotes adaptive alterations of host gut microbiome: an observational metagenomic study in mice

Genetic regulatory subnetworks and key regulating genes in rat hippocampus perturbed by prenatal malnutrition: implications for major brain disorders

Associations between maternal vitamin D status during three trimesters and cord blood 25(OH)D concentrations in newborns: a prospective Shanghai birth cohort study

Effects of early-life malnutrition on neurodevelopment and neuropsychiatric disorders and the potential mechanisms

The aberrantly expressed miR-193b-3p contributes to preeclampsia through regulating transforming growth factor-β signaling

Co-expression network analysis identified hub genes critical to triglyceride and free fatty acid metabolism as key regulators of age-related vascular dysfunction in mice

Association study of newly identified age-related macular degeneration susceptible loci SOD2, MBP, and C8orf42 in Han Chinese population

No evidence of association between variant rs2075650 in lipid metabolism-related locus APOE/TOMM40 and advanced age-related macular degeneration in Han Chinese population

DNA methylome profiling of maternal peripheral blood and placentas reveal potential fetal DNA markers for non-invasive prenatal testing

A splicing-regulatory polymorphism in DRD2 disrupts ZRANB2 binding, impairs cognitive functioning and increases risk for schizophrenia in six Han Chinese samples

A multi-model deep convolutional neural network for automatic hippocampus segmentation and classification in Alzheimer's disease

Recommendations for examining and interpreting funnel plot asymmetry in metaanalyses of randomised controlled trials

Genetic influences of dopamine transport gene on alcohol dependence: a pooled analysis of 13 studies with 2483 cases and 1753 controls

Association of BDNF and BCHE with Alzheimer's disease: Meta-analysis based on 56 genetic case-control studies of 12,563 cases and 12,622 controls

Meta-analyses of seven GIGYF2 polymorphisms with Parkinson's disease. Biomed Rep

Brain-derived neurotrophic factor gene C-270T and Val66Met functional polymorphisms and risk of schizophrenia: a moderate-scale population-based study and meta-analysis

Multi-trait analysis for genome-wide association study of five psychiatric disorders

A1166C genetic variation of the angiotensin II type I receptor gene and susceptibility to coronary heart disease: collaborative of 53 studies with 20,435 cases and 23,674 controls

Quantitative assessment of the effect of angiotensinogen gene polymorphisms on the risk of coronary heart disease. Circulation

Sex-specific Association of Circulating Ferritin Level and Risk of Type 2 diabetes: A dose-response meta-analysis of prospective studies

Causal influences of neuroticism on mental health and cardiovascular disease

Genetic evidence suggests posttraumatic stress disorder as a subtype of major depressive disorder

Smoking quantitatively increases risk for COVID-19

Genetic support of a causal relationship between iron status and type 2 diabetes: a Mendelian randomization study

Exploring the causal pathway from ischemic stroke to atrial fibrillation: a network Mendelian randomization study