Abstract. The paper is concerned with the completeness property of root functions of general boundary value problems for n × n first order systems of ordinary differential equations on a finite interval. In comparison with the recent paper [45] we substantially relax the assumptions on boundary conditions guarantying the completeness of root vectors, allowing them to be nonweakly regular and even degenerate. Emphasize that in this case the completeness property substantially depends on the values of a potential matrix at the endpoints of the interval.It is also shown that the system of root vectors of the general n × n Dirac type system subject to certain boundary conditions forms a Riesz basis with parentheses. We also show that arbitrary complete dissipative boundary value problem for Dirac type operator with a summable potential matrix admits the spectral synthesis in L 2 ([0, 1]; C n ). Finally, we apply our results to investigate completeness and the Riesz basis property of the dynamic generator of spatially non-homogenous damped Timoshenko beam model.