Distinct gut microbiome observed in insomniac
The fecal microbiota composition profiles were analyzed by the 16S rRNA gene sequencing. A total of 2,480,664 high quality sequence reads were obtained and 31,008.3± 4,463.188 sequences for each sample. There were some differences in the composition of each sample at the phylum and genus levels. At the phylum level (Figure 1A), the relative abundance of Bacteroidetes in the gut microbiota of insomniacs and healthy controls was 54.60%±16.65% and 57.93%±20.21%, respectively; followed by Firmicutes and Proteobacteria; Firmicutes accounted for 33.64%±12.93% in insomniacs and 32.02%±18.05% in healthy controls, while Proteobacteria accounted for 10.62%±12.14% and 8.43%±9.98% in insomnia patients and healthy controls, respectively. At the genus level (Figure 1B), the relative abundance of Bacteroides in the gut microbiota of insomniacs and healthy controls was 31.60%±21.32% and 33.43%±22.76% respectively, followed by Prevotella, the relative abundance of Bacteroides was 19.51%±23.77% and 21.79% ±26.89% in the gut microbiota of the patients and controls, respectively. Venn diagram chart showed the core shared OTU (at least 50% of the samples in each group contain the OTU) in insomniacs and healthy controls. There were 692 core OTUs in healthy controls, 681 core OTUs in insomniacs and 611 total OTUs were founded in both two groups. However, there was no significant difference in alpha diversity between the two groups (P > 0.05) (Figure S1). The principal coordinate analysis (PCoA) was used to compare bacterial community patterns between the two groups. The results (Figure 1D) showed that there was a significant difference in community patterns between insomniacs and healthy controls based on the unweighted UniFrac distance (Adonis, R2=0.061, P =0.0001).
Linear discriminant effect size analysis (LEfSe) was performed with the purpose of screening potential gut microbiota biomarkers related to insomnia. A LDA histogram was used to represent the predominant bacteria and the structure of the microbiota in insomniac participants and healthy control participants (Figure 2A). We screened a total of 29 biomarkers, including 2 bacterial phyla, 5 classes, 5 orders, 9 families and 8 genera, by using a Linear Discriminant Analysis (LDA) value greater than 2 as the threshold (Table 2). These results are also illustrated in the cladogram (Figure 2B). As shown in Figure 2B, the abundance of Desulfovibrio, Lactobacillus and Streptococcus were significantly enriched in the gut microbiome of insomniac participants, whereas Bifidobacterium, Gardnerella, Sneathia, Aerococcus and Atopobium were more abundant in the healthy control participants.
Table 2. Different taxonomy in LEfSe analysis
OTU
|
Name
|
Group
|
LDA value
|
P value
|
Phylum
|
Actinobacteria
|
HC
|
2.88582817
|
8.22E-09
|
Fusobacteria
|
INS
|
2.423849947
|
2.92E-10
|
Class
|
Actinobacteria
|
HC
|
2.741136487
|
5.19E-09
|
Bacilli
|
INS
|
3.207435756
|
0.001144322
|
Coriobacteriia
|
HC
|
2.291295405
|
2.92E-06
|
Erysipelotrichi
|
INS
|
2.283274477
|
0.015724472
|
Fusobacteriia
|
INS
|
2.423849948
|
2.92E-10
|
Order
|
Bifidobacteriales
|
HC
|
2.741136487
|
5.19E-09
|
Coriobacteriales
|
HC
|
2.291295404
|
2.92E-06
|
Erysipelotrichales
|
INS
|
2.283274477
|
0.015724472
|
Fusobacteriales
|
INS
|
2.423849947
|
2.92E-10
|
Lactobacillales
|
INS
|
3.25087112
|
0.000965062
|
Family
|
Aerococcaceae
|
HC
|
2.186764193
|
5.21E-08
|
Bifidobacteriaceae
|
HC
|
2.741136487
|
5.19E-09
|
Coriobacteriaceae
|
HC
|
2.291295404
|
2.92E-06
|
Erysipelotrichaceae
|
INS
|
2.283274477
|
0.015724472
|
Lachnospiraceae
|
INS
|
3.59373051
|
0.004957671
|
Lactobacillaceae
|
INS
|
2.980629362
|
0.002075563
|
Leptotrichiaceae
|
HC
|
2.291818463
|
2.01E-14
|
Streptococcaceae
|
INS
|
2.981593472
|
5.20E-05
|
Aerococcaceae
|
HC
|
2.186764193
|
5.21E-08
|
Genus
|
Aerococcus
|
HC
|
2.186764194
|
5.21E-08
|
Atopobium
|
HC
|
2.156748523
|
1.41E-05
|
Bifidobacterium
|
HC
|
2.838232315
|
3.51E-08
|
Desulfovibrio
|
INS
|
2.403106251
|
0.000188464
|
Gardnerella
|
HC
|
2.594522952
|
6.09E-10
|
Lactobacillus
|
INS
|
2.980629362
|
0.002075563
|
Sneathia
|
HC
|
2.291818463
|
2.01E-14
|
Streptococcus
|
INS
|
2.981593472
|
5.20E-05
|
Variation of serum inflammatory factors in insomniacs
During sleep, some unique substances are produced by the human body. We wanted to study these to understand difference between observations in the control group observations obtained from the insomniac group. Certain inflammatory factors have also been shown to regulate and affect the sleep in humans19. However, inflammation and infection could also alter sleep architecture, whereas a lack of sleep could also impair immune function. Therefore, we measured inflammatory factors in the serum such as interleukin-1β (IL-1β), tumor necrosis factor-α (TNF-α), and interleukin 6 (IL-6) and these measurements were obtained by using ELISA. As Figure 3 indicates, the IL-1β was significantly elevated (2.315±2.202 pg/ml in insomniac group participants VS 0.967±0.745 pg/ml in control group participants, P ˂0.0001), meanwhile TNF-α was significantly decreased in the insomniac participants (2.055±1.619 pg/ml in insomniacs VS 3.234±1.520 pg/ml in controls, P ˂0.0001). However, the IL-6 level was consistent between the patients and healthy controls (1.842±1.396 pg/ml in insomniac vs 1.805±1.541 pg/ml in controls).
Correlation between serum inflammatory factors and gut microbiota
To reveal the potential correlation between the enterobacteria and the inflammatory factors in the insomniacs, the Pearson correlation analysis was used (Figure 4). We found that 4 OTUs of Bacteroides, 2 OTUs of Bacteroides fragilis, 1 OTU of Oscillospira has a significant positive correlation with IL-1β. And IL-6 positively correlated with 5 OTUs of Prevotella copri, 2 OTUs Faecalibacterium prausnitzii, 2 OTUs of Bacteroides, and 1 OTU of Anaerostipes. At the same time, TNF-α revealed a significant positive correlation with the OTU of Prevotella copri, Parabacteroides, Oscillospira, Butyricimonas and Bacteroides.
Significant alteration of metabolomic profiles in insomniacs
Metabolic profiling of the insomniacs and healthy controls were acquired by Ultra performance liquid chromatography-tandem mass spectrometer (UHPLC-MS/MS). 5831 peaks were detected, and 5736 peaks remained after using a relative standard deviation de-noising method in positive ion mode (ES+). And in negative ion mode (ES-), 2612 peaks were detected, and 2588 peaks remained after de-noising. To discriminate the metabolic profiles between insomniacs and healthy control groups, we used principal component analysis (PCA) and orthogonal partial least-squares discriminant analysis (OPLS-DA). The results indicated a large separation of metabolome between insomniacs and healthy controls in PCA plots (Figure 5A, B) both in ES+ (R2X=0.32) and ES- (R2X=0.317). Furthermore, individuals in insomniac group were separated from and healthy control group as further evidenced by the OPLS-DA score scatter plots both in ES+ (R2X=0.169, R2Y=0.905, Q2=0.707) and ES- (R2X=0.175, R2Y=0.84, Q2=0.532). We performed a permutation test in order to further validate the OPLS-DA model. After 200 permutations, the R2 intercept was 0.71 and 0.78 in ES+ and ES- respectively, the Q2 intercept values were -0.87 and -0.73 in ES+ and ES- respectively (Figure S2A, B). In order to screen for differential metabolites, the first principal component of variable importance in the projection (VIP) was obtained (Figure S2C, D). The VIP values exceeding 1 were first selected as differential metabolites, and the Student’s t-test also be used. Then we found 25 metabolites in ES+ and 6 metabolites in ES- (VIP >1,P<0.05, Table S2). We also visualized it with a heatmap (Figure 5E). We found that, in ES+, 1-palmitoylglycerophosphocholine was identified significantly increased, however, 7-hydroxy-6-methoxy-alpha-pyrufuran, acinospesigenin A, PS(22:2(13Z,16Z)/20:2(11Z,14Z)) and PA(17:1(9Z)/13:0) were significantly decreased in the insomniacs’ serum. Meanwhile, aspartic acid, phenylalanine and phosphatidylcholine lyso 20:4 were identified to be significantly decreased in the insomniacs’ serum when under ES-. All differential metabolites were then subjected to the regulatory pathways analysis to discover the metabolic pathways exhibiting high correlations with the metabolites. According to P value and influence value, significant abnormalities were found in five metabolic pathways in insomniacs: glycerophospholipid metabolism (P <0.05), glutathione metabolism (P <0.05), nitrogen metabolism (P <0.01), alanine, aspartate and glutamate metabolism (P <0.05), and aminoacyl-tRNA biosynthesis (P < 0.05) (Figure 5F, G).
Correlation between serum metabolome and gut microbiota
To explore the functional relationship of the altered gut microbiota and metabolites in insomniacs, the correlation matrixes based on Pearson's rank correlation coefficient were formulated (Figure 6). Correlations between the varied metabolites and the gut microbiota were identified. The 1-palmitoyl lysophosphatidic acid was positively correlated with Veillonella dispar. Phenylalanine was positively correlated with Streptococcus. Phosphatidylcholine lyso 20:4 was positively correlated with [Prevotella], [Ruminococcus] gnavus, [Eubacterium] biforme, Bacteroides, Fusobacterium and Streptococcus. What’s more, L-5-oxoproline showed a significant negatively correlation with Atopobium. And some short chain polypeptides (Ala-Arg-Arg-Asn, Ala-Trp-Arg-Lys, Asn-Gly-Val, Asp-Phe,Phe-Phe, Val-Phe-Arg) were also negatively correlated with Atopobium.