Animals. Four-week-old SPF Kunming mice (n = 28; average weight, 20 g) were gained from the Slack Jingda Experimental Animal Co, Ltd. (Hunan, China). According to our previous research, male mice are used in this experiment [18]. The mice were maintained in a shielded environment at the Animal Experiment Center of the Hunan University of Chinese Medicine (Hunan, China, SYXK [Xiang] 2015-0003). Each mouse was housed individually under standard controlled laboratory conditions(12/12 dark-light cycle, temperature of 21 ± 2°C, humidity of 45 ± 10%), with ad libitum access to food and water.
Experiment design. The animal experiment was designed to intervene in intestinal microbiota homeostasis in the host and determine the role of intestinal microbiota homeostasis in host adaptive regulation of GI disorders caused by cold and humid environmental stress. After a week to habituate freely for acclimatization, all animal randomized grouping and treatments are described in Fig. 1. Firstly, to confirm the role of intestinal microbiota by the intestinal microbiota disorder, according to our previous studies [15–17] the mice were treated with antibiotics mixture to cause the intestinal microbiota disorder. The other mice were given the same dose of normal saline by gavage. The cold and humid environmental stress intervention in mice with intestinal microbiota homeostasis and intestinal microbiota disorder, forming the cold and humid environmental stress treatment group (CW-M) and microbiota disorder + cold and humid environmental stress treatment group (MD + CW), respectively. The cold and humid environmental stress intervention method: the mice were maintained in an artificial climate room from a.m. 8:00 to 12:00 with a temperature of 4 ± 0.5°C, and a humidity of 90 ± 2%, and in the remaining time, the model mice were fed in the normal shielded environment, for continual 7 days. In the meantime, the normal control group (CW-C) mice were fed in the normal shielded environment. Then, to study the metabolic mechanism of intestinal mucosal microbiota by reproducing cold and humid environmental stress animal models.
During the artificial climate room treatment, the feces of the 3 groups of mice were collected at the same time every day to determine moisture content and calculate the diarrhea index [19, 20].
Diarrhea index = l1 * L (1)
where l1 is loose stool rate (number of loose feces per mouse / total feces) and L is loose stool grade. The area of dirt formed by dilute fecal pollution filter paper is classified into 4 loose stool grades. Grade 1: pollution diameter < 1cm; Grade 2: 1cm ≤ pollution diameter < 2cm; Level 3: 2cm ≤ pollution diameter ≤ 3cm; Grade 4: pollution diameter > 3cm.
Body mass and food intakes of all experimental mice were measured every day at baseline until the end of the study, to calculate the food efficiency ratio[21, 22].
Food efficiency ratio = M1 / M * 100 (2)
where M1 is average daily weight gain (g) and M is average daily food intake(g).
After the treatment, mice were fasted for 12 h before taking samples, and the blood and intestinal tissue samples were extracted.
Surgical sampling of the intestinal mucosa.
Bacterial culturing. Fecal samples of mice were collected at an interval of 1 day until the end of the study, and samples of intestinal contents from the duodenum to the ileocecal segment were extracted, which were used for bacterial culture. Fully mix the same kind of samples in the same group, weigh 0.5 g of mixed sample, dilute with 30 mL sterile distilled water, and then shake at 120 R / min for 30 min to prepare a bacterial solution. Beef extract peptone AGAR medium was used for the bacterial culture. The bacterial solution was diluted with appropriate dilution and coated with the plate. Finally, the Petri dishes were cultured in the 37 ℃ incubator and anaerobic tank for 48 hours and counted. Repeat 3 times for each dilution. The viable bacteria counting formula was used to count, and the results were expressed in the form of colony value per gram of fecal wet weight to calculate the total number of colonies (CFU / g).
FDA microbial activity. FDA with the characteristic that can freely penetrate into and out of the complete cell membrane and be decomposed by lipase in the cell to produce polar fluorescein that can produce fluorescence was used to measure the microbial activity at 490 nm wavelength. FDA mixed acetone to prepare a 2 mg / mL FDA storage solution, which was stored in a refrigerator at 2–8 ℃ away from light. The FDA storage solution was transferred into PBS (pH 7.6) to prepare a concentration of 10 µ g / mL of FDA reaction solution. Transfer 2 ml of FDA reaction solution and 50 µ L sample diluent in a sterile test tube, and repeat 3 times. The sample was shaken at 24 ℃ for 90 min, and the reaction was terminated by adding 2 mL of acetone. The OD was measured with an ultraviolet spectrophotometer at λ490 nm. The definition of 1 microbial activity unit is the enzyme required that 1 g sample reacts for 90 min to produce 1 µ mol of fluorescein.
Mucus barrier. The fresh colon at 1 cm from the rear end of the mouse cecum was dissected and sectioned with hematoxylin-eosin (HE) staining and alcian blue - periodic acid-Schiff (AB – PAS) staining. The histological morphology of the mouse colon was observed by an optical microscope, and the sections were digitally scanned. The images were processed by Caseviewer: C.V 2.3 system and the morphological changes of colonic mucosa were compared. 10 visual fields of different samples were randomly selected to count and calculate the number of goblet cells and the thickness of the mucus layer of the intestinal mucosa.
Muc2 in colonic tissue of mice was detected by enzyme-linked immunosorbent assay (ELISA). The intestinal tissue about 3 cm from the back of the cecum of mice was cut by surgery, weighed, labeled, and quick-frozen in dry ice. After grinding frozen tissue, tissue lysate was added to prepare intestinal tissue homogenate samples. The tissue samples were detected according to the steps of the mouse Muc2 ELISA kit. Mouse Muc2 ELISA kit is provided by Jiangsu Sofia Medical Technology Co., Ltd (20200397).
Amino acid transferases and blood lipids. The mice fasted for 12 hours before sampling, and eyeball blood was collected into the procoagulant tube. The whole blood of mice was centrifuged at 3000 R / min− 1 for 10 min, and the serum was taken for standby. Amino acid transferases in serum samples were determined by the automatic biochemical instrument, such as glutamate pyruvic transa (GPT), glutamic oxaloacetic transaminase (GOT), alkaline phosphatase (ALP), and γ-glutamyltransferase (GGT). Blood lipids in serum samples were determined by the automatic biochemical instrument, such as serum total cholesterol (TC), triacylglycerol (TG), high-density leptin cholesterol (HDL-C), and low-density leptin cholesterol (LDL-C).
Full-length PCR amplication. The bacterial DNA was extracted according to the steps of the HiPure Soil DNA Mini Kit (D3142-03). The concentration of DNA was determined by NanoDrop ND-2000 (Thermo Fisher Scientific, Waltham, MA, US), and the qualified DNA was amplified by PCR. The PCR amplication primers of the bacterial full-length 16S rRNA genes:
forward: 27F (5´-AGRGTTYGATYMTGGCTCAG-3´)
reverse: 1492R (5´-RGYTACCTTGTTACGACTT-3´)
Sample-specific16-bp barcodes were incorporated into the primers for multiplex sequencing. PCR reaction procedure: a) 3 minutes at 95 ℃; b) 25 cycles×(30 sec at 95 ℃; 30 sec at 56 ℃; 60 sec at 72 ℃); c) 5 minutes at 72 ℃. PCR amplicons were puried with Agencourt AMPure Beads (Beckman Coulter, Indianapolis, IN) and quantied BY the PicoGreen dsDNA AssayKit (Invitrogen, Carlsbad, CA, USA).
PacBio HiFi sequencing. The PCR products were end-repaired and spliced to generate a library by SMRTbell Template Prep Kit 1.0-SPv3 and accurately quantified by Qubit. Agilent 2100 was used to detect the size of the inserted fragment, and the expected DNA was sequenced on the PacBio platform. Sequencing was performed using the PacBio platform with DNA/Polymerase Binding Kit 3.0 (PacBio). Software SMRT link V8 0 was used to preprocess and filter the original PacBio sequencing output data to obtain circular consensus sequences (CCS). The clean CCS sequence was filtered by removing the CCS primer and unqualified length. Finally, the intestinal mucosa samples with an average length of CCS less than 1400 dpi were excluded for subsequent analysis to ensure the accuracy of the data.
Taxonomic and biomarker analysis. Qiime2 software (qiime. ORG) was used to cluster the clean CCS sequences of all samples with 97% identity. RDP classifier Bayesian algorithm was used to classify 97% of the OTU representative sequences at similar levels (with a confidence threshold of 0.8), and OTU annotation was compared with the Silva database. Mothur software (http://www.mothur.org/wiki/) calculated alpha diversity in this study, and the difference in alpha diversity between groups of mice was analyzed by the density curve in R v4.1.0. A significant difference in beta diversity between groups were tested using a one-way analysis of variance followed by post hoc Tukey honestly significant difference at a 95% confidence level. The composition of dominant microbiota was performed by GraphPad Prism v9.0.0, and the differences between the two groups were analyzed by the Wilcoxon rank-sum test. LDA effect size (LEfSe) analysis was used to find out the biomarkers in sample division (http://huttenhower.sph.harvard.edu/galaxy). Metastat analysis was used to detect significant microbial community differences between the two groups (http://metastats.cbcb.umd.edu/detection.html).
Functional analysis of microbiota metagenome. Functional profiling was performed using PICRUSt2 (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States) v2.0.0-b with built-in EC/MetaCyc and KEGG/KO databases. The functional abundance of the sample is predicted based on the abundance of marker gene sequences in the sample. The platform inferred the metabolic pathways by using MinPath to obtain the abundance data of metabolic pathways. MetaCyc and KEGG pathway abundances were further analyzed using statistical analysis of Omic Share profiles (http://omicshare.com/tools/) and R v4.1.0 ggplot2 package. The correlation coefficient was calculated by using the cor () function in R v3.1.3. According to the P ≤ 0.05 and the correlation coefficient R ≥ 0.7, the relationship network between mucosal microbiota and amino acid metabolic function was constructed by using Cytoscape v3.9.0.
Metabonomic amino acid metabolites standard curve and limit of quantitation. Weigh an appropriate amount of 23 amino acid standards and prepare single standard mother liquor with methanol or water. The profiling of amino acid information and concentration points can be referred to (Supplementary Table S1). A liquid Chromatograph Mass Spectrometer (LC-MS) was used to detect each working standard solution and draw the standard curve. ACQUITY UPLC BEH C18 chromatographic column (2.1×100 mm, 1.7 µm) was set to 40 ℃ and a 5 µL standard solution was added. Mobile phase A-10% methanol-water (containing 0.1% formic acid), B-50% methanol-water (containing 0.1% formic acid). Gradient elution conditions were 0 ~ 6.5 min, 10 ~ 30% B; 6.5 ~ 7 min, 30 ~ 100% B; 7 ~ 8 min, 100% B; 8 ~ 8.5 min,100 ~ 10% B; 8.5 ~ 12.5 min, 10% B. Flow rate: 0 ~ 8.5 min, 0.3 mL / min; 8.5 ~ 12.5 min, 0.3 ~ 0.4 mL/min. The mass spectrometric condition was electrospray ionization (ESI) source, using positive ion ionization mode. The ion source temperature was 500 ℃, the ion source voltage was 5500 V, the collision gas was 6 psi, the air curtain gas was 30 psi, and the atomization gas and auxiliary gas were 50 psi. Multiple response monitoring (MRM) was used for scanning. The ions for quantitative analysis can be referred to (Supplementary Table S2). The concentration of the standard solution was taken as the abscissa and the area ratio of the working standard solution to the internal standard peak was taken as the ordinate to investigate the linear range and draw the standard curve (Supplementary Table S3). Eight standard solution samples reprocessed to determine intra-day and inter-day precision, repeatability, and recovery can be referred to (Supplementary Table S3).
Amino acid quantitative profiling. Transfer the mouse serum sample into a 2 ml EP tube and dilute it with 10% formic acid methanol solution ddH2O (1:1. V: V) solution. The diluted sample of 100 µ L is added to a double isotope internal standard with a concentration of 100 ppb µ L vortex oscillation for the 30s. The supernatant was filtered with a 0.22 µ mol/L membrane, and the filter solution was detected. All samples were analyzed quantitatively according to the established sample pretreatment and instrumental analysis methods. The data set was scaled by the “pheatmap” package in R v3.3.2 to export the hierarchical cluster diagram of relative quantitative values of amino acid metabolites. The metabolomics data were scaled by autoscaling, mean centering, and scaled to unit variance (UV), and then subjected to multivariate statistical analysis of partial least squares discriminant analysis (PLS-DA) and orthogonal partial least squares discriminant analysis (OPLS-DA) by using “ropls” package of R v3.3.2. The differential metabolites were screened according to the statistically significant p-value and fold change. Z-score was converted based on the content of metabolites and was used to measure the variety of metabolites at the same level.
Bioinformatics analysis. The co-occurrence network of intestinal mucosal microbiota was constructed to interpret the core bacteria in the relationship. According to Spearman’s correlation coefficient of intestinal mucosal microbiota, the positive or negative correlation between different microbial communities was investigated. Pearson correlation analysis was used to construct the relationship network of intestinal mucosal microbiota targeted amino acid metabolism in mice, to integrate the relationship between intestinal mucosal microbiota species and serum amino acid metabolites. The relationship pairs were screened in the conditions of the correlation coefficient R ≥ 0.7 and FDR < 0.05. According to the ranking of species abundance or metabolite concentration in all samples, Cytoscape v3.9.0 software was used to perform the correlation network of intestinal mucosal microbiota. The neural network algorithm was used to distribute the network nodes. The cor. test () function of R v3.1.3 was used to statistically test the correlation analysis of amino acid metabolites, and the heat map of amino acid metabolites was exported by using the “pheatmap” package of R v3.1.3 to reveal the synergy of metabolite variety. To explain the relationship between intestinal mucosal microbiota structure and amino acid metabolic environment, redundancy analysis (RDA) was used to study the multiple linear regression between microbiota response variables and explanatory variables of amino acid metabolites. The length of the arrow represents the constraint effect of different environmental variables exerted on sample distribution in two-dimension. Procrustes analysis was used to study the consistency of microbiota abundance and amino acid metabolome. RDA and Procrustes analysis were performed by using “vegan” and “ggplot2” packages of R v4.1.0. Random forest analysis was performed with the “randomForest” package of R v4.1.0. A mathematical model based on key bacteria and amino acid metabolome based on random forest variables was used to predict Muc2 content in the sampling group. Linear regression analysis was used to determine the quantitative relationship between colon Muc2 concentration and the A mathematical model.
Statistics. The continuous data conform to the normal distribution and was expressed by mean ± SEM. If the variance was homogeneous, the difference between the two groups was expressed by an independent sample t-test. One-way ANOVA was used for multi-group comparison, and the LSD test was used for multi-group pairwise comparison. If it does not conform to the normal distribution, the nonparametric rank-sum test was used. The rate of each group was expressed as a percentage (%), and the rate between groups was compared by the Chi-square test. Inspection level α= 0.05. P < 0.05 indicates significant difference; P < 0.01 indicates that the difference is highly significant.