The characteristics of participants
To explore the taxonomic composition of the gut microbiota from patients with psoriasis in comparison with healthy controls, 60 fecal samples were collected from 30 patients diagnosed with psoriasis and 30 healthy age-matched individuals for analysis. The healthy controls were composed of 20 males and 10 females, with a mean age of (43.7±13.21) years old and no history of psoriasis. The psoriasis patients were composed of 20 men and 10 women, with a mean age of (43.13±13.79) years old. The mean duration of the disease was (14.75±11.32) years. According to the medical history and physical examinations, 18 of the 30 patients with psoriasis had the symptoms of pruritus. The skin lesions were found to be more than 50% in 18 of the 30 patients (Table 1).
Table 1. Characteristics of patients with psoriasis (n = 30) and healthy individuals (n = 30) included in this study
|
Psoriasis
|
Healthy
|
P-value
|
Male/Female
|
20 / 10
|
20 / 10
|
1.000
|
Age (years)
|
43.13±13.79
|
43.7±13.21
|
0.871
|
Early-onset psoriasis ( < 40 years)
|
25 / 30
|
-
|
-
|
Duration of the disease (years)
|
14.75±11.32
|
-
|
-
|
Pustular psoriasis
|
6 / 30
|
-
|
-
|
Pruritus
|
18 / 30
|
-
|
-
|
Psoriasis area ≥ 50%
|
18 / 30
|
-
|
-
|
The biochemical examination of blood in the psoriasis group found that the mean level of white blood cell (WBC) was 7.55±0.52 ×109 L-1, of alanine transaminase (ALT) was (21.95±3.87) U L-1, of aspartate aminotransferase (AST), was (20.18±1.90) U L-1, of total cholesterol (TC), was (3.99±0.18) mmol L−1, of triglyceride (TG), was (1.62±0.19) mmol L−1, of high-density lipoprotein (HDL) was (1.02±0.05) mmol L−1, of low-density lipoprotein (LDL) was (2.68±0.15) mmol L−1, of fasting blood glucose (FBS) was (4.88±0.16) mmol L−1, and of hypersensitive c-reactive protein (hsCRP) was (11.88±4.72) mg L−1, of Immunoglobulin A (Ig A) was 2.54±0.26 g L-1, of Immunoglobulin E (Ig E) was 2.54±0.26 IU mL-1, of complement 3 (C3) was 1.09±0.03 g L-1, of interleukin-2R (IL-2R) was 678.61±106.33 U mL-1, of interleukin-6 (IL-6) was 10.91±4.59 pg mL-1, of interleukin-8 (IL-8) was 45.07±31.57 pg mL-1.
Fecal microbiota analysis
Bacterial DNA extracted from fecal samples was used for 16S rRNA sequencing analyses. 60 Samples Clean Reads generated a total of 7656 518 single-end reads with an average of 127608 reads per sample. A total of 635 OTUs were obtained, using Venn diagrams to visually show the number of shared and unique OTUs in different groups. The results showed that there were about 465 OTUs shared in the two groups, 102 and 68 OTUs in the healthy control group (N) and the psoriasis group (P), respectively (Figure 1A). Then the OTUs sequence and Greengenes database were annotated, and the result showed that predominant microbes were included (Figure 1B).
The OTU Abundance curves obtained for each sample tend to plateau, indicating that sample biodiversity was adequately covered with the applied sequencing depth (Supplement 1). The Alpha diversity box showed that there was no significant difference based on the Wilcox Test in the microbial diversity of patients with psoriasis (P), compared with the healthy controls (N). The sobs index is 140±50.824 and 150.1±45.883, P=0.34; Chao index is 162.91±56.14 and 176.75±45.51, P=0.21; ace index is 162.58±53.83 and 176.29±45.81, P=0.27; Shannon index is 2.55±0.63 and 2.56±0.65, P=0.92; simpson index is 0.18±0.11 and 0.20±0.15, P=0.81; coverage index is 1.00±0.01 and 1.00±0.01, P=0.34 in healthy controls and psoriasis patients respectively (Supplement 2 and 3).
Moreover, the β-diversity analyses revealed significant clustering (P-value < 0.001) between healthy individuals and patients with psoriasis, supporting the theory that the gut microbiota composition differed between both groups (Figure 2A and 2B).
Next, bacterial communities from both healthy individuals and psoriasis patients were analyzed. At the family and genus level, the microbial composition differed between the psoriasis patients and healthy individuals. We identified the presence of top 10 taxa in relative abundance at the family level, and the Wilcoxon test showed that Veillonellaceae and Ruminococcaceae were relatively more abundant in psoriasis patients. In contrast, Lachnospiraceae were relatively less abundant in psoriasis patients (all P<0.05). Moreover, at the genus level, Faecalibacterium and Megamonas were relatively more abundant in psoriasis patients than health control (all P<0.05) (Figure 2C and 2D). We then conducted a LEfSe comparison of the gut microbiota between the control and psoriasis groups to explore the specific bacterial taxa associated with psoriasis development. The structure and predominant bacteria of the microbiota in the control and psoriasis groups were represented in a cladogram (Figure 2E). The greatest difference in taxa from the phylum to the genus level was identified by an LDA score (Figure 2F). The two taxa with the highest score in the psoriasis group were Veillonellaceae and Megamonas, consistent with the previous one.
Correlation between gut microbiota and clinical indices in psoriasis patients
Correlation analyses were performed to investigate further identified potential associations of key microbial genera with clinical indices above (Figure 3). Exact figures of these clinical indices were showed in Supplement 4. Complement 3 (C3) showed a negative relationship with Bacteroides (P<0.001) and Escherichia (P<0.01), but it showed a positive relationship with the Prevotella group (P<0.01). Coprococcus had a great relationship with many clinical indices, and it showed a positive relationship with FBS (P<0.001) and IgA (P<0.001). Additionally, IL2R, which is abnormal in psoriasis patients, showed a positive relationship with Phascolarctobacterium (P<0.001) and a negative relationship with the Dialister group (P<0.001).
Different microbiota analysis between psoriasis vulgaris (PV) and pustular psoriasis (PP)
To compare the microbiota communities of PV (24 cases) and PP (6 cases). Based on the abundance ranks of the class level, the heat map was demonstrated through R (v3.1.1) gplots. Hierarchical clustering (Euclidean distance, complete linkage) showed that the PV group and healthy control group (N) tended to cluster together, and the PP group is being singled out (Figure 4A). Then we identified the presence of the top 10 taxa in relative abundance at the class level, and the Kruskal test showed that Gammaproteobacteria was relatively less abundant in the PP group than PV group and N group (all P-value<0.05) (Figure 4B). LefSe analysis was performed to identify differentially abundant taxa in the PP group compared with the PV group (Figure 4C). Consistent with our previous analysis, Gammaproteobacteria and Veillonellaceae had the largest LDA score indicating that Gammaproteobacteria and Veillonellaceae were consistently different between the PP group and PV group. In addition, the genus level Faecallibacterium and Anaerorhabdus were the most significant taxa in the PP group, which could help us distinguish PP from PV.