We study the pseudo-Hermitian systems with general spin-coupling point interactions and give a systematic description of the corresponding boundary conditions for PT-symmetric systems. The corresponding integrability for both bosonic and fermionic many-body systems with PT-symmetric contact interactions is investigated.Self-adjoint quantum mechanical models describing a particle moving in a local singular potential have been extensively discussed [1][2][3]. The integrability of one dimensional many-body systems with self-adjoint contact interactions has been studied according to Yang-Baxter relations [4]. The results are generalized to the case of particles with spin-coupling interactions [5].PT-symmetric quantum mechanical models have been studied from some mathematical and physical considerations [6]. In [7] the classification and spectra problem of PT-symmetric point interactions are investigated. The integrability of manybody systems with PT-symmetric interactions is clarified [8]. The δ-type spincoupling interactions with PT-symmetry is discussed in [9].In this letter we study the boundary conditions for PT-symmetric point interactions of particles with spin-coupling, and the integrability of bosonic and fermionic many-body systems with PT-symmetric, spin-coupling contact interactions characterized by these boundary conditions.One dimensional quantum mechanical models of spinless particles with point interactions at the origin can be characterized by separated or nonseparated boundary conditions imposed on the (scalar) wave function ϕ at x = 0. The family of point interactions for the Schrödinger operator −d 2 /dx 2 can be described by unitary 2 × 2 matrices via von Neumann formulas for self-adjoint extensions of symmetric operators:where ad − bc = 1, θ, a, b, c, d ∈ IR. ϕ(x) is the scalar wave function of two spinless particles with relative coordinate x.(1) also describes two particles with spin s but