Characteristics of the Participants. Finally, 9 participants from HZ(Hangzhou city in Zhejiang Province, HZ) and 10 participants from YH (Yunhe city in Zhejiang Province, YH) were included for further analysis. The detailed characteristics of all individuals after wearing masks were shown in Table 1. Among the 19 participants, 5 (26.32%) participants reported acne occurrence after wearing masks.
Table 1
Demographic of participants
Characteristics | Findings |
Age (Mean ± SD) | 3 30.21 ± 4.78 |
Region | HZ(Hangzhou) 9 |
YH(Yunhe) 10 |
Skin type (According to sebaceous gland secretion) | Combination 6 31.58% |
Oily 4 21.05% |
Dry 1 5.26% |
Neutral 7 36.84% |
Acne occurred after wearing masks | Yes 5 26.32% |
No 14 73.68% |
Daily sleep time during wearing mask | 4–8 hours 13 68.42% |
8–12 hours 6 31.58% |
Whether or not to use medical protective mask, N95 | Yes 10 52.63% |
No 9 47.37% |
The average cumulative time of wearing masks per day | 3–6 hours 11 57.89% |
≥ 6 hours 8 42.11% |
Average time of wearing mask per time | < 2 hours 7 36.84% |
2–4 hours 4 21.05% |
≥ 4 hours 8 42.11% |
Whether or not to change mask when a mask was getting damp | Yes 14 73.68% |
No 5 26.32% |
Facial Microbial Diversity Changes Before and After Wearing masks. In this study, we used principal coordinate analysis (PCoA) to assess beta diversity metrics (intersample diversity). In the PCoA plot for genus, principal coordinate 1 (PC1) represented 46.5% of the variation, and principal coordinate 2 (PC2) represented 8.4% of the variation. In the PCoA plot for pathway, principlal coordinate 1 (PC1) represented 49.9% variation, and principal coordinate 2 (PC2) represented 23.9%. Each dot in Fig. 1a and 1b represented a sample from the facial oily (triangle) or dry (circle) zone of participants before (red) or after (green) wearing masks. Here, we observed that red dots and green dots separated from each other, which indicated that wearing masks changed the facial microbiota composition and pathways.
Differences in Taxonomic Composition of Facial Microbiota Before and After Wearing Masks. As shown in Fig. 2a-2c, the phyla of Proteobacteria (25.91%), Firmicutes (22.67%), Actinobacteria (10.74%), and Bacteroidetes (7.62%) still dominated the facial microbiomes, similar to the results before wearing masks in 20196. At the genus level after wearing masks in 2020, the top nine genera were Staphylococcus (18.24%), Acinetobacter (4.21%), Streptococcus (3.92%), Corynebacterium (3.72%), Enhydrobacter (3.36%), Pseudomonas (2.60%), Deinococcus (1.76%), Xanthomonas (1.67%) and Propionibacterium (1.60%), as shown in Fig. 2d-2f.
Then the changes in taxonomic profiles in paired samples were next investigated. We compared the relative abundance (logX + 1e-6) of the respective genera before or after wearing masks using the the paired Wilcoxon rank-sum test and drew a heatmap (Fig. 3). In paired samples before and after wearing masks, Faecalibacterium, Blautia, Bacteroides, and Bifidobacterium were significantly reduced after wearing masks, while Staphylococcus, Corynebacterium, Acinetobacter, Streptococcus, Corynebacterium, Enhydrobacter, Pseudomonas, Deinococcus and Xanthomonas were markedly increased (Supplementary Figure S1 a-b, and Figure S2 a-l). Notably, Faecalibacterium, listed as the top 1 before, dropped out of the top 10, while Staphylococcus established dominance as the No. 1 after wearing masks.
Effect of variables on microbiota composition and pathways. To assess the effect of factors of wearing masks as predictors of microbiota composition and pathways, we partitioned the explained variability (R2) attributable to each factor (relative importance analysis). For microbiota composition (Fig. 4a), the contributing variables (P < 0.05) included repeated use of masks, average duration per instance of wearing the mask or total duration per day of wearing masks, the amount of sleep, whether to replace a damp mask or use medical masks. Furthermore, the physiological state of individuals, including skin type, BMI (Body mass index, BMI), and age, also influences the facial microbiota composition. As shown in Fig. 4b, the average duration per time of wearing masks, age, sleep, whether to replace a damp mask, and mask time contributed to the microbiota pathway. It was worth noting that destination, which ever was the strongest factor affecting microbial composition in our previous study6, did not play a role in either microbiota composition or pathway after wearing masks.
Focusing on the top 20 genera of facial bacterial microbiota, we analyzed the relative abundance of specific genera among groups with different wearing habits. As shown in Fig. 5, the abundance of Prevotella was higher in the group with 8–12 hours of sleeping time than in the group with 4–8 hours of sleeping time per day (P = 0.0066, Fig. 5a), especially in the oily zone (P = 0.0357, Fig. 5b). Moreover, analysis of the top 20 pathways of microbiomes revealed a decreased expression PWY-3781[aerobic_respiration_I_(cytochrome_c)] in the group of 8–12 hours with sleeping time per day (P = 0.0036, Fig. 5d) compared with the 4–8 hours group.
In addition, we also compared the relative abundance of the top 10 genera of facial bacterial microbiota between the maskne (acne due to masks) group and the normal group (non-maskne). Figure 6a-d showed that the relative abundance of Staphylococcus, Propionibacterium and Pseudomonas in the maskne group was much higher than that in the non-maskne group, while that of Paracoccus in the maskne group decreased compared to the non-maskne group. Furthermore, the abundance of Staphylococcus and Propionibacterium displayed a positive relationship between them (R2 = 0.317).