The moduli space of genus g Riemann surfaces equipped with the Teichmuller metric exhibits rich geometric, analytic, and dynamical properties. A major challenge is to understand the totally geodesic submanifolds -- these share many properties with the moduli space itself. For many decades, research focused on the one (complex) dimensional case, i.e. the captivating Teichmuller cuves, but the discovery of interesting higher-dimensional examples in recent years has led to new questions. In this talk, I will discuss joint work with Benirschke and Rached in which we study the boundary of a totally geodesic subvariety in the Deligne-Mumford compactification, showing that the boundary is itself totally geodesic.