For each discriminant $D>1$, McMullen constructed the Prym–Teichmüller curves $W_{D}(4)$ and $W_{D}(6)$ in ${\mathcal{M}}_{3}$ and ${\mathcal{M}}_{4}$, which constitute one of the few known infinite families of geometrically primitive Teichmüller curves. In the present paper, we determine for each $D$ the number and type of orbifold points on $W_{D}(6)$. These results, together with a previous result of the two authors in the genus $3$ case and with results of Lanneau–Nguyen and Möller, complete the topological characterisation of all Prym–Teichmüller curves and determine their genus. The study of orbifold points relies on the analysis of intersections of $W_{D}(6)$ with certain families of genus $4$ curves with extra automorphisms. As a side product of this study, we give an explicit construction of such families and describe their Prym–Torelli images, which turn out to be isomorphic to certain products of elliptic curves. We also give a geometric description of the flat surfaces associated to these families and describe the asymptotics of the genus of $W_{D}(6)$ for large $D$.

Fuchsian groups with a modular embedding have the richest arithmetic properties among non-arithmetic Fuchsian groups. But they are very rare, all known examples being related either to triangle groups or to Teichmüller curves. In Part I of this paper we study the arithmetic properties of the modular embedding and develop from scratch a theory of twisted modular forms for Fuchsian groups with a modular embedding, proving dimension formulas, coefficient growth estimates and differential equations. In Part II we provide a modular proof for an Apéry-like integrality statement for solutions of Picard–Fuchs equations. We illustrate the theory on a worked example, giving explicit Fourier expansions of twisted modular forms and the equation of a Teichmüller curve in a Hilbert modular surface. In Part III we show that genus two Teichmüller curves are cut out in Hilbert modular surfaces by a product of theta derivatives. We rederive most of the known properties of those Teichmüller curves from this viewpoint, without using the theory of flat surfaces. As a consequence we give the modular embeddings for all genus two Teichmüller curves and prove that the Fourier developments of their twisted modular forms are algebraic up to one transcendental scaling constant. Moreover, we prove that Bainbridge’s compactification of Hilbert modular surfaces is toroidal. The strategy to compactify can be expressed using continued fractions and resembles Hirzebruch’s in form, but every detail is different.