Background

Psoriasis and psoriatic arthritis (PsA) are inflammatory associated autoimmune disorders. MicroRNA (miR)-146a plays a crucial role in regulating inflammation. A single nucleotide polymorphism in the miR-146a gene (rs2910164), aberrantly alters its gene expression and linked with the pathogenesis of several disorders, including psoriasis and PsA. In South Africa, psoriasis and PsA are extremely rare in the indigenous African population and most common in both the Indian and Caucasian population. The aim of this study was to investigate whether the miR-146a rs2910164 contributes towards psoriasis and PsA development in South African Indian and Caucasian patients.

Methods

South African Indian (n = 84) and Caucasian (n = 32) PsA patients (total n = 116) and healthy control subjects (Indian: n = 62 and Caucasian: n = 38; total n = 100) were recruited in the study. DNA was extracted from whole blood taken from all subjects, and genotyped for the miR-146a rs2910164 using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). Data for laboratory parameters were obtained from pathology reports. The consulting rheumatologist collected all other clinical data.

Conclusion

The rs2910164 variant C-allele may play a role in the progression of PsA in the South African Indian population. The main limitation in this study was the small sample size in the case-control cohorts, with a low overall statistical power (post-hoc power analysis = 19%).

MicroRNAs (miRs) are small non-coding RNAs that control gene expression at the post-transcriptional level by negatively regulating the processing, stability, and translation of mRNA. The highly conserved “seed” region of miRs, composed of 2–7 nucleotides and located at the 5′-untranslated region (5’-UTR), binds to the 3’-UTR of their target mRNA to elicit their aforementioned functions [7]. MiRs play an invaluable role in regulating physiological processes in the body, including cell cycle progression, cell differentiation, metabolism and apoptosis [8]. When miRs are aberrantly expressed, due to single nucleotide polymorphisms (SNPs) within miR encoding genes or environmental factors (pollution, teratogens, and smoking), they can also contribute towards the pathogenesis of several inflammatory disorders [9].

MiR-146a is located on human chromosome 5q34 and plays an important role in regulating immune and inflammatory response pathways [10]. MiR-146a induction is stimulated by toll-like receptors (TLRs), interleukin (IL)-1β and tumour necrosis factor (TNF)- α. They primarily target IL receptor associated kinase 1 (IRAK1) and TNF receptor associated factor 6 (TRAF6) to modulate and prevent overstimulation of inflammatory responses in the TLR/NF-κB pathways [11]. The miR-146a G/C SNP (due to a C:U miss-pairing taking place instead of a normal G:U pairing), contributes towards the pathogenesis of several inflammatory diseases, including autoimmune disorders [12], sepsis [13], cardiovascular disease [14] and diabetes [15]. This SNP is situated within the crucial stem region of pre-miRNA-146a and affects the expression of mature miR-146a [16]. The miR-146a rs2910164 is also associated with psoriasis [17] and PsA [18], however limited data are available.

In South Africa, psoriasis and PsA is extremely rare among the indigenous African population while such cases are more common in both the Indian and Caucasian population [19]. Taking this into consideration, the present study investigated whether rs2910164 contributes towards psoriasis and PsA development in South African Indian and Caucasian patients. This was done by comparing the rs2910164 genotype and allele frequency distribution between the Caucasian and Indian PsA patients and healthy controls for any similarities or deviations.

Patient recruitment and sample collection

Blood samples were taken from South African Indian (n = 84) and Caucasian (n = 32) PsA patients (total n = 116) and healthy control subjects (Indian: n = 62 and Caucasian: n = 38; total n = 100) that were enrolled in the study after informed consent following ethical approval from the Pharma-Ethics Research Ethics Committee (ethics reference number: 13095660). The rheumatoid factor-immunoglobulin M (RF-IgM), C-reactive protein (CRP), plasma glucose, blood glycated haemoglobin (HbA1c), total cholesterol, low density lipoprotein (LDL) cholesterol, high density lipoprotein (HDL) cholesterol and Vitamin D25 levels were assessed at Lancet Laboratories (Durban, South Africa), a fully accredited South African National Laboratory. The Health Assessment Questionnaire (HAQ), an assessment tool used for measuring the overall functional health status of patients with PsA, was used to generate HAQ scores. Briefly, the HAQ Disability Index (HAQ-DI) was used to assess the level of functional ability in patients. The HAQ visual analogue (VAS) pain scale was used to assess the absence or presence of PsA related pain and its severity in patients. The HAQ VAS patient global health scale was used to assess the overall quality of life for patients where 0 = good health and 10 = poor health. HAQ score values < 0.5 and > 0.5 indicated patients had minimal functional impairments and moderate to severe functional impairments from the PsA, respectively. The physical measurements of height, weight, and abdominal and waist circumference, and patient history (age, sex, race, disease duration, radiology, HAQ scores and medications) were conducted by the consultant rheumatologist. The inclusion criteria for this study, irrespective of age and gender, were: (a) patients must be over the age of 18 years and have PsA; (b) patients must have a confirmed diagnosis of PsA and must have fulfilled the Classification Criteria for PsA (CASPAR) [20] criteria (patients with all other forms of inflammatory arthritis were excluded from this study); (c) Patients with all other forms of connective tissue disorders were excluded from this study.

DNA extraction

Genomic DNA was extracted from whole blood taken from PsA patients and controls using the FlexiGene® DNA isolation kit (Qiagen). Briefly, 750 μl cell lysis buffer was added to 300 μl whole blood to pellet out the mitochondria and cell nuclei, followed by the removal of contaminants such as proteins in the pellet by adding 150 μl denaturation buffer, which contained a chaotropic salt and protease enzyme, and incubating for 5 min at 65 °C. To this solution, 150 μl 100% isopropanol was added to precipitate out the DNA and was recovered by centrifugation. The DNA was washed in 150 μl 100% ethanol, dried at room temperature, resuspended in 15 μl hydration buffer (10 mM Tris.Cl, pH 8.5), incubated for 1 h at 65 °C and stored at -20 °C until further use. The Nanodrop 2000 spectrophotometer (Thermo Scientific) was used to determine the purity and concentration of the DNA. All DNA samples were standardised to a concentration of 10 ng/μl.

Genotyping

The miR-146a G/C rs2910164 was genotyped using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). The GoTaq® G2 Flexi DNA Polymerase PCR kit (Promega) and the CFX96 Touch™ Real-Time PCR Detection System (Bio-Rad) was used for this analysis. The 147 bp gene amplicon was amplified using 1× Green GoTaq Flexi buffer, 2.5 mM MgCl2, 200 μM of each dNTP, 0.2 Units GoTaq Flexi DNA polymerase, 20 pmol of each forward (F) and reverse (R) primer sequences, and 30 ng genomic DNA template. A no-template control was run with the positive samples to assess the overall specificity of the reaction. The forward and reverse primer sequences used were 5′-CATGGGTTGTGTCAGTGTCAGAGCT-3′, and 5′-TGCCTTCTGTCTCCAGTCTTCCAA-3′, respectively. PCR conditions were: 94 °C for 10 min (initial denaturation), followed by 30 cycles at 94 °C for 30 s (denaturation), 65 °C for 30 s (annealing) and 72 °C for 30 s (extension), and 72 °C for 7 min (final extension). The 147 bp PCR products were electrophoresed on 1.8% agarose gel containing 2 μl GelRed and visualised using the ChemiDoc™ XRS+ Imaging System (Bio-Rad). The Sac I restriction enzyme (New England BioLabs) was used to digest the PCR products at 37 °C for 16 h, electrophoresed on 3% agarose gel containing 2 μl GelRed and visualised as mentioned above. Presence of the homozygous wild-type G-allele (GG genotype) resulted in no cleavage of the 147 bp PCR product. The homozygous variant C-allele (CC genotype) yielded two fragments of 122 and 25 bp. The heterozygous GC genotype yielded three bands of 147, 122 and 25 bp. A DNA ladder was used to accurately determine the different genotypes.

Statistical analysis

The post-hoc power analysis was used to calculate the overall statistical power of the present study [21, 22]. All statistical analysis was performed using the IBM SPSS statistical software (version 24) and GraphPad Prism software (version 5.0) packages. The Student’s unpaired t-test was used to compare the characteristics of PsA patients and the control groups (Table 1). The Chi-squared (χ2) test and Fisher’s exact test were used to analyse the genotype and allele frequencies, respectively (Table 2 and Table 3). The χ2 test was also used to assess whether the genotype frequencies complied with the Hardy-Weinberg equilibrium. The Fisher’s exact test data are represented as the relative risk ratio (RR) and odds ratio (OR) at 95% confidence intervals (CI). Data were expressed as mean ± standard error (Table 1). A p value less than 0.05 was considered statistically significant.

Table 1

Clinical and demographical characteristics of PsA patients and controls

PsA patients

(n = 117)

Controls

(n = 100)

p value

Age (years)

50.34 ± 1.14

46.23 ± 1.56

0.0309*

Sex

Male, n (%)

63 (54)

35 (35)

Female, n (%)

54 (46)

65 (65)

Race

Indian, n (%)

84 (72)

62 (62)

White, n (%)

32 (27)

38 (38)

Mixed Race, n (%)b

1 (1)

0 (0)

BMI (kg/m2)

28.86 ± 0.50

27.85 ± 0.42

0.1307

Smoker

Yes, n (%)

25 (21)

No, n (%)

92 (79)

Disease duration (years)

6.43 ± 0.67

Radiology

Yes, n (%)

84 (72)

No, n (%)

33 (28)

HAQ score

0.62 ± 0.07

RF-IgM

Positive, n (%)

5 (4)

Negative, n (%)

112 (96)

Drugs

MTX, n (%)

111 (95)

SSZ, n (%)

33 (28)

LFM, n (%)

21 (18)

Biologics, n (%)a

9 (8)

CRP (mg/L)

Inclusion

18.95 ± 2.81

0.0011*

@ 6 month

9.68 ± 1.32

Plasma glucose (mmol/l)

6.17 ± 0.21

Total cholesterol (mmol/l)

5.11 ± 0.10

LDL cholesterol (mmol/l)

3.26 ± 0.12

HDL cholesterol (mmol/l)

1.14 ± 0.03

HbA1c (%)

5.85 ± 0.13

Vitamin D25 (ng/ml)

20.31 ± 1.34

Laboratory parameters are represented as mean ± standard error. Comparisons for age, BMI and CRP levels were performed using the unpaired Student’s t-test

The overall statistical power of this study (based on the post-hoc power analysis) was low (19%), however the sample size was adequate to compare the miR-146a G/C rs2910164 genotype and allele frequency distribution between the Caucasian and Indian PsA patients and healthy controls for any similarities or deviations.

In this study, we evaluated the frequency of the miR-146a G/C rs2910164 in South African Indian and Caucasian patients with PsA compared to healthy control subjects. We observed a significantly higher prevalence of the variant C-allele in Indian PsA patients compared to healthy Indian controls (35.71% vs. 22.58% respectively, p = 0.0200, OR = 1.91 95% CI 1.13–3.22). Conversely, the variant C-allele frequency between Caucasian PsA patients and healthy controls were similar. Data suggests that Indian PsA patients with the heterozygous GC and homozygous variant CC genotypes (GC + CC) are more predisposed to developing PsA compared to patients with the homozygous wild-type GG genotype.

Patients with uncontrolled psoriasis and PsA have exacerbated CRP levels [23, 24]. CRP is a biomarker of inflammation and disease severity. Psoriasis, PsA and rheumatoid arthritis are characterised by the secretion of several pro-inflammatory cytokines such as IL-2, IL-6, IL-8, IFN-γ and TNF-α. These cytokines are responsible for producing high levels of CRP and contributing to the pathophysiology of the disease. TNF-α induced secretion of IL-6 activates the production of CRP by stimulating the transcription of the CRP gene via activation of the IRAK1 and signal transducer and activator of transcription 3 (STAT3) inflammatory pathways; thereby, causing the high CRP levels observed in psoriasis, PsA and rheumatoid arthritis [25–27]. The PsA patients recruited in this study had high CRP levels (Table 1).

When miR-146a is highly expressed, it inhibits both IRAK1 and TRAF6 resulting in concomitant reductions in pro-inflammatory cytokines (IL-2, IL-6, IL-8, IFN-γ and TNF-α) expression and CRP levels [14]. The rs2910164 variant C-allele dampens the overall functionality of miR-146a, leading to an upregulation in IRAK1 and TRAF6 expression, resulting in very high cytokine production [13].

In 2014, Zhang et al. reported that Chinese patients homozygous for the wild-type GG genotype and heterozygous for the GC genotypes (GG + GC) compared to patients homozygous for the variant CC genotype had a greater risk of developing psoriasis and PsA [17]. The frequency of the wild-type G-allele was more common in psoriasis patients compared to healthy controls (48.2% versus 42.4%; p = 0.007). The frequency of the GG (21.7%) and GC (53.0%) genotypes were predominant among psoriasis patients versus controls (GG: 16.7% and GC: 51.5%) whereas the CC genotype was more common in the healthy controls versus psoriasis patients (31.8% versus 25.3%). The difference in the distribution of the rs2910164 genotypes between psoriasis patients and controls were statistically significant (p = 0.021). The combined frequency of the GG + GC genotypes were significantly higher among PSA patients compared to the controls (74.7% versus 68.2%; p = 0.018), and was associated with an increase in PsA susceptibility (adjusted OR = 1.38 95% CI 1.06–1.80) [17]. Chatzikyriakidou et al. (2010) found the frequency of the GC genotype to be higher in Greek PsA patients compared to healthy controls (41.4% versus 27.3%). The GG and CC genotypes were higher in the controls compared to PsA patients (GG: 59.1% versus 48.3% and CC: 13.6% versus 10.3%). However, no significant difference was observed in the distribution of the rs2910164 genotypes between PsA patients and controls (p = 0.394) [18].

In our study, the frequency of the GC + CC genotypes and variant C-allele were significantly higher in all PsA patients versus all healthy controls, with no significant changes in genotype distribution between patients and controls. No association between rs2910164 and PsA were noted in the Caucasian population. There was a significant difference in the distribution of the rs2910164 genotypes between Indian PsA patients and controls. Indian PsA patients had a significantly higher frequency of the GC + CC genotypes and variant C-allele versus healthy controls (Table 3).

Psoriasis and PsA patients have a higher prevalence of MetS [29]. In 2012, Langhan et al. reported that patients with mild, moderate and severe psoriasis and PsA had a 22%, 56% and 98% chance of developing MetS, respectively [30]. PsA patients recruited in our study had a mean HAQ score indicative of moderate to severe functional impairments from the disease, and displayed increased fasting plasma glucose and HbA1c levels indicative of early stages of pre-diabetes and impaired fasting glucose (Table 1). Smoking, diet and obesity are some of the environmental risk factors associated with psoriasis and PsA [1]. Cigarette smoke affects the central nervous system and immune system, and causes pro-inflammatory and anti-inflammatory cytokines to be aberrantly expressed [1]. Nicotine binds to the nicotinic acetylcholine receptors (nAChR) found in several cell types, including immune B-cells, T-cells, thymocytes and leukaemic cell lines. This can lead to defective immune and nervous system signalling processes, and can negatively regulate keratinocyte function [1, 31]. Apart from psoriasis and PsA, smoking can also play an insidious role in triggering MetS [32]. In our study, PsA patients had a mean BMI indicating overweight and obesity, and over 21% were active smokers (Table 1).

This study associated the rs2910164 with increased PsA susceptibility in the South African Indian population. A major limitation of the study is the small sample size in the case-control cohorts, with a low overall statistical power (post-hoc power analysis = 19%). The influence of rs2910164 on miR-146a expression and its role in the pathogenesis of PsA necessitates investigation in a bigger cohort.

Acknowledgments

Prof V Chandran (University of Toronto and University Health Network and Mount Sinai Hospitals, Toronto, ON, Canada) for his assistance and guidance.

Funding

None

Availability of data and materials

All data reported in this study is available from Dr AB Maharaj, Dr AA Chuturgoon and Mr P Naidoo. Also, data is available from the Discipline of Medical Biochemistry, School of Laboratory Medicine and Medical Sciences, University of KwaZulu-Natal, Durban, South Africa.

Authors’ contributions

ABM, AAC, PN, PPT and NdeV conceived and designed the study; analysed data and critically revised the draft manuscript. PN, TG, NSA, SD, TFD and PR collected all patient samples, ran the genotyping experiments, analysed the data and prepared the draft manuscript. All authors have read and approved the final version of the manuscript.

Ethics approval and consent to participate

The study received ethical approval from the Pharma-Ethics Research Ethics Committee of South Africa (Reference number: 13095660). All participants signed informed consent prior to enrolment.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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.

Department of Internal Medicine, Prince Mshiyeni Memorial Hospital and School of Clinical Medicine, College of Health Sciences, Nelson R Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa

School of Laboratory Medicine and Medical Sciences, Discipline of Medical Biochemistry and Chemical Pathology, University of KwaZulu-Natal, George Campbell Building – South Entrance, 3rd Floor, King George V Avenue, Howard College Campus, Durban, 4001, South Africa