Journal of Immunology Research
Volume 2019 (2019), Article ID 2397698, 10 pages
https://doi.org/10.1155/2019/2397698

Research Article

Association of Melatonin Pathway Gene’s Single-Nucleotide Polymorphisms with Systemic Lupus Erythematosus in a Chinese Population

Peng Wang,1,2 Lei Liu,³ Li-Fang Zhao,2,4 Chan-Na Zhao,2,4 Yan-Mei Mao,2,4 Yi-Lin Dan,2,4 Qian Wu,2,4 Xiao-Mei Li,⁵ De-Guang Wang,⁶ and Hai-Feng Pan2,4

¹Center for Genetic Epidemiology and Genomics, School of Public Health, Medical College of Soochow University, 199 Renai Road, Suzhou, Jiangsu, China
²Department of Epidemiology and Biostatistics, School of Public Health, Anhui Medical University, 81 Meishan Road, Hefei, Anhui, China
³Department of Parasitology, Medical College of Soochow University, 199 Renai Road, Suzhou, Jiangsu, China
⁴Anhui Province Key Laboratory of Major Autoimmune Diseases, 81 Meishan Road, Hefei, Anhui, China
⁵Department of Rheumatology and Immunology, Anhui Provincial Hospital, 17 Lujiang Road, Hefei, Anhui, China
⁶Department of Nephrology, The Second Affiliated Hospital of Anhui Medical University, Hefei, Anhui, China

Correspondence should be addressed to Hai-Feng Pan; panhaifeng@ahmu.edu.cn

Received 29 July 2019; Revised 21 October 2019; Accepted 1 November 2019; Published 13 November 2019

Guest Editor: Minggang Zhang

Copyright © 2019 Peng Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Objectives. This study was to investigate the association of melatonin (MTN) pathway gene’s single-nucleotide polymorphisms (SNPs) with susceptibility to systemic lupus erythematosus (SLE). Methods. We recruited 495 SLE patients and 493 healthy controls, 11 tag SNPs in MTN receptor 1a (MTNR1a), MTNR1b, and arylalkylamine N-acetyltransferase (AANAT) genes were genotyped and analyzed. Serum MTN concentration was determined by enzyme-linked immunosorbent assay (ELISA) kits. Results. Two SNPs of AANAT gene (rs8150 and rs3760138) associated with the risk of SLE; CC carriers of rs8150 had a lower risk as compared to GG (, 95% CI: 0.361, 0.799), whereas GG carrier in rs3760138 had an increased risk (, 95% CI: 1.154, 2.880) compared to TT. However, we did not find any genetic association between the other nine SNPs with SLE risk. Case-only analysis showed associations of rs2165667 and rs1562444 with arthritis, rs10830962 with malar rash, rs3760138 with immunological abnormality, and rs8150 with hematological abnormality. Furthermore, a significant difference between plasma MTN levels with different genotypes of rs1562444 was observed. Haplotype analyses revealed that haplotype of CCTAT, CTAGT, and GGG was significantly associated with the increased risk in SLE susceptibility, but TCTAT and CTG appeared to be a protective haplotype. Conclusions. The present study supported the genetic association of MTN pathway genes with SLE susceptibility and specific clinical manifestations, suggesting the potential role of MTN pathway genes in the pathogenesis and development of SLE.

1. Introduction

Systemic lupus erythematosus (SLE) is a chronic and inflammatory autoimmune disease characterized by antinuclear autoantibody production and a multitude of immune-complex deposition, which is involved in multisystem; such as the skin, kidney, and brain; and caused organs/tissue destruction [1, 2]. Until now, the etiology of SLE is still not fully understood. A large body of literatures have suggested that the interactions between exogenous (infectious triggers, ultraviolet radiation, and dietary habit) and endogenous sources (hereditary susceptibility, endocrine disorders and disturbed status of oxidative metabolism, and autoimmune responsiveness sympathetic nervous system) are responsible for the pathogenesis and development of SLE [3–9].

Melatonin (MTN), as one of the major neuroendocrine hormones, is mainly produced and secreted by the pineal gland; it allows to regulate the circadian day-night rhythm and seasonal biorhythms and is also a key player in the neuroendocrine-immune pathway [10, 11]. Several studies have suggested the important role of MTN on the regulation of the immune system; it can skew the immune response by repressing the production of several proinflammatory cytokines (such as tumor necrosis factor- (TNF-) α, interleukin- (IL-) 1β, and IL-6), as well as blocking the DNA-binding activity of NF-κB both in vitro and in vivo, exerting an anti-inflammatory effect [11–16].

Previous studies have investigated the underlying role of MTN in SLE. Lechner et al. observed an increased MTN level in MRL/MP-fas lupus-prone mice (represents an animal model for human autoimmune diseases, which spontaneously develops lupus-like glomerulonephritis, systemic vasculitis, arthritis, and sialadenitis) during the light phase; under the regulation of MTN, the levels of autoantibodies were reduced, and the histological changes were improved in female lupus MRL/MP-fas lupus-prone mice [17]. MTN is able to inhibit IgM, anti-dsDNA, and anti-histone antibodies, thus, decreases the levels of IL-6 and IL-13 and increases the IL-12 levels [18]. In patients with SLE, a lower daily MTN level was observed as compared to healthy controls, and this decreased daily MTN level inversely correlated with the systemic lupus erythematosus disease activity index (SLEDAI) [19]. In addition, study has also revealed the seasonal pattern of MTN levels in SLE, with an elevated daily plasma MTN levels in December than in June [20]. These findings provided the possible evidence that MTN might play a potential role in the pathogenesis of SLE.

The synthesis and function of MTN mainly depend on three MTN pathway genes, MTN receptors 1a/1b (MTNR1a/MTNR1b) are largely responsible for mediating the downstream effects of MTN, and arylalkylamine N-acetyltransferase (AANAT) is the major enzyme in MTN synthesis [21, 22]. The genetic association of MTN genes with some diseases has been demonstrated in several studies, including SLE, multiple sclerosis (MS), breast cancer, and major depression [23–26]. Nevertheless, associations between genetic variation in MTN pathway genes and SLE susceptibility have not been determined.

In the present study, we conducted a case–control study to comprehensively evaluate the role of common genetic variation in the MTNR1a, MTNR1b, and AANAT genes to SLE susceptibility in a Chinese population.

2. Materials and Methods

2.1. Study Subjects and Methods

This case–control genotyping study recruited a total of 988 subjects (495 SLE patients and 493 healthy controls). The sample size and power calculation of the study were computed by power and sample size program, where both the minor allele frequency (MAF) and statistical significant level were set as 0.05, odds ratio (OR) was 1.5, case and control ratio was 1 : 1 and when the statistical power was 0.8, the computed sample size for cases was 451. In addition, we also evaluated the statistical power for the 496 included cases, and the results showed that the computed statistical power was 0.833.

Patients with SLE were recruited from the Department of Rheumatology and Immunology at Anhui Provincial Hospital, The First Affiliated Hospital of Anhui Medical University. The diagnosis of SLE was established by the presence of four or more 1997 revised American College of Rheumatology (ACR) classification criteria [27]. Patients with viral infections and any history of cancer, pregnancy, and recurrent spontaneous abortions were excluded based on reviews of their appropriate history. The normal controls consisted of age, gender, and ethnicity-matched healthy individuals who belonged to the same geographical area as that of cases; normal controls were excluded if they had a family history of SLE or any other autoimmune disease and history of any chronic or lifestyle diseases like depression, hypothyroidism, hypertension, diabetes mellitus, and tuberculosis (TB). Demographics, clinical features, and related laboratory results were obtained from hospital medical records and then reviewed by experienced physicians.

The Ethical Committee of Anhui Medical University (Hefei, Anhui, China) approved this study. All the study subjects provided informed consent to participate in this study.

All studies on humans described in the present manuscript were carried out with the approval of the responsible ethics committee and in accordance with the national law and the Declaration of Helsinki 1975 (in its current, revised form).

2.2. MTN Pathway Gene’s SNP Selection and Genotyping

Ensembl Gene Browser 37 (GRCh37) (http://grch37.ensembl.org/index.html) (Ensembl Archive Release 90) was implemented to acquire the genetic and location information of MTNR1a, MTNR1b, and AANAT genes [28], and linkage pedigree file (PED) and marker information file were downloaded. Then, the downloaded files were used to select the tag SNPs in Haploview 4.2 software (Broad Institute, Cambridge, MA, USA), with a MAF above 0.05 in Chinese Han population (CHB) of Beijing, linkage disequilibrium (LD) with an threshold of 0.8. A total of 46 tag SNPs (23 MTNR1a tag SNPs, 13 MTNR1b tag SNPs, and 10 AANAT tag SNPs) were captured for further evaluation. The function prediction for 46 tag SNPs was assessed by the online bioinformatics tools (https://snpinfo.niehs.nih.gov/snpinfo/snpfunc.html) [29]; the basic information of these tag SNPs is shown in Supplementary Table 1. In addition, the published literatures about the MTN pathway gene’s polymorphisms were also carefully reviewed. Based on one or more of the following criteria: assumed functionality (located in the regulatory regions, for example, 3-untranslated regions (UTR), 5-UTR, or amino acid change), and , as well as previous studies reported SNPs), the five tag SNPs (rs10030173, rs2119882, rs2165667, rs4861722, and rs6847693) in MTNR1a gene, three tag SNPs (rs1562444, rs10830962, and rs3781637) in MTNR1b gene, and three tag SNPs (rs8150, rs3760138, and rs12942767) in AANAT gene were finally chosen for further genotyping.

Genomic DNA was extracted from peripheral venous blood of patients and healthy controls using a QIAGEN kit (QIAGEN, Hilden, Germany) based on the manufacturer’s instructions, and quantification and concentration of DNA was determined using NanoDrop 2000 Spectrophotometer (Thermo Fisher Scientific, USA). Qualified sample requirements were shown as follows: concentration greater than 50 ng/μl, total amount greater than 600 ng, and no obvious degradation.

Genotyping was performed in using improved multiple ligase detection reaction (iMLDR), with technical support from the Center for Genetic & Genomic Analysis, Genesky Biotechnologies Inc., Shanghai. A multiplex PCR-ligase detection reaction method was used in the iMLDR. For each SNP, the alleles were distinguished by different fluorescent labels of allele-specific oligonucleotide probe pairs. Different SNPs were further distinguished by different extended lengths at the 3 end. Two negative controls were set: one with double-distilled water as template and the other with DNA sample without primers while keeping all other conditions the same in one plate. Duplicate tests were designed, and the results were consistent. A random sample accounting for ~3% of the total DNA samples was directly sequenced using BigDye Terminator version 3.1 and an ABI3730XL automated sequencer (Applied Biosystems) to confirm the results of iMLDR.

2.3. Plasma MTN Determination

Blood samples were collected from 5 ml of whole blood of all study subjects and then stored at -80°C until assayed. Plasma MTN concentration was determined by enzyme-linked immunosorbent assay (ELISA) kits (the lower detection limit was 0.1 pg/ml) which were purchased from Anhui Xinle Biotechnology Co. Ltd.; the results of MTN were expressed as picograms per milliliter. The interassay and intraassay variation coefficients of the ELISA kit of our study were 6.9% and 7.6%, respectively.

2.4. Statistical Analysis

The allelic and genotypic association analyses between SLE patients and healthy controls were performed in using the Chi-square or Fisher’s exact test. Logistic regression analyses were utilized to calculate odds ratios (ORs), and 95% confidence intervals (CIs) for the association between genotype and SLE susceptibility, additive, dominant, recessive, and allelic models were also considered. A nonparametric test was used to compare the difference in plasma MTN levels among patients with different genotypes. Statistical analysis was implemented with the use of the SPSS (IBM Corp. Released 2015. IBM SPSS Statistics for Windows, Version 23.0. Armonk, NY: IBM Corp.).

Hardy–Weinberg equilibrium (HWE) among controls was assessed by comparing the observed-to-expected genotype frequencies using the Chi-square test. Online software SHEsis was used for haplotype analyses of each MTN pathway genes; all the haplotypes with a were ignored in the analysis [30]. All results with a two-tailed were considered to be statistically significant. The Bonferroni correction was used for multiple testing.

3. Results

3.1. Characteristics of the Study Population

The present study recruited 988 subjects, with 495 SLE patients and 493 healthy controls. In SLE patients, there were 57 males and 438 females with a median age of 37.00 (28.00, 46.00) years and the median disease duration was 4.10 (range from 1.05 to 9.06) years, while there were 55 males and 434 females with the median age of 38.00 (30.00, 47.00) in healthy controls. No significant differences in gender and age distribution were observed between SLE patients and normal controls (Table 1). The major clinical manifestations of SLE were immunological abnormality (73.1%), hematological abnormality (68.3%), arthritis (49.5%), malar rash (45.3%), photosensitivity (39.0%), and renal abnormality (37.2%). In control groups, the presence of observed genotype frequencies of all included tag SNPs was distributed in compliance with the HWE (all ).

Table 1: Demographic characteristics and clinical features of patients with SLE and control subjects.

Parameters
Patients with SLE ()
Healthy controls ()

Demographic characteristics

    Age (year)
37.00 (28.00, 46.00)
38.00 (30.00, 47.00)

    Female, (%)
438 (88.5)
434 (88.0)

    Male, (%)
58 (11.5)
59 (12.0)

    Disease duration (year)
4.10 (1.05, 9.06)
NA

    BMI (kg/m²)
21.48 (19.72, 23.43)
NA

    SLEDAI
11.40 ± 9.07
NA

Disease manifestations

NA

    Malar rash, (%)
224 (45.3)
NA

    Discoid rash, (%)
94 (19.0)
NA

    Photosensitivity, (%)
193 (39.0)
NA

    Oral ulcers, (%)
119 (24.0)
NA

    Arthritis, (%)
245 (49.5)
NA

    Serositis, (%)
45 (9.1)
NA

    Renal abnormality, (%)
184 (37.2)
NA

    Neurological abnormality, (%)
21 (4.2)
NA

    Hematological abnormality, (%)
338 (68.3)
NA

    Immunological abnormality, (%)
362 (73.1)
NA

BMI: Body mass index; : number; SLE: systemic lupus erythematosus; SLEDAI: systemic lupus erythematosus disease activity index. Data with normal distribution were described as mean and standard deviation (SD). Data with skewed distribution were described as median and interquartile range (IQR).

3.2. Association of MTNR1a/b and AANAT Gene’s Polymorphisms with Susceptibility to SLE

There were no significant differences in allele and genotype distribution of 8 tag SNPs in MTNR1a/b genes between SLE and healthy controls (all ) (Table 2). However, when analyzing the allele and genotype frequency of 3 tag SNPs in AANAT genes, the results showed that two SNPs of rs8150 and rs3760138 were associated with the risk of SLE susceptibility, where CC carriers of rs8150 had a lower risk as compared to GG (, 95% CI: 0.361, 0.799) (), whereas GG carrier of rs3760138 had an increased risk as compared to TT (, 95% CI: 1.154, 2.880) (), but we did not observe other positive findings regarding the SNPs of rs12942767 (Table 3).

Table 2: Genotype frequencies of MTNR1a/b SNPs in SLE patients and healthy controls.

SNPs
Analyzed model

SLE
Control
OR (95% CI)
value

rs10030173 (MTNR1a)
Genotypes
CC
107
93
1.308 (0.916, 1.867)
0.139

CT
242
234
1.176 (0.883, 1.565)
0.267

TT
146
166
1.000

Additive model
CC
107
93
1.308 (0.916, 1.867)
0.139

TT
146
166
1.000

rs2119882 (MTNR1a)
Genotypes
CC
71
78
0.792 (0.542, 1.156)
0.226

CT
217
235
0.803 (0.612, 1.054)
0.114

TT
207
180
1.000

Additive model
CC
71
78
0.792 (0.542, 1.156)
0.226

TT
207
180
1.000

rs2165667 (MTNR1a)
Genotypes
AA
230
200
1.267 (0.849, 1.890)
0.246

AT
206
228
0.995 (0.667, 1.484)
0.982

TT
59
65
1.000

Additive model
AA
230
200
1.267 (0.849, 1.890)
0.246

TT
59
65
1.000

rs4861722 (MTNR1a)
Genotypes
AA
25
28
0.978 (0.548, 1.743)
0.939

GA
165
189
1.243 (0.707, 2.183)
0.450

GG
304
274
1.000

Additive model
AA
25
28
0.978 (0.548, 1.743)
0.939

GG
304
274
1.000

rs6847693 (MTNR1a)
Genotypes
CC
182
167
0.868 (0.658, 1.145)
0.318

CT
229
242
0.918 (0.635, 1.326)
0.647

TT
84
84
1.000

Additive model
CC
182
167
0.868 (0.658, 1.145)
0.318

TT
84
84
1.000

rs1562444 (MTNR1b)
Genotypes
AA
254
262
1.039 (0.659, 1.636)
0.870

GA
199
186
1.146 (0.720, 1.826)
0.565

GG
42
45
1.000

Additive model
AA
254
262
1.039 (0.659, 1.636)
0.870

GG
42
45
1.000

rs10830962 (MTNR1b)
Genotypes
CC
176
145
1.349 (0.951, 1.913)
0.093

CG
220
236
1.036 (0.746, 1.438)
0.834

GG
99
110
1.000

Additive model
CC
176
145
1.349 (0.951, 1.913)
0.093

GG
99
110
1.000

rs3781637 (MTNR1b)
Genotypes
CC
11
11
1.009 (0.420, 2.424)
0.984

CT
111
110
1.003 (0.429, 2.341)
0.995

TT
373
372
1.000

Additive model
CC
11
11
1.009 (0.420, 2.424)
0.984

TT
373
372
1.000

SLE: systemic lupus erythematosus; SNPs: single-nucleotide polymorphisms; OR: odds ratio; MTNR1a: melatonin receptor 1a; MTNR1b: melatonin receptor 1b.The values are not corrected for multiple testings, Bonferroni corrected .

Table 3: Genotype frequencies of AANAT SNPs in SLE patients and healthy controls.

SNPs
Analyzed model

SLE
Control
OR (95% CI)
value

rs8150
Genotypes
CC
176
202
0.537 (0.361, 0.799)
0.002

GC
232
236
1.128 (0.860, 1.480)
0.384

GG
86
53
1.000

Additive model
CC
176
202
0.537 (0.361, 0.799)
0.002

GG
86
53
1.000

rs3760138
Genotypes
GG
271
223
1.823 (1.154, 2.880)
0.010

GT
188
216
1.306 (0.820, 2.078)
0.261

TT
36
54
1.000

Additive model
GG
271
223
1.823 (1.154, 2.880)
0.010

TT
36
54
1.000

rs12942767
Genotypes
GG
450
447
—
1.000

GA
44
46
—
1.000

AA
1
0
1.000

Additive model
GG
450
447
—
1.000

AA
1
0
1.000

SLE: systemic lupus erythematosus; SNPs: single-nucleotide polymorphisms; OR: odds ratio; AANAT: arylalkylamine N-acetyltransferase. The values are not corrected for multiple testings, Bonferroni corrected .

3.3. Association of MTNR1a/b and AANAT Gene’s Polymorphisms with Clinical Features in Patients with SLE

Case-only analysis was conducted to further explore the genetic association of MTNR1a/b and AANAT gene’s polymorphisms with specific clinical features of SLE. In MTNR1a/b genes, a significantly increased AA genotype frequency of rs2165667 (MTNR1a) and A allele frequency of rs1562444 (MTNR1b) were found in patients with arthritis than those without (both ). In addition, the frequency of CC/CG genotype in rs10830962 (MTNR1b) was significantly lower in patients with malar rash than in those without (). In terms of the genetic association of AANAT gene with clinical features of SLE, there was a higher GG/GT/TT genotype distribution of rs3760138 in patients with positive immunological abnormality than those with negative (); the C/G allele and CC/CG/GG genotype frequency of rs8150 appeared to have a significantly increased risk in patients with positive hematological abnormality compared with those with negative (, , respectively). Nevertheless, no other positive findings were revealed regarding the MTNR1a/b and AANAT gene’s polymorphisms with SLE clinical features (Table 4).

Table 4: The positive findings on association of clinical characteristics with genotype and allele frequencies in MTNR1a/b and AANAT genes.

Gene (SNPs)
Allele
(M/m)
Clinical features
Group
Genotypes ()
value
Alleles ()
value

MM
Mm
mm

rs2165667 (MTNR1a)
A/T
Arthritis
Positive
125
87
33
0.024
337
153
0.318

Negative
105
119
26

329
171

rs10830962 (MTNR1b)
C/G
Malar rash
Positive
77
90
57
0.018
244
204
0.055

Negative
99
130
42

328
214

rs1562444 (MTNR1b)
A/G
Arthritis
Positive
136
94
15
0.072
366
124
0.024

Negative
118
105
27

341
159

rs3760138 (AANAT)
G/T
Immunological abnormality
Positive
198
131
33
0.024
527
197
0.264

Negative
73
57
3

203
63

rs8150 (AANAT)
C/G
Hematological abnormality
Positive
131
155
51
0.039
417
257
0.010

Negative
45
77
35

167
147

SNPs: single-nucleotide polymorphisms; MTNR1a: melatonin receptor 1a; MTNR1b: melatonin receptor 1b; AANAT: arylalkylamine N-acetyltransferase.

3.4. Association of Plasma MTN Concentrations with Genotypes in Patients with SLE

The results indicated that, in patients with SLE, there was a significant difference of MTN level among genotype of AA, AG, and GG in rs1562444 (MTNR1b) (), which GG genotype showed an elevated MTN concentration than in AA and AG genotype (20.57 pg/ml vs 14.08 pg/ml vs 10.36 pg/ml). However, no significant differences of plasma MTN concentrations were observed in other SNPs between patients (Table 5).

Table 5: Association of plasma MTN levels with genotype in MTNR1a/b and AANAT.

SNPs
Genotypes
Number
Plasma MTN levels (pg/ml)
value

(₂₅, ₇₅)

rs10030173 (MTNR1a)
CC
19
14.22 (9.23, 20.57)
0.946

CT
45
13.45 (9.91, 19.06)

TT
19
12.39 (10.27, 18.69)

rs2119882 (MTNR1a)
CC
9
12.19 (9.54, 25.08)
0.312

CT
44
13.79 (10.34, 20.00)

TT
30
10.82 (9.16, 16.33)

rs2165667 (MTNR1a)
AA
32
10.69 (8.41, 16.12)
0.079

AT
44
13.24 (10.39, 19.91)

TT
7
19.38 (10.49, 28.21)

rs4861722 (MTNR1a)
AA
6
15.47 (8.24, 34.58)
0.640

GA
30
11.74 (10.19, 16.35)

GG
47
13.71 (9.91, 21.46)

rs6847693 (MTNR1a)
CC
25
10.62 (9.08, 16.41)
0.333

CT
47
13.45 (10.36, 20.20)

TT
11
18.37 (8.59, 21.95)

CT
1
11.07

rs1562444 (MTNR1b)
AA
41
14.08 (10.63, 20.08)
0.001

AG
35
10.36 (8.94, 13.86)

GG
7
20.57 (18.37, 28.21)

rs10830962 (MTNR1b)
CC
25
13.71 (10.32, 22.10)
0308

CG
49
11.58 (9.80, 18.30)

GG
9
15.57 (10.05, 33.40)

rs3781637 (MTNR1b)
CC
3
18.37 (13.71, 28.21)
0.392

CT
16
13.12 (8.57, 20.27)

TT
64
12.55 (10.25, 17.83)

rs8150 (AANAT)
CC
36
13.68 (10.28, 23.20)
0.518

GC
34
11.88 (9.62, 18.45)

GG
13
13.04 (10.44, 14.82)

rs3760138 (AANAT)
GG
45
12.19 (9.78, 17.09)
0.802

GT
31
13.45 (10.05, 21.95)

TT
7
14.27 (12.01, 21.46)

rs12942767 (AANAT)
GA
5
13.03 (9.88, 16.83)
0.878

GG
78
13.07 (9.92, 19.59)

SNPs: single-nucleotide polymorphisms; : median; MTN: melatonin; MTNR1a: melatonin receptor 1a; MTNR1b: melatonin receptor 1b; AANAT: arylalkylamine N-acetyltransferase.

3.5. Haplotype Analyses

The haplotype of tag SNPs in MTNR1a/b and AANAT genes was constructed by using SHEsis software. Haplotype analyses implied that haplotype CCTAT (MTNR1a), CTAGT (MTNR1a), and GGG (AANAT) were significantly associated with the increased risk in SLE susceptibility, but TCTAT (MTNR1a) and CTG (AANAT) appeared to be a protective haplotype (all ). However, no positive findings of other haplotypes were observed (Tables 6–8).

Table 6: Haplotype analysis of five SNPs in MTNR1a gene in SLE patients and healthy controls.

Haplotype
SLE ( (%))
Controls ( (%))

OR (95% CI)

rs10030173, rs2119882, rs2165667, rs4861722, rs6847693

CCTAT
44.05 (0.045)
27.21 (0.028)
4.135
0.042031
1.650 (1.014, 2.686)

CTAGC
341.42 (0.346)
331.12 (0.337)
0.241
0.623479
1.049 (0.866, 1.270)

CTAGT
38.30 (0.039)
19.92 (0.020)
5.983
0.014466
1.962 (1.133, 3.400)

TCTAT
135.01 (0.137)
180.70 (0.184)
8.071
0.004512
0.702 (0.549, 0.897)

TCTGT
117.73 (0.119)
118.09 (0.120)
0.001
0.974785
0.996 (0.757, 1.309)

TTAGC
217.53 (0.220)
216.35 (0.220)
0.003
0.954914
1.006 (0.811, 1.249)

SLE: systemic lupus erythematosus; SNPs: single-nucleotide polymorphisms; MTNR1a: melatonin receptor 1a; OR: odds ratio. Total , . All the haplotypes with a were ignored in the analysis.

Table 7: Haplotype analysis of three SNPs in MTNR1b gene in SLE patients and healthy controls.

Haplotype
SLE ( (%))
Control ( (%))

OR (95% CI)

rs1562444, rs10830962, rs3781637

CAT
296.17 (0.299)
260.80 (0.265)
2.573
0.108701
1.175 (0.965, 1.430)

CGC
130.72 (0.132)
127.39 (0.129)
0.014
0.906006
1.016 (0.782, 1.320)

CGT
144.94 (0.146)
137.77 (0.140)
0.123
0.725622
1.046 (0.813, 1.346)

GAT
408.55 (0.413)
442.60 (0.450)
3.208
0.073303
0.849 (0.710, 1.016)

SLE: systemic lupus erythematosus; SNPs: single-nucleotide polymorphisms; MTNR1b: melatonin receptor 1b; OR: odds ratio. Total , . All the haplotypes with a were ignored in the analysis.

Table 8: Haplotype analysis of three SNPs in AANAT gene in SLE patients and healthy controls.

Haplotype
SLE ( (%))
Control ( (%))

OR (95% CI)

rs8150, rs3760138, rs12942767

CGG
339.49 (0.344)
324.11 (0.329)
0.559
0.454842
1.074 (0.891, 1.295)

CTG
244.50 (0.247)
315.27 (0.320)
12.475
0.000415
0.701 (0.576, 0.854)

GGA
43.12 (0.044)
44.37 (0.045)
0.018
0.894335
0.971 (0.633, 1.492)

GGG
347.39 (0.352)
291.50 (0.296)
7.336
0.006775
1.299 (1.075, 1.571)

SLE: systemic lupus erythematosus; SNPs: single-nucleotide polymorphisms; AANAT: arylalkylamine N-acetyltransferase; OR: odds ratio. Total , . All the haplotypes with a were ignored in the analysis.

4. Discussion

The neuroendocrine-immune system is regarded as a fundamental network supporting the health state that could play an important role in the development of autoimmune disorders [31, 32]. MTN, as one of the pineal gland-driven hormones, secreted in a circadian rhythm and regulated by photoperiod with the highest peak at midnight and lowest level after sunrise, acting mainly as a regulator for sleeping rhythm [33–36]. In recent years, the bidirectional associations between the pineal gland and the immune system have been suggested to depend on the immune-modulating effect of MTN and the pineal regulation by different lymphokines [37]. MTN might have a direct effect on immune-competent cells, thus, fulfilling critical roles on the development and progression of autoimmune diseases [38–40].

Our previous study has evaluated the level of MTN in patients with SLE as compared to healthy controls. Although there was no significant difference in MTN concentration between those two groups, we observed a slightly lower level of MTN in SLE patients than in healthy controls; in addition, an inverse correlation of MTN concentration with IgM was also revealed [41]. Similarly, Robeva et al. revealed that there was a decreased daily MTN level in women with SLE and found an inverse relationship between daily MTN concentrations and disease activity [19]. In the subarctic region, the study has demonstrated the presence of seasonal variations in daily MTN level, where the increased level of MTN was discovered in December than that in June [20]. These evidences suggest that the change of MTN concentration may be involved in the pathogenesis of SLE.

MTNR1a/b and AANAT play an important role during the pathway of MTN from biosynthesis to its functioning; the former MTNR1a/b genes are largely responsible for mediating the downstream effects of MTN, while the latter AANAT gene is the major enzyme in MTN synthesis. Several previous studies have investigated the underlying role of MTN pathway genes in a number of human diseases. In SLE, primary study, performed by Tanev et al., has demonstrated no significant differences in allelic and genotype distribution of MTNR1b gene (rs1562444, rs10830962, and rs10830963) polymorphisms between 109 patients with SLE and 101 controls, yet, in SLE patients, C/C genotype of rs10830963 in MTNR1b gene was related to increased prevalence of leucopenia compared to C/G and G/G genotype; the rs10830963 G/G carriers had a lower number of lupus criteria than in those with C/C genotype [25]. The MTNR1b rs10830962 and rs10830963 polymorphisms have been predominantly investigated in the context of metabolic disorders, of which rs10830962 and rs10830963 G alleles were reported to associate with reduced insulin secretion, increased fasting plasma glucose concentrations, and increased risk for diabetes in different populations [42–44]. In autoimmune diseases of multiple sclerosis (MS), there were no significant allelic associations of SNPs rs4753426 and rs10830963 in MTNR1b gene with susceptibility to MS, but the rs10830963-rs4753426 G-T haplotype associated with the risk of MS in the progressive MS group [23]. Deming et al. analyzed the MTN pathway gene’s polymorphisms in human breast cancer patients, and they supported that AA genotype of MTNR1b rs10765576 was associated with a decreased risk of breast cancer, the GG genotype in premenopausal women correlated with an increased risk for breast cancer, and however, in postmenopausal women, the GG genotype were related with a decreased risk of breast cancer; they did not observe any significant breast cancer associations for variants in the AANAT gene [24]. In patients with major depression, the two SNPs of AANAT (rs3760138 and rs4238989) were reported to be associated with an increased contribution to major depression [26].

In the present study, the tag SNPs of rs8150 and rs3760138 in AANAT gene were associated with genetic susceptibility to SLE, but no genetic association regarding the other nine tag SNPs with SLE susceptibility was found. Case-only analysis indicated that AA genotype frequency in rs2165667 (MTNR1a) AA genotype and rs1562444 (MTNR1b) A/G allele frequency were at increased risk for arthritis and rs10830962 (MTNR1b) CC/CG genotype was at decreased risk for malar rash. In AANAT gene, rs3760138 GG/GT/TT genotype associated with positive immunological abnormality than those with negative, and rs8150 CC/CG/GG genotype and its C/G allele appeared to have an increased risk for hematological abnormality. Moreover, we also found that there was a significant difference of MTN concentration samong the genotype of AA, AG, and GG in rs1562444 (MTNR1b), where the GG genotype showed an elevated MTN concentration than in AA and AG genotype. We might hypothesize that the rs1562444 variant polymorphisms lead to the aberrant expression of MTN in patients with SLE. Later, the haplotype of MTNR1a/b and AANAT was identified; MTNR1a gene haplotype of CCTAT and CTAGT, and AANAT gene haplotype of GGG showed an increased risk in SLE susceptibility, but haplotype of TCTAT (MTNR1a) and CTG (AANAT) appeared to have a protective role.

The present study investigated the genetic association of SNPs in MTN pathway genes (MTNR1a/b and AANAT) with SLE susceptibility. However, there are some limitations in our study. First, the current study might be due to inherent selection biases such as a relative small sample size, the limited number of variables accounted, and the lack of information regarding body mass index (BMI) in healthy controls. Second, the potential confounding factors, such as type of treatments and concomitant infections, may have an effect on the level of MTN. Furthermore, although our study represented significant genetic variations regarding AANAT gene in Chinese SLE patients, the detailed mechanism about potential effect of the AANAT gene variation on SLE is scarce.

In conclusion, our study demonstrated that, in the Chinese population, the genetic polymorphism of MTN pathway genes associated with the susceptibility to SLE, as well as with specific clinical manifestations, suggesting that the MTN pathway genes might be involved in the pathogenesis and development of SLE. However, further large sample size studies in other population are needed to further reveal the significance of MTN pathway gene’s polymorphisms in SLE. In addition, related mechanism researches are necessary to better understand the function of the MTN pathway gene SNP in different immune cell types and to evaluate its correlation with clinical features.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Ethical Approval

This study was approved by the Ethical Committee of Anhui Medical University (Hefei, Anhui, China). All the study subjects provided informed consent to participate in this study. All studies on humans described in the present manuscript were carried out with the approval of the responsible ethics committee and in accordance with national law and the Declaration of Helsinki 1975 (in its current, revised form).

Conflicts of Interest

The authors confirm that there are no conflicts of interest.

Acknowledgments

The authors thank the study participants and the staff involved in the collection of blood samples. This work was supported by grants from the National Natural Science Foundation of China (81573222).

Supplementary Materials

Supplementary Table 1: characteristics of the 46 Tag SNPs. (Supplementary Materials)

References

R. Cervera, M. A. Khamashta, and G. R. Hughes, “The Euro-lupus project: epidemiology of systemic lupus erythematosus in Europe,” Lupus, vol. 18, no. 10, pp. 869–874, 2009.

G. C. Tsokos, “Systemic lupus erythematosus,” The New England Journal of Medicine, vol. 365, no. 22, pp. 2110–2121, 2011.

J. C. Crispin, C. M. Hedrich, and G. C. Tsokos, “Gene-function studies in systemic lupus erythematosus,” Nature Reviews Rheumatology, vol. 9, no. 8, pp. 476–484, 2013.

L. M. Olsson, A. C. Johansson, B. Gullstrand et al., “A single nucleotide polymorphism in the NCF1 gene leading to reduced oxidative burst is associated with systemic lupus erythematosus,” Annals of the Rheumatic Diseases, vol. 76, no. 9, pp. 1607–1613, 2017.

L. Wen, C. Zhu, Z. Zhu et al., “Exome-wide association study identifies four novel loci for systemic lupus erythematosus in Han Chinese population,” Annals of the Rheumatic Diseases, vol. 77, no. 3, p. 417, 2018.

J. Li, G. C. Wu, T. P. Zhang et al., “Association of long noncoding RNAs expression levels and their gene polymorphisms with systemic lupus erythematosus,” Scientific Reports, vol. 7, no. 1, article 15119, 2017.

K. Tizaoui, S. H. Kim, G. H. Jeong et al., “Association of PTPN22 1858C/T polymorphism with autoimmune diseases: a systematic review and Bayesian approach,” Journal of Clinical Medicine, vol. 8, no. 3, p. 347, 2019.

K. Lee, A. Kronbichler, D. Pereira Vasconcelos et al., “Genetic variants in antineutrophil cytoplasmic antibody-associated vasculitis: a Bayesian approach and systematic review,” Journal of Clinical Medicine, vol. 8, no. 2, p. 266, 2019.

W. D. Xu, L. Fu, X. Y. Liu et al., “Association between TL1A gene polymorphisms and systemic lupus erythematosus in a Chinese Han population,” Journal of Cellular Physiology, vol. 234, no. 12, pp. 22543–22553, 2019.

D. X. Tan, L. C. Manchester, R. Hardeland et al., “Melatonin: a hormone, a tissue factor, an autocoid, a paracoid, and an antioxidant vitamin,” Journal of Pineal Research, vol. 34, no. 1, pp. 75–78, 2003.

D. Pozo, M. Delgado, J. M. Fernandez-Santos et al., “Expression of the Mel1a-melatonin receptor mRNA in T and B subsets of lymphocytes from rat thymus and spleen,” FASEB Journal, vol. 11, no. 6, pp. 466–473, 1997.

S. Garcia-Maurino, D. Pozo, J. R. Calvo, and J. M. Guerrero, “Correlation between nuclear melatonin receptor expression and enhanced cytokine production in human lymphocytic and monocytic cell lines,” Journal of Pineal Research, vol. 29, no. 3, pp. 129–137, 2000.

S. Garcia-Maurino, M. G. Gonzalez-Haba, J. R. Calvo, R. Goberna, and J. M. Guerrero, “Involvement of nuclear binding sites for melatonin in the regulation of IL-2 and IL-6 production by human blood mononuclear cells,” Journal of Neuroimmunology, vol. 92, no. 1-2, pp. 76–84, 1998.

M. Aparicio-Soto, C. Alarcon-de-la-Lastra, A. Cardeno, S. Sanchez-Fidalgo, and M. Sanchez-Hidalgo, “Melatonin modulates microsomal PGE synthase 1 and NF-E2-related factor-2-regulated antioxidant enzyme expression in LPS-induced murine peritoneal macrophages,” British Journal of Pharmacology, vol. 171, no. 1, pp. 134–144, 2014.

D. Acuña-Castroviejo, G. Escames, C. Venegas et al., “Extrapineal melatonin: sources, regulation, and potential functions,” Cellular and Molecular Life Sciences, vol. 71, no. 16, pp. 2997–3025, 2014.

V. Raghavendra, V. Singh, A. V. Shaji, H. Vohra, S. K. Kulkarni, and J. N. Agrewala, “Melatonin provides signal 3 to unprimed CD4⁺ T cells but failed to stimulate LPS primed B cells,” Clinical and Experimental Immunology, vol. 124, no. 3, pp. 414–422, 2001.

O. Lechner, H. Dietrich, A. Oliveira dos Santos et al., “Altered Circadian Rhythms of the Stress Hormone and Melatonin Response in Lupus-prone MRL/MP-fasIpr Mice,” Journal of Autoimmunity, vol. 14, no. 4, pp. 325–333, 2000.

L. L. Zhou, W. Wei, J. F. Si, and D. P. Yuan, “Regulatory effect of melatonin on cytokine disturbances in the pristane- induced lupus mice,” Mediators of Inflammation, vol. 2010, Article ID 951210, 7 pages, 2010.

R. Robeva, D. Tanev, G. Kirilov et al., “Decreased daily melatonin levels in women with systemic lupus erythematosus - a short report,” Balkan Medical Journal, vol. 30, no. 3, pp. 273–276, 2013.

H. J. Haga, J. G. Brun, O. P. Rekvig, and L. Wetterberg, “Seasonal variations in activity of systemic lupus erythematosus in a subarctic region,” Lupus, vol. 8, no. 4, pp. 269–273, 1999.

M. Gonzalez-Arto, T. R. Hamilton, M. Gallego et al., “Evidence of melatonin synthesis in the ram reproductive tract,” Andrology, vol. 4, no. 1, pp. 163–171, 2016.

E. Velarde, J. M. Cerda-Reverter, A. L. Alonso-Gomez, E. Sanchez, E. Isorna, and M. J. Delgado, “Melatonin-synthesizing enzymes in pineal, retina, liver, and gut of the goldfish (Carassius): mRNA expression pattern and regulation of daily rhythms by lighting conditions,” Chronobiology International, vol. 27, no. 6, pp. 1178–1201, 2010.

R. Natarajan, E. Einarsdottir, A. Riutta et al., “Melatonin pathway genes are associated with progressive subtypes and disability status in multiple sclerosis among Finnish patients,” Journal of Neuroimmunology, vol. 250, no. 1-2, pp. 106–110, 2012.

S. L. Deming, W. Lu, A. Beeghly-Fadiel et al., “Melatonin pathway genes and breast cancer risk among Chinese women,” Breast Cancer Research and Treatment, vol. 132, no. 2, pp. 693–699, 2012.

D. Tanev, R. Robeva, S. Andonova et al., “Melatonin receptor 1b polymorphisms in women with systemic lupus erythematosus,” Acta Reumatologica Portuguesa, vol. 41, no. 1, pp. 62–67, 2016.

V. Soria, È. Martínez-Amorós, G. Escaramís et al., “Resequencing and association analysis of arylalkylamine N‐acetyltransferase (AANAT) gene and its contribution to major depression susceptibility,” Journal of Pineal Research, vol. 49, no. 1, pp. 35–44, 2010.

M. C. Hochberg, “Updating the American College of Rheumatology revised criteria for the classification of systemic lupus erythematosus,” Arthritis and Rheumatism, vol. 40, no. 9, p. 1725, 1997.

B. L. Aken, P. Achuthan, W. Akanni et al., “Ensembl 2017,” Nucleic Acids Research, vol. 45, no. D1, pp. D635–D642, 2017.

Z. Xu and J. A. Taylor, “SNPinfo: integrating GWAS and candidate gene information into functional SNP selection for genetic association studies,” Nucleic Acids Research, vol. 37, pp. W600–W605, 2009.

Z. Li, Z. Zhang, Z. He et al., “A partition-ligation-combination-subdivision EM algorithm for haplotype inference with multiallelic markers: update of the SHEsis (http://analysis.bio-x.cn),” Cell Research, vol. 19, no. 4, pp. 519–523, 2009.

N. Deckx, W. P. Lee, Z. N. Berneman, and N. Cools, “Neuroendocrine immunoregulation in multiple sclerosis,” Clinical and Developmental Immunology, vol. 2013, Article ID 705232, 23 pages, 2013.

S. ThyagaRajan and H. P. Priyanka, “Bidirectional communication between the neuroendocrine system and the immune system: relevance to health and diseases,” Annals of Neurosciences, vol. 19, no. 1, pp. 40–46, 2012.

M. Cutolo, A. Sulli, C. Pizzorni et al., “Circadian rhythms: glucocorticoids and arthritis,” Annals of the New York Academy of Sciences, vol. 1069, pp. 289–299, 2006.

M. Cutolo, B. Villaggio, K. Otsa, O. Aakre, A. Sulli, and B. Seriolo, “Altered circadian rhythms in rheumatoid arthritis patients play a role in the disease's symptoms,” Autoimmunity Reviews, vol. 4, no. 8, pp. 497–502, 2005.

M. Cutolo, K. Otsa, O. Aakre, and A. Sulli, “Nocturnal hormones and clinical rhythms in rheumatoid arthritis,” Annals of the New York Academy of Sciences, vol. 1051, pp. 372–381, 2005.

M. Cutolo and A. T. Masi, “Circadian rhythms and arthritis,” Rheumatic Diseases Clinics of North America, vol. 31, no. 1, pp. 115–129, 2005.

M. Szczepanik, “Melatonin and its influence on immune system,” Journal of Physiology and Pharmacology, vol. 58, Supplement 6, pp. 115–124, 2007.

P. Medrano-Campillo, H. Sarmiento-Soto, N. Álvarez-Sánchez et al., “Evaluation of the immunomodulatory effect of melatonin on the T-cell response in peripheral blood from systemic lupus erythematosus patients,” Journal of Pineal Research, vol. 58, no. 2, pp. 219–226, 2015.

P. Plaimee, M. Khamphio, N. Weerapreeyakul, S. Barusrux, and N. P. Johns, “Immunomodulatory effect of melatonin in SK-LU-1 human lung adenocarcinoma cells co-cultured with peripheral blood mononuclear cells,” Cell Proliferation, vol. 47, no. 5, pp. 406–415, 2014.

E. M. Kuklina, “Melatonin as potential inducer of Th17 cell differentiation,” Medical Hypotheses, vol. 83, no. 3, pp. 404–406, 2014.

P. Wang, H. M. Li, Y. F. Zou, J. H. Tao, and H. F. Pan, “Plasma melatonin levels do not differ in SLE patients,” Zeitschrift für Rheumatologie, vol. 77, pp. 66–70, 2016.

E. Reiling, E. van ’t Riet, M. J. Groenewoud et al., “Combined effects of single-nucleotide polymorphisms in GCK, GCKR, G6PC2, and MTNR1B on fasting plasma glucose and type 2 diabetes risk,” Diabetologia, vol. 52, no. 9, pp. 1866–1870, 2009.

H. Staiger, F. Machicao, S. A. Schäfer et al., “Polymorphisms within the novel type 2 diabetes risk locus MTNR1B determine β-cell function,” PLoS One, vol. 3, no. 12, article e3962, 2008.

V. Lyssenko, C. L. Nagorny, M. R. Erdos et al., “Common variant in MTNR1B associated with increased risk of type 2 diabetes and impaired early insulin secretion,” Nature Genetics, vol. 41, no. 1, pp. 82–88, 2009.