  The notion of higher-order topological phases can have interesting
generalizations to systems with subsystem symmetries that exhibit fractonic
dynamics for charged excitations. In this work, we systematically study the
higher-order topological phases protected by a combination of subsystem
symmetries and ordinary global symmetries in two and three-dimensional
interacting boson systems, with some interacting fermionic examples.
