phyloseq object

OTU table should be prepared, metadata and taxonomic table details

Example: phyloseq

The diet swap data set represents a study with African and African American groups undergoing a two-week diet swap. For details, see https://www.nature.com/articles/ncomms7342.

phyloseq-class experiment-level object
otu_table()   OTU Table:         [ 130 taxa and 222 samples ]
sample_data() Sample Data:       [ 222 samples by 8 sample variables ]
tax_table()   Taxonomy Table:    [ 130 taxa by 3 taxonomic ranks ]

Looking into the data

  • Number of taxonomic groups: 130

  • Number of samples: 222

  • Number of metadata fields: 8:

  • Metadata variables are: subject, sex, nationality, group, sample, timepoint, timepoint.within.group, bmi_group

Unique taxa at Phylum:

phyloseq::get_taxa_unique(phy, “Phylum”)
Actinobacteria
Firmicutes
Proteobacteria
Verrucomicrobia
Bacteroidetes
Spirochaetes
Fusobacteria
Cyanobacteria

Unique taxa at Genus:

phyloseq::get_taxa_unique(phy, “Genus”)
Actinomycetaceae
Aerococcus
Aeromonas
Akkermansia
Alcaligenes faecalis et rel.
Allistipes et rel.
Anaerobiospirillum
Anaerofustis
Anaerostipes caccae et rel.
Anaerotruncus colihominis et rel.
Anaerovorax odorimutans et rel.
Aneurinibacillus
Aquabacterium
Asteroleplasma et rel.
Atopobium
Bacillus
Bacteroides fragilis et rel.
Bacteroides intestinalis et rel.
Bacteroides ovatus et rel.
Bacteroides plebeius et rel.
Bacteroides splachnicus et rel.
Bacteroides stercoris et rel.
Bacteroides uniformis et rel.
Bacteroides vulgatus et rel.
Bifidobacterium
Bilophila et rel.
Brachyspira
Bryantella formatexigens et rel.
Bulleidia moorei et rel.
Burkholderia
Butyrivibrio crossotus et rel.
Campylobacter
Catenibacterium mitsuokai et rel.
Clostridium (sensu stricto)
Clostridium cellulosi et rel.
Clostridium colinum et rel.
Clostridium difficile et rel.
Clostridium felsineum et rel.
Clostridium leptum et rel.
Clostridium nexile et rel.
Clostridium orbiscindens et rel.
Clostridium ramosum et rel.
Clostridium sphenoides et rel.
Clostridium stercorarium et rel.
Clostridium symbiosum et rel.
Clostridium thermocellum et rel.
Collinsella
Coprobacillus catenaformis et rel.
Coprococcus eutactus et rel.
Corynebacterium
Desulfovibrio et rel.
Dialister
Dorea formicigenerans et rel.
Eggerthella lenta et rel.
Enterobacter aerogenes et rel.
Enterococcus
Escherichia coli et rel.
Eubacterium biforme et rel.
Eubacterium cylindroides et rel.
Eubacterium hallii et rel.
Eubacterium limosum et rel.
Eubacterium rectale et rel.
Eubacterium siraeum et rel.
Eubacterium ventriosum et rel.
Faecalibacterium prausnitzii et rel.
Fusobacteria
Gemella
Granulicatella
Haemophilus
Helicobacter
Klebisiella pneumoniae et rel.
Lachnobacillus bovis et rel.
Lachnospira pectinoschiza et rel.
Lactobacillus catenaformis et rel.
Lactobacillus gasseri et rel.
Lactobacillus plantarum et rel.
Lactobacillus salivarius et rel.
Lactococcus
Leminorella
Megamonas hypermegale et rel.
Megasphaera elsdenii et rel.
Methylobacterium
Micrococcaceae
Mitsuokella multiacida et rel.
Moraxellaceae
Novosphingobium
Oceanospirillum
Oscillospira guillermondii et rel.
Outgrouping clostridium cluster XIVa
Oxalobacter formigenes et rel.
Papillibacter cinnamivorans et rel.
Parabacteroides distasonis et rel.
Peptococcus niger et rel.
Peptostreptococcus anaerobius et rel.
Peptostreptococcus micros et rel.
Phascolarctobacterium faecium et rel.
Prevotella melaninogenica et rel.
Prevotella oralis et rel.
Prevotella ruminicola et rel.
Prevotella tannerae et rel.
Propionibacterium
Proteus et rel.
Pseudomonas
Roseburia intestinalis et rel.
Ruminococcus bromii et rel.
Ruminococcus callidus et rel.
Ruminococcus gnavus et rel.
Ruminococcus lactaris et rel.
Ruminococcus obeum et rel.
Serratia
Sporobacter termitidis et rel.
Staphylococcus
Streptococcus bovis et rel.
Streptococcus intermedius et rel.
Streptococcus mitis et rel.
Subdoligranulum variable at rel.
Sutterella wadsworthia et rel.
Tannerella et rel.
Uncultured Bacteroidetes
Uncultured Chroococcales
Uncultured Clostridiales I
Uncultured Clostridiales II
Uncultured Mollicutes
Uncultured Selenomonadaceae
Veillonella
Vibrio
Weissella et rel.
Wissella et rel.
Xanthomonadaceae
Yersinia et rel.

Useful plots

The taxa ranks in the previous phyloseq dataset are: Phylum, Family, Genus

Important to remember taxa abundances are compositional that are need appropriate transform to be informative, here I will use the relative abundances via compositional transform.

Abundance at Phylum level

Abundance at Family level

Abundance at Genus level

Heatmaps for the core taxa with certain prevalences & detections

Core heatmaps

Acknowledgement to "Leo Lahti, Sudarshan Shetty et al. (2017). Tools for microbiome analysis in R. Version . URL: http://microbiome.github.com/microbiome

LS0tDQp0aXRsZTogIm92ZXJ2aWV3Ig0KYXV0aG9yOiAiV2lzYW0iDQpkYXRlOiAiYHIgU3lzLkRhdGUoKWAiDQpvdXRwdXQ6IA0KICBodG1sX2RvY3VtZW50Og0KICAgIHRvYzogdHJ1ZQ0KICAgIHRvY19kZXB0aDogNA0KICAgIHRvY19mbG9hdDogdHJ1ZQ0KICAgIGNvZGVfZG93bmxvYWQ6IHRydWUNCi0tLQ0KDQpgYGB7ciBwYWNrYWdlcywgZWNobz1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgZXJyb3I9RkFMU0UsIHdhcm5pbmc9RkFMU0V9DQpsaWJyYXJ5KGtuaXRyKQ0KbGlicmFyeShwaHlsb3NlcSkNCmxpYnJhcnkodmVnYW4pDQpsaWJyYXJ5KGdncGxvdDIpDQpsaWJyYXJ5KHZlZ2FuKQ0KbGlicmFyeSh0aWN0b2MpDQpsaWJyYXJ5KHBlcm11dGUpDQpsaWJyYXJ5KGtuaXRyKQ0KbGlicmFyeShsYXR0aWNlKQ0Kb3B0aW9ucyhtYXgucHJpbnQ9Ijc1IikNCiAga25pdHI6Om9wdHNfY2h1bmskc2V0KGZpZy53aWR0aD04LA0KICAgICAgICAgICAgICAgICAgICAgICAgZmlnLmhlaWdodD02LA0KICAgICAgICAgICAgICAgICAgICAgICAgZXZhbD1UUlVFLA0KICAgICAgICAgICAgICAgICAgICAgICAgY2FjaGU9VFJVRSwNCiAgICAgICAgICAgICAgICAgICAgICAgIGVjaG89VFJVRSwNCiAgICAgICAgICAgICAgICAgICAgICAgIHByb21wdD1GQUxTRSwNCiAgICAgICAgICAgICAgICAgICAgICAgIHRpZHk9VFJVRSwNCiAgICAgICAgICAgICAgICAgICAgICAgIGNvbW1lbnQ9TkEsDQogICAgICAgICAgICAgICAgICAgICAgICBtZXNzYWdlPUZBTFNFLA0KICAgICAgICAgICAgICAgICAgICAgICAgd2FybmluZz1GQUxTRSkNCm9wdHNfa25pdCRzZXQod2lkdGg9NzUpDQoNCmBgYA0KDQojIHBoeWxvc2VxIG9iamVjdA0KDQpPVFUgdGFibGUgc2hvdWxkIGJlIHByZXBhcmVkLCBtZXRhZGF0YSBhbmQgdGF4b25vbWljIHRhYmxlIFtkZXRhaWxzXShodHRwczovL2pvZXk3MTEuZ2l0aHViLmlvL3BoeWxvc2VxL2ltcG9ydC1kYXRhLmh0bWwjZXhhbXBsZXMpDQoNCiMjIEV4YW1wbGU6IHBoeWxvc2VxDQoNClRoZSBkaWV0IHN3YXAgZGF0YSBzZXQgcmVwcmVzZW50cyBhIHN0dWR5IHdpdGggQWZyaWNhbiBhbmQgQWZyaWNhbiBBbWVyaWNhbiBncm91cHMgdW5kZXJnb2luZyBhIHR3by13ZWVrIGRpZXQgc3dhcC4gRm9yIGRldGFpbHMsIHNlZSBodHRwczovL3d3dy5uYXR1cmUuY29tL2FydGljbGVzL25jb21tczczNDIuDQoNCg0KYGBge3IgZGF0YXNldCwgZWNobz1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgZXJyb3I9RkFMU0UsIHdhcm5pbmc9RkFMU0V9DQpsaWJyYXJ5KG1pY3JvYmlvbWUpDQpkYXRhKCJkaWV0c3dhcCIpDQpwaHk8LSBkaWV0c3dhcA0KcGh5DQpgYGANCg0KIyMjIExvb2tpbmcgaW50byB0aGUgZGF0YQ0KDQogKiBOdW1iZXIgb2YgdGF4b25vbWljIGdyb3VwczogYHIgbnRheGEocGh5KWANCiAqIE51bWJlciBvZiBzYW1wbGVzOiBgciBuc2FtcGxlcyhwaHkpYA0KICogTnVtYmVyIG9mIG1ldGFkYXRhIGZpZWxkczogYHIgbmNvbChnZXRfdmFyaWFibGUocGh5KSlgOiAgDQoNCiAqIE1ldGFkYXRhIHZhcmlhYmxlcyBhcmU6IGByIHNhbXBsZV92YXJpYWJsZXMocGh5KWANCg0KDQojIyMjIFVuaXF1ZSB0YXhhIGF0IFBoeWx1bToNCg0KYGBge3IgcmFua3NfcGh5bHVtLCBlY2hvPUZBTFNFLCBtZXNzYWdlPUZBTFNFLCBlcnJvcj1GQUxTRX0NCmQgPC0gYXMuZGF0YS5mcmFtZShwaHlsb3NlcTo6Z2V0X3RheGFfdW5pcXVlKHBoeSwgJ1BoeWx1bScpKQ0Ka2FibGUoZCkNCmBgYA0KDQojIyMjIFVuaXF1ZSB0YXhhIGF0IEdlbnVzOg0KDQpgYGB7ciByYW5rcywgZWNobz1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgZXJyb3I9RkFMU0V9DQpkIDwtIGFzLmRhdGEuZnJhbWUocGh5bG9zZXE6OmdldF90YXhhX3VuaXF1ZShwaHksICdHZW51cycpKQ0Ka2FibGUoZCkNCg0KYGBgDQoNCg0KDQojIFVzZWZ1bCBwbG90cw0KDQpUaGUgdGF4YSByYW5rcyBpbiB0aGUgcHJldmlvdXMgcGh5bG9zZXEgZGF0YXNldCBhcmU6IGByIHJhbmtfbmFtZXMocGh5KWANCg0KKipJbXBvcnRhbnQgdG8gcmVtZW1iZXIqKiB0YXhhIGFidW5kYW5jZXMgYXJlIGNvbXBvc2l0aW9uYWwgdGhhdCBhcmUgbmVlZCBhcHByb3ByaWF0ZSB0cmFuc2Zvcm0gdG8gYmUgaW5mb3JtYXRpdmUsIGhlcmUgSSB3aWxsIHVzZSB0aGUgcmVsYXRpdmUgYWJ1bmRhbmNlcyB2aWEgKipjb21wb3NpdGlvbmFsKiogdHJhbnNmb3JtLg0KDQojIyBBYnVuZGFuY2UgYXQgUGh5bHVtIGxldmVsDQoNCmBgYHtyIHBoeWx1bSwgZWNobz1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgZmlnLndpZHRoPTEwLCBmaWcuaGVpZ2h0PTgsIGZpZy5zaG93PSJrZWVwIiwgb3V0LndpZHRoPSI1MCUifQ0KcGhfY29tcCA8LSB0cmFuc2Zvcm1fc2FtcGxlX2NvdW50cyhwaHksIGZ1bmN0aW9uKE9UVSkgT1RVL3N1bShPVFUpKQ0KDQpwbG90X2JhcihwaF9jb21wLCAnUGh5bHVtJywgZmlsbCA9ICdibWlfZ3JvdXAnKSArIHlsYWIoIlJlbGF0aXZlIGFidW5kYW5jZSIpDQpgYGANCg0KIyMgQWJ1bmRhbmNlIGF0IEZhbWlseSBsZXZlbA0KDQpgYGB7ciBGYW1pbHksIGVjaG89RkFMU0UsIG1lc3NhZ2U9RkFMU0UsIGZpZy53aWR0aD0xMCwgZmlnLmhlaWdodD04LCBmaWcuc2hvdz0ia2VlcCIsIG91dC53aWR0aD0iNTAlIn0NCnBsb3RfYmFyKHBoX2NvbXAsICdGYW1pbHknLCBmaWxsID0gJ2JtaV9ncm91cCcpK3lsYWIoIlJlbGF0aXZlIGFidW5kYW5jZSIpDQpgYGANCg0KIyMgQWJ1bmRhbmNlIGF0IEdlbnVzIGxldmVsDQoNCmBgYHtyIEdlbnVzLCBlY2hvPUZBTFNFLCBtZXNzYWdlPUZBTFNFLCBlcnJvcj1GQUxTRX0NCm15X3ggPSB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChzaXplID0gNSkpDQpwPC1wbG90X2JhcihwaF9jb21wLCAnR2VudXMnLCBmaWxsID0gJ2JtaV9ncm91cCcpK3lsYWIoIlJlbGF0aXZlIGFidW5kYW5jZSIpDQpwK215X3gNCmBgYA0KDQojIyBIZWF0bWFwcyBmb3IgdGhlIGNvcmUgdGF4YSB3aXRoIGNlcnRhaW4gKipwcmV2YWxlbmNlcyoqICYgKipkZXRlY3Rpb25zKiogDQoNCg0KQ29yZSBoZWF0bWFwcw0KDQoqKkFja25vd2xlZGdlbWVudCoqIHRvICJMZW8gTGFodGksIFN1ZGFyc2hhbiBTaGV0dHkgZXQgYWwuICgyMDE3KS4gVG9vbHMgZm9yIG1pY3JvYmlvbWUgYW5hbHlzaXMgaW4gUi4gVmVyc2lvbiAuIFVSTDogaHR0cDovL21pY3JvYmlvbWUuZ2l0aHViLmNvbS9taWNyb2Jpb21lDQoNCg0KDQpgYGB7ciBoZWF0bWFwcywgZWNobz1GQUxTRSwgbWVzc2FnZT1GQUxTRSwgZXJyb3I9RkFMU0UsIHJlc3VsdHM9J2hpZGUnfQ0KbGlicmFyeShSQ29sb3JCcmV3ZXIpDQpsaWJyYXJ5KHJlc2hhcGUpDQpwIDwtbWljcm9iaW9tZTo6cGxvdF9jb3JlKHBoX2NvbXAsIHBsb3QudHlwZSA9ICJoZWF0bWFwIiwgY29sb3VycyA9IGdyYXkoc2VxKDAsMSxsZW5ndGg9NSkpLCBwcmV2YWxlbmNlcyA9IHNlcSguMDUsIDEsIC4wNSksIGRldGVjdGlvbnMgPSAxMF5zZXEobG9nMTAoMWUtMyksIGxvZzEwKC4yKSwgbGVuZ3RoID0gMTApKSArDQogICAgbGFicyh4ID0gIkRldGVjdGlvbiBUaHJlc2hvbGQgKFJlbGF0aXZlIEFidW5kYW5jZSAoJSkpIikNCm15X3k8LXRoZW1lKGF4aXMudGV4dC55ID0gZWxlbWVudF90ZXh0KHNpemUgPSA0KSkNCmBgYA0KDQoNCmBgYHtyIGhlYXRtYXBzX3Bsb3QsIGVjaG89RkFMU0UsIG1lc3NhZ2U9RkFMU0UsIGVycm9yPUZBTFNFfQ0KcCArIHRoZW1lKGF4aXMudGl0bGUueCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTQpLA0KICAgICAgICAgIGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KHNpemUgPSA1KSwNCiAgICAgICAgICAgICAgICAgYXhpcy50ZXh0LnkgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDQpLA0KICAgICAgICAgICAgICAgICAgICAgIGxlZ2VuZC50aXRsZT1lbGVtZW50X3RleHQoc2l6ZT0xMCksbGVnZW5kLnRleHQ9ZWxlbWVudF90ZXh0KHNpemU9OSkpDQpgYGANCg0KDQo=