  We study matrix models involving Pfaffian interactions as generalizations of
the standard $\beta = 1$ and $\beta = 4$ matrix models. We present the Pfaffian
formulas for the partition function and the characteristic polynomial averages.
We also explore the matrix chain with the Pfaffian interaction, which realizes
the BCD-type quiver matrix models.
