By using the K-free complex bosons and the K-free complex fermions, we construct the Formula omitted supersymmetric Formula omitted algebra which is the matrix generalization of previous Formula omitted supersymmetric Formula omitted algebra. By twisting this Formula omitted supersymmetric Formula omitted algebra, we obtain the Formula omitted supersymmetric Formula omitted algebra which is the matrix generalization of known Formula omitted supersymmetric topological Formula omitted algebra. From this two-dimensional symmetry algebra, we propose the operator product expansion (OPE) between the soft graviton and gravitino (as a first example), at nonzero deformation parameter, in the supersymmetric Einstein-Yang-Mills theory explicitly. Other six OPEs between the graviton, gravitino, gluon and gluino can be determined completely. At vanishing deformation parameter, we reproduce the known result of Fotopoulos, Stieberger, Taylor and Zhu on the above OPEs via celestial holography.