- Research article
- Open Access
Gastrointestinal microbial populations can distinguish pediatric and adolescent Acute Lymphoblastic Leukemia (ALL) at the time of disease diagnosis
BMC Genomicsvolume 17, Article number: 635 (2016)
An estimated 15,000 children and adolescents under the age of 19 years are diagnosed with leukemia, lymphoma and other tumors in the USA every year. All children and adolescent acute leukemia patients will undergo chemotherapy as part of their treatment regimen. Fortunately, survival rates for most pediatric cancers have improved at a remarkable pace over the past three decades, and the overall survival rate is greater than 90 % today. However, significant differences in survival rate have been found in different age groups (94 % in 1–9.99 years, 82 % in ≥10 years and 76 % in ≥15 years). ALL accounts for about three out of four cases of childhood leukemia. Intensive chemotherapy treatment coupled with prophylactic or therapeutic antibiotic use could potentially have a long-term effect on the resident gastrointestinal (GI) microbiome. The composition of GI microbiome and its changes upon chemotherapy in pediatric and adolescent leukemia patients is poorly understood. In this study, using 16S rRNA marker gene sequences we profile the GI microbial communities of pediatric and adolescent acute leukemia patients before and after chemotherapy treatment and compare with the microbiota of their healthy siblings.
Our study cohort consisted of 51 participants, made up of matched pediatric and adolescent patients with ALL and a healthy sibling. We elucidated and compared the GI microbiota profiles of patients and their healthy sibling controls via analysis of 16S rRNA gene sequencing data. We assessed the GI microbiota composition in pediatric and adolescent patients with ALL during the course of chemotherapy by comparing stool samples taken before chemotherapy with stool samples collected at varying time points during the chemotherapeutic treatment. The microbiota profiles of both patients and control sibling groups are dominated by members of Bacteroides, Prevotella, and Faecalibacterium. At the genus level, both groups share many taxa in common, but the microbiota diversity of the patient group is significantly lower than that of the control group. It was possible to distinguish between the patient and control groups based on their microbiota profiles. The top taxa include Anaerostipes, Coprococcus, Roseburia, and Ruminococcus2 with relatively higher abundance in the control group. The observed microbiota changes are likely the result of several factors including a direct influence of therapeutic compounds on the gut flora and an indirect effect of chemotherapy on the immune system, which, in turn, affects the microbiota.
This study provides significant information on GI microbiota populations in immunocompromised children and opens up the potential for developing novel diagnostics based on stool tests and therapies to improve the dysbiotic condition of the microbiota at the time of diagnosis and in the earliest stages of chemotherapy.
Acute lymphoblastic leukemia (ALL) is a malignant disease of the bone marrow in which early lymphoid precursors proliferate and replace the normal hematopoietic cells of the marrow . The lymphoid progenitor cells in patients with ALL are affected by the disease, leading to an impaired immune system typically observed at the time of diagnosis. ALL is the most common type of leukemia in children in the United States accounting for 26 % of all cancers in children up to 14 years of age and for 75 % of all pediatric leukemia cases . All children and adolescents diagnosed with ALL will undergo chemotherapy as part of their treatment plan, and their health can be severely compromised due to the treatment with chemotherapeutic drugs. Chemotherapy can damage healthy cells in the lining of the digestive system and gastrointestinal (GI) disturbances are often induced in response to chemotherapy. Chemotherapeutic and antibiotic treatment has a detrimental impact on the host microbial ecosystem, which is essential for host mucosal protection . The cytotoxic effects of these treatments lead to additional immunosuppression, which entails episodes of febrile neutropenia and potentially life-threatening bloodstream infections. Further, use of prophylactic and therapeutic antibiotics disrupts the GI microbiome’s ecological balance . Medium-term gastrointestinal health outcomes of chemotherapy and disruptions of the intestinal microbial ecology include vomiting, diarrhea, constipation and Clostridium difficile infection-associated diarrhea (CDAD). To varying degrees, survivors of children and adolescents treated for cancer experience a range of long-term growth and developmental, organ function, fertility and reproduction and psychological adverse outcomes. The growing population of survivors of pediatric cancer reflects a highly vulnerable group of individuals who will probably experience adverse health-related and quality-of-life outcomes during their subsequent lifetimes, as a result of their curative cancer treatment [5–8]. There is increasing evidence that the composition of the GI microbiome may affect, and is modulated by, the human immune system . Perturbed GI microbiomes have been associated with decreased immune competence, detrimental metabolic changes (e.g. obesity and malnutrition), susceptibility to GI infections and inflammatory syndromes [10, 11]. Likewise, several microbiome studies have examined the impact of antibiotics on the microbiota of individuals using 16S rRNA gene sequencing alone [12–16] or combined with metatranscriptomics . These studies have demonstrated that administration of antibiotics can perturb this microflora temporarily and in certain cases permanently. Antibiotic-mediated alteration of the gut microbiome can converts the global metabolic profile to one that favors C. difficile germination and growth . Furthermore, chemotherapeutics have a detrimental effect on the intestinal microbial composition, coinciding in time with the development of chemotherapy-induced mucositis . There is, however, a limited understanding of the GI microbiome composition of pediatric and adolescent leukemia patients and the impact, if any, of contemporary chemotherapeutic treatments. The rationale for this study is to define and compare the GI microbiota composition of pediatric and adolescent leukemia patients with their healthy sibling controls via analysis of high-throughput sequencing data. And also, to assess the changes in microbiota structure of pediatric and adolescent leukemia patients during chemotherapy by comparing the samples taken before and after chemotherapy at varying time points during the chemotherapeutic treatment. Our results provide significant information on GI microbiota composition in immunocompromised children and indicate that the baseline microbiota of immunocompromised children was substantially different from their healthy siblings. It creates the potential for the better management of GI and systemic complications associated with immunodeficiency and other disease conditions of this type. Furthermore, characterizing the GI microbiota dynamics following chemotherapy treatment will address what alterations happen to the GI microbiota during and following chemotherapy regimens and can correlate with response to treatment.
A cohort study was designed to assess the impact of chemotherapy on the GI microbiota of pediatric and adolescent patients diagnosed with acute B-cell leukemia. The study cohort consisted of 51 participants, made up of 23 matched patients and a healthy sibling and five unmatched patients. Five patients who did not have enrolled healthy siblings were also included in the cohort. Three subjects not complete the study: two withdrew and one subject was deceased. Subject demographics by age and gender are shown in Table 1. All study participants were enrolled in the Hyundai Cancer Institute, Children’s Hospital Orange County (CHOC Children’s), California, USA. Human subject protocol and consent forms were established, and approved by the Institutional Review Boards at CHOC Children’s and the J. Craig Venter Institute (JCVI). Stool samples were collected at the completion of each treatment stage during the patient's stay at the hospital, referred to as “sampling visits”. Samples marked “visit 1” were collected at the time of diagnosis before any chemotherapy was administered, and thus provided baseline microbiota for each patient. As such, patient’s samples were collected before chemotherapy, during induction chemotherapy (chemotherapy given to induce a remission), consolidation chemotherapy (chemotherapy given once a remission is achieved) and during maintenance therapy (chemotherapy given in lower doses to assist in prolonging a remission). All healthy sibling controls were sampled once, aligning with the time period before chemotherapy began on the patient (Additional file 1: Table S1), however, four siblings samples were collected at two time points, which were excluded from the analysis. All patients with ALL enrolled in the study received antibiotic prophylaxis with sulfamethoxazole and trimethoprim during treatment and steroid prophylaxis at the induction stage. Incidental use of antibiotics and occurrence of infections in the month before each visit were recorded. Additional file 1: Table S1 provides details of the sampling visits for each patient over the period of enrollment.
Samples were collected using the Human Microbiome Project (HMP) collection protocol section 7.3.3. with no modifications. Stool specimens were collected and transported to CHOC Children’s for deoxyribonucleic acid (DNA) extraction .
Bacterial DNA was extracted from the stool samples using the PowerSoil® DNA Isolation Kit from MO BIO Laboratories, Inc. (catalog no: 12888) and by using the protocol as described in Yooseph et al., .
Library construction and sequencing
DNA was amplified using primers that targeted the V1-V3 regions of the 16S rRNA gene . These primers included the i5 and i7 adaptor sequences for Illumina MiSeq sequencing as well as unique 8 bp indices incorporated onto both primers such that each sample received its own unique barcode pair. This method of incorporating the adaptors and index sequences onto the primers at the polymerase chain reaction (PCR) stage provided minimal loss of sequence data when compared to previous library construction methods that would ligate the adaptors to every amplicon after amplification. This method also allows generating sequence reads which were all in the same 5′-3′ orientation. Using approximately 100 ng of extracted DNA, the amplicons were generated with Platinum Taq polymerase (ThermoFisher, catalog no: 11304-011) and by using the following cycling conditions: 95 °C for 5 min for an initial denaturing step followed by 95 °C for 30 s, 55 °C for 30 s, 72 °C for 30 s for a total of 35 cycles followed by a final extension step of 72 °C for 7 min then stored at 4 °C. Once the PCR for each sample was completed, the amplicons were purified using the QIAquick PCR purification kit (QIAGEN, catalog no: 28104), quantified using Tecan fluorometric methods (Tecan Group, Männedorf, Switzerland), normalized, and then pooled in preparation for cluster generation followed by Illumina MiSeq sequencing using the dual index 2x300 bp format (Roche, Branford, CT) following the manufacturer’s protocol.
16S rRNA sequence data processing
After primer trimming, the paired-end reads were quality trimmed using the DynamicTrim program (available in the SolexaQA suite ). Subsequently, mothur (v.1.34.4)  was used to merge overlapping forward and reverse reads to generate contig sequences from the paired-end reads. Chimeric sequences were identified using the UCHIME  implementation in mothur; these sequences were removed from the downstream analysis. The resulting sequence set was clustered at different sequence identity thresholds (90, 95 and 97 %) using CD-HIT  to generate Operational Taxonomic Units (OTUs); OTU representatives were assigned taxonomy using the mothur implementation of the Ribosomal Database Project (RDP) Classifier . While the RDP classifier has a goal of generating genus-level taxonomic assignments, not all sequences could, however, be confidently assigned taxonomy to the genus level (using the bootstrap confidence threshold of 80 %); we denote these sequences by appending the tag unclassified to the end of their taxonomic assignment (one of phylum, class, order, or family levels).
Electronic medical records
Study-specific participant data was collected at each specimen collection time point. Also, multiple time-stamped data points on each participant starting from the initial visit were gathered in the electronic medical record (EMR) at CHOC Children’s. Events recorded in the EMR were obtained for each patient identified by anonymized participant identification numbers.
16S rRNA data analysis
The Shannon diversity index was calculated from the OTUs in the samples to assess the alpha diversity of the microbial communities they represent. This was done using mothur . Calculations were performed on sub-sampled sequences to account for differences in sequencing depth across the samples. To test whether the microbiota diversity difference between the patients and controls was statistically significant, we applied the Wilcoxon Rank Sum test in R, the statistical programming language,  to calculate the p-value for the comparison. The mean age difference between patient and control groups was evaluated using the t-test. For the average microbiota diversity between different time points, we used a paired (one sample) t-test.
The signature associated with the patient and control groups were identified using Random Forests (RF) , as implemented in the RF package (version 4.6–10) that is available in the R programming language and environment for statistical computing (version 3.2.2) . For this analysis, OTUs (at 97 % identity threshold) with the same taxonomic classification were combined into a single bin, thereby generating a set of taxon bins where each bin had a unique taxonomy; we note here that the best taxonomic resolution for a taxon bin is at the genus level and that some of the bins may have less resolved taxonomy (that is, at one of phylum, class, order, or family levels). The sequence counts in these taxon bins were used to calculate their relative abundance in each sample, and these abundances were used as input features for the RF analysis. Gender, antibiotic use, and alpha diversity (Shannon index) were also included as input features.
Results and discussion
The patient cohort included 28 children and adolescent participants with ALL, age 3–23 years (median age 9.9 years), and consisted of 19 male and nine female members. Similarly, the control cohort included 23 sibling controls age 2–25 years (median age 9.3 years) with 16 females and seven males (Table 1). The patient cohort showed a high prevalence of ALL in males; it has been demonstrated previously that childhood cancer occurs more frequently in males than in females . The mean age difference between patient and control groups was not statistically significant based on a t-test. A total of 180 fecal samples were collected and used over a period of one year for 16S rRNA sequencing and analysis (Additional file 1: Table S1). There were an average 16,270 annotated sequences per sample, (82,969 – maximum, and 10,498 – median). Four patients were positive for C. difficile. The low incidence of C. difficile in this immunocompromised patient population may reflect the efficiency of antibiotic prophylaxis and other mitigation strategies employed by the clinical care team.
GI microbiota composition
The microbiota profiles (Visit 1) of both patient and control groups were dominated by members of the genera Bacteroides, Prevotella, and Faecalibacterium, with these having mean abundances of 62.2, 7.3, and 6.4 % respectively, in the patient group, and 40.2, 12.2, and 8.3 % respectively, in the control group. The rank orders of the remaining lower abundant taxa were different in the two groups (Fig. 1). However, at the genus level, both patient and control groups share many taxa in common (Fig. 1).
The OTUs identified from the samples were used to compute the alpha diversity of the microbial communities (using the Shannon index). These calculations reveal the microbiota diversity (using OTUs at 97 % identity threshold) of the Patient group (Visit 1) to be lower than that of the Control group (Visit 1), with this difference being statistically significant (p-value 0.0012, Wilcoxon Rank Sum test) (Fig. 2). The Patient group was further partitioned into those who reported taking antibiotics in the one month period prior to Visit 1 (Patient_A group) and those who did not take antibiotics in this period (Patient_NA group). The microbiota diversity in each of these groups was compared with that of the Control group, and also found to be different using the Wilcoxon Rank Sum test (p-value for Control vs. Patient_A: 0.0027; p-value for Control vs. Patient_NA: 0.028). These observations also hold for diversity calculations based on OTUs generated using 95 and 90 % identity thresholds as well (Additional file 2: Figure S1).
It was possible to distinguish between the patient and control groups (Visit 1) based on their microbiota profiles. For this assessment, we used RF to identify features that could discriminate between the two groups. We observed a high classification accuracy, as measured by Area Under the Curve (AUC), of 87.9 % (Fig. 3). The input features (taxon abundances and sample metadata) were ranked by order of their importance in the classification (using their MeanDecreaseGini value). Based on this, the top taxa include Anaerostipes, Coprococcus, Roseburia, and Ruminococcus2 (all of these being in higher relative abundance in the Control group), while metadata on antibiotic use and alpha diversity are also important.
The protective role of commensal intestinal bacteria in human disease is increasingly being appreciated. However, microbiome studies that investigate the role of bacteria in human disease have focused primarily on inflammatory bowel disease, which is caused by a chronic inflammatory process. Previous studies have shown that the intestinal microbiota in patients with inflammatory bowel disease (IBD) is characterized by a contraction of Firmicutes and Bacteroidetes and an expansion of Proteobacteria . Lachnospiraceae (which comprises the Clostridium XIVa and IV groups within the order Clostridiales) and Roseburia (butyrate-producing bacterium) were greatly reduced, in acute leukemia patients compared to a healthy sibling (Fig. 1), but the Bacteroides are increased in patients compare to sibling controls. Bacteria producing butyrate play a major role in the composition of the mucus layer, as butyrate is an important energy source for intestinal epithelial cells and plays a role in the maintenance of colonic homeostasis . Several intestinal bacteria produce short chain fatty acids (SCFAs), with butyrate being the most thoroughly investigated. Butyrate is produced by F. prausnitzii and Clostridium XIVa and has been shown to have profound anti-inflammatory effects [32–36]. The observed reduction in Lachnospiraceae and Roseburia in acute leukemia patients may increase the risk of developing chemotherapy induced mucositis and other GI complications in childhood leukemia patients.
Microbiota diversity changes during chemotherapy and maintenance therapy
We assessed the changes in microbiota diversity for the patients after chemotherapy (using OTUs calculated at 97 % identity threshold). As part of this, we evaluated differences in diversity for three cases: (a) diversity at Visit 2 (after chemotherapy) versus diversity at Visit 1 (before chemotherapy), (b) diversity at the final visit versus diversity at Visit 1, and (c) the average diversity for visits after Visit 1 versus diversity at Visit 1. For each of these comparisons, we used a paired (one sample) t-test to assess whether the difference in diversity was 0 or greater than 0. As shown in Fig. 4, while the diversity increase was not significant for case (a) with a p-value of 0.318, the increase for subsequent visits were significant, with p-values of 0.00026 and 0.00643 for cases (b) and (c) respectively.
It is likely that these observed changes in the GI microbiota diversity during therapy are the result of several factors including a direct influence of some of the therapeutic compounds or combination of chemotherapy and steroid prophylaxis on the gut flora, an indirect effect of chemotherapy on the immune system, which, in turn, affects the microbiome, or dietary choices during therapy. The lymphoid progenitor cells in ALL patients are affected during the disease, leading to the impaired immune system typically observed at the time of diagnosis [21, 37]. It is known that microbiome community structure is determined by both host and environmental factors . If any one of these factors is greatly perturbed, a drastic composition shift in the composition of the microbiome can be expected.
Several microbiome studies have demonstrated that administration of antibiotics can perturb this microflora temporarily and in particular cases permanently [12–15, 39]. Antibiotic-induced shifts can increase susceptibility to C. difficile infection . Similarly, mucosal barrier injury, characterized by both inflammation and cell loss in the epithelial barrier lining of the gastrointestinal tract, is one of the most debilitating side effects of radiotherapy and chemotherapy treatment [39, 40]. However, the composition of GI microbiota in pediatric and adolescent leukemia patients and the microbiota changes after contemporary chemotherapy has not been investigated. Our study is the first, to the best of our knowledge, to address a population of pediatric and adolescent patients with acute leukemia and to compare these patients with their sibling controls. This characterization of the GI microbiota in pediatric and adolescent patients with acute leukemia has provided significant information on GI microbiota populations in immunocompromised individuals and opens up the potential for developing novel diagnostics based on stool tests as well as developing therapies to improve the dysbiotic condition of the microbiota at the time of diagnosis and in the earliest stages of chemotherapy. This creates the potential for the better management of GI and systemic complications associated with immunodeficiency and other disease conditions of this type. Furthermore, characterizing the GI microbiota dynamics and immune response following chemotherapy treatment will address what alterations happen to the GI microbiota during and following chemotherapy regimens and may be correlated with response to treatment. We anticipate expanding this study to get more detailed information on microbial profiles that are associated with or that lead to the development of various infections including C. difficile and persistent diarrhea. Further, we may find that the composition of the microbiota may ultimately be used as an indicator as to how well a patient may respond to different chemotherapy treatments.
ALL, acute lymphoblastic leukemia; GI, gastrointestinal; HMP, Human Microbiome Project; OTUs, Operational Taxonomic Units; PCR, polymerase chain reaction; RDP, Ribosomal Database Project; RF, Random Forests; rRNA, ribosomal ribonucleic acid.
Bates D, Maechler M, Bolker B, Walker S. lme4: Linear mixed-effects models using Eigen and S4. R package version 1.1-5. In: R package version 02–3. 2014.
Burke PW, Douer D. Acute lymphoblastic leukemia in adolescents and young adults. Acta Haematol. 2014;132(3–4):264–73.
Croswell A, Amir E, Teggatz P, Barman M, Salzman NH. Prolonged impact of antibiotics on intestinal microbial ecology and susceptibility to enteric Salmonella infection. Infect Immun. 2009;77(7):2741–53.
Holler E, Butzhammer P, Schmid K, Hundsrucker C, Koestler J, Peter K, Zhu W, Sporrer D, Hehlgans T, Kreutz M, et al. Metagenomic analysis of the stool microbiome in patients receiving allogeneic stem cell transplantation: loss of diversity is associated with use of systemic antibiotics and more pronounced in gastrointestinal graft-versus-host disease. Biol Blood Marrow Transplant. 2014;20(5):640–5.
Hudson MM, Mertens AC, Yasui Y, Hobbie W, Chen H, Gurney JG, Yeazel M, Recklitis CJ, Marina N, Robison LR, et al. Health status of adult long-term survivors of childhood cancer: a report from the Childhood Cancer Survivor Study. JAMA. 2003;290(12):1583–92.
Hudson MM, Ness KK, Gurney JG, Mulrooney DA, Chemaitilly W, Krull KR, Green DM, Armstrong GT, Nottage KA, Jones KE, et al. Clinical ascertainment of health outcomes among adults treated for childhood cancer. JAMA. 2013;309(22):2371–81.
Oeffinger KC, Mertens AC, Sklar CA, Kawashima T, Hudson MM, Meadows AT, Friedman DL, Marina N, Hobbie W, Kadan-Lottick NS, et al. Chronic health conditions in adult survivors of childhood cancer. N Engl J Med. 2006;355(15):1572–82.
Tai E, Buchanan N, Townsend J, Fairley T, Moore A, Richardson LC. Health status of adolescent and young adult cancer survivors. Cancer. 2012;118(19):4884–91.
Macpherson AJ, Gatto D, Sainsbury E, Harriman GR, Hengartner H, Zinkernagel RM. A primitive T cell-independent mechanism of intestinal mucosal IgA responses to commensal bacteria. Science. 2000;288(5474):2222–6.
Vijay-Kumar M, Aitken JD, Carvalho FA, Cullender TC, Mwangi S, Srinivasan S, Sitaraman SV, Knight R, Ley RE, Gewirtz AT. Metabolic syndrome and altered gut microbiota in mice lacking Toll-like receptor 5. Science. 2010;328(5975):228–31.
Guinane CM, Cotter PD. Role of the gut microbiota in health and chronic gastrointestinal disease: understanding a hidden metabolic organ. Ther Adv Gastroenterol. 2013;6(4):295–308.
Dethlefsen L, Huse S, Sogin ML, Relman DA. The pervasive effects of an antibiotic on the human gut microbiota, as revealed by deep 16S rRNA sequencing. PLoS Biol. 2008;6(11):e280.
Dethlefsen L, Relman DA. Incomplete recovery and individualized responses of the human distal gut microbiota to repeated antibiotic perturbation. Proc Natl Acad Sci U S A. 2011;108 Suppl 1:4554–61.
Hernandez E, Bargiela R, Diez MS, Friedrichs A, Perez-Cobas AE, Gosalbes MJ, Knecht H, Martinez-Martinez M, Seifert J, von Bergen M, et al. Functional consequences of microbial shifts in the human gastrointestinal tract linked to antibiotic treatment and obesity. Gut Microbes. 2013;4(4):306–15.
Perez-Cobas AE, Artacho A, Knecht H, Ferrus ML, Friedrichs A, Ott SJ, Moya A, Latorre A, Gosalbes MJ. Differential effects of antibiotic therapy on the structure and function of human gut microbiota. PLoS One. 2013;8(11):e80201.
Pallav K, Dowd SE, Villafuerte J, Yang X, Kabbani T, Hansen J, Dennis M, Leffler DA, Newburg DS, Kelly CP. Effects of polysaccharopeptide from Trametes versicolor and amoxicillin on the gut microbiome of healthy volunteers: a randomized clinical trial. Gut Microbes. 2014;5(4):458–67.
Maurice CF, Haiser HJ, Turnbaugh PJ. Xenobiotics shape the physiology and gene expression of the active human gut microbiome. Cell. 2013;152(1–2):39–50.
Theriot CM, Koenigsknecht MJ, Carlson Jr PE, Hatton GE, Nelson AM, Li B, Huffnagle GB, J ZL, Young VB. Antibiotic-induced shifts in the mouse gut microbiome and metabolome increase susceptibility to Clostridium difficile infection. Nat Commun. 2014;5:3114.
van Vliet MJ, Harmsen HJ, de Bont ES, Tissing WJ. The role of intestinal microbiota in the development and severity of chemotherapy-induced mucositis. PLoS Pathog. 2010;6(5):e1000879.
McInnes P, Cutting M. Core microbiome sampling protocol A HMP Protocol # 07–001. Manual of Procedures. Hum Microbiome Proj. 2010;11:1–109.
Kitchingman GR, Rovigatti U, Mauer AM, Melvin S, Murphy SB, Stass S. Rearrangement of immunoglobulin heavy chain genes in T cell acute lymphoblastic leukemia. Blood. 1985;65(3):725–9.
Jeraldo P, Chia N, Goldenfeld N. On the suitability of short reads of 16S rRNA for phylogeny-based analyses in environmental surveys. Environ Microbiol. 2011;13(11):3000–9.
Cox MP, Peterson DA, Biggs PJ. SolexaQA: At-a-glance quality assessment of Illumina second-generation sequencing data. BMC Bioinf. 2010;11:485.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200.
Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22(13):1658–9.
Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261–7.
R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2008. ISBN 3-900051-07-0.
Breiman L. Machine Learning. 2001;45:5. doi:10.1023/A:1010933404324.
Pearce MS, Parker L. Childhood cancer registrations in the developing world: still more boys than girls. Int J Cancer. 2001;91(3):402–6.
Frank DN, St Amand AL, Feldman RA, Boedeker EC, Harpaz N, Pace NR. Molecular-phylogenetic characterization of microbial community imbalances in human inflammatory bowel diseases. Proc Natl Acad Sci U S A. 2007;104(34):13780–5.
Hamer HM, Jonkers D, Venema K, Vanhoutvin S, Troost FJ, Brummer RJ. Review article: the role of butyrate on colonic function. Aliment Pharmacol Ther. 2008;27(2):104–19.
Barcenilla A, Pryde SE, Martin JC, Duncan SH, Stewart CS, Henderson C, Flint HJ. Phylogenetic relationships of butyrate-producing bacteria from the human gut. Appl Environ Microbiol. 2000;66(4):1654–61.
Duncan SH, Hold GL, Harmsen HJ, Stewart CS, Flint HJ. Growth requirements and fermentation products of Fusobacterium prausnitzii, and a proposal to reclassify it as Faecalibacterium prausnitzii gen. nov., comb. nov. Int J Syst Evol Microbiol. 2002;52(Pt 6):2141–6.
Hamer HM, Jonkers DM, Bast A, Vanhoutvin SA, Fischer MA, Kodde A, Troost FJ, Venema K, Brummer RJ. Butyrate modulates oxidative stress in the colonic mucosa of healthy humans. Clin Nutr. 2009;28(1):88–93.
Nancey S, Bienvenu J, Coffin B, Andre F, Descos L, Flourie B. Butyrate strongly inhibits in vitro stimulated release of cytokines in blood. Dig Dis Sci. 2002;47(4):921–8.
Zhang XL, Komada Y, Chipeta J, Li QS, Inaba H, Azuma E, Yamamoto H, Sakurai M. Intracellular cytokine profile of T cells from children with acute lymphoblastic leukemia. Cancer Immunol Immunother. 2000;49(3):165–72.
Human Microbiome Project C. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486(7402):207–14.
Vangay P, Ward T, Gerber JS, Knights D. Antibiotics, pediatric dysbiosis, and disease. Cell Host Microbe. 2015;17(5):553–64.
Sonis ST. Pathobiology of mucositis. Semin Oncol Nurs. 2004;20(1):11–5.
This project was funded by Hyundai Motor America and Hyundai Hope on Wheels and the J. Craig Venter Institute.
Availability of data and materials
The sequencing data (16S rRNA) utilized in the study is available via the NIH BioProject ID PRJNA318000.
Managed study participant’s consent and sample collection: LS, KBZ. Participated in experiments: SVR, KJM, MT, KEN. Performed computational analysis: SY, DMH, AT. Analyzed the data and provide input to the manuscript: SRV, SY, DMH, KBZ, MT, AT, SKH, RP, LS, KEN. Wrote the paper: SY, SVR, KEN, RP. Conceived, designed, and supervised: KEN, LS. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Informed consent (and assent, where appropriate) and authorization to use, create and disclose health information for research was obtained from and documented for each research participant enrolled to study in the Hyundai Cancer Institute at CHOC Children’s.
Ethics approval and consent to participate
This study was reviewed and approved by the In House Institutional Review Board at CHOC Children’s (protocol number 120555).
List of fecal samples collected for each patient over the period of enrollment for 16S rRNA gene sequencing and analysis. (XLSX 17 kb)
Box-plots of the alpha diversity of OTUs at 95 % and 90 % identity threshold of the Control and Patient groups. The Patient group is further partitioned into the group taking antibiotics 1-month period Visit 1 (Patient_A) and the group not taking antibiotics (Patient_NA). (A) Alpha diversity for 95 % OTUs, the Y-axis denotes alpha diversity (Shannon Index values). The Wilcoxon Rank Sum test p-value = 0.00253 for Control vs Patient, p-value = 0.00328 for Control vs Patient_A and p-value = 0.05969 for Control vs Patient_NA. (B) Alpha diversity for 90 % OTUs, the Y-axis denotes alpha diversity (Shannon Index values). The Wilcoxon Rank Sum test p-value = 0.00156 for Control vs Patient, p-value = 0.00328 for Control vs Patient_A and p-value = 0.03119 for Control vs Patient_NA. In both cases (A, B) the Patient group has a lower microbiota diversity (statistically significant) compared to the Control group (p-value < 0.0026). The diversities of the Patient_A and Patient_NA groups are also significantly lower (p-value < 0.05) than the Control group. (PDF 203 kb)