Abstract
Despite the impact of childhood diarrhea on morbidity and mortality, our understanding of its sequelae has been significantly hampered by the lack of studies that examine samples across the entire intestinal tract. Infant rhesus macaques are naturally susceptible to human enteric pathogens and recapitulate the hallmarks of diarrheal disease such as intestinal inflammation and growth faltering. Here, we examined intestinal biopsies, lamina propria leukocytes, luminal contents, and fecal samples from healthy infants and those experiencing growth faltering with distant acute or chronic active diarrhea. We show that growth faltering in the presence or absence of active diarrhea is associated with a heightened systemic and mucosal pro-inflammatory state centered in the colon. Moreover, polyclonal stimulation of colonic lamina propria leukocytes resulted in a dampened cytokine response, indicative of immune exhaustion. We also detected a functional and taxonomic shift in the luminal microbiome across multiple gut sites including the migration of Streptococcus and Prevotella species between the small and large intestine, suggesting a decompartmentalization of gut microbial communities. Our studies provide valuable insight into the outcomes of diarrheal diseases and growth faltering not attainable in humans and lays the groundwork to test interventions in a controlled and reproducible setting.
Similar content being viewed by others
Introduction
Despite significant improvements in healthcare delivery and sanitation, over 500,000 children under the age of 5 die of diarrheal disease every year1. In addition to high mortality rates, diarrhea early in life can have a lasting effect on host physiology2. Specifically, diarrhea results in a dangerous feedback loop of malnutrition3,4,5,6 and intestinal damage3,7,8, which ultimately leads to growth faltering, poor responses to oral vaccination9,10,11,12, and diminished resilience to subsequent enteric infections13. This cycle can occur even in the absence of overt diarrhea14,15. Indeed, growth faltering children in developing countries have been shown to have a significant enteric pathogen burden without acute disease14,16,17. These subclinical infections are believed to perpetuate intestinal damage and restrict healthy growth18. Moreover, unknown and emerging etiological agents are thought to play an underappreciated role in childhood diarrheal disease and growth faltering. Specifically, multiple large-scale international studies have found that less than 50% of diarrheal cases can be attributed to a specific pathogen19,20. This pattern is likely a combination of high pathogen burden in the absence of symptoms19,20,21 and the use of clinical tests that limit detection of novel pathogens not routinely detected22,23. More recent studies have begun to identify a set of non-traditional pathogens, especially in the genus Campylobacter, that are missed by clinical tests but highly abundant in developing countries24,25,26. Additional studies are needed to determine the role of these newly identified pathogens in the cycle of diarrheal disease, malnutrition, and growth faltering.
In addition, few studies have investigated diarrheal disease-induced damage within the gut itself. One study reported increased gut permeability and inflammation markers in duodenal biopsies obtained from Gambian infants with growth faltering compared to those obtained from healthy children living in the United Kingdom9. While this study reported an association between growth faltering and damage to the upper small intestine, the involvement of other gut sites remains unknown. More importantly, the interpretation of these observations is complicated by the lack of comparison to samples obtained from asymptomatic children living under the same conditions. Transcriptional analysis of fecal samples obtained from infants with environmental enteric dysfunction revealed increased expression of immune genes27; however, the cellular source of these expression changes remains unclear.
Alterations in mucosal immune cell frequency and function across the GI tract due to growth faltering/diarrheal disease have yet to be identified due to the challenges associated with the collection of intestinal biopsies from healthy infants and those experiencing growth faltering while controlling for key variables such as geographic location, diet, and socioeconomic status. A deeper understanding of the pleiotropic impact of diarrhea on GI immunological and microbial community composition is key to successfully identify sub-clinical cases and developing targeted treatments against growth faltering and malnutrition.
Malnutrition, growth faltering, and diarrhea have all been associated with dysbiosis of the gut microbiome at a taxonomic level using fecal samples28,29,30, as well as the persistence of an immature gut microbiome state31,32. Moreover, growth faltering is associated with the decompartmentalization of microbial communities as illustrated by the recent report of potential infiltration of nasopharyngeal microbes into the GI tract33. Interestingly, levels of enteric pathogens such as Escherichia coli and Campylobacter spp. are elevated in the fecal microbiome of children with malnutrition and growth faltering regardless of diarrhea16,25,34,35. However, the vast majority of clinical findings associated with diarrheal disease and growth faltering are centered around the small intestine, while the majority of microbiome data is generated from fecal samples. This dichotomy may conceal clinically relevant information and no studies to date have interrogated microbial changes within the multiple intestinal tract sections at the taxonomic or functional level.
Answering these questions would be greatly facilitated by the availability of an animal model that was susceptible to human enteric pathogens and displayed the hallmarks of growth faltering. Rhesus macaques are the gold standard model for the study of several infectious diseases due to their genetic and physiological similarity to humans. Rhesus macaques housed at the Oregon National Primate Research Center (ONPRC) and California National Primate Research Center (CNPRC) naturally develop the diarrheal disease at a predictable yearly rate with Campylobacter and Shigella being the primary pathogens isolated at clinic admission36,37,38. In addition, infant macaques that experience diarrhea suffer growth faltering37. Rhesus macaques also display a high pathogen burden even in the absence of diarrheal disease, as described in infants19. Moreover, we recently reported that the gut microbiome of infant rhesus macaques closely resembles that of children in the developing world where diarrheal disease and growth faltering is most prevalent29.
In this study, we leveraged the infant macaque model to characterize the impact of growth faltering with distant or active chronic diarrhea on the mucosal immune system and microbial communities throughout the intestinal tract. Our data revealed immunological and microbial dysregulation that was more pronounced in the colon. Specifically, we detected alterations in the frequency of immune cell populations and their response to polyclonal stimulation, as well as shifts in microbial communities along the entire GI tract, as well as the migration of microbes between gut sites suggestive of decompartmentalization. Shotgun metagenomics of luminal contents revealed key differences in the functional potential of luminal communities and led to the assembly of novel Campylobacter strains in animals with growth faltering. Together, our findings suggest that the pathology of growth faltering extends beyond the small intestine and may be more prominent in the large intestine.
Results
Cohort description and systemic markers of growth faltering and diarrhea
In this study, a total of 17 6–12-month old infant macaques were classified into three groups; Healthy growth (HG), Growth Faltering with acute distance or no history of diarrhea (GF), and Growth Faltering with active chronic diarrhea (GF-DX) based on growth rate and clinical status at the time of necropsy (Fig. 1a, Supp. Table 1). Briefly, HG (n = 7) grew at or above the average growth trajectory with only 1 of the 7 infants experiencing diarrhea ~3.3 months before necropsy (Supp. Table 1). In the GF (n = 5) group, 3 of 5 of the infants had experienced clinical diarrhea but had not had an episode of diarrhea for at least one 1-month before necropsy (Fig. 1a, Supp. Table 1). GF-DX (n = 5) all had a history of active chronic diarrhea and had diarrhea at the time of necropsy (Fig. 1a, Supp. Table 1). The clinical history of infants that experienced diarrhea was previously described29.
We next measured the levels of several systemic markers associated with barrier dysfunction and inflammation. We observed significantly higher levels of circulating soluble CD14 in GF-DX infants, with no differences between the GF and HG groups (Fig. 1b). No differences in circulating free or IgM bound LPS, or internal fatty acid-binding protein (iFABP) were observed between the three groups (Supp. Fig. 1a–c). Circulating levels of cytokines, chemokines, and growth factors were measured using multiplex ELISA (Luminex). A principal component analysis (PCA) revealed that the circulating immune mediator profile of both GF and GF-DX infants is distinct from that of the HG group (Fig. 1c). Levels of pro-inflammatory cytokines (Eotaxin, IFNα, IFNγ, IL-1β, IL-2, IL-15, and IL-12p70), as well as chemokines (BLC, IP-10, and PDF-2), were higher in GF compared to HG infants (Fig. 1d). On the other hand, levels of BDNF, MIP-1b, and VEGF-A were higher in GF-DX infants compared to HG infants (Fig. 1d). In addition, the levels of IL-8, IL-23, TNFα, VEGF-d were significantly different between the three groups by Kruskal Wallis non-parametric ANOVA, but not significantly elevated in a single group.
Growth faltering is associated with altered immune cell frequencies
We next measured the frequency of several immune subsets in ileal and colonic LPLs using flow cytometry. The overall frequencies of T or B-lymphocytes within colonic LPLs did not differ significantly between the three groups (Fig. 2a). However, the frequency of colonic CD4+ effector memory (EM) T cells was significantly increased while that of CD4+ transitional effector memory (TEM) T cells was decreased in GF-DX animals (Fig. 2b) indicative of accelerated differentiation of colonic CD4 T cells in GF-DX infants (Fig. 2b). No differences were detected in the relative frequencies of naïve or memory colonic CD8+ T cell subsets (Fig. 2c). In the ileum, we observed a significant decrease in the abundance of total CD4+ T cells in GF-DX infants (Supp. Fig. 1d), but no differences were observed in ileal CD8+ T cells or B cells (Supp. Fig. 1d). No differences were observed in the frequency of ileal naïve and memory CD4 or CD8 subsets (Supp. Fig. 1e, f).
No differences were detected in the frequency of dendritic cells (DCs), macrophages, or NK cells at either gut site (Fig. 2d, Supp. Fig. 1g). The frequency of colonic myeloid dendritic cells (mDC) was decreased in GF animals and trended lower in GF-DX animals (Fig. 2e), while the frequency of activated colonic CD16+ macrophages increased in GF animals (Fig. 2f). No differences in the frequencies of these subsets were detected between the three groups in the ileum (Supp. Fig. 1h, i).
Growth faltering with chronic active diarrhea results in dampened colonic LPL response to stimulation
To explore the impact of growth faltering/diarrhea on the function of LPLs, cells were stimulated with PMA/ionomycin (PMAi) for 16 h and the production of 33 immune mediators was measured using Luminex. No significant differences in the production of immune mediators by ileal LPLs were detected between the three groups of infants with or without PMAi stimulation as illustrated by the overlap of samples in the PCA plot (Supp. Fig. 2a) and the absolute amounts measured in supernatants (Supp. Fig. 2b). However, PMAi stimulation resulted in the increased production of CD-40L and CCL5 in only the HG group (Supp. Fig. 2b). In contrast, immune mediator production by colonic LPLs in response to PMAi was quite distinct between the three groups (Fig. 3a). Specifically, while colonic LPLs from healthy and GF animals formed distinct clusters post-stimulation, those from GF-DX animals overlapped with the unstimulated samples (Fig. 3a). Indeed, the concentration of pro-inflammatory cytokines (e.g., TNFα, IL-17, IFNγ), T cell-specific factors (e.g., IL-2, IL-4, IL5, CD-40L), and the growth factor GM-CSF were significantly reduced compared to those observed from the HG (Fig. 3b). To determine if this defect in immune activation was limited to the gut, PBMC were stimulated with PMAi for 6 h, and cytokine production measured by intracellular cytokine staining. Due to the paucity of available samples, the GF and GF-DX groups were combined for this experiment. As described for LPLs, the frequency of CD4 T cells secreting IL-17 and TNFα was significantly reduced in GF/GF-DX animals (Fig. 3c, d).
We further explored the impact of growth faltering/diarrhea on colonic LPLs response to PMAi using RNA-seq. PCA shows that the transcriptional profile of the 3-groups was distinct regardless of stimulation and that LPLs from GF-DX animals are distinct from those from HG and GF animals (Fig. 4a). Pairwise differential gene expression (DEG) showed a robust transcriptional response to PMAi stimulation with 387 and 1274 DEGs in HG and GF infants, respectively (Fig. 4b). In contrast, only 8 DEGs were detected in colonic LPLs from GF-DX animals in response to stimulation (Fig. 4b). Functional enrichment using Metascape revealed that PMAi-induced DEGs in HG and GF LPLs enriched to gene ontology (GO) terms associated with a response to immune stimulation such as “cell activation involved in an immune response”, “regulation of cytokine production” and “negative regulation of immune system processes” (Fig. 4c). In agreement with our cytokine measurements, IL-17A, CD40LG, IL2, and IFNg were highly expressed in colonic LPLs from HG animals and to a lesser extent from GF animals after stimulation (Fig. 4d). Moreover, IL10 and CD74 were upregulated in response to stimulation in HG and GF infants but downregulated in infants with acute diarrhea (Fig. 4d).
Gut transcriptional response to GF and GF-DX
To elucidate the effects of growth faltering/diarrhea on gut tissue homeostasis, we compared the transcriptomes of ileal and colonic tissue biopsies obtained from the three groups. PCA showed that the gut site exerted the most significant impact on the overall transcriptional profiles (Fig. 5a). DEG analysis between groups within colonic biopsies resulted in the identification of 87 DEGs between GF-DX and HG infants, 76 DEGs between GF-DX and GF infants, and only 3 DEGs between GF and HG infants (Fig. 5b, Supp. Table 2).
Functional enrichment using MetaScape showed that colonic DEGs unique to GF-DX vs. HG comparison enriched to GO terms associated with inflammation such as “regulation of inflammatory response”, “cytokine secretion”, and “lymphocyte-mediated immunity” (Fig. 5d). DEGs that enriched these GO terms were highly expressed in the colon of GF-DX infants (Fig. 5e) including TLR2, TLR8, LYZ, IDO1, CXCL10, and CCL5 (Fig. 5e). DEGs unique to the GF-DX vs. GF comparison enriched to the GO-term “release of Cytochrome C from the mitochondria” a key step in apoptotic cell death (Fig. 5d). Finally, we used IMMQUANT to infer cell populations based on transcriptional profiles. In agreement with our Flow cytometry experiment, this analysis predicted an increase in macrophages in GF infants, as well as a decrease in the frequency of early B cells and pro-B cells within the colons of GF and GF-DX infants, respectively relative to HG infants (Supp. Fig. 2c).
In the ileum, only 2 DEGs were identified between GF-DX and HG, 78 DEGs between GF-DX and GF infants, and 16 DEGs between GF and HG infants with 7 common DEGs (Fig. 5c). DEGs unique to GF-DX vs. GF infants enriched to immunological GO terms including “lymphocyte activation” and “cytokine-mediated signaling pathway” (Fig. 5d). DEGs enriching to these GO terms were highly expressed in the ileal biopsies of GF infants and included; IL10RA, ICOS, CD27, CXCR5, and CXCL13 (Fig. 5f). IMMQUANT analysis of ileal expression profiles predicts a decrease in effector memory CD8+ T cells in GF infants along with a decrease in granulocytes and an increase in plasmacytoid dendritic cells in GF-DX infants (Supp. Fig. 2d).
Since our RNA seq data suggested a significant increase in tissue inflammation in the colon with minimal changes in the ileum, we next evaluated inflammatory changes through histopathology. Examination of H&E stained slides from the mi-ileum revealed no gross differences between all three groups (Fig. 5g). Limited villous blunting and uniform intraepithelial lymphocyte infiltration was noted in all three groups (Fig. 5g). In contrast, analysis of the descending colon revealed significant damage in both the GF and GF-DX animal including crypt abscessations, inflammation in the submucosa, crypt hyperplasia, and neutrophil infiltration, all of which were most pronounced in the GF-DX animal (Fig. 5g).
Growth faltering and diarrhea are associated with dysbiosis and loss of compartmentalization of the intestinal microbiome
We next interrogated the association of growth faltering/diarrhea with the taxonomic landscape of the luminal microbiome across multiple gut sites (duodenum, jejunum, ileum, traverse colon, descending colon) and fecal samples obtained from the same HG, GF, and GF-DX infants using 16S rRNA gene amplicon sequencing. Gut site, rather than host status or individual, was the major determinant of community composition for weighted UniFrac dissimilarity (54% of total variation) with microbial communities of the small intestine (duodenum, jejunum, and ileum) clustering separately from those of the colon and fecal materials (Supp. Fig. 3a, b). Nevertheless, we observed an increase in phylogenetic diversity in jejunal and ileal communities from GF-DX infants compared to those from HG infants (Supp. Fig. 3c). Regardless of host status, the luminal microbiome of the small intestine was dominated by Streptococcus, Pasteurellaceae, Actinobacillus, and Veillonella; while Prevotella was the most abundant taxa in the large intestine and fecal samples (Fig. 6a).
Next, we determined which bacterial taxa were differentially abundant between healthy and faltering infants at each gut site. Only amplicon sequence variants (ASVs) detected in at least two samples at greater than 0.01% relative abundance were included in this analysis (Fig. 6b, Supp. Table 3). A total of 75 ASVs were differentially abundant across all gut sites (1 Duodenum, 8 Jejunum, 17 Ileum, 23 Transverse Colon, 13 Descending Colon, 13 Fecal) (Supp. Table 3). Notable differences include a higher abundance of multiple Pasteurellaceae, Lactobacillus, and Veillonella ASVs in the small intestine luminal contents from HG infants (Fig. 6b, Supp. Table 3). Similarly, small intestine luminal contents of GF infants were enriched in multiple Pasteurellaceae ASVs, while those from GF-DX infants had a higher abundance of ASVs assigned to Prevotella, Catenibacterium, and Lactobacillus salivarius (Fig. 6b, Supp. Table 3). The large intestine and fecal microbiome of HG infants were enriched in multiple Prevotella and Ruminococcaceae ASVs, as well as Faecalibacterium and Bracyspira (Fig. 6b, Supp. Table 3). On the other hand, the colonic microbiome of GF-DX infants was enriched in Lactobacillus salivarius along with multiple different Prevotella and Alloprevotella ASVs (Fig. 6b, Supp. Table 3). Finally, the fecal microbiome of GF infants was enriched in ASVs assigned to Campylobacter hyointestinalis, and Escherichia-Shigella (Fig. 6b, Supp. Table 3).
Growth faltering and diarrheal diseases are hypothesized to lead to a decompartmentalization of the gut microbiome; however, this pattern has only been inferred and has not been observed directly using paired samples. Fecal samples collected from GF/GF-DX animals had significantly lower observed ASVs and more intra-group community variability at all sights except the duodenum (Supp. Fig. 3d,e). We next plotted the abundance of taxa enriched in the luminal microbiomes of GF/GF-DX compared to HG infants along the GI tract (Fig. 6c–e, Supp. Fig. 3f–h, Supp. Table 4). This analysis showed that Streptococcus ASV 124, which is normally relegated to the small intestine luminal microbiome in HG animals, was abundant in the large intestine and fecal microbiome of GF/GF-DX (Fig. 6c). Similarly, the relative abundance of several ASVs normally relegated to the large intestine in HG animals (Prevotella ASV 147 and 23, Cantenibacterium ASV 156, and Lachnospiraceae ASV 24) was higher in the small intestine of GF/GF-DX infants (Fig. 6d,e, Supp. Fig. 3g,h). Finally, Alloprevotella ASV 336 was found almost exclusively in the lumen of GF/GF-DX infants and was significantly more abundant in both colon sections (Supp. Fig. 3f).
We next explored associations between immune dysregulation and microbial dysbiosis in GF and GF-DX animals. In the traverse colon, we observed 25 significant (p < 0.05) correlations between immune cell populations and the luminal microbiome (Fig. 6f). These correlations included positive correlations between Prevotella and CD4+ TEM T cells, Faecalibacterium, and total CD20+ B cells, along with Campylobacter hyointestinalis and CD4+ EM T cells (Fig. 6f). On the other hand, Pasteurellaceae and Ruminococcaceae were negatively correlated with total DC and CD8+ T cells, respectively (Fig. 6f). In the Ileum, we observed 8 significant correlations between local immune cell populations and the microbiome (Supp. Fig. 3i). These included positive correlations between total Macrophages and Lactobacillus, and a negative correlation between Alloprevotella and total CD8+ T cells (Supp. Fig. 3i).
Growth faltering and diarrhea are associated with functional shifts in the luminal microbiome
To enhance the resolution at which microbial communities are defined, we utilized shotgun metagenomics. More reads mapped to the host genome in the ileal luminal contents compared to colonic luminal contents due to the lower microbial biomass found in the small intestine (Supp. Fig. 4a, Supp Table 5). At the Phyla level, colonic luminal contents had significantly more Bacteroidetes, Spirochetes, Actinobacteria, and Euryarchaeota than ileal luminal contents (Supp. Fig. 4b, Supp Table 6). Reads that mapped to microbial taxa or functional genes in the ileum were distinct from those obtained from the colon regardless of host health status (Supp. Fig. 4c, d). A total of 46 species were differentially abundant between the ileum and the colon (Supp Table 6). Notably, Lactobacillus amaylovorus and Lactobacillus mucosae were enriched in the ileum while Prevotella copri, Helicobacter macacae, and Treponema succinifaciens were enriched in the colon (Supp Table 6).
Next, we compared the abundance of taxa between groups. The ileal luminal contents from HG infants were enriched in Lactobacillus reuteri while those of GF-DX infants had a higher abundance of Lactobacillus salivarius (Supp. Fig. 4e). Larger differences between the three groups were observed in the colonic luminal contents. Microbial communities from HG infants had a higher abundance of the commensals Helicobacter macacae and Ruminococcus obeum (Supp. Fig. 4g). In contrast, microbial communities of GF and GF-DX are enriched in potential pathobionts notably Fusobacterium nucleatum and Helicobacter cinaedi in GF infants (Supp. Fig. 4g) and Lactobacillus salivarius, Escherichia coli, and Enterococcus faecalis in GF-DX infants (Supp. Fig. 4g).
Although the overall functional potential of the ileal contents did not significantly differ between the three groups as measured by Bray Curtis dissimilarity (permanova p > 0.05) (Supp. Fig. 4d), the ileal microbiome of HG infants was enriched in the pathway for stachyose degradation, folate transformation and galactose degradation (Supp. Fig. 4f). The ileal microbiome of GF infants was enriched in pathways for L-histidine degradation III and D-beta-fucofuranose biosynthesis, while that of GF-DX infants were enriched for L-methionine and L-ornithine biosynthesis (Supp. Fig. 4f). Similarly, the overall functional potential of colonic luminal contents did not differ significantly between HG and GF infants as measured by Bray Curtis dissimilarity (permanova p > 0.05) (Supp. Fig. 4d). Nevertheless, the colonic microbiome of HG infants was enriched for pyruvate fermentation to butanoate, pyruvate fermentation to hexanol, and glutaryl-CoA degradation pathways (Supp. Fig. 4h), while those of GF infants were enriched in the citrulline biosynthesis pathway, and those of infants with GF-DX in pathways for thiazole biosynthesis, fatty acid degradation, and ornithine degradation (Supp. Fig. 4h).
We were able to assemble 23 putative microbial genomes from colonic luminal contents with >80% completeness and <2% contamination as measured by checkM39. including 9 Prevotella, 6 Helicobacter, and 3 Campylobacter genomes (Supp. Table 7). The three assembled Campylobacter genomes were distinct from each other and most closely related to Campylobacter garcilis, Campylobacter coli/jejuni, and the recently proposed species, Campylobacter infanis25. (Fig. 7a). Importantly, all Campylobacter genomes were only assembled from GF and GF-DX infants, and significantly more reads mapped to these assembled genomes in samples from GF and GF-DX infants compared to HG infants (Fig. 7b). Helicobacter genomes were assembled from all three groups (Fig. 7c) and no differences were noted in the percentage of reads mapping to these genomes between the three groups (Fig. 7d). In addition, we assembled a diverse set of Prevotella genomes primarily from GF infants (Fig. 7e); however, no differences in the abundance of reads that mapped to these genomes were observed between the groups (Fig. 7f).
Discussion
In addition to high mortality, diarrheal diseases are a major driver of malnutrition, as well as growth and cognitive stunting4,6,13,40,41. One of the major knowledge gaps in the field is the lack of studies pertaining to host-microbial interactions within the gut. These studies are difficult to conduct in humans due to challenges associated with obtaining samples along the entire GI tract. Rhesus macaques are naturally susceptible to human enteric pathogens that result in diarrhea and growth faltering36,37. In this study, we leveraged this translational model to elucidate the immunological and microbial changes associated with growth faltering and a history of diarrhea throughout the GI tract in addition to fecal samples.
In contrast to what has been reported in humans42,43 we did not observe changes in systemic levels of molecules associated with barrier dysfunction (iFABP) or bacterial translocation (LAL, IgM bound endotoxin) in GF and GF/DX infants. However, our analysis revealed higher levels of systemic inflammatory markers in both groups of growth faltering infants. In addition, we found a stark increase in the levels of BDNF which was detectable in 2/6 HG infants, 3/4 GF infants, and 4/4 GF-DX infants. BDNF can be produced by multiple cell types in the gut44,45 and is important for the regulation of tight junctions46. Moreover, systemic levels of BDNF have been positively correlated with visceral hypersensitivity (pain) in irritable bowel syndrome (IBS)45,47. Our data suggest that plasma levels of BDNF could be used as an alternative marker of intestinal injury. Interestingly, plasma levels of IL-1β, IFNα, IFNγ, and IL-15 were higher in GF infants, indicating that distant history of diarrhea can have long-lasting ramifications on the host immune status. High levels of systemic IL-1β have been linked to poor neurodevelopment in Bangladeshi children48, while IFNγ has been associated with increased microbial translocation in growth faltering infants49. This enhanced systemic inflammation was accompanied by a dampened IL-17 and TNFα response of peripheral CD4 T cells to polyclonal stimulation, potentially an indication of an immunoregulatory phenotype.
Significant differences in immune cell frequencies and response to polyclonal stimulation between HG and GF-DX animals were observed in colonic, but not ileal, LPLs. Specifically, we observed a reduction in the relative frequency of colonic mDCs, a major antigen-presenting cell subset50, in both faltering groups. This loss could be one of the factors contributing to the reduced immune response to oral Rotavirus and Polio vaccines observed in children with growth faltering and diarrheal diseases51,52. In addition, we detected a higher relative frequency of the non-classical colonic CD16+ macrophages in GF infants, indicative of higher inflammation even in the absence of active diarrhea. Indeed, CD16+ macrophages in the gut have been shown to mediate fibrosis in inflammatory bowel disease53. We also observed an increase in the percentage of terminally differentiated colonic memory CD4+ T, suggesting an accelerated differentiation. CD4+ T cells are the most abundant LPL population and play important roles in tissue maintenance and anti-microbial responses54,55,56. In healthy infants, tissue-resident T cells are typically naive, and regulatory and resident effector memory cells do not develop until adulthood57,58. High pathogen load or repeated exposure to antigen could drive the observed shifts in memory CD4 T cell populations. Terminally differentiated memory T cells are usually associated with increased inflammatory responses. However, colonic LPLs from GF-DX animals generated a dampened response to polyclonal stimulation relative to colonic LPLs from HG infants failing to produce canonical cytokines IL-2, IL-17, and IFNγ. This dampened response was also evident at the transcriptional level. Interestingly colonic LPLs from GF infants generated a larger transcriptional response to polyclonal stimulation in line with increased plasma levels of inflammatory cytokines. Together these observations suggest that GF with a limited diarrhea history is associated with exacerbated inflammation while GF with chronic diarrhea is associated with immune exhaustion. However, further measurement of exhaustion markers on LPL such as PD-1 and CTLA-4 in the colon is needed to confirm this.
Our transcriptomic analysis of ileal and colonic biopsies revealed the downregulation of genes important for immune activation and highly expressed by gut-associated B cells59 in GF-DX compared to GF infants. This suggests that GF in the absence of active diarrhea is associated with increased B cell activity in the ileum. B cells have been shown to increase in abundance and aggregate in the small intestine of growth faltering infants in the Gambia49. This chronic immune activation may contribute to systemic inflammation and inhibit the host from adequately responding to mucosal vaccines/infection. In contrast, the transcriptional landscape in the colon was indicative of enhanced immune activation in GF-DX infants compared to GF infants. Many of the genes upregulated in the colon of GF-DX animals were associated with T cell and macrophage activation in agreement with our flow cytometry data, as well as granulocytes (not included in our immunophenotyping). Finally, histological examination of tissue biopsies from the mid-ileum and traverse colon confirmed that GF is associated with a heightened inflammatory landscape in the colon. Taken together, our phenotypic, functional, transcriptional, and histological data strongly suggest that GF, especially when coupled with chronic active diarrhea, primarily impacts the colon rather than the small intestine. While more recent human studies have focused on the small intestine, earlier studies identified mucosal abnormalities in rectal biopsies of volunteers in southern India which they termed tropical colonopathy60.
This study provided the first insight into the microbial communities within luminal contents along the intestinal tract of infant macaques. As previously reported for adult macaques61, the luminal microbiome of the small and large intestines are distinct. Phylogenetic diversity in the small intestinal was higher in GF-DX infants. On the other hand, and as described in human infants with, malnutrition32 and growth faltering31, the alpha diversity in the fecal samples of GF/GF-DX macaques was reduced compared to HG macaques. Furthermore, inter-sample variability was higher in microbial communities from GF/GF-DX animals, indicative of microbiome instability or a more immature state as has been shown in the fecal microbiome of malnourished children32.
In addition, multiple bacterial taxa normally sequestered to either the small or large intestine in HG infants were detected in several gut sites in GF/GF-DX infants, suggesting the loss of strict compartmentalization between gut-sites with growth faltering. Specifically, the relative abundance of Lactobacillus salivarius, a core member of the oral microbiome in humans62, were increased in the lumen of both the small and large intestine of GF-DX infants. L. salivarius has been reported to be enriched in the fecal samples of growth-stunted children33. Moreover, Streptococcus, also part of the core oral microbiome62 and normally relegated to the small intestine in healthy animals, was detected in colonic luminal contents. This observation is similar to recent reports describing the detection of nasopharyngeal taxa to the gut in children experiencing growth stunting33. This loss of strict compartmentalization could be due to a reduced ability of the host immune system to tightly control microbial growth.
We also observed differences in the functional potential in both the ileal and colonic luminal microbiome. Specifically, we observed a reduced capacity of short-chain fatty acid (SCFA) production in GF and GF-DX infants in both ileal and colonic microbial communities. The microbial production of SCFA has been closely linked to the activation of mucosal Th1, Th17, and Tregs63. However, SCFA has also been reported to skew T cells towards anti-inflammatory IL-10 production and play a role in preventing colitis and other inflammatory conditions64. In addition, the SCFA butyrate can modulate the function of innate immune cells by reducing the expression and nuclear translocation of NF-κB in intestinal macrophages65,66. Therefore, the loss of SCFA production potential in GF and GF-DX animals may contribute to the heightened inflammatory state we observed.
We also observed correlations between intestinal immune cell frequencies and the relative abundance of microbes. In the Ileum this included a positive correlation between Lactobacillus abundance and Ileal Macrophages. Lactobacillus strains have also been shown to activate macrophages in the small intestine67,68. The increased abundance of Lactobacillus salivarius in the small and large intestine of GF-DX infants coupled with increased frequency of activated macrophages suggest an infiltration of oral microbes into the small intestine resulting in activation of the innate immune system. In the colon, we found an inverse relationship between Prevotella and Campylobacter hyointestinalis with colonic CD4+ T cell subsets. While gut Prevotella has been linked to systemic inflammation and reduced systemic CD4+ T cells in the western world, particularly in HIV infection69,70,71. It is important to note that Prevotella is the most abundant microbe in the gut of humans with a non-westernized lifestyle and is largely believed to be beneficial in this context72,73. However, the positive correlation of CD4+ EM T cells with Campylobacter hyointestinalis may indicate recurrent infection and the expansion of a pathogen-specific T cell subset.
We detected an increased abundance of Campylobacter hyointestinalis and Escherichia-Shigella in fecal samples collected from GF infants, indicative of an increased pathogen burden, as reported in multiple human studies19,74,75. Diversity within the Campylobacter genus is often missed in 16S amplicon surveys due to their highly conserved 16S rRNA gene76; therefore, we resorted to metagenomic genome assembly. We constructed 3 Campylobacter genomes only from the colonic luminal contents of faltering infants. Moreover, a higher number of short reads mapped to these three genomes in colonic luminal contents of faltering infants. One of our assembled Campylobacter genomes, as well as some of our previously described29. Campylobacter assemblies from infant rhesus macaque feces were most closely related to the newly described “Candidatus Campylobacter infanis”25. Non-coli/jejuni Campylobacter has recently been implicated in diarrheal diseases and childhood mortality in the developing world25,26,77.
In contrast, we did not observe differences in the abundance of our assembled Helicobacter genomes. These data differ from our recent report of a decreased abundance of Helicobacter in the fecal microbiome of infants that experienced diarrhea and received antibiotics. This difference could be potentially be explained by the fact that Helicobacter is tightly associated with gut mucosa and greatly reduced in the lumen61. We also assembled several new Prevotella genomes mostly from GF animals. These new genomes differed from the cluster of Prevotella genomes we recently identified from the fecal microbiome of 1-month old infants that would later experience diarrhea29. Interestingly, the frequency of reads mapping to these genomes did not differ between the three groups even though 90% of the genomes were assembled from GF/GF-DX animals.
In summary, this study provides the first comprehensive look at the association of growth faltering and diarrheal disease with luminal microbiomes and local immune dysfunction that would not be possible in humans. Together, our data suggest that the colon is the major site of immune dysfunction during active diarrhea. We also found that immune dysfunction in GF and GF-DX is associated with a decompartmentalization of the luminal microbiome. Data presented in this manuscript highlight the need for novel multifaceted approaches that can address the complex etiology of GF and diarrheal disease with the goal of reducing pathogen burden, restoring gut integrity, and rejuvenating the intestinal immune system. Previous interventions aimed at reducing childhood mortality due to growth faltering and enteric disease78 including the prophylactic administration of azithromycin77,79, micronutrients80, probiotics (Lactobacillus species)81. However, many of these interventions have resulted in little to no improvement in growth or other primary endpoints78. Potential new strategies could include broad-acting phage therapy against Campylobacter, probiotic treatment with commensals that are relevant to children who live in the developing world, SCFA supplementation to reduce inflammation and improve gut immune fitness, and immune-modulatory drugs that can improve gut integrity. Infant macaques provide a highly relevant translational animal model to test these novel interventions and evaluate their impact on tissue inflammation, immune function, and microbial communities across all sections of the gut.
Methods
Cohort description
Samples were obtained from a total of 17 6–12-month old infant macaques (Supp. Table 1). All rhesus macaque studies were overseen and approved by the OHSU/ONPRC and University of California-Davis/CNPRC Institutional Animal Care and Use Committees’ (IACUC) per the National Institutes of Health guide for the care and use of laboratory animals. Animals were housed per the standards established by the US Federal Animal Welfare Act and The Guide for the Care and Use of Laboratory Animals. All animals were tested for simian viruses (Simian Immunodeficiency Virus, Simian Retrovirus 2, Macacine herpesvirus 1, and Simian T lymphotropic virus) and received a tuberculin test semi-annually.
Measurement of systemic markers
Levels of Endotoxin-core antibodies IgM (Hycult Biotech, Uden, Netherlands), soluble CD14 (sCD14) (R&D Systems, Minneapolis, MN, United States), and internal fatty acid-binding protein (iFABP) (MyBioSource, San Diego, CA, United States) in plasma samples were measured using an enzyme-immunoassay technique (ELISA) following the manufacturer’s protocol. Plasma levels of free bacterial endotoxin were measured using a chromogenic Limulus amebocyte lysate assay (Hycult Biotech, Uden, Netherlands) following the manufacturer’s protocol. Plasma levels of circulating immune factors were measured using a Magnetic 33-plex assay (R&D Systems, Minneapolis, MN, United States) and a MAGPIX instrument (Luminex, Austin, TX, United States). Bubble plots were generated using ggplot2 R package with concentration (pg/ml).
Isolation of lamina propria leukocytes (LPL)
Tissue biopsies were collected from the ileum and descending colon into Roswell Park Memorial Institute medium supplemented with 10% fetal bovine serum, streptomycin/penicillin, and L-glutamine (RP10). Lamina propria leukocytes (LPLs) were isolated using collagenase IX/DNAse digestion followed by centrifugation over a Percoll gradient (GE Healthcare, Waukesha, WI, USA) as previously described82. Cells were cryopreserved in FetalPlex Animal Serum Complex (Gemini Bio-Products, West Sacramento, CA, USA) and DMSO.
Flow cytometry
1 × 106 LPLs were stained using antibodies against CD3, CD4, CD8b, CD95, CD28, CCR7, and CD20 to delineate naïve and memory CD4 and CD8 T cells, as well as CD20 B cells. The cells were subsequently fixed and permeabilized before the addition of anti-Ki-67 to measure proliferation. T cells were then divided into Naive (CD28+, CD95−; Naive), central memory (CCR7+ CD28+ CD95+; CM), transitional effector memory (CCR7− CD28+ CD95+; TEM), and Effector memory (CCR7− CD28− CD95+; EM), and. The second tube of 1 × 106 LPLs was stained using antibodies against: CD3, CD20, HLA-DR, CD14, CD11c, CD123, CD16, and CD8a to delineate macrophages (CD3− CD20− CD14 + HLA-DR+), dendritic cells (DC, CD3− CD20-CD14− HLA-DR+), and natural killer (NK, CD3-CD20-CD14− HLA-DR− CD8a+) cell subsets. Macrophages were further subdivided into classical (CD16−) and non-classical macrophages (CD16+). DCs were further subdivided into myeloid DC (mDC, CD123− CD11c+) and plasmacytoid DC (pDC, CD123+CD11c−;). All flow cytometry samples were acquired using Attune NxT (Life Technologies, Carlsbad, CA, United States) and analyzed using FlowJo (TreeStar, Ashland, OR, United States).
LPL stimulation
LPLs (5 × 105/well) were thawed and incubated for 16 h in the absence (unstimulated) or presence (stimulated) of 5 ng/ml PMA and 1 μg/ml Ionomycin (Sigma-Aldrich, St. Louis, MO, United States) at 37 °C in a humidified incubator (5% CO2). At the end of the incubation, cells were centrifuged for 5 min at 2000 rpm. Supernatants were collected to measure the production of immune mediators using a Magnetic 33-plex assay (R&D Systems, Minneapolis, MN, United States) and a MAGPIX instrument (Luminex, Austin, TX, United States). Absolute change of immune mediators was calculated by subtracting the basal level of production from the level post-stimulation in a pairwise fashion. Bubble plots were generated using ggplot2 R package with concentration (pg/ml).
RNA-sequencing and data analysis
Total RNA was isolated from ileum and colon biopsies, as well as ileal and colonic LPLs following overnight culture in the presence/absence of PMA/Ionomycin using Qiagen miRNeasy Kit. RNA concentration and quality were verified using Agilent 2100 Bioanalyzer. RNA-seq libraries were constructed using the TruSeq Stranded Total RNA Kit (Illumina, San Diego, CA, United States). rRNA-depleted RNA was fragmented and converted to double-stranded cDNA. Adapters were ligated, and ~300 base pair fragments were amplified by PCR and selected by size exclusion using Ampure beads XP (Beckman Coulter, Brea, CA). Each library was labeled with a unique barcode for multiplexing. To ensure proper sizing, quantitation, and quality before sequencing, libraries were analyzed on the Agilent 2100 Bioanalyzer. Multiplexed libraries were sequenced (100 base pair single-end) using the HiSeq 4000 platform.
Data analysis was performed with the RNA sequencing workflow module of the systemPiperR package available on Bioconductor83. Quality filtration and trimming were performed using Trimmomatic with an average phred score cutoff of 30 and a minimum length of 75 bp. Reads were mapped with the splice-aware aligner suite Bowtie2/Tophat2 against the Macaca mulatta genome (Macaca_mulatta.MMUL_1.78). Raw count values were generated with summarizeOverlaps function. Only reads overlapping exonic regions of genes were counted, discarding reads mapping to ambiguous regions of exons from overlapping genes. Analysis of differentially expressed genes (DEGs) was performed using edgeR package83. DEGs were defined as those with a fold change of ≥2 and a Benjamini-Hochberg–controlled false discovery rate of <0.05. Functional enrichment analysis was performed to identify significant Gene Ontology (GO) terms and using Metascape84. Immquant was used to infer alteration in immune cell frequency based on tissue transcriptional profiles85.
Intestinal pathology imaging
Tissue samples from each animal were collected from the mid ileum, and descending colon. The tissues were placed in 10% neutral buffered formalin at room temperature for 48–72 h and processed in a Tissue Tek VIP 5 Tissue Processor. The blocks were paraffin-embedded, cut, and stained with hematoxylin and eosin. Representative images were captured at ×10 magnification by a blinded pathologist.
16S amplicon sequencing
Total DNA was extracted from ~0.25 g of luminal contents using the PowerSoil DNA Isolation Kit (MO BIO Laboratories, Carlsbad, CA, USA). Next, we amplified the hypervariable V4-V5 region of the 16S rRNA gene using PCR primers (515 F/926 R with the forward primers including a 12-bp barcode)86. PCR reactions contained 12.5 μl GoTaq mastermix, 9.5 μl nuclease-free H20, 1 μl template DNA, and 1 μl 10 μM primer mix. Thermal cycling parameters were 94 °C for 5 min, 35 cycles of 94 °C for 20 s, 50 °C for 20 s, 72 °C for 30 s, followed by 72 °C for 5 min. PCR products were purified using a MinElute 96 UF PCR Purification Kit (Qiagen, Valencia, CA, USA). Libraries were sequenced (2 ×300 bases) using an Illumina MiSeq.
Raw FASTQ 16S rRNA gene amplicon sequences were uploaded and processed using the QIIME2 analysis pipeline87. Briefly, sequences were demultiplexed and the quality filtered using DADA288, which filters chimeric sequences and generates an amplicon sequence variant (ASV) table equivalent to an operational taxonomic unit (OTU) table at 100% sequence similarity. Sequence variants were then aligned using MAFFT89. and a phylogenetic tree was constructed using FastTree290. Taxonomy was assigned to sequence variants using q2-feature-classifier91. against the SILVA database (release 119)92. To prevent sequencing depth bias samples were rarified to 13,000 sequences per sample before alpha and beta diversity analysis. QIIME 2 was also used to generate the following alpha diversity metrics: richness (as observed taxonomic units), Shannon evenness, and phylogenetic diversity. Beta diversity was estimated in QIIME 2 using weighted and unweighted UniFrac distances93.
Shotgun metagenomic library preparation and analysis
Shotgun metagenomic libraries were prepared from 50 ng of gDNA using the Nextera library prep (Illumina, La Jolla CA) per Illumina’s recommended protocol and sequenced on an Illumina HiSeq 4000 2 × 100. Raw demultiplexed reads were quality filtered using Trimmomatic94, and potential host reads were removed by aligning trimmed reads to the Macaca mulata genome (Mmul 8.0.1) using BowTie295. Trimmed and decontaminated reads were then annotated using the HUMAnN296. pipeline using default settings with the UniRef50 database and assigned to Metacyc pathways. Functional annotations were normalized using copies per million (CPM) reads before statistical analysis97,98,99. Species-level taxonomy was assigned to quality-controlled short reads using Metaphlan2100. Genome assemblies were generated for each sample. Trimmed and decontaminated reads were assembled into contigs using meta-SPAdes with default parameters101. and binned into putative genomes using MetaBat102. Genome completeness/contamination was tested using CheckM39, and all bins with a completeness >80% and contamination <2% were annotated using PATRIC103. Taxonomy of draft genomes was determined using PATRICs’ similar genome finder.
Statistical analysis
All statistical analyses were conducted using PRISM (V8) or the R package Vegan104. Bray Curtis dissimilarity matrices were constructed for both species-level relative abundance and normalized gene annotations using the vegdist function in the R package Vegan. Principal coordinate analysis (PcoA) was conducted using the R function cmdscale. PERMANOVAs were performed using the Vegan function ADONIS. Unpaired t-test, 1-way, and 2-way ANOVA were implemented using PRISM8 were noted to generate p-values and utilizing the corresponding post-hoc-test when the initial ANOVA was significant. The LEfSe algorithm was used to identify differentially abundant taxa and pathways between groups with logarithmic Linear discriminant analysis (LDA) score cutoff of 299.
Spearman rank correlations were used to investigate the relationship between Lamina propria immune cell populations and the taxonomy of the luminal contents collected from the same site. The “rcorr” function from the R package Hmisc to calculate Spreaman rho and p-values105. These values were then plotted using the R package corrplot106.
Data availability
The sequencing datasets generated for this study are available in the NCBI SRA repository, under the bioproject ID no. PRJNA729051.
References
Liu, L. et al. Global, regional, and national causes of under-5 mortality in 2000-15: an updated systematic analysis with implications for the Sustainable Development Goals. Lancet 388, 3027–3035 (2016).
Clasen, T. F. et al. Interventions to improve water quality for preventing diarrhoea. Cochrane Database Syst. Rev. 10, CD004794 (2015).
Ferdous, F. et al. Severity of diarrhea and malnutrition among under five-year-old children in rural Bangladesh. Am. J. Trop. Med. Hyg. 89, 223–228 (2013).
Mata, L. Diarrheal disease as a cause of malnutrition. Am. J. Trop. Med. Hyg. 47, 16–27 (1992).
Baqui, A. H. et al. Malnutrition, cell-mediated immune deficiency, and diarrhea: a community-based longitudinal study in rural Bangladeshi children. Am. J. Epidemiol. 137, 355–365 (1993).
Yip, R. & Sharp, T. W. Acute malnutrition and high childhood mortality related to diarrhea. Lessons from the 1991 Kurdish refugee crisis. JAMA 270, 587–590 (1993).
Brown, K. H., Khatun, M. & Ahmed, G. Relationship of the xylose absorption status of children in Bangladesh to their absorption of macronutrients from local diets. Am. J. Clin. Nutr. 34, 1540–1547 (1981).
Manary, M. J. et al. Perturbed zinc homeostasis in rural 3-5-y-old Malawian children is associated with abnormalities in intestinal permeability attributed to tropical enteropathy. Pediatr. Res. 67, 671–675 (2010).
Campbell, D. I., Elia, M. & Lunn, P. G. Growth faltering in rural Gambian infants is associated with impaired small intestinal barrier function, leading to endotoxemia and systemic inflammation. J. Nutr. 133, 1332–1338 (2003).
Ngure, F. M. et al. Water, sanitation, and hygiene (WASH), environmental enteropathy, nutrition, and early child development: making the links. Ann. N. Y. Acad. Sci. 1308, 118–128 (2014).
John, C. C., Black, M. M. & Nelson, C. A. Neurodevelopment: the impact of nutrition and inflammation during early to middle childhood in low-resource settings. Pediatrics 139, S59–S71 (2017).
Czerkinsky, C. & Holmgren, J. Vaccines against enteric infections for the developing world. Philos. Trans. R. Soc. Lond. B 370, 20150142 (2015).
Moore, S. R. et al. Prolonged episodes of acute diarrhea reduce growth and increase risk of persistent diarrhea in children. Gastroenterology 139, 1156–1164 (2010).
Guerrant, R. L., Oria, R. B., Moore, S. R., Oria, M. O. & Lima, A. A. Malnutrition as an enteric infectious disease with long-term effects on child development. Nutr. Rev. 66, 487–505 (2008).
McCormick, B. J. J. & Lang, D. R. Diarrheal disease and enteric infections in LMIC communities: how big is the problem? Trop. Dis. Travel Med. Vaccines 2, 11 (2016).
Rouhani, S. et al. Gut microbiota features associated with Campylobacter burden and postnatal linear growth deficits in a Peruvian birth cohort. Clin. Infect. Dis. 71, 1000–1007 (2019).
Rogawski, E. T. et al. Use of quantitative molecular diagnostic methods to investigate the effect of enteropathogen infections on linear growth in children in low-resource settings: longitudinal analysis of results from the MAL-ED cohort study. Lancet Glob. Health 6, e1319–e1328 (2018).
Schaible, U. E. & Kaufmann, S. H. Malnutrition and infection: complex mechanisms and global impacts. PLoS Med. 4, e115 (2007).
Platts-Mills, J. A. et al. Pathogen-specific burdens of community diarrhoea in developing countries: a multisite birth cohort study (MAL-ED). Lancet Glob. Health 3, e564–e575 (2015).
Kotloff, K. L. et al. Burden and aetiology of diarrhoeal disease in infants and young children in developing countries (the Global Enteric Multicenter Study, GEMS): a prospective, case-control study. Lancet 382, 209–222 (2013).
Lee, G. et al. Symptomatic and asymptomatic Campylobacter infections associated with reduced growth in Peruvian children. PLoS Negl. Trop. Dis. 7, e2036 (2013).
Kotloff, K. L. et al. The Global Enteric Multicenter Study (GEMS) of diarrheal disease in infants and young children in developing countries: epidemiologic and clinical methods of the case/control study. Clin. Infect. Dis. 55, S232–S245 (2012).
Peeling, R. W., Smith, P. G. & Bossuyt, P. M. A guide for diagnostic evaluations. Nat. Rev. Microbiol. 4, S2–S6 (2006).
Platts-Mills, J. A. et al. Detection of Campylobacter in stool and determination of significance by culture, enzyme immunoassay, and PCR in developing countries. J. Clin. Microbiol. 52, 1074–1080 (2014).
Bian, X. et al. Campylobacter abundance in breastfed infants and identification of a new species in the global enterics multicenter study. mSphere 5, e00735-19 (2020).
Francois, R. et al. The other Campylobacters: not innocent bystanders in endemic diarrhea and dysentery in children in low-income settings. PLoS Negl. Trop. Dis. 12, e0006200 (2018).
Yu, J. et al. Environmental enteric dysfunction includes a broad spectrum of inflammatory responses and epithelial repair processes. Cell Mol. Gastroenterol. Hepatol. 2, 158–174 e151 (2016).
Rouhani, S. et al. Diarrhea as a potential cause and consequence of reduced gut microbial diversity among undernourished children in Peru. Clin. Infect. Dis. 71, 989–999 (2020).
Rhoades, N. et al. Maturation of the infant rhesus macaque gut microbiome and its role in the development of diarrheal disease. Genome Biol. 20, 173 (2019).
Pop, M. et al. Diarrhea in young children from low-income countries leads to large-scale alterations in intestinal microbiota composition. Genome Biol. 15, R76 (2014).
Gough, E. K. et al. Linear growth faltering in infants is associated with Acidaminococcus sp. and community-level changes in the gut microbiota. Microbiome 3, 24 (2015).
Subramanian, S. et al. Persistent gut microbiota immaturity in malnourished Bangladeshi children. Nature 510, 417–421 (2014).
Vonaesch, P. et al. Stunted childhood growth is associated with decompartmentalization of the gastrointestinal tract and overgrowth of oropharyngeal taxa. Proc. Natl Acad. Sci. USA 115, E8489–E8498 (2018).
Rouhani, S. et al. Diarrhea as a potential cause and consequence of reduced gut microbial diversity among undernourished children in Peru. Clin. Infect. Dis. 71, 989–999 (2019).
Gallardo, P. et al. Distinctive gut microbiota is associated with diarrheagenic Escherichia coli infections in Chilean children. Front. Cell Infect. Microbiol. 7, 424 (2017).
Prongay, K., Park, B. & Murphy, S. J. Risk factor analysis may provide clues to diarrhea prevention in outdoor-housed rhesus macaques (Macaca mulatta). Am. J. Primatol. 75, 872–882 (2013).
Haertel, A. J., Prongay, K., Gao, L., Gottlieb, D. H. & Park, B. Standard growth and diarrhea-associated growth faltering in captive infant rhesus macaques (Macaca mulatta). Am. J. Primatol. 80, e22923 (2018).
Quintel, B. K. et al. Vaccine-mediated protection against Campylobacter-associated enteric disease. Sci. Adv. 6, eaba4511 (2020).
Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P. & Tyson, G. W. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055 (2015).
Briend, A. Is diarrhoea a major cause of malnutrition among the under-fives in developing countries? A review of available evidence. Eur. J. Clin. Nutr. 44, 611–628 (1990).
Checkley, W. et al. Multi-country analysis of the effects of diarrhoea on childhood stunting. Int. J. Epidemiol. 37, 816–830 (2008).
Zhang, R. X. et al. Primary tumor location as a predictor of the benefit of palliative resection for colorectal cancer with unresectable metastasis. World J. Surg. Oncol. 15, 138 (2017).
Kaonga, P. et al. Direct biomarkers of microbial translocation correlate with immune activation in adult Zambians with environmental enteropathy and Hepatosplenic Schistosomiasis. Am. J. Trop. Med Hyg. 97, 1603–1610 (2017).
Lommatzsch, M. et al. Abundant production of brain-derived neurotrophic factor by adult visceral epithelia. Implications for paracrine and target-derived Neurotrophic functions. Am. J. Pathol. 155, 1183–1193 (1999).
Wang, P. et al. BDNF contributes to IBS-like colonic hypersensitivity via activating the enteroglia-nerve unit. Sci. Rep. 6, 20320 (2016).
Yu, Y. B. et al. BDNF modulates intestinal barrier integrity through regulating the expression of tight junction proteins. Neurogastroenterol. Motil 29, e12967 (2017).
Joo, Y. E. Increased expression of brain-derived neurotrophic factor in irritable bowel syndrome and its correlation with abdominal pain (Gut 2012;61:685-694). J. Neurogastroenterol. Motil. 19, 109–111 (2013).
Jiang, N. M. et al. Febrile illness and pro-inflammatory cytokines are associated with lower neurodevelopmental scores in Bangladeshi infants living in poverty. BMC Pediatr. 14, 50 (2014).
Marie, C., Ali, A., Chandwe, K., Petri, W. A. Jr & Kelly, P. Pathophysiology of environmental enteric dysfunction and its impact on oral vaccine efficacy. Mucosal Immunol. 11, 1290–1298 (2018).
Bell, S. J. et al. Migration and maturation of human colonic dendritic cells. J. Immunol. 166, 4958–4967 (2001).
Taniuchi, M. et al. Impact of enterovirus and other enteric pathogens on oral polio and rotavirus vaccine performance in Bangladeshi infants. Vaccine 34, 3068–3075 (2016).
Naylor, C. et al. Environmental enteropathy, oral vaccine failure and growth faltering in infants in Bangladesh. EBioMedicine 2, 1759–1766 (2015).
Salvador, P. et al. CD16+ macrophages mediate fibrosis in inflammatory bowel disease. J. Crohns Colitis 12, 589–599 (2018).
Tauschmann, M. et al. Distribution of CD4(pos) -, CD8(pos) - and regulatory T cells in the upper and lower gastrointestinal tract in healthy young subjects. PLoS ONE 8, e80362 (2013).
Gad, M., Brimnes, J. & Claesson, M. H. CD4+ T regulatory cells from the colonic lamina propria of normal mice inhibit proliferation of enterobacteria-reactive, disease-inducing Th1-cells from SCID mice with colitis. Clin. Exp. Immunol. 131, 34–40 (2003).
Yu, X. et al. Intestinal Lamina Propria CD4(+) T cells promote bactericidal activity of macrophages via Galectin-9 and Tim-3 Interaction during Salmonella enterica Serovar Typhimurium Infection. Infect. Immun. 86, e00769-17 (2018).
Thome, J. J. et al. Early-life compartmentalization of human T cell differentiation and regulatory function in mucosal and lymphoid tissues. Nat. Med. 22, 72–77 (2016).
Zens, K. D., Connors, T. & Farber, D. L. Tissue compartmentalization of T cell responses during early life. Semin. Immunopathol. 39, 593–604 (2017).
Immunological Genome, P. ImmGen at 15. Nat. Immunol. 21, 700–703 (2020).
Mathan, M. M. & Mathan, V. I. Rectal mucosal morphologic abnormalities in normal subjects in southern India: a tropical colonopathy? Gut 26, 710–717 (1985).
Yasuda, K. et al. Biogeography of the intestinal mucosal and lumenal microbiome in the rhesus macaque. Cell Host Microbe 17, 385–391 (2015).
Griffen, A. L. et al. CORE: a phylogenetically-curated 16S rDNA database of the core oral microbiome. PLoS ONE 6, e19051 (2011).
Park, J. et al. Short-chain fatty acids induce both effector and regulatory T cells by suppression of histone deacetylases and regulation of the mTOR-S6K pathway. Mucosal Immunol. 8, 80–93 (2015).
Sun, M. et al. Microbiota-derived short-chain fatty acids promote Th1 cell IL-10 production to maintain intestinal homeostasis. Nat. Commun. 9, 3555 (2018).
Luhrs, H. et al. Butyrate inhibits NF-kappaB activation in lamina propria macrophages of patients with ulcerative colitis. Scand. J. Gastroenterol. 37, 458–466 (2002).
Segain, J. P. et al. Butyrate inhibits inflammatory responses through NFkappaB inhibition: implications for Crohn’s disease. Gut 47, 397–403 (2000).
Rocha-Ramirez, L. M. et al. Probiotic Lactobacillus strains stimulate the inflammatory response and activate human macrophages. J. Immunol. Res. 2017, 4607491 (2017).
Udayan, S. et al. Macrophage cytokine responses to commensal Gram-positive Lactobacillus salivarius strains are TLR2-independent and Myd88-dependent. Sci. Rep. 11, 5896 (2021).
Pinacchio, C. et al. High abundance of genus Prevotella is associated with dysregulation of IFN-I and T cell response in HIV-1-infected patients. AIDS 34, 1467–1473 (2020).
Rhoades, N. et al. Altered immunity and microbial dysbiosis in aged individuals with long-term controlled HIV infection. Front. Immunol. 10, 463 (2019).
Iljazovic, A. et al. Perturbation of the gut microbiome by Prevotella spp. enhances host susceptibility to mucosal inflammation. Mucosal Immunol. 14, 113–124 (2021)
Tett, A. et al. The Prevotella copri complex comprises four distinct clades underrepresented in westernized populations. Cell Host Microbe 26, 666–679 e667 (2019).
Gupta, V. K., Paul, S. & Dutta, C. Geography, ethnicity or subsistence-specific variations in human microbiome composition and diversity. Front. Microbiol. 8, 1162 (2017).
Amour, C. et al. Epidemiology and impact of campylobacter infection in children in 8 low-resource settings: results from the MAL-ED study. Clin. Infect. Dis. 63, 1171–1179 (2016).
Lindenbaum, J., Harmon, J. W. & Gerson, C. D. Subclinical malabsorption in developing countries. Am. J. Clin. Nutr. 25, 1056–1061 (1972).
Korczak, B. M. et al. Genetic relatedness within the genus Campylobacter inferred from rpoB sequences. Int. J. Syst. Evol. Microbiol. 56, 937–945 (2006).
Doan, T. et al. Gut microbiome alteration in MORDOR I: a community-randomized trial of mass azithromycin distribution. Nat. Med. 25, 1370–1376 (2019).
Tickell, K. D., Atlas, H. E. & Walson, J. L. Environmental enteric dysfunction: a review of potential mechanisms, consequences and management strategies. BMC Med. 17, 181 (2019).
Keenan, J. D. et al. Azithromycin to reduce childhood mortality in sub-Saharan Africa. N. Engl. J. Med. 378, 1583–1592 (2018).
Gehrig, J. L. et al. Effects of microbiota-directed foods in gnotobiotic animals and undernourished children. Science 365, eaau4732 (2019).
Galpin, L. et al. Effect of Lactobacillus GG on intestinal integrity in Malawian children at risk of tropical enteropathy. Am. J. Clin. Nutr. 82, 1040–1045 (2005).
Barr, T. et al. Chronic ethanol consumption alters lamina propria leukocyte response to stimulation in a region-dependent manner. FASEB J. 33, 7767–7777 (2019).
Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140 (2010).
Tripathi, S. et al. Meta- and orthogonal integration of Influenza “OMICs” data defines a role for UBR4 in virus budding. Cell Host Microbe 18, 723–735 (2015).
Frishberg, A., Brodt, A., Steuerman, Y. & Gat-Viks, I. ImmQuant: a user-friendly tool for inferring immune cell-type composition from gene-expression data. Bioinformatics 32, 3842–3843 (2016).
Parada, A. E., Needham, D. M. & Fuhrman, J. A. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ. Microbiol. 18, 1403–1414 (2016).
Bolyen, E. et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol 37, 852–857 (2019).
Callahan, B. J. et al. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583 (2016).
Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780 (2013).
Price, M. N., Dehal, P. S. & Arkin, A. P. FastTree 2-approximately maximum-likelihood trees for large alignments. PLoS ONE 5, e9490 (2010).
Bokulich, N. A. et al. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome 6, 90 (2018).
Quast, C. et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596 (2013).
Lozupone, C., Lladser, M. E., Knights, D., Stombaugh, J. & Knight, R. UniFrac: an effective distance metric for microbial community comparison. ISME J. 5, 169–172 (2011).
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).
Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).
Abubucker, S. et al. Metabolic reconstruction for metagenomic data and its application to the human microbiome. PLoS Comput. Biol. 8, e1002358 (2012).
Suzek, B. E., Wang, Y., Huang, H., McGarvey, P. B. & Wu, C. H. UniRef clusters: a comprehensive and scalable alternative for improving sequence similarity searches. Bioinformatics 31, 926–932 (2015).
Karp, P. D., Riley, M., Paley, S. M. & Pellegrini-Toole, A. The MetaCyc database. Nucleic Acids Res. 30, 59–61 (2002).
Segata, N. et al. Metagenomic biomarker discovery and explanation. Genome Biol. 12, R60 (2011).
Segata, N. et al. Metagenomic microbial community profiling using unique clade-specific marker genes. Nat. Methods 9, 811–814 (2012).
Bankevich, A. et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput Biol. 19, 455–477 (2012).
Kang, D. D., Froula, J., Egan, R. & Wang, Z. MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities. PeerJ 3, e1165 (2015).
Wattam, A. R. et al. Assembly, annotation, and comparative genomics in PATRIC, the all bacterial bioinformatics resource center. Methods Mol. Biol. 1704, 79–101 (2018).
Oksanen, J. et al. R package “vegan”: Community Ecology Package. (Version 2.5-7), https://CRAN.R-project.org/package=vegan (2020).
Harrell, F. et al. R package “Hmisc”: Harrell Miscellaneous. (Version 4.5.0), https://github.com/harrelfe/Hmisc/ (2021).
Wei, T., Simko, V. R package “corrplot”: Visualization of a Correlation Matrix. (Version 0.89), https://github.com/taiyun/corrplot (2021).
Acknowledgements
This work was supported by the Bill and Melinda Gates Foundation (OPP1149233 and OPP1159806). N.R. was supported by NIH/NIAID T32 training grants (AI141346 and AI007319).
Author information
Authors and Affiliations
Contributions
I.M. and M.S. designed the study. S.H., K.P., A.H., and L.G. carried out the sample and clinical data collection and veterinary care. N.R. generated data and conducted computational/statistical analysis. N.R., M.S., and I.M. wrote the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Rhoades, N.S., Hendrickson, S.M., Prongay, K. et al. Growth faltering regardless of chronic diarrhea is associated with mucosal immune dysfunction and microbial dysbiosis in the gut lumen. Mucosal Immunol 14, 1113–1126 (2021). https://doi.org/10.1038/s41385-021-00418-2
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41385-021-00418-2
This article is cited by
-
Infant diarrheal disease in rhesus macaques impedes microbiome maturation and is linked to uncultured Campylobacter species
Communications Biology (2024)
-
Environmental enteric dysfunction: gut and microbiota adaptation in pregnancy and infancy
Nature Reviews Gastroenterology & Hepatology (2023)
-
Gut microbiota bridges the iron homeostasis and host health
Science China Life Sciences (2023)
-
Water, sanitation, hygiene practices, health and nutritional status among children before and during the COVID-19 pandemic: longitudinal evidence from remote areas of Dailekh and Achham districts in Nepal
BMC Public Health (2022)
-
Mucin O-glycan-microbiota axis orchestrates gut homeostasis in a diarrheal pig model
Microbiome (2022)