- 1Obstetrics and Gynaecology Department, Peking University First Hospital, Peking University, Beijing, China
- 2School of Basic Medical Sciences, Peking University, Ministry of Education, Beijing, China
Vulvar lichen sclerosus (VLS) is a chronic inflammatory skin disease that brings life-long and psychological distress to patients. It remains unclear whether this condition is related to changes in the skin microbial community. The aim of this study was to evaluate the compositional characteristics of the vulvar skin microbiota between VLS patients and asymptomatic postmenopausal women. We included 60 cases of postmenopausal patients in the outpatient vulvar clinic of Peking University First Hospital from August 2020 to October 2020. Thirty-one patients were diagnosed with VLS by vulvar skin biopsy (VLS group), while 29 women were asymptomatic volunteers (control group). DNA was extracted from vulvar skin swabs of the VLS and control groups. The V3-V4 fragments of 16S rRNA were targeted for high-throughput sequencing and gene sequence analysis. The sequencing results were analysed by α diversity, β diversity, species composition, LEfSe analysis to compare the compositional differences of the vulvar skin microbiota between the two groups. Our study revealed that at the phylum level, patients with VLS had a lower relative abundance of Firmicutes (p < 0.0001) and a higher relative abundance of Proteobacteria than the control group (p < 0.0001). At the genus level, Lactobacillus spp. accounted for the largest proportion of the microflora in the asymptomatic controls, while the proportion of Prevotella spp. in the VLS group was the highest. In the VLS group, the relative abundance of Finegoldia spp., Ralstonia spp., Peptoniphilus spp., Anaerococcus spp., Campylobacter spp., Providencia spp. Kelbsiella spp., Ezakiella spp., and Escherichia-Shigella spp. was significantly increased compared with the control group. Although there was no significant difference in the α diversity of the vulvar skin microbiota, the β diversity differed significantly between the two groups.
Introduction
Vulvar lichen sclerosus (VLS) is a chronic inflammatory skin disease of the vulva that may lead to vulvar scarring and sexual dysfunction. VLS affects premenarchal children and postmenopausal women, with an incidence rate ranging from 1 in 70 to 1 in 1,000 (Fruchter et al., 2017). However, the actual incidence may be underestimated. The pathogenesis of VLS remains uncertain. Autoimmune factors, genetic background, hormonal factors, and local factors have been associated with the pathogenesis of VLS. With the development of microbiome sequencing technology, it has been demonstrated that changes in the composition of the human microbiota may contribute to the development of chronic inflammatory diseases. Symbiotic bacteria can beneficially regulate host immunity, reducing the risk of infection-induced autoimmune disease and/or inflammation, while maladaptive dysbiosis of the microbiome may lead to disruption of immune homeostasis and disease. It has been reported that colonization and alteration of the skin microbiota are associated with a variety of skin diseases, such as dermatitis, rosacea, and psoriasis (Holmes et al., 2013; Nakatsuji et al., 2017; Thio et al., 2018; Wang et al., 2020). Currently, microbiome profiling has been used to explore the microbial composition of the vulvar skin and the role of the microbiome in vulvar health and disease (JayaramA et al., 2014; Pagan et al., 2021; Murina et al., 2020). However, few studies have focused on the correlation between vulvar skin microbiota and VLS. The aim of our study was to evaluate the compositional characteristics of the vulvar skin microbiota between VLS patients and asymptomatic women using 16S rRNA sequencing.
2 Materials and Methods
Patients and Sample Collection
Thirty-one VLS patients and 29 asymptomatic volunteers were recruited from the outpatient clinic of Peking University First Hospital from August 2020 to October 2020. The inclusion criteria for the VLS group were postmenopausal patients pathologically diagnosed with VLS by vulvar skin punch biopsy. The inclusion criteria for the control group were postmenopausal women without any vulvar symptoms or vulvar skin lesions. Women with acute vulvar infection, vaginitis, cervicitis or pelvic inflammation, use of topical glucocorticoids, oral or topical antibiotics within the last month, history of malignant tumours, radiotherapy/chemotherapy, or systemic skin diseases were excluded. Basic information (including body mass index, race and dietary habits) and medical history were recorded. The Cattaneo Clinical Score was used to evaluate the severity of VLS (Cattaneo A et al., 1991). Vaginal secretion wet pap was taken routinely to rule out bacterial vaginosis. Skin samples were collected from the labium majus and labium minus with a sterile swab in the clinic room. The sampling area was approximately 3 cm × 3 cm. Samples were stored in 2 ml sterile Eppendorf tubes at −80°C for DNA extraction.
Sequence Processing and Data Analysis
Genomic DNA was extracted using the bacterial Genome DNA extraction kit [Hi-Swab DNA Kit, TIANGEN ® BIOTECH (BEIJING) CO., LTD.]. The primer pairs used to amplify the hypervariable 16S rRNA regions V3-V4 were 341F (5′-CCTAYGGGRBGCASCAG-3′) upstream and 806R (5′-GGACTACNNGGGTATCTAAT-3′) downstream. The HiSeq 2,500 Illumina sequencing platform was used for high-throughput sequencing. High-quality sequences (length> 100 bp, mass fractions> 20) were clustered at a 97% similarity level with the Silva database into an operational taxon (OTU). The Shannon diversity index was used to compare the α diversity of the microbial communities between the two groups. We used the ACE and Chao1 indices to compare the richness of the microbiota and the Pielou index to compare the evenness. Richness and evenness were both measurements of α diversity. The Mann-Whitney test, principal coordinate analysis (PCoA) and ANOSIM were used to evaluate the differences in abundance and β diversity between the two groups. Quantitative Insights into Microbial Ecology (QIIME2) and the R package vegan were used for data processing and graph drawing.
Results
The Demographic Information of Vulvar Lichen Sclerosus Patients and Controls
The demographic information of the VLS patients and controls is shown in Table 1. All the patients and controls enrolled were Asian. The median Cattaneo Clinical Score of VLS patients was 8. There was no significant difference in age, allergy history, dietary habits, or daily hygiene habits between the two groups (p = 0.3103, 0.3022, 0.9999, and 0.7765 respectively), except for body mass index (BMI) (p = 0.0110).
Data Results and Phylogenetic Analyses (OTUs)
A total of 5,578,399 sequences were generated with an average of 92,973 sequences per sample. After quality filtering, a total of 5,537,018 high-quality sequences were obtained for downstream analysis. These sequences were divided into 10,868 OTUs, including VLS 4415 and controls 6,453. There were 679 overlapping OTUs between the two groups. Statistical data of identified OTU numbers at all levels in the VLS group and control group were shown in Table 2 and Figure 1.
Changes in the Structure of the Vulvar Microbiota in Patients With Vulvar Lichen Sclerosus
We found that there was no difference in the α diversity as measured by Shannon’s index (M-W test: p = 0.8882) (Figure 2A). The ACE and Chao1 richness indices in the VLS group were significantly lower than those in the control group (M-W test: p = 0.0069, p = 0.0033) (Figures 2B,C). The Pielou index for measurements of evenness was not significantly different between the two groups (M-W test: p = 0.2769) (Figure 2D). In terms of β diversity, which was used to compare the composition of taxa in the two cohorts, weighted UniFrac PCoA distinguished the VLS group from the asymptomatic control group (Figure 3A). The ANOSIM analysis showed R > 0, p = 0.001, indicating significant differences between the two groups (Figure 3B).
FIGURE 2. (A) Shannon diversity index (α-diversity) of the vulvar microbiota between the VLS and control groups. ACE index (B) and Chao1 index (C) (richness) of vulvar microbiota between the VLS and control groups. (D) Pielou index (D) (evenness) of vulvar microbiota between the VLS and control groups. *p values less than 0.05 were considered statistically significant.
FIGURE 3. (A) The results of principal coordinate analysis (PCoA) based on weighted UniFrac showed that the bacterial microbiome clustered separately in the VLS and control groups. (B) The plot of ANOSIM analysis measured the β diversity of vulvar microbiota between the VLS and control groups. (The Between column is the merged information of the two groups.)
Study on Changes in the Vulvar Microbiota in Vulvar Lichen Sclerosus Patients
According to the Silva analysis, the samples were divided into 45 phyla and 752 genera. The distribution of phyla and genera indicated differences between VLS patients and controls. The four most common phyla in both asymptomatic controls and VLS skin were Firmicutes, Bacteroidetes, Actinobacteria and Proteobacteria (Figure 4). Compared with asymptomatic controls, patients with VLS had an increased relative abundance of Proteobacteria (VLS vs. control, 19.22 vs. 2.66%, p < 0.0001) and a decreased relative abundance of Firmicutes (VLS vs. control, 43.47 vs. 68.67%, p < 0.0001) (Table 3). At the genus level, 9 key functional bacterial genera (Finegoldia spp., Ralstonia spp., Peptoniphilus spp., Anaerococcus spp., Campylobacter spp., Providencia spp. Kelbsiella spp., Ezakiella spp., and Escherichia-Shigella spp.) increased significantly between the VLS group and the control group; however, only 3 genera (Lactobacillus spp., Atopobium spp. and Gardnerella spp.) decreased significantly in VLS patients. Lactobacillus was predominant in the control group, and its relative abundance was significantly decreased in the VLS patients (54.17 vs. 9.73%, p < 0.0001). Prevotella spp. (belonging to Bacteroidetes) was predominant in the VLS, but there was no significant difference when compared with the control group (17.69 vs. 8.96%, p = 0.3630), which was consistent with the Bacteroidetes phyla (Table 4). Species’ relative abundance at the genus level of the VLS and control groups were shown in Figure 5. We also used the LEfSe algorithm to identify the specific taxa that were variably distributed in the groups. By using the default LDA cut-off of+/−4.0, we found that 16 taxa were overrepresented and 14 taxa were underrepresented in the VLS group compared to the control group. The key categories with significant differences between the VLS group and the control group are shown in Figure 6.
TABLE 3. The mean relative abundance of microbiota at phylum level in the VLS group and control group.
TABLE 4. The mean relative abundance of the microbiota at genus level in the VLS group and control group.
Discussion
Skin microecology is a skin microbiota composed of millions of bacteria, fungi, and viruses (Byrd et al., 2018). Like the microbes in the gut, skin microbes play a crucial role in the regulation of the immune system. The composition of the microbial community depended on the physiological state of the skin, and the relative abundance of bacteria changed with different humidity and sebaceous microenvironments. The vulva microbiota is diverse, and the normal vulva microbiota includes the characteristics of the vagina, urethra, and colon (Brown et al., 2007; Pagan et al., 2021). Previous studies on the vulvar skin microbiota, from cultivation to culture-independent molecular techniques, have shown that Lactobacillus and Corynebacterium spp., as well as Staphylococcus aureus and Prevotella, were the main bacterial components of healthy vulvar skin. (Elsner and Maibach, 1990a; Costello et al., 2009; Bruning et al., 2020). Wolf et al. (2017) transplanted vulvar skin of VLS patients onto the thigh and returned to normal skin, while normal skin was transplanted onto the vulvar area of VLS patients, and the skin turned into VLS. This suggested that the local microenvironment could be associated with the pathogenesis of VLS. Therefore, the natural microbiome of the vulvar skin may play a significant role in immunizing and maintaining vulva health, as well as preventing invasive organisms from gaining a foothold.
Lactobacillus is the permanent probiotic of the vaginal mucosa, which constitutes a majority of the vaginal flora. As the portal of the female reproductive tract, the flora composition of the vulva is influenced by the flora of the vagina. Therefore, lactobacilli dominate the vulva flora of healthy women (Brown et al., 2007; Costello et al., 2009; Murina et al., 2020). Researchers have used 16S rRNA sequencing to study the vulvar flora in premenopausal nonpregnant women. Samples were taken from 10 Japanese women before and on the second day of menstruation. Lactobacillus was dominant in all 7 subjects before menstruation and during the menstruation period (Shiraishi et al., 2011). In our study, we also found that the relative abundance of Lactobacillus in the asymptomatic control group was 54.17%, while the proportion of other bacterial groups was small. However, in the study of Chattopadhyay et al. (2021), sequencing and analysis of the vulva skin flora of 5 VLS child patients and 3 healthy children found that at the genus level, Porphyromonas spp., Prevotella, Corynebacterium and Peptoniphilus were dominant in the two groups, which may be due to the difference in vulva flora between women and children.
There was limited evidence on the skin flora of VLS patients. Watchorn et al. found a significant increase in the relative abundance of Prevotella compared to the normal control group in a study of male sclerosing lichens (20 vs. 4%) (Watchorn et al., 2021). Chattopahyay et al. (2021) found that the relative abundance of Prevotella in the vulvar flora of VLS children was significantly higher than that of healthy children. However, we found that there was no significant difference in the relative abundance of Prevotella between the VLS group and the control group (17.69 vs. 8.96%, p = 0.3630), which indicated that the vulvar skin microbiota of VLS patients might vary with age and postmenopausal status.
In addition to Prevotella, the abundance of other anaerobes in the vulva skin flora of VLS patients was also higher than that of healthy women. Guet-Revillet et al. (2014) found that strict anaerobes, such as Prevotella, Peptonophilus, and Porphyromonas gingivalis, were present in 87% of chronic suppurative bacteria, suggesting that these bacteria may be involved in chronic inflammation. In the same way, these anaerobic bacteria may also be involved in the formation of chronic skin inflammation and the development of VLS through similar immune regulation. Our study found that Firmicutes, Bacteroidetes, Actinobacteria, and Proteobacteria accounted for more than 97% of the total sequence in the vulvar skin swab samples in both groups. Compared to the control group, the number of Proteobacteria in the VLS group was significantly increased (19.22 vs. 2.66% p = 0.000997). Previous studies have shown that Corynebacterium predominate on wet skin, while the predominant colonization mixture is Proteobacteria and Flavobacteria on dry skin (Grice et al., 2009). Our result is consistent with the dry vulvar skin condition in patients with VLS.
Other potential factors influencing the vulvar skin microbiota include menstrual cycle, body mass index, pH, temperature, dryness, ethnicity, dietary habits, antibiotic use, and other factors, which play a crucial role in maintaining and stabilizing the microbiota. Vongsa et al. (2019) found that Lactobacillus spp. were more prevalent on the vulva of women with normal BMI, whereas Corynebacterium spp. and Anaerococcus spp. were more prevalent on the vulva of women with a high BMI. Therefore, it should be noted that our results may be influenced by BMI, as there were significant differences in BMI between the VLS patients and asymptomatic controls. It has been speculated that the vulva pH value may be somewhere between the skin value (estimated at pH 4.7) and the vagina value (average pH 3.5), which ranges from 3.8 to 4.2 in the menstrual cycle (Chen et al., 2017). A range of factors may affect the vulva pH value, including endogenous factors (e.g., humidity, sweat and vaginal discharge, menstrual cycle, urine and faecal pollution, anatomical folds, genetics, and age) and exogenous factors (such as soap, detergent, cosmetics, lubricants and spermicide, with tight clothes or sanitary napkin shade, shaving, and leather products). Long-term dryness of vulvar skin can significantly reduce its pH (Elsner and Maibach, 1990b), and pH can subsequently influence the growth of various microorganisms.
Therefore, our study has the following shortcomings. Due to the limitations of 16S rRNA, this study was not able to assess whether vulvar skin fungi and viruses changed in VLS patients. Second, questions have been raised regarding whether skin biopsy samples of the deep epidermis, dermis and glands might provide additional information on the skin microbiota (Grice et al., 2009). Therefore, swab sampling may have limitations in the research results. Finally, due to the limitations of the algorithm, our evaluation could not reach the species level, and metagenomics is required for further analysis in future work.
Conclusion
This study compared the compositional characteristics of the vulvar skin microbiota in patients with VLS and asymptomatic menopausal women. Although there was no significant difference in the α diversity, the β diversity differed significantly between the two groups. The predominant bacteria on the vulvar skin of asymptomatic women were Lactobacillus spp., while the proportion of Prevotella spp. was most prevalent in the VLS group. The relative abundance of Finegoldia spp., Ralstonia spp., Peptoniphilus spp., Anaerococcus spp., Campylobacter spp., Providencia spp. Kelbsiella spp., Ezakiella spp., and Escherichia-Shigella spp. was significantly increased in the VLS group. Metabolomics and transcriptomics are required to explain the relationship more completely between VLS and the vulvar skin microflora, which may provide a theoretical basis for clarifying the pathogenesis of VLS and new concepts for its treatment.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA791515.
Ethics Statement
The studies involving human participants were reviewed and approved by the Peking University First Hospital Institutional Review Board. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
JH and HW designed the study. XL analysed the data. CX performed the experimental validation. XL, YZ and YZ searched the literature and wrote the manuscript. All authors read and approved the manuscript.
Funding
This study was supported by Peking University First Hospital Youth Physician Foundation (No. 2015QN031).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We acknowledge HW for his guidance on this study.
References
Brown, C. J., Wong, M., Davis, C. C., Kanti, A., Zhou, X., and Forney, L. J. (2007). Preliminary Characterization of the normal Microbiota of the Human Vulva Using Cultivation-independent Methods. J. Med. Microbiol. 56 (Pt 2), 271–276. doi:10.1099/jmm.0.46607-0
Bruning, E., Chen, Y., McCue, K. A., Rubino, J. R., Wilkinson, J. E., and Brown, A. D. G. (2020). A 28 Day Clinical Assessment of a Lactic Acid-Containing Antimicrobial Intimate Gel Wash Formulation on Skin Tolerance and Impact on the Vulvar Microbiome. Antibiotics 9 (2), 55. doi:10.3390/antibiotics9020055
Byrd, A. L., Belkaid, Y., and Segre, J. A. (2018). The Human Skin Microbiome. Nat. Rev. Microbiol. 16 (3), 143–155. doi:10.1038/nrmicro.2017.157
Cattaneo, A., Bracco, G. L., Maestrini, G., Carli, P., Taddei, G. L., Colafranceschi, M., et al. (1991). Lichen Sclerosus and Squamous Hyperplasia of the Vulva. A Clinical Study of Medical Treatment. J. Reprod. Med. 36 (4), 301–305.
Chattopadhyay, S., Arnold, J. D., Malayil, L., Hittle, L., Mongodin, E. F., Marathe, K. S., et al. (2021). Potential Role of the Skin and Gut Microbiota in Premenarchal Vulvar Lichen Sclerosus: A Pilot Case-Control Study. PloS one 16 (1). e0245243. doi:10.1371/journal.pone.0245243
Chen, Y., Bruning, E., Rubino, J., and Eder, S. E. (2017). Role of Female Intimate hygiene in Vulvovaginal Health: Global hygiene Practices and Product Usage. Womens Health (Lond Engl. 13 (3), 58–67. doi:10.1177/1745505717731011
Costello, E. K., Lauber, C. L., Hamady, M., Fierer, N., Gordon, J. I., and Knight, R. (2009). Bacterial Community Variation in Human Body Habitats across Space and Time. Science 326 (5960), 1694–1697. doi:10.1126/science.1177486
Elsner, P., and Maibach, H. I. (1990a). Microbiology of Specialized Skin: the Vulva. Semin. Dermatol. 9 (4), 300–304.
Elsner, P., and Maibach, H. I. (1990b). The Effect of Prolonged Drying on Transepidermal Water Loss, Capacitance and pH of Human Vulvar and Forearm Skin. Acta Derm Venereol. 70 (2), 105–109.
Fruchter, R., Melnick, L., and Pomeranz, M. K. (2017). Lichenoid Vulvar Disease: A Review. Int. J. women's Dermatol. 3 (1), 58–64. doi:10.1016/j.ijwd.2017.02.017
Grice, E. A., Kong, H. H., Conlan, S., Deming, C. B., Davis, J., Young, A. C., et al. (2009). Topographical and Temporal Diversity of the Human Skin Microbiome. Science 324 (5931), 1190–1192. doi:10.1126/science.1171700
Guet-Revillet, H., Coignard-Biehler, H., Jais, J.-P., Quesne, G., Frapy, E., Poirée, S., et al. (2014). Bacterial Pathogens Associated with Hidradenitis Suppurativa, France. Emerg. Infect. Dis. 20 (12), 1990–1998. doi:10.3201/eid2012.140064
Holmes, A. D. (2013). Potential Role of Microorganisms in the Pathogenesis of Rosacea. J. Am. Acad. Dermatol. 69 (6), 1025–1032. doi:10.1016/j.jaad.2013.08.006
Jayaram, A., Witkin, S. S., Zhou, X., Brown, C. J., Rey, G. E., Linhares, I. M., et al. (2014). The Bacterial Microbiome in Paired Vaginal and Vestibular Samples from Women with Vulvar Vestibulitis Syndrome. Pathog. Dis. 72 (3), a–n. doi:10.1111/2049-632X.12197
Murina, F., Caimi, C., Di Pierro, F., Di Francesco, S., and Cetin, I. (2020). Features of the Vaginal and Vestibular Microbioma in Patients with Vestibulodynia: A Case-Control Study. J. Low Genit Tract Dis. 24 (3), 290–294. doi:10.1097/LGT.0000000000000523
Nakatsuji, T., Chen, T. H., Narala, S., Chun, K. A., Two, A. M., Yun, T., et al. (2017). Antimicrobials from Human Skin Commensal Bacteria Protect against Staphylococcus aureus and Are Deficient in Atopic Dermatitis. Sci. Transl Med. 9 (378), eaah4680. doi:10.1126/scitranslmed.aah4680
Pagan, L., Ederveen, R. A. M., Huisman, B. W., Schoones, J. W., Zwittink, R. D., Schuren, F. H. J., et al. (2021). The Human Vulvar Microbiome: A Systematic Review. Microorganisms, 9(12), 2568. doi:10.3390/microorganisms9122568
Shiraishi, T., Fukuda, K., Morotomi, N., Imamura, Y., Mishima, J., Imai, S., et al. (2011). Influence of Menstruation on the Microbiota of Healthy Women's Labia Minora as Analyzed Using a 16S rRNA Gene-Based Clone Library Method. Jpn. J. Infect. Dis. 64 (1), 76–80.
Thio, H. B. (2018). The Microbiome in Psoriasis and Psoriatic Arthritis: The Skin Perspective. J. Rheumatol. 94, 30–31. Supplement. doi:10.3899/jrheum.180133
Vongsa, R., Hoffman, D., Shepard, K., and Koenig, D. (2019). Comparative Study of Vulva and Abdominal Skin Microbiota of Healthy Females with High and Average BMI. BMC Microbiol. 19 (1), 16. doi:10.1186/s12866-019-1391-0
Wang, R., Farhat, M., Na, J., Li, R., and Wu, Y. (2020). Bacterial and Fungal Microbiome Characterization in Patients with Rosacea and Healthy Controls. Br. J. Dermatol. 183 (6), 1112–1114. doi:10.1111/bjd.19315
Watchorn, R. E., Munckhof, E. H. A., Quint, K. D., Eliahoo, J., Koning, M. N. C., Quint, W. G. V., et al. (2021). Balanopreputial Sac and Urine Microbiota in Patients with Male Genital Lichen Sclerosus. Int. J. Dermatol. 60 (2), 201–207. doi:10.1111/ijd.15252
Keywords: Vulvar lichen sclerosus, Skin microbiota, 16S rRNA, Lactobacillus, Alpha diversity, Beta diversity
Citation: Liu X, Zhuo Y, Zhou Y, Hu J, Wen H and Xiao C (2022) Analysis of the Vulvar Skin Microbiota in Asymptomatic Women and Patients With Vulvar Lichen Sclerosus Based on 16S rRNA Sequencing. Front. Cell Dev. Biol. 10:842031. doi: 10.3389/fcell.2022.842031
Received: 23 December 2021; Accepted: 11 March 2022;
Published: 05 April 2022.
Edited by:
Mario Antonio Bianchet, Johns Hopkins Medicine, United StatesReviewed by:
Mehmet Demirci, Kırklareli University, TurkeyAlan J. Wolfe, Loyola University Chicago, United States
Copyright © 2022 Liu, Zhuo, Zhou, Hu, Wen and Xiao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jun Hu, hujunjoyce@126.com
†These authors have contributed equally to this work and share second authorship