Abstract
Commensal microbial–host interaction is crucial for host metabolism, growth, development, and immunity. However, research on microbial–host immunity in large animal models has been limited. This study was conducted to investigate the effects of the commensal microbiota on immune function in two model groups: germ-free (GF) and specific-pathogen-free (SPF) piglets. The weight and organ index of the spleen of the GF piglet were larger than those in the SPF piglet (P < 0.05). The histological structure of the red pulp area and mean area of germinal centers were larger in the SPF piglet than in the GF piglet (P < 0.05), whereas the areas of staining of B cells and T cells in the spleen and mesenteric lymph nodes (MLNs) were lower in the GF piglet (P < 0.05). We identified immune-related genes in the spleen and MLNs using RNA sequencing, and used real-time quantitative PCR to analyze the expression of core genes identified in gene set enrichment analysis. The expression levels of genes in the transforming growth factor-β/SMAD3 signaling pathway, Toll-like receptor 2/MyD88/nuclear factor-κB signaling pathway, and pro-inflammatory factor genes IL-6 and TNF-α in the spleen and MLNs were higher in the SPF piglet and in splenic lymphocytes compared with those in the GF and control group, respectively, under treatment with acetic acid, propionic acid, butyric acid, lipopolysaccharide (LPS), or concanavalin A (ConA). The abundances of plasma cells, CD8++ T cells, follicular helper T cells, and resting natural killer cells in the spleen and MLNs were significantly greater in the SPF piglet than in the GF piglet (P < 0.05). In conclusion, the commensal microbiota influenced the immune tissue structure, abundances of immune cells, and expression of immune-related pathways, indicating the importance of the commensal microbiota for spleen and MLNs development and function. In our study, GF piglet was used as the research model, eliminating the interference of microbiota in the experiment, and providing a suitable and efficient large animal research model for exploring the mechanism of “microbial-host” interactions.
1 Introduction
The commensal microbiota colonizes various tissues of the body, including the gastrointestinal tract, oral cavity, respiratory tract, urogenital tract, and skin, providing immunomodulatory signals during the development, differentiation, and activation of immune cells, which are important for the maintenance of immune homeostasis (Lloyd and Marsland, 2017; ). Hosts have coevolved with their commensal microbiota in a complex mechanism that affects innate and adaptive immunity (; ). Adaptive immunity mainly comprises B cell-mediated humoral immunity and T cell-mediated cellular immunity (McCoy et al., 2017). The commensal microbiota directly regulates the differentiation and activation of B and T cells (), affecting the secretion of immunoglobulin (Ig)A, IgG, IgM, IgD, and IgE by plasma cells, forming the first barrier, and protecting the body from harm via infection (). As the first line of host defense, the innate immune system recognizes microbial pathogens with host germline-encoded pattern recognition receptors (PRRs) that bind unique pathogen-associated molecular patterns (PAMPs) (Takeuchi and Akira, 2010).
The spleen, as the largest secondary lymphoid tissue in the body of animals, is divided into different compartments by B and T lymphocytes: follicles containing B cells and periarteriolar lymph sheaths containing T cells (Steiniger et al., 2011; ). The most developed lymphoid regions in mesenteric lymph nodes (MLNs) comprise three main areas: follicular region, paracortex, and cortex (Zhang et al., 2012). Dendritic cells (DCs) in the spleen can reach various tissues of the body through the bloodstream, taking up antigens in tissues. After activation of Toll-like receptors (TLRs), cells involved in immune functions present antigens through white pulp to trigger initial activation of T cells, which can induce immune tolerance (Zarember and Godowski, 2002; Yi et al., 2015). Microbiome-derived TLR ligands and metabolites act directly on enterocytes and intestinal immune cells, but can also travel through the systemic circulation to modulate immunity in remote tissues, such as the spleen (; Shao et al., 2016). Tight junctions between the intestinal epithelial cells and mucus produced by goblet cells and hematopoietic cells, such as mononuclear cells and innate lymphoid cells (ILCs), including natural killer (NK) cells, are important effector cells of the innate immune system (; ). NK cells express T-bet and produce tumor necrosis factor (TNF)-α (). The commensal microbiota has been identified as an important regulator of immune cell activation and inflammation, acting through nuclear factor kappa B (NF-κB) and inflammasome activation (Xing et al., 2021). As a key enforcer of immune tolerance and a suppressor of inflammation, transforming growth factor-β (TGF-β) restricts multiple functions of the innate and adaptive immune systems, with Smad3 as an intracellular signaling protein in the TGF-β pathway (; Massagué and Sheppard, 2023). TGF-β not only enhances the B cell response (Reboldi et al., 2016) but also suppresses T cell proliferation and activation through differentiation of regulatory T cells (Tregs) (). Furthermore, specific intestinal microbes, possibly acting through metabolites such as short-chain fatty acids (SCFAs), contribute to immune homeostasis by modulating TGF-β production (Pang et al., 2018).
Germ-free (GF) mice have been used to study the effects of the commensal microbiota on the host immune system (; ; ). However, rodent models are limited by a number of important physiological and metabolic differences from humans (). In addition, the results of our recent experiments in Peyer’s patches (PPs) of GF piglets (Zhang et al., 2023) are not entirely consistent with those in PPs of GF mice (), suggesting that it may be difficult to replicate immune-related research conclusions based on rodent models in large animal models. GF and specific-pathogen-free (SPF) piglets are ideal large animal models for researching microbial–host interactions (Zhou et al., 2021). In this study, we explored the effects of the commensal microbiota on histological structure, gene expression, immune cell composition, and immune-related signal pathways in the host spleen and MLNs in GF and SPF piglet models. Our findings further elucidate the influence of the commensal microbiota on the development of the piglet immune system and support the use of the piglet immunological model as an alternative to widely used rodent models.
2 Materials and methods
2.1 GF and SPF piglet preparation
All animals were fed at the experimental base of Chongqing Academy of Animal Sciences. The preparation of Bama strain GF piglets and GF environment testing were based on Zhang et al. (2023). The GF piglets were reared in sterile isolators under GF conditions and were handfed Co60 γ-irradiated sterile 4.8%-fat cow milk powder diluted with sterile water (Anyou Group, Jiangsu, China). Once per week, the GF environments were checked for aerobic and anaerobic bacterial contamination (Sun et al., 2018). Of six GF piglets, three were assigned to the GF piglet; the other three were raised in rearing isolators in a barrier environment which the cleanliness could reach 10,000 and required independent flow of clean air with a high ventilation frequency (20 ∼ 60 times/h) (Ogden et al., 2016) and assigned to the SPF piglet. The pathogen detection report of SPF piglet was attached Supplementary Figure 1. The Ethics Committee of the Chongqing Academy of Animal Science reviewed the relevant ethical issues and approved the study (permit number XKY-No. 20210606).
2.2 Sample collection
At 42 days of age, all GF and SPF piglets were weighed, euthanized under isoflurane anesthesia, and exsanguinated in a sterile environment, followed by careful removal and weighing of the spleens. The mass of the MLNs in the jejunum was too slight to be weighed. A portion of the fresh spleen and MLNs were fixed with 4% (w/v) paraformaldehyde (Beyotime, China) for histomorphologic analysis. The remainder of the spleen and MLNs samples were snap-frozen in liquid nitrogen, and then stored at −80°C for further analysis.
2.3 Hematoxylin-eosin and immunohistochemistry staining
For hematoxylin-eosin (HE) staining, fixed tissues were dehydrated and embedded in paraffin. Tissue blocks were sliced, and the tissue sections were deparaffinized and stained with HE (Beyotime, China). B- and T-cell distributions were examined by immunohistochemistry (IHC) staining using CD3+ and CD20+ pig-specific makers, respectively. IHC staining was performed following the method described by Zhang et al. (2023). Semi-quantitative evaluation of IHC staining using an Image-Pro Plus (version 6.0) followed the method introduced in a previous report (Wang et al., 2009). Briefly, after B cell and T cell zones were identified based on the structure of the splenic central arteriole and the packing density of the cell populations (van Krieken et al., 1985), parameters including area and integrated optical density (IOD) were measured. The mean optical staining densities of the B-cell and T-cell zones were calculated as: sum of IOD/sum of area.
2.4 Total RNA extraction, library preparation, and sequencing
Total RNA was isolated from the frozen spleen and MLNs using TRIzol reagent (Takara Bio, Japan), in accordance with the manufacturer’s instructions. Total RNA samples with a ratio of absorbance at 260/280 nm ranging from 1.8 to 2.0 and an RNA integrity number >8 were used for library construction. We pooled the samples from each tissue by group and prepared a total of 12 sequencing libraries: three replicates each of GF versus SPF in spleen and MLNs. All libraries were constructed and sequenced on the DNBSEQ-T7 seq platform (Novogene, China), and 150-base pair paired-end reads were obtained. High-quality data (clean data) were generated from raw data through rigorous quality control, removing poly-N and low-quality reads, including those with ≥10% N, >10 nucleotides aligned to the adapter with ≤10% mismatches allowed, and >50% of bases with Phred quality <5. The RNA-seq data have been deposited in the GSA with the accession number CRA015219.
2.5 Expression analysis of mRNA
Levels of mRNA were quantified using Kallisto software (v 0.44.0), and the transcripts per kilobase of exon model per million mapped reads (TPM) value of each sample was calculated. An mRNA with a TPM value >1 and at least three repeats within one group was considered expressed. Differentially expressed (DE) genes were identified by edgeR (v 3.10) using cutoffs of P < 0.05 and | log2(Fold Change) | > 1(). Gene functional enrichment analysis and gene set enrichment analysis (GSEA) were performed at the Metascape1 (Zhou et al., 2019) and Omicshare2 websites, respectively. We uploaded each sample matrix to the CIBERSORTx database,3 to calculate the abundances of each immune cell type (Yao et al., 2023), selecting samples with a P value < 0.05.
2.6 Splenic lymphocyte isolation and culture
Spleens from 42-day-old piglets were aseptically removed and bubbled in pre-cooled phosphate-buffered saline (PBS) (Boster Bio, China) containing 5% penicillin/streptomycin for 3–5 min, followed by removal of the connective tissues around the spleen. After three washes in pre-cooled PBS, the spleen was cut into small pieces. Steps for isolating spleen lymphocytes using the Porcine Spleen Lymphocyte Isolation Kit (TBD Sciences, China) followed the protocol described by Ren et al. (2020). Splenic lymphocytes were suspended in RPMI-1640 complete medium and counted using trypan blue staining exclusion. Experiments were conducted on splenic lymphocytes reaching a concentration of 3.75 × 106 cells/ml with a viability >95%. Splenic lymphocytes were cultured in RPMI-1640 medium (Procell Life Science & Technology Co., Ltd., China) with 10% fetal bovine serum (Procell Life Science & Technology Co., Ltd., China) and 1% penicillin/streptomycin in a 37°C incubator with 5% CO2. The cells were divided into five treatment groups: control (treated with Phosphate Buffered Saline) (Servicebio, China), SCFAs (acetic acid, propionic acid, and butyric acid), lipopolysaccharide (LPS) (Beyotime, China), and concanavalin A (ConA) (Beyotime, China). SCFAs concentrations were based on Zhang et al. (2023), LPS and ConA concentrations were based on .
2.7 Cell proliferation assays
Cell proliferation was measured using a Cell Counting Kit (CCK)-8 Assay Kit (Beyotime, China) following the manufacturer’s instructions. Splenic lymphocytes were seeded in a 96-well plate (2 × 105 cells/well), divided into three replicates for each of the six treatment groups. After applying each treatment, cells were cultured and assayed at 1, 3, 6, and 12 h by the addition of 10 μl/well CCK-8. Optical density (OD) was measured at 450 nm using a microplate reader (Thermo Fisher Scientific, Spain).
2.8 Real-time quantitative PCR
Total RNA was extracted from frozen spleen, MLNs samples, and splenic lymphocytes using HiPure Total RNA Mini Kits (Magen, China), in accordance with the manufacturer’s instructions. Synthesis of cDNA was carried out as described by Pistol et al. (2015), using a PrimeScript RT reagent Kit with gDNA Eraser (Takara Bio, Japan) following the manufacturer’s instructions. Primers were designed through the National Center for Biotechnology Information database and synthesized commercially by Tsingke (China). Real-time quantitative PCR (RT-qPCR) was performed on an ABI Prism 7000 detection system in a two-step protocol with SYBR Green (Applied Bio-Systems, Foster City, CA, USA). Each reaction was performed in a volume containing 1 μl cDNA, 5 ml SYBR Premix Ex Taq (2×), 0.2 μl ROX reference dye (50×), 0.4 μl of each forward and reverse primer, and 3 μl PCR-grade water. The following PCR reaction was performed in a Quant Studio 6 Flex (Thermo Fisher Scientific, USA): 95°C for 30 s, followed by 40 cycles of 95°C for 5 s and 60°C for 34 s, and ending with a melting curve analysis (65–95°C). The relative expression of mRNA was calculated using the 2–ΔΔCT method and expressed as fold-change relative to the GF or cell culture control group. GADPH was used to normalize the gene expression. The specific primer sequences in this study are listed in Supplementary Table 1.
2.9 Statistical analysis
All statistical analyses were performed using Graphpad Prism 9.5. Unless otherwise stated, numerical data are presented as the mean ± standard deviation (SD). For single comparisons between two groups, the unpaired t-test was used. P-values < 0.05 and <0.01 were considered to indicate statistical significance.
3 Results
3.1 Effects of the commensal microbiota on histological structure in the piglet spleen and MLNs
The spleen and MLNs play important roles in innate and adaptive immunity. We established GF and SPF piglets to explore the histological differences in spleen and MLNs between the two models. The spleen weight and organ index were significantly lower in the GF piglet than in the SPF piglet (P < 0.05; Figures 1A, B). Although red pulp structures and lymphatic follicles in the spleen were obvious in both piglets (Figure 1C), the red pulp area was larger in the SPF piglet (Figure 1E). The mature lymphoid regions in MLNs were as follows: the follicular region, in which many germinal centers were located; the paracortex, a lymphatic region located between the cortex and the medulla; and the sinusoid region, comprising the lymphatic cord and sinus (Figure 1D). The mean area of the germinal center was larger in the SPF piglet than in the GF piglet (Figure 1F). These results indicated that the commensal microbiota promotes the development and maturation of the spleen and MLNs in piglets.
FIGURE 1
3.2 Effects of the commensal microbiota on the distribution of B cells and T cells in the piglet spleen and MLNs
The spleen and MLNs contain a variety of immune cells that co-participate in B cell-mediated humoral immunity and T cell-mediated cellular immunity (McCoy et al., 2017). IHC staining of porcine-specific B and T cell markers revealed that B and T cells were mainly distributed in the lymphatic follicles of the spleen and in the cortical part of the MLNs (Figures 2A, B). The proportions of T cells were slightly higher in the spleen of the SPF piglet than in the GF piglet (P = 0.0843) and were significantly higher in the MLNs of the SPF piglet than in those of the GF piglet (P < 0.05; Figure 2C). Furthermore, the proportions of B cells were significantly higher in the spleen and MLNs of the SPF piglet than in those of the GF piglet (P < 0.05; Figure 2D), indicating that B and T cell development in the spleen and MLNs are influenced by colonization of the commensal microbiota.
FIGURE 2
3.3 Effects of the commensal microbiota on mRNA expression in the piglet spleen and MLNs
RNA sequencing (RNA-seq) was used to explore the gene expression patterns of spleen and MLNs in the presence or absence of the commensal microbiota. We obtained ∼218.10 G clean data from 12 libraries (18.175 G per library) (Supplementary Table 2), with a total of 12,434 mRNAs identified. Principal component analysis (PCA) and hierarchical clustering showed different expression patterns between the spleens and MLNs of the GF and SPF piglets (Figures 3A, B). A total of 809 (600 upregulated and 209 downregulated) and 346 (120 upregulated and 226 downregulated) DE mRNAs were obtained in the spleen and MLNs, respectively (Figure 3C and Supplementary Table 3). We had obtained some highly upregulated genes, such as TNF receptor superfamily (TNFRSF) was found in spleen and MLNs upregulated genes, namely TNFRSF19 and TNFRSF17, which function to activate signaling pathways for development, organogenesis, cell death, and survival (Zhao et al., 2018). In addition, USP1 associated with innate immunity (Zhai et al., 2024), was also found in the upregulated genes, and IRF7 expressed in lymphocytes consistented with the above results (). Gene Ontology (GO) enrichment showed that many immune-related pathways, such as cellular response to external stimulus and defense response to virus were enriched in the spleen, while adaptive immune response and antiviral innate immune response were enriched in MLNs (Figures 3D, E and Supplementary Table 4). These results indicated that colonization of the commensal microbiota enhances the expression of immune-related signaling pathways and promotes the function of the immune system.
FIGURE 3
3.4 Effects of the commensal microbiota on innate and adaptive immune-related signaling pathways in the piglet spleen and MLNs
The colonization of the commensal microbiota significantly stimulated the immune status of the spleen and MLNs. Therefore, we used the GO database to download genes and products related to inflammatory factors (Supplementary Table 5). In the spleen and MLNs, the majority of genes in the gene set of the SPF piglets was significantly upregulated (Supplementary Figure 2). However, gene set analysis only focused on specific related pathway. Therefore, we used GSEA to analyze the overall gene expression data of the spleen and MLNs, focusing on the synergistic changes of the entire gene set, and exploring the impact of the commensal microbiota on the overall immune-related signaling pathways in the spleen and MLNs. Compared with the GF piglet, the antigen processing and presentation of exogenous antigen-related signaling pathways (GO:0048002, GO:0019884, and GO:0003823) were significantly activated in the spleen of the SPF piglet (Figure 4A). The SPF piglet also showed not only activation of the innate and adaptive signaling pathways (GO:0045088 and GO:0002821), but also positive regulation of the lymphocyte-mediated immunity signaling pathway (G0:0002708) in MLNs compared with the GF piglet (Figure 4B). In addition, RT-qPCR was used to verify the vital factors regulating inflammatory signaling pathways and inflammation-related factors triggered by the antigen response shown in the results of GSEA. We detected the expression of an immune-suppressive cytokine (TGF-β)-related signaling pathway and found that the expression levels of TGF-β and SMAD3 were higher in the spleen and MLNs of the SPF piglet compared with those in the GF piglet (Figures 4C, D). Additionally, the expression levels of TLR2, NF-κB, MyD88, and proinflammatory factor genes IL-6 and TNF-α were also higher in the spleen and MLNs of the SPF piglet compared with those in the GF piglet (Figures 4C, D). These results revealed that the commensal microbiota impacts the specific immune-related signaling pathways of the spleen and MLNs, enhancing the expression of immune response pathways, including the TGF-β signaling pathway, and inflammatory signaling pathways (TLR2/MyD88/NF-κB).
FIGURE 4
3.5 Effects of the commensal microbiota on immune cellular abundance in the spleen and MLNs
The commensal microbiota regulates the function and abundance of multiple immune cell types (Scott and Mann, 2020). CIBERSORTx is a suitable tool for the assessment of cellular abundance and cell type-specific gene expression patterns from bulk tissue transcriptome profiles (Steen et al., 2020). We identified 16 and 17 immune cell types in the spleen and MLNs, respectively, using CIBERSORTx (Figures 5A, B). The abundances of plasma cells, CD8+ T cells, follicular helper T cells, and resting NK cells were significantly increased in the spleen and MLNs of the SPF piglet compared with those in the GF piglet (P < 0.05; Figure 5C). By contrast, the abundance of memory B cells was lower in the spleen and MLNs of the SPF piglet than in the GF piglet (Figure 5D). In addition, compared with the spleen, the abundances of resting NK cells and monocytes in MLNs were higher and lower, respectively (Figures 5E, F). These results indicated that the commensal microbiota alters the abundance of immune cells by influencing the cell type-specific gene expression patterns of spleen and MLNs. In addition, the commensal microbiota had different effects on the abundances of resting NK cells and monocytes.
FIGURE 5
3.6 Immune-related pathways verified by splenic lymphocyte culture
To verify the immune-related pathways underlying the regulation of tissue immunity function, we established an in vitro model in splenic lymphocytes (Supplementary Figure 3). Cell proliferation was significantly increased under treatment with acetic acid, propionic acid, butyric acid, LPS, or ConA for 1 h compared with control (P < 0.01; Figure 6A). However, proliferation dropped thereafter, with no significant differences between the treatment and control groups 3, 6, and 12 h of treatment (Figure 6A). Next, we used RT-qPCR to analyze the expression of immune-related pathway genes at 1 h post-treatment, including TGF-β/SMAD3 and TLR2/MyD88/NF-κB pathway genes, as well as inflammatory factor genes IL-6 and TNF-α. The expression levels of TGF-β were significantly increased under all treatments compared with control (P < 0.05; Figure 6B), while the expression of SMAD3 showed an increasing trend under all treatments (Figure 6B). The expression levels of TLR2, MyD88, NF-κB, IL-6, and TNF-α were also significantly increased in all treatment groups compared with control (Figure 6B). These results revealed that the tendency for expression of specific immune-signal pathways in the spleen and MLNs in vivo were consistent with those in lymphocytes in vitro.
FIGURE 6
4 Discussion
Removal of the commensal microbiota leads to defects in the immune system and incomplete immune cell development (). At present, the impacts of the commensal microbiota on the immune system are mostly being studied in GF mice or pseudo-GF mice treated with antibiotics (). Nevertheless, the porcine immune system has several peculiarities that separate it from that of mice. For instance, pigs have a high frequency of peripheral CD4+ and CD8+ effector/memory T cells (). Furthermore, compared with mice, pigs are more anatomically and immunologically similar to humans (Wang and Donovan, 2015). GF and SPF piglets provide ideal animal models for investigating the impacts of the commensal microbiota on the immune function of the spleen and MLNs.
The spleen in piglets with colonization of the commensal microbiota exhibited significantly higher weight and organ index in this study. The same results were also confirmed in the spleen and MLNs of antibiotic-treated mice (; Thackray et al., 2018). Furthermore, the presence of the commensal microbiota promoted the development of red pulp and germinal centers in the spleen and MLNs, respectively. In the absence of the commensal microbiota, the same immune structure was present, but was underdeveloped. This finding suggested that the commensal microbiota, instead of impacting organogenesis, influenced immune system development (). Furthermore, the commensal microbiota increased the distribution areas of specific B cells and T cells in the spleen. These results were confirmed in piglet PPs (Zhang et al., 2023), indicating that the commensal microbiota is crucial for the development of immune tissues and the differentiation and maturation of B and T cells (Wang et al., 2019).
Many immune-related DE mRNAs in the spleen and MLNs between GF and SPF piglets were observed in this study, which is consistent with previous findings in mice (Wang et al., 2023b). Conventional pigs and GF pigs also exhibit many immune-related DE mRNAs in the spleen (Sun et al., 2018). Interestingly, our functional enrichment analysis showed that spleen-specific and MLNs-specific DE mRNAs were involved in response to virus and antiviral innate immune response. These findings are consistent with previous studies showing that the commensal microbiota influences resistance to external stimuli aspects of the jejunum, colon, spleen, and oral mucosa (Sun et al., 2018). Additionally, our team confirmed the importance of the commensal microbiota in the phenotypic characteristics and gene expression of PPs in piglets (Zhang et al., 2023).
We used the GO database of immunomarker genes as a reference for the GSEA to identify specific signaling pathways that influence the immune-related functions of the spleen and MLNs. This revealed that antigen processing and presentation of exogenous, innate immunity and adaptive immunity had far greater enrichment in the spleen and MLNs under colonization of the commensal microbiota. Compared with GF animals, the presence of the commensal microbiota significantly stimulated the host immune system, thereby achieving its purpose of self-protection (; ). Inflammatory signaling pathway expression and activation of inflammatory factors indicate that homeostasis is maintained through an active immunoregulatory process (). In the presence of the commensal microbiota, RT-qPCR showed increased expression of the TGF-β/SMAD3 signaling pathway in this study. The dominant function of TGF-β is to regulate peripheral immune homeostasis, particularly in the microbe-rich, antigen-rich environment of the gut. The SMAD pathway regulates the production of IgA by B cells, maintains the protective mucosal barrier, and promotes the balanced differentiation of CD4+ T cells into inflammatory T helper 17 cells and suppressive FOXP3-expressing Tregs (Malhotra and Kang, 2013). Under treatment with SCFAs, LPS, and ConA, the expression levels of TLR2, MyD88, and NF-κB increased in the spleen and MLNs of the SPF piglet, as well as in splenic lymphocytes. Resident enteric bacteria activated intestinal IL-10-producing B cells through the TLR2/MyD88/PI3K pathway to reduce colonic T cell activation and maintain mucosal homeostasis (Mishima et al., 2019). Consistently, TLRs significantly enhance the ability of DCs to present antigen through the TLR/MyD88/NF-κB signaling pathway (). The release of proinflammatory cytokines is one way that the body resists infection of tissues or cells by pathogenic bacteria. An appropriate amount of proinflammatory cytokines regulates the immune response to some extent, leading to resistance or clearing of pathogens, promoting the repair of damaged tissues, inducing tumor cell apoptosis, and so forth (Rostami et al., 2013). Consistent with other studies, we confirmed the impact of SCFAs on the metabolism of immune cells, including immunosuppressive and inflammatory cells (Wang et al., 2023a). These data in piglets support the notion that the colonization of the commensal microbiota increases the expression of proinflammatory factors that, together with timely expression of inflammatory suppressive factors that inhibit their overexpression, maintain the immune homeostasis of the body.
The effect of the commensal microbiota on the abundance of immune cells has been investigated in multiple organs in mice and humans (Papotto et al., 2021). CIBERSORTx is an emerging technology that identifies cell types based on relative subsets of known RNA transcripts without physical cell sorting (Wang et al., 2021), and has been validated in multiple tumor types (). Using CIBERSORTx, we were surprised to find that the commensal microbiota significantly affected the abundances of innate and adaptive immune cells in the spleen and MLNs, especially those of memory B cells, plasma cells, multiple T cell types, and resting NK cells. It also revealed fewer B cells in the ileum and blood in the presence of the commensal microbiota or viruses in SPF piglets compared with those in GF piglets (Sinkora et al., 2014; Potockova et al., 2015). Stimulation of commensal microbes also led to significantly greater abundance of plasma cells in the SPF piglet than in the GF piglet, consistent with the decrease in IgA-producing plasma cells in GF mice (Liu and Rhoads, 2013). Additionally, commensal microbes promoted the abundances of related T cell types, consistent with previously reported results (). Moreover, colonization of the commensal microbiota led to the appearance of effector CD4+CD8+ αβ T helper and CD2+CD8– γδ T cells (Sinkora et al., 2011), while colonization with Staphylococcus aureus was negatively associated with the abundance of memory B cells (). Although commensal microbes are known to affect the abundance of lymphocytes (Sinkora and Butler, 2016), the reasons for the significant differences in memory B cells require further exploration. Direct contact between MLNs and the commensal microbiota resulted in significantly greater abundance of NKs in the MLNs compared to that in the spleen of the SPF piglet. The abundance of NK cells in different ecological niches of the small intestine also varied (Potockova et al., 2015). Similarly, the abundance of monocytes, as peripheral immune cells, was significantly greater in the spleen than in the MLNs. These findings suggested that colonization of the commensal microbiota had varying degrees of impact on immune cells in different parts of the piglet body, as previously shown in mice (Van Praet et al., 2015).
CIBERSORTx algorithm was accurate and feasible to analyze the proportion of immune cells in tissues using RNA-seq data, which has been repeatedly mentioned in many journals (Newman et al., 2019; Steen et al., 2020; Wang et al., 2021). However, the use of bulk RNA-seq data for deconvolution analysis in this study had the following shortcomings. First, because the CIBERSORTx database features only 22 specific types of immune cells, we were unable to obtain the exact types and abundance of all immune cells. Second, the resolution of bulk RNA-seq was insufficient in that it only captures the average expression level of immune tissues. Third, spleen and MLNs exhibit tissue heterogeneities in terms of cell subpopulations and functions associated with physiological processes that RNA-seq cannot fully elucidate. In future studies, we aim to explore this issue further by combining single-cell RNA-seq and spatial transcriptome technologies to explore the effects of the commensal microbiota on the precise cellular composition of distinct tissue regions and spatial gene expression of the immune system.
5 Conclusion
The presence of the commensal microbiota promoted the development and maturation of the spleen and MLNs, impacting the abundance of various immune cell types and the proliferation and differentiation of lymphocytes. This work also explored microbial–host interactions and provided evidence to support the use of the piglet immunological model as an alternative to widely used rodent models.
Statements
Data availability statement
The original contributions presented in the study are publicly available. The RNA-seq data have been deposited in the Genome Sequence Archive (GSA) with the accession number CRA015219.
Ethics statement
The animal studies were approved by the Chongqing Academy of Animal Science. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.
Author contributions
YL: Writing – original draft, Writing – review & editing. JZ: Writing – original draft, Writing – review & editing. GY: Validation, Writing – review & editing. CT: Data curation, Writing – review & editing. XKL: Data curation, Writing – review & editing. LL: Formal analysis, Writing – review & editing. KL: Formal analysis, Writing – review & editing. JS: Validation, Writing – review & editing. YD: Validation, Writing – review & editing. XWL: Formal analysis, Writing – review & editing. ML: Formal analysis, Writing – review & editing. LG: Supervision, Writing – review & editing. JM: Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was funded by the National Natural Science Foundation of China (32072687 and 32302712), the National Center of Technology Innovation for Pigs (NCTIP-XD/B13), the Science Foundation of the Sichuan Province (2022YFQ0022), Sichuan International Science and Technology Innovation Cooperation/Hong Kong, Macao, and Taiwan Science and Technology Innovation Cooperation Project (2021YFH0033), Ministry of Agriculture and Rural Affairs (no. A120202), and the Major Science and Technology Projects of Tibet Autonomous Region (no. XZ202101ZD0005N).
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1398631/full#supplementary-material
References
1
AhmedS.TravisS. D.íaz-BahamondeF. V. D.PorterD. D. L.HenryS. N.MykinsJ.et al (2021). Early influences of microbiota on white matter development in germ-free piglets.Front. Cell Neurosci.15:807170. 10.3389/fncel.2021.807170
2
AlexanderK. L.TarganS. R.ElsonC. O. (2014). Microbiota activation and regulation of innate and adaptive immunity.Immunol. Rev.260206–220. 10.1111/imr.12180
3
AntoniniM.Lo ConteM.SoriniC.FalconeM. (2019). How the interplay between the commensal microbiota, gut barrier integrity, and mucosal immunity regulates brain autoimmunity.Front. Immunol.10:1937. 10.3389/fimmu.2019.01937
4
BayerF.AscherS.PontarolloG.ReinhardtC. (2019). Antibiotic treatment protocols and germ-free mouse models in vascular research.Front. Immunol.10:2174. 10.3389/fimmu.2019.02174
5
BergenW. G.WuG. (2009). Intestinal nitrogen recycling and utilization in health and disease.J. Nutr.139821–825. 10.3945/jn.109.104497
6
ChoJ.HanS. C.Ho HwangJ.SongJ. (2023). Characterization of immune development of fetal and early-life of minipigs.Int. Immunopharmacol.120:110310. 10.1016/j.intimp.2023.110310
7
ClarkeE. C.BradfuteS. B. (2020). The use of mice lacking type I or both type I and type II interferon responses in research on hemorrhagic fever viruses. Part 1: Potential effects on adaptive immunity and response to vaccination.Antiviral. Res.174:104703. 10.1016/j.antiviral.2019.104703
8
DurandA.Audemard-VergerA.GuichardV.MattiuzR.DelpouxA.HamonP.et al (2018). Profiling the lymphoid-resident T cell pool reveals modulation by age and microbiota.Nat. Commun.9:68. 10.1038/s41467-017-02458-4
9
EberlG.ColonnaM.Di SantoJ. P.McKenzieA. N. (2015). Innate lymphoid cells. Innate lymphoid cells: A new paradigm in immunology.Science348:aaa6566. 10.1126/science.aaa6566
10
EndersM.FrankenL.PhilippM. S.KesslerN.BaumgartA. K.EichlerM.et al (2020). Splenic red pulp macrophages cross-prime early effector CTL that provide rapid defense against viral infections.J. Immunol.20487–100. 10.4049/jimmunol.1900021
11
Ganal-VonarburgS. C.DuerrC. U. (2020). The interaction of intestinal microbiota and innate lymphoid cells in health and disease throughout life.Immunology15939–51. 10.1111/imm.13138
12
Garcia-FloresV.RomeroR.FurcronA. E.LevensonD.GalazJ.ZouC.et al (2020). Prenatal maternal stress causes preterm birth and affects neonatal adaptive immunity in mice.Front. Immunol.11:254. 10.3389/fimmu.2020.00254
13
GensollenT.IyerS. S.KasperD. L.BlumbergR. S. (2016). How colonization by microbiota in early life shapes the immune system.Science352539–544. 10.1126/science.aad9378
14
GeukingM. B.BurkhardR. (2020). Microbial modulation of intestinal T helper cell responses and implications for disease and therapy.Mucosal Immunol.13855–866. 10.1038/s41385-020-00335-w
15
GolaA.DorringtonM. G.SperanzaE.SalaC.ShihR. M.RadtkeA. J.et al (2021). Commensal-driven immune zonation of the liver promotes host defence.Nature589131–136. 10.1038/s41586-020-2977-2
16
GommermanJ. L.RojasO. L.FritzJ. H. (2014). Re-thinking the functions of IgA(+) plasma cells.Gut Microbes5652–662. 10.4161/19490976.2014.969977
17
GuW.MadridD. M. C.JoyceS.DriverJ. P. (2022). A single-cell analysis of thymopoiesis and thymic iNKT cell development in pigs.Cell Rep.40:111050. 10.1016/j.celrep.2022.111050
18
HeY.BenedictisJ. N. DeCaimentF.van BredaS. G. J. (2023). Analysis of cell-specific transcriptional responses in human colon tissue using CIBERSORTx.Sci. Rep.13:18281. 10.1038/s41598-023-45582-6
19
HooperL. V.MacphersonA. J. (2010). Immune adaptations that maintain homeostasis with the intestinal microbiota.Nat. Rev. Immunol.10159–169. 10.1038/nri2710
20
HuY.CongX.ChenL.QiJ.WuX.ZhouM.et al (2016). Synergy of TLR3 and 7 ligands significantly enhances function of DCs to present inactivated PRRSV antigen through TRIF/MyD88-NF-κB signaling pathway.Sci. Rep.6:23977. 10.1038/srep23977
21
HuangT.WeiB.VelazquezP.BornemanJ.BraunJ. (2005). Commensal microbiota alter the abundance and TCR responsiveness of splenic naïve CD4+ T lymphocytes.Clin. Immunol.117221–230. 10.1016/j.clim.2005.09.012
22
IvanovI. I.TuganbaevT.SkellyA. N.HondaK. (2022). T cell responses to the microbiota.Annu. Rev. Immunol.40559–587. 10.1146/annurev-immunol-101320-011829
23
JinC.LagoudasG. K.ZhaoC.BullmanS.BhutkarA.HuB.et al (2019). Commensal microbiota promote lung cancer development via γδ T cells.Cell176:998-1013.e16. 10.1016/j.cell.2018.12.040
24
KandasamyS.ChatthaK. S.VlasovaA. N.RajashekaraG.SaifL. J. (2014). Lactobacilli and Bifidobacteria enhance mucosal B cell responses and differentially modulate systemic antibody responses to an oral human rotavirus vaccine in a neonatal gnotobiotic pig disease model.Gut Microbes5639–651. 10.4161/19490976.2014.969972
25
KimM.KimC. H. (2017). Regulation of humoral immunity by gut microbial products.Gut Microbes8392–399. 10.1080/19490976.2017.1299311
26
LeesC. W.BarrettJ. C.ParkesM.SatsangiJ. (2011). New IBD genetics: Common pathways with other diseases.Gut601739–1753. 10.1136/gut.2009.199679
27
LiC.LeeA.GrigoryanL.ArunachalamP. S.ScottM. K. D.TrisalM.et al (2022). Mechanisms of innate and adaptive immunity to the Pfizer-BioNTech BNT162b2 vaccine.Nat. Immunol.23543–555. 10.1038/s41590-022-01163-9
28
LiM. O.SanjabiS.FlavellR. A. (2006). Transforming growth factor-beta controls development, homeostasis, and tolerance of T cells by regulatory T cell-dependent and -independent mechanisms.Immunity25455–471. 10.1016/j.immuni.2006.07.011
29
LiS.ZhuZ.YangF.CaoW.YangJ.MaC.et al (2021). Porcine epidemic diarrhea virus membrane protein interacted with IRF7 to inhibit type I IFN production during viral infection.J. Immunol.2062909–2923. 10.4049/jimmunol.2001186
30
LiY.ZhongS.JinZ.YeG.ZhangT.LiuZ.et al (2023). Peyer’s patch-involved gut microbiota facilitates anti-HBV immunity in mice.Virus Res.331:199129. 10.1016/j.virusres.2023.199129
31
LiuY.RhoadsJ. M. (2013). Communication between B-cells and microbiota for the maintenance of intestinal homeostasis.Antibodies2535–553. 10.3390/antib2040535
32
LloydC. M.MarslandB. J. (2017). Lung homeostasis: Influence of age, microbes, and the immune system.Immunity46549–561. 10.1016/j.immuni.2017.04.005
33
MalhotraN.KangJ. (2013). SMAD regulatory networks construct a balanced immune system.Immunology1391–10. 10.1111/imm.12076
34
MassaguéJ.SheppardD. (2023). TGF-β signaling in health and disease.Cell1864007–4037. 10.1016/j.cell.2023.07.036
35
McCoyK. D.RonchiF.GeukingM. B. (2017). Host-microbiota interactions and adaptive immunity.Immunol. Rev.27963–69. 10.1111/imr.12575
36
MishimaY.OkaA.LiuB.HerzogJ. W.EunC. S.FanT. J.et al (2019). Microbiota maintain colonic homeostasis by activating TLR2/MyD88/PI3K signaling in IL-10-producing regulatory B cells.J. Clin. Invest.1293702–3716. 10.1172/jci93820
37
NewmanA. M.SteenC. B.LiuC. L.GentlesA. J.ChaudhuriA. A.SchererF.et al (2019). Determining cell type abundance and expression from bulk tissues with digital cytometry.Nat. Biotechnol.37773–782. 10.1038/s41587-019-0114-2
38
OgdenB. E.Pang WilliamW.AguiT.LeeB. H. (2016). Laboratory animal laws, regulations, guidelines and standards in China Mainland, Japan, and Korea.Ilar J.57301–311. 10.1093/ilar/ilw018
39
PangX.TangY. J.RenX. H.ChenQ. M.TangY. L.LiangX. H. (2018). Microbiota, Epithelium, Inflammation, and TGF-β Signaling: An Intricate Interaction in Oncogenesis.Front. Microbiol.9:1353. 10.3389/fmicb.2018.01353
40
PapottoP. H.YilmazB.Silva-SantosB. (2021). Crosstalk between γδ T cells and the microbiota.Nat. Microbiol.61110–1117. 10.1038/s41564-021-00948-2
41
PistolG. C.BraicuC.MotiuM.GrasM. A.MarinD. E.StancuM.et al (2015). Zearalenone mycotoxin affects immune mediators, MAPK signalling molecules, nuclear receptors and genome-wide gene expression in pig spleen.PLoS One10:e0127503. 10.1371/journal.pone.0127503
42
PotockovaH.SinkorovaJ.KarovaK.SinkoraM. (2015). The distribution of lymphoid cells in the small intestine of germ-free and conventional piglets.Dev. Comp. Immunol.5199–107. 10.1016/j.dci.2015.02.014
43
ReboldiA.ArnonT. I.RoddaL. B.AtakilitA.SheppardD.CysterJ. G. (2016). IgA production requires B cell interaction with subepithelial dendritic cells in Peyer’s patches.Science352:aaf4822. 10.1126/science.aaf4822
44
RenZ.GuoC.HeH.ZuoZ.HuY.YuS.et al (2020). Effects of deoxynivalenol on mitochondrial dynamics and autophagy in pig spleen lymphocytes.Food Chem Toxicol.140:111357. 10.1016/j.fct.2020.111357
45
RostamiK.Rostami NejadM.AsadzadehH.ZaliM. (2013). PTU-199 proinflammatory cytokine (IL-8) in microscopic enteritis.Gut62A130–A130. 10.1136/gutjnl-2013-304907.289
46
ScottN. A.MannE. R. (2020). Regulation of mononuclear phagocyte function by the microbiota at mucosal sites.Immunology15926–38. 10.1111/imm.13155
47
ShaoL.FischerD. D.KandasamyS.SaifL. J.VlasovaA. N. (2016). Tissue-specific mRNA expression profiles of porcine Toll-like receptors at different ages in germ-free and conventional pigs.Vet. Immunol. Immunopathol.1717–16. 10.1016/j.vetimm.2016.01.008
48
SinkoraM.ButlerJ. E. (2016). Progress in the use of swine in developmental immunology of B and T lymphocytes.Dev. Comp. Immunol.581–17. 10.1016/j.dci.2015.12.003
49
SinkoraM.ButlerJ. E.LagerK. M.PotockovaH.SinkorovaJ. (2014). The comparative profile of lymphoid cells and the T and B cell spectratype of germ-free piglets infected with viruses SIV, PRRSV or PCV2.Vet. Res.45:91. 10.1186/s13567-014-0091-x
50
SinkoraM.StepanovaK.ButlerJ. E.FrancisD.Santiago-MateoK.PotockovaH.et al (2011). Ileal Peyer’s patches are not necessary for systemic B cell development and maintenance and do not contribute significantly to the overall B cell pool in swine.J. Immunol.1875150–5161. 10.4049/jimmunol.1101879
51
SteenC. B.LiuC. L.AlizadehA. A.NewmanA. M. (2020). Profiling cell type abundance and expression in bulk tissues with CIBERSORTx.Methods Mol. Biol.2117135–157. 10.1007/978-1-0716-0301-7_7
52
SteinigerB.BetteM.SchwarzbachH. (2011). The open microcirculation in human spleens: A three-dimensional approach.J. Histochem. Cytochem.59639–648. 10.1369/0022155411408315
53
SunJ.ZhongH.DuL.LiX.DingY.CaoH.et al (2018). Gene expression profiles of germ-free and conventional piglets from the same litter.Sci. Rep.8:10745. 10.1038/s41598-018-29093-3
54
TakeuchiO.AkiraS. (2010). Pattern recognition receptors and inflammation.Cell140805–820. 10.1016/j.cell.2010.01.022
55
ThackrayL. B.HandleyS. A.GormanM. J.PoddarS.BagadiaP.et al (2018). Oral antibiotic treatment of mice exacerbates the disease severity of multiple flavivirus infections.Cell Rep.22:3440–3453.e6. 10.1016/j.celrep.2018.03.001
56
van KriekenJ. H.Te VeldeJ.KleiverdaK.Leenheers-BinnendijkL.van de VeldeC. J. (1985). The human spleen; a histological study in splenectomy specimens embedded in methylmethacrylate.Histopathology9571–585. 10.1111/j.1365-2559.1985.tb02841.x
57
Van PraetJ. T.DonovanE.VanasscheI.DrennanM. B.WindelsF.DendoovenA.et al (2015). Commensal microbiota influence systemic autoimmune responses.Embo J.34466–474. 10.15252/embj.201489966
58
WangC.-J.ZhouZ.-G.HolmqvistA.ZhangH.LiY.AdellG.et al (2009). Survivin expression quantified by Image Pro-Plus compared with visual assessment.Appl. Immunohistochem. Mol. Morphol.17530–535. 10.1097/PAI.0b013e3181a13bf2
59
WangJ.HuangM.DuY.ChenH.LiZ.ZhaiT.et al (2023a). Lactobacillus rhamnosus GG regulates host IFN-I through the RIG-I signalling pathway to inhibit herpes simplex virus type 2 infection.Probiotics Antimicrob. Proteins12602–623. 10.1007/s12602-023-10137-8
60
WangJ.ZhuN.SuX.GaoY.YangR. (2023b). Gut-microbiota-derived metabolites maintain gut and systemic immune homeostasis.Cells12:793. 10.3390/cells12050793
61
WangL.ZhuL.QinS. (2019). Gut microbiota modulation on intestinal mucosal adaptive immunity.J. Immunol. Res.2019:4735040. 10.1155/2019/4735040
62
WangM.DonovanS. M. (2015). Human microbiota-associated swine: Current progress and future opportunities.Ilar J.5663–73. 10.1093/ilar/ilv006
63
WangY. H.HouH. A.LinC. C.KuoY. Y.YaoC. Y.HsuC. L.et al (2021). A CIBERSORTx-based immune cell scoring system could independently predict the prognosis of patients with myelodysplastic syndromes.Blood Adv.54535–4548. 10.1182/bloodadvances.2021005141
64
XingC.WangM.AjibadeA. A.TanP.FuC.ChenL.et al (2021). Microbiota regulate innate immune signaling and protective immunity against cancer.Cell Host Microbe29:959–974.e7. 10.1016/j.chom.2021.03.016
65
YaoS.HuangZ.WeiC.WangY.XiaoH.ChenS.et al (2023). CD79A work as a potential target for the prognosis of patients with OSCC: Analysis of immune cell infiltration in oral squamous cell carcinoma based on the CIBERSORTx deconvolution algorithm.BMC Oral Health23:411. 10.1186/s12903-023-02936-w
66
YiT.LiJ.ChenH.WuJ.AnJ.XuY.et al (2015). Splenic dendritic cells survey red blood cells for missing self-CD47 to trigger adaptive immune responses.Immunity43764–775. 10.1016/j.immuni.2015.08.021
67
ZaremberK. A.GodowskiP. J. (2002). Tissue expression of human toll-like receptors and differential regulation of Toll-like receptor mRNAs in leukocytes in response to microbes, their products, and cytokines.J. Immunol.168554–561. 10.4049/jimmunol.168.2.554
68
ZhaiY.DuY.YuanH.FanS.ChenX.WangJ.et al (2024). Ubiquitin-specific proteinase 1 stabilizes PRRSV nonstructural protein Nsp1β to promote viral replication by regulating K48 ubiquitination.J. Virol.98:e0168623. 10.1128/jvi.01686-23
69
ZhangJ.ShenY.YangG.SunJ.TangC.LiangH.et al (2023). Commensal microbiota modulates phenotypic characteristics and gene expression in piglet Peyer’s patches.Front. Physiol.14:1084332. 10.3389/fphys.2023.1084332
70
ZhangW.NasuT.HosakaY. Z.YasudaM. (2012). Comparative studies on the distribution and population of immunocompetent cells in bovine hemal node, lymph node and spleen.J Vet Med Sci.74405–411. 10.1292/jvms.11-0405
71
ZhaoX.WangY.GuoJ.WangJ. (2018). Correlation analyses of CpG island methylation of cluster of differentiation 4 protein with gene expression and T lymphocyte subpopulation traits.Asian Aust. J. Anim. Sci.311141–1149. 10.5713/ajas.17.0805
72
ZhouH.SunJ.YuB.LiuZ.ChenH.HeJ.et al (2021). Gut microbiota absence and transplantation affect growth and intestinal functions: An investigation in a germ-free pig model.Anim. Nutr.7295–304. 10.1016/j.aninu.2020.11.012
73
ZhouY.ZhouB.PacheL.ChangM.KhodabakhshiA. H.TanaseichukO.et al (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets.Nat. Commun.10:1523. 10.1038/s41467-019-09234-6
Summary
Keywords
commensal microbiota, immunity, germ-free piglet, spleen, mesenteric lymph node
Citation
Liu Y, Zhang J, Yang G, Tang C, Li X, Lu L, Long K, Sun J, Ding Y, Li X, Li M, Ge L and Ma J (2024) Effects of the commensal microbiota on spleen and mesenteric lymph node immune function: investigation in a germ-free piglet model. Front. Microbiol. 15:1398631. doi: 10.3389/fmicb.2024.1398631
Received
10 March 2024
Accepted
27 May 2024
Published
12 June 2024
Volume
15 - 2024
Edited by
Silvia Turroni, University of Bologna, Italy
Reviewed by
David Wylensek, Institute of Medical Microbiology, Germany
Krzysztof Flisikowski, Technical University of Munich, Germany
Updates
Copyright
© 2024 Liu, Zhang, Yang, Tang, Li, Lu, Long, Sun, Ding, Li, Li, Ge and Ma.
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: Liangpeng Ge, geliangpeng1982@163.comJideng Ma, jideng.ma@sicau.edu.cn
†These authors share first authorship
Disclaimer
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.