A non-negative matrix factorization framework for identifying modular patterns in metagenomic profile data