The ability to detect and count certain substructures in graphs is important for solving many tasks on graph-structured data, especially in the contexts of computational chemistry and biology as well as social network analysis. Inspired by this, we propose to study the expressive power of graph neural networks (GNNs) via their ability to count attributed graph substructures, extending recent works that examine their power in graph isomorphism testing. We distinguish between two types of substructure counting: matching-count and containment-count, and establish mostly negative answers for a wide class of GNN architectures. Specifically, we prove that Message Passing Neural Networks (MPNNs), Weisfeiler-Lehman (WL) and 2-Invariant Graph Networks (2-IGNs) cannot perform matching-count of substructures consisting of 3 or more nodes, while they can perform containment-count of star-shaped substructures. We also provide partial results for k-WL and k-IGNs. We then conduct experiments that support several of the theoretical results, and demonstrate that local relational pooling strategies inspired by Murphy et al. (2019) are more effective for substructure counting. In addition, we prove that WL and 2-IGNs are equivalent in distinguishing non-isomorphic graphs, partly answering an open problem raised in Maron et al. (2019).