Association of novel polymorphisms in TMEM39A gene with systemic lupus erythematosus in a Chinese Han population
© The Author(s). 2017
Received: 20 July 2016
Accepted: 13 April 2017
Published: 20 April 2017
This study aimed to assess the association between 14 single nucleotide polymorphisms (SNPs) in six genes (IRF8, TMEM39A, IKZF3, ORMDL3, GSDMB, and ZPBP2) and systemic lupus erythematosus (SLE) in a Chinese Han population sample.
We carried out a case-control study of 415 patients with SLE and 470 healthy controls without autoimmune disease or cancer. DNA for genetic analysis was isolated from the blood of all subjects using standard phenol-chloroform method. TagSNPs were identified using genotype data from the panel (Han Chinese in Beijing) of the HapMap Project and were selected using the Haploview program. Genotyping assay was conducted using the Sequenom MassARRAY iPLEX Gold platform. The frequencies of the alleles and genotypes were calculated and analyzed. Association studies and haplotype analysis were also performed.
The genotypic frequencies of rs12493175 and rs13062955 were significantly different between the SLE patients and the healthy controls. Compared with the common homozygous genotype, the CT and CT + TT genotypes in rs12493175 and the AC and AC + AA genotypes in rs13062955 was observed to significantly reduce the risk of SLE. The haplotype analysis of TMEM39A polymorphisms showed that the CGTA haplotype frequency was significantly low in the SLE patients.
Our findings identified three novel associations in SNPs located in the TMEM39A gene associated with SLE susceptibility in a Chinese Han population.
KeywordsSingle nucleotide polymorphism Systemic lupus erythematosus Susceptibility
Systemic lupus erythematosus (SLE) is typically characterized by the dysregulation of T cell response and B cell activation which usually causes the formation of immune complexes in multiple organs and tissues . Although the pathogenesis of SLE is largely unknown to date, it most likely involves environmental and genetic factors. Several candidate-gene studies and genome-wide association (GWA) scans have successfully discovered multiple susceptibility genes that fall into key pathways implicating immune complex clearance, immune signal transduction and interferon pathways contributing to the development of SLE [2, 3]. However, much of the heritable risk needs to be identified.
Recently, multiethnic approach was utilized to find that three SLE risk loci exceeded the genome-wide significance threshold, including interferon regulatory factor 8 (IRF8), transmembrane protein 39A (TMEM39A), and 17q21 between IKAROS family of zinc finger 3 (IKZF3) and zona pellucida binding protein 2 (ZPBP2) . The 17q21 region was originally associated with asthma in family linkage study [5–7]. More single nucleotide polymorphisms (SNPs) in the 17q21 region have been identified as being associated with the susceptibility to autoimmune diseases, including rheumatoid arthritis, ankylosing spondylitis and SLE [4, 8–10]. IRF8 is a family member of transcription factors that play a critical role in the regulation of cell apoptosis and immune response . It is required for promoting type I interferon responses which can induce the overexpression of genes reported in SLE, and several variants within IRF8 could influence binding to the regulatory elements [4, 12]. Although very limited biological data of TMEM39A is published so far, its polymorphisms have been found to be associated with multiple sclerosis and SLE [4, 13–15]. Additionally, several studies found genetic variants in orosomucoid like 3 (ORMDL3) and gasdermin B (GSDMB) were associated with the risk of autoimmune disease [16, 17]. On the basis of these studies, we hypothesized that certain novel variants in the loci described previously may contribute to the susceptibility to SLE.
In this study, we selected the six candidate genes, namely, IRF8, TMEM39A, IKZF3, ORMDL3, GSDMB, and ZPBP2, and screened for the putatively functional tagSNPs. We aimed to determine the association between the polymorphisms and susceptibility to SLE in a Chinese population.
General characteristics of the study population
SLE patients (n = 415)
Healthy controls (n = 470)
P > 0.05
Age (mean ± SD), (years)
38.2 ± 11.4
35.8 ± 13.6
P > 0.05
Fever, n (%)
Baldness, n (%)
Light sensitivity, n (%)
Facial erythema, n (%)
Oral ulcer, n (%)
Arthritis, n (%)
Lupus nephritis, n (%)
TagSNPs are representative SNPs which can capture most of the genetic variation in a region of the genome on the basis that they are in high linkage disequilibrium (LD) with other SNPs . TagSNPs genotyped in this study were selected by analyzing the genotype data of Chinese Han population from HapMap dbSNP (http://www.hapmap.ncbi.nlm.nih.gov) using LD-based tagSNP selection with a pairwise algorithm LDSelect, available in the Tagger function implemented in Haploview version 4.2 (http://www.broadinstitute.org/mpg/haploview) [20, 21].
Common SNPs captured using the selected 14 tagSNPs in the six candidate genes based on the HapMap population data for Chinese in Beijing (Release 27)
rs9635726, rs3816470, rs9303277, rs10445308, rs9909593
rs11078927, rs1008723, rs4795400, rs9303281, rs2305480, rs7219923, rs869402, rs2305479, rs7224129, rs11078926, rs2290400, rs7216389
rs4795403, rs4795402, rs3744246, rs4795404
rs4378650, rs8076131, rs12603332
rs11557467, rs12936231, rs11078925, rs1054609, rs11557466, rs11870965, rs10852936, rs9907088
rs188602, rs170033, rs2270502, rs381139
rs4843860, rs12926854, rs4843861
rs12492859, rs13094625, rs13081197, rs13078312, rs12493326, rs16829853, rs13081067, rs2282170, rs13062955, rs12492315, rs12493175, rs13096213, rs12496277, rs12492609
rs7629750, rs2282171, rs3772136, rs9846088, rs4687859, rs9872589, rs3195852
rs17281647, rs2282175, rs1132202
Details of the primers used in the polymorphism genotyping by MassArray
Forward and reverse primer
Reverse transcriptase–PCR of candidate gene mRNA levels
To examine the relation between the associated polymorphisms and the gene mRNA levels, forty patients stratified by polymorphic genotypes were randomly selected. The relative expression levels of twenty patients with common homogenous genotype carriers were set to a unity, and the relative expression levels of twenty patients with heterogeneous and rare homogenous genotypes were expressed relative to those of the common homogenous genotype carriers. Total RNA prepared from peripheral blood mononuclear cells were reversely transcribed using the TaqMan reverse transcription reagents (Applied Biosystems, Weiterstadt, Germany). RT-PCR was carried out using the ABI Universal Master Mix on an ABI PRISM 7000 Sequence Detection System.
Data were managed and stored using the SPSS software 16.0. Allele and genotype frequencies were compared between patient and control groups by the chi-square (χ 2) test. The quality of the genotype data was evaluated by Hardy-Weinberg equilibrium (HWE) in the case and control subjects using Fisher’s exact test (P > 0.05). The association between each polymorphism and risk of SLE was estimated by logistic regression and was expressed as odds ratio (OR) with 95% confidence intervals (95% CI). The haplotypes were assigned using the online software platform SHEsis (http://www.analysis.bio-x.cn). The haplotype construction element is based on the standard Full-Precise-Iteration (FPI) algorithm . All tests were two-tailed, and P values < 0.05 were considered as statistically significant.
The details of the identified genetic single-nucleotide variants
CADD = 3.564
CADD = 2.429
CADD = 8.630
CADD = 4.988
CADD = 0.720
The five SNPs call rates in patients and control individuals and HWE p-values
Call rate (%)
Genotype and allele association analysis of three tagSNPs
OR (95% CI)
2.95 × 10−4
The genotypic frequencies of rs12493175 and rs13062955 located in TMEM39A gene were significantly different between the SLE patients and the healthy controls. Compared with the common homozygous genotype, the CT and CT + TT genotypes in rs12493175 (p.adjust = 0.005, odds ratio (OR) 0.58, 95% CI 0.42 to 0.79; p.adjust = 0.007, OR 0.62, 95% CI 0.46 to 0.84, respectively) and the AC and AC + AA genotypes in rs13062955 (p.adjust = 0.002, OR 0.55, 95% CI 0.40 to 0.76; p.adjust = 0.007, OR 0.61, 95% CI 0.45 to 0.84, respectively) was observed to significantly reduce the risk of SLE. On the other hand, the difference in the frequency of rs2282175 was only marginal. The CT and CT + TT genotypes in rs2282175 was observed to modestly increase the risk of SLE (p.adjust = 0.054, OR 1.63, 95% CI 1.06 to 2.38; p.adjust = 0.054, OR 1.56, 95% CI 1.05 to 2.27, respectively). However, we did not find any other tagSNP associated with SLE risk in the genes of IRF8, IKZF3, ORMDL3, GSDMB and ZPBP2. We also evaluated the relation between the associated polymorphisms and the gene mRNA levels in peripheral blood mononuclear cells from 40 patients. Nevertheless, we failed to find any correlation between them (data not shown).
Frequencies of the haplotypes formed by TMEM39A rs2282175, rs12493175 and rs13062955 SNPs
OR (95% CI)
IRF8, TMEM39A and IKZF3-ZPBP2 were previously identified as susceptibility loci for SLE in the multiracial replication study , Besides, ORMDL3 and GSDMB were found to have susceptibility loci for autoimmune diseases [16, 17]. Thus we hypothesized certain novel associations in SNPs located in these genes could be identified in Chinese populations. To test this hypothesis, we selected 14 tagSNPs in these candidate genes to determine the association between the polymorphisms and SLE susceptibility in a Chinese Han population. Our findings showed that TMEM39A rs2282175, rs12493175, and rs13062955 were associated with SLE risk.
To date, almost no biological data of TMEM39A have been reported and only two SNPs in TMEM39A were identified as being associated with the susceptibility of autoimmune diseases. TMEM39A rs1132200 have been found to be associated with susceptibilities to multiple sclerosis and SLE in multiracial replication study [4, 13, 14]. but the recent studies showed that TMEM39A rs12494314, instead of rs1132200, was associated with SLE susceptibility in the Chinese population [15, 24]. In our current study, we identified three novel associations in SNPs located in TMEM39A as being associated with SLE susceptibility. The genotypic frequencies of rs12493175 and rs13062955 were significantly different between the SLE patients and the healthy controls, while the difference in the frequency of rs2282175 was only marginal. Among these polymorphisms, rs12493175 T-allele and rs13062955 A-allele were found to be associated with a reduced SLE risk, suggesting a protective factor to SLE. In contrast, rs2282175 T-allele was found to be associated with an increased SLE risk, suggesting a susceptibility factor to SLE. Haplotype analysis for TMEM39A SNPs revealed that the haplotype CGTA conferred a reduced risk of SLE. It is possible that the haplotype CGTA provides protection to SLE, resulting from the rs12493175 T and rs13062955 A alleles. It is worth noting that rs2282175 is located in the region of 5' upstream in TMEM39A and predicted to be a binding site of certain transcription factor. It is speculated that the C → T allele change at the rs2282175 site may influence the DNA binding ability of transcription factor c-Rel, which was predicted according to the different variants using the search tool of Alibaba 2.1 (http://www.gene-regulation.com/pub/programs/alibaba2). Although we did not find any relation between mRNA expression and the polymorphisms, it may be required to explore the possible biological significance of the SNPs in different cell subsets.
Several limitations in the current study should also be noted. First, due to the restricted number of study subjects and limited analysis capacity, we did not analyze the SNPs with rare MAF in the Chinese Han population, including those reported as risk loci for SLE in other ancestries, and mainly focused on the SNPs with predicted functional effect. Further research on the role of the rare variants should be carried out in a larger number of samples. As different populations have different genetic backgrounds, it is still necessary to perform the genetic analysis of multiracial study. Second, we still could not determine the causality of SLE-associated SNPs. For those variants in the large strong LD region, such as chromosome 17q21, it is difficult to determine which SNP is the true functional locus that contributes to SLE susceptibility independently. Better understanding whether the SNP is functionally relevant will require mechanistic and fine-mapping experiments. Third, our study did not assess the SNP-SNP interaction (epistasis). For the genetically complex disease, multiple interacting loci could contribute to SLE susceptibility. Additionally, as a heterogenetic disease, the contribution of genetic and environmental factors is very important for the disease . Therefore, interaction between genetic and environmental factors is required to further clarify the pathogenesis of SLE.
This study identified three novel associations in SNPs located in the TMEM39A gene associated with SLE susceptibility in a Chinese Han population. Functional study and further independent large-scale study in other racial populations are still needed to confirm our results.
IKAROS family of zinc finger 3
Interferon regulatory factor 8
Minor allele frequency
Orosomucoid like 3
Systemic lupus erythematosus
Single nucleotide polymorphism
Transmembrane protein 39A
Zona pellucida binding protein 2
We especially thank all SLE patients who participated for making this study possible.
Availability of data and material
All data generated or analysed during this study are included in this published article.
XC contributed to the collection, analysis of the clinical data, and manuscript preparation. WH and XL contributed to the analysis of the SNP results. LW and YJ coordinated the study and helped draft the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The study was approved by the Human Ethics Committee of China Medical University and the informed consents were obtained from all donors.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Rother N, van der Vlag J. Disturbed T cell signaling and altered Th17 and regulatory T cell subsets in the pathogenesis of systemic lupus erythematosus. Front Immunol. 2015;6:610.View ArticlePubMedPubMed CentralGoogle Scholar
- Moser KL, Kelly JA, Lessard CJ, Harley JB. Recent insights into the genetic basis of systemic lupus erythematosus. Genes Immun. 2009;10(5):373–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Deng Y, Tsao BP. Genetic susceptibility to systemic lupus erythematosus in the genomic era. Nat Rev Rheumatol. 2010;6(12):683–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Lessard CJ, Adrianto I, Ice JA, Wiley GB, Kelly JA, Glenn SB, Adler AJ, Li H, Rasmussen A, Williams AH, et al. Identification of IRF8, TMEM39A, and IKZF3-ZPBP2 as susceptibility loci for systemic lupus erythematosus in a large-scale multiracial replication study. Am J Hum Genet. 2012;90(4):648–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Dizier MH, Besse-Schmittler C, Guilloud-Bataille M, Annesi-Maesano I, Boussaha M, Bousquet J, Charpin D, Degioanni A, Gormand F, Grimfeld A, et al. Genome screen for asthma and related phenotypes in the French EGEA study. Am J Respir Crit Care Med. 2000;162(5):1812–8.View ArticlePubMedGoogle Scholar
- Marinho S, Custovic A, Marsden P, Smith JA, Simpson A. 17q12-21 variants are associated with asthma and interact with active smoking in an adult population from the United Kingdom. Ann Allergy Asthma Immunol. 2012;108(6):402–11. e409.View ArticlePubMedGoogle Scholar
- Fang Q, Zhao H, Wang A, Gong Y, Liu Q. Association of genetic variants in chromosome 17q21 and adult-onset asthma in a Chinese Han population. BMC Med Genet. 2011;12:133.View ArticlePubMedPubMed CentralGoogle Scholar
- Qiu R, Zhang H, Zhao H, Li J, Guo C, Gong Y, Liu Q. Genetic variants on 17q21 are associated with ankylosing spondylitis susceptibility and severity in a Chinese Han population. Scand J Rheumatol. 2013;42(6):469–72.View ArticlePubMedGoogle Scholar
- Kurreeman FA, Stahl EA, Okada Y, Liao K, Diogo D, Raychaudhuri S, Freudenberg J, Kochi Y, Patsopoulos NA, Gupta N, et al. Use of a multiethnic approach to identify rheumatoid- arthritis-susceptibility loci, 1p36 and 17q12. Am J Hum Genet. 2012;90(3):524–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Cai X, Qiao Y, Diao C, Xu X, Chen Y, Du S, Liu X, Liu N, Yu S, Chen D, et al. Association between polymorphisms of the IKZF3 gene and systemic lupus erythematosus in a Chinese Han population. PLoS One. 2014;9(10), e108661.View ArticlePubMedPubMed CentralGoogle Scholar
- Tamura T, Yanai H, Savitsky D, Taniguchi T. The IRF family transcription factors in immunity and oncogenesis. Annu Rev Immunol. 2008;26:535–84.View ArticlePubMedGoogle Scholar
- Baechler EC, Batliwalla FM, Karypis G, Gaffney PM, Ortmann WA, Espe KJ, Shark KB, Grande WJ, Hughes KM, Kapur V, et al. Interferon-inducible gene expression signature in peripheral blood cells of patients with severe lupus. Proc Natl Acad Sci U S A. 2003;100(5):2610–5.View ArticlePubMedPubMed CentralGoogle Scholar
- (IMSGC) IMSGC. Comprehensive follow-up of the first genome-wide association study of multiple sclerosis identifies KIF21B and TMEM39A as susceptibility loci. Hum Mol Genet. 2010;19(5):953–62.View ArticleGoogle Scholar
- Varade J, Comabella M, Ortiz MA, Arroyo R, Fernandez O, Pinto-Medel MJ, Fedetz M, Izquierdo G, Lucas M, Gomez CL, et al. Replication study of 10 genes showing evidence for association with multiple sclerosis: validation of TMEM39A, IL12B and CBLB [correction of CLBL] genes. Mult Scler. 2012;18(7):959–65.View ArticlePubMedGoogle Scholar
- Sheng YJ, Xu JH, Wu YG, Zuo XB, Gao JP, Lin Y, Zhu ZW, Wen LL, Yang C, Liu L, et al. Association analyses confirm five susceptibility loci for systemic lupus erythematosus in the Han Chinese population. Arthritis Res Ther. 2015;17:85.View ArticlePubMedPubMed CentralGoogle Scholar
- Verlaan DJ, Berlivet S, Hunninghake GM, Madore A-M, Larivière M, Moussette S, Grundberg E, Kwan T, Ouimet M, Ge B. Allele-specific chromatin remodeling in the ZPBP2/GSDMB/ORMDL3 locus associated with the risk of asthma and autoimmune disease. Am J Hum Genet. 2009;85(3):377–93.View ArticlePubMedPubMed CentralGoogle Scholar
- Li X, Ampleford EJ, Howard TD, Moore WC, Torgerson DG, Li H, Busse WW, Castro M, Erzurum SC, Israel E. Genome-wide association studies of asthma indicate opposite immunopathogenesis direction from autoimmune diseases. J Allergy Clin Immunol. 2012;130(4):861–868. e867.View ArticlePubMedPubMed CentralGoogle Scholar
- Yu C, Gershwin ME, Chang C. Diagnostic criteria for systemic lupus erythematosus: a critical review. J Autoimmun. 2014;48–49:10–3.View ArticlePubMedGoogle Scholar
- Kathiresan S, Gabriel SB, Yang Q, Lochner AL, Larson MG, Levy D, Tofler GH, Hirschhorn JN, O’Donnell CJ. Comprehensive survey of common genetic variation at the plasminogen activator inhibitor-1 locus and relations to circulating plasminogen activator inhibitor-1 levels. Circulation. 2005;112(12):1728–35.View ArticlePubMedGoogle Scholar
- Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21(2):263–5.View ArticlePubMedGoogle Scholar
- Carlson CS, Eberle MA, Rieder MJ, Yi Q, Kruglyak L, Nickerson DA. Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium. Am J Hum Genet. 2004;74(1):106–20.View ArticlePubMedGoogle Scholar
- Yi X, Li J, Yu S, Zhang A, Xu J, Yi J, Zou J, Nie X, Huang J, Wang J. A new PCR-based mass spectrometry system for high-risk HPV, part I: methods. Am J Clin Pathol. 2011;136(6):913–9.View ArticlePubMedGoogle Scholar
- Yong Y, Lin H. SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Res. 2005;15(2):97–8.View ArticleGoogle Scholar
- You Y, Zhai ZF, Chen FR, Chen W, Hao F. Autoimmune risk loci of IL12RB2, IKZF1, XKR6, TMEM39A and CSK in Chinese patients with systemic lupus erythematosus. Tissue Antigens. 2015;85(3):200–3.View ArticlePubMedGoogle Scholar
- Sparks JA, Costenbader KH. Genetics, environment, and gene-environment interactions in the development of systemic rheumatic diseases. Rheum Dis Clin North Am. 2014;40(4):637–57.View ArticlePubMedPubMed CentralGoogle Scholar