In this paper, we study the nonlinear Klein-Gordon systems arising from relativistic physics and quantum field theories {u(tt) - u(xx) + bu + epsilon v +f(t, x, u) = 0, v(tt) - v(xx) + bv + epsilon u + g(t, x, v) = 0, where u, v satisfy the Dirichlet boundary conditions on spatial interval [0, pi], b > 0 and f, g are 2 pi-periodic in t. We are concerned with the existence, regularity and asymptotic behavior of time-periodic solutions to the linearly coupled problem as e goes to 0. Firstly, under some superlinear growth and monotonicity assumptions on f and g, we obtain the solutions (u(epsilon), v epsilon) with time period 2 pi for the problem as the linear coupling constant e is sufficiently small, by constructing critical points of an indefinite functional via variational methods. Secondly, we give a precise characterization for the asymptotic behavior of these solutions, and show that, as epsilon -> 0, (u(epsilon), v(epsilon)) converge to the solutions of the wave equations without the coupling terms. Finally, by careful analysis which is quite different from the elliptic regularity theory, we obtain some interesting results concerning the higher regularity of the periodic solutions.