Differential gene abundance analysis with 16s seq data

Hello,

Is it possible to do differential gene „abundance“ (I don’t want to say expression because it was not sequenced) analysis from two conditions of 16s RNA seq data (Microbiome analysis)? So basically, by 16s seq we can assign sequences to bacterial taxonomies and then calculate the relative bacterial composition in a given tissue. Now my question is, if we know the genome of the species and if the genes are annotated, then wouldn’t it be possible to derive from that data that say in condition 2 the pathways for post-translational acetylation are increased (at least in relative terms)? Because the pathways consist of genes, not every bacterial species has these genes, thus we could do a GSEA or GO.
I hope that wasn’t too confusing.
Is there such an analysis? Thank you in advance.