Zhiqiang Wan1†Xufei Zhang2†Xianhao Jia2†Yuhua Qin3Ning Sun1Jinge Xin4Yan Zeng1
Bo Jing1
Jing Fang1
Kangcheng Pan1
Dong Zeng1
Yang Bai4
Hesong Wang5*
Hailin Ma2*
Xueqin Ni1*- 1Animal Microecology Institute, College of Veterinary, Sichuan Agricultural University, Chengdu, China
- 2Plateau Brain Science Research Center, Tibet University, Lhasa, China
- 3College of Basic Medical Sciences, Chengdu University of Traditional Chinese Medicine, Chengdu, China
- 4Guangdong Provincial Key Laboratory of Gastroenterology, Department of Gastroenterology, Institute of Gastroenterology of Guangdong Province, Nanfang Hospital, Southern Medical University, Guangzhou, China
- 5Guangzhou Beneco Biotechnology Co. Ltd., Guangzhou, China
Background: Intestinal microbiota plays an important role in maintaining the microecological balance of the gastrointestinal tract in various animals. Disturbances in the intestinal microbiota may lead to the proliferation of potentially pathogenic bacteria that become the dominant species, leading to intestinal immune disorders, intestinal inflammation, and other intestinal diseases. Numerous studies have been confirmed that high-altitude exposure affects the normal function of the intestine and the composition of the intestinal microbiota. However, it is still necessary to reveal the changes in intestinal microbiota in high-altitude exposure environments, and clarify the relationship between the proliferation of potentially pathogenic bacteria and intestinal injury in this environment. In addition, explored probiotics that may have preventive effects against intestinal diseases.
Methods and results: C57BL/6 mice were randomly divided into three groups, a high-altitude group (HA), control group (C), and high-altitude probiotic group (HAP). The HA and HAP groups were subjected to hypoxia modeling for 14 days in a low-pressure oxygen chamber with daily gavage of 0.2 mL of normal saline (HA) and Lactobacillus johnsonii YH1136 bacterial fluid (HAP), while the control group was fed normally. L. johnsonii YH1136 was isolated from feces of a healthy Tibetan girl in Baingoin county, the Nagqu region of the Tibet Autonomous Region, at an altitude of 5000 meters. Our observations revealed that gavage of YH1136 was effective in improving the damage to the intestinal barrier caused by high-altitude exposure to hypoxic environments and helped to reduce the likelihood of pathogenic bacteria infection through the intestinal barrier. It also positively regulates the intestinal microbiota to the extent of Lactobacillus being the dominant microbiome and reducing the number of pathogenic bacteria. By analyzing the expression profile of ileal microRNAs and correlation analysis with intestinal microbiota, we found that Staphylococcus and Corynebacterium1 cooperated with miR-196a-1-3p and miR-3060-3p, respectively, to play a regulatory role in the process of high-altitude hypoxia-induced intestinal injury.
Conclusion: These findings revealed the beneficial effect of L. johnsonii YH1136 in preventing potential endogenous pathogenic bacteria-induced intestinal dysfunction in high-altitude environments. The mechanism may be related to the regulation of intestinal injury from the perspective of the gut microbiota as well as miRNAs.
Introduction
Probiotics can maintain the stability of intestinal microbiota in the body, which is of great significance to the health of humans and animals (1). In recent years, there has been an increasing understanding of the use of probiotics for the prevention and treatment of intestinal diseases caused by pathogenic bacteria or intestinal immune dysregulation; however, there is a lack of understanding that exists in high-altitude exposure environments. The features of high-altitude exposure are low air pressure and hypoxia, which profoundly affect the physiological function and status of organs and are considered serious problems (2, 3).Studies have shown that intestinal microbiota disorders are closely related to impaired intestinal function, and acute high-altitude exposure leads to intestinal microbiota disorders that affect intestinal barrier function, which in turn may cause bacterial translocation and, consequently, intestinal inflammation (4–6). These are also considered to be important potential factors in the progression of inflammatory bowel disease (IBD) (6–8), but the mechanisms by which high-altitude hypoxic environments affect intestinal function are not yet clear. Disturbances in the intestinal microbiota can lead to a predominance of pathogenic bacteria in the gut, resulting in intestinal infections. Pan et al. (9) reported that the gut microbiota in a hypoxic environment appeared different, characterized by an increased abundance of the genera Parabacteroides, Alistipes and Lactococcus and an increased ratio of Bacteroides to Prevotella. In addition, disruption of the gut microbiota increases susceptibility to infection (10), and the intestinal microbiota usually plays an important role in immune regulation in maintaining host immune homeostasis and maintaining intestinal health together with the intestinal barrier (11). Previous studies have revealed that a special living environment makes the gut microbiota of animals in high-altitude areas unique facilitating their adaptation to high-altitude environments (12, 13). This suggests that altitude-related intestinal issues can be investigated from the perspective of the intestinal microbiota. A novel approach to coping with intestinal injury induced by a high-altitude hypoxic environment may be the use of probiotics to modulate the composition of the intestinal microbiota.
However, few studies have examined the effects of probiotics in plateau environment. Based on previous studies, we hypothesized that sudden exposure to a plateau environment disrupts the intestinal microbial balance, where pathogenic bacteria that become dominant species may disrupt intestinal immunity leading to intestinal inflammation, and that probiotic supplementation may have a regulatory effect on this process. We intend to use L. johnsonii YH1136 in further studies. Lactobacillus spp. is a widely used probiotic. We isolated Lactobacillus johnsonii YH1136 (CCTCC M 20221116) from feces of a healthy Tibetan girl in Baingoin county, the Nagqu region of the Tibet Autonomous Region, at an altitude of 5000 meters. Based on our previous studies, Lactobacillus johnsonii shows promising results by improving intestinal microbiota in preventing and treating nonalcoholic fatty liver disease, modulating stress and memory impairment based on the gut-brain axis, and plays an important regulatory role in processes involved in renal and intestinal injury induced by high fluoride exposure (10, 14–19). Therefore, this study evaluated whether L.johnsonii YH1136 is effective in mitigating damage to intestinal integrity due to hypoxic environments in rats exposed to high altitudes. In addition to the intestinal microbiota, microRNAs (miRNAs) also play regulatory role during intestinal injury. It has been reported that miRNAs can regulate proteins involved in intestinal barrier integrity and permeability, and are considered biomarkers of some intestinal diseases (20–22). Therefore, we detected and analyzed the expression of miRNAs in the ileum and the composition characteristics of intestinal microbiota concurrently, in order to describe the changes in the intestinal environment after high-altitude hypoxic stress and explore the preventive mechanism of L. johnsonii YH1136.
Materials and methods
Animal and experiment design
A total of 36 C57BL/6 mice of the same age and similar growth conditions (8 weeks old, purchased from Chengdu Dashuo Institute of Biology, China) were randomly selected and, divided into three groups and given a week to adapt to their new environment. The three groups were the control (C), high altitude (HA), and high-altitude probiotic (HAP) groups. The air pressure level for group C was maintained at 94.5 kPa, while the HA and HAP groups were exposed to a low-pressure oxygen chamber simulating an altitude of 3500-4000m for 14 days, with the air pressure in the chamber set to 60-65 kPa. Mice in the HAP group were given 0.2mL L. johnsonii YH1136 by gavage every day, while mice in the C and HA groups were given 0.2mL physiological saline (pH 7.0) by gavage as a substitute. The temperature and humidity of the low-pressure oxygen chamber were identical to those outside the chamber. The chamber was opened for one hour per day to supply food (Chengdu Dashuo Institute of Biology, China) and water. All animal experiments were conducted in accordance with the guidelines for the feeding and use of experimental animals (approval no: SYXKchuan2019-187) approved by the Committee of Sichuan Agriculture University. All animals were raised in the animal room of the Animal Microecology Research Institute of Sichuan Agriculture University, with a 12-hour diurnal cycle. The experiment lasted 14 days.
Sample collection
On day 15, six mice from each group were randomly selected and subjected to cervical dislocation. The ileal segment of mice (approximately 0.5 cm) was collected, washed with cold diethylpyrocarbonate water, and stored in liquid nitrogen until further miRNA analysis. The serum was collected by centrifuging the blood of mice and storing it at -80°C. The ileum was collected from each group, fixed in 4% paraformaldehyde solution, and stored at 4°C for immunohistochemical analysis.
16S rRNA gene amplicon sequencing and analysis
PCR amplification of the bacterial 16S rRNA genes V3–V4 region was performed using the forward primer 338F (5’-ACTCCTACGGGAGGCAGCAG-3’) and reverse primer 806R (5’-GGACTACHVGGGTWTCTAAT-3’). The PCR components contained 5μL of buffer (5×), 0.25μL of Fast Pfu DNA Polymerase (5U/μL), 2μL (2.5 mM) of dNTPs, 1μL (10 μM) of each forward and reverse primer, 1 μL of DNA template, and 14.75 μL of ddH2O. Thermal cycling consisted of initial denaturation at 98°C for 5 min, followed by 25 cycles consisting of denaturation at 98°C for 30 s, annealing at 53°C for 30 s, and extension at 72°C for 45 s, with a final extension of 5 min at 72°C. PCR amplicons were purified using Vazyme VAHTSTM DNA Clean Beads (Vazyme, Nanjing, China) and quantified using the Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, USA). After the individual quantification step, amplicons were pooled in equal amounts, and pair-end 2× 250 bp sequencing was performed using the Illlumina MiSeq platform with MiSeq Reagent Kit v3 at Shanghai Personal Biotechnology Co., Ltd (Shanghai, China). Bioinformatics analysis was performed using QIIME2 2020.11 (https://docs.qiime2.org/2020.11/tutorials/) and the (V3.1.2)
Serum biochemical test and statistical analysis
A mouse-specific ELISA kit was used to detect the levels of serum diamine oxidase (DAO) and D-lactate. Statistical analysis of DAO and D-lactate levels was performed by one-way analysis of variance (ANOVA) using IBM SPSS Statistics version 27, to detect the normality and statistical significance of the data. The data was considered to confirm to a normal distribution if the significance value obtained from the Shapiro-Wilk test was > 0.05. One-way ANOVA was used for analysis, and least significant difference (LSD) and Dunnett T3 tests were used for multiple comparisons. Statistical significance was set at P <0.05.
Immunohistochemical and immunofluorescence analysis
The fixed tissues were flushed with water, dehydrated, made transparent, waxed, embedded, and sectioned to prepare paraffin sections. The sections were placed in citric acid antigen retrieval buffer (pH 6.0) in a microwave oven and incubated after boiling for antigen retrieval. The sections were thereafter placed in a 3% H2O2 solution and incubated for 25 min in the dark at room temperature, after which the sections were placed in Phosphate-buffered saline (PBS, pH 7.4) and shaken three times for 5 min each on a shaker. The tissue was thereafter uniformly covered and blocked with 3% bovine serum albumin (BSA) for 30 min at room temperature. The blocking solution was gently discarded, and rabbit polyclonal anti-occludin (1:200), rabbit polyclonal anti-claudin-1 (1:200), rabbit polyclonal anti-ZO-1 (1:200), rabbit polyclonal anti-IL-1β (1:200), rabbit polyclonal anti-TGF-β (1:200), and rabbit polyclonal anti-TNF-α (1:200) antibodies were sequentially added to the sections and incubated overnight at 4°C. The antibody solutions (occludin, Claudin-1 and ZO-1) were removed by washing with PBS, and Horseradish (HRP)-labeled goat antirabbit IgG (1:200) was added dropwise and incubated for 50 min at room temperature. Finally, the cells and nuclei were stained using the immunohistochemical DAB chromogen kit and hematoxylin, respectively. Sections for inflammatory factors (occludin, claudin, and ZO-1) were counterstained for nuclei with DAPI; an autofluorescence quencher was added, and an anti-fade mounting medium was used to block the slides. The prepared sections were placed under a fluorescence microscope for observation and to obtain images. All above reagents were purchased from Servicebio Technology Co., Ltd. (Wuhan, China).
Small RNA library preparation and sequencing
Total RNA was extracted from each group using a Total RNA Isolation Kit (Omega Inc., Norcross, Georgia, United States) according to the manufacturer’s instructions. The library was constructed and sequenced as follows. First, an Agilent 2100 Bioanalyzer was used to detect the RNA concentration and integrity. The TruSeq Small RNA Sample Prep Kit was used to construct a small RNA library. Next, the library was enriched by PCR amplification, and the sequencing adaptor and index parts were added. The library was thereafer purified by gel electrophoresis. An Agilent 2100 Bioanalyzer was used to test the quality of the library with an Agilent High Sensitivity Kit, where the qualified library should have a single peak. The the library was quantified using a Quant-iT PicoGreen dsDNA Assay Kit, and finally, an appropriate loading amount was selected and sequenced on an Illumina platform.
Analysis and identification of microRNA
The raw data was processed to remove the adaptor sequences and filter the unqualified data. We then merged identical sequences and recorded the richness of the remaining sequences. Based on the mouse reference genome, the processed data were annotated and compared to obtain relevant information on small RNA. Small RNA includes miRNAs, piwi-interacting RNA, and non-coding RNAs (ncRNAs). Subsequent analysis was performed mainly on the miRNAs. We downloaded the sequences of mouse precursors and matures miRNA form miRBase, and then compared unique reads with the downloaded sequences to annotate the sequencing results. According to the comparison results, we used mireap to analyze the sequences that have not yet been annotated, and predicted novo miRNAs.
Analysis of known miRNA expression
The read counts of the measured miRNAs were calculated according to the sequence number of mature miRNAs in mice. Because a miRNA may have multiple precursors that produce potentially identical mature miRNAs, we chose to name the abundance of the first occurrence of the same name miRNAs as the abundance of the miRNAs for subsequent analysis.
miRNA target prediction and functional analysis
miRNAs mainly bind to target sites through complementary pairing. We used the miRanda software to predict the target genes of the differentially expressed miRNAs. Using the 3- ‘UTR’ of mouse mRNA as the target sequence, the number of target genes and target sites of differentially expressed miRNAs were predicted. We then used the topGO software for enrichment analysis, where we used Gene Ontology (GO) terms to annotate the different miRNA target genes and calculate the number of miRNA target genes of each term in order to determine which of the different miRNA target genes were significantly enriched in the whole genomic context (P-value < 0.05 indicated significant enrichment). Thus, the main biological functions of differentially expressed miRNA target genes were determined.
The validation of differential expression of miRNA by RT-qPCR
We used the stem-loop RT-qPCR method (23) to verify the differential expression of miRNAs and randomly selected differentially expressed miRNAs from each of the two groups were balanced. Total RNA was extracted using a Total RNA Isolation Kit (Omega Bio-Tek, Norcross, GA, USA) according to the manufacturer’s instructions. The purity and concentration of the RNA were assesed using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific Inc., Waltham, MA, USA), and cDNA was generated by reverse transcription using the M-MLV 4 First-Strand cDNA Synthesis Kit (Biomed, Beijing, China). The expression of miRNAs from the cDNA was detected using SYBR Green real-time fluorescence quantitative PCR. Stem-loop RT-qPCR is a common method used for detecting miRNA expression. It is characterized by the use of special primers with a stem structure to reverse transcribe RNA in the reverse transcription stage. Specific primers were designed to combine the products for subsequent fluorescence quantitative detection experiments. The reaction volume was 10 µL, which including 5 µL SYBR Green Master Mix (Bio-Rad),1 µL forward primer, 1µL reverse primer, and 100ng cDNA. The remained was supplemented to 10 µL with RNase-free water. The reaction was conducted on a LightCycler ®96 System (Roche, Germany). The programs for real time qPCR reactions were as follows: pre-incubation at 95°C for 300 s, 40 cycles of 95°C denaturation for 10s and 60°C annealing for 30s. Final melting curve analysis was performed to determine the purity of the PCR products. The specific primers for RT-PCR is shown in Table S1 and the Primers for qPCR is shown in Table S2.
Correlation analysis between gut microbiota and miRNAs
We quantified changes in community networks underlying microbial associations among the three groups using Netshift analysis, identifying microbes that may have a driving role for some of them.
To identify the crucial microbial species for high-altitude exposure and probiotic treatment, indicator species analysis was applied to calculate the indicator value of an Amplicon Sequence Variant (ASV) positive specificity and fidelity with one or more group, which was performed with 104 permutations and identified as significant at P < 0.05. Then, Sparcc correlation was used to calculate microbial genera in each group or ASV in all group associations, and a highly significant correlation (ρ > 0.8 and P < 0.01) was considered significant. The greedy optimization of modularity algorithm was applied to select and visualize the community modules in the co-occurrence network. Further, Netshift analysis was applied to explore the changes in the microbial network in the gut community among normal mice, hypoxia-exposed mice, and YH1136-treated mice, from which the genus was found to the main driving force. After these analyses, the crucial microbial genera and key ASVs were integrated to identify the species that could have the greatest effects on the obvious prevention of inflammatory injury induced by long-term hypoxic exposure by gut probiotics. Finally, Spearman’s correlation was used to calculate the degree of association (ρ > 0.8 and P < 0.01) between the selected ASVs and the miRNAs that were significantly different between groups and to visualize their networks and heatmaps.
Results
Lactobacillus johnsonii YH1136 alleviates intestinal barrier damage and inflammation induced by high-altitude hypoxic stress
The results of serum enzyme-linked immunosorbent assay (ELISA) showed that the levels of D-lactate and DAO in the HA group were significantly higher than those in groups C and HAP, but there was no significant difference between groups HAP and C (Figures 1A, B). Immunohistochemical analysis (Figure 1C) of the three tight junction proteins in the HA group decreased (the number of brown positive cells decreased), and supplementation with L. johnsonii YH1136 partially alleviated this decreasing trend. DAPI staining showed the expression profiles of IL-1β, TNF-α, and TGF-β, which are inflammatory factors in cells (Figure 2). These results suggest that the hypoxic environment at high altitudes may lead to intestinal inflammation and improve intestinal mucosal permeability or damage the intestinal barrier in mice, and that intragastric administration of L. johnsonii YH1136 could alleviate these effects to a certain extent.
Figure 1 Evaluation and analysis of intestinal permeability (A) Serum D-Lactic content. (B) Serum DAO content. Data are presented with the meas ± standard deviation (n=6). The marked * between the histograms indicates that the difference is significant. * means P<0.05, ** means P<0.01, *** means P<0.001. (C) Immunohistochemical analysis of the expression of mouse ileal tight junction proteins (Claudin-1, Occludin and ZO-1). The positive cells of the three proteins were stained brown. ns, not significant.
Figure 2 Analysis of inflammation factors in ileum DAPI was used to stain the nucleus, and different fluorescent genes were connected through secondary antibodies to show the expression of the target gene. The target genes were showed in red.
Lactobacillus johnsonii YH1136 modulates gut microbiota disturbance induced by high-altitude hypoxic stress
Figures 3A, B show the composition analysis of the intestinal microbiota at the phylum and genus levels, respectively. The results of the phylum level analysis (Figure 3A) showed that Firmicutes were most abundant in the three groups. The abundance of Firmicutes was the highest in group C and lower in the remaining two groups. There was a significant increase in the abundance of Proteobacteria in the HA group compared to that in group C, whereas the abundance of Proteobacteria in the HAP group reverted to that of group C. In addition, Bacteroidetes showed distinct differences among the three groups: a substantial decrease in abundance in the HA group and an increase in abundance in the HAP group compared to that in group C. These results suggest that high-altitude exposure may be beneficial to the colonization of Proteobacteria and adverse to colonization by Bacteroides, which was reversed by L.johnsonii YH1136 supplementation. In the analysis of species composition at the genus level (Figure 3B), changes in Staphylococcus spp. and Lactobacillus spp. were of the greatest concern. In the high-altitude exposure environment, Lactobacillus spp. almost completely disappeared in group HA, replaced by Staphylococcus spp., which is a common opportunistic pathogen in mammals (24), and its abundance sharply increased from less than 5% in group C to > 60% in group HA. In the HAP group, Lactobacillus abundance was effectively reversed after supplementation with L.johnsonii YH1136.
Figure 3 Effects of hypoxia environment and Lactobacillus johnsonii YH1136 supplementation on the composition, structure and diversity of mice ileal micriobiome (A) Phylum level species composition. (B) Genus level species composition. (C) Shannon diversity analysis. (D) Chao1 Index (E) PCoA (Bray_Curtis) analysis. (F) Indicator species analysis, the species sensitive to different treatments are divided into one plant different node colors represent different Phylum levels.
As shown in Figures 3C, D, the Shannon diversity and Chao1 index of groups HA and HAP were higher than those of the control group, and there was no significant difference between the other groups. Principal coordinate analysis (Figure 3E) showed that the distance between group HAP and group C was relatively close, indicating that the structure of intestinal microbiota in group HAP was the most similar to that in group C; the distance between groups HA and C was the farthest, indicating that the similarity of intestinal microbiota structure between group HA and group C was the worst. Indicator species analysis is usually used to explore species differences in intestinal microbiota between groups and the common and representative species of each group. Based on the results shown in Figure 3F, we observed that the important indicator species belonged to Firmicutes, Proteobacteria, Bacteroidetes, and Actinobacteria under the different treatments. In addition, there were four, fifteen, and five special indicator species in C, HA and HAP groups belonging to Firmicutes, Actinobacteria and Proteobacteria, respectively.
LEfSe is an algorithm for the discovery of biomarkers that identify genomic features (genes, pathways, or taxa) that characterize the differences between two or more biological conditions (25). Compared with the HAP group (Figures 4A, B), the significant discriminative taxa in group HA were Psychrobacter, Moraxellaceae, Pseudomonadales, Gammaproteobacteria, Staphylococcus, Staphylococcaceae, and Bacillus; and in the HAP group, Lactobacillaceae, Lactobacillus, Lactobacillales, Erysipelottrichales, Erysipelotrichaceae, Erysipelotrichia, and Dubosiella. Compared with the HA group (Figures 4C, D), the significant discriminative taxa in group C were Lactobacillaceae, Lactobacillus, Lactobacillales, Bacteroidales, Muribaculaceae, Bacteroidia, Bacteroidetes, uncultured-Bacteroidalesbacterium, Erysipelotrichales, Erysipelotrichaceae, Erysipelotrichia, and Dubosiella.
Figure 4 Lefse analysis of ileal microbiota in mice The linar discriminant analysis (LDA) score (A, C) and cladogram (B, D) were generated from LDA effect size. Taxa with LDA values larger than 4.5 are shown in figure.
Effects of high-altitude exposure and Lactobacillus johnsonii YH1136 supplementation on ileal miRNA expression
The sequencing of all groups yielded more than 15000000 raw reads on average, and the clean read ratio was > 80% in every sample (Table 1). We removesd the adapters and filtered the low-quality reads to obtain the clean reads. Clean reads with a sequence length of 18-36nt were counted, and the identical sequences in a single sample were merged to remove the repetitive reads, which are referred to as unique reads (Figure S1). There were 185830 unique reads shared by the three groups. Next, we counted the sequence abundance, and the results showed that the length of most unique reads ranged from 18 to 23nt (Figure S2). We compared unique reads with the reference genome to obtain the distribution of small RNA in the genome (Figure S3). The results showed that small RNAs were distributed in all chromosomes, mainly concentrated in the positive chain of the chromosome 1 reference genome and the negative chain of the chromosome 17 reference genome.
Table 1 The ileum small RNA sequencing reads of mice.
By comparing precursor and mature miRNAs of mice in the miRBase database, we annotated the unique reads, and obtained a statistical table of miRNA annotation results (Table 2). We then used mireap to perform a new miRNA prediction analysis on sequences that were not annotated; the results are shown in Table 3. In this way, we identified 147, 149, and 127 miRNAs in the HA, C and HAP groups, respectively. We also annotated the classification and proportion of all small RNAs (Figure 5A), analyzed the characteristics of all known miRNAs (Figure 5B), and performed principal component analysis (PCA) on all samples (Figure 5C). Before differential expression analysis, the correlation between miRNA expression levels among the samples was tested to ensure the reliability and rationality of the test (Figure 5D). According to the miRNA expression data in each sample, the differential expression of miRNA was analyzed by DEseq (version 1.18.0), and the different conserved miRNAs were screened according to the multiple differences in expression (|log2foldchange| > 1) and the significance of expression differences (p-value < 0.05). A volcano map of differentially expressed miRNAs was drawn in the R statistical environment using the ggplot2 package (Figures 5F–H). A total of 44 differentially expressed miRNAs were identified among the three groups. Cluster analysis was used to determine the expression patterns of differentially expressed miRNAs under different experimental conditions (Figure 5E). miRNAs with high expression correlations among samples were classified into one category.
Table 2 The ileum small RNA sequencing reads of mice.
Table 3 Prediction of novel miRNAs.
Figure 5 miRNA sequencing analysis (A) Unique reads of Small RNA classification statistics. (B) miRNA first base preference. (C) Principal component analysis of samples in each group. (D) Correlation analysis of miRNA expression levels among samples. (E) Cluster analysis of differentially expressed miRNAs. (F–H) The volcano plot of differentially expressed miRNAs.
Target gene prediction of differentially expressed miRNA and RT-qPCR verification
We used miRanda to predict the target genes of differentially expressed miRNA sequences with the 3 ‘UTR’ sequence of the mRNA of the species as the target sequence. TopGo was used for GO enrichment analysis to determine the main biological functions of the differentially expressed miRNA target genes. We used the Kyoto Encyclopedia of Gene and Genomes (KEGG) database to enrich and analyze the functions of the miRNA target genes and their related pathways. Figures 6 shows the top 20 GO terms (Figures 6A–C) and KEGG pathways (Figures 6D–F) with the lowest false discovery rate (FDR) value.
Figure 6 GO and KEGG Pathway enrichment analysis of target genes of differentially expressed miRNAs According to the enrichment results of Go, and the degree of enrichment is measured by rich factor, FDR value and the number of miRNA target genes enriched on this Go term and pathway. The larger the rich factor, the greater the degree of enrichment. The general value range of FDR is 0-1. The closer it is to zero, the more significant the enrichment is. (A–C) The GO analysis. (D–F) The KEGG Pathway enrichment analysis. The horizontal axis is rich factor, which refers to the ratio of the number of differential miRNA target genes enriched into GO term or pathway to the number of differential miRNA target genes annotated, with a larger value indicating more significant enrichment. The vertical axis is GO term or pathway.
GO analysis showed that the differentially expressed miRNA target genes in the HA and HAP groups were mainly enriched in system development, signal transduction, and multicellular organism development, while the differentially expressed miRNA target genes in groups C and HA were mainly enriched in signal regulation, cellular communication regulation, and biomass regulation. The differentially expressed miRNAs target gene in the groups C and HAP weere mainly enriched in signaling receiver activity, cellular protein modification processes, and protein modification processes.
In addition, KEGG pathway analysis showed that the significantly enriched pathways of the target genes of differentially expressed miRNAs were the Wnt signaling pathway, miRNAs in cancer, parathyroid hormone synthesis, secretion and action, lysine degradation, phosphatidylinositol signaling system, inositol phosphate metabolism, lysine degradation, thyroid hormone signaling pathway, and the mTOR signaling pathway.
We used RT-qPCR to analyze miRNA expression and compared it with the sequencing results to determine the accuracy of the results. Four miRNAs with significant differential expression were randomly selected from HA vs. HAP and HA vs. C for quantitative fluorescence analysis. Primer design is presented in Tables S1, S2, and the results are shown in Figure 7, with U6 snRNA as the housekeeping gene. Fold change was calculated using the △△Ct method and was compared with the results of small RNA sequencing. The results showed that the fold-change trends of the DE miRNAs obtained by the two methods were consistent between groups, and the values were also relatively close, indicating that the sequencing results were reliable.
Figure 7 Validation of differentilly expressed of miRNA. RT-qPCR validation experiment of randomly selected differentially expressed miRNAs. Four miRNAs including mmu-miR-3743a, mmu-miR-20a-3p, mmu-miR-122-5p and mmu-miR-3060-3p were selected from HAP group and HA group respectively, and four miRNAs including mmu-miR-3102-3p, mmu-miR-802-3p, mmu-miR-34b-5p and mmu-miR-5114 were selected from HA group and C group to complete the validation test of RT-qPCR.
The correlation analysis between miRNAs and gut microbiota
In the results of network properties, we noticed an evident decrease in the sub-network average path length, number of nodes, edges, and exclusive edges, as well as an obvious increase in the sub-network density of the microbial community of HA mice compared with the control mice (Figures 8A–E). Similarly, with the intake of L.johnsonii YH1136, the HAP group reduced the total number of sub-networks of nodes, edges, and exclusive edges; instead, the density and, average path length of the microbial community increased (Figures 8F–J).
Figure 8 Netshift and co-occurrence network analysis (A-E) Changes of network properties between group C and group HA. (F-J) Changes of network properties between group HA and group HAP. (K) Common sub-networks analysis of group C and group HA. (L) Common sub-networks analysis of group HA and group HAP. (M) Each node in the co-occurrence network represented an ASV closely related species, which was grouped into the same module, and different colors were used to highlight the response pattern of the ASV module.
To explore the driving microbes associated with network changes, we focused on the importance of different microbial genera using the DelBet index and NESH scores. Therefore, nodes (genus) with a positive DelBet index and higher NESH scores in the network shift were observed. The “driver microbes” (top five) that altered the symbiotic network of Genus in the group HA compared to the group C (Table S3) are Globicatella, Staphylococcus, Agathobacter, Alloprevotella and Aminobacterium. Moreover, compared with group HA (Table S4), the “driving microbes” (top five) that changed the genus associtation network in the HAP group were Corynebacterium 1, Desulfovibrio, Bifidobacterium, Pseudochrobactrum, and Prevotella 9.
By integrating the results of the crucial microbial genus and key ASV (Table S5), Staphylococcus (ASV14, ASV16 and ASV20) and Corynebacterium 1 (ASV73, ASV60, ASV61 and ASV80) were identified as the driving microbes.
To identify the core modules of the microbial community, we identified five important communities in the co-occurrence network, including indicator species, which corresponded to the response patterns in the species community under hypoxic exposure and L. johnsonii YH1136 treatment (Figure 8M). We noticed that modules 1 and 2 communities need to be focused on, as most of the species in these module communities had an obvious response to altitude hypoxic environment, but this response disappeared after the supplementation with L. johnsonii YH1136. Therefore, the indicator species classified as HA responsive are the key species that we will focus on in further analysis.
As shown in Figure 9, we analyzed the association between differentially expressed miRNAs and intestinal microbiota using the Spearman correlation analysis, and performed a Netshift assay. The results of the Netshift analysis revealed that Staphylococcus, which had higher NESH scores in the plateau environment, was a notable gut microbiota driving-species. Additionally, it showed a strong negative correlation with Lactobacillus. In the correlation analysis between miRNAs and gut microbiota, miR-196a-1-3p in the HA group had a strong negative correlation with Staphylococcus, while miR-3060-3p in the HAP group had a strong positive correlation with Corynebacterium1. TargetScan and mirpath were used (26) to predict the target genes and signaling pathways of these two miRNAs (Table S6). The results showed that miR-196a-1-3p was enriched in ECM-receptor interaction and miR-3060-3p was enriched in the Hippo signaling pathway, which might be the regulatory pathways of high-altitude hypoxic environments and probiotics.
Figure 9 Correlation analysis between miRNAs and intestinal microbiota (A, D) Correlations between miRNA in C group and intestinal microbiota. (B, E) Correlations between miRNA in HA groups and intestinal microbiota. (C, F) Correlations between miRNA in HAP groups and intestinal microbiota. (G) Correlation analysis of driving microbes (Staphylococcus and Corynebacterium 1). The colors range form blue to red corresponds to negative correlation and positive correlation, respectively. The ‘*’ means P value < 0.05, ‘**’ means P value < 0.01.
Discussion
High-altitude areas have lower air pressure and oxygen levels compared to low-altitude areas, which has an impact on the physiological functions of animals and humans, particularly on intestinal microbiota and function (2, 3, 27), inducing gastrointestinal diseases and intestinal microbiota disorders (6, 28, 29). In the past decades of research, changes in the gut microbiota have been increasingly associated with the development of a variety of gastro-intestinal diseases such as inflammatory diseases (8, 30–32). However, the changes in the intestinal microbiota and the regulatory mechanism of probiotics in high-saltitude exposure environment are still unclear. Through 16S high-throughput sequencing analysis, we can understand the changes in intestinal microbiota exposed to high altitude environments, which is helpful in clarifying the resilience of animals and humans to high-altitude exposure environments. In our study, under the simulated hypoxic environments, damage to the intestinal mucosal barrier occurred, with increased permeability and reduced tight junction protein structure. The results of DAPI staining showed that the ileum exhibited inflammatory symptoms in high-altitude exposure environments, and that supplementation with probiotic L.johnsonii YH1136 by gavage could alleviate these asymptoms. D-lactate, a metabolite of intestinal bacteria, and DAO produced by intestinal mucosal cells can enter the blood when the intestinal barrier is injured, so the detection of these two indicators in serum can reflect damage to the intestinal barrier (33, 34). In our experiment, probiotic supplementation helped reduce the increased serum D-lactate and DAO content caused by high-altitude modeling, indirectly demonstrating that probiotics played a protective role. This is similar to the results of Khanna et al. (29), on acute hypoxic exposure in rats, which also found varying degrees of intestinal mucosal damage in the form of increased mucosal permeability, and impaired intestinal villus structure. In another study, symbiotic supplementation alleviated intestinal barrier damage induced by hypobaric hypoxic exposure and reduced intestinal inflammation (35).
We found that the composition of the intestinal microbiota in the ileum of mice changed significantly, with Lactobacillaceae predominating in group C, whereas Staphylococcaceae predominated in the HA group. Meanwhile, the intestinal microbiota of the group HAP gavaged with L.johnsonii YH1136 showed a composition similar to that of group C. Staphylococcus aureus belongigns to the genus Staphylococcus is mostly the pathogenic bacteria, which has been found to be closely related to intestinal infection in lot of studies (36, 37), can cause inflammatory diseases of the intestine, and is a common potentially pathogenic bacterium of the intestine (38, 39). Lactobacillus salivarius (40), Lactobacillus fermentum (40), Lactobacillus rhamnosus (36), and Lactobacillus kiferi (41) have been confirmed in several studies to have direct or indirect antibacterial activity against Staphylococcus aureus. In our study, the abundance of Staphylococcus was significantly increased by simulated high-altitude, which suggesting that the number of pathogenic bacteria in Staphylococcus may also be significantly higher, in other words, a high-altitude hypoxic environment may cause higher numbers of pathogenic bacteria to become the dominant species. But this hypothesis remains to be determined by subsequent analysis. A series of 16S sequencing analysis results revealed that supplementation with L. johnsonii YH1136 reversed this trend, returning the gut microbiota at the ileum genus level to the normal Lactobacillus-dominated structure.
With the in-depth study of probiotics, there have been several reports on the prevention and treatment of intestinal structural and functional injury (42–44). Small RNA sequencing technology has been continuously developing, and miRNAs are increasingly used in the study of various molecular mechanisms (45, 46). miRNAs are small non-coding RNA, with a length is usually 19-25 nucleotides, which can regulate the post-transcriptional silencing of target genes (47). miRNAs can be involved in regulating a variety of cell activities, including cell differentiation, growth, development, and apoptosis (48). In addition, it has also been demonstrated that hypoxia can alter miRNA expression. Kulshreshtha et al. (49) demonstrated that a specific spectrum of miRNAs (including miR-23, miR-24, miR-26, miR-27, miR-103, miR-107, miR-181, miR-210, and miR-213) was induced in response to hypoxia, some of which were induced through a hypoxia-inducible-factor (HIF)-dependent mechanism. Bhandari et al. (50) observed hypoxia related to miRNA imbalance in cancer and verified its role in hypoxic regulation. Li et al. (51) showed that miR-21-5p has a protective effect by promoting bone matrix differentiation under hypoxic conditions. Probiotics have also been shown to regulate miRNA expression (52–54).
Muenchau et al. (55), analyzed the miRNA expression profile in human intestinal epithelial cells (IECs) in a hypoxic environment and found that miRNA-320a has a direct role in regulating IEC barrier function, suggesting proposing that a hypoxic environment in the intestinal lumen may not only mediate the regulation of tight junctions and adhesion protein expression through HIF, but also affect intestinal barrier function through miRNA-based regulation of cell-cell formation. In this study, by analyzing the sequencing results of miRNAs, we identified more than 12 600 miRNAs from all samples and 161 novo miRNAs in all groups. Subsequently, differentially expressed miRNAs in the three groups were identified. A total of 44 differentially expressed miRNAs (DEmiRNAs) were obtained, of which 13 were between HA and HAP groups, 8 were up-expressed and 5 were down-regulated. There were 22 up-regulated and 11 down-regulated DEmiRNAs between groups C and HA. There were 9 DE miRNAs between groups C and HAP, with four up-regulated and five down-regulated. These DEmiRNAs may play a role in the regulation of intestinal structure and function under hypoxia induced by L. johnsonii YH1136. Shi et al. (56). showed that hypoxia could inhibit the expression of miR-499-5p, and that overexpression of miR-499-5p could inhibit apoptosis induced by hypoxia/reoxygenation by targeting SOX6. In this study, the sequencing results showed that mmu-miR-499a-5p was significantly down-regulated in group HA compared to group C, which was consistent with the results of Shi et al. miR-499-5p was significantly up-regulated in the HAP group compared with the HA group, suggesting that L. johnsonii YH1136 may exert beneficial effects by reducing apoptosis through miR-499-5p-SOX6. MiR-122-5p was also significantly differential expressed between the HAP and HA groups, with higher expression in both groups (compared with other low expression miRNAs), which was significantly down regulated in the HAP group. Hou et al. (57) analyzed the miRNA expression profile of plasma exosomes and reported that some miRNAs, including miR-122-5p, were positively correlated with disease activity in intestinal Behçet’s syndrome.
miR-196a-1-3p and miR-3060-3p were considered important regulatory miRNAs through correlation analysis, and their target genes and potentially enriched pathways were further predicted. miR-196a-1-3p is enriched in ECM-receptor interactions, while the structural components of the ECM actively participate in intestinal inflammation (58). miR-3060-3p is enriched in the Hippo signaling pathway. A study by Yu et al. (59) showed that miR-590-5p could alleviate intestinal inflammation by targeting YAP, a component of the Hippo signaling pathway. This may explain the regulation of intestinal inflammation and injury by miR-3060-3p in response to high-altitude conditions. However, whether miR-3060-3p has similar effects requires further exploration and verification.
We further studied the role of L. johnsonii YH1136 supplementation in regulating miRNAs involved in the processes of intestinal injury caused by a hypoxic environment and used miRanda to predict target genes. Among them, we focused on the target genes of miRNAs differentially expresseed between groups HA and HAP. These target genes are mainly invloved in system development, signal transduction, and multicellular organism development. The Wnt signaling pathway was significantly enriched in all groups. Genetic studies have shown that the Wnt/β- catenin signaling pathway plays an important role in the proliferation of intestinal epithelial cells (60). Li et al. (61) found that the miR-155/HBP1 axis activates Wnt/β-catenin signaling pathway and induces intestinal fibrosis. In our microbiota-miRNA correlation analysis, Lactobacillus showed a strong positive or negative correlation with a variety of miRNAs. Therefore, we speculated that L. johnsonii YH1136 significantly affects the intestinal barrier by up or down regulating various miRNAs. MicroRNAs have been confirmed to remodel the composition of intestinal microbiota in several studies (62, 63), while probiotics can regulate the expression of certain miRNAs. This suggests that probiotics may not only directly affect the composition of the intestinal microbiota but may also remodel it by affecting specific miRNAs expression. In addition to the miR-196a-1-3p and miR-3060-3p, miR-485-5p, miR-7219-3p, and miR-103-2-5p were also of interest, and their target gene were predicted using the miRDB database. After screening, we found that a target gene of miR-485-5p, Pak1, the p21 activated kinase, is associated with multiple diseases (64, 65) and may play a regulatory role in intestinal inflammation induced by a high-altitude hypoxic environment. Overexpression studies of PAK1, which has been shown to have higher expression in inflammatory bowel disease, showed higher cell proliferation and reduced apoptosis when overexpressed (66). In addition, the target gene of miR-7219-3p, Adamdec1, which protects the intestine from chemical and bacterial insults that may contribute to the development of Crohn’s disease (67), was also confirmed to be downregulated by probiotics in one study (68). Likewise, Cyld, a target gene of miR-103-2-5p, plays an important regulatory role in intestinal inflammation (69). More in-depth studies are required to determine pathways that plays key regulatory roles in intestinal diseases. Rodríguez-Nogales et al. (54), pointed out that the probiotic Saccharomyces boulardii could play a protective role in DSS-induced colitis in mice by affecting the expression of miRNAs and the composition of the intestinal microbiota. This is similar to the protective effect of L. johnsonii YH1136 on the intestinal barrier in hypoxic environments.
Conclusion
High-altitude exposure can lead to disorders in the intestinal microbiota and promote intestinal barrier impairment in mice, triggering intestinal inflammation, and endogenous pathogenic bacteria induced intestinal dysfunction. Supplementation with L. johnsonii YH1136 via gavage could be effective in alleviating these symptoms by reshaping the intestinal microbiota, reducing the abundance of pathogenic bacteria and regulating the expression of target genes in conjunction with miRNAs in the ileum. Though a number of positive results of L. johnsonii YH1136 were found in the present experiment, more studies in mice exposed at different altitudes should be undertaken to fully support the protective effect of such a potential probiotic strain against high altitude exposure. The specific regulatory mechanisms of miRNAs in the process of high-altitude exposure leading to intestinal barrier damage should also be demonstrated.
Data availability statement
The data presented in the study are deposited in the NCBI SRA database, accession number PRJNA865817 (MicroRNA sequencing raw data) and PRJNA866881(16S rRNA amplicon sequencing raw data).
Ethics statement
The animal study was reviewed and approved by the Committee of Sichuan Agriculture University. Animal experiments were conducted in accordance with the guidelines for the feeding and use of experiment animals (approval No: SYXKchuan2019-187) approved by Committee of Sichuan agriculture University.
Author contributions
ZW, XZ and XJ conceived the research. ZW, XZ, NS, XJ and JX performed the experiments. ZW, NS and XZ wrote the manuscript. XJ, YQ and NS analyzed the experimental data. All authors contributed to the article and approved the submitted version.
Funding
This work is supported by the Sichuan Science and Technology Program (2021YJ0166), Scientific Development funds for Local Region from the Chinese Government in 2022 (XZ202201YD0018C), Reformation and Development Funds for Local Region Universities from the Chinese Government in 2021 (zf21003002) and Guangdong Basic and Applied Basic Research Foundation (Grant No: 2020A1515110693).
Conflict of interest
HW was employed by the Guangzhou Beneco Biotechnology Co., Ltd.
The remaining 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/fimmu.2022.1007737/full#supplementary-material
Abbreviations
L., Lactobacillus; ASV, Amplicon Sequence Variants; IBD, Inflammation Bowel Disease; miRNA, MicroRNA; DEPC, Diethyl pyrocarbonate; DAO, Diamine oxidase; ANOVA, Analysis of Variance; LSD, Least Significant Difference; PBS, Phosphate Buffered Saline; BSA, Bovine Serum Albumin; HRP, Horseradish Peroxidase; DAB, Diaminobenzidine; DAPI, 4’,6-diamidino-2-phenylindole; PCR, Polymerase Chain Reaction; RT-qPCR, Real Time Quantitative PCR; ASV, Amplicon Sequence Variants; TNF-α, Tumor Necrosis Factor-α; TGF-β, Transforming Growth Factor-β; IL-1β (Interleukin-1β).
References
1. Quigley EMM. Prebiotics and probiotics in digestive health. Clin Gastroenterol Hepatol (2019) 17(2):333–44. doi: 10.1016/j.cgh.2018.09.028
2. Jia Z, Zhao X, Liu X, Zhao L, Jia Q, Shi J, et al. Impacts of the plateau environment on the gut microbiota and blood clinical indexes in han and Tibetan individuals. mSystems (2020) 5(1):e00660-19. doi: 10.1128/mSystems.00660-19
3. Wen J, Zhang JQ, Nie ZL, Zhong Y, Sun H. Evolutionary diversifications of plants on the qinghai-Tibetan plateau. Front Genet (2014) 5:4. doi: 10.3389/fgene.2014.00004
4. Sánchez de Medina F, Romero-Calvo I, Mascaraque C, Martínez-Augustin O. Intestinal inflammation and mucosal barrier function. Inflammation Bowel Dis (2014) 20(12):2394–404. doi: 10.1097/mib.0000000000000204
5. Al-Sadi R, Nighot P, Nighot M, Haque M, Rawat M, Ma TY. Lactobacillus acidophilus induces a strain-specific and toll-like receptor 2-dependent enhancement of intestinal epithelial tight junction barrier and protection against intestinal inflammation. Am J Pathol (2021) 191(5):872–84. doi: 10.1016/j.ajpath.2021.02.003
6. Wang Y, Shi Y, Li W, Wang S, Zheng J, Xu G, et al. Gut microbiota imbalance mediates intestinal barrier damage in high-altitude exposed mice. FEBS J (2022) 289(16):4850–68. doi: 10.1111/febs.16409
7. McKenna ZJ, Gorini Pereira F, Gillum TL, Amorim FT, Deyhle MR, Mermier CM. High-altitude exposures and intestinal barrier dysfunction. Am J Physiol Regul Integr Comp Physiol (2022) 322(3):R192–r203. doi: 10.1152/ajpregu.00270.2021
8. Ananthakrishnan AN, Bernstein CN, Iliopoulos D, Macpherson A, Neurath MF, Ali RAR, et al. Environmental triggers in IBD: a review of progress and evidence. Nat Rev Gastroenterol Hepatol (2018) 15(1):39–49. doi: 10.1038/nrgastro.2017.136
9. Pan Z, Hu Y, Huang Z, Han N, Li Y, Zhuang X, et al. Alterations in gut microbiota and metabolites associated with altitude-induced cardiac hypertrophy in rats during hypobaric hypoxia challenge. Sci China Life Sci (2022), 1–21. doi: 10.1007/s11427-021-2056-1
10. Wang H, He S, Xin J, Zhang T, Sun N, Li L, et al. Psychoactive effects of lactobacillus johnsonii against restraint stress-induced memory dysfunction in mice through modulating intestinal inflammation and permeability-a study based on the gut-brain axis hypothesis. Front Pharmacol (2021) 12:662148. doi: 10.3389/fphar.2021.662148
11. Cristofori F, Dargenio VN, Dargenio C, Miniello VL, Barone M, Francavilla R. Anti-inflammatory and immunomodulatory effects of probiotics in gut inflammation: A door to the body. Front Immunol (2021) 12:578386. doi: 10.3389/fimmu.2021.578386
12. Ma Y, Ma S, Chang L, Wang H, Ga Q, Ma L, et al. Gut microbiota adaptation to high altitude in indigenous animals. Biochem Biophys Res Commun (2019) 516(1):120–6. doi: 10.1016/j.bbrc.2019.05.085
13. Zhang Z, Xu D, Wang L, Hao J, Wang J, Zhou X, et al. Convergent evolution of rumen microbiomes in high-altitude mammals. Curr Biol (2016) 26(14):1873–9. doi: 10.1016/j.cub.2016.05.012
14. Wang H, Sun Y, Xin J, Zhang T, Sun N, Ni X, et al. Lactobacillus johnsonii BS15 prevents psychological stress-induced memory dysfunction in mice by modulating the gut-brain axis. Front Microbiol (2020) 11:1941. doi: 10.3389/fmicb.2020.01941
15. Xin J, Zeng D, Wang H, Ni X, Yi D, Pan K, et al. Preventing non-alcoholic fatty liver disease through lactobacillus johnsonii BS15 by attenuating inflammation and mitochondrial injury and improving gut environment in obese mice. Appl Microbiol Biotechnol (2014) 98(15):6817–29. doi: 10.1007/s00253-014-5752-1
16. Xin J, Sun N, Wang H, Ma H, Wu B, Li L, et al. Preventive effects of lactobacillus johnsonii on the renal injury of mice induced by high fluoride exposure: Insights from colonic microbiota and co-occurrence network analysis. Ecotoxicol Environ Saf (2021) 228:113006. doi: 10.1016/j.ecoenv.2021.113006
17. Sun N, Zhu B, Xin J, Li L, Gan B, Cao X, et al. Psychoactive effects of lactobacillus johnsonii BS15 on preventing memory dysfunction induced by acute ethanol exposure through modulating intestinal microenvironment and improving alcohol metabolic level. Front Microbiol (2022) 13:847468. doi: 10.3389/fmicb.2022.847468
18. Sun N, Ni X, Wang H, Xin J, Zhao Y, Pan K, et al. Probiotic lactobacillus johnsonii BS15 prevents memory dysfunction induced by chronic high-fluorine intake through modulating intestinal environment and improving gut development. Probiotics Antimicrob Proteins (2020) 12(4):1420–38. doi: 10.1007/s12602-020-09644-9
19. Xin J, Wang H, Sun N, Bughio S, Zeng D, Li L, et al. Probiotic alleviate fluoride-induced memory impairment by reconstructing gut microbiota in mice. Ecotoxicol Environ Saf (2021) 215:112108. doi: 10.1016/j.ecoenv.2021.112108
20. Liu L, Yao J, Li Z, Zu G, Feng D, Li Y, et al. miR-381-3p knockdown improves intestinal epithelial proliferation and barrier function after intestinal ischemia/reperfusion injury by targeting nurr1. Cell Death Dis (2018) 9(3):411. doi: 10.1038/s41419-018-0450-z
21. Al-Sadi R, Engers J, Abdulqadir R. Talk about micromanaging! role of microRNAs in intestinal barrier function. Am J Physiol Gastrointest Liver Physiol (2020) 319(2):G170–g174. doi: 10.1152/ajpgi.00214.2020
22. Ikemura K, Iwamoto T, Okuda M. MicroRNAs as regulators of drug transporters, drug-metabolizing enzymes, and tight junctions: implication for intestinal barrier function. Pharmacol Ther (2014) 143(2):217–24. doi: 10.1016/j.pharmthera.2014.03.002
23. Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, et al. Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res (2005) 33(20):e179. doi: 10.1093/nar/gni178
24. Lindsay JA. Staphylococci: Evolving genomes. Microbiol Spectr (2019) 7(6):485–98. doi: 10.1128/microbiolspec.GPP3-0071-2019
25. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol (2011) 12(6):R60. doi: 10.1186/gb-2011-12-6-r60
26. Vlachos IS, Zagganas K, Paraskevopoulou MD, Georgakilas G, Karagkouni D, Vergoulis T, et al. DIANA-miRPath v3.0: deciphering microRNA function with experimental support. Nucleic Acids Res (2015) 43(W1):W460–6. doi: 10.1093/nar/gkv403
27. Tian L, Jia Z, Xu Z, Shi J, Zhao X, He K. Transcriptional landscape in rat intestines under hypobaric hypoxia. PeerJ (2021) 9:e11823. doi: 10.7717/peerj.11823
28. Ji Q, Zhang Y, Zhou Y, Gamah M, Yuan Z, Liu J, et al. Effects of hypoxic exposure on immune responses of intestinal mucosa to citrobacter colitis in mice. BioMed Pharmacother (2020) 129:110477. doi: 10.1016/j.biopha.2020.110477
29. Khanna K, Mishra KP, Chanda S, Eslavath MR, Ganju L, Kumar B, et al. Effects of acute exposure to hypobaric hypoxia on mucosal barrier injury and the gastrointestinal immune axis in rats. High Alt Med Biol (2019) 20(1):35–44. doi: 10.1089/ham.2018.0031
30. Lee M, Chang EB. Inflammatory bowel diseases (IBD) and the microbiome-searching the crime scene for clues. Gastroenterology (2021) 160(2):524–37. doi: 10.1053/j.gastro.2020.09.056
31. Harris KG, Chang EB. The intestinal microbiota in the pathogenesis of inflammatory bowel diseases: new insights into complex disease. Clin Sci (Lond) (2018) 132(18):2013–28. doi: 10.1042/cs20171110
32. Chassaing B, Van de Wiele T, De Bodt J, Marzorati M, Gewirtz AT. Dietary emulsifiers directly alter human microbiota composition and gene expression ex vivo potentiating intestinal inflammation. Gut (2017) 66(8):1414–27. doi: 10.1136/gutjnl-2016-313099
33. Song WB, Lv YH, Zhang ZS, Li YN, Xiao LP, Yu XP, et al. Soluble intercellular adhesion molecule-1, d-lactate and diamine oxidase in patients with inflammatory bowel disease. World J Gastroenterol (2009) 15(31):3916–9. doi: 10.3748/wjg.15.3916
34. Cai J, Chen H, Weng M, Jiang S, Gao J. Diagnostic and clinical significance of serum levels of d-lactate and diamine oxidase in patients with crohn's disease. Gastroenterol Res Pract (2019) 2019:8536952. doi: 10.1155/2019/8536952
35. Khanna K, Mishra KP, Chanda S, Ganju L, Singh SB, Kumar B. Effect of synbiotics on amelioration of intestinal inflammation under hypobaric hypoxia. High Alt Med Biol (2021) 22(1):32–44. doi: 10.1089/ham.2020.0062
36. Liu G, Pang B, Li N, Jin H, Li J, Wu W, et al. Therapeutic effect of lactobacillus rhamnosus SHA113 on intestinal infection by multi-drug-resistant staphylococcus aureus and its underlying mechanisms. Food Funct (2020) 11(7):6226–39. doi: 10.1039/d0fo00969e
37. Larcombe S, Jiang JH, Hutton ML, Abud HE, Peleg AY, Lyras D. A mouse model of staphylococcus aureus small intestinal infection. J Med Microbiol (2020) 69(2):290–7. doi: 10.1099/jmm.0.001163
38. Kang SS, Noh SY, Park OJ, Yun CH, Han SH. Staphylococcus aureus induces IL-8 expression through its lipoproteins in the human intestinal epithelial cell, caco-2. Cytokine (2015) 75(1):174–80. doi: 10.1016/j.cyto.2015.04.017
39. Azimirad M, Dezfulian A, Alebouyeh M, Bahreiny Esfehani R, Shahrokh S, Zali MR. Infection with enterotoxigenic staphylococcus aureus as a concern in patients with gastroenteritis. J Glob Antimicrob Resist (2017) 9:111–4. doi: 10.1016/j.jgar.2017.02.009
40. Kang MS, Lim HS, Oh JS, Lim YJ, Wuertz-Kozak K, Harro JM, et al. Antimicrobial activity of lactobacillus salivarius and lactobacillus fermentum against staphylococcus aureus. Pathog Dis (2017) 75(2). doi: 10.1093/femspd/ftx009
41. Xiu L, Sheng S, Hu Z, Liu Y, Li J, Zhang H, et al. Exopolysaccharides from lactobacillus kiferi as adjuvant enhanced the immuno-protective against staphylococcus aureus infection. Int J Biol Macromol (2020) 161:10–23. doi: 10.1016/j.ijbiomac.2020.06.005
42. Riehl TE, Alvarado D, Ee X, Zuckerman A, Foster L, Kapoor V, et al. Lactobacillus rhamnosus GG protects the intestinal epithelium from radiation injury through release of lipoteichoic acid, macrophage activation and the migration of mesenchymal stem cells. Gut (2019) 68(6):1003–13. doi: 10.1136/gutjnl-2018-316226
43. Wu H, Xie S, Miao J, Li Y, Wang Z, Wang M, et al. Lactobacillus reuteri maintains intestinal epithelial regeneration and repairs damaged intestinal mucosa. Gut Microbes (2020) 11(4):997–1014. doi: 10.1080/19490976.2020.1734423
44. Wu Y, Wang B, Xu H, Tang L, Li Y, Gong L, et al. Probiotic bacillus attenuates oxidative stress- induced intestinal injury via p38-mediated autophagy. Front Microbiol (2019) 10:2185. doi: 10.3389/fmicb.2019.02185
45. Shan TD, Ouyang H, Yu T, Li JY, Huang CZ, Yang HS, et al. miRNA-30e regulates abnormal differentiation of small intestinal epithelial cells in diabetic mice by downregulating Dll4 expression. Cell Prolif (2016) 49(1):102–14. doi: 10.1111/cpr.12230
46. Li Y, Zhu L, Chen P, Wang Y, Yang G, Zhou G, et al. MALAT1 maintains the intestinal mucosal homeostasis in crohn's disease via the miR-146b-5p-CLDN11/NUMB pathway. J Crohns Colitis (2021) 15(9):1542–57. doi: 10.1093/ecco-jcc/jjab040
47. Lu TX, Rothenberg ME. MicroRNA. J Allergy Clin Immunol (2018) 141(4):1202–7. doi: 10.1016/j.jaci.2017.08.034
48. Saliminejad K, Khorram Khorshid HR, Soleymani Fard S, Ghaffari SH. An overview of microRNAs: Biology, functions, therapeutics, and analysis methods. J Cell Physiol (2019) 234(5):5451–65. doi: 10.1002/jcp.27486
49. Kulshreshtha R, Ferracin M, Wojcik SE, Garzon R, Alder H, Agosto-Perez FJ, et al. A microRNA signature of hypoxia. Mol Cell Biol (2007) 27(5):1859–67. doi: 10.1128/mcb.01395-06
50. Bhandari V, Hoey C, Liu LY, Lalonde E, Ray J, Livingstone J, et al. Molecular landmarks of tumor hypoxia across cancer types. Nat Genet (2019) 51(2):308–18. doi: 10.1038/s41588-018-0318-2
51. Li L, Jiang D. Hypoxia-responsive miRNA-21-5p inhibits Runx2 suppression by targeting SMAD7 in MC3T3-E1 cells. J Cell Biochem (2019) 120(10):16867–75. doi: 10.1002/jcb.28944
52. Davoodvandi A, Marzban H, Goleij P, Sahebkar A, Morshedi K, Rezaei S, et al. Effects of therapeutic probiotics on modulation of microRNAs. Cell Commun Signal (2021) 19(1):4. doi: 10.1186/s12964-020-00668-w
53. Zhao Y, Zeng Y, Zeng D, Wang H, Zhou M, Sun N, et al. Probiotics and MicroRNA: Their roles in the host-microbe interactions. Front Microbiol (2020) 11:604462. doi: 10.3389/fmicb.2020.604462
54. Rodríguez-Nogales A, Algieri F, Garrido-Mesa J, Vezza T, Utrilla MP, Chueca N, et al. Intestinal anti-inflammatory effect of the probiotic saccharomyces boulardii in DSS-induced colitis in mice: Impact on microRNAs expression and gut microbiota composition. J Nutr Biochem (2018) 61:129–39. doi: 10.1016/j.jnutbio.2018.08.005
55. Muenchau S, Deutsch R, de Castro IJ, Hielscher T, Heber N, Niesler B, et al. Hypoxic environment promotes barrier formation in human intestinal epithelial cells through regulation of MicroRNA 320a expression. Mol Cell Biol (2019) 39(14):e00553-18. doi: 10.1128/mcb.00553-18
56. Shi Y, Han Y, Niu L, Li J, Chen Y. MiR-499 inhibited hypoxia/reoxygenation induced cardiomyocytes injury by targeting SOX6. Biotechnol Lett (2019) 41(6-7):837–47. doi: 10.1007/s10529-019-02685-3
57. Hou CC, Bao HF, Shen Y, Ye JF, Ma HF, Guan JL. Expression of miRNAs derived from plasma exosomes in patients with intestinal behçet's syndrome. Clin Exp Rheumatol (2021). doi: 10.55563/clinexprheumatol/6xgxzk
58. Shimshoni E, Yablecovitch D, Baram L, Dotan I, Sagi I. ECM remodelling in IBD: innocent bystander or partner in crime? the emerging role of extracellular molecular events in sustaining intestinal inflammation. Gut (2015) 64(3):367–72. doi: 10.1136/gutjnl-2014-308048
59. Yu M, Luo Y, Cong Z, Mu Y, Qiu Y, Zhong M. MicroRNA-590-5p inhibits intestinal inflammation by targeting YAP. J Crohns Colitis (2018) 12(8):993–1004. doi: 10.1093/ecco-jcc/jjy046
60. Richmond CA, Breault DT. Regulation of gene expression in the intestinal epithelium. Prog Mol Biol Transl Sci (2010) 96:207–29. doi: 10.1016/b978-0-12-381280-3.00009-9
61. Li N, Ouyang Y, Xu X, Yuan Z, Liu C, Zhu Z. MiR-155 promotes colitis-associated intestinal fibrosis by targeting HBP1/Wnt/β-catenin signalling pathway. J Cell Mol Med (2021) 25(10):4765–75. doi: 10.1111/jcmm.16445
62. Liu S, da Cunha AP, Rezende RM, Cialic R, Wei Z, Bry L, et al. The host shapes the gut microbiota via fecal MicroRNA. Cell Host Microbe (2016) 19(1):32–43. doi: 10.1016/j.chom.2015.12.005
63. Teng Y, Ren Y, Sayed M, Hu X, Lei C, Kumar A, et al. Plant-derived exosomal MicroRNAs shape the gut microbiota. Cell Host Microbe (2018) 24(5):637–52. doi: 10.1016/j.chom.2018.10.001
64. Cheng WL, Zhang Q, Li B, Cao JL, Jiao L, Chao SP, et al. PAK1 silencing attenuated proinflammatory macrophage activation and foam cell formation by increasing PPARγ expression. Oxid Med Cell Longev (2021) 2021:6957900. doi: 10.1155/2021/6957900
65. Dammann K, Khare V, Harpain F, Lang M, Kurtovic A, Mesteri I, et al. PAK1 promotes intestinal tumor initiation. Cancer Prev Res (Phila) (2015) 8(11):1093–101. doi: 10.1158/1940-6207.Capr-15-0205-t
66. Khare V, Dammann K, Asboth M, Krnjic A, Jambrich M, Gasche C. Overexpression of PAK1 promotes cell survival in inflammatory bowel diseases and colitis-associated cancer. Inflammation Bowel Dis (2015) 21(2):287–96. doi: 10.1097/mib.0000000000000281
67. O'Shea NR, Chew TS, Dunne J, Marnane R, Nedjat-Shokouhi B, Smith PJ, et al. Critical role of the disintegrin metalloprotease ADAM-like decysin-1 [ADAMDEC1] for intestinal immunity and inflammation. J Crohns Colitis (2016) 10(12):1417–27. doi: 10.1093/ecco-jcc/jjw111
68. Plaza-Díaz J, Robles-Sánchez C, Abadía-Molina F, Morón-Calvente V, Sáez-Lara MJ, Ruiz-Bravo A, et al. Adamdec1, ednrb and Ptgs1/Cox1, inflammation genes upregulated in the intestinal mucosa of obese rats, are downregulated by three probiotic strains. Sci Rep (2017) 7(1):1939. doi: 10.1038/s41598-017-02203-3
Keywords: high-altitude exposure, probiotic, Lactobacillus johnsonii, intestinal microbiota, miRNA
Citation: Wan Z, Zhang X, Jia X, Qin Y, Sun N, Xin J, Zeng Y, Jing B, Fang J, Pan K, Zeng D, Bai Y, Wang H, Ma H and Ni X (2022) Lactobacillus johnsonii YH1136 plays a protective role against endogenous pathogenic bacteria induced intestinal dysfunction by reconstructing gut microbiota in mice exposed at high altitude. Front. Immunol. 13:1007737. doi: 10.3389/fimmu.2022.1007737
Received: 30 July 2022; Accepted: 06 September 2022;
Published: 10 October 2022.
Edited by:
Wentao Yang, Jilin Agriculture University, ChinaReviewed by:
Muhammad Shoaib, Baku State University, AzerbaijanXiao-Xuan Zhang, Qingdao Agricultural University, China
Copyright © 2022 Wan, Zhang, Jia, Qin, Sun, Xin, Zeng, Jing, Fang, Pan, Zeng, Bai, Wang, Ma and Ni. 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: Xueqin Ni, eHVlcWlubmlAZm94bWFpbC5jb20=; Hailin Ma, RGF2aWRfbWE3OUAxNjMuY29t; Hesong Wang, c3VuYXNhMTAzMEBmb3htYWlsLmNvbQ==
†These authors have contributed equally to this work