- Department of Human Parasitology, School of Basic Medicine, Hubei University of Medicine, Shiyan, China
Cerebral malaria (CM), as one of the most common complications in severe malaria, has threatened millions of individuals’ neurological health and even their lives. Macrophage migration inhibitory factor (MIF), a pleiotropic proinflammatory factor in humans, seems to be a risk factor for death in patients with CM, but its functional mechanism remains unclear. To verify whether affecting the intestinal microbes of the host was one of the mechanisms by which MIF regulates CM, C57BL/6 mice, including WT + PbA, MIF-KO + PbA and their uninfected controls, were sent for 16S rRNA-based sequencing targeting the V4 region of the intestinal microbiota through the Illumina MiSeq platform. The results showed that OTU clustering, alpha and beta diversity in the four groups involved had evident variation. The relative abundance at different taxonomic levels, especially the dominant intestinal flora, was obviously changed. The LEfSe analysis screened out several biomarkers, including significantly reduced Ligilactobacillus (Lactobacillus murinus) in WPbA mice compared to the WT group and Akkermansia (Akkermansia_muciniphila) in KPbA mice compared to the WPbA group. For MIF KO groups, mice infected with PbA or uninfected showed significant enrichment of producers of short-chain fatty acids, including Dubosiella and Faecalibaculum (Faecalibaculum rodentium) in KPbA, and Lachnospiraceae_NK4A136_group and Firmicutes_bacterium_M10-2 in KO. This study not only further proved the gut microbiota changes in C57BL/6 mice caused by PbA infection, but also found that MIF deletion directly affected the changes in the gut microbiota of C57BL/6 mice before and after PbA infection. This finding reveals a potential mechanism by which MIF regulates CM. Combining MIF with potential microbial biomarkers will provide a promising idea to develop combined drugs for improving CM in the future.
Introduction
Malaria caused by Plasmodium infection is a widespread infectious disease in tropical and subtropical regions. According to the latest report by the World Health Organization (WHO; WHO, 2021), the estimated malaria cases reached up to 241 million worldwide in 2020 and were distributed in 85 malaria-endemic areas. Among these, approximately 627,000 deaths occurred. Most of them (95%) appeared in Africa, followed by Southeast Asia (2%), and 77% of patients were children under 5 years of age. Severe malaria is multisyndromic. Its most common manifestations are severe malaria anemia (SMA) and cerebral malaria (CM). Without timely and effective treatment, the mortality rate is fairly high (WHO, 2021). Of these, CM is a serious neurological complication in malaria patients and is mainly caused by infection with Plasmodium falciparum. Even with effective antimalarial therapy and intensive care, CM survivors may also develop long-term neurological deficits (Carter et al., 2005; Boivin et al., 2007; John et al., 2008). Therefore, actively improving the disease process of CM is a beneficial attempt to improve the survival rate and quality of life of patients, especially in high-endemic areas of P. falciparum.
Several studies have shown that human macrophage migration inhibitory factor (MIF), a pleiotropic proinflammatory factor, is closely related to the occurrence and development of various parasitic diseases (Rosado Jde and Rodriguez-Sosa, 2011; Bozza et al., 2012). Studies have found that MIF may be implicated in the pathogenesis of CM (Clark and Awburn, 2002; Clark et al., 2003; Jain et al., 2009). Of these, a study of clinical cases suggested that elevated MIF levels seemed to be a risk factor for death in patients with CM (Jain et al., 2009). These results suggested that MIF may be a potential molecular target for improving CM and even other severe malaria. Conversely, another study found that MIF had a protective role in CM (Awandare et al., 2007), implying that MIF probably has both protective and adverse effects against malaria. However, as a key rationale for the development of MIF-targeted therapeutic regimens, the role and mechanism of MIF in the regulation of malaria, especially CM, remain unclear. In vivo, MIF is abundantly expressed in gastric, intestinal, and colonic epithelia, as well as in the epithelium of human intestinal xenografts (Maaser et al., 2002). Additionally, MIF is involved in the maintenance of physiological microbiota diversity and immune surveillance, thereby enabling normal intestinal barrier functions (Vujicic et al., 2018). The absence of MIF increased the intestinal permeability (Vujičić et al., 2020). These results all indicated the important role of MIF in maintaining host intestinal function. Several studies have shown that malaria may lead to host intestinal lesions (Taniguchi et al., 2015; Shimada et al., 2019). Especially in patients infected with P. falciparum, frequent gastrointestinal symptoms, such as abdominal pain and diarrhea, and changes in the composition of the gut microbiome were observed (Taniguchi et al., 2015). Therefore, a speculation was boldly proposed that perhaps affecting the patient’s intestinal function, such as changing the microbial community of the host intestine, acts as one of the mechanisms by which MIF regulates CM.
Another important reason that the study mainly focused on the regulation of MIF on gut microbiome of CM individuals is as follows. In recent years, regulating the development of diseases by artificially interfering with gut microbes (e.g., adding probiotic yogurt to the daily diet) seems to be an emerging and promising strategy for the prevention and treatment of various diseases (Kataoka, 2016; Clemente et al., 2018; Grochowska et al., 2018). A study by Yilmaz et al. found that specific bacterial components in the human intestinal flora could activate natural defense mechanisms, such as inducing the host to produce high levels of circulating natural antibodies that can recognize and kill Plasmodium, thereby efficiently protecting the host and preventing malaria spread (Yilmaz et al., 2014). Therefore, gut microbiota modulation probably has outstanding potential as a novel approach to prevent or improve severe malaria (Waide and Schmidt, 2020). Additionally, Villarino et al. found that the severity of malaria in mice was profoundly affected by the composition of the gut microbiota, in which, susceptible mice that accepted laboratory-cultured yogurt made from the species Lactobacillus and Bifidobacterium displayed a decreased parasite burden (Villarino et al., 2016). Previous studies showed that the MIF level could be changed by probiotic treatment, indicating an interaction between intestinal microbiota and MIF. Importantly, commonalities of intestinal microbial communities between African children and Plasmodium-susceptible mice suggest that the possibility of probiotics modulating the mouse gut microbiota to control severe malaria may also be suitable for humans (Villarino et al., 2016). Therefore, by using C57BL/6 and PbA to build a murine CM (MCM) model that can mimic human CM (HCM), exploring the intestinal microbial changes caused by Plasmodium infection and the factors that regulate this change (such as MIF) can help to identify biomarkers associated with severe outcomes of CM.
In brief, the study focused on the regulation of MIF on the host gut microbiota in CM models based on WT and MIF knockout C57BL/6 mice infected with PbA to reveal a potential mechanism that MIF regulates CM. The results may bear the expectation that a synergistic treatment of MIF and intestinal flora can more effectively improve the severity of malaria in malaria treatment targeting MIF in the future.
Materials and methods
Mice
Eight- to 10-week-old male C57BL/6 mice (22–30 g weight) were used in the present study. C57BL/6 (wild type, WT) mice were purchased from HNSJA Co., Ltd., Changsha, China. Mif−/− mice with a C57BL/6 genetic background were constructed with the assistance of a professional institute (Gempharmatech Co., Ltd., Jiangsu, China). All animals were maintained in appropriate living conditions (25°C ± 3°C) and were entitled to freely consume a standard diet and pure water daily. Protocols for the use of animals were approved by the Institutional Animal Care and Use Committee of the Hubei University of Medicine (HBMU-S20160414).
Parasite and establishment of murine cerebral malaria
Plasmodium berghei ANKA (PbA) was kindly provided by Professor Wenyue Xu (Third Military Medical University, China) and stored in liquid nitrogen before usage. The parasites were revived by intraperitoneal inoculation in WT C57BL/6 mice, which did not act as donor mice until at least the third passage.
Before the sequencing of intestinal flora, the infection and model of murine cerebral malaria (MCM) were first established and confirmed. In brief, WT (n = 10) and MIF KO (n = 10) mice were intraperitoneally inoculated with approximately 1 × 106 infected red blood cells (iRBCs) from the donor mice in sterile saline (0.9% NaCl) with a volume of 200 μl. From the first day post inoculation (dpi), parasitemia was determined using thin blood smears made from mouse tails and stained with Giemsa. For each mouse, infected red blood cells (iRBCs) in at least 1,000 RBCs were counted under light microscopy. Moreover, assessment of blood–brain barrier permeability was performed as previously reported with some modifications (Taniguchi et al., 2015; Huang et al., 2019). In brief, five mice were randomly taken from both the WT + PbA and KO + PbA groups at 6 dpi and two uninfected mice were injected intravenously with 200 μl of 1% Evans blue dye (Sigma-Aldrich, St Louis, United States) and allowed to circulate in the body for 30 min. Then the mice were anesthetized and sent for heart perfusion with 20 ml of sterile saline (0.9% NaCl). After image documentation, whole brain tissues were placed in formamide for 48 h at 37°C to extract Evans blue. Then, aliquots (100 μl) of the solution of each mouse were analyzed on a microtiter plate for three repeats. Absorbance was measured at 630 nm on an ELISA reader and compared to a preformed standard curve following the document.
After the approach of MCM establishment was confirmed, mice that were randomly assigned to the WT + PbA (WPbA, n = 12) and KO + PbA (KPbA, n = 12) groups were infected with PbA under the same conditions for sequencing of intestinal flora. For the WT control (WT, n = 8) and KO control (KO, n = 7) groups, mice were only injected with sterile saline (0.9% NaCl) in the same volume.
Preparation of gut contents
Gut content samples were collected at 6 dpi, when manifestations of MCM have appeared in most mice, but they had not yet widely died. Only the surviving mice were sacrificed to extrude feces from the whole gut into a sterile EP tube. After being immediately frozen in liquid nitrogen, all gut content samples were transferred to a professional sequencing company (Novogene Bioinformatics Technology Co., Ltd. in Tianjin, China) with dry ice.
DNA extraction, PCR amplification, library construction, and sequencing
Genomic DNA of the fecal samples was extracted using the CTAB or SDS method. After detecting the purity and concentration of DNA by agarose gel electrophoresis, an appropriate amount of sample DNA was diluted to 1 ng/μl with sterile water to be the PCR template. Then, specific primers based on 16S rRNA (515F and 806R, 392 bp), Phusion® High-Fidelity PCR Master Mix with GC Buffer and high-efficiency high-fidelity enzyme from New England Biolabs were used for PCR to ensure the efficiency and accuracy of amplification. PCR products were verified by 2% agarose gel electrophoresis, and only the qualified PCR products were purified by magnetic beads and quantified.
The mixed PCR products in equal amounts were finally purified with the Qiagen Gel Extraction Kit (Qiagen, Germany). Sequencing libraries were constructed with a TruSeq® DNA PCR-Free Sample Preparation Kit (Illumina, United States) and quantified by Qubit and Q-PCR. After the library was qualified, NovaSeq6000 was used for on-machine sequencing.
Sequencing data processing and statistical analysis
According to the barcode sequence and PCR primers, the data of each sample were split from the offline data. After truncating the barcode and primer sequences, FLASH was used to splice the reads of each sample to obtain the original tag data (raw tags; Magoč and Salzberg, 2011). The raw tags were subjected to strict filtering (Bokulich et al., 2013) to obtain high-quality tag data (clean tags). Referring to the tag quality control process of QIIME (Caporaso et al., 2010), the chimera sequences were detected by comparison with the species annotation database (Rognes et al., 2016) and removed (Haas et al., 2011) to obtain the final effective data (effective tags).
The Uparse algorithm was used to cluster all the effective tags of all samples (Edgar, 2013). By default, the sequences were clustered into operational taxonomic units (OTUs) based on 97% identity and the sequences with the highest frequency in the OTUs were screened out as representative sequences. Species annotation analysis was performed using the Mothur method and the SSUrRNA database (Quast et al., 2013) of SILVA138 (with a threshold of 0.8 ~ 1; Wang et al., 2007). Meanwhile, the community composition of each sample was counted at all taxonomic levels: kingdom, phylum, class, order, family, genus and species. A rapid multiple sequence alignment was performed using MUSCLE software to obtain the phylogenetic relationship of all representative sequences of OTUs (Edgar, 2004). Finally, the data of each sample were homogenized, and the sample with the least amount of data was used as the standard for homogenization. The subsequent alpha diversity analysis and beta diversity analysis were based on the homogenized data.
The observed_species, Chao1, Shannon and PD_whole_tree were all calculated in QIIME software (Version 1.9.1). The rarefaction curves and species accumulation curves were all plotted using R software (Version 2.15.3). Differences among groups regarding the alpha diversity index were also obtained in R software using Tukey’s test and the Wilcoxon test.
The UniFrac distance was calculated by QIIME software (Version 1.9.1), and the unweighted pair-group method with arithmetic mean (UPGMA) sample clustering tree was constructed. PCoA graphs were drawn using R software (Version 2.15.3). PCoA analyses were completed with the ade4, ggplot2 and vegan packages in R software. Differences among groups regarding the beta diversity index were also obtained in R software using the Tukey’s test and the Wilcoxon test.
LEfSe analysis was performed using LEfSe software, and the default setting of the LDA score was 4. A permutation test between groups under each classification level (phylum, class, order, family, genus, species) was performed with Metastats analysis by using R software to obtain the p-value, and then the p-value was corrected with the Benjamini and Hochberg False Discovery Rate method to acquire the q-value (White et al., 2009). Anosim and MRPP analyses were finished with the anosim and mrpp functions of the R packages, respectively. Species with significant differences between groups were also analyzed by T-test and plotted with R software.
Results
Confirmation of infection and murine cerebral malaria
As shown in Figure 1A, an average of 23.23% and 25.33% iRBCs were detected in mice from the WPbA and KPbA groups at 6 dpi, respectively. Parasitemia was slightly higher in the KPbA group than in the WPbA group from 4 to 6 dpi. Importantly, the leakage of Evans blue dye into the brains indicated that the collapse of the blood–brain barrier and the development of MCM in both PbA-infected groups occurred at 6 dpi (Figure 1B). Significantly more severe damage to blood–brain barrier permeability appeared in the KPbA mice (Figure 1C). Of course, the closely related clinical symptoms of MCM (Carroll et al., 2010) including ruffled fur, ataxic gait, hunching, limb paralysis, and even convulsions or coma, were also observed at 6 dpi. These results indicated that the infection was sustained and MCM has appeared in the WPbA and KPbA groups at 6 dpi. The absence of MIF seems to play an adverse role in MCM in the present study.
Figure 1. Differences in the pathology between WT and KO mice infected with PbA. (A) Parasitemia was monitored in both the WPbA and KPbA groups, daily. Parasitemia values represent the means ± SD from 10 mice in each group. (B) Mice were injected intravenously with Evans blue at 6 dpi. Uninfected mice were used as controls. (C) The amounts of EB per g of brain from infected mice are shown (n = 5). Statistical analysis was performed using the T-test. *p < 0.05.
By using the same approach for MCM establishment as above, mice prepared for analysis of intestinal flora have also showed MCM signs (except for one and three deaths in WPbA and KPbA before sequencing, respectively) at 6 dpi, including piloerection, mild ataxia in gait and decreased exploratory behavior, strength and reflexes and self-preservation, which are usually followed by severe CM and death, probably at any time of the next several days according to previous experiences and documents (Engwerda et al., 2005; McQuillan et al., 2011; Schmutzhard et al., 2012).
Overview and quality assessment of sequencing data
As shown in Table 1, a total of 2,948,760 raw sequences were generated from all samples after sequencing. Among them, the number of filtered effective tags was 2,219,911 (average effective% = 76.90%), ranging from 58,864 to 68,731 with an average of 63,426 ± 2,883. An average of 98.12 ± 0.13 in Q30 (from 97.79 to 98.36) guaranteed the accuracy of sequencing. Effective tags were clustered with ≥97% sequence identity as the operational taxonomic units (OTUs), and a total of 29,836 OTUs were finally obtained, with an average number of 852 ± 39 (from 752 to 917). The number of observed species was 28,131 in total, and the mean number of all samples was 804 ± 38 (from 699 to 871). Interestingly, the mean number of OTUs in the WPbA group was higher than that in the WT group. However, compared with the KO group, the average number of OTUs was decreased in the KPbA group. The quantity of shared and unique OTUs among groups or different samples within a group were displayed using Venn diagrams or flower figures, respectively (Supplementary Figures 1A,B). Briefly, both PbA-infected groups had a reduction in the amount of shared OTUs relative to their corresponding control groups, with a greater reduction in KPbA mice. The rarefaction curves showed that the species diversity in each sample approached saturation levels, indicating that the sequencing depth had been progressively justified (Supplementary Figure 2A). The species accumulation boxplot tended to be flattened, indicating that the number of samples for sequencing was sufficient (Supplementary Figure 2B).
Table 1. Operational taxonomic unit (OTU)-based diversity indexes of gut microbiota during PbA infection in mice.
Overview of the relative abundance and clustering of species
According to the annotation results, the top 10 phyla, families and genera with the largest abundance in each sample or group were used to generate the column accumulation charts of relative abundance. As shown in Figure 2A, the most abundant phylum was Firmicutes, with an average proportion of 41.71%, followed by Bacteroidota (37.19%). However, for different groups, the proportion of dominant phyla was different. Compared with the WT (35.85%) group, Firmicutes had no significant change in WPbA (36.04%) mice but was obviously increased in KO (48.68%) and KPbA (48.43%) mice. Bacteroidota was slightly decreased in WPbA (40.53%) compared to the WT (43.37%) group and was apparently decreased in the KO (28.09%) and KPbA (34.69%) groups. The remaining seven of the top 10 phyla (excluding unidentified bacteria, 6.78%) were all changed among the four groups and characterized by limited relative abundance, even less than an average proportion of 2.5%.
Figure 2. Variation in microbial community components at different taxonomic levels. (A) Proportions of bacterial communities at the phylum level. (B) Proportions of bacterial communities at the family level. (C) Proportions of bacterial communities at the genus level. Each column represents the constitution and proportion of the top 10 most abundant microbial communities in each individual or group. (D) Taxonomic heatmap of the top 35 most abundant microbial communities at the genus level. The left pattern represents the clustering trees at the genus level. The above profile displays clustering trees among different groups.
As shown in Figure 2B, the most abundant family was Muribaculaceae, with an average proportion of 29.59%, followed by Lachnospiraceae (14.89%), Erysipelotrichaceae (6.75%), and Lactobacillaceae (5.52%). Muribaculaceae and Lactobacillaceae were obviously reduced in the KO (21.92%, 7.97%) and KPbA groups (27.74%, 3.04%) compared with the WT (35.75%, 8.46%) and WPbA (31.51%, 3.84%) groups. In contrast, Lachnospiraceae and Erysipelotrichaceae were overall increased in the KO (16.27%, 8.22%) and KPbA groups (16.06%, 11.98%) compared with the WT (12.37%, 4.10%) and WPbA (14.88%, 3.45%) groups.
As shown in Figure 2C, the most abundant genus was Lachnospiraceae_NK4A136_group, with an average proportion of 5.06%, followed by Ligilactobacillus (4.43%), Dubosiella (3.91%) and Bacteroides (3.75%). Ligilactobacillus and Bacteroides were obviously reduced in the KO (6.18%, 2.82%) and KPbA groups (1.79%, 3.77%) compared with the WT (7.59%, 3.42%) and WPbA (3.19%, 4.56%) groups, respectively. Conversely, Lachnospiraceae_NK4A136_group, Dubosiella and Faecalibaculum were increased overall in the KO (6.36%, 4.89%, and 2.01%) and KPbA groups (5.52%, 5.88%, and 5.17%) compared with the WT (3.44%, 3.15%, and 0.80%) and WPbA (5.05%, 2.21%, and 0.81%) groups, respectively. It is interesting that Akkermansia was increased by approximately threefold in the WPbA group compared to the WT group, but decreased by twofold in the KPbA group compared to the KO group. Additionally, the top 35 genera in abundance were clustered at both the species and sample levels, and the heatmap of the results was shown in Figure 2D. Overall, mice in the WT and WPbA groups tended to cluster together, and the same applied to the samples in the KO and KPbA groups. Compared to the WT group, higher amounts of Oscillibacter, Muribaculum, Akkermansia, Escherichia-Shigella, Parabacteroides and Bacteroides, and decreased amounts of Dubosiella, Desulfovibrio and Lactobacillus were enriched in the WPbA group. However, for the KO group, the most abundant species clustered were Intestinimonas, Bifidobacterium, Mycoplasma, and Lactobacillus. Compared to the KO group, the KPbA group was characterized by higher amounts of [Eubacterium]_fissicatena_group, [Eubacterium] _ ruminantium_group, Faecalibaculum and Roseburia, and decreased amounts of [Eubacterium]_siraeum_group, Alistipes, Odoribacter and Candidatus_Saccharimonas.
Diversity of bacterial community within groups
The richness and diversity of microbes within communities were assessed by alpha diversity analysis, including observed_species, Chao1, Shannon and PD_whole_tree (Table 1; Figure 3). Compared with the WT group, the average number of OTUs and the median of observed_species, Chao1 and Shannon in the WPbA group were all increased. In contrast, compared with the KO group, the above indicators in the KPbA group were all decreased. Slightly different, for PD_whole_tree, a gradual downward trend in the four groups WT, WPbA, KO, and KPbA was observed.
Figure 3. The richness and diversity of microbes within communities by alpha analysis. (A) Observed_species, (B) Chao1, (C) Shannon, and (D) PD_whole_tree. Differences among groups regarding the alpha diversity index were obtained in R software using the Wilcoxon test. Significant differences are indicated via **p < 0.01 and ***p < 0.001.
Observed_species (p = 0.0074, p = 0.0006, respectively; Wilcox rank test, respectively) and Chao1 (p = 0.0018, p = 0.0001; Wilcox rank sum test, respectively) were all significantly decreased in the KPbA group when compared with the WPbA and KO groups. However, the median values of these two indexes were all increased in the WPbA and KO groups compared with the WT group, but there was no significant difference (Figures 3A,B). The Shannon value was significantly elevated in the WPbA group (p = 0.0017; Wilcoxon rank sum test) compared with the WT group, but conversely, it was significantly decreased in the KPbA group (p = 0.0032; Wilcoxon rank-sum test) compared with the WPbA group (Figure 3C). In contrast, PD_whole_tree was significantly decreased in the KO group (p = 0.0024; Wilcoxon rank-sum test) compared to the WT group and significantly reduced in the KPbA group (p = 0.0009; Wilcoxon rank-sum test) compared with the WPbA group (Figure 3D).
Comparative analysis of the bacterial communities among groups
Beta diversity was used for comparative analysis of the microbial community composition among different groups. Unweighted UniFrac analysis showed that medians were significantly lower in both WPbA (p = 0.0041; Wilcox rank-sum test) and KO (p = 0.0004; Wilcox rank-sum test) mice than in WT mice. However, it was significantly higher in the KPbA mice (p = 0.0001, p = 0.0008, respectively; Wilcoxon rank sum test) than in the KO and WPbA groups (Figure 4A). In contrast, the median values among the four groups showed no significant difference by weighted UniFrac analysis (Figure 4B).
Figure 4. Comparative analysis of the microbial community composition among different groups by beta diversity. (A) Unweighted UniFrac, (B) weighted UniFrac, (C) principal coordinates analysis (PCoA) based on unweighted, and (D) weighted UniFrac distance matrixes for intestinal contents from four groups involved. Significant differences are indicated via **p < 0.01.
To assess the structural similarity of species composition, PCoA based on unweighted UniFrac and weighted UniFrac was performed on each sample (Figures 4C,D). The results presented a relatively consistent clustering trend by the two methods. Anosim and MRPP are two similar analytical methods that can be used for significance testing of the differences in microbial community structure between groups. The results showed that differences among all groups were significant and greater than those within group samples (Table 2), indicating that both MIF gene knockout and PbA infection could cause significant changes in the gut microbial community. Additionally, the 35 samples were subjected to UPGMA cluster analysis with a weighted UniFrac distance matrix (Supplementary Figure 3). The results revealed that the WT and WPbA groups or the KO and KPbA groups were clustered into two separate branches at the phylum level.
Table 2. Within- vs. among-group dissimilarities analysis via analysis of similarities (ANOSIM) and multi response permutation procedure (MRPP).
Potential biomarker discovery
The LDA effect size (LEfSe) analysis was able to find statistically significant biomarkers between groups. As shown in Figure 5, potential biomarkers existed at different taxonomic levels. At the genus and species level, the biomarker of WT mice was Ligilactobacillus (Lactobacillus murinus) compared with WPbA but was Faecalibaculum (Faecalibaculum rodentium) for the KPbA compared with the KO groups (Figures 5A,B). Compared to WT mice, Lachnospiraceae_NK4A136_group (genus) and Firmicutes_bacterium_M10-2 (species) were highlighted as significantly different biomarkers in the KO group (Figure 5C). Interestingly, when comparing WPbA with KPbA, the biomarkers for WPbA were Akkermansia (Akkermansia_muciniphila), and the biomarkers for KPbA included Dubosiella and Faecalibaculum (Faecalibaculum rodentium; Figure 5D).
Figure 5. The differentially abundant bacterial taxa identified by LEfSe among different groups. (A) Linear discriminant analysis (LDA) score distribution between the WT and WPbA groups. (B) LDA score distribution between the KO and KPbA groups. (C) LDA score distribution between the WT and KO groups. (D) LDA score distribution between the WPbA and KPbA groups. The histogram of LDA value distribution only showed the species whose LDA score was greater than the default setting (threshold > 4), which represented the biomarkers with significant differences between groups.
Discussion
The elevated level of MIF, an important immune molecule in humans, was reported to be a possible risk factor for death in patients with CM (Jain et al., 2009), making it a potential molecular target for improving severe malaria. In contrast, another study revealed that MIF had a protective role in CM (Awandare et al., 2007), implying that MIF probably has both protective and adverse effects against malaria. In general, its actions and mechanisms seem very complicated and remain unclear. Altering the microbial community of patients may be one of the pathways by which MIF affects CM. The results of gut microbiota sequencing showed that the OTU clustering, as well as the alpha and beta diversity in the four groups of the present study all had obvious changes.
For the top two dominant phyla, Firmicutes was elevated by ~12% in both KO and KPbA mice, but Bacteroidota was decreased by 15.28% in KO and 5.8% in KPbA compared to the WT and WPbA groups, respectively. A previous study reported that immune protection in the intestine depends on the development of Bacteroidetes (Sequeira et al., 2020). Therefore, its reduction in the KO and KPbA groups may hint that the deficiency of MIF would weaken, at least partially, microbial immune protection in the intestine of the PbA-infected and uninfected hosts. Additionally, the Firmicutes/Bacteroidetes (F/B) ratios in KO and KPbA mice were approximately 2.1 and 1.6 times as much as those in the WT and WPbA groups, respectively. Studies have shown that F/B ratios were significantly higher in prostate enlargement patients (Takezawa et al., 2021), spontaneously hypertensive rats (Yang et al., 2015), autistic individuals and spinal-onset amyotrophic lateral sclerosis patients (Strati et al., 2017). These results seemed to suggest that the increase in F/B was probably a potential risk marker for gut dysbiosis in MIF-deficient mice. Muribaculaceae, as the dominant family in the Bacteroidota phylum, had a similar change trend. Muribaculaceae was decreased by 13.83% in the KO and 3.77% in the KPbA groups compared to the WT and WPbA groups, respectively. Studies have shown that members of Muribaculaceae represent major foragers of mucin monosaccharide (Pereira et al., 2020), which forms a natural barrier in the human intestine to avoid intimate contact between virulent microbes and the intestinal wall. Therefore, MIF deficiency appears to affect uninfected and PbA-infected individuals by reducing the abundance of intestinal mucin-degrading bacteria. Conversely, Erysipelotrichaceae was increased by 4.12% in the KO and 8.53% in the KPbA groups compared to the WT and WPbA groups, respectively. Members of this family appear to be highly immunogenic, play an important role in gastrointestinal inflammation-related diseases, and cause host lipid-related metabolic disturbances (Kaakoush, 2015). Characterizing important changes in this family in KPbA mice may ultimately provide promising microbial targets to combat the metabolic disorders caused by malaria. As one of the most common probiotic families (Zawistowska-Rojek et al., 2022), Lactobacillaceae was significantly lower in the WPbA and KPbA groups than in their corresponding control groups, but was less affected by MIF knockout, indicating that it is mainly PbA infection that caused the decrease in Lactobacillaceae.
The LEfSe analysis in this study identified a number of biomarkers with significant differences between groups at different taxonomic levels. Here, we focus on the species level. When compared with WT, the relative abundance of Lactobacillus murinus (belonging to the genus Ligilactobacillus), one of the natural commensal gut bacteria in the intestine of healthy individuals (Lebovitz and Theus, 2019), was dramatically decreased in WPbA. Several studies have revealed the positive functions of L. murinus in antimicrobial production (Nardi et al., 2005), protection of the intestinal barrier (Delucchi et al., 2017) and neonatal necrotizing enterocolitis (Isani et al., 2018), indicating that L. murinus was probably also a promising probiotic candidate for CM improvement. When compared with the KO group, the biomarker was Faecalibaculum (Faecalibaculum rodentium) for the KPbA group. Similarly, when compared with WPbA, the biomarkers for KPbA included significantly increased Faecalibaculum (Faecalibaculum rodentium) and Dubosiell and decreased Akkermansia (Akkermansia_muciniphila). Interestingly, these three genera/species microbiota were all beneficial for the production of short-chain fatty acids (SCFAs; Shen et al., 2021; Ye et al., 2021), which, as one of the safeguards for maintaining gut health, plays an important role in strengthening epithelial barrier integrity and inhibiting inflammation. However, a previous study found that oral administration of F. rodentium in antibiotic-treated resilient Ephx2 KO mice led to depression-like behaviors (Wang et al., 2021). Another study reported that Dubosiella was positively correlated with glycolipid metabolism disorders (Miao et al., 2022). Akkermansia spp., as important mucus-degrading bacteria, have significant anti-inflammatory activity and are considered as promising probiotics (Romano et al., 2021; Ye et al., 2021). Meanwhile, significant increases in Lachnospiraceae_NK4A136_group (genus) and Firmicutes_bacterium_M10-2 (species) were also highlighted in the KO group compared with the WT group. The Lachnospiraceae_NK4A136_group represents a gut bacterium naturally inhabiting healthy individuals (Stadlbauer et al., 2020). Interestingly, it is also a producer of butyrate, which is an important SCFA (Gheorghe et al., 2022) and has a positive role in maintaining the gut barrier integrity of mice (Hu et al., 2019). However, the function of Firmicutes_bacterium_M10-2 is still poorly understood in both humans and rodents. In brief, two interesting findings were discovered in the LEfSe analysis: 1) The probiotic Ligilactobacillus (Lactobacillus murinus) was significantly reduced in WPbA mice compared to the WT group, and a candidate probiotic Akkermansia (Akkermansia_muciniphila) was significantly reduced in the KPbA mice compared to the WPbA group. 2) MIF knockout mice, uninfected or infected with PbA, showed significant enrichment of producers of short-chain fatty acids; however, the predominant biomarkers were different between the two groups. The above results indicated that PbA infection could lead to an apparent decrease in the abundance of different probiotics or potential probiotics in WT and MIF KO mice. However, whether MIF deficiency is likely to play a beneficial role in malaria by upregulating the abundance of producers of short-chain fatty acids or otherwise exert adverse effects through other pathways, as previously reported or in an unknown way, still needs further exploration. Most likely, the role of different probiotics or potential probiotics changes in food habits of individuals as well as the pathogenesis of parasitic infection in human beings, which may culminate with short-chain fatty acids effects through different pathways.
In addition, Taniguchi et al. found that PbA infection resulted in significant changes in host intestinal pathology and microbiota (Taniguchi et al., 2015). Guan et al. found that the gut microbiota of C57BL/6 mice had been rebuilt on day 4 after PbA infection (Guan et al., 2021). However, different changes in the gut microbiome and significant biomarkers were reported in their studies. Studies have even found that mice from different commercial suppliers showed differences in gut bacterial communities and affected their susceptibility to Plasmodium (Ivanov et al., 2008, 2009). This may suggest that differences in gut microbiota are an important factor in whether different individuals develop severe malaria in the clinic. Therefore, the development of adjuvant therapy based on probiotic or harmful biomarkers screened in a single study will have great limitations. In other words, the bacterial flora changes and advantageous biomarkers reported in different studies may be different, but the consistent changes and shared significant biomarkers in different studies, in our opinion, are probably the key factors that affect the severity of different malaria patients and act as the most potential targets for improving severe malaria, including CM. By conducting comprehensive analysis on the limited documents about the influence of intestinal microbes on CM (Taniguchi et al., 2015; Guan et al., 2021) and biomarkers in the present study, it seems that the genus Lactobacillus that was frequently screened out probably be most responsible for CM.
This study revealed interesting changes in the gut microbiota of MIF KO and WT C57BL/6 mice before and after PbA infection. However, this study also had some limitations. Result in the present study showed an adverse role of MIF deficiency in mice with CM. Nevertheless, it is unknown how MIF affects the severity and pathogenesis of CM after changing the intestinal flora. A dynamic observation of the gut microbiota and the corresponding pathological, metabolic and immune changes is lacking. In future studies, we should also fully explore whether the addition of the identified potential probiotic biomarkers or the elimination of harmful bacteria with specific antibiotics can reduce the severity of CM and their microorganism-dependent mechanisms.
In summary, this study further proved the gut microbiota changes in C57BL/6 mice caused by PbA infection and found that MIF deletion directly affected the changes in the gut microbiota of C57BL/6 mice before and after PbA infection. This finding reveals a potential mechanism by which MIF regulates CM. The obtained potential biomarkers combined with MIF will also provide a promising idea to develop combined drugs for improving CM in the future.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA852657.
Ethics statement
The animal study was reviewed and approved by Institutional Animal Care and Use Committee of the Hubei University of Medicine (HBMU-S20160414).
Author contributions
WG, YZ, YX, and JL conceived the study, participated in its design, performed the data analysis and interpretation, and drafted the manuscript. WG, YZ, SY, KG, SC, XH, HS, and YX carried out the experiments. WG, YZ, YX, and JL participated in analyzing and interpreting the data. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by grants from the Post-Graduates of Hubei University of Medicine (2018QDJZR17 and 2016QDJZR04), foundation from the Research Project of Hubei Provincial Department of Education (Q20192102 and Q20172102), and partially supported by the Principle Investigator Program of Hubei University of Medicine (HBMUPI202101).
Acknowledgments
The authors sincerely thank Wenyue Xu (Third Military Medical University) for the Plasmodium berghei ANKA strain (PbA) parasite gift.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.978644/full#supplementary-material
References
Awandare, G. A., Kremsner, P. G., Hittner, J. B., Keller, C. C., Clark, I. A., Weinberg, J. B., et al. (2007). Higher production of peripheral blood macrophage migration inhibitory factor in healthy children with a history of mild malaria relative to children with a history of severe malaria. Am. J. Trop. Med. Hyg. 76, 1033–1036. doi: 10.4269/ajtmh.2007.76.1033
Boivin, M. J., Bangirana, P., Byarugaba, J., Opoka, R. O., Idro, R., Jurek, A. M., et al. (2007). Cognitive impairment after cerebral malaria in children: a prospective study. Pediatrics 119, e360–e366. doi: 10.1542/peds.2006-2027
Bokulich, N. A., Subramanian, S., Faith, J. J., Gevers, D., Gordon, J. I., Knight, R., et al. (2013). Quality-filtering vastly improves diversity estimates from illumina amplicon sequencing. Nat. Methods 10, 57–59. doi: 10.1038/nmeth.2276
Bozza, M. T., Martins, Y. C., Carneiro, L. A., and Paiva, C. N. (2012). Macrophage migration inhibitory factor in protozoan infections. J. Parasitol. Res. 2012:413052. doi: 10.1155/2012/413052
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Carroll, R. W., Wainwright, M. S., Kim, K.-Y., Kidambi, T., Gómez, N. D., Taylor, T., et al. (2010). A rapid murine coma and behavior scale for quantitative assessment of murine cerebral malaria. PLoS One 5:e13124. doi: 10.1371/journal.pone.0013124
Carter, J. A., Mung'ala-Odera, V., Neville, B. G., Murira, G., Mturi, N., Musumba, C., et al. (2005). Persistent neurocognitive impairments associated with severe falciparum malaria in Kenyan children. J. Neurol. Neurosurg. Psychiatry 76, 476–481. doi: 10.1136/jnnp.2004.043893
Clark, I., and Awburn, M. (2002). Migration inhibitory factor in the cerebral and systemic endothelium in sepsis and malaria. Crit. Care Med. 30, S263–S267. doi: 10.1097/00003246-200205001-00015
Clark, I. A., Awburn, M. M., Whitten, R. O., Harper, C. G., Liomba, N. G., Molyneux, M. E., et al. (2003). Tissue distribution of migration inhibitory factor and inducible nitric oxide synthase in falciparum malaria and sepsis in African children. Malar. J. 2:6. doi: 10.1186/1475-2875-2-6
Clemente, J. C., Manasson, J., and Scher, J. U. (2018). The role of the gut microbiome in systemic inflammatory disease. BMJ 360:j5145. doi: 10.1136/bmj.j5145
Delucchi, L., Fraga, M., and Zunino, P. (2017). Effect of the probiotic Lactobacillus murinus lbp2 on clinical parameters of dogs with distemper-associated diarrhea. Can. J. Vet. Res. 81, 118–121.
Edgar, R. C. (2004). Muscle: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797. doi: 10.1093/nar/gkh340
Edgar, R. C. (2013). Uparse: highly accurate otu sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Engwerda, C., Belnoue, E., Grüner, A. C., and Rénia, L. (2005). Experimental models of cerebral malaria. Curr. Top. Microbiol. Immunol. 297, 103–143.
Gheorghe, A. S., Negru, Ș. M., Preda, M., Mihăilă, R. I., Komporaly, I. A., Dumitrescu, E. A., et al. (2022). Biochemical and metabolical pathways associated with microbiota-derived butyrate in colorectal cancer and omega-3 fatty acids implications: a narrative review. Nutrients 14:1152. doi: 10.3390/nu14061152
Grochowska, M., Wojnar, M., and Radkowski, M. (2018). The gut microbiota in neuropsychiatric disorders. Acta Neurobiol. Exp. 78, 69–81. doi: 10.21307/ane-2018-008
Guan, W., Yang, S., Zhao, Y., Cheng, W., Song, X., Yao, Y., et al. (2021). Observation of the gut microbiota profile in c57bl/6 mice induced by Plasmodium berghei ANKA infection. Front. Cell. Infect. Microbiol. 11:680383. doi: 10.3389/fcimb.2021.680383
Haas, B. J., Gevers, D., Earl, A. M., Feldgarden, M., Ward, D. V., Giannoukos, G., et al. (2011). Chimeric 16S rRNA sequence formation and detection in sanger and 454-pyrosequenced PCR amplicons. Genome Res. 21, 494–504. doi: 10.1101/gr.112730.110
Hu, S., Wang, J., Xu, Y., Yang, H., Wang, J., Xue, C., et al. (2019). Anti-inflammation effects of fucosylated chondroitin sulphate from Acaudina molpadioides by altering gut microbiota in obese mice. Food Funct. 10, 1736–1746. doi: 10.1039/c8fo02364f
Huang, B., Huang, S., Chen, X., Liu, X. B., Wu, Q., Wang, Y., et al. (2019). Activation of mast cells promote Plasmodium berghei ANKA infection in murine model. Front. Cell. Infect. Microbiol. 9:322. doi: 10.3389/fcimb.2019.00322
Isani, M., Bell, B. A., Delaplain, P. T., Bowling, J. D., Golden, J. M., Elizee, M., et al. (2018). Lactobacillus murinus hf12 colonizes neonatal gut and protects rats from necrotizing enterocolitis. PLoS One 13:e0196710. doi: 10.1371/journal.pone.0196710
Ivanov, I. I., Atarashi, K., Manel, N., Brodie, E. L., Shima, T., Karaoz, U., et al. (2009). Induction of intestinal th17 cells by segmented filamentous bacteria. Cell 139, 485–498. doi: 10.1016/j.cell.2009.09.033
Ivanov, I. I., Frutos, R. D. L., Manel, N., Yoshinaga, K., Rifkin, D. B., Sartor, R. B., et al. (2008). Specific microbiota direct the differentiation of il-17-producing t-helper cells in the mucosa of the small intestine. Cell Host Microbe 4, 337–349. doi: 10.1016/j.chom.2008.09.009
Jain, V., McClintock, S., Nagpal, A. C., Dash, A. P., Stiles, J. K., Udhayakumar, V., et al. (2009). Macrophage migration inhibitory factor is associated with mortality in cerebral malaria patients in India. BMC. Res. Notes 2:36. doi: 10.1186/1756-0500-2-36
John, C. C., Bangirana, P., Byarugaba, J., Opoka, R. O., Idro, R., Jurek, A. M., et al. (2008). Cerebral malaria in children is associated with long-term cognitive impairment. Pediatrics 122, e92–e99. doi: 10.1542/peds.2007-3709
Kaakoush, N. O. (2015). Insights into the role of erysipelotrichaceae in the human host. Front. Cell. Infect. Microbiol. 5:84. doi: 10.3389/fcimb.2015.00084
Kataoka, K. (2016). The intestinal microbiota and its role in human health and disease. J. Med. Invest. 63, 27–37. doi: 10.2152/jmi.63.27
Lebovitz, Y., and Theus, M. H. (2019). Molecular phenotyping and genomic characterization of a novel neuroactive bacterium strain, Lactobacillus murinus hu-1. Front. Pharmacol. 10:1162. doi: 10.3389/fphar.2019.01162
Maaser, C., Eckmann, L., Paesold, G., Kim, H. S., and Kagnoff, M. F. (2002). Ubiquitous production of macrophage migration inhibitory factor by human gastric and intestinal epithelium. Gastroenterology 122, 667–680. doi: 10.1053/gast.2002.31891
Magoč, T., and Salzberg, S. L. (2011). Flash: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
McQuillan, J. A., Mitchell, A. J., Ho, Y. F., Combes, V., Ball, H. J., Golenser, J., et al. (2011). Coincident parasite and cd8 t cell sequestration is required for development of experimental cerebral malaria. Int. J. Parasitol. 41, 155–163. doi: 10.1016/j.ijpara.2010.08.003
Miao, M., Wang, Q., Wang, X., Fan, C., Luan, T., Yan, L., et al. (2022). The protective effects of inulin-type fructans against high-fat/sucrose diet-induced gestational diabetes mice in association with gut microbiota regulation. Front. Microbiol. 13:832151. doi: 10.3389/fmicb.2022.832151
Nardi, R. M., Santoro, M. M., Oliveira, J. S., Pimenta, A. M., Ferraz, V. P., Benchetrit, L. C., et al. (2005). Purification and molecular characterization of antibacterial compounds produced by Lactobacillus murinus strain l1. J. Appl. Microbiol. 99, 649–656. doi: 10.1111/j.1365-2672.2005.02632.x
Pereira, F. C., Wasmund, K., Cobankovic, I., Jehmlich, N., Herbold, C. W., Lee, K. S., et al. (2020). Rational design of a microbial consortium of mucosal sugar utilizers reduces Clostridiodes difficile colonization. Nat. Commun. 11:5104. doi: 10.1038/s41467-020-18928-1
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The silva ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Rognes, T., Flouri, T., Nichols, B., Quince, C., and Mahé, F. (2016). Vsearch: a versatile open source tool for metagenomics. PeerJ 4:e2584. doi: 10.7717/peerj.2584
Romano, S., Savva, G. M., Bedarf, J. R., Charles, I. G., Hildebrand, F., and Narbad, A. (2021). Meta-analysis of the Parkinson's disease gut microbiome suggests alterations linked to intestinal inflammation. NPJ Parkinsons Dis. 7:27. doi: 10.1038/s41531-021-00156-z
Rosado Jde, D., and Rodriguez-Sosa, M. (2011). Macrophage migration inhibitory factor (MIF): A key player in protozoan infections. Int. J. Biol. Sci. 7, 1239–1256. doi: 10.7150/ijbs.7.1239
Schmutzhard, J., Kositz, C. H., Glueckert, R., Schmutzhard, E., Schrott-Fischer, A., and Lackner, P. (2012). Apoptosis of the fibrocytes type 1 in the spiral ligament and blood labyrinth barrier disturbance cause hearing impairment in murine cerebral malaria. Malar. J. 11:30. doi: 10.1186/1475-2875-11-30
Sequeira, R. P., McDonald, J. A. K., Marchesi, J. R., and Clarke, T. B. (2020). Commensal bacteroidetes protect against klebsiella pneumoniae colonization and transmission through il-36 signalling. Nat. Microbiol. 5, 304–313. doi: 10.1038/s41564-019-0640-1
Shen, T., Yue, Y., He, T., Huang, C., Qu, B., Lv, W., et al. (2021). The association between the gut microbiota and Parkinson's disease, a meta-analysis. Front. Aging Neurosci. 13:636545. doi: 10.3389/fnagi.2021.636545
Shimada, M., Hirose, Y., Shimizu, K., Yamamoto, D. S., Hayakawa, E. H., and Matsuoka, H. (2019). Upper gastrointestinal pathophysiology due to mouse malaria Plasmodium berghei ANKA infection. Trop. Med. Health 47:18. doi: 10.1186/s41182-019-0146-9
Stadlbauer, V., Engertsberger, L., Komarova, I., Feldbacher, N., Leber, B., Pichler, G., et al. (2020). Dysbiosis, gut barrier dysfunction and inflammation in dementia: a pilot study. BMC Geriatr. 20:248. doi: 10.1186/s12877-020-01644-2
Strati, F., Cavalieri, D., Albanese, D., De Felice, C., Donati, C., Hayek, J., et al. (2017). New evidences on the altered gut microbiota in autism spectrum disorders. Microbiome 5:24. doi: 10.1186/s40168-017-0242-1
Takezawa, K., Fujita, K., Matsushita, M., Motooka, D., Hatano, K., Banno, E., et al. (2021). The firmicutes/bacteroidetes ratio of the human gut microbiota is associated with prostate enlargement. Prostate 81, 1287–1293. doi: 10.1002/pros.24223
Taniguchi, T., Miyauchi, E., Nakamura, S., Hirai, M., Suzue, K., Imai, T., et al. (2015). Plasmodium berghei ANKA causes intestinal malaria associated with dysbiosis. Sci. Rep. 5:15699. doi: 10.1038/srep15699
Villarino, N. F., LeCleir, G. R., Denny, J. E., Dearth, S. P., Harding, C. L., Sloan, S. S., et al. (2016). Composition of the gut microbiota modulates the severity of malaria. Proc. Natl. Acad. Sci. U. S. A. 113, 2235–2240. doi: 10.1073/pnas.1504887113
Vujičić, M., Despotović, S., Saksida, T., and Stojanović, I. (2020). The effect of macrophage migration inhibitory factor on intestinal permeability: Fitc-dextran serum measurement and transmission electron microscopy. Methods Mol. Biol. 2080, 193–201. doi: 10.1007/978-1-4939-9936-1_17
Vujicic, M., Saksida, T., Despotovic, S., Bajic, S. S., Lalić, I., Koprivica, I., et al. (2018). The role of macrophage migration inhibitory factor in the function of intestinal barrier. Sci. Rep. 8:6337. doi: 10.1038/s41598-018-24706-3
Waide, M. L., and Schmidt, N. W. (2020). The gut microbiome, immunity, and plasmodium severity. Curr. Opin. Microbiol. 58, 56–61. doi: 10.1016/j.mib.2020.08.006
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/aem.00062-07
Wang, S., Ishima, T., Qu, Y., Shan, J., Chang, L., Wei, Y., et al. (2021). Ingestion of Faecalibaculum rodentium causes depression-like phenotypes in resilient ephx2 knock-out mice: A role of brain-gut-microbiota axis via the subdiaphragmatic vagus nerve. J. Affect. Disord. 292, 565–573. doi: 10.1016/j.jad.2021.06.006
White, J. R., Nagarajan, N., and Pop, M. (2009). Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput. Biol. 5:e1000352. doi: 10.1371/journal.pcbi.1000352
Yang, T., Santisteban, M. M., Rodriguez, V., Li, E., Ahmari, N., Carvajal, J. M., et al. (2015). Gut dysbiosis is linked to hypertension. Hypertension 65, 1331–1340. doi: 10.1161/hypertensionaha.115.05315
Ye, X., Liu, Y., Hu, J., Gao, Y., Ma, Y., and Wen, D. (2021). Chlorogenic acid-induced gut microbiota improves metabolic endotoxemia. Front. Endocrinol. 12:762691. doi: 10.3389/fendo.2021.762691
Yilmaz, B., Portugal, S., Tran, T. M., Gozzelino, R., Ramos, S., Gomes, J., et al. (2014). Gut microbiota elicits a protective immune response against malaria transmission. Cell 159, 1277–1289. doi: 10.1016/j.cell.2014.10.053
Keywords: Plasmodium berghei, macrophage migration inhibitory factor, cerebral malaria, gut microbiota, 16S rRNA
Citation: Xie Y, Guan W, Zhao Y, Yan S, Guo K, Chen S, Hu X, Shi H and Li J (2022) Deficiency of migration inhibitory factor influences the gut microbiota of C57BL/6 mice infected with Plasmodium berghei ANKA. Front. Microbiol. 13:978644. doi: 10.3389/fmicb.2022.978644
Edited by:
George Grant, University of Aberdeen, United KingdomReviewed by:
Natarajan Gopalan, Central University of Tamil Nadu, IndiaHeny Arwati, Airlangga University, Indonesia
Fátima Brant, Federal University of Minas Gerais, Brazil
Copyright © 2022 Xie, Guan, Zhao, Yan, Guo, Chen, Hu, Shi and Li. 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: Jian Li, eXhsaWppYW5AMTYzLmNvbQ==
†These authors have contributed equally to this work