Phylogeny-corrected identification of microbial gene families relevant to human gut colonization

PLoS Computational Biology

Volumn 14, Issue 8, 2018, Pages

  Bradley, Patrick H ^a   Nayfach, Stephen ^a,b   Pollard, Katherine S ^a,c,d  

^a GLADSTONE INSTITUTES   (United States)

^b LAWRENCE BERKELEY NATIONAL LABORATORY   (United States)

^c UNIVERSITY OF CALIFORNIA   (United States)

^d CHAN ZUCKERBERG BIOHUB   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

BACTERIA; HEALTH; PHYLOGENETIC; PROBIOTICS;

COLONISATION; DISEASE STATE; GENE FAMILIES; HEALTHY HUMANS; HUMAN GUTS; METAGENOMES; MICROBIOME; PHYLOGENETICS; PROBIOTICS; TREE OF LIFE;

GENES;

ARTICLE; BACTERIUM CONJUGATION; BACTEROIDETES; CROHN DISEASE; FALSE POSITIVE RESULT; FIMBRIA; FIRMICUTES; INTESTINE FLORA; MARKER GENE; METAGENOME; METAGENOMICS; MICROBIAL COLONIZATION; MICROBIAL GENE; NONHUMAN; PHYLOGENY; POPULATION ABUNDANCE; PROTEOBACTERIA; STATISTICAL MODEL; BACTERIUM; GENE EXPRESSION REGULATION; GENETICS; HUMAN; MICROFLORA; PHYSIOLOGY; PROCEDURES;

BACTERIA; GASTROINTESTINAL MICROBIOME; GENE EXPRESSION REGULATION, BACTERIAL; GENES, MICROBIAL; HUMANS; METAGENOME; METAGENOMICS; MICROBIOTA; PHYLOGENY;

EID: 85053085817     PISSN: 1553734X     EISSN: 15537358     Source Type: Journal    
DOI: 10.1371/journal.pcbi.1006242     Document Type: Article

References (107)

1. Slack E, Hapfelmeier S, Stecher B, Velykoredko Y, Stoel M, Lawson MAE, et al. Innate and adaptive immunity cooperate flexibly to maintain host-microbiota mutualism. Science. 2009;325(5940):617–620. doi: 10.1126/science.1172747 pmid: 19644121
2. Atarashi K, Tanoue T, Shima T, Imaoka A, Kuwahara T, Momose Y, et al. Induction of colonic regulatory T cells by indigenous Clostridium species. Science. 2011;331(6015):337–341. doi: 10.1126/science.1198469 pmid: 21205640
3. Mazmanian SK, Round JL, Kasper DL, A microbial symbiosis factor prevents intestinal inflammatory disease. Nature. 2008;453(7195):620–625. doi: 10.1038/nature07008 pmid: 18509436
4. Sassone-Corsi M, Raffatellu M, No vacancy: how beneficial microbes cooperate with immunity to provide colonization resistance to pathogens. Journal of Immunology. 2015;194(9):4081–7. doi: 10.4049/jimmunol.1403169
5. Yano JM, Yu K, Donaldson GP, Shastri GG, Ann P, Ma L, et al. Indigenous bacteria from the gut microbiota regulate host serotonin biosynthesis. Cell. 2015;161(2):264–76. doi: 10.1016/j.cell.2015.02.047 pmid: 25860609
6. Peng L, Li ZR, Green RS, Holzman IR, Lin J, Butyrate enhances the intestinal barrier by facilitating tight junction assembly via activation of AMP-activated protein kinase in Caco-2 cell monolayers. Journal of Nutrition. 2009;139(9):1619–1625. doi: 10.3945/jn.109.104638 pmid: 19625695
7. Reber SO, Siebler PH, Donner NC, Morton JT, Smith DG, Kopelman JM, et al. Immunization with a heat-killed preparation of the environmental bacterium Mycobacterium vaccae promotes stress resilience in mice. Proceedings of the National Academy of Sciences. 2016;113(22):E3130–E3139. doi: 10.1073/pnas.1600324113
8. Garrett WS, Gallini CA, Yatsunenko T, Michaud M, DuBois A, Delaney ML, et al. Enterobacteriaceae act in concert with the gut microbiota to induce spontaneous and maternally transmitted colitis. Cell Host & Microbe. 2010;8(3):292–300. doi: 10.1016/j.chom.2010.08.004
9. Kullberg MC, Ward JM, Gorelick PL, Caspar P, Hieny S, Cheever A, et al. Helicobacter hepaticus triggers colitis in specific-pathogen-free interleukin-10 (IL-10)-deficient mice through an IL-12- and gamma interferon-dependent mechanism. Infection and Immunity. 1998;66(11):5157–66. pmid: 9784517
10. Kostic A, Chun E, Robertson L, Glickman J, Gallini C, Michaud M, et al. Fusobacterium nucleatum potentiates intestinal tumorigenesis and modulates the tumor-immune microenvironment. Cell Host & Microbe. 2013;14(2):207–215. doi: 10.1016/j.chom.2013.07.007
11. Bartlett JG, Clostridium difficile-associated enteric disease. Current Infectious Disease Reports. 2002;4(6):477–483. doi: 10.1007/s11908-002-0032-0 pmid: 12433321
12. van Nood E, Vrieze A, Nieuwdorp M, Fuentes S, Zoetendal EG, de Vos WM, et al. Duodenal infusion of donor feces for recurrent Clostridium difficile. New England Journal of Medicine. 2013;368(5):407–415. doi: 10.1056/NEJMoa1205037 pmid: 23323867
13. Weingarden A, González A, Vázquez-Baeza Y, Weiss S, Humphry G, Berg-Lyons D, et al. Dynamic changes in short- and long-term bacterial composition following fecal microbiota transplantation for recurrent Clostridium difficile infection. Microbiome. 2015;3(1):10. doi: 10.1186/s40168-015-0070-0 pmid: 25825673
14. Khanna S, Vazquez-Baeza Y, González A, Weiss S, Schmidt B, Muñiz-Pedrogo DA, et al. Changes in microbial ecology after fecal microbiota transplantation for recurrent C. difficile infection affected by underlying inflammatory bowel disease. Microbiome. 2017;5(1):55. doi: 10.1186/s40168-017-0269-3 pmid: 28506317
15. Carvalho F, Koren O, Goodrich J, Johansson MV, Nalbantoglu I, Aitken J, et al. Transient inability to manage Proteobacteria promotes chronic gut inflammation in TLR5-deficient mice. Cell Host & Microbe. 2012;12(2):139–152. doi: 10.1016/j.chom.2012.07.004
16. Chassaing B, Koren O, Carvalho FA, Ley RE, Gewirtz AT, AIEC pathobiont instigates chronic colitis in susceptible hosts by altering microbiota composition. Gut. 2014;63(7):1069–1080. doi: 10.1136/gutjnl-2013-304909 pmid: 23896971
17. Vilhjálmsson BJ, Nordborg M, The nature of confounding in genome-wide association studies. Nature Reviews Genetics. 2012;14(1):1–2. doi: 10.1038/nrg3382 pmid: 23165185
18. Kim PJ, Price ND, Genetic co-occurrence network across sequenced microbes. PLoS Computational Biology. 2011;7(12):e1002340. doi: 10.1371/journal.pcbi.1002340 pmid: 22219725
19. Porter SS, Chang PL, Conow CA, Dunham JP, Friesen ML, Association mapping reveals novel serpentine adaptation gene clusters in a population of symbiotic Mesorhizobium. The ISME Journal. 2017;11(1):248–262. doi: 10.1038/ismej.2016.88 pmid: 27420027
20. Collins C, Didelot X, A phylogenetic method to perform genome-wide association studies in microbes that accounts for population structure and recombination. PLoS Computational Biology. 2018;14(2):e1005958. doi: 10.1371/journal.pcbi.1005958 pmid: 29401456
21. Zhao N, Chen J, Carroll I, Ringel-Kulka T, Epstein M, Zhou H, et al. Testing in microbiome-profiling studies with MiRKAT, the microbiome regression-based kernel association test. The American Journal of Human Genetics. 2015;96(5):797–807. doi: 10.1016/j.ajhg.2015.04.003 pmid: 25957468
22. Silverman JD, Washburne AD, Mukherjee S, David LA, A phylogenetic transform enhances analysis of compositional microbiota data. eLife. 2017;6. doi: 10.7554/eLife.21887
23. Tang ZZ, Chen G, Alekseyenko AV, Li H, A general framework for association analysis of microbial communities on a taxonomic tree. Bioinformatics. 2016;33(9):btw804. doi: 10.1093/bioinformatics/btw804
24. Lozupone CA, Hamady M, Cantarel BL, Coutinho PM, Henrissat B, Gordon JI, et al. The convergence of carbohydrate active gene repertoires in human gut microbes. Proceedings of the National Academy of Sciences. 2008;105(39):15076–15081. doi: 10.1073/pnas.0807339105
25. Felsenstein J, Phylogenies and the comparative method. The American Naturalist. 1985;. doi: 10.1086/284325
26. Grafen A, The phylogenetic regression. Philosophical Transactions of the Royal Society of London Series B, Biological Sciences. 1989;326(1233):119–57. doi: 10.1098/rstb.1989.0106 pmid: 2575770
27. Dunn CW, Zapata F, Munro C, Siebert S, Hejnol A, Pairwise comparisons across species are problematic when analyzing functional genomic data. Proceedings of the National Academy of Sciences. 2018;115(3):E409–E417. doi: 10.1073/pnas.1707515115
28. Levy A, Salas Gonzalez I, Mittelviefhaus M, Clingenpeel S, Herrera Paredes S, Miao J, et al. Genomic features of bacterial adaptation to plants. Nature Genetics. 2018;50(1):138–150. doi: 10.1038/s41588-017-0012-9 pmid: 29255260
29. Ord TJ, Martins EP, Tracing the origins of signal diversity in anole lizards: phylogenetic approaches to inferring the evolution of complex behaviour. Animal Behaviour. 2006;71(6):1411–1429. doi: 10.1016/j.anbehav.2005.12.003
30. Zaneveld JRR, Parfrey LW, Van Treuren W, Lozupone C, Clemente JC, Knights D, et al. Combined phylogenetic and genomic approaches for the high-throughput study of microbial habitat adaptation. Trends in Microbiology. 2011;19(10):472–82. doi: 10.1016/j.tim.2011.07.006 pmid: 21872475
31. Washburne AD, Morton JT, Sanders J, McDonald D, Zhu Q, Oliverio AM, et al. Methods for phylogenetic analysis of microbiome data. Nature Microbiology. 2018;3(6):652–661. doi: 10.1038/s41564-018-0156-0 pmid: 29795540
32. Nayfach S, Rodriguez-Mueller B, Garud N, Pollard KS, An integrated metagenomics pipeline for strain profiling reveals novel patterns of bacterial transmission and biogeography. Genome Research. 2016;26(11):1612–1625. doi: 10.1101/gr.201863.115 pmid: 27803195
33. Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457(7228):480–4. doi: 10.1038/nature07540 pmid: 19043404
34. Ives AR, Garland T, Phylogenetic logistic regression for binary dependent variables. Systematic Biology. 2010;59(1):9–26. doi: 10.1093/sysbio/syp074 pmid: 20525617
35. Human Microbiome Project ConsortiumStructure, function and diversity of the healthy human microbiome. Nature. 2012;486(7402):207–14. doi: 10.1038/nature11234 pmid: 22699609
36. Gevers D, Kugathasan S, Denson L, 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(3):382–392. doi: 10.1016/j.chom.2014.02.005
37. Meyer F, Overbeek R, Rodriguez A, FIGfams: yet another set of protein families. Nucleic Acids Research. 2009;37(20):6643–54. doi: 10.1093/nar/gkp698 pmid: 19762480
38. Sakamoto M, Ohkuma M, Bacteroides reticulotermitis sp. nov., isolated from the gut of a subterranean termite (Reticulitermes speratus). International Journal of Systematic and Evolutionary Microbiology. 2013;63(Pt 2):691–695. doi: 10.1099/ijs.0.040931-0 pmid: 22544795
39. Browne HP, Forster SC, Anonye BO, Kumar N, Neville BA, Stares MD, et al. Culturing of “unculturable” human microbiota reveals novel taxa and extensive sporulation. Nature. 2016;533(7604):543–546. doi: 10.1038/nature17645 pmid: 27144353
40. Swick MC, Koehler TM, Driks A, Surviving between hosts: sporulation and transmission. Microbiology Spectrum. 2016;4(4). doi: 10.1128/microbiolspec.VMBF-0029-2015 pmid: 27726794
41. De Biase D, Pennacchietti E, Glutamate decarboxylase-dependent acid resistance in orally acquired bacteria: function, distribution and biomedical implications of the gadBC operon. Molecular Microbiology. 2012;86(4):770–786. doi: 10.1111/mmi.12020 pmid: 22995042
42. Srinivasa Rao PS, Lim TM, Leung KY, Functional genomics approach to the identification of virulence genes involved in Edwardsiella tarda pathogenesis. Infection and Immunity. 2003;71(3):1343–51. doi: 10.1128/IAI.71.3.1343-1351.2003 pmid: 12595451
43. Cotter PD, Gahan CG, Hill C, A glutamate decarboxylase system protects Listeria monocytogenes in gastric fluid. Molecular Microbiology. 2001;40(2):465–75. doi: 10.1046/j.1365-2958.2001.02398.x pmid: 11309128
44. Wargo MJ, Meadows JA, Carnitine in bacterial physiology and metabolism. Microbiology. 2015;161(6):1161–1174. doi: 10.1099/mic.0.000080 pmid: 25787873
45. Staley C, Weingarden AR, Khoruts A, Sadowsky MJ, Interaction of gut microbiota with bile acid metabolism and its influence on disease states. Applied Microbiology and Biotechnology. 2017;101(1):47–64. doi: 10.1007/s00253-016-8006-6 pmid: 27888332
46. Marques JC, Oh IK, Ly DC, Lamosa P, Ventura MR, Miller ST, et al. LsrF, a coenzyme A-dependent thiolase, catalyzes the terminal step in processing the quorum sensing signal autoinducer-2. Proceedings of the National Academy of Sciences. 2014;111(39):14235–14240. doi: 10.1073/pnas.1408691111
47. Thompson J, Oliveira R, Djukovic A, Ubeda C, Xavier K, Manipulation of the quorum sensing signal ai-2 affects the antibiotic-treated gut microbiota. Cell Reports. 2015;10(11):1861–1871. doi: 10.1016/j.celrep.2015.02.049 pmid: 25801025
48. Wu M, McNulty NP, Rodionov DA, Khoroshkin MS, Griffin NW, Cheng J, et al. Genetic determinants of in vivo fitness and diet responsiveness in multiple human gut Bacteroides. Science. 2015;350(6256):aac5992–aac5992. doi: 10.1126/science.aac5992 pmid: 26430127
49. Wang J, Jia H, Metagenome-wide association studies: fine-mining the microbiome. Nature Reviews Microbiology. 2016;14(8):508–522. doi: 10.1038/nrmicro.2016.83 pmid: 27396567
50. Shin NR, Whon TW, Bae JW, Proteobacteria: microbial signature of dysbiosis in gut microbiota. Trends in Biotechnology. 2015;33(9):496–503. doi: 10.1016/j.tibtech.2015.06.011 pmid: 26210164
51. Lynch SV, Pedersen O, The human intestinal microbiome in health and disease. The New England Journal of Medicine. 2016;375(24):2369–2379. doi: 10.1056/NEJMra1600266 pmid: 27974040
52. Nielsen HB, Almeida M, Juncker AS, Rasmussen S, Li J, Sunagawa S, et al. Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes. Nature Biotechnology. 2014;32(8):822–828. doi: 10.1038/nbt.2939 pmid: 24997787
53. Li J, Jia H, Cai X, Zhong H, Feng Q, Sunagawa S, et al. An integrated catalog of reference genes in the human gut microbiome. Nature Biotechnology. 2014;32(8):834–841. doi: 10.1038/nbt.2942 pmid: 24997786
54. Glasser AL, Boudeau J, Barnich N, Perruchot MH, Colombel JF, Darfeuille-Michaud A, Adherent invasive Escherichia coli strains from patients with Crohn’s disease survive and replicate within macrophages without inducing host cell death. Infection and Immunity. 2001;69(9):5529–37. doi: 10.1128/IAI.69.9.5529-5537.2001 pmid: 11500426
55. Barnich N, Boudeau J, Claret L, Darfeuille-Michaud A, Regulatory and functional co-operation of flagella and type 1 pili in adhesive and invasive abilities of AIEC strain LF82 isolated from a patient with Crohn’s disease. Molecular Microbiology. 2003;48(3):781–794. doi: 10.1046/j.1365-2958.2003.03468.x pmid: 12694621
56. Small CLN, Reid-Yu SA, McPhee JB, Coombes BK, Persistent infection with Crohn’s disease-associated adherent-invasive Escherichia coli leads to chronic inflammation and intestinal fibrosis. Nature Communications. 2013;4:1957. doi: 10.1038/ncomms2957 pmid: 23748852
57. Lawley TD, Klimke WA, Gubbins MJ, Frost LS, F factor conjugation is a true type IV secretion system. FEMS Microbiology Letters. 2003;224(1):1–15. doi: 10.1016/S0378-1097(03)00430-0 pmid: 12855161
58. Stecher B, Denzler R, Maier L, Bernet F, Sanders MJ, Pickard DJ, et al. Gut inflammation can boost horizontal gene transfer between pathogenic and commensal Enterobacteriaceae. Proceedings of the National Academy of Sciences. 2012;109(4):1269–1274. doi: 10.1073/pnas.1113246109
59. Ives AR, Helmus MR, Generalized linear mixed models for phylogenetic analyses of community structure. Ecological Monographs. 2011;81(3):511–525. doi: 10.1890/10-1264.1
60. Konstantinidis KT, Tiedje JM, Genomic insights that advance the species definition for prokaryotes. Proceedings of the National Academy of Sciences. 2005;102(7):2567–72. doi: 10.1073/pnas.0409727102
61. Jain C, Rodriguez-R LM, Phillippy AM, Konstantinidis KT, Aluru S. High-throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. bioRxiv. 2017; p. 225342.
62. Rosselló-Móra R, Amann R, Past and future species definitions for Bacteria and Archaea. Systematic and Applied Microbiology. 2015;38(4):209–216. doi: 10.1016/j.syapm.2015.02.001 pmid: 25747618
63. Cadillo-Quiroz H, Didelot X, Held NL, Herrera A, Darling A, Reno ML, et al. Patterns of gene flow define species of thermophilic archaea. PLoS Biology. 2012;10(2):e1001265. doi: 10.1371/journal.pbio.1001265 pmid: 22363207
64. Richter M, Rosselló-Móra R, Shifting the genomic gold standard for the prokaryotic species definition. Proceedings of the National Academy of Sciences. 2009;106(45):19126–31. doi: 10.1073/pnas.0906412106
65. Winter DJ. rentrez: An R package for the NCBI eUtils API. PeerJ Preprints. 2017
66. Edgar RC, Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–2461. doi: 10.1093/bioinformatics/btq461 pmid: 20709691
67. Wattam AR, Abraham D, Dalay O, Disz TL, Driscoll T, Gabbard JL, et al. PATRIC, the bacterial bioinformatics database and analysis resource. Nucleic Acids Research. 2014;42(Database issue):D581–91. doi: 10.1093/nar/gkt1099 pmid: 24225323
68. Price MN, Dehal PS, Arkin AP, FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS ONE. 2010;5(3):e9490. doi: 10.1371/journal.pone.0009490 pmid: 20224823
69. Edgar RC, MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research. 2004;32(5):1792–1797. doi: 10.1093/nar/gkh340 pmid: 15034147
70. Paradis E, Claude J, Strimmer K, APE: Analyses of Phylogenetics and Evolution in R language. Bioinformatics. 2004;20(2):289–290. doi: 10.1093/bioinformatics/btg412 pmid: 14734327
71. Garamszegi LZ, Gonzalez-Voyer A, Working with the tree of life in comparative studies: how to build and tailor phylogenies to interspecific datasets. In: Modern phylogenetic comparative methods and their application in evolutionary biology. Berlin, Heidelberg: SpringerBerlin Heidelberg; 2014. p. 19–48. Available from: http://link.springer.com/10.1007/978-3-662-43550-2_2.
72. Karlsson FH, Tremaroli V, Nookaew I, Bergström G, Behre CJ, Fagerberg B, et al. Gut metagenome in European women with normal, impaired and diabetic glucose control. Nature. 2013;498(7452):99–103. doi: 10.1038/nature12198 pmid: 23719380
73. Qin J, Li Y, Cai Z, Li S, Zhu J, Zhang F, et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature. 2012;490(7418):55–60. doi: 10.1038/nature11450 pmid: 23023125
74. Zhu Y, Stephens RM, Meltzer PS, Davis SR, SRAdb: query and use public next-generation sequencing data from within R. BMC Bioinformatics. 2013;14(1):19. doi: 10.1186/1471-2105-14-19 pmid: 23323543
75. Leinonen R, Akhtar R, Birney E, Bower L, Cerdeno-Tarraga A, Cheng Y, et al. The European Nucleotide Archive. Nucleic Acids Research. 2011;39(Database):D28–D31. doi: 10.1093/nar/gkq967 pmid: 20972220
76. si Tung Ho L, Ané C, A linear-time algorithm for Gaussian and non-Gaussian trait evolution models. Systematic Biology. 2014;63(3):397–408. doi: 10.1093/sysbio/syu005
77. Storey JD, Bass AJ, Dabney A, Robinson D. qvalue: Q-value estimation for false discovery rate control; 2015. Available from: http://github.com/jdstorey/qvalue.
78. Storey JD, Tibshirani R, Statistical significance for genomewide studies. Proceedings of the National Academy of Sciences. 2003;100(16):9440–5. doi: 10.1073/pnas.1530509100
79. Kappelman MD, Rifas-Shiman SL, Kleinman K, Ollendorf D, Bousvaros A, Grand RJ, et al. The prevalence and geographic distribution of Crohn’s disease and ulcerative colitis in the United States. Clinical Gastroenterology and Hepatology. 2007;5(12):1424–1429. doi: 10.1016/j.cgh.2007.07.012 pmid: 17904915
80. Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Teamnlme: Linear and Nonlinear Mixed Effects Models; 2018. Available from: https://CRAN.R-project.org/package=nlme.
81. Delignette-Muller ML, Dutang C, fitdistrplus: An R Package for Fitting Distributions. Journal of Statistical Software. 2015;64(4):1–34. doi: 10.18637/jss.v064.i04
82. Gloor GB, Macklaim JM, Pawlowsky-Glahn V, Egozcue JJ, Microbiome datasets are compositional: and this is not optional. Frontiers in Microbiology. 2017;8:2224. doi: 10.3389/fmicb.2017.02224 pmid: 29187837
83. Fernandes AD, Reid JN, Macklaim JM, McMurrough TA, Edgell DR, Gloor GB, Unifying the analysis of high-throughput sequencing datasets: characterizing RNA-seq, 16S rRNA gene sequencing and selective growth experiments by compositional data analysis. Microbiome. 2014;2(1):15. doi: 10.1186/2049-2618-2-15 pmid: 24910773
84. Aitchison J, The statistical analysis of compositional data. Journal of the Royal Statistical Society, Series B (Methodological). 1982;44(2):139–177.
85. Friedman J, Alm EJ, Inferring correlation networks from genomic survey data. PLoS Comput Biol. 2012;8. doi: 10.1371/journal.pcbi.1002687
86. Weiss S, Van Treuren W, Lozupone C, Faust K, Friedman J, Deng Y, et al. Correlation detection strategies in microbial data sets vary widely in sensitivity and precision. ISME J. 2016;10(7). doi: 10.1038/ismej.2015.235 pmid: 26905627
87. Xu L, Paterson AD, Turpin W, Xu W, Assessment and selection of competing models for zero-inflated microbiome data. PLoS ONE. 2015;10(7):e0129606. doi: 10.1371/journal.pone.0129606 pmid: 26148172
88. Overbeek R, Begley T, Butler RM, Choudhuri JV, Chuang HY, Cohoon M, et al. The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucleic Acids Research. 2005;33(17):5691–5702. doi: 10.1093/nar/gki866 pmid: 16214803
89. Hochberg Y, Benjamini Y, Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B (Methodological). 1995;1:289–300.
90. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences of the United States of America. 2005;102(43):15545–50. doi: 10.1073/pnas.0506580102 pmid: 16199517
91. R Core TeamR: A language and environment for statistical computing; 2017. Available from: https://www.R-project.org/.
92. Allaire J, Xie Y, McPherson J, Luraschi J, Ushey K, Atkins A, et al. rmarkdown: Dynamic Documents for R; 2018. Available from: https://CRAN.R-project.org/package=rmarkdown.
93. McMurdie PJ, Holmes S, phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8(4):e61217. doi: 10.1371/journal.pone.0061217 pmid: 23630581
94. Solymos P, Zawadzki Z. pbapply: Adding Progress Bar to’*apply’ Functions; 2018. Available from: https://CRAN.R-project.org/package=pbapply.
95. Lang DT, the CRAN Team. XML: Tools for Parsing and Generating XML Within R and S-Plus; 2018. Available from: https://CRAN.R-project.org/package=XML.
96. Bache SM, Wickham H. magrittr: A Forward-Pipe Operator for R; 2014. Available from: https://CRAN.R-project.org/package=magrittr.
97. Revell LJ, phytools: An R package for phylogenetic comparative biology (and other things). Methods in Ecology and Evolution. 2012;3:217–223. doi: 10.1111/j.2041-210X.2011.00169.x
98. Yu G, Smith D, Zhu H, Guan Y, Lam TTY, ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods in Ecology and Evolution. 2017;8:28–36. doi: 10.1111/2041-210X.12628
99. Auguie B. gridExtra: miscellaneous functions for “grid” graphics; 2017. Available from: https://CRAN.R-project.org/package=gridExtra.
100. Venables WN, Ripley BD. Modern Applied Statistics with S. 4th ed. New York: Springer; 2002. Available from: http://www.stats.ox.ac.uk/pub/MASS4.
101. Mersmann O, Trautmann H, Steuer D, Bornkamp B. truncnorm: truncated normal distribution; 2018. Available from: https://CRAN.R-project.org/package=truncnorm.
102. Paradis E, Claude J, Strimmer K, APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–290. doi: 10.1093/bioinformatics/btg412 pmid: 14734327
103. Harmon LJ, Weir JT, Brock CD, Glor RE, Challenger W, GEIGER: investigating evolutionary radiations. Bioinformatics. 2008;24:129–131. doi: 10.1093/bioinformatics/btm538 pmid: 18006550
104. Martin AD, Quinn KM, Park JH, MCMCpack: Markov Chain Monte Carlo in R. Journal of Statistical Software. 2011;42(9):22. doi: 10.18637/jss.v042.i09
105. Xie Y, Stodden V, Leisch F, Peng RD, knitr: A comprehensive tool for reproducible research in R. In: Implementing reproducible computational research. Chapman and Hall/CRC; 2014. Available from: http://www.crcpress.com/product/isbn/9781466561595.
106. Dowle M, Srinivasan A. data.table: Extension of ‘data.frame‘; 2018. Available from: https://CRAN.R-project.org/package=data.table.
107. Wickham H, Francois R, Henry L, Müller K. dplyr: a grammar of data manipulation; 2017. Available from: https://CRAN.R-project.org/package=dplyr.