  We use the framework of supersymmetric transformations in the construction of
coupled systems of Dirac fermions. Its energy operator is a composite of the
generators of the associated superalgebra, and the two coupled Dirac fermions
acquire different Fermi velocities. We discuss in detail the peculiar spectral
properties of the new system. The emergent phenomena like level crossing or
generation of bound states in the continuum (BICs) are illustrated in two
explicit examples.
