Skip to main content

Correlation between gut microbiota composition, enteric infections and linear growth impairment: a case–control study in childhood stunting in Pidie, Aceh, Indonesia



Gut microbiota is pivotal in maintaining children's health and well-being. The ingestion of enteric pathogens and dysbiosis lead to Environmental Enteric Dysfunction (EED), which is essential in stunting pathogenesis. The roles of gut microbiome and enteric infections have not been explored comprehensively in relation to childhood stunting in Indonesia. This study aimed to determine the correlation between gut microbiota composition, enteric infections, and growth biomarker, Insulin-like Growth Factor 1 (IGF-1), in stunted children from Pidie, Aceh, Indonesia.


This study was a case–control study involving 42 subjects aged 24 to 59 months, comprising 21 stunted children for the case and 21 normal children for the control group. The IGF-1 serum level was quantified using ELISA. The gut microbiome profiling was conducted using 16S rDNA amplicon sequencing. The expression of enteric pathogens virulence genes was determined using quantitative PCR (qPCR) assay. The correlations of observed variables were analysed using suitable statistical analyses.


The result showed that the IGF-1 sera levels in stunted were lower than those in normal children (p ≤ 0.001). The abundance of Firmicutes (50%) was higher than Bacteroidetes (34%) in stunted children. The gut microbiome profile of stunted children showed enriched genera such as Blautia, Dorea, Collinsella, Streptococcus, Clostridium sensu stricto 13, Asteroleplasma and Anaerostipes. Meanwhile the depleted genera comprised Prevotella, Lactococcus, Butyrivibrio, Muribaculaceae, Alloprevotella, Akkermansia, Enterococcus, Terrisporobacter and Turicibacter. The abundance of water biological contaminants such as Aeromonas, Stappiaceae, and Synechococcus was also higher in stunted children compared to normal children. The virulence genes expression of Enteroaggregative Escherichia coli (aaiC), Enterotoxigenic E. coli (estA), Enteropathogenic E. coli (eaeA), Shigella/Enteroinvasive E. coli (ipaH3) and Salmonella enterica (ompC) in stunted was higher than in normal children (p ≤ 0.001), which negatively correlated to height and level of IGF-1.


The present study showed the distinctive gut microbiome profile of stunted and normal children from Pidie, Aceh, Indonesia. The gut microbiota of stunted children revealed dysbiosis, comprised several pro-inflammatory, metabolic abnormalities and high-fat/low-fiber diet-related taxa, and expressed virulence genes of enteric pathogens. These findings provide evidence that it is imperative to restore dysbiosis and preserve the balance of gut microbiota to support linear growth in children.


Stunting or linear growth faltering is the most prevalent form of undernutrition. To date, 149.2 million children under five years are affected by stunting worldwide, and more than half of those children live in Asia. The Covid-19 pandemic may substantially increase the number due to limitations in accessing adequate nutrition and healthcare services [1]. Indonesia is one of the Asian countries still struggling to eradicate stunting. In 2022, Indonesia was fourth on the list of the countries with the highest prevalence of stunting, reaching 21.6% [2].

Stunting is more than just growth impairment. This undernourished state also deteriorates intellectual and cognitive ability. The quality of health will also decline due to future metabolic complications such as diabetes and other cardiovascular diseases. These will affect the quality of human resources and become significant threats to the nation's future. The devastating impacts of stunting can last for a lifetime and affect the next generation [1].

Previous studies have revealed the role of gut microbiota in the pathogenesis of stunting. The gut microbiota imbalance, namely dysbiosis, was the underlying mechanism caused by the chronic ingestion of enteric pathogens. The existence of enteric pathogens activated the gut immune system and led to local and systemic inflammation. These conditions disturb gut physiology by reducing the function of the gut barrier, increasing permeability, and impairing the gut structure. These states were defined as Environmental Enteric Dysfunction (EED). EED causes the impairment of nutrient absorption, resulting in nutrient deficiency [3].

A worldwide study revealed that childhood diarrhoea was one of the major causes of stunting in developing countries [4]. World Health Organization has proposed enteric infections and diarrhoeal diseases as contributing factors to stunting. Most studies in Indonesia revealed a strong correlation between stunting in children under 5 years of age with diarrhoeal episodes [5]. Previous studies revealed the common enteric pathogens associated with diarrhoea and linear growth impairment in children, such as Enteroaggregative Escherichia coli /EAEC [6], Enterotoxigenic E. coli/ETEC [7, 8], Enteropathogenic E. coli/EPEC [9], Enteroinvasive E. coli/EIEC [10], Shigella [6], and Salmonella [11].

Insulin-like Growth Factor 1 (IGF-1) is the most common growth biomarker used in studying linear growth impairment. The previous studies revealed the correlation between the low concentration of IGF-1, the increase of inflammation markers, and the presence of enteric infections in stunted children [12,13,14,15,16]. The secretion of IGF-1 was also a result of gut microbiota modulation activity, which was confirmed by previous animal studies [17,18,19]. The role of IGF-1 in linear growth in humans was indicated by the activity of IGF-1in supporting bone growth through the Growth Hormone (GH)-IGF-1 axis [20,21,22]. The low level of IGF-1 was also associated with several non-communicable diseases, such as diabetes [23]. Therefore, IGF-1 could be developed as a potential predictor or marker for long-term health status in children [24].

Stunting is widely distributed in Indonesia and Aceh is one of the provinces with a high prevalence of stunting, which is 33.2% in 2021 [2]. Not only does it have to deal with stunting, but Aceh also struggles with other health-related issues. A high prevalence of diarrhoea among toddlers was also reported in Aceh, reaching 13.8% in 2019. This number was higher than the national prevalence, which was 11%. Aceh was also the province with the lowest basic immunization coverage in Indonesia, reaching 50.9% [25]. The location of this study, Pidie district, was one of the areas in Aceh with high stunting cases [2]. It has been designated as one of the stunting loci in Indonesia [26].

The vast implication of the stunting impacts on health and development of this nation demands some progressive measures. Several established stunting eradication programs have successfully lowered the prevalence of stunting. The nutritional approaches have been mainly focused on the supplementation of macronutrients such as carbohydrates and protein and micronutrients such as zinc, vitamin A, and iron [27]. The success of supplementation relies significantly on the gut absorption function, which declines due to dysbiosis and EED in the case of stunting.

This study aimed to determine the correlation between gut microbiota composition, enteric infections, and growth biomarker, Insulin-like Growth Factor 1 (IGF-1), in stunted children from Pidie, Aceh, Indonesia. To our knowledge, the microbiome study of stunted children in Aceh has never been published. The results of this study are essential in capturing certain types of gut microbiota that play a pivotal role in linear growth. It provides beneficial information since gut microbiota exhibits dynamic features related to gender, age, diet, geographical regions, socio-demographic conditions, and health status. These findings may support the more comprehensive approaches to stunting prevention and eradication.


Characteristics of subjects

A total of 42 subjects consisting of 21 stunted and 21 normal children aged 24 to 59 months were enrolled in this study (Table 1). All subjects in the case group were categorized as stunted based on WHO growth standards [28]. The birthweight and BMI between observed groups was not statistically different. Most stunted children experienced bloody or watery diarrhoea (86%) and had incomplete basic immunization status (71%). The IGF-1 sera level in stunted children was lower than in normal children (Table 1).

Table 1 Characteristics and IGF-1 concentration of observed groups

Gut microbiome composition in stunted and normal children

The human gut microbiome comprises a complex microbial community that changes temporally and spatially, and in relation to health status [29]. The 16S rDNA amplicon sequencing is the most extensively used approach in determining gut microbiota composition [30]. A total of 14,423 Amplicon Sequence Variant (ASV) were observed in 16S rDNA amplicon sequencing and proceeded to bioinformatics analyses. Phylum-level abundance analysis showed a predominant portion of Firmicutes compared to the Bacteroidetes in stunted children (Fig. 1a), resulting in an increasing Firmicutes to Bacteroidetes (F/B) ratio. Differential abundance analysis of gut microbiome (Linear Discriminant Analysis (LDA) score > 3) in stunted children revealed the significant abundance of Streptococcus, Blautia, Dorea, Collinsella, Clostridium sensu stricto 13, Aeromonas, Stappiaceae, and Synechococcus CC9902 (Fig. 1b). Meanwhile, the important taxa in normal children comprised of Prevotella, Butyrivibrio, Muribaculaceae, and Alloprevotella (Fig. 1b). Weighted UniFrac (Fig. 1c) and Unweighted UniFrac (Additional file 1: Fig. S1) models of beta diversity showed the different structures of the gut microbiome community among observed groups (p = 0.042 and p = 0.001). Chao1 and Shannon index represented the gut microbiome diversity among samples/individuals. The rarefaction curve showed sufficient sequencing depth in this study (Additional file 1: Fig. S2). Chao1 value indicated the species richness, which was lower in stunted children than in normal children. Shannon index showed more complex species diversity, and it was not statistically different between observed groups (Fig. 1d).

Fig. 1
figure 1

Gut microbiota composition in normal and stunted children. a Relative abundance of bacterial phyla. b Significant different genera in stunted and normal children (LDA score > 3) c Weighted UniFrac model of beta diversity. d Chao1 and Shannon index of alpha diversity. The statistically significant difference was indicated by an asterisk (ns not significant; ***p ≤ 0.001)

Spearman analysis revealed the depleted and enriched taxa in stunted children (Table 2). Several depleted taxa were Prevotella, Lactococcus, Butyrivibrio, Muribaculaceae, Alloprevotella, Akkermansia, Enterococcus, Terrisporobacter and Turicibacter. In the other hand, Blautia, Dorea, Collinsella, Asteroleplasma, Streptococcus, Lactobacillus, Anaerostipes, Aeromonas, Stappiaceae, and Synechococcus CC9902 were found to be enriched in stunted children.

Table 2 Genus-level differences between stunted and normal children

Difference expression of virulence genes among observed groups

The result of 16S rDNA amplicon sequencing revealed the presence of several taxa related to enteric infections, such as Escherichia-Shigella and Salmonella (Additional file 2). Pathogenic bacteria express virulence factors to invade the host cells and retain intracellular survival, persistence, and spreading [31]. The present study used qPCR assay to determine the expression level of virulence genes from several enteric pathogens by measuring the relative quantification. The measurement revealed the higher expression of virulence genes of enteric pathogens in stunted children compared to normal children as follows: aaiC expression of EAEC (p < 0.001); eaeA expression of EPEC (p < 0.001); estA expression of ETEC (p < 0.001); ipaH3 expression of Shigella/EIEC (p < 0.001); and ompC expression of Salmonella enterica (p < 0.001)(Fig. 2).

Fig. 2
figure 2

Virulence genes expression in stunted and normal children. Statistical analysis was performed using Mann–Whitney U Test (p < 0.05). The statistically significant difference was indicated by an asterisk (***p ≤ 0.001)

Correlation between virulence gene expression and linear growth

The ingestion of enteric pathogens lead to dysbiosis and linear growth impairment [32]. The correlation matrix was used to visualize the correlation between the virulence genes expression, IGF-1, and height (Fig. 3). IGF-1 level was significantly associated with height (p = 0.007). There was a negative correlation between virulence genes expression and IGF-1 sera level, although it was not statistically significant. The expression of estA, gene encoding ETEC heat-stable enterotoxin, was negatively correlated to height (p = 0.020). There was a positive correlation among each virulence gene of enteric pathogen observed in this study as follows: the aaiC expression was positively correlated to ompC and eaeA expression (p < 0.001); the ompC expression was positively correlated to eaeA (p = 0.001), ipaH3 (p = 0.008) and estA expression (p = 0.006); the eaeA expression was positively correlated to ipaH3 (p = 0.03), and estA expression (p = 0.03); and the ipaH3 expression was positively correlated to estA (p = 0.001) (Fig. 3).

Fig. 3
figure 3

Correlation between the expression of virulence genes and level of IGF-1. The blue colour indicated a positive correlation, and the red showed a negative correlation. The analysis was performed using Spearman-rank test and statistically significant correlation was marked by an asterisk (*p ≤ 0.05; **p ≤ 0.01; ***p ≤ 0.001)


In accordance with the previous reports, this study showed a significant correlation between the history of diarrhoea and stunting [33,34,35]. Diarrhoea could be one of the clinical manifestations of EED, representing enteric pathogens exposure, although the exposure could be present with or without diarrhoea. Diarrhoea causes electrolyte imbalance and water loss, leading possibly to nutrition depletion [3]. High energy expenditure, malabsorption, and lower appetite due to diarrhoea contribute to malnutrition [36]. A dissimilar result from another study in Indonesia showed no association between a history of diarrhoea and stunting since it was only an acute episode; therefore, there was no significant effect on linear growth [37]. Subclinical infections have raised a concern about their contribution to linear growth retardation. The absence of recognized episodes of diarrhoea (three or more formless stools per day) and asymptomatic carriages of known enteric pathogens might lead to subclinical gut inflammation and impaired linear growth [38]. A study showed that atypical EPEC could retain gut colonization due to the absence of the bfpA gene. The infections commonly emerge as mild and subclinical, although some cases report more severe clinical symptoms [39].

Basic immunization in early life is fundamental for children to reach optimal growth. A secondary analysis study from 13 provinces in Indonesia demonstrated a positive correlation between incomplete immunization status and stunting in children aged 2 to 5 years [40]. A study in Thai children also revealed a similar outcome [41]. An extensive study in Africa showed that early-life vaccination decreased the odds of stunting and anemia. The unvaccinated children might be prone to recurring infections leading to prolonged inflammation [42]. A longitudinal study in Zimbabwe revealed that the prolonged inflammation in stunted infants was negatively correlated to IGF-1 levels and erythropoietin [12].

This study also revealed lower levels of IGF-1 sera in stunted children than in normal children. A study of stunted children aged 4 to 10 in Egypt also revealed a declined level of IGF-1 compared to the controls [15]. Another study in Indonesia showed a positive association between low IGF-1 levels and stunting in children with transfusion-dependent thalassemia [16]. The pro-inflammatory markers in stunting reduced the circulating IGF-1, impairing linear growth [12,13,14]. Children with higher levels of IGF-1 during the first year of life were less likely to develop stunting [14]. However, a different result showed that IGF-1 levels were independently associated with stunting in Bangladeshi children aged 12 to 18 months [43].

The present study used IGF-1 as a growth biomarker in stunting since it was more sensitive to nutritional deficiency than another biomarker, namely IGF Binding Protein 3 (IGFBP3) [44]. The IGF-1 was more stable than GH and also unaffected by diurnal fluctuations and pre-analytical factors such as food intake, exercise, and stress proximately before blood sampling [45]. Free form of IGF-1 in serum, as was used in this study, demonstrated more prominent physiological activity and clinical relevance than total IGF-1 (bound and free form) [46]. A concern has been raised due to the low concentration of free IGF-1 serum in several catabolic conditions such as malnutrition, anorexia nervosa, and poorly controlled type 1 diabetes, as well as in hypothyroidism [45]. The concern was no longer applied since these conditions were not found in all subjects recruited. Further research is needed to elaborate on the relationship between IGF-1 levels and stunting.

The increased F/B ratio in this study was also found in other reports, such as in West Java [47] and Chimalhuacán, Mexico. The higher abundance of Firmicutes than Bacteroidetes indicated dysbiosis and inflammation, and it was strongly correlated to high sugar-low fiber intake and lipid metabolism impairment, resulting in childhood obesity [48].

The high abundance of several taxa in stunted children also indicated dysbiosis, inflammation, and an imbalanced diet. Blautia was associated with high levels of long-chain triglycerides and a pre-diabetic state in patients with Type I Diabetes Mellitus. Blautia and Dorea represented a high fat intake and significantly correlated to inflammation and the impairment of the gut barrier [49]. Blautia and Clostridium sensu stricto 13 were positively correlated to serum lipids [50]. Streptococcus was a pro-inflammatory microbiota and commonly found in systemic inflammation states. Collinsella was associated with low-fiber intake, hyperglycemia state, and decreased intestinal integrity [51]. The presence of metabolic abnormalities-related taxa may become a predictor of future metabolic complications in stunted children.

Interestingly, this study found a higher abundance of water biological contaminants in stunted children than in normal children, such as Aeromonas, Stappiaceae, and Synechococcus CC9902. Aeromonas is commonly known as fish pathogens and naturally reside in aquatic environments. This organism has gained a warning on human health due to its ability to colonize and infect humans. Acute gastroenteritis was a common clinical manifestation of Aeromonas infection [52]. A study revealed Aeromonas as the significant pathogens in moderate to severe diarrhoea in children from Bangladesh and Pakistan, and it was also associated to the degree of stunting in children aged 24 to 59 months. Aeromonas-associated diarrhoea was co-occurred with dysentery caused by Shigella spp. [53]. Stappiaceae is an alphaproteobacterium that resides in marine environments, especially in mangrove sediments. It is commonly found in shrimps and oysters [54]. Synechococcus, one of the cyanobacteria, has been known as a seafood contaminant. The contamination occurred due to the accumulation of cyanotoxins on the seaweed, crustaceans, or fish. A high accumulation of toxins could harm humans [55]. The cyanotoxin-producing Synechococcus spp. has been reported on freshwater environments in Singapore. This organism produced cylindrospermopsin (CYN) and anatoxin-a [56]. The CYN of Cyanobacteria showed effects on liver, small intestine and colon. The lipopolysaccharides of Cyanobacteria affected the intestinal and immune cells, although the impacts on humans were still debatable [57]. Since the study site was located near the coastal area and most subjects consumed marine products as a primary source of proteins, the contamination may be caused by the ingestion of contaminated products. Further investigation regarding the sources of contamination and daily food processing and handling is needed to elaborate more on the relationship between those biological contaminants and stunting. Whether the contamination is site-specific is also interesting to be explored comprehensively.

Prevotella, Butyrivibrio, and Muribaculaceae were predominantly found in normal children. This microbiota was associated with low fat, high-fiber diet, carbohydrate intake, and increased propionate production. The Alloprevotella genus has been known for its anti-inflammatory effect on the gut and ability to produce butyrate [58].

Beta diversity metrics, Weighted and Unweighted UniFrac, demonstrated that nutritional status significantly influenced the microbiota diversity between observed groups. The present study showed a lower species richness (Chao1 alpha diversity metric) in stunted children than in normal children. This finding was also revealed in another study in India [59]. However, a dissimilar result from previous studies showed higher alpha diversity in stunted children [47] and no difference between stunted and normal children [60, 61]. Since it represents strong functional and ecological stability, higher alpha diversity is mainly found in a healthy or disease-improvement state [29]. Low-diversity gut microbiota environment was more susceptible to enteric infections [62]. However, different results in alpha diversity were relevant to the dynamic features of gut microbiota and their contributing factors.

The present study revealed a number of depleted and enriched genera in stunted children. However, dissimilar results showed the depleted abundance of Blautia and Collinsella and enriched abundance of Akkermansia [61], Butyrivibrio and Alloprevotella [47] in stunted children. Both studies revealed a higher abundance of Copprococcus. Meanwhile, the present study did not show a significant difference in the abundance of Coprococcus between observed groups. The exploration of the significant taxa based on their function revealed the beneficial and harmful effects of depleted and enriched genera in stunted children. Enterococcus, a commensal microbiota that was depleted in stunted children, played an important role in maintaining colonic homeostasis and decreasing the severity of infectious diarrhoea in children. It has been developed as a probiotic, although there were also some reports regarding hospital-associated infections [63]. Terrisporobacter and Turicibacter were enriched genera in healthy persons with normal glycemic levels [64]. Asteroleplasma, an enriched genus in stunted children, was also enriched in Type 2 Diabetes Mellitus and positively correlated to the increased HbA1c level [65]. Lactobacillus, a well-known probiotic genus, has been associated with intestinal inflammation, such as Crohn’s disease and Rheumatoid Arthritis. Whether Lactobacillus was directly involved in pathogenesis or adapted to survive the pro-inflammatory gut environment was poorly understood [66]. The present study showed an enriched abundance of Anaerostipes in stunted children. A similar result was also revealed in Pakistani stunted children, and Anaerostipes was negatively correlated to height [66]. However, other studies revealed the beneficial roles of Anaerostipes, especially regarding its ability to produce butyrate from the degradation of carbohydrates and commonly depleted in multiple diseases [67, 68]. Gut microbiota composition is dynamic and fluctuating in response to external factors and internal biological processes. Defining the ‘abnormal’ microbiome is challenging due to the temporal dynamics of gut microbiota [69], as well as distinct variations of gender, age, geographical, social-demographic, and medical contexts. Therefore, it was indispensable to investigate the ‘healthy’ microbiome in a proper setting [70].

This study discovered a higher expression of enteric pathogen virulence genes in stunted children than in normal children. To our knowledge, the virulence genes expression study in stunted children has never been reported. In the present study, the virulence genes not only existed but were also assuredly expressed, which later were translated into proteins and played critical roles in facilitating enteric infections. Gut inflammation contributed to the loss of microbial density and diversity, therefore supporting the growth of the Enterobacteriaceae family due to the lack of competition with other endogenous microbiota. The accumulation of Sulphur-containing by-product tetrathionate and nitrate from gut inflammation supported the rapid outgrowth of Salmonella typhimurium and E. coli [71]. The depletion of Enterococcus and Alloprevotella, which played essential roles in maintaining gut homeostasis and exhibiting anti-inflammatory effects, might promote susceptibility to enteric infections. However, this should be further explored.

Previous studies showed that the presence of EAEC harbouring aaiC gene was negatively correlated to linear growth in Bangladesh and Peru [72, 73]. It was explained that the correlations were site-specific [73]. EAEC and ETEC were the most prevalent causes of diarrhoea in children less than five years of age in China. Hot and warm temperatures contributed to the high prevalence [74].

The present study demonstrated a higher expression of ETEC’s heat-stabile toxin-encoded gene (estA) in stunted children than in normal children. ETEC was one of four diarrhoeal pathogens associated with moderate to severe diarrhoea cases among children under five in Asia and Africa. ETEC-associated stunting was reported in low-income and lower-middle-income countries, and it increased mortality [75]. A previous study showed that asymptomatic heat-stable toxin ETEC infections were significantly associated with childhood stunting [76].

This study only used eaeA to confirm the existence of EPEC among samples; therefore, it could not be differentiated whether the infections were typical or atypical. EPEC has been known as the most prevalent E. coli pathotype in children under five with diarrhoea [76]. A study in Bangladesh indicated that EPEC was associated with stunting and underweight [10].

Shigella/EIEC was the predominant cause of linear growth impairment in children with diarrhoea [77]. Shigella infections showed the strongest negative association with linear growth in children from developing countries [6] and increased mortality rate in stunted children [75].

This study used the ompC gene to reveal the existence of Salmonella enterica (S. enterica) among fecal specimens; however, the exact serovar was undetermined. A study showed that Salmonella was detected in fecal samples of stunted children in India using the culture method, although it was not statistically significant [78]. Another study revealed that non-typhoidal Salmonella infections were significantly associated with wasting but not stunting [11].

The present study highlighted the negative correlation between the expression of virulence genes, IGF-1 level, and height. IGF-1 through GH/IGF-1 axis is essential in stimulating the proliferation and differentiation of chondrocytes in the epiphyseal plate. IGF-1 also induces osteoblast proliferation, collagen secretion, and bone matrix mineralization, thus supporting linear growth [20,21,22]. The modulation by gut microbiota might take part in IGF-1 secretion, as previously reported in animal studies [17,18,19]. A previous study reported the correlation between diarrhoea due to enteric infections, IGF-1 concentration, systemic inflammation biomarker, and linear growth impairment in stunted children aged 6 to 24 months from Brazil [13].

Exposure to enteric pathogens caused EED. The changes in intestinal structures and functions and systemic inflammation have led to impairment of nutrient adsorption, growth hormone resistance, and disruption of bone growth and remodelling, thus resulting in linear growth faltering [32]. Diarrhoea magnified these conditions, especially by promoting nutrition depletion. Previous studies revealed the association between enteric infections and inflammation. ETEC, Shigella/EIEC, and EAEC were associated with increased EED scores in stunted children [73]. Previous studies revealed the association between enteric infections and the increase of inflammation markers, as follows: EAEC, EPEC, and ETEC were positively associated with the level of fecal calprotectin [10, 75]; Shigella, EPEC, and ETEC were positively associated with the levels of myeloperoxidase, and EIEC was positively associated with the levels of lactoferrin, interleukin (IL) 8, and IL-1b [79]. Acute Salmonella infection increased intestinal permeability. Chronic infection elevated the level of myeloperoxidase and activated the pro-inflammatory pathways, thus increasing the susceptibility to intestinal inflammation [80]. However, this study was unable to determine whether subjects exhibited the EED since the EED parameters were not evaluated.

Coinfection of enteric pathogens has been reported in pediatric diarrhoeal cases [81], and it might increase virulence. The existence of certain pathogens might support or facilitate others. A previous study showed that EAEC was coinfected with ETEC and EPEC [10]. The concurrence infections of enteric pathogens in children under five with diarrhoea predominantly happened in under 2-year-old children compared to 2–5 years of age [82]. The present study showed the existence of multiple enteric pathogen virulence genes. Whether the infections occurred in the same or different periods was undeterminable in this study.

The present study revealed a significant association between a history of diarrhoea and stunting. Furthermore, the expression of enteric pathogen virulence genes was higher in stunted children than in normal children. However, the lack of detailed information regarding types, onset, duration, diagnosis, and treatments of diarrhoea, as well as other clinical symptoms, was one of the limitations of this study. Recall bias from parents/guardians, as well as undocumented event or incomplete information regarding history of diarrhoea in Maternal and Child Report, were the major causes of the shortcomings. There was also no comprehensive survey regarding dietary intake (nutritional survey), which restricted the exploration of the association between dietary pattern and gut microbiota composition in stunted and normal children. An extensive study involving more subjects will provide more conclusive results regarding the gut microbiome profile in stunted and normal children and their role in linear growth.


The present study showed that a history of diarrhoea and incomplete basic immunization status were strongly associated with stunting. It also discovered the distinctive gut microbiome composition in stunted and normal children from Pidie, Aceh, Indonesia. Gut microbiome profile revealed dysbiosis, a higher abundance of inflammation and metabolic abnormalities-related taxa, and unhealthy diets (high-fat/low-fiber diet). This study highlighted a higher abundance of water biological contaminants in gut microbiome of stunted children. The virulence genes expression study indicated a more vulnerable gut environment to enteric pathogens in stunted children, which might lead to linear growth retardation. These findings corroborate the importance of restoring dysbiosis and preserving the balanced state of gut microbiota in supporting linear growth in children. The existing stunting eradication programs need to be emphasized on several critical points. Prevention and adequate treatments of childhood diarrhoea and improvement of immunization coverage will provide better protection. Proper food and drinking water processing and adequate hygiene and sanitation are essential in reducing exposure to enteric pathogens. It must be persuasively and consistently educated to mothers and other food providers. Consuming healthy fermented food, adequate dietary fibers, and a proportional intake of fats, carbohydrates, and proteins may be beneficial in altering and shaping the balance of gut microbiota composition. Future research is encouraged to elaborate more on the relationship between water biological contaminants and linear growth impairment and explore the interactions/co-occurrences between enteric pathogens and other microbiota in stunted and normal children.


Study design and subject recruitment

This study was a case–control study involving 42 participants recruited using a purposive sampling method in Pidie District, one of the stunting pocket areas in Aceh Province, Indonesia (5°18′16.5"N 96°00′35.3" E). The sample size estimation was determined using the Lameshow formula for case–control study [83]. It was calculated based on the proportion of stunting in children under 5 years in Aceh Province, Indonesia (21.1%), as reported in Indonesia Basic Health Research in 2019 [84]. The protocol of this study was approved by the Research Ethics Committee of the Faculty of Medicine, Universitas Syiah Kuala (dossier No.050/EA/FK-RSUDZA/2021). Written informed consent was acquired from the parents or guardians of the participants in the presence of a third person. Study participants for the case (n = 21) and control groups (n = 21) were recruited based on the following criteria: weaned toddlers, 24–59 months of age, and no history of antibiotic treatment within the previous four weeks. Participants were of any gender or ethnicity and had no acute or chronic illnesses, symptoms, or congenital diseases. The length/height for age and weight for age for the case group was between -3 standard deviation (SD) and -2 SD of the median value of WHO Child Growth Standards. On the other hand, the value of length/height for age and weight for age for the control group was above -2 SD and below + 2 SD of the median value of standards, which meant that the participants were within the normal range and other types of malnutrition were excluded. The selected normal children for the control group were of the same age, sex, and living areas as the stunted children, with a 1:1 comparison. The anthropometric measurements of all participants were conducted in Puskesmas (community health center) facility from January to May 2021. The measurement was performed by trained health workers using standardized methods and equipment and extrapolated to stunting versus normal based on the WHO child growth standard [28]. Kartu Menuju Sehat (Health Card) and Buku Kesehatan Ibu dan Anak (Maternal and Child Report) was used to confirm the health status of the subjects, including the record of growth, health and immunization status. History of diarrhoea was collected using direct interview with parents/guardians. It referred to previous recurrent incidents of acute diarrhoea with three or more bloody or watery stools daily and happened for < 7 days with or without hospitalization. Basic immunization status was indicated as the fulfilment of mandatory vaccinations (Expanded Programs of Immunization/EPI) of participants from 0 to 23 months consisted of Bacille Calmette-Guérin (BCG), Hepatitis B, diphtheria, tetanus, pertussis-Haemophilus influenzae type B-Hepatitis B (DTP-Hib-HepB), Polio, and Measles.

IGF-1 measurement using ELISA

The sera isolated from the whole blood were stored at − 20 °C and used for the immunochemistry analysis to measure the level of human IGF-1 (Human IGF-1 ELISA Kit, Sigma-Aldrich) using the sandwich ELISA method. The intensity of the colour was measured by spectrophotometer (Multiskan Go, Thermo Scientific) at 450 nm. The amounts of IGF-1 were derived by interpolation from the reference/standard curve generated from the known concentrations of IGF-1 in the same assay. The results of unknown samples were calculated using the four-parameter logistic function [85].

16S rDNA amplicon sequencing

About 1 g of the faecal specimen for each sample was mixed into 9 ml DNA-RNA shield solution (Zymo Research, USA) and stored at -80 °C until further used for DNA and RNA extraction. The DNA from fecal specimens was extracted using Quick-DNA Fecal/Soil Microbe Miniprep Kit (Zymo Research). The genomic DNA for metagenomics analysis was ensured to meet several requirements such as amount (≥ 200 ng), concentration (≥ 20 ng/µl), volume (≥ 12 ul), and purity (OD 260/280 = 1.8–2.0, with no degradation or contamination). The amount of DNA was quantified using Nanodrop Spectrophotometer (Maestrogen). The degradation of DNA was assessed using agarose gel electrophoresis. The DNA extraction was performed by PT. Genetika Science Indonesia, a metagenomics services provider. The gut microbiota profiling was determined using NGS (Illumina NovaSeq 6000, Novogene), which targets the V3-V4 region of bacterial 16S rDNA. The V3-V4 region was amplified using forward primer 341F 5’-CCTAYGGGRBGCASCAG-3’ and reverse primer 806R 5’- GGACTACNNGGGTATCTAAT-3’ [86, 87]. The metagenomics analysis workflow consisted of amplification using Polymerase Chain Reaction (PCR), library preparation, sequencing, and Bioinformatics analysis. Quality control was performed in all steps to ensure the quality and validity of data. Bioinformatics analyses were conducted using the microbiome bioinformatics platform QIIME2 version 2022.11 ( in processing the reads, picking, and classifying the ASV (SILVA version 138) as well as measuring microbiome diversity. The prevalence cut-off was set to 10%, which meant that ASVs found in fewer than 10% of total samples were removed.

Virulence genes expression analysis

The RNA from fecal specimens was extracted using RNAeasy Power Microbiome Kit (Qiagen). The process was optimized to fulfil the requirement, such as purity (OD 260/280 was 1.8–2.0 and A260/230 was above 2) and sufficient concentration for expression study. The RNA was used as a template for cDNA synthesis (ReverTraAce with gDNA Remover, Toyobo). Each reaction with 10 µl of total volume consisted of 150 ng RNA, 2 ul 4 × DN Master Mix (contained DNAse I), and 2 ul 5 × RT Master Mix II. The cDNA synthesis program consisted of denaturation of RNA at 65 °C for 5 min, followed by genomic DNA removal at 37 °C for 5 min, cDNA synthesis at 37 °C for 15 min, and final heat at 98 °C for 5 min.

The expression study was conducted by quantitative real-time PCR (CFX96 Touch Real-Time PCR, BioRad) using DNA-binding fluorescent dye/SYBR Green (SensiFAST SYBR® No-ROX, Bioline). The amplifications were carried out in triplicate. Each reaction (20 µl of total volume) consisted of 10 µl of 2 × SensiFAST SYBR® No-ROX Mix (Bioline), forward and reverse primers, 1 µl cDNA, and nuclease-free water. The amplification program consisted of polymerase activation at 95 °C for 2 min, followed by 40 cycles of denaturation at 95 °C for 5 s, annealing for 10 s, and elongation at 72 °C for 20 s. The bacterial 16S rDNA gene (rrs) was used as a reference gene and it has shown a stabile performance in observed groups (Additional file 1: Table S1). Since the efficiency values of qPCR reaction using the set of primers were not equal to 2, the relative quantification was calculated using the Pfaffl method with efficiency-corrected calculation [88]. The list of primers used in this study was shown in Additional file 1: Table S2.

Statistical analyses

The difference between groups in anthropometric values and alpha diversity Shannon index were analysed using an independent T-test. The categorical data were analysed using Fisher Exact Test. The level of IGF-1 and Alpha diversity Chao1 were analysed using Mann–Whitney U Test due to the non-normal distribution of data. Beta diversity analysis was performed using PERMANOVA with Weighted and Unweighted UniFrac models. The differential abundance analysis of statistically significant ASV between groups (normal vs. stunted) was performed using Linear Discriminant Analysis Effect Size (LEfSe). The association of taxonomic abundance and nutritional status was calculated using the Mann–Whitney U Test and corrected using False Discovery Rate (FDR) by Benjamin-Hochberg with a significant threshold of 10%. The correlation between virulence gene expression and nutritional status was evaluated using Mann–Whitney U Test. The correlation matrix between virulence genes expression, IGF-1 sera level, and height was constructed using the Spearman-rank test. All statistical analyses were performed using the R software package version 4.3.0 [89]. The results of statistical analyses were visualized using R studio version 2023.03.0 [90] and GraphPad Prism version 8.4.3.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article and its additional files. Further request for information can be directed to the corresponding author.



Insulin-like Growth Factor 1


Enzyme Linked Immunosorbent Assay


Ribosomal Deoxyribonucleic acid


Quantitative Polymerase Chain Reaction


Enteroaggregative Escherichia coli


Enterotoxigenic Escherichia coli


Enteropathogenic Escherichia coli


Enteroinvasive Escherichia coli


Environmental Enteric Dysfunction


Body Mass Index


Bacille Calmette-Guérin


Diphtheria, tetanus, pertussis-Haemophilus influenzae type B-Hepatitis B


Linear discriminant analysis effect size


Linear Discriminant Analysis


Growth hormone




Standard seviation


World Health Organization


Amplicon Sequence Variant


Optical density


False discovery rate


  1. UNICEF/WHO/World Bank Group. Levels and trends in child malnutrition. World Health Organization. 2021. Accessed 15 Jan 2022.

  2. Ministry of Health. Indonesia Nutritional Survey 2022. Jakarta, Indonesia; 2023. Accessed 20 Jun 2023.

  3. Budge S, Parker AH, Hutchings PT, Garbutt C. Environmental enteric dysfunction and child stunting. Nutr Rev. 2019;77:240–53.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Danaei G, Andrews KG, Sudfeld CR, Fink G, McCoy DC, Peet E, et al. Risk factors for childhood stunting in 137 developing countries: a comparative risk assessment analysis at global, regional, and country levels. PLoS Med. 2016;13:1–18.

    Article  Google Scholar 

  5. Beal T, Tumilowicz A, Sutrisna A, Izwardy D, Neufeld LM. A review of child stunting determinants in Indonesia. Matern Child Nutr. 2018;14:1–10.

    Article  Google Scholar 

  6. Rogawski ET, Liu J, Platts-Mills JA, Kabir F, Lertsethtakarn P, Siguas M, 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 Heal. 2018;6:e1319–28.

    Article  Google Scholar 

  7. Liu J, Gratz J, Amour C, Nshama R, Walongo T, Maro A, et al. Optimization of quantitative PCR methods for enteropathogen detection. PLoS ONE. 2016;11:1–11.

    Google Scholar 

  8. Bagamian KH, Anderson JD, Muhib F, Cumming O, Laytner LA, Wierzba TF, et al. Heterogeneity in enterotoxigenic Escherichia coli and shigella infections in children under 5 years of age from 11 African countries: a subnational approach quantifying risk, mortality, morbidity, and stunting. Lancet Glob Heal. 2020;8:e101–12.

    Article  Google Scholar 

  9. Nasrin D, Blackwelder WC, Sommerfelt H, Wu Y, Farag TH, Panchalingam S, et al. Pathogens associated with linear growth faltering in children with diarrhea and impact of antibiotic treatment: the global enteric multicenter study. J Infect Dis. 2021;224:S848–55.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Gazi MA, Alam MA, Fahim SM, Wahid BZ, Khan SS, Islam MO, et al. Infection With Escherichia coli pathotypes is associated with biomarkers of gut enteropathy and nutritional status among malnourished children in Bangladesh. Front Cell Infect Microbiol. 2022;12:1–9.

    Article  Google Scholar 

  11. Das R, Haque MA, Chisti MJ, Faruque ASG, Ahmed T. Association between non-typhoidal salmonella infection and growth in children under 5 years of age: Analyzing data from the global enteric multicenter study. Nutrients. 2021;13:1–11.

    Article  CAS  Google Scholar 

  12. Prendergast AJ, Rukobo S, Chasekwa B, Mutasa K, Ntozini R, Mbuya MNN, et al. Stunting is characterized by chronic inflammation in zimbabwean infants. PLoS ONE. 2014;9:e86928.

    Article  PubMed  PubMed Central  Google Scholar 

  13. DeBoer MD, Scharf RJ, Leite AM, Férrer A, Havt A, Pinkerton R, et al. Systemic inflammation, growth factors, and linear growth in the setting of infection and malnutrition. Nutrition. 2017;33:248–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Syed S, Manji KP, McDonald CM, Kisenge R, Aboud S, Sudfeld C, et al. Biomarkers of systemic inflammation and growth in early infancy are associated with stunting in young Tanzanian children. Nutrients. 2018;10:1–14.

    Article  Google Scholar 

  15. Abdou S, El-Boghdady N, Abd El-Maksoud A, Khairy S, El-Sawalhi M. Evaluation of insulin-like growth factor-1, total ghrelin, and insulin resistance in nutritionally stunted Egyptian children. Bull Fac Pharmacy, Cairo Univ. 2019;57:55–65.

    Article  Google Scholar 

  16. Pratiwi IGAPE, Irawan R, Ugrasena IDG, Faizi M. Vitamin D, insulin-like growth factor-1, and stunting in children with transfusion-dependent thalassemia. Paediatr Indones Indones. 2022;62:98–103.

    Article  Google Scholar 

  17. Yan J, Herzog JW, Tsang K, Brennan CA, Bower MA, Garrett WS. Gut microbiota induce IGF-1 and promote bone formation and growth. PNAS. 2016.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Schwarger M, Makki K, Storelli G, Machuca-Gayet I, Strutkova D, Hermanova P, et al. Lactobacillus plantarum strain maintains growth of infant mice during chronic undernutrition. Microbiome. 2016;351:854–7.

    Google Scholar 

  19. Poinsot P, Schwarzer M, Peretti N, Leulier F. 40 years of IGF1: The emerging connections between IGF1, the intestinal microbiome, lactobacillus strains and bone growth. J Mol Endocrinol. 2018;61:T103–13.

    Article  CAS  PubMed  Google Scholar 

  20. Laron Z. IGF-1 and insulin as growth hormones. In: Bock, G and Goode, J, editors. Biology of IGF-1: Its Interaction with Insulin in Health and Malignant States. Novartis Foundation. 2008;262:56–77.

  21. Blum WF, Alherbish A, Alsagheir A, El Awwa A, Kaplan W, Koledova E, et al. The growth hormone-insulin-like growth factor-I axis in the diagnosis and treatment of growth disorders. Endocr Connect. 2018;7:R212–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Dixit M, Poudel SB, Yakar S. Effects of GH/IGF axis on bone and cartilage. Mol Cell Endocrinol. 2021;519: 111052.

    Article  CAS  PubMed  Google Scholar 

  23. Bourdon C, Lelijveld N, Thompson D, Dalvi PS, Gonzales GB, Wang D, et al. Metabolomics in plasma of Malawian children 7 years after surviving severe acute malnutrition: “ChroSAM” a cohort study. EBioMedicine. 2019. hal. 464–72.

  24. Pham TPT, Raoult D, Million M. IGF1 levels in children with severe acute malnutrition after nutritional recovery: a good predictor for children’s long-term health status. EBioMedicine. 2019;45:9–10.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Ministry of Health. Indonesia Health Profile 2019. 2020. Accessed 10 Sep 2021.

  26. Ministry of National Development Planning of the Republic of Indonesia. Decree of Ministry of National Development Planning of the Republic of Indonesia, Number 42 year of 2020 regarding the designation of Stunting Loci in Indonesia. 2020. Accessed 15 May 2023.

  27. Ministry of National Development Planning of the Republic of Indonesia. National Strategic Plan in Stunting Eradication. 2019. Accessed 12 Jan 2020.28.

  28. World Health Organization. WHO Child Growth Standards. 2006.

  29. Li Z, Zhou J, Liang H, Ye L, Lan L, Lu F, et al. Differences in alpha diversity of gut microbiota in neurological diseases. Front Neurosci. 2022;16:1–13.

    Google Scholar 

  30. Kameoka S, Motooka D, Watanabe S, Kubo R, Jung N, Midorikawa Y, et al. Benchmark of 16S rRNA gene amplicon sequencing using Japanese gut microbiome data from the V1–V2 and V3–V4 primer sets. BMC Genomics. 2021;22:1–10.

    Article  Google Scholar 

  31. Denzer L, Schroten H, Schwerk C. From gene to protein—how bacterial virulence factors manipulate host gene expression during infection. Int J Mol Sci. 2020;21:1–37.

    Article  Google Scholar 

  32. Mbuya MNN, Humphrey JH. Preventing environmental enteric dysfunction through improved water, sanitation and hygiene: an opportunity for stunting reduction in developing countries. Matern Child Nutr. 2016;12:106–20.

    Article  PubMed  Google Scholar 

  33. Ezeh OK, Abir T, Zainol NR, Mamun A Al, Milton AH, Haque MR, et al. Trends of stunting prevalence and its associated factors among Nigerian children aged 0–59 months residing in the northern Nigeria, 2008–2018. Nutrients. 2021;13:2008–18.

  34. Arini D, Nursalam N, Mahmudah M, Faradilah I. The incidence of stunting, the frequency/duration of diarrhea and Acute respiratory infection in toddlers. J Public health Res. 2020;9:117–20.

    Article  Google Scholar 

  35. Nasrin D, Liang Y, Powell H, Casanova IG, Sow SO, Hossain MJ, et al. Moderate-to-severe diarrhea and stunting among children younger than 5 years: findings from the vaccine impact on diarrhea in Africa (VIDA) Study. Clin Infect Dis. 2023;76:S41–8.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Wasihun AG, Dejene TA, Teferi M, Marugán J, Negash L, Yemane D, et al. Risk factors for diarrhoea and malnutrition among children under the age of 5 years in the Tigray Region of Northern Ethiopia. PLoS ONE. 2018;13:32–9.

    Article  Google Scholar 

  37. Setyani LI, Anwar K. The association between feeding practices and history of diarrhea with nutritional status of toddlers. Indones J Public Heal Nutr. 2022;3:43–9.

    Google Scholar 

  38. Rogawski ET, Guerrant RL. The burden of enteropathy and subclinical infections. Pediatr Clin North Am. 2017;64:1–17.

    Article  Google Scholar 

  39. Hu J, Torres AG. Enteropathogenic Escherichia coli: foe or innocent bystander? Clin Microbiol Infect. 2015;21:729–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Fajariyah RN, Hidajah AC. Correlation between immunization status and mother’s height, and stunting in children 2–5 years in Indonesia. J Berk Epidemiol. 2020;8:89.

    Article  Google Scholar 

  41. Shinsugi C, Mizumoto A. Associations of nutritional status with full immunization coverage and safe hygiene practices among Thai children aged 12–59 months. Nutrients. 2022;14:34.

    Article  CAS  Google Scholar 

  42. Berendsen MLT, Smits J, Netea MG, van der Ven A. Non-specific effects of vaccines and stunting: timing may be essential. EBioMedicine. 2016;8:341–8.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Hossain M, Nahar B, Haque A, Mondal D, Mahfuz M, Naila NN, et al. Serum adipokines, growth factors, and cytokines are independently associated with stunting in Bangladeshi children. Nutrients. 2019;11:1–16.

    Article  CAS  Google Scholar 

  44. Gupta N, Lustig RH, Kohn MA, McCracken M, Vittinghoff E. Sex differences in statural growth impairment in Crohn’s disease: role of IGF-1. Inflamm Bowel Dis. 2011;17:2318–25.

    Article  PubMed  Google Scholar 

  45. Schilbach K, Olsson DS, Boguszewski MCS, Bidlingmaier M, Johannsson G, Jørgensen JOL. Biomarkers of GH action in children and adults. Growth Horm IGF Res. 2018;40:1–8.

    Article  CAS  PubMed  Google Scholar 

  46. Gómez JM, Maravall FJ, Gómez N, Navarro MA, Casamitjana R, Soler J. Interactions between serum leptin, the insulin-like growth factor-I system, and sex, age, anthropometric and body composition variables in a healthy population randomly selected. Clin Endocrinol (Oxf). 2003;58:213–9.

    Article  PubMed  Google Scholar 

  47. Surono IS, Widiyanti D, Kusumo PD, Venema K. Gut microbiota profile of Indonesian stunted children and children with normal nutritional status. PLoS ONE. 2021;16:1–18.

    Google Scholar 

  48. Méndez-Salazar EO, Ortiz-López MG, Granados-Silvestre MDLÁ, Palacios-González B, Menjivar M. Altered gut microbiota and compositional changes in firmicutes and proteobacteria in Mexican undernourished and obese children. Front Microbiol. 2018;9:1–11.

    Google Scholar 

  49. Vacca M, Celano G, Calabrese FM, Portincasa P, Gobbetti M, De Angelis M. The controversial role of human gut lachnospiraceae. Microorganisms. 2020;8:1–25.

    Article  Google Scholar 

  50. Zeng Q, Li D, He Y, Li Y, Yang Z, Zhao X, et al. Discrepant gut microbiota markers for the classification of obesity-related metabolic abnormalities. Sci Rep. 2019;9:1–10.

    Article  Google Scholar 

  51. Li BY, Xu XY, Gan RY, Sun QC, Meng JM, Shang A, et al. Targeting gut microbiota for the prevention and management of diabetes mellitus by dietary natural products. Foods. 2019;8:440.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Pessoa RBG, de Oliveira WF, dos Correia MTS, Fontes A, Coelho LCBB. Aeromonas and human health disorders: clinical approaches. Front Microbiol. 2022;13:868890.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Qamar FN, Nisar MI, Quadri F, Shakoor S, Sow SO, Nasrin D, et al. Aeromonas-associated diarrhea in children under 5 years: the gems experience. Am J Trop Med Hyg. 2016;95:774–80.

    Article  PubMed  PubMed Central  Google Scholar 

  54. de Santana CO, Spealman P, Melo V, Gresham D, de Jesus T, Oliveira E, et al. Large-scale differences in diversity and functional adaptations of prokaryotic communities from conserved and anthropogenically impacted mangrove sediments in a tropical estuary. PeerJ. 2021;9:e12229.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Abeysiriwardena NM, Gascoigne SJL, Anandappa A. Algal bloom expansion increases cyanotoxin risk in food. Yale J Biol Med. 2018;91:129–42.

    CAS  PubMed  PubMed Central  Google Scholar 

  56. Gin KYH, Sim ZY, Goh KC, Kok JWK, Te SH, Tran NH, et al. Novel cyanotoxin-producing Synechococcus in tropical lakes. Water Res. 2021;192: 116828.

    Article  CAS  PubMed  Google Scholar 

  57. Kubickova B, Babica P, Hilscherová K, Šindlerová L. Effects of cyanobacterial toxins on the human gastrointestinal tract and the mucosal innate immune system. Environ Sci Eur. 2019;31:1–27.

    Article  CAS  Google Scholar 

  58. Pietrucci D, Teofani A, Milanesi M, Fosso B, Putignani L, Messina F, et al. Machine learning data analysis highlights the role of parasutterella and alloprevotella in autism spectrum disorders. Biomedicines. 2022;10:1–22.

    Article  Google Scholar 

  59. Dinh DM, Ramadass B, Kattula D, Sarkar R, Braunstein P, Tai A, et al. Longitudinal analysis of the intestinal microbiota in persistently stunted young children in south India. PLoS ONE. 2016;11:1–17.

    Article  Google Scholar 

  60. Zambruni M, Ochoa TJ, Somasunderam A, Cabada MM, Morales ML, Mitreva M, et al. Stunting is preceded by intestinal mucosal damage and microbiome changes and is associated with systemic inflammation in a cohort of Peruvian infants. Am J Trop Med Hyg. 2019;101:1009–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Gatya M, Fibri DLN, Utami T, Suroto DA, Rahayu ES. Gut microbiota composition in undernourished children associated with diet and sociodemographic factors: a case–control study in Indonesia. Microorganisms. 2022;10:1748.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Velly H, Britton RA, Preidis GA. Mechanisms of cross-talk between the diet, the intestinal microbiome, and the undernourished host. Gut Microbes. 2017;8:98–112.

    Article  PubMed  Google Scholar 

  63. Dubin K, Pamer EG. Enterococci and their interactions with the intestinal microbiome. Microbiol Spectr. 2018;611:309–30.

    Google Scholar 

  64. Radwan S, Gilfillan D, Eklund B, Radwan HM, El Menofy NG, Lee J, et al. A comparative study of the gut microbiome in Egyptian patients with Type I and Type II diabetes. PLoS ONE. 2020;15:1–17.

    Article  Google Scholar 

  65. De D, Nayak T, Chowdhury S, Dhal PK. Insights of host physiological parameters and gut microbiome of Indian type 2 diabetic patients visualized via metagenomics and machine learning approaches. Front Microbiol. 2022;13:1–15.

    Article  Google Scholar 

  66. Heeney DD, Gareau MG, Marco ML. Intestinal Lactobacillus in health and disease, a driver or just along for the ride? Curr Opin Biotechnol. 2018;176:139–48.

    Google Scholar 

  67. Tanno H, Fujii T, Ose R, Hirano K, Tochio T, Endo A. Characterization of fructooligosaccharide-degrading enzymes in human commensal Bifidobacterium longum and Anaerostipes caccae. Biochem Biophys Res Commun. 2019;518:294–8.

    Article  CAS  PubMed  Google Scholar 

  68. Louis P, Duncan SH, Sheridan PO, Walker AW, Flint HJ. Microbial lactate utilisation and the stability of the gut microbiome. Gut Microbiome. 2022;3:1–16.

    Article  Google Scholar 

  69. Olsson LM, Boulund F, Nilsson S, Khan MT, Gummesson A, Fagerberg L, et al. Dynamics of the normal gut microbiota: a longitudinal one-year population study in Sweden. Cell Host Microbe. 2022;30:726-739.e3.

    Article  CAS  PubMed  Google Scholar 

  70. Robertson RC, Manges AR, Finlay BB, Prendergast AJ. The human microbiome and child growth – first 1000 days and beyond. Trends Microbiol. 2019;27:131–47.

    Article  CAS  PubMed  Google Scholar 

  71. Charlie GB, Eric GP. Microbiota-mediated colonization resistance against intestinal pathogens. Nat Rev Immunol. 2014;13:790–801.

    Google Scholar 

  72. Acosta GJ, Vigo NI, Durand D, Riveros M, Arango S, Zambruni M, et al. Diarrheagenic Escherichia coli: prevalence and pathotype distribution in children from peruvian rural communities. Am J Trop Med Hyg. 2016;95:574–9.

    Article  PubMed  PubMed Central  Google Scholar 

  73. Das R, Palit P, Haque MA, Mahfuz M, Faruque ASG, Ahmed T. Site specific incidence rate of virulence related genes of enteroaggregative Escherichia coli and association with enteric inflammation and growth in children. Sci Rep. 2021;11:1–10.

    Article  Google Scholar 

  74. Zhou HL, Bessey T, Wang SM, Mo ZJ, Barclay L, Wang JX, et al. Burden and etiology of moderate and severe diarrhea in children less than 5 years of age living in north and south of China: prospective, population-based surveillance. Gut Pathog. 2021;13:1–11.

    Article  Google Scholar 

  75. Anderson JD, Bagamian KH, Muhib F, Amaya MP, Laytner LA, Wierzba T, et al. Burden of enterotoxigenic Escherichia coli and shigella non-fatal diarrhoeal infections in 79 low-income and lower middle-income countries: a modelling analysis. Lancet Glob Heal. 2019;7:e321–30.

    Article  Google Scholar 

  76. Haque MA, Nasrin S, Palit P, Das R, Wahid BZ, Gazi MA, et al. Site-specific analysis of the incidence rate of Enterotoxigenic Escherichia coli infection elucidates an association with childhood stunting, wasting, and being underweight: a secondary analysis of the MAL-ED birth cohort. Am J Trop Med Hyg. 2023;108:1192–200.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Bona M, Medeiros PH, Santos AK, Freitas T, Prata M, Veras H, et al. Virulence-related genes are associated with clinical and nutritional outcomes of Shigella/Enteroinvasive Escherichia coli pathotype infection in children from Brazilian semiarid region: a community case-control study. Int J Med Microbiol. 2019;309:151–8.

    Article  CAS  PubMed  Google Scholar 

  78. Vonaesch P, Tondeur L, Breurec S, Bata P, Nguyen LBL, Frank T, et al. Factors associated with stunting in healthy children aged 5 years and less living in Bangui (RCA). PLoS ONE. 2017.

    Article  PubMed  PubMed Central  Google Scholar 

  79. George CM, Burrowes V, Perin J, Oldja L, Biswas S, Sack D, et al. Enteric Infections in young children are associated with environmental enteropathy and impaired growth. Trop Med Int Heal. 2018;23:26–33.

    Article  CAS  Google Scholar 

  80. Zha L, Garrett S, Sun J. Salmonella infection in chronic inflammation and gastrointestinal cancer. Diseases. 2019;7:28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Ledwaba SE, Kabue JP, Barnard TG, Traore AN, Potgieter N. Enteric pathogen co-infections in the paediatric population from rural communities in the Vhembe District, South Africa. South Afr J Child Heal. 2018;12:170–4.

    Article  Google Scholar 

  82. Shrivastava AK, Kumar S, Mohakud NK, Suar M, Sahu PS. Multiple etiologies of infectious diarrhea and concurrent infections in a pediatric outpatient-based screening study in Odisha. India Gut Pathog. 2017;9:1–12.

    CAS  Google Scholar 

  83. Ogston SA, Lemeshow S, Hosmer DW, Klar J, Lwanga SK. Adequacy of sample size in Health Studies. Biometrics. 1991;47:347.

    Article  Google Scholar 

  84. Ministry of Health. Indonesia Basic Health Research 2019. 2019. Accesses 15 Jan 2020.

  85. Sigma Aldrich. Technical Bulletin EB-127TC. 2010. Accessed 20 Jun 2020.

  86. Novogene. Primers applied to Amplicon Metagenomic Sequencing. 2020. Accessed 26 May 2020.

  87. Chen Z, Hui PC, Hui M, Yeoh YK, Wong PY, Chan MCW, et al. Impact of preservation method and 16S rRNA hypervariable region on gut microbiota profiling. mSystems. 2019.

    Article  PubMed  PubMed Central  Google Scholar 

  88. Pfaffl MW. Real-time PCR. Dorak MT, editor. New York: Taylor and Francis Group; 2007.

  89. R Core Team. A language and environment for statistical computing. R Foundation for Statistical Computing. 2023. Accessed 15 Jan 2023.

  90. RStudio Team. RStudio: Integrated development for R. Boston: RStudio. PBC; 2020.

Download references


We would like to thank the late Prof. Dr. Debbie Sofie Retnoningrum who had contributed tremendously to the finalization of the research concepts.


This work was funded by Lembaga Pengelola Dana Pendidikan (LPDP) through BUDI-DN scholarship for a Doctoral Program. The Bacterial 16S rDNA amplicon sequencing was also partially funded by ASAHI Glass Foundation through ASAHI Overseas Grant 2022 on behalf of Dr. Lucy Sasongko; however, the funder has no role in manuscript preparation and publication.

Author information

Authors and Affiliations



TR, CR, and AA contributed to conceptualization and study design. TR managed funding acquisition, sample collection, laboratory works and analyses, and the first draft of the manuscript. LS supervised the sample collection and contributed to funding acquisition. CR and AA supervised the laboratory works and analyses. LS, CR, and AA performed manuscript review and editing. All authors read and approve the final manuscript.

Corresponding author

Correspondence to Lucy Sasongko.

Ethics declarations

Ethics approval and consent to participate

The protocol of this study was approved by Research Ethics Committee of Faculty of Medicine, Universitas Syiah Kuala (dossier No.050/EA/FK-RSUDZA/2021).

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Fig. S1.

Unweighted UniFrac model for beta diversity. The metric was visualized using QIIME2 emperor. Fig. S2. Rarefaction curve of alpha diversity. The minimum sample depth was used as rarefaction depth (51,636 reads). The curve were plotted using Ampvis2 package in R and visualized in RStudio. Table S1. The reference gene (rrs) expression among observed groups. Table S2. Primers used in this study.

Additional file 2

. Taxonomy table of classified taxa from 16S rDNA amplicon sequencing. The table was constructed and converted using QIIME2. Enteric infections-related taxa were bolded.

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rinanda, T., Riani, C., Artarini, A. et al. Correlation between gut microbiota composition, enteric infections and linear growth impairment: a case–control study in childhood stunting in Pidie, Aceh, Indonesia. Gut Pathog 15, 54 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: