Continuous pre- and post-transplant exposure to a disease-associated gut microbiome promotes hyper-acute graft-versus-host disease in wild-type mice

References

• Zeiser R, Blazar BR. Acute graft-versus-host disease — biologic process, prevention, and therapy. N Engl J Med. 2017;377:2167–2179.
   Google Scholar
• Castilla-Llorente C, Martin PJ, McDonald GB, Storer BE, Appelbaum FR, Deeg HJ, Mielcarek M, Shulman H, Storb R, Nash RA. Prognostic factors and outcomes of severe gastrointestinal GVHD after allogeneic hematopoietic cell transplantation. Bone Marrow Transplant. 2014;49:966. doi:10.1038/bmt.2014.69.
   Google Scholar
• Jones JM, Wilson R, Bealmear PM. Mortality and gross pathology of secondary disease in germfree mouse radiation chimeras. Radiat Res. 1971;45(3):577–588. doi:10.2307/3573066.
   Google Scholar
• Bekkum D, Roodenburg J, Heidt PJ, Waaij D. Mitigation of secondary disease of allogeneic mouse radiation chimeras by modification of the intestinal microflora. JNCI: J National Cancer Inst. 1974;52(2):401–404. doi:10.1093/jnci/52.2.401.
   Google Scholar
• Hill GR, Crawford JM, Cooke KR, Brinson YS, Pan L, Ferrara JLM. Total body irradiation and acute graft-versus-host disease: the role of gastrointestinal damage and inflammatory cytokines. Blood. 1997;90(8):3204–3213. doi:10.1182/blood.V90.8.3204.
   Google Scholar
• Kennedy GA, Varelias A, Vuckovic S, Le Texier L, Gartlan KH, Zhang P, Thomas G, Anderson L, Boyle G, Cloonan N, et al. Addition of interleukin-6 inhibition with tocilizumab to standard graft-versus-host disease prophylaxis after allogeneic stem-cell transplantation: a phase 1/2 trial. Lancet Oncol. 2014;15(13):1451–1459. doi:10.1016/S1470-2045(14)71017-4.
   Google Scholar
• Cooke KR, Gerbitz A, Crawford JM, Teshima T, Hill GR, Tesolin A, Rossignol DP, Ferrara JL. LPS antagonism reduces graft-versus-host disease and preserves graft-versus-leukemia activity after experimental bone marrow transplantation. J Clin Invest. 2001;107(12):1581–1589. doi:10.1172/JCI12156.
   Google Scholar
• Teshima T, Ordemann R, Reddy P, Gagin S, Liu C, Cooke KR, Ferrara JL. Acute graft-versus-host disease does not require alloantigen expression on host epithelium. Nat Med. 2002;8(6):575–581. doi:10.1038/nm0602-575.
   Google Scholar
• Koyama M, Hill GR. Alloantigen presentation and graft-versus-host disease: fuel for the fire. Blood. 2016;127(24):2963–2970. doi:10.1182/blood-2016-02-697250.
   Google Scholar
• Whangbo J, Ritz J, Bhatt A. Antibiotic-mediated modification of the intestinal microbiome in allogeneic hematopoietic stem cell transplantation. Bone Marrow Transplant. 2016;52:183–190. doi:10.1038/bmt.2016.206.
   Google Scholar
• Khoruts A, Hippen KL, Lemire AM, Holtan SG, Knights D, Young JAH. Toward revision of antimicrobial therapies in hematopoietic stem cell transplantation: target the pathogens, but protect the indigenous microbiota. Transl Res. 2017;179:116–125. doi:10.1016/j.trsl.2016.07.013.
   Google Scholar
• Koyama M, Mukhopadhyay P, Schuster IS, Henden AS, Hülsdünker J, Varelias A, Vetizou M, Kuns RD, Robb RJ, Zhang P, et al. MHC class II antigen presentation by the intestinal epithelium initiates graft-versus-host disease and is influenced by the microbiota. Immunity. 2019;51(5):885–898
   Google Scholar
• Jenq RR, Ubeda C, Taur Y, Menezes CC, Khanin R, Dudakov JA, Liu C, West ML, Singer NV, Equinda MJ, et al. Regulation of intestinal inflammation by microbiota following allogeneic bone marrow transplantation. J Exp Med. 2012;209(5):903–911. doi:10.1084/jem.20112408.
   Google Scholar
• Eriguchi Y, Takashima S, Oka H, Shimoji S, Nakamura K, Uryu H, Shimoda S, Iwasaki H, Shimono N, Ayabe T, et al. Graft-versus-host disease disrupts intestinal microbial ecology by inhibiting Paneth cell production of α-defensins. Blood. 2012;120(1):223–231. doi:10.1182/blood-2011-12-401166.
   Google Scholar
• Heimesaat MM, Nogai A, Bereswill S, Plickert R, Fischer A, Loddenkemper C, Steinhoff U, Tchaptchet S, Thiel E, Freudenberg MA, et al. MyD88/TLR9 mediated immunopathology and gut microbiota dynamics in a novel murine model of intestinal graft-versus-host disease. Gut. 2010;59(8):1079–1087. doi:10.1136/gut.2009.197434.
   Google Scholar
• 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–645. doi:10.1016/j.bbmt.2014.01.030.
   Google Scholar
• Taur Y, Jenq RR, Perales M-A, Littmann ER, Morjaria S, Ling L, No D, Gobourne A, Viale A, PB D, et al. The effects of intestinal tract bacterial diversity on mortality following allogeneic hematopoietic stem cell transplantation. Blood. 2014;124(7):1174–1182. doi:10.1182/blood-2014-02-554725.
   Google Scholar
• Jenq RR, Taur Y, Devlin SM, Ponce DM, Goldberg JD, Ahr KF, Littmann ER, Ling L, Gobourne AC, Miller LC, et al. Intestinal Blautia is associated with reduced death from graft-versus-host disease. Biol Blood Marrow Transplant. 2015;21(8):1373–1383. doi:10.1016/j.bbmt.2015.04.016.
   Google Scholar
• Biagi E, Zama D, Nastasi C, Consolandi C, Fiori J, Rampelli S, Turroni S, Centanni M, Severgnini M, Peano C, et al. Gut microbiota trajectory in pediatric patients undergoing hematopoietic SCT. Bone Marrow Transplant. 2015;50(7):992–998. doi:10.1038/bmt.2015.16.
   Google Scholar
• Doki N, Suyama M, Sasajima S, Ota J, Igarashi A, Mimura I, Morita H, Fujioka Y, Sugiyama D, Nishikawa H, et al. Clinical impact of pre-transplant gut microbial diversity on outcomes of allogeneic hematopoietic stem cell transplantation. Ann Hematol. 2017;96(9):1517–1523. doi:10.1007/s00277-017-3069-8.
   Google Scholar
• Varelias A, Ormerod KL, Bunting MD, Koyama M, Gartlan KH, Kuns RD, Lachner N, Locke KR, Lim CY, Henden AS, et al. Acute graft-versus-host disease is regulated by an IL-17–sensitive microbiome. Blood. 2017;129(15):2172–2185. doi:10.1182/blood-2016-08-732628.
   Google Scholar
• Lagkouvardos I, Pukall R, Abt B, Foesel BU, Meier-Kolthoff JP, Kumar N, Bresciani A, Martínez I, Just S, Ziegler C, et al. The Mouse Intestinal Bacterial Collection (miBC) provides host-specific insight into cultured diversity and functional potential of the gut microbiota. Nat Microbiol. 2016;1:16131. doi:10.1038/nmicrobiol.2016.131.
   Google Scholar
• Ormerod KL, Wood DLA, Lachner N, Gellatly SL, Daly JN, Parsons JD, Dal’Molin CGO, Palfreyman RW, Nielsen LK, Cooper MA, et al. Genomic characterization of the uncultured Bacteroidales family S24-7 inhabiting the guts of homeothermic animals. Microbiome. 2016;4(1):1–17. doi:10.1186/s40168-016-0181-2.
   Google Scholar
• Desai MS, Seekatz AM, Koropatkin NM, Kamada N, Hickey CA, Wolter M, Pudlo NA, Kitamoto S, Terrapon N, Muller A, et al. A dietary fiber-deprived gut microbiota degrades the colonic mucus barrier and enhances pathogen susceptibility. Cell. 2016;167(5):1339–1353.e1321. doi:10.1016/j.cell.2016.10.043.
   Google Scholar
• Potempa J, Pike RN. Corruption of innate immunity by bacterial proteases. J Innate Immun. 2009;1(2):70–87. doi:10.1159/000181144.
   Google Scholar
• Thornton RF, Murphy EC, Kagawa TF, O’Toole PW, Cooney JC. The effect of environmental conditions on expression of Bacteroides fragilis and Bacteroides thetaiotaomicron C10 protease genes. BMC Microbiol. 2012;12(1):190. doi:10.1186/1471-2180-12-190.
   Google Scholar
• Capestany CA, Kuboniwa M, Jung I-Y, Park Y, Tribble GD, Lamont RJ. Role of the Porphyromonas gingivalis InlJ protein in homotypic and heterotypic biofilm development. Infect Immun. 2006;74(5):3002–3005. doi:10.1128/IAI.74.5.3002-3005.2006.
   Google Scholar
• Sakai E, Naito M, Sato K, Hotokezaka H, Kadowaki T, Kamaguchi A, Yamamoto K, Okamoto K, Nakayama K. Construction of recombinant hemagglutinin derived from the gingipain-encoding gene of Porphyromonas gingivalis, identification of its target protein on erythrocytes, and inhibition of hemagglutination by an interdomain regional peptide. J Bacteriol. 2007;189(11):3977–3986. doi:10.1128/JB.01691-06.
   Google Scholar
• Frank AC. Molecular host mimicry and manipulation in bacterial symbionts. FEMS Microbiol Lett. 2019;366:4. doi:10.1093/femsle/fnz038.
   Google Scholar
• Kaysen A, Heintz-Buschart A, Muller EEL, Narayanasamy S, Wampach L, Laczny CC, Graf N, Simon A, Franke K, Bittenbring J, et al. Integrated meta-omic analyses of the gastrointestinal tract microbiome in patients undergoing allogeneic hematopoietic stem cell transplantation. Transl Res. 2017;186:79–94.e71. doi:10.1016/j.trsl.2017.06.008.
   Google Scholar
• Simms-Waldrip TR, Sunkersett G, Coughlin LA, Savani MR, Arana C, Kim J, Kim M, Zhan X, Greenberg DE, Xie Y, et al. Antibiotic-induced depletion of anti-inflammatory Clostridia is associated with the development of graft-versus-host disease in pediatric stem cell transplantation patients. Biol Blood Marrow Transplant. 2017;23(5):820–829. doi:10.1016/j.bbmt.2017.02.004.
   Google Scholar
• Nesta B, Spraggon G, Alteri C, Gomes Moriel D, Rosini R, Veggi D, Smith S, Bertoldi I, Pastorello I, Ferlenghi I, et al. FdeC, a novel broadly conserved Escherichia coli adhesin eliciting protection against urinary tract infections. mBio. 2012;3(2):e00010–00012. doi:10.1128/mBio.00010-12.
   Google Scholar
• Schroeder BO, Birchenough GMH, Ståhlman M, Arike L, Johansson MEV, Hansson GC, Bäckhed F. Bifidobacteria or fiber protects against diet-induced microbiota-mediated colonic mucus deterioration. Cell Host Microbe. 2018;23(1):27–40.e27. doi:10.1016/j.chom.2017.11.004.
   Google Scholar
• Borton MA, Sabag-Daigle A, Wu J, Solden LM, O’Banion BS, Daly RA, Wolfe RA, Gonzalez JF, Wysocki VH, Ahmer BMM, et al. Chemical and pathogen-induced inflammation disrupt the murine intestinal microbiome. Microbiome. 2017;5(1):47. doi:10.1186/s40168-017-0264-8.
   Google Scholar
• Lagkouvardos I, Lesker TR, Hitch TCA, Gálvez EJC, Smit N, Neuhaus K, Wang J, Baines JF, Abt B, Stecher B, et al. Sequence and cultivation study of Muribaculaceae reveals novel species, host preference, and functional potential of this yet undescribed family. Microbiome. 2019;7(1):28. doi:10.1186/s40168-019-0637-2.
   Google Scholar
• Rigottier-Gois L. Dysbiosis in inflammatory bowel diseases: the oxygen hypothesis. Isme J. 2013;7:1256. doi:10.1038/ismej.2013.80.
   Google Scholar
• Kung AL, Robinson C, George D, Karamehmet E, Garvin JH, Sosna J, Pinkney K, Foca MD, Bhatia M, Satwani P, et al. Acute gastrointestinal graft-vs-host disease is associated with increased enteric bacterial bloodstream infection density in pediatric allogeneic hematopoietic cell transplant recipients. Clin Infect Dis. 2015;61(3):350–357. doi:10.1093/cid/civ285.
   Google Scholar
• Sano H, Hilinski JA, Qayed M, Applegate K, Newton JG, Watkins B, Chiang KY, Horan J. Early blood stream infection following allogeneic hematopoietic stem cell transplantation is a risk factor for acute grade III-IV GVHD in children and adolescents. Pediatr Blood Cancer. 2018;65:2. doi:10.1002/pbc.26821.
   Google Scholar
• Beelen DW, Elmaagacli A, Müller K-D, Hirche H, Schaefer UW. Influence of intestinal bacterial decontamination using metronidazole and ciprofloxacin or ciprofloxacin alone on the development of acute graft-versus-host disease after marrow transplantation in patients with hematologic malignancies: final results and long-term follow-up of an open-label prospective randomized trial. Blood. 1999;93:3267–3275.
   Google Scholar
• Beelen DW, Haralambie E, Brandt H, Linzenmeier G, Muller KD, Quabeck K, Sayer HG, Graeven U, Mahmoud HK, Schaefer UW. Evidence that sustained growth suppression of intestinal anaerobic bacteria reduces the risk of acute graft-versus-host disease after sibling marrow transplantation. Blood. 1992;80(10):2668–2676. doi:10.1182/blood.V80.10.2668.2668.
   Google Scholar
• Routy B, Letendre C, Enot D, Chénard-Poirier M, Mehraj V, Séguin NC, Guenda K, Gagnon K, Woerther P-L, Ghez D, et al. The influence of gut-decontamination prophylactic antibiotics on acute graft-versus-host disease and survival following allogeneic hematopoietic stem cell transplantation. OncoImmunology. 2017;6(1):e1258506. doi:10.1080/2162402X.2016.1258506.
   Google Scholar
• Franzosa EA, Huang K, Meadow JF, Gevers D, Lemon KP, Bohannan BJ, Huttenhower C. Identifying personal microbiomes using metagenomic codes. Proc Natl Acad Sci U S A. 2015;112(22):E2930–2938. doi:10.1073/pnas.1423854112.
   Google Scholar
• Dun CAJ, de Bont ESJM, van Vliet MJ, Kamps WA, Tissing WJE, Harmsen HJM, Meessen NEL. Chemotherapy treatment in pediatric patients with acute myeloid leukemia receiving antimicrobial prophylaxis leads to a relative increase of colonization with potentially pathogenic bacteria in the gut. Clin Infect Dis. 2009;49(2):262–270. doi:10.1086/599346.
   Google Scholar
• Montassier E, Batard E, Massart S, Gastinne T, Carton T, Caillon J, Le Fresne S, Caroff N, Hardouin JB, Moreau P, et al. 16S rRNA gene pyrosequencing reveals shift in patient faecal microbiota during high-dose chemotherapy as conditioning regimen for bone marrow transplantation. Microb Ecol. 2014;67(3):690–699. doi:10.1007/s00248-013-0355-4.
   Google Scholar
• Montassier E, Gastinne T, Vangay P, Al-Ghalith GA, Bruley Des Varannes S, Massart S, Moreau P, Potel G, de La Cochetière MF, Batard E, et al. Chemotherapy-driven dysbiosis in the intestinal microbiome. Aliment Pharmacol Ther. 2015;42(5):515–528. doi:10.1111/apt.2015.42.issue-5.
   Google Scholar
• Shono Y, Docampo MD, Peled JU, Perobelli SM, Velardi E, Tsai JJ, Slingerland AE, Smith OM, Young LF, Gupta J, et al. Increased GVHD-related mortality with broad-spectrum antibiotic use after allogeneic hematopoietic stem cell transplantation in human patients and mice. Sci Transl Med. 2016;8(339):339ra371–339ra371. doi:10.1126/scitranslmed.aaf2311.
   Google Scholar
• Taur Y, Xavier JB, Lipuma L, Ubeda C, Goldberg J, Gobourne A, Lee YJ, Dubin KA, Socci ND, Viale A, et al. Intestinal domination and the risk of bacteremia in patients undergoing allogeneic hematopoietic stem cell transplantation. Clin Infect Dis. 2012;55(7):905–914. doi:10.1093/cid/cis580.
   Google Scholar
• DeFilipp Z, Peled JU, Li S, Mahabamunuge J, Dagher Z, Slingerland AE, Del Rio C, Valles B, Kempner ME, Smith M, et al. Third-party fecal microbiota transplantation following allo-HCT reconstitutes microbiome diversity. Blood Adv. 2018;2(7):745–753. doi:10.1182/bloodadvances.2018017731.
   Google Scholar
• Kaito S, Toya T, Yoshifuji K, Kurosawa S, Inamoto K, Takeshita K, Suda W, Kakihana K, Honda K, Hattori M, et al. Fecal microbiota transplantation with frozen capsules for a patient with refractory acute gut graft-versus-host disease. Blood Adv. 2018;2(22):3097–3101. doi:10.1182/bloodadvances.2018024968.
   Google Scholar
• Taur Y, Coyte K, Schluter J, Robilotti E, Figueroa C, Gjonbalaj M, Littmann ER, Ling L, Miller L, Gyaltshen Y, et al. Reconstitution of the gut microbiota of antibiotic-treated patients by autologous fecal microbiota transplant. Sci Transl Med. 2018;10(460):eaap9489. doi:10.1126/scitranslmed.aap9489.
   Google Scholar
• Stecher B. The roles of inflammation, nutrient availability and the commensal microbiota in enteric pathogen infection. Microbiol Spectr. 2015;3:3.
   Google Scholar
• Varelias A, Bunting MD, Ormerod KL, Koyama M, Olver SD, Straube J, Kuns RD, Robb RJ, Henden AS, Cooper L, et al. Recipient mucosal-associated invariant T cells control GVHD within the colon. J Clin Invest. 2018;128(5):1919–1936. doi:10.1172/JCI91646.
   Google Scholar
• Cooke KR, Kobzik L, Martin TR, Brewer J, Delmonte J Jr., Crawford JM, Ferrara JL. An experimental model of idiopathic pneumonia syndrome after bone marrow transplantation: I. The Roles of Minor H Antigens and Endotoxin. Blood. 1996;88:3230–3239.
   Google Scholar
• Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25(14):1754–1760. doi:10.1093/bioinformatics/btp324.
   Google Scholar
• Nurk S, Meleshko D, Korobeynikov A, Pevzner PA. 2017. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 27(5):824–834.
   Google Scholar
• Kang DD, Froula J, Egan R, Wang Z. MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities. PeerJ. 2015;3:e1165. doi:10.7717/peerj.1165.
   Google Scholar
• Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25(7):1043–1055. doi:10.1101/gr.186072.114.
   Google Scholar
• Olm MR, Brown CT, Brooks B, Banfield JF. dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. Isme J. 2017;11(12):2864–2868. doi:10.1038/ismej.2017.126.
   Google Scholar
• Parks DH, Chuvochina M, Waite DW, Rinke C, Skarshewski A, Chaumeil P-A, Hugenholtz P. A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life. Nat Biotechnol. 2018;36:996. doi:10.1038/nbt.4229.
   Google Scholar
• Pedersen BS, Quinlan AR. Mosdepth: quick coverage calculation for genomes and exomes. Bioinformatics. 2017;34(5):867–868.
   Google Scholar
• Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–336. doi:10.1038/nmeth.f.303.
   Google Scholar
• Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. doi:10.1186/s13059-014-0550-8.
   Google Scholar
• Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30(14):2068–2069. doi:10.1093/bioinformatics/btu153.
   Google Scholar
• Mistry J, Finn RD, Eddy SR, Bateman A, Punta M. Challenges in homology search: HMMER3 and convergent evolution of coiled-coil regions. Nucleic Acids Res. 2013;41(12):e121. doi:10.1093/nar/gkt263.
   Google Scholar
• Yin Y, Mao X, Yang J, Chen X, Mao F, Xu Y. dbCAN: a web resource for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2012;40(WebServer issue):W445–W451. doi:10.1093/nar/gks479.
   Google Scholar
• Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, Heger A, Hetherington K, Holm L, Mistry J, et al. Pfam: the protein families database. Nucleic Acids Res. 2014;42(D1):D222–D230. doi:10.1093/nar/gkt1223.
   Google Scholar
• Haft DH, Selengut JD, White O. The TIGRFAMs database of protein families. Nucleic Acids Res. 2003;31(1):371–373. doi:10.1093/nar/gkg128.
   Google Scholar
• Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–410. doi:10.1016/S0022-2836(05)80360-2.
   Google Scholar
• Suzek BE, Huang HZ, McGarvey P, Mazumder R, Wu CH. UniRef: comprehensive and non-redundant UniProt reference clusters. Bioinformatics. 2007;23(10):1282–1288. doi:10.1093/bioinformatics/btm098.
   Google Scholar
• Nguyen L-T, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2014;32(1):268–274. doi:10.1093/molbev/msu300.
   Google Scholar
• Wang H-C, Minh BQ, Susko E, Roger AJ. Modeling site heterogeneity with posterior mean site frequency profiles accelerates accurate phylogenomic estimation. Syst Biol. 2017;67(2):216–235. doi:10.1093/sysbio/syx068.
   Google Scholar
• Lechner M, Findeiß S, Steiner L, Marz M, Stadler PF, Prohaska SJ. Proteinortho: detection of (Co-)orthologs in large-scale analysis. BMC Bioinformatics. 2011;12(1):1–9. doi:10.1186/1471-2105-12-124.
   Google Scholar
• Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14:587. doi:10.1038/nmeth.4285.
   Google Scholar
• Ludwig W, Strunk O, Westram R, Richter L, Meier H, Kumar Y, Buchner A, Lai T, Steppi S, Jobb G, et al. ARB: a software environment for sequence data. Nucleic Acids Res. 2004;32(4):1363–1371. doi:10.1093/nar/gkh293.
   Google Scholar
• Chen L, Xiong Z, Sun L, Yang J, Jin Q. VFDB 2012 update: toward the genetic diversity and molecular evolution of bacterial virulence factors. Nucleic Acids Res. 2012;40(D1):D641–D645. doi:10.1093/nar/gkr989.
   Google Scholar
• Li W, Godzik A. CD-HIT: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22:1658–1659.
   Google Scholar
• Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’Hara RB, Simpson GL, Solymos P, Stevens MHH, Wagner H 2015. vegan: community ecology package. R package version 2.3-1. http://CRAN.R-project.org/package=vegan.
   Google Scholar
• Paulson JN, Stine OC, Bravo HC, Pop M. Differential abundance analysis for microbial marker-gene surveys. Nat Methods. 2013;10(12):1200–1202. doi:10.1038/nmeth.2658.
   Google Scholar
• Rohart F, Gautier B, Singh A, Lê Cao K-A. mixOmics: an R package for ‘omics feature selection and multiple data integration. PLoS Comput Biol. 2017;13:e1005752.
   Google Scholar
• Revelle W 2019. psych: procedures for psychological, psychometric, and personality research. R package version 1. 8.12. https://CRAN.R-project.org/package=psych.
   Google Scholar
• Kolde R 2015. pheatmap: pretty heatmaps. R package version 1.0.7. http://CRAN.R-project.org/package=pheatmap.
   Google Scholar
• Wickham H, Chang W, Henry L, Pedersen TL, Takahashi K, Wilke C, Woo K 2019. ggplot2: create elegant data visualisations using the grammar of graphics. R package version 3.1.1. https://CRAN.R-project.org/package=ggplot2.
   Google Scholar
• Schloerke B, Crowley J, Cook D, Briatte F, Marbach M, Thoen E, Elberg A, Larmarange J 2018. GGally: extension to ‘ggplot2ʹ. R package version 1.4.0. https://CRAN.R-project.org/package=GGally.
   Google Scholar