We introduce and solve exactly a family of invariant 2 × 2 random matrices, depending on one parameter η, and we show that rotational invariance and real Dyson index β are not incompatible properties. The probability density for the entries contains a weight function and a multiple trace-trace interaction term, which corresponds to the representation of the Vandermonde-squared coupling on the basis of power sums. As a result, the effective Dyson index β eff of the ensemble can take any real value in an interval. Two weight functions (Gaussian and non-Gaussian) are explored in detail and the connections with β-ensembles of Dumitriu-Edelman and the so-called Poisson-Wigner crossover for the level spacing are respectively highlighted. A curious spectral twinning between ensembles of different symmetry classes is unveiled: as a consequence, the identification between symmetry group (orthogonal, unitary or symplectic) and the exponent of the Vandermonde (β = 1, 2, 4) is shown to be potentially deceptive. The proposed technical tool more generically allows for designing actual matrix models which (i) are rotationally invariant; (ii) have a real Dyson index β eff ; (iii) have a pre-assigned confining potential or alternatively levelspacing profile. The analytical results have been checked through numerical simulations with an excellent agreement. Eventually, we discuss possible generalizations and further directions of research.