In this paper we consider time-changed models of population evolution Xf(t) = X(Hf(t)), where X is a counting process and Hf is a subordinator with Laplace exponent f. In the case where X is a pure birth process, we study the form of the distribution, the intertimes between successive jumps, and the condition of explosion (also in the case of killed subordinators). We also investigate the case where X represents a death process (linear or sublinear) and study the extinction probabilities as a function of the initial population size n0. Finally, the subordinated linear birth–death process is considered. Special attention is devoted to the case where birth and death rates coincide; the sojourn times are also analysed.