- Open Access
Microbiota profile in new-onset pediatric Crohn’s disease: data from a non-Western population
Gut Pathogensvolume 10, Article number: 49 (2018)
The role of microbiota in Crohn’s disease (CD) is increasingly recognized. However, most of the reports are from Western populations. Considering the possible variation from other populations, the aim of this study was to describe the microbiota profile in children with CD in Saudi Arabia, a non-Western developing country population.
Significantly more abundant genera in children with CD included Fusobacterium, Peptostreptococcus, Psychrobacter, and Acinetobacter; whereas the most significantly-depleted genera included Roseburia, Clostridium, Ruminococcus, Ruminoclostridium, Intestinibacter, Mitsuokella, Megasphaera, Streptococcus, Lactobacillus, Turicibacter, and Paludibacter. Alpha diversity was significantly reduced in stool (p = 0.03) but not in mucosa (p = 0.31). Beta diversity showed significant difference in community composition between control and CD samples (p = 0.03).
In this developing country, we found a pattern of microbiota in children with CD similar to Western literature, suggesting a role of recent dietary lifestyle changes in this population on microbiota structure.
Crohn’s disease (CD) is the most common phenotype of inflammatory bowel disease (IBD). Although the incidence is highest in Western populations, increasing time trend is reported worldwide in adults and children [1,2,3]. Despite extensive research, the causes of all phenotypes of IBD remain unknown. However, a multifactorial theory is most likely. In a genetically susceptible individual, environmental factors trigger uncontrolled inflammation . Diet and microbiota are the most likely environmental factors, acting as triggers. It was found that there was an increased risk of CD with a high intake of polyunsaturated fatty acids, omega-6 fatty acids, saturated fats, and meat, but there was a decreased risk with a high intake of dietary fiber, omega-3 fatty acids, vegetables, and fruits [5,6,7,8,9,10]. The role of microbiota in IBD in general, and CD in particular, has been increasingly recognized. Several studies documented reduced diversity of the microbial community and defined associations of certain taxa with CD, suggesting a role of beneficial and harmful microbes [11,12,13,14].
Almost all the literature on microbiota in IBD is regarding Western populations, who have well-defined environmental, cultural, and dietary lifestyles, which are different from populations in developing countries. Since IBD is newly recognized in these populations, these are commonly referred to as “new populations” with IBD. The study of the characteristics of IBD in these populations may increase our understanding of this condition.
In previous reports, we have defined the incidence and clinical profile of pediatric CD in Saudi Arabia [15, 16]. The objective of this report is to describe the bacterial microbiota profile in a cohort of newly-diagnosed treatment-naïve children in Saudi Arabia, a non-Western developing country population in comparison with other populations reported in relevant literature.
All children were Saudi nationals (17 with CD and 18 controls). The median (range) age was 15 (7.3–17.8) years for the children with CD and 16.3 (3.9–18.6) years for controls. Gender distribution indicated that 11/17 (65%) of the CD and 12/18 (67%) of the control subjects were males.
Dietary history revealed similar prevalence of breastfeeding in children with CD (13/17; 76%) and controls (13/18; 72%). The quality of family food intake by children estimated by mothers was very good to good in 6/17 (35%) and 9/18 (50%) in children with CD and controls respectively. Daily fast food consumption was 10/17 (58%) and 9/18 (50%), daily sweetened gaseous drink in 11/17 (65%) and 9/18 (50%) in children with CD and controls respectively. Daily fruit consumption was 1/17 (6%) and 2/18 (11%) and vegetable 3/17 (18%) and 10/18 (56%) in children with CD and controls respectively. In summary, compared to controls, children with CD tended to eat less family food, more fast food and sweetened gaseous drinks, but less vegetable and fruit consumption.
The clinical presentation in the 17 children with CD included abdominal pain in 15 (88%), diarrhea in 14 (82%), weight loss in 13 (76%), blood in stools in 8 (47%), and perianal disease in 3 (18%) children. At presentation, laboratory tests revealed anemia in 8 (47%), thrombocytosis in 6 (35%), elevated erythrocyte sedimentation rate (ESR) in 9 (53%), and high C reactive protein (CRP) in 6 (35%) children. At diagnosis, CD locations were ileal (L1) in 4 (24%) and Ileocolonic (L3) in 13 (76%) children while CD behavior was non-constricting non-penetrating (B1) in 13 (75%) and constricting (B2) in 4 (24%) children.
The clinical presentation of the 18 children classified as controls included recurrent abdominal pain in 9 (50%) children, finally diagnosed as functional abdominal pain, diarrhea in 5 (28%) children, finally diagnosed as nonspecific diarrhea, and 4 (22%) children with rectal bleeding, finally diagnosed as juvenile polyps. Hemoglobin, platelets, ESR, and CRP were normal in all controls.
The CD-associated taxa for the family, genus, and species phylogenetic levels in stool and mucosa are presented in Tables 1 and 2 respectively. Significantly-abundant taxa in children with CD compared with controls at the genus level included Fusobacterium, Peptostreptococcus, Psychrobacter, and Acinetobacter, and at the species level, Fusobacterium nucleatum, Bacteroides clarus, and Psychrobacter pulmonis. The most significantly-depleted genera in children with CD included Roseburia, Clostridium, Ruminococcus, Ruminoclostridium, Intestinibacter, Mitsuokella, Megasphaera, Streptococcus, Lactobacillus, Turicibacter, and Paludibacter; whereas, significantly-depleted species included Roseburia inulivorans, Clostridium disporicum, Blaucia ruminococcus spp., Eubacterium seraeum, Intestinibacter bartelitii, Eubacterium spp., Streptococcus salivarius, Turicibacter spp., Bacteroides xylanolyticus, Clostridium perfringens, and Bifidobacterium catenulatum. It is to be noted that no significantly more abundant taxa were found in mucosal samples. By contrast, a large number of taxa were depleted from stool and mucosa samples as detailed in Tables 1 and 2. The direction of change of most taxa (gain or loss) is depicted in Fig. 1 and the rank abundance distribution of the 20 most abundant genera in stool and mucosa samples is illustrated in Fig. 2.
Alpha diversity, as measured by the Shannon index, is shown in Fig. 3, indicating significantly-reduced alpha diversity in the stool of children with CD compared to controls (p = 0.03); whereas, the difference in CD mucosa was not significant (p = 0.32).
Beta diversity as evaluated by the Bray–Curtis distance indicated a statistically significant community dissimilarity between control and CD samples in stool (p = 0.03).
Abundance and diversity in inflamed and un inflamed mucosa
Comparison of species abundance between inflamed and uninflamed mucosa of children with CD for the four most different species (Bacteroides nordii, Escherichia coli, Eisenbergiella tayi, Bacteroides caccae), indicated no significant difference (p > 0.9). In addition, there was no significant difference in alpha diversity between inflamed and uninflamed mucosa in children with CD (p = 0.31).
The national incidence of pediatric IBD in Saudi Arabia (0.47/105), including CD (0.27/105), has been reported recently, indicating a lower but steady increase in incidence similar to that in the Western literature . In addition, the clinical, laboratory, endoscopic, and histopathologic characteristics have been reported, indicating similar patterns to descriptions from Western countries [16, 17]. In Saudi Arabia, marked socioeconomic improvement which led to improved education, nutrition, and health care, was accompanied by changes in lifestyle. For example, in a larger report on Saudi children with CD, only 10% and 30% consumed fruits daily and twice weekly respectively; whereas about 50% and 30% consumed fast food and sweetened soft drinks daily and twice weekly respectively . This dietary lifestyle pattern indicating less consumption of fruits and high consumption of fast food and sweetened soft drinks is similar to descriptions in the Western literature [19, 20].
IBD in Saudi Arabia evolved from occurring rarely, to a commonly-diagnosed condition with increasing incidence, suggesting a role of recent changes in environment including dietary lifestyle. It is well known that dietary components acting directly or through alteration of intestinal microbiota have a significant role in triggering inflammation [21,22,23,24].
This report, to our knowledge, is the first description of microbiota profile in newly-diagnosed treatment-naïve children with CD from a non-Western population. We identified a large number of taxa in the CD fecal and mucosal samples from order to species levels in the phylogenetic tree. However, it should be noted that taxonomic species designations based upon 16 s are tentative assignments and caution is advised in interpretations related to species classification. In view of the variation of the microbiome between newly-diagnosed and treated patients with IBD , comparison will be mainly with the few most similar reports of newly-diagnosed and treatment-naïve children with CD.
All of the 20 most abundant taxa from the order to the species levels in this study, have been reported in the Western literature [4, 14, 26, 27]. However, the significance of associations of some taxa with CD, contrasts with reports from Western populations. For example, Enterobacteriaceae, reported significantly more abundant in CD , was not found significantly associated with CD in our study. Similarly, Faecalibacterium prausnitzii, which is reported as significantly-depleted bacteria with possible anti-inflammatory properties , was not found to be significantly associated with CD in our samples. Variations in the significance of associations exists between studies even from the same Western populations. For example, at least two reports found Faecalibacterium prausnitzii to be more abundant in mucosal samples of children with CD, contrasting with other reports of the protective role of this bacterium in patients with CD [29, 30]. These observations are in line with the well-known variability of microbiota both within and between subjects in the same population. Finally, it is important to note that associations described in this study as well as in the literature, do not imply functional or causal effects. Specifically, it is still unclear whether changes in microbiota in children with CD were the cause or the result of inflammation.
In this study, alpha diversity was reduced in mucosa and stools of children with CD relative to controls, a finding similar to that in the literature . However, this reduction was statistically significant only in stools (p = 0.03) and not in tissue samples (p = 0.32). Beta diversity in our cohort indicated an overall statistically-significant distance difference between CD and control samples (p = 0.03), a finding similar to reports from Western populations .
Lack of difference between inflamed and uninflamed CD mucosa
In this study, the lack of significant difference in the four most abundant bacterial species and alpha diversity between sites with inflamed and uninflamed mucosa are consistent with most reports and suggest an unlikely role of bacteria in the pathogenesis of the patchy distribution of lesions in CD [32, 33]. However, one report suggested that uninflamed tissue forms an intermediate bacterial population between controls and inflamed tissues , and another reported a significant difference in microbial community structure between inflamed and uninflamed mucosal sites, but there was great variation between individuals, suggesting no obvious bacterial signature that was positively associated with the inflamed gut . It appears, therefore, that our findings of no significant difference in the bacterial community between inflamed and uninflamed mucosa in children with CD are consistent with most reports.
In this study from a developing country population, the finding of a microbiota profile similar to that in Western populations was unexpected in view of different culture and lifestyle. However, recent changes to a more westernized dietary lifestyle, affecting microbiota structure, explains at least in part this similarity.
The most important limitation is the small sample size. However, Crohn’s disease is evolving in this part of the world and the characterization of the microbiota associated with new onset Crohn’s is unique.
In this developing country, we found a pattern of microbiota in children with CD similar to Western literature, suggesting an effect of recent dietary lifestyle changes on microbiota structure. This report suggests a possible role of dietary lifestyle related to alteration of microbiota and the increasing incidence of CD in the Saudi population.
The study population
In addition to controls, the study population included all children diagnosed with CD according to standard guidelines . The children referred for colonoscopy were enrolled prospectively. Two hospitals participated in the study. These were King Khalid University Hospital, King Saud University (a free-access primary and tertiary care hospital), and Al Mofarreh Polyclinics (a private gastroenterology institution). Demographic information, socioeconomic family status, nutritional history, drug history, history of the present illness, past medical and surgical history, including any medications, physical examination, laboratory, imaging, endoscopic, and histopathological findings were recorded at presentation. Controls were enrolled if they had no evidence of IBD or other causes of inflammation proven by endoscopy and histopathology.
Sample collection, storage, and processing
In view of the known variations in the microbial community along the gastrointestinal (GI) tract [37,38,39,40,41], mucosal samples were collected from the ileum and different colonic sites to minimize the effects of these variations. Samples were collected from 17 children with CD and 18 controls. A total of 44 tissue samples from the children with CD (8 from the Ileum, 6 from each of the cecum, ascending, transverse, descending, sigmoid colon and rectum) and 14 from controls (6 from the Ileum, 3 transverse colon, 2 sigmoid colon and 3 from the rectum). For logistic reasons, mucosal samples were not taken from the ileum and all colonic segments of each subject. Similarly, not all subjects gave stool samples. A total of 20 stool samples were collected from children with CD (10) and controls (10) before bowel preparation (75%), or from the first stool passed after the start of bowel preparation (25%), to minimize washout effects . All samples were collected in cryovials (no fixatives or stabilizers), immediately placed in ice, transported to the laboratory, and stored at − 80 °C within 5 to 20 min. The average storage duration before analysis was 3 years. At the time of microbiota analysis, all samples were shipped in dry ice by express mail to the USA (MR DNA, Shallowater, TX, USA). The samples were received frozen in about 36 h.
DNA extraction and sequencing methods
DNA was extracted using the Mobio Powersoil Kit as per the manufacturer’s instructions (MOBIO, Carlsbad, CA, USA). Amplicon sequencing (bTEFAP®) was performed at MR DNA (Shallowater, TX, USA) and used for bacterial analysis . The primers 515F GTGCCAGCMGCCGCGGTAA and 806R GGACTACHVGGGTWTCTAAT were used to evaluate the microbial ecology of swabs on the Illumina MiSeq with methods based upon the bTEFAP®. A single-step 28-cycle polymerase chain reaction (PCR) with the HotStarTaq Plus Master Mix Kit (Qiagen, Valencia, CA, USA) was employed under the following conditions: 94 °C for 3 min, followed by 28 cycles of 94 °C for 30 s, 53 °C for 40 s, and 72 °C for 1 min; after this, a final elongation step at 72 °C for 5 min was performed. Following PCR, all amplicon products from different samples were mixed in equal concentrations and purified using Agencourt Ampure Beads (Agencourt Bioscience Corporation, MA, USA). Samples were sequenced utilizing Illumina MiSeq Chemistry following the manufacturer’s protocols. The Q25 sequence data derived from the sequencing process was processed using a standard analysis pipeline (http://www.mrdnalab.com; MR DNA, Shallowater, TX, USA). Paired sequences were merged and depleted of barcodes and primers, then short sequences < 150 base pairs (bp) were removed, sequences with ambiguous base calls were removed, and sequences with homopolymer runs exceeding 6 bp were removed. Sequences were then denoised and chimeras were removed. Operational taxonomic units (OTUs) were defined after removal of singleton sequences, clustering at 3% divergence (97% similarity) [43,44,45,46]. OTUs were then taxonomically classified using the Nucleotide Basic Local Alignment Search Tool against a 16 s National Center for Biotechnology Information (NCBI)-derived database (http://www.ncbi.nlm.nih.gov, http://rdp.cme.msu.edu), and compiled into each taxonomic level into both ‘counts’ and ‘percentage’ files. Counts files contain the actual number of sequences and percentage files contain the relative (proportional) percentage of sequences within each sample, which map to the designated taxonomic classification.
The analysis was performed using Python and R software [47, 48]. To increase statistical power, taxa with low representation in the samples were excluded from the analysis. Specifically, we excluded samples below 1000 reads as well as taxa absent from more than 50% of both CD and control samples. Custom functions implementing the permutation test were written to detect the taxa whose abundances were significantly different between two sample categories, e.g., CD and control or inflamed and uninflamed. When more than one sample was available from the same patient for the analysis, the log-relative abundances from these samples were averaged. We performed all statistical analysis on log-transformed data after adding pseudo counts of 1 read for each taxonomic group.
To understand which members of the bacterial microbiota might contribute to CD, we examined the difference in microbial abundance between CD stool and controls, CD mucosa and controls, and inflamed and uninflamed mucosa in CD. First, we compared uninflamed and inflamed mucosa in children with CD and found no significant changes in microbial abundances (p > 0.9). Given this lack of difference, all CD mucosal samples were included in the analysis. Associations were determined based on the difference in the mean log-relative abundance. Statistical significance was assessed via a permutation test (exact Fisher’s test) followed by a correction for multiple hypothesis testing. Specifically, the permutation test yielded raw, uncorrected p-values, which were corrected to q-values following the Benjamini–Hochberg procedure to account for the false discovery rate (FDR) . Although less significantly-associated taxa may be biologically important, we reported only statistically-significant associations when the corrected FDR-corrected p-value was < 0.05.
Diversity analysis was used to study the richness of taxa and the evenness of habitat composition, as well as the community dissimilarity between samples. Alpha diversity, a measure of taxa richness was evaluated by the Shannon index. This measure quantifies the number of taxa and their representative proportions in the habitat; a high alpha diversity indicates that there is a high number of taxa with similar abundance. The difference in diversity between CD mucosa and controls, stool CD and controls, or inflamed and uninflamed mucosa were analyzed. The sample-wise difference in community composition (Beta diversity) was quantified by the Bray–Curtis dissimilarity, which accounts for both patterns of presence-absence of taxa and changes in their relative abundance between samples.
The beta diversity separations were analyzed by the ANOSIM or analysis of (dis)similarity. The ANOSIM statistic compares the mean of ranked dissimilarities between groups to the mean of ranked dissimilarities within groups. The significance of the statistic was determined by an exact permutation test.
Kappelman MD, Moore KR, Allen JK, Cook SF. Recent trends in the prevalence of Crohn’s disease and ulcerative colitis in a commercially insured US population. Dig Dis Sci. 2009;58:519–25.
Molodecky NA, Soon IS, Rabi DM, Ghali WA, Ferris M, Chernoff G, et al. Increasing incidence and prevalence of the inflammatory bowel diseases with time, based on systematic review. Gastroenterology. 2012;142:46–54.
Benchimol EI, Fortinsky KJ, Gozdyra P, den Heuvel M, Limbergen JV, Griffiths AM. Epidemiology of pediatric inflammatory bowel disease: a systematic review of international trends. Inflamm Bowel Dis. 2011;17:423–39.
Green C, Elliott L, Beaudoin C, Bernstein CN. A population-based ecologic study of inflammatory bowel disease: searching for etiologic clues. Am J Epidemiol. 2006;164:615–23.
Hou JK, Abraham B, El-Serag H. Dietary intake and risk of developing inflammatory bowel disease: a systematic review of the literature. Am J Gastroenterol. 2011;106:563–73.
Ananthakrishnan AN, Khalili H, Konijeti GG, Higuchi LM, de Silva P, Korzenik JR, et al. A prospective study of long-term intake of dietary fiber and risk of Crohn’s disease and ulcerative colitis. Gastroenterology. 2013;145:970–7.
Li F, Liu X, Wang W, Zhang D. Consumption of vegetables and fruit and the risk of inflammatory bowel disease: a meta-analysis. Eur J Gastroenterol Hepatol. 2015;27:623–30.
Ananthakrishnan AN, Khalili H, Konijeti GG, Higuchi LM, de Silva P, Fuchs CS, et al. Long-term intake of dietary fat and risk of ulcerative colitis and Crohn’s disease. Gut. 2014;63:776–84.
Chan SSM, Luben R, Olsen A, Tjonneland A, Kaaks R, Lindgren S, et al. Association between high dietary intake of the n3 polyunsaturated fatty acid docosahexaenoic acid and reduced risk of Crohn’s disease. Aliment Pharmacol Ther. 2014;39:834–42.
Lewis JD, Abreu MT. Diet as a trigger or therapy for inflammatory bowel diseases. Gastroenterology. 2017;152:398–414.
Shaw KA, Bertha M, Hofmekler T, Chopra P, Vatanen T, Srivatsa A, et al. Dysbiosis, inflammation, and response to treatment: a longitudinal study of pediatric subjects with newly diagnosed inflammatory bowel disease. Genome Med. 2016;8:75.
Maukonen J, Kolho KL, Paasela M, Honkanen J, Klemetti P, Vaarala O, et al. Altered fecal microbiota in paediatric inflammatory bowel disease. J Crohns Colitis. 2015;9:1088–95.
Kolho KL, Korpela K, Jaakkola T, Pichai MV, Zoetendal EG, Salonen A, et al. Fecal microbiota in pediatric inflammatory bowel disease and its relation to inflammation. Am J Gastroenterol. 2015;110:921–30.
Gevers D, Kugathasan S, Denson LA, Vázquez-Baeza Y, Van Treuren W, Ren B, et al. The treatment-naive microbiome in new-onset Crohn’s disease. Cell Host Microbe. 2014;15:382–92.
El Mouzan MI, Saadah O, Al-Saleem K, Al Edreesi M, Hasosah M, Alanazi A, et al. Incidence of pediatric inflammatory bowel disease in Saudi Arabia: a multicenter national study. Inflamm Bowel Dis. 2014;20:1085–90.
Saadah OI, El Mouzan M, Al Mofarreh M, Al Mehaidib A, Al Edreesi M, Hasosah M, Al-Hussaini A, et al. Characteristics of pediatric Crohn’s disease in Saudi children: a multicenter national study. Gastroenterol Res Pract. 2016;2016:7403129.
Al-Hussaini A, El Mouzan M, Hasosah M, Al-Mehaidib A, ALSaleem K, Saadah OI, et al. Clinical pattern of early-onset inflammatory bowel disease in Saudi Arabia: a multicenter national study. Inflamm Bowel Dis. 2016;22:1961–70.
El Mouzan M, Al Mofarreh M, Al Sarkhy A, Assiri A, Hamed Y. Pre-Illness diet as risk factor in pediatric inflammatory bowel disease in Saudi Arabia. Saudi J Gastroenterol. 2017;23:287–90.
Grimm KA, Kim SA, Yaroch AL, Scanlon KS. Fruit and vegetable intake during infancy and early childhood. Pediatrics. 2014;134(Suppl 1):S63–9.
Savige GS, Ball K, Worsley A, Crawford D. Food intake patterns among Australian adolescents. Asia Pac J Clin Nutr. 2007;16:738–47.
Laitinen K, Collado MC, Isolauri E. Early nutritional environment: focus on health effects of microbiota and probiotics. Benef Microbes. 2010;1:383–90.
Kau AL, Ahern PP, Griffin NW, Goodman AL, Gordon JI. Human nutrition, the gut microbiome and the immune system. Nature. 2011;474:327–36.
Clemente JC, Ursell LK, Parfrey LW, Knight R. The impact of the gut microbiota on human health: an integrative view. Cell. 2012;148:1258–70.
Rapozo D, Bernardazzi C, Pereira de Souza H. Diet and microbiota in inflammatory bowel disease: the gut in disharmony. World J Gastroenterol. 2017;23:2124–40.
Bibiloni R, Mangold M, Madsen KL, Fedorak RN, Tannock GW. The bacteriology of biopsies differs between newly diagnosed, untreated, Crohn’s disease and ulcerative colitis patients. J Med Microbiol. 2006;55:1141–9.
Morgan XC, Tickle TL, Sokol H, Gevers D, Devaney KL, Ward DV, et al. Dysfunction of the intestinal microbiome in inflammatory bowel disease and treatment. Genome Biol. 2012;13:R79.
Papa E, Docktor M, Smillie C, Weber S, Preheim SP, Gevers D, et al. Non-invasive mapping of the gastrointestinal microbiota identifies children with inflammatory bowel disease. PLoS ONE. 2012;7:e39242.
Sokol H, Pigneur B, Watterlot L, Lakhdari O, Bermúdez-Humarán LG, Gratadoux JJ, et al. Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. PNAS. 2008;105:16731–6.
Hansen R, Russell RK, Reiff C, Louis P, McIntosh F, Berry SH, et al. Microbiota of de-novo pediatric IBD: increased Faecalibacterium prausnitzii and reduced bacterial diversity in Crohn’s but not in ulcerative colitis. Am J Gastroenterol. 2012;107:1913–22.
Assa A, Butcher J, Li J, Sherman PM, Muise AM, Stintzi A, et al. Mucosa-associated ileal microbiota in new onset pediatric Crohn’s disease. Inflamm Bowel Dis. 2016;22:1533–9.
Manichanh C, Rigottier-Gois L, Bonnaud E, Gloux K, Pelletier E, Frangeul L, et al. Reduced diversity of faecal microbiota in Crohn’s disease revealed by a metagenomic approach. Gut. 2006;55:205–511.
Seksik P, Lepage P, de la Cochetière MF, Bourreille A, Sutren M, Galmiche JP, et al. Search for localized dysbiosis in Crohn’s disease ulcerations by temporal temperature gradient gel electrophoresis of 16S rRNA. J Clin Microbiol. 2005;43:4654–8.
Vasquez N, Mangin I, Lepage P, Seksik P, Duong JP, Blum S, et al. Patchy distribution of mucosal lesions in ileal Crohn’s disease is not linked to differences in the dominant mucosa-associated bacteria: a study using fluorescence in situ hybridization and temporal temperature gradient gel electrophoresis. Inflamm Bowel Dis. 2007;13:684–92.
Sepehri S, Kotlowski R, Bernstein CN, Krause DO. Microbial diversity of inflamed and noninflamed gut biopsy tissues in inflammatory bowel disease. Inflamm Bowel Dis. 2007;13:675–83.
Walker AW, Sanderson JD, Churcher C, Parkes GC, Hudspith BN, Rayment N, et al. High-throughput clone library analysis of the mucosa-associated microbiota reveals dysbiosis and differences between inflamed and noninflamed regions of the intestine in inflammatory bowel disease. BMC Microbiol. 2011;11:7.
Levine A, Koletzko S, Turner D, Escher JC, Cucchiara S, de Ridder L, et al. European society of pediatric gastroenterology, hepatology, and nutrition. ESPGHAN revised Porto criteria for the diagnosis of inflammatory bowel disease in children and adolescents. J Pediatr Gastroenterol Nutr. 2014;58:795–806.
Swidsinski A, Ladhoff A, Pernthaler A, Swidsinski S, Loening-Baucke V, Ortner M, et al. Mucosal flora in inflammatory bowel disease. Gastroenterology. 2002;122:44–54.
Hu S, WangY Lichtenstein L, TaoY Musch MW, Jabri B, et al. Regional differences in colonic mucosa-associated microbiota determine the physiological expression of host heashock proteins. Am J Physiol Gastrointest Liver Physiol. 2010;299:G1266–75.
WangY Antonopoulos D, Zhu X, Harrell L, Hanan I, Alverdy J, et al. Laser capture microdissection and metagenomic analysis of intact mucosa-associated microbial communities of human colon. Appl Microbiol Biotechnol. 2010;88:1333–42.
WangY Devkota S, Musch MW, Jabri B, Nagler C, Antonopoulos DA, et al. Regional mucosa—associated microbiota determine physiological expression of TLR2 and TLR4 in murin colon. PLoS ONE. 2010;5:e13607.
Shobar RM, Velineni S, Keshavarzian A, Swanson G, DeMeo MT, Melson JE, et al. The effects of bowel preparation on microbiota-related metrics differ in health and in inflammatory bowel disease and for the mucosal and luminal microbiota compartments. Clin Transl Gastroenterol. 2016;7:e143.
Dowd SE, Callaway TR, Wolcott RD, Sun Y, McKeehan T, Hagevoort RG, et al. Evaluation of the bacterial diversity in the feces of cattle using 16S rDNA bacterial tag-encoded FLX amplicon pyrosequencing [bTEFAP]. BMC Microbiol. 2008;8:125.
Dowd SE, Sun Y, Wolcott RD, Domingo A, Carroll JA. Bacterial tag encoded FLX amplicon pyrosequencing [bTEFAP] for microbiome studies: bacterial diversity in the ileum of newly weaned Salmonella-infected pigs. Foodborne Pathog Dis. 2008;5:459–72.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.
Eren AM, Zozaya M, Taylor CM, Dowd SE, Martin DH, Ferris MJ. Exploring the diversity of Gardnerella virginals in the genitourinary tract microbiota of monogamous couples through subtle nucleotide variation. PLoS ONE. 2011;6:e26732.
Swanson KS, Dowd SE, Suchodolski JS, Middelbos IS, Vester BM, Barry KA, et al. Phylogenetic and gene-centric metagenomics of the canine intestinal microbiome reveals similarities with humans and mice. ISME J. 2011;5:639–49.
Core R, Team R. A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2015.
McMurdie PJ, Holmes S. An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8:e61217.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B. 1995;1995:289–300.
MIE designed the study and wrote the manuscript. HSW participated in study design, reviewed and approved the final manuscript. AMA, AAS, and MAM, participated in patient enrollment and collection of data and approved the final manuscript. KSK and RM performed bioinformatics, statistical analysis and approved the final manuscript. SED performed microbial DNA extraction and sequencing. All authors read and approved the final manuscript.
The authors extend their appreciations to the Deanship of Scientific Research at King Saud University in Riyadh, Kingdom of Saudi Arabia for funding this work through Research Group No (RGP-1436-007). Harland Winter is supported in part by a philanthropic grant from Martin Schlaff and the B. Hasso Family Foundation. This work was also partially supported by a grant from Simons Foundation (No. 409704 to Kirill Korolev). Rajita Menon was partially supported by a Hariri Graduate Fellowship and a Hariri Research Award (to Kirill Korolev and Daniel Serge) both from Boston University. Data analysis was carried out on Shared Computing Cluster at Boston University.
The authors declare that they have no competing interests.
Availability of data and materials
The data that support the findings of this study are available upon reasonable request from the principal investigator (Professor Mohammad El Mouzan) and with permission of the Deanship of Scientific Research, King Saud University, Riyadh, The Kingdom of Saudi Arabia.
Consent for publication
Ethics approval and consent to participate
The protocol of this study was approved by the Institutional Review Board of the College of Medicine, King Saud University in Riyadh, Kingdom of Saudi Arabia (No: 10/2647/IRB). Signed informed consent was obtained from the guardians and/or patients before enrollment in the study.
This study was funded by the Deanship of Scientific Research, King Saud University, Riyadh, Kingdom of Saudi Arabia. Grant No. RGP-1436-007.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.