- 1Department of Anorectal Surgery, Qilu Hospital of Shandong University (Qingdao), Shandong, China
- 2State Key Laboratory of Microbial Technology, Microbial Technology Institute, Shandong University, Jinan, Qingdao, China
Hemorrhoids are a common ailment that can cause significant disruptions to one’s daily life. While some researchers have speculated about a potential link between hemorrhoid development and gut microbes, there is currently insufficient evidence to support this claim. In this study, we collected samples from 60 hemorrhoid patients and analyzed the composition and characteristics of microbiomes in hemorrhoids. PCoA results revealed distinct differences between the microbiomes of hemorrhoids, skin-originated microbiomes, and gut microbes, highlighting the complex nature of hemorrhoidal microbiomes. The distribution characteristics of Staphylococcus suggest that the skin microbiome influences the microbiome of hemorrhoids. Additionally, we observed higher levels of Prevotella in two cases of thrombosed hemorrhoids compared to non-thrombosed hemorrhoids. This finding suggests that Prevotella may play a crucial role in the development of thrombosed hemorrhoids.
Introduction
Hemorrhoids are the most common condition in the field of proctology, with a prevalence ranging from 4 to 38.93% (Johanson and Sonnenberg, 1990; Riss et al., 2012). This condition typically involves the displacement of the anal cushion, leading to discomfort and impacting the daily lives of those affected (Loder et al., 1994). Hemorrhoids can be classified into internal hemorrhoids and external hemorrhoids based on their location and clinical manifestations (Adler, 1905). In addition, thrombosed hemorrhoids are a specific type of hemorrhoids characterized by the formation of blood clots within the hemorrhoidal blood vessels. This condition results in swelling, redness, and pain in the hemorrhoidal tissues (Zhang et al., 2019).
In the medical field, there is no consensus on the etiology of hemorrhoids. Currently, there are several main theories on the causes of hemorrhoids: (1) the varicose vein theory, (2) the vascular hyperplasia theory, and (3) the theory of sliding anal canal lining, which is widely accepted (Palumbo et al., 2023). The earliest theory proposed that hemorrhoids may be related to varicose veins, such as varicose swelling caused by increased local pressure from hard stools (Burkitt, 1972; Chawla and Dilawari, 1989). Some scholars believe that the formation of hemorrhoids may be a result of vascular proliferation during the inflammatory response and tissue remodeling process (Chung et al., 2004; Aigner et al., 2009). The latest theory suggests that the pathological sliding of the anal cushions is the main cause of hemorrhoid formation (Thomson, 1975; Haas et al., 1984; Loder et al., 1994; Lohsiriwat, 2012). The above three theories are based on the characteristics of human tissues. In addition to the reasons mentioned above, dietary habits such as a low-fiber diet, lack of regular exercise, and low water intake are also believed to be associated with the formation of hemorrhoids (Feyen et al., 2022).
In the past decade, with the rapid advancement of microbiology, people have gradually realized that the role of microorganisms is far more powerful than we imagined. The traditional concept of the gut-brain axis has also gradually transformed into the gut-brain-microbiome (GBM) interactions. Many researchers believe that there is a close relationship between microorganisms and the development of hemorrhoids, particularly the microorganisms residing in the gut. Current research shows that the intestinal microbiota plays a crucial role in the development of the host immune system, digestion, and other pathways. Moreover, an increasing number of studies have demonstrated a correlation between fecal microbiota and chronic constipation (Jeffery et al., 2012; Parthasarathy et al., 2016).
With the popularity of high-throughput sequencing technology, the 16S ribosomal RNA (rRNA) gene amplicon sequencing has also been widely used in various gastrointestinal diseases (Takahashi and Andoh, 2016; Nishino et al., 2018). Compared with traditional microbiological culture, sequencing technology can accurately and realistically reflect the composition and proportion of microbial communities with high throughput (Lundberg et al., 2013). Furthermore, based on the viewpoint that the majority of microorganisms in different environments are unculturable, 16S rRNA can accurately represent the microbial composition in a given environment (Steen et al., 2019). From the above considerations, in the study on the correlation between the clinical occurrence of hemorrhoids and microbiota, in order to find the relationship between microbe and hemorrhoids, 16S rDNA amplicons sequencing was used to analyze the samples of 60 patients with hemorrhoids.
Materials and methods
Sample collection and preparation
The samples were taken from various regions of patients with hemorrhoids in the anorectal department of Qilu Hospital of Shandong University (Qingdao). When sampling, use a sterile cotton swab along with an anal reamer to carefully collect an adequate amount of tissue fluid from the relevant sites. The swabs are then quickly placed in sterile Eppendorf tubes containing 15% glycerin, stored at a low temperature, and transported to the lab.
16S rDNA amplicon sequencing and bioinformatics analysis
Total DNA was extracted from swab samples using the DNeasy® PowerSoil® Kit (Qiagen, Germany). The extracted DNA was checked using the NanoDrop One spectrophotometer (NanoDrop Technologies, Wilmington, DE) and the Qubit 3.0 Fluorometer (Life Technologies, Carlsbad, CA, United States). Amplification of the 16S V3-V4 region was performed using a custom primer with a barcode (341F: 5′-CCTACGGGNGGCWGCAG-3′, 805R: 5′-GACTACHVGGGTATCTAATCC-3′). Utilizing the library building kit from Vazyme, specifically the VAHTS® Universal DNA Library Prep Kit for Illumina. 16S rDNA amplicon sequencing was performed on the Illumina NovaSeq 6000 sequencer (Illumina, Inc., San Diego, CA, USA). Using QIIME2 v2021.4.0, the primer sequences were removed (Bokulich et al., 2018). Using the QIIME2 plugin vsearch (with a default sequence threshold of 97%), Operational Taxonomic Units (OTUs) classification was determined. The SILVA ribosomal RNA database (v138.1) was used for annotation.
Data analysis and statistics
Unweighted UniFrac distance was calculated using the phyloseq package in the R platform. Hierarchical clustering was performed using the UPGMA method. MaAsLin2 (Microbiome Multivariable Association with Linear Models version 2.0) was used to identify biomarkers of microbiomes using linear models. ANOSIM and NMDS were performed using the vegan package on the R platform.
All statistical analyses were performed using Jamovi 2.3.21.0 or Prism 9.4.1 (681).
Results and discussion
Sample collection and experimental design
Samples were collected from 60 hemorrhoid-suffering patients to investigate the microbiomic signatures of hemorrhoids and the potential origins of microbes associated with hemorrhoids. These patients include 22 males and 38 females, with an average age of 43.75. Standard diagnostic procedures were performed to confirm the disease, and samples were taken before and after surgery (for hemorrhoids).
Four samples were taken from each patient, including swabs of feces, rectum surface, hemorrhoid surface before surgery, and hemorrhoid after surgery (Figure 1). These sampling sites were designed to identify potential microbiomic biomarkers for each site, especially for hemorrhoids, and to evaluate the potential origins of the hemorrhoidal microbes.
Bacterial compositions of collected samples
The bacterial compositions of the collected samples were determined using 16S rDNA amplicon sequencing. On the phylum level, Firmicutes, Proteobacteria, and Bacteroidota are the three most abundant bacterial taxa, consistent with prior research on human gastrointestinal microbiomes (Figure 2A). On the genus level, Prevotella, Escherichia-Shigella, Bacteroides, Faecalibacterium, and Streptococcus, are present at high abundances. These genera are also frequently observed in gastrointestinal tract samples (Figure 2B).
Figure 2. Microbial community structures of the four sampling locations. (A), on phylum level; (B), on genus level.
Analysis of alpha diversities revealed that the microbiomes of the hemorrhoid surface and rectum surface exhibit similar alpha diversities (Figure 3; Table 1). Similarly, feces and hemorrhoid microbiomes exhibit similar α diversities (Figure 3; Table 1). However, both feces and hemorrhoid microbiomes exhibit lower richness (as indicated by Chao1 indices and ACE estimators) and diversity (as measured by Shannon indices and Simpson indices) levels. Therefore, feces and hemorrhoid carry microbiomes with lower diversities, which can be explained by the lack of exposure of hemorrhoid and feces to the environment. This is particularly evident when considering that the skin has been removed from the hemorrhoids. However, it is surprising to discover that hemorrhoid microbiomes are as diverse as fecal microbiomes, which could be one of the most diverse bacterial communities inside the human body. This reveals the complexity of microbes in hemorrhoids.
An analysis of similarities between microbiomes was performed to explore the relationships among microbiomes from various sampling sites. Weighted UniFrac distances were calculated between each microbiome using a maximum likelihood phylogenetic tree constructed with FastTree. PCoA analysis revealed that the microbiomes at the four sampling sites are significantly different (Figure 4, ANOSIM p = 0.001). However, the microbiomes are “paired” in nature, indicating that the four microbiomes from four sampling sites of the same subject form a distinct group of samples. Therefore, simply grouping microbiomes of the same sampling site from different samples is statistically flawed. With this consideration, a second approach was subsequently applied to assess the distance between microbiomes of every pair of sampling sites and to conduct paired t-tests to identify significant differences in the distances. Clear differences in distances between each pair of sampling sites can be observed, as shown in Figure 5A. Subsequent statistical comparison (Figure 5A) showed that the distances between the microbiomes of each pair of sampling sites follow a clear hierarchy: hemorrhoid surface and hemorrhoid = rectum surface to hemorrhoid = feces to hemorrhoid > hemorrhoid surface to feces > hemorrhoid surface to rectum surface > rectum surface to feces.
Figure 5. Weighted UniFrac distances between different pairs of sampling sites (A) and their significances (B).
An interesting observation in this distance hierarchy is that the hemorrhoid microbiome is equally distant from other microbiomes, which are closer to each other than to the hemorrhoid microbiome. Although PCoA of averaged distances cannot clearly reflect this observation because it is essentially a dimension reduction technique and cannot accurately represent true distances on a two-dimensional plot (Figure 6), as shown in the three-dimensional illustration of the four microbiomes in Figure 6, it is evident that hemorrhoid microbiomes are distinct from all other three sampling sites. It can therefore be suggested that hemorrhoids harbor a significantly different microbiome compared to all other sampling sites. Considering that all three other sampling sites are part of the gastrointestinal tract and closely interact with the gastrointestinal microbiomes, this finding strongly suggests that hemorrhoid microbiomes may have sources of influence beyond the gastrointestinal microbiomes.
Bacterial signatures of hemorrhoid microbiomes
To understand the signatures of hemorrhoid microbiomes, a linear model-based MaAsLin2 algorithm was used to identify biomarkers for hemorrhoid microbiomes (Mallick et al., 2021). This method outperforms LEfSe, a similar algorithm, by reducing false positive discoveries. Bacterial genera with an average abundance of over 0.1% and a discovery-adjusted p-value (FDR) of <0.05 were reported. No significant microbiomic signature was found in hemorrhoid microbiomes compared to rectum surface microbiomes. However, six genera were found to distinguish the microbiomes of hemorrhoid surfaces, and three genera were found to distinguish the microbiomes of feces compared to hemorrhoid microbiomes. Distinct patterns were found while analyzing the microbiomic signatures of hemorrhoid surface and feces compared to hemorrhoid microbiomes. All six genera identified had higher abundances in hemorrhoid surface samples than in hemorrhoid samples (Table 2), while two out of three bacterial groups identified in the hemorrhoid-feces comparison had higher abundances in hemorrhoid samples (Table 3).
The higher abundances of bacterial genera on the surface of hemorrhoids compared to samples from inside the hemorrhoids suggest a bacterial gradient from the surface to the interior of hemorrhoids. This is interestingly true for both feces and rectum surface samples: all six genera are present in significantly higher abundances in hemorrhoid surface samples when analyzed with MaAsLin2. For Actinomyces, the average abundances in the rectum surface and feces are 1.11 × 10−3 and 1.17 × 10−3, respectively, significantly lower than the abundance of 5.24 × 10−3 in hemorrhoid surface (FDR = 8.32 × 10−6 and 1.33 × 10−7, respectively). For Corynebacterium, the average abundances in the rectum surface and feces are 2.13 × 10−3 and 2.89 × 10−3, respectively, significantly lower than the abundance of 0.018 in the hemorrhoid surface (FDR = 2.31 × 10−5 and 4.13 × 10−9, respectively). For Finegoldia, the average abundances in the rectum surface and feces are 0.012 and 5.56 × 10−3, respectively, significantly lower than 0.046 in the hemorrhoid surface (FDR = 4.04 × 10−3 and 4.29 × 10−6, respectively). For Anaerococcus, the average abundances in the rectum surface and feces are 9.93 × 10−3 and 6.62 × 10−3, respectively, significantly lower than the abundance of 0.032 on the hemorrhoid surface (FDR = 9.55 × 10−3 and 1.12 × 10−5, respectively). For Peptoniphilus, the average abundances in the rectum surface and feces are 0.011 and 0.011, respectively, which are significantly lower than the abundance of 0.034 in the hemorrhoid surface (FDR = 0.011 and 1.63 × 10−4, respectively). For Staphylococcus, the average abundances in the rectum surface and feces are 2.18 × 10−3 and 3.39 × 10−3, respectively, significantly lower than the abundance of 0.019 on the hemorrhoid surface (FDR = 1.29 × 10−5 and 1.63 × 10−4, respectively). Staphylococcus is a genus particularly associated with the skin and mucosa. The high presence of microbiomes on hemorrhoid surfaces suggests that they originate from the skin.
The hemorrhoid microbiome differs from the fecal microbiome in three genera, including Staphylococcus. A higher abundance of Staphylococcus was found in hemorrhoids (7.68 × 10−3) compared to feces (3.39 × 10−3) and the rectum surface (2.18 × 10−3). This, in combination with the higher presence of Staphylococcus in hemorrhoid samples and the skin origin of this group, appears to indicate a skin origin of Staphylococcus in hemorrhoid. Although not studied in this work, previous attempts to analyze fecal microbiomes of healthy adults (Feng et al., 2023; Liu et al., 2023) and children (Yuan et al., 2023) were unable to identify characteristic skin bacteria in their feces. Therefore, a proposal can be made that skin microbiomes at least partly influence the microbiomes in hemorrhoids. It must be stressed here that we did not carry out similar comparisons in healthy subjects. Therefore, this proposal was not validated using healthy controls to determine the absence of skin-oriented microbes in them. Hence, conclusions need to be carefully evaluated.
Microbes in thrombosed hemorrhoid
Two hemorrhoid samples we were able to obtain are from patients suffering from thrombosed hemorrhoids. This is a type of hemorrhoid that can lead to acute pain and swelling (Halverson, 2007). There has been no report or proposal connecting this form of hemorrhoids to microbial infections. The microbiomes associated with the two thrombosed hemorrhoid samples were examined and compared with other hemorrhoid samples. Although the microbiomes in thrombosed hemorrhoid samples Z31 and Z60 appear vastly different from averaged hemorrhoid samples, the most striking observation is that both thrombosed hemorrhoid samples contain a very high level of Prevotella (Figure 7; Z31, 46.37%; Z60, 35.45%; averaged hemorrhoid samples, 9.80%). Due to the rarity of patients with thrombosed hemorrhoids, we were unable to collect additional samples for a more comprehensive investigation. Nevertheless, this observation of a high Prevotella level implies that Prevotella may play a role in the development of thrombosed hemorrhoids, and that the blood clots within thrombosed hemorrhoids do contain microbes. Indeed, preliminary screening successfully isolated bacteria from thrombosed hemorrhoids. However, Prevotella strains were not cultivated due to our lack of awareness regarding the potential enrichment of this genus. Additionally, we did not utilize media that could enrich this genus or employ anaerobic cultivating conditions.
Conclusion
In conclusion, the microbiomic signatures of hemorrhoids and related samples were investigated in 60 patients. The microbiomes of hemorrhoids are significantly different from those found on hemorrhoid surfaces, rectum surfaces, and feces. Further analysis and comparison of microbiomic biomarkers suggest that microbiomes in hemorrhoids are influenced significantly by the surfaces of the hemorrhoids. Specifically, Staphylococcus, which are bacteria associated with the skin or mucosa, are more abundant on hemorrhoid surfaces, less abundant in hemorrhoids, and even less so on rectum surfaces and fecal samples. This implies a potential transmission of skin microbes into hemorrhoids. Further investigation of thrombosed hemorrhoid samples suggests a higher Prevotella content, indicating the presence of microbes in thrombosed hemorrhoids. Therefore, Prevotella may play a role in the development of thrombosed hemorrhoids.
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 in the article/supplementary material.
Ethics statement
The studies involving humans were approved by Medical Ethics Committee of Qilu Hospital of Shandong University (Qingdao) under approval number KYLL-2023045. The studies were conducted in accordance with the local legislation and institutional requirements. The 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
YQW: Data curation, Writing – original draft. WYS: Formal analysis, Writing – original draft. ZL: Writing – original draft & Software. YHW: Data curation, Writing – original draft. LL: Writing – original draft. HX: Writing – review & editing. MW: Funding acquisition, Writing – review & editing. WLS: Funding acquisition, Resources, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Foundation of Qingdao Key Health Discipline Development Fund under grant number QDZDZK-2022098, the National Key Research and Development Program of China under grant number 2022YFE0199800; Key R&D Program of Shandong Province under grant number 2020CXGC011305; and the National Natural Science Foundation of China under grant number 82271658.
Acknowledgments
We would like to thank Wuhan Benagen Technology Co., Ltd. for conducting the16S rDNA amplicon sequencing.
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.
References
Adler, L. H. (1905). Internal and external hemorrhoids—with special reference to their treatment. J. Am. Med. Assoc. XLIV, 193–196. doi: 10.1001/jama.1905.92500300025001g
Aigner, F., Gruber, H., Conrad, F., Eder, J., Wedel, T., Zelger, B., et al. (2009). Revised morphology and hemodynamics of the anorectal vascular plexus: impact on the course of hemorrhoidal disease. Int. J. Color. Dis. 24, 105–113. doi: 10.1007/s00384-008-0572-3
Bokulich, N. A., Kaehler, B. D., Rideout, J. R., Dillon, M., Bolyen, E., Knight, R., et al. (2018). Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome 6:90. doi: 10.1186/s40168-018-0470-z
Burkitt, D. P. (1972). Varicose veins, deep vein-thrombosis, and hemorrhoids-epidemiology and suggested etiology. Brit. Med. J. 2, 556–561. doi: 10.1136/bmj.2.5813.556
Chawla, Y. K., and Dilawari, J. B. (1989). Anorectal varices, hemorrhoids, and portal-hypertension. Lancet 1:725
Chung, Y. C., Hou, Y. C., and Pan, A. C. H. (2004). Endoglin (CD105) expression in the development of haemorrhoids. Europ. J. Clin. Invest. 34, 107–112. doi: 10.1111/j.1365-2362.2004.01305.x
Feng, J., Wu, Y., Dai, P., Wang, D., Liu, L., and Chai, B. (2023). Gut microbial signatures of patients with primary hepatocellular carcinoma and their healthy first-degree relatives. J. Appl. Microbiol. 134:lxad221. doi: 10.1093/jambio/lxad221
Feyen, L., Freyhardt, P., Schott, P., Wullstein, C., Dimitriou-Zarra, A., and Katoh, M. (2022). Embolization of hemorrhoidal arteries: a novel endovascular therapy option for the treatment of hemorrhoidal disease. RöFo 194, 266–271. doi: 10.1055/a-1662-5487
Haas, P. A., Fox, T. A. J., and Haas, G. P. (1984). The pathogenesis of hemorrhoids. Dis. Colon Rectum 27, 442–450. doi: 10.1007/bf02555533
Jeffery, I. B., O'Toole, P. W., Ohman, L., Claesson, M. J., Deane, J., Quigley, E. M. M., et al. (2012). An irritable bowel syndrome subtype defined by species-specific alterations in faecal microbiota. Gut 61, 997–1006. doi: 10.1136/gutjnl-2011-301501
Johanson, J. F., and Sonnenberg, A. (1990). The prevalence of hemorrhoids and chronic constipation. Gastroenterology 98, 380–386. doi: 10.1016/0016-5085(90)90828-o
Liu, P., Yang, J., Chen, Y., Zhu, Y., Tang, Y., Xu, X., et al. (2023). Alterations of gut microbiota and metabolome in early chronic kidney disease patients complicated with hyperuricemia. Heliyon 9:e20328. doi: 10.1016/j.heliyon.2023.e20328
Loder, P. B., Kamm, M. A., Nicholls, R. J., and Phillips, R. K. S. (1994). Hemorrhoids-pathology, pathophysiology and etiology. Brit. J. Surg. 81, 946–954. doi: 10.1002/bjs.1800810707
Lohsiriwat, V. (2012). Hemorrhoids: from basic pathophysiology to clinical management. World J. Gastroenterol. 18, 2009–2017. doi: 10.3748/wjg.v18.i17.2009
Lundberg, D. S., Yourstone, S., Mieczkowski, P., Jones, C. D., and Dangl, J. L. (2013). Practical innovations for high-throughput amplicon sequencing. Nat. Methods 10, 999–1002. doi: 10.1038/nmeth.2634
Mallick, H., Rahnavard, A., McIver, L. J., Ma, S., Zhang, Y., Nguyen, L. H., et al. (2021). Multivariable association discovery in population-scale meta-omics studies. PLoS Comput. Biol. 17:e1009442. doi: 10.1371/journal.pcbi.1009442
Nishino, K., Nishida, A., Inoue, R., Kawada, Y., Ohno, M., Sakai, S., et al. (2018). Analysis of endoscopic brush samples identified mucosa-associated dysbiosis in inflammatory bowel disease. J. Gastroenterol. 53, 95–106. doi: 10.1007/s00535-017-1384-4
Palumbo, V. D., Tutino, R., Messina, M., Santarelli, M., Nigro, C., Lo Secco, G., et al. (2023). Altered gut microbic flora and haemorrhoids: could they have a possible relationship? J. Clin. Med. 12:2198. doi: 10.3390/jcm12062198
Parthasarathy, G., Chen, J., Chen, X. F., Chia, N., O'Connor, H. M., Wolf, P. G., et al. (2016). Relationship between microbiota of the colonic mucosa vs feces and symptoms, colonic transit, and methane production in female patients with chronic constipation. Gastroenterology 150, 367–379.e1. doi: 10.1053/j.gastro.2015.10.005
Riss, S., Weiser, F. A., Schwameis, K., Riss, T., Mittlbock, M., Steiner, G., et al. (2012). The prevalence of hemorrhoids in adults. Int. J. Color. Dis. 27, 215–220. doi: 10.1007/s00384-011-1316-3
Steen, A. D., Crits-Christoph, A., Carini, P., DeAngelis, K. M., Fierer, N., Lloyd, K. G., et al. (2019). High proportions of bacteria and archaea across most biomes remain uncultured. ISME J. 13, 3126–3130. doi: 10.1038/s41396-019-0484-y
Takahashi, K., and Andoh, A. (2016). Reduced abundance of butyrate-producing bacteria species in the fecal microbial community in Crohn's disease. J. Gastroenterol. Hepatol. 93, 174–178. doi: 10.1159/000444541
Thomson, W. H. (1975). The nature and cause of haemorrhoids. Proc. R. Soc. Med. 68, 574–575. doi: 10.1177/003591577506800916
Yuan, Y., Lu, Y., Zhang, Z., Cheng, W., Yan, K., Zheng, Y., et al. (2023). Characteristics of the Cajal interstitial cells and intestinal microbiota in children with refractory constipation. Microb. Pathog. 184:106373. doi: 10.1016/j.micpath.2023.106373
Keywords: hemorrhoids, microbial community structure, bacterial infection, microbiomic, thrombosed hemorrhoid
Citation: Wang Y, Su W, Liu Z, Wang Y, Li L, Xu H, Wang M and Shen W (2024) The microbiomic signature of hemorrhoids and comparison with associated microbiomes. Front. Microbiol. 15:1329976. doi: 10.3389/fmicb.2024.1329976
Edited by:
Thi Thu Hao Van, RMIT University, AustraliaReviewed by:
Jinlyung Choi, PrognomiQ Inc., United StatesVincenzo Davide Palumbo, Euro-Mediterranean Institute of Science and Technology (IEMEST), Italy
Copyright © 2024 Wang, Su, Liu, Wang, Li, Xu, Wang and Shen. 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: Mingyu Wang, wangmingyu@sdu.edu.cn; Wenlong Shen, swl028073@qlyyqd.com
†These authors have contributed equally to this work