```r
OTU = otu_table(as.matrix(abund_table), taxa_are_rows = FALSE)
TAX = tax_table(as.matrix(OTU_taxonomy))
SAM = sample_data(meta_table)
OTU_tree <- compute.brlen(OTU_tree, method = "Grafen")
physeq <- merge_phyloseq(phyloseq(OTU, TAX), SAM, OTU_tree)
```