Abstract
Aims:
The toll-like receptor (TLR) genes were shown to be involved in the pathogenesis of rheumatoid arthritis (RA). We aimed to investigate the genetic associations between the TLR-1, -2, -4, and -6 genes polymorphisms with RA susceptibility in a Chinese Han population.
Methods:
Six polymorphisms [TLR-1 (rs5743610, rs5743618), -2 (rs5743708), -4 (rs4986790, rs4986791), and -6 (rs5743810)] in TLRs genes were genotyped in 360 patients with RA and 560 matched healthy controls by direct sequencing. The odds ratios (ORs) and 95% confidence intervals (CIs) were evaluated using a standard logistic regression analysis.
Results:
No significant associations were observed under the allelic, dominant, or recessive models for TLR-1 rs5743610, TLR-2 rs5743708, TLR-4 rs4986790 and rs4986791, and TLR-6 rs5743810 polymorphisms and RA risk (all p > 0.05). However, significant associations were detected under the allelic, dominant, and recessive models for TLR-1 rs5743618 and RA risk (allelic: OR [95% CI] = 2.21 [1.73-2.81], p < 0.0001; dominant: OR [95% CI] = 2.33 [1.75-3.09], p < 0.0001; recessive models: OR [95% CI] = 3.70 [1.85-7.41], p = 0.0002). In addition, TLR6 rs5743810 was found to be associated with the rheumatoid factor (RF)− and anticyclic citrullinated peptide (anti-CCP)− antibody in RA group (RF: OR [95% CI] = 2.29 [1.42-3.69], p = 0.0007; anti-CCP: OR [95% CI] = 2.33 [1.39-3.89], p = 0.001).
Conclusions:
The allelic, dominant, and recessive models of TLR1 rs5743618 might be associated with RA susceptibility. Also, the TLR6 rs5743810 marker may be associated with RF and the anti-CCP antibody of RA in the Chinese Han population.
Introduction
Rheumatoid arthritis (RA) is a chronic, aggressive, and systemic autoimmune disease (Gallacher, 2015). The incidence of RA varies by region and race, and the prevalence of RA in China is about 0.3-0.5%. Also, the RA prevalence in females is three times than that in males (Yu et al., 2018). Epidemiological studies have found that the incidence of RA is characterized by familial aggregation and co-twin disease, and the heritability of RA in identical twins is 53-65% (Alarcón-Segovia et al., 2005; Kuo et al., 2017). Although the etiology of RA is undetermined, a number of genetic components have been established (Karlson et al., 2013).
Toll-like receptors (TLRs) are the specific type I transmembrane receptors and pathogenic pattern recognition receptors in the natural immune system, and play an important role in inflammation, cell signal transduction, and apoptosis (Netea et al., 2006; Ospelt and Gay, 2010). It has been found that TLRs, especially TLR2 and TLR4, on the cell surface play a key role in the development of RA (Ruckdeschel et al., 2004). Studies showed that the expressions of TLR-2, -3, -4, and -7 in synovial tissues, synovial fibroblasts, peripheral blood monocytes, and CD14+ macrophages in synovial fluid of RA patients significantly increased (Iwahashi et al., 2004; Thwaites et al., 2014; Huang et al., 2007).
In addition, according to Spitzer et al., TLRs could interact and modulate with each other. TLR1 and TLR4 can form heterodimer complexes. Their extracellular segments can bind to TLR4 and block the connection between TLR4 and ligand, thus blocking the signal transduction of TLR4 (Spitzer et al., 2002). Furthermore, Bulut et al. (2001) suggested that the synergistic effect of TLR2 and TLR6 may be involved in the identification of all TLR2 ligands. Further experiments also confirmed that the fused TLR2 and TLR6 could inhibit the activation of the ligands of TLR2 and TLR4 (Spitzer et al., 2002).
One of the possible mechanisms driving TLRs dysfunction in RA is genetic variants in genes encoding these receptors. Recently, increasing number of studies has been carried out on the TLRs gene polymorphisms and the susceptibility of RA. Most of studies have focused on the genetic association between TLR2 and TLR4 genes polymorphisms and RA susceptibility. Two well-known missense mutations “Asp299Gly (rs4986790)” and “Thr399Ile (rs4986791)” in TLR4 gene have been shown to be correlated with RA risk. However, such associations were not shown in all the investigated populations (Kilding et al., 2003; Radstake et al., 2004; Zheng et al., 2010). With regard to TLR2 Arg753Gln (rs5743708), it has been shown not to be the risk factor for arthritis in a Spanish population (Sanchez et al., 2004). As for TLR1 and TLR6, no correlation of TLR1 Arg80Thr (rs5743610), Ser602Ile (rs5743618), and TLR6 Ser249Pro (rs5743810) polymorphisms in RA was observed (Jaen et al., 2009).
However, evaluation of the associations between the TLRs, especially the TLR1, TLR2, and TLR6 gene polymorphisms and RA susceptibility in Chinese Han population has been very limited. Thus, by utilizing a case-control method, we aim to investigate whether polymorphisms of the TLR-1, -2, -4, and -6 gene polymorphisms contribute to the development of RA in a Chinese Han population.
Materials and Methods
Study subjects
The experimental protocol was evaluated and approved by the Ethics Committee of The First Affiliated Hospital of Soochow University (China) (EC/19021, December 10, 2019). Written informed consent for genetic analysis was obtained from all individuals. Three hundred sixty patients who fulfill the American College of Rheumatology 1987 criteria for RA were enrolled as case group in the present study. A control group consisted of 560 healthy individuals, who carried out physical examination from the First Affiliated Hospital of Soochow University, were also included. The patients and controls were Chinese Han origin and matched for age and sex. Detailed information of patients and controls was concluded in Table 1.
Clinical Characteristics of Rheumatoid Arthritis Patients and Healthy Controls
Disease activity score 28 (DAS28): a score for evaluation of RA activity by assessing the state of 28 joints; anti-CCP: anticyclic citrullinated peptide.
CRP, C-reactive protein; ESR, erythrocyte sedimentation rate; RA, rheumatoid arthritis; RF, rheumatoid factor; SD, standard deviation;.
Genotyping
Genomic DNA was extracted by standard phenol chloroform method. Genotyping detection was performed on six polymorphisms of TLR-1 (rs5743610, rs5743618), TLR-2 (rs5743708), TLR-4 (rs4986790, rs4986791), and TLR-6 (rs5743810) genes by direct sequencing.
Statistical analysis
The SPSS 20.0 were applied for statistical analysis of the data. The Hardy-Weinberg equilibrium (HWE) of the single nucleotide polymorphisms (SNPs) in control group was tested with χ2. The chi-squared test was used for the difference of genotype and allele frequency between control and RA groups. The corresponding odds ratios (ORs) and 95% confidence intervals (CIs) of RA risk of individuals with genetic models, including allelic, dominant, and recessive models of selected SNPs, were calculated using logistic regression analysis.
In addition, the genetic association between the alleles of selected SNPs and autoantibody status (rheumatoid factor [RF] and anticyclic citrullinated peptide [anti-CCP]) were also performed. Multiple comparison correction was performed to adjust p values using Benja Mini-Hochberg (BH) method based on False Discovery Rate (FDR) standard. Statistical significance was set at p < 0.05. Calculation power was obtained at the 0.05 level of significance, assuming an OR of 1.5 (small effect size) by using the G*Power software (www.gpower.hhu.de).
Results
A total of 360 patients with RA and 560 healthy controls were included in the present study. The genotype distribution was in HWE in the control group for all the SNPs (p > 0.05) (Table 2). Moreover, our study obtained >74.5% power at the 5% significant level (two-tailed) (Table 2).
Distribution of the Genotypes of TLR1,2,4,6 Gene Polymorphisms in Rheumatoid Arthritis Cases and Controls
p Value were calculated using Fisher's exact test.
The FDR adjusting was carried out to adjust the p value.
“—”: not significant; padj, padjusted.
95% CI, 95% confidence interval; FDR, False Discovery Rate; HWE, Hardy-Weinberg equilibrium; OR, odds ratio; SNP, single nucleotide polymorphism; TLR, toll-like receptor.
The associations between TLRs genes polymorphisms and RA risk
Significant associations were detected between the recessive model of the TLR1 rs5743610 (OR [95% CI] = 1.60 [1.04-2.46], p = 0.03), TLR6 rs5743810 (OR [95% CI] = 2.21 [1.15-4.27], p = 0.02), and RA risk. However, the significant associations disappeared after the FDR correction (rs5743610: padj = 0.12; rs5743810: padj = 0.06) (Table 2).
Significant difference was found between the distributions of the allelic and dominant models of the TLR2 rs5743708 (allelic model: OR [95% CI] = 1.38 [1.05-1.81], p = 0.02; dominant model: OR [95% CI] = 1.39 [1.02-1.87], p = 0.03) in RA and control groups. However, these significant associations did not exist after the FDR correction (p > 0.05) (Table 2).
In addition, the frequencies of the allelic (OR [95% CI] = 2.21 [1.73-2.81], p < 0.0001), dominant (OR [95% CI] = 2.33 [1.75-3.09], p < 0.0001), and recessive models (OR [95% CI] = 3.70 [1.85-7.41], p = 0.0002) of TLR1 rs5743618 were significantly different in RA group from that in the control group, even after the FDR correction (p < 0.05).
Furthermore, no TLR4 rs4986790 and rs4986791 mutations were detected in both RA patients and control subjects (Table 2).
The associations between risk allele in TLRs genes polymorphisms and clinical characteristics
Correlation between the TLR-1, -2, -4, and -6 risk alleles and the status of RF and anti-CCP in patients with RA was evaluated by regression analysis.
As shown in Table 3, 360 patients with RF and anti-CCP information available were divided into two groups (RF status: RF-positive [RF+] RA group and RF-negative [RF−] RA group; anti-CCP status: anti-CCP-positive [anti-CCP+] RA group and anti-CCP-negative [anti-CCP−] RA group). Results revealed that no significant association was detected between TLR1 rs5743610, rs5743618, TLR2 rs5743708, TLR4 rs4986790, rs4986791, and RF status, as well as anti-CCP status (p > 0.05). Moreover, the distribution of the TLR6 rs5743810 in RF+ (OR [95% CI] = 2.29 [1.42-3.69], p = 0.0007) and anti-CCP+ (OR [95% CI] = 2.33 [1.39-3.89], p = 0.001) RA groups were significantly different from those in RF− and anti-CCP− RA groups. These significant associations existed even after the FDR correction (p < 0.05) (Table 3).
TLR1,2,4,6 Gene Alleles Frequencies and Autoantibody Levels in Patients with Rheumatoid Arthritis
p Value were calculated using Fisher's exact test.
—, not calculated; .
Discussion
TLRs play an important role in both innate and adaptive immune responses (Majewska and Szczepanik, 2006). In innate immunity, TLRs directly trigger the intracellular bactericidal mechanism or induce the production of immune inflammatory factors based on the recognition of pathogenic microorganisms to expand the nonspecific defense (Cho et al., 2017). In adaptive immune response, TLRs can induce dendritic cells to mature, activate them to secrete cytokines and chemokines, and induce their expression of costimulatory molecules (Sun et al., 2018). In addition, TLRs could activate naive B cells, affect the intensity and quality of the memory T cell response, and initiate the CD8+T cell response against soluble protein antigen (Massonnet et al., 2009).
Thus, TLRs are widely involved in various pathways of specific immune response and in the occurrence and development of chronic immune inflammatory diseases. Furthermore, increased expression of TLRs was found in peripheral blood mononuclear cells, synovial tissues, synovial fluid, synovial macrophages, and synovial fibroblasts in RA patients. TLR-2, -3, -4, -5, -7, -8, and -9 were also demonstrated to involve in the development of RA (Ospelt et al., 2008; Sacre et al., 2016).
TLR2 and TLR4 are the most studied molecules in the TLR family. The TLR2 gene is located in the 4q31.3 region and works in concert with at least two other heterodimers such as TLR1 and TLR6. Previous studies have found that the positive rate and average fluorescence intensity of TLR2 in peripheral blood mononuclear cells in patients with RA were significantly higher than those in healthy controls suggesting that TLR2 plays an important role in the occurrence and development of RA (Yang et al., 2011).
However, little is known about the relationship of TLR2 rs5743708 and the susceptibility of RA. TLR2 rs5743708 is located at amino acid position 753 of intron. The A/G mutation causes arginine to become glutamine, which is a nonsynonymous mutation. In the present study, we found no association between the TLR2 rs5743708 and RA risk, as well as the RF and anti-CCP status, which was similar with the results reported by Sanchez et al. (2004) in a Spanish population.
The TLR4 gene is located in the 4q33.1 region. Studies on TLR4 gene polymorphisms and RA are currently focused on two missense mutation sites “Asp299Gly” and “Thr399Ile.” The Asp299Gly (TLR4/896) polymorphism is an amino acid variation, in which aspartic acid at 299 of TLR4 protein is replaced by glycine (Asp/Gly) due to missense mutation of “A/G” at 896th position of TLR4 gene. The Thr399Ile (TLR4/1196) refers to the substitution of cytosine (C) at position 1196 for thymine (T), resulting in substitution of threonine at position 399 by isoleucine (Thr/Ile). Lines of studies have reported the association between TLR4 rs4986790 and rs4986791 and RA in multiple populations. However, the results were inconsistent in different ethnicity groups. Radstake et al. found that TLR4 rs4986790 was associated with RA susceptibility, but not with RA disease activity and prognosis in a Dutch population.
However, no association was found between the TLR4 rs4986790 and rs4986791 and RA in other Caucasian populations, including French (Jaen et al., 2009), Finn (Kuuliala et al., 2006), Spanish (Sanchez et al., 2004), and British (Kilding et al., 2003). The results of the relationship of TLR4 rs4986790 and rs4986791 and RA in Asian populations were more complex. In the present study, both the TLR4 rs4986790 and rs4986791 were found to be not polymorphismic, which was similar with the previous results reported by Zheng et al. (2010), Yuan et al. (2010), and Kang and Lee (2007).
In addition, no significant association between the TLR4 rs4986790 and rs4986791 polymorphisms and RA risk was identified by meta-analysis (Lee and Song, 2013; Tizaoui et al., 2015). Thus, we may conclude that the TLR2 rs5743708, TLR4 rs4986790, and rs4986791 did not contribute to the pathogenesis of RA. These results need to be confirmed in larger number of case-control studies.
TLR1 and TLR6 are the receptors of TLR2, which could induce different transduction pathways (Elass et al., 2005; Massari et al., 2006). TLR6, TLR1 and TLR10 gene are clustered in a 54-kb region on chromosome 4p14 and encode proteins that share a high degree of homology in their overall amino acid sequences (Mikacenic et al., 2013). Genetic variations within TLR6-TLR1-TLR10 gene cluster have been reported to be associated with sarcoidosis (Veltkamp et al., 2011), prostate cancer (Stevens et al., 2010), Meniere's disease (Requena et al., 2013), Legionnaires' disease (Misch et al., 2013), and Crohn's disease (Morgan et al., 2012).
Eleven TLR1 polymorphisms have been reported in a population-based study in Sweden, and three were found to have an association with increased risk of prostate cancer (Chen et al., 2007). A 5′-UTR SNP of TLR1 gene, rs5743565, was demonstrated to have a protective effect against Graves' disease (Xiao et al., 2014). For other two TLR1 polymorphisms, rs5743610 and rs5743618, the TLR1 rs5743618 is a common nonsynonymous SNP lying just at the junction of the transmembrane and cytoplasmic domain of TLR1, resulting in a reduced response to heterodimer antagonists. In addition, the TLR1 rs5743610 and rs5743618 were reported to be associated with complicated skin and skin structure infections (cSSSIs) (Stappers et al., 2014) and Crohn's disease (Marie et al., 2006).
However, rare study has been conducted on the genetic associations between the TLR1 gene polymorphisms and RA risk. To our knowledge, this is the first study of the genetic association between the allelic, dominant, and recessive models of TLR1 rs5743610 and rs5743618 and RA susceptibility was conducted and demonstrated that the TLR1 rs5743618 was a susceptible factor for RA. Further studies with larger number of subjects are needed to investigate the genetic association between TLR1 rs5743618 and RA risk for the relatively small sample size in the present study.
Several TLR6 SNPs have been reported to be related with multiple diseases, including tuberculosis (Zhang et al., 2013), mild malaria (Leoratti et al., 2008), cSSSIs (Stappers et al., 2014), asthma (Hoffjan et al., 2005), and invasive aspergillosis (Kesh et al., 2005). The TLR6 rs5743810 T allele was found to reduce nuclear factor kappa B signaling which led to an altered level of IL-6 production (Shey et al., 2010). However, no study was conducted on the genetic association between TLR6 polymorphisms and RA. Notable, significant associations were first observed between the TLR6 rs5743810 and RA status, as well as anti-CCP status, which was different from the results reported by Jaen et al. (2009). The inconsistency may be due to the different genetic background between the two populations. The function of this SNP is still hardly known. We hypothesized that TLR6 rs5743810 might be associated with disease severity by increasing autoantibody production. To confirm this result, functional study is necessary in the future.
In summary, our results suggested that the allelic, dominant, and recessive models of TLR1 rs5743618 might be associated with the pathogenesis of RA. The TLR6 rs5743810 might be associated with clinical characteristics (RF and anti-CCP antibody) in RA in Chinese Han population. Further studies of the TLR1 rs5743618 and TLR6 rs5743810 in ethnically well-defined populations as well as functional studies of this variation are warranted to evaluate their role in the pathogenesis of RA.
Footnotes
Author Disclosure Statement
No competing financial interests exist.
Funding Information
This study was funded by the Foundation of Changshu department of Science and Technology (2017BY28020).
