Associations between infant fungal and bacterial dysbiosis and childhood atopic wheeze in a nonindustrialized setting

Marie-Claire Arrieta, Andrea Arévalo, Leah Stiemsma, Pedro Dimitriu, Martha E Chico, Sofia Loor, Maritza Vaca, Rozlyn C T Boutin, Evan Morien, Mingliang Jin, Stuart E Turvey, Jens Walter, Laura Wegener Parfrey, Philip J Cooper, Brett Finlay, Marie-Claire Arrieta, Andrea Arévalo, Leah Stiemsma, Pedro Dimitriu, Martha E Chico, Sofia Loor, Maritza Vaca, Rozlyn C T Boutin, Evan Morien, Mingliang Jin, Stuart E Turvey, Jens Walter, Laura Wegener Parfrey, Philip J Cooper, Brett Finlay

Abstract

Background: Asthma is the most prevalent chronic disease of childhood. Recently, we identified a critical window early in the life of both mice and Canadian infants during which gut microbial changes (dysbiosis) affect asthma development. Given geographic differences in human gut microbiota worldwide, we studied the effects of gut microbial dysbiosis on atopic wheeze in a population living in a distinct developing world environment.

Objective: We sought to determine whether microbial alterations in early infancy are associated with the development of atopic wheeze in a nonindustrialized setting.

Methods: We conducted a case-control study nested within a birth cohort from rural Ecuador in which we identified 27 children with atopic wheeze and 70 healthy control subjects at 5 years of age. We analyzed bacterial and eukaryotic gut microbiota in stool samples collected at 3 months of age using 16S and 18S sequencing. Bacterial metagenomes were predicted from 16S rRNA data by using Phylogenetic Investigation of Communities by Reconstruction of Unobserved States and categorized by function with Kyoto Encyclopedia of Genes and Genomes ontology. Concentrations of fecal short-chain fatty acids were determined by using gas chromatography.

Results: As previously observed in Canadian infants, microbial dysbiosis at 3 months of age was associated with later development of atopic wheeze. However, the dysbiosis in Ecuadorian babies involved different bacterial taxa, was more pronounced, and also involved several fungal taxa. Predicted metagenomic analysis emphasized significant dysbiosis-associated differences in genes involved in carbohydrate and taurine metabolism. Levels of the fecal short-chain fatty acids acetate and caproate were reduced and increased, respectively, in the 3-month stool samples of children who went on to have atopic wheeze.

Conclusions: Our findings support the importance of fungal and bacterial microbiota during the first 100 days of life on the development of atopic wheeze and provide additional support for considering modulation of the gut microbiome as a primary asthma prevention strategy.

Keywords: Asthma; atopy; gut microbiome; mycobiome; nonindustrialized setting; short-chain fatty acids; wheeze.

Copyright © 2017 The Authors. Published by Elsevier Inc. All rights reserved.

Figures

Graphical abstract
Graphical abstract
Fig 1
Fig 1
Bacterial dysbiosis associated with development of atopic wheeze. A, Relative abundance of bacterial genera within the top 100 OTUs among the 2 phenotypes at 3 months. Rectangle colors correspond to the bacterial genera in the legend. Rectangles represent merged OTUs belonging to the same genus. NA, OTUs without genus-level taxonomic assignment. B, Log fold change of OTUs that are significantly (false discovery rate < 0.05) abundant between AWs and control subjects calculated by using DESeq2. Error bars denote SDs. Taxonomic identification and OTU numbers are denoted on the y-axis.
Fig 2
Fig 2
Fungal dysbiosis associated with development of atopic wheeze. A, Relative abundance of fungal genera within the top 100 OTUs among the 2 phenotypes at 3 months. Rectangle colors correspond to fungal genera in the legend. Rectangles represent merged OTUs belonging to the same genus or higher taxa if genus information was not obtained in the analysis. NA, OTUs without genus-level taxonomic assignment. B, Log fold change of OTUs significantly (false discovery rate < 0.05) abundant between AWs and control subjects calculated by using DESeq2. Error bars denote SDs. Taxonomic identification and OTU numbers are denoted on the y-axis.
Fig 3
Fig 3
Fungus-specific qPCR and Illumina sequence data. A, qPCR quantification (standard curve method) of P kudriavzevii in feces of 3-month-old infants from the ECUAVIDA study (control subjects, 68; AWs, 27). B, Relative abundance of fungal reads over the total number of eukaryotic reads from 18S data set. Only samples with more than 1000 sequences after filtering were included in the 18S high-throughput sequencing analysis (control subjects, 23; AWs, 10). C, qPCR quantification (standard curve method) of fungus-specific 18S copies in feces of 3-month-old infants from the ECUAVIDA study (control subjects, 68; AWs, 27). In Fig 3, A and C, graphs are based on a logarithmic scale. All graphs show Mann-Whitney tests: *P < .05 and ***P < .001.
Fig 4
Fig 4
Concentration of the 6 most abundant SCFAs in feces of infants with atopic wheeze and control infants at 3 months of age measured by using gas chromatography and normalized to feces wet weight. Acetate (A), butyrate (B), propionate (C), isobutyrate (D), valerate (E), and caproate (F). (AWs, 27; control subjects, 70; *P < .05 and **P < .01, Mann-Whitney test).
Fig E1
Fig E1
Standard curves for fungus-specific 18S (A) and P kudriavzevii(B) qPCR reactions. Standard curves were generated by using 1:5 dilutions of a 1 ng/μL stock of 18S amplicons obtained from PCR reactions done with fungus-specific 18S primers (FR1 and FF390) and purified P kudriavzevii template DNA (see the Methods section). DNA concentrations are denoted in logarithmic scale (y-axis) in relation to PCR cycle number (x-axis). Exponential trend lines and equations, as well as R2 values, are included in each graph.
Fig E2
Fig E2
Relative abundance of the 100 most abundant bacterial OTUs in feces of Canadian versus Ecuadorian 3-month-old infants. Rectangle colors correspond to the bacterial genera in the legend. Rectangles represent merged OTUs belonging to the same genus. NA, OTUs without genus-level taxonomic assignment.
Fig E3
Fig E3
Both α (A and C; Chao1 index) and β (B and D; principal coordinates analysis [Bray-Curtis]) bacterial (Fig E3, A and B) and fungal (Fig E3, C and D) diversity in fecal samples from 3-month-old infants who presented with atopic wheeze at 5 years (AW) versus healthy control subjects (CTRL). Sequencing results from all 16S samples and from 33 18S samples (control subjects, 23; AWs, 10) were retained after applying a cutoff of 1000 sequences per sample. The atopic wheeze phenotype did not explain significant changes in these measurements (P = .323 and .676 [Mann-Whitney test and permutational multivariate ANOVA] for Fig E3, A and C, respectively).
Fig E4
Fig E4
Heat map of biweight correlations between the top 100 bacterial (x-axis) and top 100 fungal (y-axis) OTUs in 33 fecal samples (control subjects, 23; AWs, 10) from 3-month-old infants recruited in the ECUAVIDA cohort study. Significant correlation values are denoted with a cross (P < .05, false discovery rate).
Fig E5
Fig E5
Phylogenetic placement of P kudriavzevii OTUs within Pichia species sequences from the SILVA 123 reference database. 18S ribosomal DNA sequences labeled as Pichia species or Pichiaceae in the 99% SILVA reference database were aligned and then used to construct a reference tree with RAxML. OTUs from this study assigned to Pichia species were placed in this tree with the RAxML evolutionary placement algorithm. All OTUs cluster within a clade of P kudriavzevii.
Fig E6
Fig E6
PICRUSt-predicted Kyoto Encyclopedia of Genes and Genomes functional categories with significant differences in relative abundance between control subjects and AWs (Welch t test, q values are shown in figure).

References

    1. Asher I., Pearce N. Global burden of asthma among children. Int J Tuberc Lung Dis. 2014;18:1269–1278.
    1. Martinez F.D., Vercelli D. Asthma. Lancet. 2013;382:1360–1372.
    1. Arrieta M.C., Stiemsma L.T., Amenyogbe N., Brown E.M., Finlay B. The intestinal microbiome in early life: health and disease. Front Immunol. 2014;5:427.
    1. Russell S.L., Gold M.J., Hartmann M., Willing B.P., Thorson L., Wlodarska M. Early life antibiotic-driven changes in microbiota enhance susceptibility to allergic asthma. EMBO Rep. 2012;13:440–447.
    1. Russell S.L., Gold M.J., Reynolds L.A., Willing B.P., Dimitriu P., Thorson L. Perinatal antibiotic-induced shifts in gut microbiota have differential effects on inflammatory lung diseases. J Allergy Clin Immunol. 2015;135:100–109.
    1. Russell S.L., Gold M.J., Willing B.P., Thorson L., McNagny K.M., Finlay B.B. Perinatal antibiotic treatment affects murine microbiota, immune responses and allergic asthma. Gut Microbes. 2013;4:158–164.
    1. Arrieta M.C., Stiemsma L.T., Dimitriu P.A., Thorson L., Russell S., Yurist-Doutsch S. Early infancy microbial and metabolic alterations affect risk of childhood asthma. Sci Transl Med. 2015;7:307ra152.
    1. Stiemsma L., Arrieta M.C., Dimitriu P., Cheng J., Thorson L., Lefebvre D. Shifts in Lachnospira and Clostridium sp. in the 3-month stool microbiome are associated with preschool-age asthma. Clin Sci (Lond) 2016;130:2199–2207.
    1. Stiemsma L.T., Turvey S.E. Asthma and the microbiome: defining the critical window in early life. Allergy Asthma Clin Immunol. 2017;13:3.
    1. Takaro T.K., Scott J.A., Allen R.W., Anand S.S., Becker A.B., Befus A.D. The Canadian Healthy Infant Longitudinal Development (CHILD) birth cohort study: assessment of environmental exposures. J Expo Sci Environ Epidemiol. 2015;25:580–592.
    1. Fujimura K.E., Sitarik A.R., Havstad S., Lin D.L., Levan S., Fadrosh D. Neonatal gut microbiota associates with childhood multisensitized atopy and T cell differentiation. Nat Med. 2016;22:1187–1191.
    1. The International Study of Asthma and Allergies in Childhood (ISAAC) Steering Committee. Worldwide variation in prevalence of symptoms of asthma, allergic rhinoconjunctivitis, and atopic eczema: ISAAC. Lancet. 1998;351:1225–1232.
    1. Mallol J., Crane J., von Mutius E., Odhiambo J., Keil U., Stewart A. The International Study of Asthma and Allergies in Childhood (ISAAC) Phase Three: a global synthesis. Allergol Immunopathol (Madr) 2013;41:73–85.
    1. Cooper P.J., Rodrigues L.C., Cruz A.A., Barreto M.L. Asthma in Latin America: a public heath challenge and research opportunity. Allergy. 2009;64:5–17.
    1. Yatsunenko T., Rey F.E., Manary M.J., Trehan I., Dominguez-Bello M.G., Contreras M. Human gut microbiome viewed across age and geography. Nature. 2012;486:222–227.
    1. Escobar J.S., Klotz B., Valdes B.E., Agudelo G.M. The gut microbiota of Colombians differs from that of Americans, Europeans and Asians. BMC Microbiol. 2014;14:311.
    1. Forno E., Gogna M., Cepeda A., Yanez A., Sole D., Cooper P. Asthma in Latin America. Thorax. 2015;70:898–905.
    1. Cooper P.J., Chico M.E., Platts-Mills T.A., Rodrigues L.C., Strachan D.P., Barreto M.L. Cohort profile: the Ecuador Life (ECUAVIDA) study in Esmeraldas Province, Ecuador. Int J Epidemiol. 2015;44:1517–1527.
    1. World Health Organization Parasitic Disease Programme . World Health Organization; Geneva: 1985. Diagnostic techniques for intestinal parasitic infections (IPI) applicable to primary health care (PHC) services.
    1. Kozich J.J., Westcott S.L., Baxter N.T., Highlander S.K., Schloss P.D. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013;79:5112–5120.
    1. Comeau A.M., Douglas G.M., Langille M.G. Microbiome Helper: a custom and streamlined workflow for microbiome research. mSystems. 2017;2
    1. Schloss P.D., Westcott S.L., Ryabin T., Hall J.R., Hartmann M., Hollister E.B. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–7541.
    1. Kuczynski J., Stombaugh J., Walters W.A., Gonzalez A., Caporaso J.G., Knight R. Using QIIME to analyze 16S rRNA gene sequences from microbial communities. Curr Protoc Bioinformatics. 2011 Chapter 10:Unit 10.7.
    1. Eren A.M., Morrison H.G., Lescault P.J., Reveillaud J., Vineis J.H., Sogin M.L. Minimum entropy decomposition: unsupervised oligotyping for sensitive partitioning of high-throughput marker gene sequences. ISME J. 2015;9:968–979.
    1. Eren A.M., Maignien L., Sul W.J., Murphy L.G., Grim S.L., Morrison H.G. Oligotyping: differentiating between closely related microbial taxa using 16S rRNA gene data. Methods Ecol Evol. 2013;4:1111–1119.
    1. DeSantis T.Z., Hugenholtz P., Larsen N., Rojas M., Brodie E.L., Keller K. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–5072.
    1. Yilmaz P., Parfrey L.W., Yarza P., Gerken J., Pruesse E., Quast C. The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. Nucleic Acids Res. 2014;42(Database issue):D643–D648.
    1. Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–1313.
    1. Caporaso J.G., Bittinger K., Bushman F.D., DeSantis T.Z., Andersen G.L., Knight R. PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics. 2010;26:266–267.
    1. Berger S.A., Krompass D., Stamatakis A. Performance, accuracy, and Web server for evolutionary placement of short sequence reads under maximum likelihood. Syst Biol. 2011;60:291–302.
    1. Chemidlin Prevost-Boure N., Christen R., Dequiedt S., Mougel C., Lelievre M., Jolivet C. Validation and application of a PCR primer set to quantify fungal communities in the soil environment by real-time quantitative PCR. PLoS One. 2011;6:e24166.
    1. Langille M.G., Zaneveld J., Caporaso J.G., McDonald D., Knights D., Reyes J.A. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31:814–821.
    1. Parks D.H., Tyson G.W., Hugenholtz P., Beiko R.G. STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics. 2014;30:3123–3124.
    1. Love M.I., Huber W., Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
    1. Campbell J.M., Fahey G.C., Jr., Wolf B.W. Selected indigestible oligosaccharides affect large bowel mass, cecal and fecal short-chain fatty acids, pH and microflora in rats. J Nutr. 1997;127:130–136.
    1. McMurdie P.J., Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8:e61217.
    1. Neuwirth E. 2011. RColorBrewer: ColorBrewer palettes. Available at: .
    1. Oksanen J., Blanchet F.G., Kindt R., Legendre P., Minchin P.R., O'Hara R.B. 2013. vegan: Community Ecology Package. Available at: .
    1. Ploner A. 2012. Heatplus: Heatmaps with row and/or column covariates and colored clusters. Available at: .
    1. Pollard K., Gilbert H.N., Ge Y., Taylor S., Dudoit S. Springer; 2005. Multiple testing procedures: R multtst Package and Applications to Genomics, in Bioinformatics and Computational Biology Solutions Using R and Bioconductor.
    1. Wickham H. 2009. ggplot2: elegant graphics for data analysis. Available at:
    1. Wickham H. The split-apply-combine strategy for data analysis. J Stat Softw. 2011;40:1–29.
    1. Morgan X.C., Tickle T.L., Sokol H., Gevers D., Devaney K.L., Ward D.V. Dysfunction of the intestinal microbiome in inflammatory bowel disease and treatment. Genome Biol. 2012;13:R79.
    1. Martinez I., Stegen J.C., Maldonado-Gomez M.X., Eren A.M., Siba P.M., Greenhill A.R. The gut microbiota of rural Papua New Guineans: composition, diversity patterns, and ecological processes. Cell Rep. 2015;11:527–538.
    1. Cooper P.J., Rodrigues L.C., Barreto M.L. Influence of poverty and infection on asthma in Latin America. Curr Opin Allergy Clin Immunol. 2012;12:171–178.
    1. Trompette A., Gollwitzer E.S., Yadava K., Sichelstiel A.K., Sprenger N., Ngom-Bru C. Gut microbiota metabolism of dietary fiber influences allergic airway disease and hematopoiesis. Nat Med. 2014;20:159–166.
    1. Fallani M., Young D., Scott J., Norin E., Amarri S., Adam R. Intestinal microbiota of 6-week-old infants across Europe: geographic influence beyond delivery mode, breast-feeding, and antibiotics. J Pediatr Gastroenterol Nutr. 2010;51:77–84.
    1. Arpaia N., Campbell C., Fan X., Dikiy S., van der Veeken J., deRoos P. Metabolites produced by commensal bacteria promote peripheral regulatory T-cell generation. Nature. 2013;504:451–455.
    1. Underhill D.M., Iliev I.D. The mycobiota: interactions between commensal fungi and the host immune system. Nat Rev Immunol. 2014;14:405–416.
    1. Wheeler M.L., Limon J.J., Bar A.S., Leal C.A., Gargus M., Tang J. Immunological consequences of intestinal fungal dysbiosis. Cell Host Microbe. 2016;19:865–873.
    1. Noverr M.C., Falkowski N.R., McDonald R.A., McKenzie A.N., Huffnagle G.B. Development of allergic airway disease in mice following antibiotic therapy and fungal microbiota increase: role of host genetics, antigen, and interleukin-13. Infect Immun. 2005;73:30–38.
    1. Noverr M.C., Noggle R.M., Toews G.B., Huffnagle G.B. Role of antibiotics and fungal microbiota in driving pulmonary allergic responses. Infect Immun. 2004;72:4996–5003.
    1. Erb-Downward J.R., Noverr M.C. Characterization of prostaglandin E2 production by Candida albicans. Infect Immun. 2007;75:3498–3505.
    1. Kim Y.G., Udayanga K.G., Totsuka N., Weinberg J.B., Nunez G., Shibuya A. Gut dysbiosis promotes M2 macrophage polarization and allergic airway inflammation via fungi-induced PGE(2) Cell Host Microbe. 2014;15:95–102.
    1. Noverr M.C., Phare S.M., Toews G.B., Coffey M.J., Huffnagle G.B. Pathogenic yeasts Cryptococcus neoformans and Candida albicans produce immunomodulatory prostaglandins. Infect Immun. 2001;69:2957–2963.
    1. Kurtzman C.P., Robnett C.J., Basehoar-Powers E. Phylogenetic relationships among species of Pichia, Issatchenkia and Williopsis determined from multigene sequence analysis, and the proposal of Barnettozyma gen. nov., Lindnera gen. nov. and Wickerhamomyces gen. nov. FEMS Yeast Res. 2008;8:939–954.
    1. Madeeha I.R., Ikram A., Imran M. A preliminary insight of correlation between human fecal microbial diversity and blood lipid profile. Int J Food Sci Nutr. 2016;67:865–871.

Source: PubMed

赞助者和合作者

医疗条件

药物干预

3
订阅