Can we use deseq2 for 16s metagenome analysis? If yes is there any specific pipeline for it?

You can take a look at the Bioconductor package phyloseq, and the associated references and workflows from those authors.