Many-body localized systems in which interactions and disorder come together defy the expectations of quantum statistical mechanics: In contrast to ergodic systems, they do not thermalize when undergoing nonequilibrium dynamics. What is less clear, however, is how topological features interplay with many-body localized phases as well as the nature of the transition between a topological and a trivial state within the latter. In this work, we numerically address these questions, using a combination of extensive tensor network calculations, specifically DMRG-X, as well as exact diagonalization, leading to a comprehensive characterization of Hamiltonian spectra and eigenstate entanglement properties.