ORIGINAL RESEARCH article

Front. Nutr., 15 July 2022

Sec. Nutritional Immunology

Volume 9 - 2022 | https://doi.org/10.3389/fnut.2022.944390

Research on the Mechanism of HRP Relieving IPEC-J2 Cells Immunological Stress Based on Transcriptome Sequencing Analysis

  • 1. College of Animal Science and Veterinary Medicine, Heilongjiang Bayi Agricultural University, Daqing, China

  • 2. Shanxi Animal Husbandry and Veterinary School, Taiyuan, China

  • 3. College of Food Science, Heilongjiang Bayi Agricultural University, Daqing, China

  • 4. College of Animal Science and Technology, Jilin Agricultural University, Changchun, China

Article metrics

View details

29

Citations

3,1k

Views

1,4k

Downloads

Abstract

Early weaning increased the economic benefits of piglets. However, early weaning damages the intestinal barrier of piglets and causes immunological stress. The mechanism by which Hippophae rhamnoides polysaccharide (HRP) alleviates lipopolysaccharide (LPS)-induced intestinal porcine epithelial cells (IPEC-J2) inflammatory damage was investigated using proteomics in our previous studies. In this study we employed RNA-sequencing (RNA-seq) to determine the level and function of differentially expressed genes (DEGs) and further explore the mechanism of the HRP anti-inflammatory and immune process. The differential expression analysis indicated that 3622, 1216, and 2100 DEGs in the IPEC-J2 cells were identified in C vs. L, L vs. H6-L, and C vs. H6-L, respectively. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis foundsix identified pathways related to the immune system. Additionally, we used the Science, Technology, Engineering, and Math (STEM) program to categorize the 3,134 DEGs that were differentially expressed in H2-L, H4-L and H6-L into eight possible expression profiles, in which 612 were clustered into two profiles. The accuracy and consistency of RNA-seq data were validated by the results of qRT-PCR of the nuclear factor of kappa light polypeptide gene enhancer in B-cells 2 (NFKB2), MAP kinase interacting serine/threonine kinase 2 (MKNK2), mitogen-activated protein kinase kinase 1 (MAP2K1), mitogen-activated protein kinase kinase kinase 8 (MAP3K8), Ras-related protein R-Ras (RRAS), TNF receptor-associated factor 1 (TRAF1), NF-kappa-B inhibitor alpha (NFKBIA), interleukin 8 (IL8), tumor necrosis factor, alpha-induced protein 3 (TNFAIP3), and transforming growth factor beta-1 (TGFB1). Transcriptome sequencing also indicated that HRP reduced the expression levels of related DEGs and inhibited the activation of the mitogen-activated protein kinase (MAPK)/nuclear factor kappa-B (NF-κB) signaling pathway. Our findings indicate that the application of HRP in piglet diets during the early weaning period can improve intestinal epithelial function and integrity, and relieve intestinal damage, and improve piglet health.

Introdution

With the acceleration of the large-scale and intensified process of pig production, early weaning techniques for piglets have been gradually implemented (1). Weaned piglets are affected by stresses such as nutrition, immunity, and environment, and piglets have early weaning syndromes such as reduced feed utilization, poor growth, and diarrhea (2, 3). Piglet weaning is accompanied by the occurrence of intestinal inflammation, which causes a series of negative reactions in the incompletely developed intestinal tracts of piglets, such as intestinal mucosal injury, intestinal villi damage, and intestinal wall injury, Intestinal inflammation also affects the digestive and absorptive function of the intestinal tract in piglets and leads to sluggish growth, diarrhea, and even death, causing great economic losses to the swine industry (4). The intestinal mucosa is the host's first line of defense against pathogenic microorganisms. Intestinal epithelial cells (IECs) are an important part of the intestinal mucosal barrier (5, 6). Injury to IECs is an important pathological basis for intestinal dysfunction. IECs can produce severe immune stress in a variety of physiological, pathological, dietary, or environmental conditions, and the accumulation of excessive inflammatory cytokines can damage intestinal epithelial cells, resulting in intestinal dysfunction. Therefore, the key to reducing the diarrhea rate and improving the production performance of weaned piglets is to reduce weaning stress and protect the intestinal structure and function of piglets (7). Seeking immunostimulants to promote growth and reduce the intestinal inflammation caused by weaning has become a hot scientific issue.

In recent years, plant extracts as immunomodulators have received worldwide attention due to their nutritional and medicinal potential (8, 9). Hippophae rhamnoides L. (sea buckthorn) is a traditional medicinal plant (10). H.rhamnoides extracts have shown antioxidant, anti-inflammatory, and anti-viral effects by decreasing cytotoxicity and reactive oxygen species (ROS) generation (11, 12). A recent study showed that Hippophae rhamnoides polysaccharide (HRP) protected mice livers from CCl4-induced damage through its anti-inflammatory effect and the modulation of the balance between anti-inflammatory cytokines inimmune cells (13). The protective mechanism of chitosan oligosaccharide against LPS-induced inflammatory responses in IPEC-J2 and in mice with DSS dextran sulfate sodium-induced colitis is reported (14). In our previous study, the role of several key regulatory genes and proteins involved in HRP immunoregulation was examined (15, 16). However, the composition and mechanisms of the underlying global regulatory networks at the transcriptome levels are still poorly understood. To confirm the pre-protective effect of HRP on lipopolysaccharide (LPS)-induced intestinal porcine epithelial cells (IPEC-J2) in terms of anti-inflammatory or immunoregulatoryproperties, the pathways enriched by differentially expressed genes (DEGs) should be biologically validated. We employed RNA-sequencing (RNA-seq) to determine the abundance and function of genes (17, 18) and lay the foundation for further study of the anti-inflammatory immune response of HRP in order to provide a theoretical basis and technical support for the treatment and prevention of diseases such as diarrhea caused by intestinal damage in early-weaned piglets.

Materials and Methods

Materials

HRP, ≥98% (HPLC) purchased from Nanjing Zelang Biological Technology Co., Ltd., China. HRP was extracted by water decoction and alcohol precipitation as described previously (19, 20). The protein in the filtrate was removed primarily by the Sevage method, HRP was purified on a Sephadex G150 gel (Pharmacia) and the purity was also verified by HPLC. HRP consists of 70% carbohydrate and 14.2% uronic acid. The polysaccharide is composed of mannose, arabinose, glucose, galactose and rhamnose with a ratio of 2.02:1.02:4.24:1:9.22 as the indication of chromatographic analysis of HRP (19).

Cell Culture and Treatment

The IPEC-J2 cell line source and cells culture method was as described in our previous study (21). HRP was pretreated for 24 h. Subsequently, cells were exposed to LPS for 16 h. The cells were collected, and quickly frozen in liquid nitrogen, and stored at −80°C for future transcriptome analysis. This study included five treatments C represents the control group IPEC-J2 cells without treatment. L represents the treatment group of IPEC-J2 cells induced by 10 μg/mL LPS. H2-L, H4-L, and H6-L represent the pre-treatment of IPEC-J2 cells with 200, 400, and 600 μg/mL HRP, respectively, followed by treatment with 10 μg/mL LPS. The HRP and LPS concentrations used were based on the results of our previous articles (16, 17).

Total RNA Extraction, cDNA Library Construction, and RNA-seq

According to the manufacturer's instructions, TRIzol reagent was used to extract RNA from the IPEC-J2 cells in different groups. cDNA library construction and RNA-seq were as described by Shao et al. (17).

Bioinformatics Analysis

High quality clean reads were obtained by removing reads containing adapters, more than 10% unknown nucleotides (N), and low-quality reads containing more than 50% low quality (Q-value ≤ 20) bases. The short reads alignment tool Bowtie2 (22) was used to map reads to the ribosome RNA (rRNA) database. The rRNA mapped reads were removed. The remaining reads were further used in the assembly and analysis of the transcriptome. Gene abundances were quantified by the RSEM software (23). Transcript abundance was normalized to fragments per kilobase of exon model per million mapped reads (FPKM). DEGs between groups were analyzed using the DEGs R package (http://www.rproject.org/). Genes with fold change (FC) >1.2 and a false discovery rate (FDR) <0.02 were considered significant. Significant enrichment of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways determined the most important biochemical metabolic pathways and signal transduction pathways that the protein participated in. The KEGG website (http://www.kegg.jp/kegg/) was used to query the immune-related signal pathways involved in significantly differently changed genes, and heat map analysis was performed on the significantly differentially expressed genes involved in immune-related signal pathways. For each treatment group, it is included three replicates to increase reliability of the data.

Verification of Transcriptome Data Using qRT-PCR

To validate the accuracy of the RNA-seq data, the DEG expressions for the nuclear factor of kappa light polypeptide gene enhancer in B-cells 2 (NFKB2), MAP kinase interacting serine/threonine kinase 2 (MKNK2), mitogen-activated protein kinase kinase 1 (MAP2K1), mitogen-activated protein kinase kinase kinase 8 (MAP3K8), Ras-related protein R-Ras (RRAS), TNF receptor-associated factor 1 (TRAF1), NF-kappa-B inhibitor alpha (NFKBIA), interleukin 8 (IL8), tumor necrosis factor, alpha-induced protein 3 (TNFAIP3), and transforming growth factor beta-1 (TGFB1) were determined by qRT-PCR analyses as described previously (16). We focused on the 600 μg/mL concentration of H. rhamnoides polysaccharide to pre-treat the IPEC-J2 cells in the following results (16, 21). The primers for the selected genes are listed in Supplementary Table S1.

Statistical Analysis

For the qPCR analysis, all data analyses were performed by using SPSS software (SPSS version 20.0, Chicago, IL, USA). Tukey's multiple range test was used to compare the mean values (P < 0.05) to indicate significant differences and the results were expressed as mean ± SD.

Results

Transcriptome Profiles

The aggregate of 356,912,016 raw reads were accumulated by using the High-throughput RNA sequencing to do the paired-end sequencing of the fifteen constructed libraries. After quality control, assessment of contaminated rRNA and low-quality sequences, 356,207,308 qualified Illumina reads were obtained. Approximately 97.44% of the clean reads were mapped to the reference genome and then used for further gene expression analysis (Table 1). We focused on C, L, and H6-L samples in the following experiment.

Table 1

SampleGroupRaw readsClean readsBF_Q30 (%)AF_Q30
(%)
BF_GC
(%)
AF_GC
(%)
Total_MappedMapping_
rate (%)
Unique_MappedMultiple_Mapped
C-1C2,08,95,1362,08,62,49893.8793.9957.0657.062,03,31,96197.781,97,11,0716,20,890
C-22,10,28,5802,09,87,07293.4393.5455.1455.132,03,48,75197.451,97,28,8266,19,925
C-32,25,18,5322,247,268893.4893.6155.3255.322,17,40,75397.472,10,55,1866,85,567
L-1L2,03,00,2922026257293.6593.7956.1756.161,95,66,82297.521,88,75,3436,91,479
L-22,37,73,9662,37,31,95493.6893.8156.0956.092,30,56,06497.612,22,76,1957,79,869
L-32,36,25,6402,35,82,59893.5193.6355.9755.962,29,31,15397.662,21,51,5007,79,653
H2-L-1H2-L2,11,43,1642,11,06,06693.994.0256.4956.492,04,55,40097.631,97,62,8466,92,554
H2-L-22,24,45,9062,23,96,56493.1993.3256.1456.142162994897.262,09,15,0537,14,895
H2-L-32,46,04,6182,45,58,73693.6693.7856.0456.032,37,74,06697.562,29,82,4877,91,579
H4-L-12,87,80,5402,87,16,00893.1393.2756.3856.382,77,37,54097.332,68,17,0899,20,451
H4-L-2H4-L2,45,62,6642,45,06,63092.9593.156.6256.632,37,01,50497.272,29,11,1807,90,324
H4-L-32,77,16,8202,76,60,21093.4993.6155.6855.682,67,26,02697.322,58,29,8668,96,160
H6-L-12,68,24,3102,67,75,97493.4493.5656.1356.132,58,68,45997.42,50,02,3068,66,153
H6-L-2H6-L2,55,24,1642,54,70,01093.4593.656.2756.272,45,00,61497.062,36,43,2988,57,316
H6-L-32,31,67,6842,31,17,72893.3793.555.9355.922,23,36,84597.252,15,99,2897,37,556

Summary statistics for sequence quality and alignment information of IPEC-J2 cells sample in every group.

C1–C3 represent the control-group IPEC-J2 cells without treatment, L1-L3 represent the replicates of treatment IPEC-J2 cells induced by LPS with 10 μg/mL, H2-L-1- H2-L-3 represent the replicates of pre-treatment IPEC-J2 cells with 200 μg/mL HRP and followed by co-treatment with 10 μg/mL LPS, H4-L-1- H4-L-3 represent the replicates of pre-treatment IPEC-J2 cells with 400 μg/mL HRP and followed by co-treatment with 10 μg/mL LPS, and H6-L-1- H6-L-3 represent the replicates of pre-treatment IPEC-J2 cells with 600 μg/mL HRP and followed by co-treatment with 10 μg/mL LPS.

DEGs were identified using digital gene expression tags. The expression analysis indicated that 3,622, 1,216, and 2,100 DEGs in the IPEC-J2 cells were identified in C vs. L, L vs. H6-L, and C vs. H6-L, respectively (FC >1.2 and FDR <0.02) The most DEGs were identified from C to L (Figures 1A–D and Supplementary Tables S2S4). These results demonstrated that number of inflammatory cytokines were accentuated after LPS induction. However, number of inflammatory cytokines in HRP pre-treated groups were significantly reduced. The volcano plot shows that HRP plays an important role in the immune regulation of cellular inflammatory damage.

Figure 1

Functional Enrichment of DEGs

DEGs were categorized into three Gene Ontology (GO) groups: biological process, cellular component, and molecular function. In the biological process group, many DEGs were categorized as the cellular process, single-organism process, and metabolic process. In the cellular component group, many DEGs were categorized as the cell, cell part and organelle. In the molecular function group, many DEGs were categorized as the binding, catalytic activity, and transporter activity (Figures 2A–C and Supplementary Tables S5S7).

Figure 2

As shown in Figures 3A–C and Supplementary Tables S8S10, in C vs. L, 1,442 DEGs were mapped into 322 KEGG pathways. The key pathways were Apoptosis (ko04210) and the NF-kappa B signaling pathway (ko04064). In L vs. H6-L, 531 DEGs were mapped into 307 KEGG pathways. The key pathways were MAPK (ko04010), NF-kappa B (ko04064), TNF (ko04668), Apoptosis (ko04210), HIF-1 (ko04066), and the TGF-beta signaling pathway (ko04350). In C vs. H6-L, 892 DEGs were mapped into 318 KEGG pathways. The key pathway was the IL-17 signaling pathway (ko04657).

Figure 3

Heat Map of the DEGs Related to Immune Pathways in L vs. H6-L

In Figure 4 and Table 2, the sample gene expression of the C and L, L and H6-L, C and H6-L groups is shown using a pseudo color scale, with high levels of expression shown in red and low levels shown in blue. The results showed that most of the gene expression levels were down-regulated after HRP pre-treatment compared with the LPS group. The Log2-fold change and p-value statistics of DEGs related to the main immune pathways in the L-vs.-H6-L group are shown in Table 3. This data indicated that genes were significantly expressed in six immune pathways.

Figure 4

Table 2

Gene IDGene
name
Description
LOC100153829NFKB2Nuclear factor of kappa light polypeptide gene enhancer in B-cells 2
LOC100517077MKNK2MAP kinase interacting serine/threonine kinase 2
LOC100233191MAP2K1Mitogen-activated protein kinase kinase 1
LOC396648HSPA1sHeat shock 70kda protein 1/2/6/8
LOC100158082PPP3CSerine/threonine-protein phosphatase 2B catalytic subunit
LOC100152997GADD45Growth arrest and dna-damage-inducible protein
LOC100622217MAP3K8Mitogen-activated protein kinase kinase kinase 8
LOC100240743DDIT3DNA damage-inducible transcript 3
LOC493184HSPB1Heat shock protein beta-1
LOC100513084CSF1Macrophage colony-stimulating factor 1
LOC100516558RRASRas-related protein R-Ras
LOC110259297RRAS2Ras-related protein R-Ras2
LOC396985PLAUurokinase plasminogen activator
LOC100627029TRAF1TNF receptor-associated factor 1
LOC406188NFKBIANF-kappa-B inhibitor alpha
LOC100125555LY96Lymphocyte antigen 96
LOC100623776TRIFToll-like receptor adapter molecule 1
LOC396750ICAM1Intercellular adhesion molecule 1
LOC100153830MALT1Mucosa-associated lymphoid tissue lymphoma translocation protein 1
LOC396880IL8Interleukin 8
LOC100518917PARPPoly
LOC100622156TNFAIP3Tumor necrosis factor, alpha-induced protein 3
LOC100622590BIRC2_3Baculoviral IAP repeat-containing protein 2/3
LOC100515635CREB5Cyclic amp-responsive element-binding protein 5
LOC396915EDN1Endothelin-1
LOC397454ITPR1Inositol 1,4,5-triphosphate receptor type 1
LOC397278PMAIP1Phorbol-12-myristate-13-acetate-induced protein 1
LOC407245LDHL-lactate dehydrogenase
LOC396823GAPDHGlyceraldehyde 3-phosphate dehydrogenase
LOC445512HMOX1Heme oxygenase 1
LOC407608PGKPhosphoglycerate kinase
LOC100512013ALDOFructose-bisphosphate aldolase, class I
LOC100126864IFNGR2Interferon gamma receptor 2
LOC397062TFRCTransferrin receptor
LOC396862TIMP1Metallopeptidase inhibitor 1
LOC100623582NBL1Neuroblastoma suppressor of tumorigenicity 1
LOC397490INHBBInhibin beta B chain
LOC100518554INHBEInhibin beta e chain
MSTRG.7730PDGFBPlatelet-derived growth factor subunit b
LOC100517416DUSP10Dual specificity phosphatase 10
LOC397260SMAD2_3Mothers against decapentaplegic homolog 2/3
LOC100152069SMAD6Mothers against decapentaplegic homolog 6
LOC100522839ID1DNA-binding protein inhibitor ID1
LOC100522677E2F4_5Transcription factor e2f4/5
LOC396691BMPR1BBone morphogenetic protein receptor type-1b
LOC654298ID2DNA-binding protein inhibitor ID2
LOC100626978ID3DNA-binding protein inhibitor ID3
LOC100144508ID4DNA-binding protein inhibitor ID4
LOC100521305SMAD7Mothers against decapentaplegic homolog 7
LOC100516289PFKFB36-phosphofructo-2-kinase / fructose-2,6-biphosphatase 3
LOC100156226EP300E1A/CREB-binding protein
LOC100192445PDPK13-phosphoinositide dependent protein kinase-1
LOC397027CAPN1Calpain-1
LOC493186SOCS3Suppressor of cytokine signaling 3
LOC100621398MRASRas-related protein M-Ras
LOC100620969EREGEpiregulin
LOC397078TGFB1Transforming growth factor beta-1
LOC100524751RIPK1Receptor-interacting serine/threonine-protein kinase 1
LOC100157144DUSPDual specificity MAP kinase phosphatase
MSTRG.7893NR4A1Nuclear receptor subfamily 4 group A member 1
LOC100524021CACNB3Voltage-dependent calcium channel beta-3
LOC102166564PDGFRBPlatelet-derived growth factor receptor beta
LOC100514115FGFR3Fibroblast growth factor receptor 3
LOC654328METProto-oncogene tyrosine-protein kinase Met

Significant differentially expressed genes related to the main immune pathways in the L-vs-H6-L comparison group.

Table 3

Gene IDGene
name
Log2 fold changeP-value
C-vs-LL-vs-H6-LC-vs-LL-vs-H6-L
LOC100153829NFKB21.16−0.439.63E-392.65E-18
LOC100517077MKNK21.16−0.355.10E-652.88E-39
LOC100233191MAP2K10.61−0.452.23E-123.99E-12
LOC396648HSPA1s0.82−0.592.61E-197.87E-114
LOC100158082PPP3C0.48−0.295.85E-070.00045
LOC100152997GADD452.28−2.189.91E-068.74E-06
LOC100622217MAP3K81.77−0.761.82E-225.36E-09
LOC100240743DDIT31.61−0.762.11E-238.35E-09
LOC4 93184HSPB11.98−0.654.03E-801.82E-68
LOC100513084CSF12.07−1.018.41E-281.90E-11
LOC100516558RRAS1.40−0.424.76E-525.56E-15
LOC110259297RRAS20.61−0.294.70E-123.17E-13
LOC396985PLAU1.28−0.571.14E-581.10E-70
LOC100627029TRAF12.28−0.882.27E-110.00050
LOC406188NFKBIA0.97−0.774.03E-215.22E-21
LOC100125555LY961.79−0.842.07E-161.29E-07
LOC100623776TRIF0.02−0.480.284300.00119
LOC396750ICAM10.82−0.339.75E-181.23E-11
LOC100153830MALT11.49−0.781.83E-171.33E-08
LOC396880IL81.71−0.984.04E-372.07E-18
LOC100518917PARP0.40−0.350.000130.00022
LOC100622156TNFAIP31.58−0.994.94E-374.43E-28
LOC100622590BIRC2_30.50−0.610.001081.80E-05
LOC100515635CREB52.53−1.275.95E-174.61E-07
LOC396915EDN10.11−0.680.029096.20E-11
LOC397454ITPR10.81−0.391.26E-192.07E-08
LOC397278PMAIP10.90−0.691.34E-091.72E-07
LOC407245LDH0.24−0.321.53E-073.24E-28
LOC396823GAPDH0.71−0.283.78E-198.99E-41
LOC445512HMOX12.38−1.376.68E-532.62E-28
LOC407608PGK0.59−0.444.28E-234.92E-58
LOC100512013ALDO0.98−0.603.18E-172.34E-12
LOC100126864IFNGR21.39−0.509.80E-315.24E-10
LOC397062TFRC0.03−0.610.077178.04E-21
LOC396862TIMP10.81−0.281.02E-164.82E-06
LOC100623582NBL11.29−0.341.58E-220.00014
LOC397490INHBB1.95−0.997.16E-106.23E-05
LOC100518554INHBE0.80−0.374.79E-152.46E-08
MSTRG.7730PDGFB1.40−0.401.78E-261.14E-05
LOC100517416DUSP102.11−1.005.37E-251.01E-10
LOC397260SMAD2_3−0.460.395.08E-063.01E-17
LOC100152069SMAD6−1.570.991.65E-136.05E-07
LOC100522839ID1−1.891.642.86E-633.67E-217
LOC100522677E2F4_5−0.900.387.12E-112.86E-06
LOC396691BMPR1B−1.080.412.31E-142.66E-05
LOC654298ID2−1.341.033.63E-132.60E-12
LOC100626978ID3−1.070.812.45E-195.13E-26
LOC100144508ID4−1.750.668.12E-180.00041
LOC100521305SMAD7−0.970.637.04E-079.90E-05
LOC100516289PFKFB3−0.800.683.84E-062.42E-11
LOC100156226EP300−1.150.523.85E-193.15E-08
LOC100192445PDPK1−0.730.313.17E-070.00189
LOC397027CAPN1−0.660.362.47E-061.36E-13
LOC493186SOCS3−1.110.940.000581.46E-06
LOC100621398MRAS2.52−1.472.76E-060.00043
LOC100620969EREG−0.761.812.39E-060.00046
LOC397078TGFB1−0.560.432.25E-055.57E-13
LOC100524751RIPK1−0.810.351.44E-080.00057
LOC100157144DUSP−0.260.350.411876.08E-05
MSTRG.7893NR4A1−0.290.850.392077.41-17
LOC100524021CACNB3−0.270.360.393430.00266
LOC102166564PDGFRB−1.201.130.000168.65E-06
LOC100514115FGFR3−1.200.578.41E-213.70E-35
LOC654328MET−1.040.485.19E-232.23E-17

Log2-fold change and p-value statistics of significant differentially expressed genes related to the main immune pathways in the L-vs-H6-L group.

Cluster Analysis of DEGs Among the Three HRP Concentrations

To determine the gene expression trajectories, we used the STEM program to categorize the 3,134 DEGs that were differentially expressed in H2-L, H4-L, and H6-L into eight possible expression profiles (P < 0.05) (Figures 5A,B, 6A–H, and Supplementary Table S11), in which 612 were clustered into two profiles (P ≤ 0.05), including two up-regulated patterns (Profile 6 and Profile 7). Profile 6 and 7 contained 452 and 160 DEGs, respectively. The consistent up-regulation of genes of Profile 6 in H2-L, H4-L, and H6-L indicated that DEGs may contribute to stimulatory functions during the polysaccharide anti-inflammatory function. The up-regulated genes of Profile 7 between only H2-L and H4-L revealed that these DEGs played a key role in the anti-inflammatory process. There was no significant difference in the up-regulated genes between H4-L and H6-L.

Figure 5

Figure 6

KEGG Pathway Enrichment Analysis of Differentially Expressed Genes Among the Three HRP Concentrations

A total of 11.1% (940/8,499) of the DEGs could be annotated. As shown in Table 4, the metabolic pathways (ko01100), Cytokine-cytokine receptor interaction (ko04060), Neuroactive ligand-receptor interaction (ko04080), Pathways in cancer (ko05200), Olfactory transduction (ko04740), PI3K-Akt signaling pathway (ko04151), Calcium signaling pathway (ko04020), Human papillomavirus infection (ko05165), Jak-STAT signaling pathway (ko04630), and Systemic lupus erythematosus (ko05322) pathways were significantly enriched. The 13 genes among 125 DEGs (22.40%) in Profile 6, and five genes accounting for 11.63% of 43 DEGs in Profile 7 were annotated to the cytokine-cytokine receptor interaction. The three genes among 125 DEGs (2.40%) in Profile 6, and two genes accounting for 4.65% of 43 DEGs in Profile 7 were annotated to the PI3K-Akt signaling pathway. The seven genes among 125 DEGs (5.60%) in Profile 6, and two genes accounting for 4.65% of 43 DEGs in Profile 7 were annotated to the Jak-STAT signaling pathway.

Table 4

pathwaysNo. of DEGs with pathway annotationPathway ID
All profiles
(940)
Profile
6 (125)
Profile
7 (43)
Metabolic pathways162 (17.23%)28 (22.40%)5 (11.63%)ko01100
Cytokine-cytokine receptor interaction81 (8.62%)13 (10.40%)4 (9.30%)ko04060
Neuroactive ligand-receptor interaction73 (7.77%)9 (7.20%)1 (2.33%)ko04080
Pathways in cancer69 (7.34%)5 (4.00%)4 (9.30%)ko05200
Olfactory transduction67 (7.13%)8 (6.40%)3 (6.98%)ko04740
PI3K-Akt signaling pathway48 (5.11%)3 (2.40%)2 (4.65%)ko04151
Calcium signaling pathway42 (4.47%)5 (4.00%)1 (2.33%)ko04020
Human papillomavirus infection42 (4.47%)2 (1.60%)4 (9.30%)ko05165
Jak-STAT signaling pathway40 (4.26%)7 (5.60%)2 (4.65%)ko04630
Systemic lupus erythematosus37 (3.94%)5 (4.00%)1 (2.33%)ko05322

10 top KEGG pathways with high representation of the DEGs.

Validation of RNA-seq Data Using qRT-PCR

Ten genes (NFKB2, MKNK, MAP2K1, MAP3K8, RRAS, TRAF1, NFKBIA, IL8, TNFAIP3, and TGFB1) were selected for qRT-PCR analysis.The results showed a strong correlation between the RNA sequencing data and the qRT-PCR data (Figure 7). This suggested that the expression results generated by RNA sequencing were reliable.

Figure 7

Discussion

In the current piglet production system, early weaning is an important means to improve the efficiency of pig production (24, 25). However, early weaning is very stressful to piglets, and can easily cause piglet immune stress and affect the healthy growth of piglets. How to alleviate the immune stress of piglets has become an area of concerning the piglet industry (26, 27). Hippophae rhamnoides extracts are widely used to enhancing immunity in both healthy and diseased animals (28). Polysaccharide is the main active ingredient of H.rhamnoides. HRP has been shown to have immunomodulatory effects (10). No research has been conducted on the molecular mechanism of HRP in piglets by transcriptome sequencing. Variation in gene expression may provide a key to uncovering the mechanisms of diseases. Transcriptome sequencing, also called RNA-seq, provides a new technique to quantify whole-genome expression profiling in any organism. It promises digital transcriptome profiling with high resolution and is rapidly replacing microarray technology (29, 30). Studies have used RNA-seq technology to explore the protective mechanism of IPEC-J2 cells stimulated by LPS after astragalus polysaccharide (APS) pretreatment. APS relieves cell damage by inhibiting the activation of the MAPK and NF-κB inflammatory pathways, thereby reducing intestinal inflammation (31). Cluster analysis revealed 3134 DEGs that were differentially expressed in H2-L, H4-L, and H6-L into eight possible expression profiles, in which 168 were clustered into two profiles. The up-regulated genes of Profile 6, only between H2-L and H4-L, revealed that these DEGs played a key role in the anti-inflammatory process. The consistent up-regulation of genes of the Profile 7 in H2-L, H4-L, and H6-L indicated that DEGs may contribute to stimulatory functions during the polysaccharide anti-inflammatory process. KEGG enrichment analysis found that the six identified pathways were related to the immune system. Among the six identified pathways related to the immune system, the MAPK signaling pathway and NF-κB signaling pathway play important immunomodulatory roles in our study. Finally, we selected 10 DEGs (NFKB2, MKNK2, MAP2K1 (MEK1), MAP3K8, RRAS, TRAF1, NFKBIA, IL8, TNFAIP3, and TGFB1) related to the main immune pathways to validate the RNA-Seq data using qRT-PCR.

The downstream signal transduction pathways mediated by LPS mainly include the NF-κB signal transduction pathway and the MAPK signal transduction pathway (32). Extracellular regulated protein kinases (ERK), c-Jun N-terminal kinase (JNK), and p38 mitogen-activated protein kinase (p38 MAPK) belong to the three subtypes of the MAPK signaling pathway (33). The Ras/Mitogen-activated protein kinase kinase (MEK)/ERK pathway is one of the most important signal transduction pathways among the MAPK pathways. The Ras/MEK/ERK pathway involves the regulation of a variety of physiological functions of cells and plays a key role in the pathogenesis and pathophysiology of various diseases (34). The activation of ERK is a key step in transferring signals from surface receptors to the nucleus. The activation of ERK induced by LPS leads to the secretion of large amounts of tumor necrosis factor-alpha (TNF-α), as well as interleukin-6 (IL-6) and IL-8, and increases the expression of inducible nitric oxide synthase and nitric oxide. Son of Sevenless (SOS) binds to RRAS-Guanosine diphosphate (GDP), prompting guanosine triphosphate (GTP) to replace GDP on RRAS and activate RRAS protein, then activate MEK and ERK sequentially. In recent years, some initial reports on targeting SOS to inhibit the activation of RRAS, thereby inhibiting the activation of the MAPK signaling pathway, have also achieved satisfactory results (35). The activation of ERK promotes the secretion of MKNK. MKNK is an important downstream protein kinase of ERK that has an important immunomodulatory function. MKNK dysfunction can inhibit the inflammatory signal of upstream ERK and affect downstream eIf4E and CREB and other effector proteins, thereby preventing cell inflammation. MAP3K8 is essential for the activation of the intracellular MAPK/ERK pathway induced by LPS in cells (36). Therefore, MAP3K8 is a critical factor for the production of pro-inflammatory cytokines during immune responses (37). Therefore, in our study, the results showed that MKNK2, MAP2K1, MAP3K8, and RRAS gene expression levels were down-regulated after HRP pre-treatment compared with the LPS group. The reduction of MKNK2, MAP2K1, MAP3K8, and RRAS gene levels plays an important immuno-regulatory role in HRP alleviating LPS-induced cell damage.

LPS, a trigger of inflammation, can activate the NF-κB signaling pathway (38). NFKBIA is a specific inhibitor of NF-κB that binds to NF-κB at a resting state to cause NF-κB to enter an inactive state. Phosphorylated NFKBIA is separated from NF-κB, and NF-κB is activated. Activated NF-κB migrates to the nucleus, where NF-κB nuclear transcription factor can up-regulate the levels of inflammation-related genes TNF, IL-6, and IL-8 and down-regulate the level of TGFB1 (39). Studies have proven that LPS promotes the degradation of NFKBIA, activates the DNA binding ability of NF-κB, and regulates the gene expression level of cytokines (40). Moreover, TRAFs are key regulatory proteins in NF-κB signaling pathways. TRAF1 enhances the activation of TNF-R2 induced by NF-κB (41), therefore promoting the release of a large number of inflammatory cytokines. A previous study showed that TRAF1 is over-expressed in a variety of lymphoma and leukemia cell lines and is a crucial mediator of diverse oncogenic signaling in the development of lymphoid malignancies. TNFAIP3 is a cytokine-induced protein that inhibits apoptosis and activates NF-κB (42). The main function of TNFAIP3 is to inhibit the activity of NF-κB and inhibit TNF-mediated apoptosis, thereby having an important impact on immune regulation and inflammatory processes (43, 44). The release of NFKB2, TRAF1, NFKBIA, IL8, and TNFAIP3, as important regulatory genes of NF-κB, can activate the upstream pathway NF-κB. Jayashankar et al. found that the intervention of supercritical carbon dioxide extract from seabuckthorn leaves can inhibit the expression levels of TNF-α and IL-6 after LPS-induced inflammatory damage, inhibit the activation of the MAPK/NF-κB signaling pathway, reduce inflammation, and play a immunomodulatory role (45). Our study indicated that NFKB2, TRAF1, NFKBIA, IL8, and TNFAIP3were increased and TGFB1 was reduced after LPS induction. However, after HRP pre-treatment, the gene expression level showed the opposite trend and they played an important immune-regulatory role in HRP alleviating LPS-induced cell damage, which provided more targets and prevention directions for theoretical and basic research on intestinal health. Studies have shown that APS may block radiation-induced bystander effects (RIBE) in bone mesenchymal stem cells (BMSCs) induced by the irradiated A549 through regulating the MAPK/NF-κB pathway (46).

Conclusions

This study was the first using a RNA-Seq technique to establish a dynamic transcriptomic profile of three stages (C, L, and H6-L) related to pre-treatment with HRP followed by challenge with LPS in IPEC-J2 cells. Subsequently, bioinformatics analysis (GO, KEGG, and series cluster) helped us to identify key regulatory genes (IL8 and NFKB2, among others.) related IPEC-J2 cellular immune regulation. Transcriptome analysis also showed that HRP protected IPEC-J2 cells from LPS-induced inflammation and decreased the expression of inflammatory cytokines by mainly inhibiting the MAPK/NF-κB signaling pathway. This study does not only provide the useful transcriptomic reference for HRP to effectively protect LPS-induced inflammatory damage in IPEC-J2 cells, but also provides a benchmark for the discovery of biomarkers related to HRP immune regulation.

Funding

This project was supported by grant from the National Key Research and Development Program of China (2017YFD0500506) and the Personnel Foundation of Heilongjiang Bayi Agricultural University (NO. XYB202015).

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.

Statements

Data availability statement

The data presented in the study are deposited in the following repository: https://www.ncbi.nlm.nih.gov/Traces/study/, accession number PRJNA854604.

Author contributions

ML, HS, and LZ contributed to conception, design of study, drafting the manuscript, and critical revision. YZ conducted acquisition of data. LC conducted analysis of data. All authors read and approved the final manuscript.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.944390/full#supplementary-material

References

  • 1.

    HeQTangHRenPKongXWangY. Dietary supplementation with l-arginine partially counteracts serum metabonome induced by weaning stress in piglets. J Prot Res. (2011) 10:521421. 10.1021/pr200688u

  • 2.

    ShanahanF. Probiotics in inflamatory bowel disease. Gut. (2001) 48:609. 10.1136/gut.48.5.609

  • 3.

    HoqueMASkerrattLFRahmanMARabiul Alam BegABMDebnathNC. Factors limiting traditional household duck production in Bangladesh. Trop Ani Health Prod. (2010) 42:157987. 10.1007/s11250-010-9609-z

  • 4.

    LakatosLRednikA. Astroenterologic sub-acute care unit at the hospital department of internal medicine letter. Orvosi Hetilap. (1997) 138:1668.

  • 5.

    BerschneiderM. Development of a normal cultured small intestinal epithelial cell line which transport na and cl. Gastroenterology. (1989) 96:A41.

  • 6.

    HermesRGManzanillaEGSusanaMMartín-OrúeJoséFPérezet al. Influence of dietary ingredients on in vitro inflammatory response of intestinal porcine epithelial cells challenged by an enterotoxigenic escherichia coli (k88). Compar Immunol Microbiol Infect Dis. (2011) 34:47988. 10.1016/j.cimid.2011.08.006

  • 7.

    PanLQinGZhaoYWangJLiuFCheD. Effects of soybean agglutinin on mechanical barrier function and tight junction protein expression in intestinal epithelial cells from piglets. Int J Molecul. (2013) 14:21689704. 10.3390/ijms141121689

  • 8.

    MlotC. Antidotes for antibiotic use on the farm. BioScience. (2000) 50, 955960. 10.1641/0006-3568(2000)050(0955:AFAUOT)2.0.CO;2

  • 9.

    VanCottJLKobayashiTYamamotoMPillaiSMcGheeJRKiyonoH. Induction of pneumococcal polysaccharide-specific mucosal immune responses by oral immunization. Vaccine. (1996) 14:3928. 10.1016/0264-410X(95)00198-A

  • 10.

    ZhangASunHWangX. Recent advances in natural products from plants for treatment of liver diseases. Eur J Med Chem. (2013) 63:5707. 10.1016/j.ejmech.2012.12.062

  • 11.

    WangXLiuJZhangXZhaoSZouKXieJ. Seabuckthorn berry polysaccharide extracts protect against acetaminophen induced hepatotoxicity in mice via activating the Nrf-2/HO-1-SOD-2 signaling pathway. Phytomedicine. (2018) 38:907. 10.1016/j.phymed.2017.11.007

  • 12.

    SuryakumarGGuptaA. Medicinal and therapeutic potential of sea buckthorn (hippophae rhamnoides l.). J Ethnopharmacol. (2011) 138:26878. 10.1016/j.jep.2011.09.024

  • 13.

    Zhang W ZhangXZouKXieJZhaoSLiuJ. Seabuckthorn berry polysaccharide protects against carbon tetrachloride-induced hepatotoxicity in mice via anti-oxidative and anti-inflammatory activities. Food Func. (2017) 8:313038. 10.1039/C7FO00399D

  • 14.

    ShiLFangBYongYLiXGongDLiJ. Chitosan oligosaccharide-mediated attenuation of LPS-induced inflammation in IPEC-J2 cells is related to the TLR4/NF-κB signaling pathway. Carbohyd Poly. (2019) 219:26979. 10.1016/j.carbpol.2019.05.036

  • 15.

    ZhaoLLiMSuSGengTSunH. Hippophae rhamnoides Linn polysaccharide enhances antioxidant enzyme activity, cytokine level and related mRNA expression in intestinal porcine epithelial cells. Canadian J Ani Sci. (2020) 100:193204. 10.1139/cjas-2019-0134

  • 16.

    ZhaoLLiMSunKSuSGengTSunH. Hippophae rhamnoides polysaccharides protect ipec-j2 cells from lps-induced inflammation, apoptosis and barrier dysfunction in vitro via inhibiting tlr4/nf-κb signaling pathway. Int J Biol Macromol. (2020) 155:120215. 10.1016/j.ijbiomac.2019.11.088

  • 17.

    ShaoDHuYWangQTongHShiS. Transcriptome sequencing reveals genetic mechanisms of reproduction performance stimulated by dietary daidzein in laying breeder Hens Theriogenol. (2019) 142:12030. 10.1016/j.theriogenology.2019.09.040

  • 18.

    GaoYLiSBaoXLuoCYangHWangJ. Transcriptional and proteomic analysis revealed a synergistic effect of aflatoxin m1 and ochratoxin a mycotoxins on the intestinal epithelial integrity of differentiated human Caco-2 cells. J Prot Res. (2018) 17:312842. 10.1021/acs.jproteome.8b00241

  • 19.

    LiuHZhangWDongS. Protective effects of sea buckthorn polysaccharide extracts against LPS/d-GalN-induced acute liver failure in mice via suppressing TLR4-NF-κB signaling. J Ethnopharmacol. (2015) 176(Complete):6978. 10.1016/j.jep.2015.10.029

  • 20.

    NiWGaoTWangHDuYLiJLiC. Anti-fatigue activity of polysaccharides from the fruits of four Tibetan plateau indigenous medicinal plants. J Ethnopharmacol. (2013) 150:529535. 10.1016/j.jep.2013.08.055

  • 21.

    ZhaoLGengTSunKSuSZhaoYBaoN. Proteomic analysis reveals the molecular mechanism of Hippophae rhamnoides polysaccharide intervention in LPS-induced inflammation of IPEC-J2 cells in piglets. Int J Biol Macromol. (2020) 164:3294304. 10.1016/j.ijbiomac.2020.08.235

  • 22.

    LangmeadBSalzbergSL. Fast gapped-read alignment with Bowtie 2. Nat Meth. (2012) 9:3579. 10.1038/nmeth.1923

  • 23.

    LiBDeweyCN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. (2011) 12:1. 10.1186/1471-2105-12-323

  • 24.

    BombaLMinutiAMoisáSJTrevisiEEufemiELizierM. Gut response induced by weaning in piglet features marked changes in immune and inflammatory response. Func Int Genom. (2014) 14:65771. 10.1007/s10142-014-0396-x

  • 25.

    KhafipourEMunyakaPMNyachotiCMKrauseDORodriguez-LecompteJC. Effect of crowding stress and escherichia coli k88+ challenge in nursery pigs supplemented with anti-escherichia coli k88+ probiotics. J Animal Sci. (2014) 92:201729. 10.2527/jas.2013-7043

  • 26.

    PluskeJRHampsonDJWilliamsIH. Factors influencing the structure and function of the small intestine in the weaned pig: a review. Livestock Prod Sci. (1997) 51:21536. 10.1016/S0301-6226(97)00057-2

  • 27.

    XiaoT. Innate immune recognition of nucleic acids. Immunol Res. (2009) 43:98108. 10.1007/s12026-008-8053-x

  • 28.

    TanwarHShweta SinghDSinghSBGanjuL. Anti-inflammatory activity of the functional groups present in Hippophae rhamnoides (Seabuckthorn) leaf extract. Inflammopharmacology. (2017) 26:291301. 10.1007/s10787-017-0345-0

  • 29.

    HerwigR. Predictive network modelling with toxicogenomics data. Toxicol Lett. (2014) 229:S4S21. 10.1016/j.toxlet.2014.06.043

  • 30.

    DavoliRZambonelliPHedeegardJHornshojHNanni CostaLStellaA. Transcriptome analysis of skeletal muscle tissue to identify genes involved in pre-slaughter stress response in pigs. Italian J Anil Sci. (2009) 8:6971. 10.4081/ijas.2009.s2.69

  • 31.

    DongNLiXXueCZhangLWangCXuX. Astragalus polysaccharides alleviates lps-induced inflammation via the nf-κb/mapk signaling pathway. J Cell Physiol. (2020) 235:116. 10.1002/jcp.29452

  • 32.

    WuCXSunHLiuQGuoHGongJP. Lps induces hmgb1 relocation and release by activating the nf-κb-cbp signal transduction pathway in the murine macrophage-like cell line raw264.7. J Surg Res. (2012) 175:88100. 10.1016/j.jss.2011.02.026

  • 33.

    CostaAPLopesMWRiegerDKBarbosaSGRGonçalvesFMXikotaJC. Differential activation of mitogen-activated protein kinases, ERK 1/2, 38MA.PK and JNK p54/p46 during postnatal development of rat hippocampus. Neurochem Res. (2015) 41:11609. 10.1007/s11064-015-1810-z

  • 34.

    OkanoJISnyderLCRustgiAK. Paclitaxel induces prolonged activation of the ras/mek/erk pathway independently of activating the programmed cell death machinery. Gastroenterology. (2011) 120:A6612. 10.1016/S0016-5085(01)83290-X

  • 35.

    ZhaoHLiYYFuciniRVRossSEPessinJEKoretzkyGA. T cell receptor-induced phosphorylation of sos requires activity of cd45, lck, and protein kinase c, but not erk. J Biol Chem. (1997) 272:2162534. 10.1074/jbc.272.34.21625

  • 36.

    SantagSSiegelFWengnerAMLangeCPetersenK. Abstract 341: preclinical mode of action and anti-tumor efficacy of the selective mknk1 inhibitor bay 11,43,269 in nsclc models. Cancer Res. (2016) 76(14 Supplement):341. 10.1158/1538-7445.AM2016-341

  • 37.

    WeiRYangQHanBLiYYaoKYanX. Microrna-375 inhibits colorectal cancer cells proliferation by downregulating jak2/stat3 and map3k8/erk signaling pathways. Oncotarget. (2017) 8:1663341. 10.18632/oncotarget.15114

  • 38.

    ShishodiaSKoulDAggarwalBB. Cyclooxygenase (cox)-2 inhibitor celecoxib abrogates tnf-induced nf-κb activation through inhibition of activation of iκbα kinase and akt in human non-small cell lung carcinoma: correlation with suppression of cox-2 synthesis. J Immunol. (2004) 173:201122. 10.4049/jimmunol.173.3.2011

  • 39.

    KohY. Inhibition of SRC Tyrosine kinases suppresses activation of nuclear factor-κB, and Serine and Tyrosine Phosphorylation of IκB-α in Lipopolysaccharide-stimulated Raw 264.7 macrophages. J Toxicol Environ Health, Part A. (2005) 68:164362. 10.1080/15287390500192114

  • 40.

    LiMYSunLNiuXTChenXMTianJXKongYDet al. Astaxanthin protects lipopolysaccharide-induced inflammatory response in Channa argus through inhibiting NF-κB and MAPKs signaling pathways. Fish Shellfish Immunol. (2019) 86:2806. 10.1016/j.fsi.2018.11.011

  • 41.

    TangXZhangLWeiW. Roles of TRAFs in NF-κB signaling pathways mediated by BAFF. Immunol Lett. (2018) 196:1138. 10.1016/j.imlet.2018.01.010

  • 42.

    RheeLMurphySFKolodziejLEGrimmWAWeberCRLodolceJP. Expression of TNFAIP3 in intestinal epithelial cells protects from DSS- but not TNBS-induced colitis. Am J Physiol-Gastroint Liver Physiol. (2012) 303:G2207. 10.1152/ajpgi.00077.2012

  • 43.

    GuiJYueYChenRXuWXiongS. A20 (TNFAIP3) alleviates CVB3-induced myocarditis via inhibiting NF-κB signaling. PLoS ONE. (2012) 7:e46515. 10.1371/journal.pone.0046515

  • 44.

    VarfolomeevEGoncharovTMaeckerHZobelKKomuvesLGDeshayesKet al. Cellular inhibitors of apoptosis are global regulators of NF- B and MAPK activation by members of the TNF family of receptors. Sci Signal. (2012) 5:ra22. 10.1126/scisignal.2001878

  • 45.

    JayashankarBMishraKPKumarMSYUdayasankarKMisraKGanjuLet al. A supercritical CO2 extract from seabuckthorn leaves inhibits pro-inflammatory mediators via inhibition of mitogen activated protein kinase p38 and transcription factor nuclear factor-κB. Int Immunopharmacol. (2012) 13:4617. 10.1016/j.intimp.2012.05.011

  • 46.

    ZhangLLuoYLuZHeJLiuY. Astragalus polysaccharide inhibits ionizing radiation-induced bystander effects by regulating mapk/nf-kb signaling pathway in bone mesenchymal stem cells (bmscs). Med Mon Int Med J Exp Clin Res. (2018) 24:4649. 10.12659/MSM.909153

Summary

Keywords

IPEC-J2 cells, HRP, anti-inflammatory, MAPK/NF-κB signaling pathway, transcriptome

Citation

Li M, Chen L, Zhao Y, Sun H and Zhao L (2022) Research on the Mechanism of HRP Relieving IPEC-J2 Cells Immunological Stress Based on Transcriptome Sequencing Analysis. Front. Nutr. 9:944390. doi: 10.3389/fnut.2022.944390

Received

15 May 2022

Accepted

13 June 2022

Published

15 July 2022

Volume

9 - 2022

Edited by

Baojun Xu, United International College, China

Reviewed by

Utoomporn Surayot, Chiang Mai University, Thailand; Susmita Barman, University of Nebraska Medical Center, United States

Updates

Copyright

*Correspondence: Lei Zhao Hui Sun

This article was submitted to Nutritional Immunology, a section of the journal Frontiers in Nutrition

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics