A series of Barrovian sequence samples ranging from garnet to sillimanite zones were investigated to infer their porphyroblast-forming reactions and mineral inclusion histories. Quartz is overgrown and partly consumed during garnet formation and remains as inclusion-rich layers in porphyroblasts of the garnet zone. Staurolite crystals in the staurolite zone display sharp transitions between inclusion-rich and inclusion-free areas, suggesting two stages of growth with a different role of quartz in each. The inclusion-rich domains formed similarly to those in garnet by simple overgrowth and resorption of matrix minerals, with thermodynamic constraints suggesting that this staurolite-forming reaction required the presence of chloritoid that is now absent from the examined samples. The participation of garnet was limited in staurolite formation, with chloritoid breakdown supplying sufficient material to form the large amounts (c. 25 vol%) of staurolite found in the rock. This reaction produces an excess of SiO2, which leaves the crystal domain as SiO2aq and thus caused the formation of the inclusion-free areas in the staurolite and precipitation of quartz in the matrix. In the sillimanite zone, staurolite is consumed forming new garnet. The newly formed garnet has less quartz inclusions than its core due to a proportionally greater consumption of quartz by the second garnet-forming reaction than by the initial, garnet-grade reactions. Textural and thermodynamic data both suggest that inclusions in these porphyroblasts represent leftovers of a preferentially overgrown matrix than co-products of the porphyroblast-forming reaction.