We present the third-order analytic solution of the matter density fluctuation in the proper-time hypersurface of nonrelativistic matter flows by solving the nonlinear general relativistic equations. The proper-time hypersurface provides a coordinate system that a local observer can set up without knowledge beyond its neighborhood, along with physical connections to the local Newtonian descriptions in the relativistic context. The initial condition of our analytic solution is set up by the curvature perturbation in the comoving gauge, clarifying its impact on the nonlinear evolution. We compute the effective non-Gaussian parameters due to the nonlinearity in the relativistic equations. With proper coordinate rescaling, we show that the equivalence principle is respected and the relativistic effect vanishes in the large-scale limit. 98.80.Jk,98.62.Py Rapid developments in large-scale galaxy surveys over the past decades have enabled the precision measurements of galaxy clustering, which can be used to probe the nature of dark energy and the perturbation generation mechanism in the early Universe [1]. In parallel, the recent theoretical development ([2-6]; see [7] for review) has revealed that the subtle relativistic effects are present in galaxy clustering, providing new opportunities to extract additional and critical information about the gravity on large scales and the initial conditions for structure formation. In particular, the relativistic formalism has been extended [8-10] to the second-order in perturbations for the computation of higher-order statistics such as the bispectrum.One of the critical elements in the relativistic formalism is galaxy bias, which relates the galaxy number density to the underlying matter distribution. Beyond the linear order in perturbations, however, galaxy bias poses a nontrivial problem due to the gauge issues in general relativity. It was shown [11] that the proper-time hypersurface of nonrelativistic matter flows provides a physical description of the local observer, moving with dark matter and baryons that will collapse to form galaxies. This physical justification has led us to study the nonlinear relativistic effects of the matter density fluctuation in the proper-time hypersurface of nonrelativistic matter flows. This is timely and interesting, for there has been recent conflict in literature -It is argued in [12,13] by performing the second-order relativistic calculations that the nonlinear evolution of gravity generates the local-type non-Gaussianity, which would in turn exhibit a prominent signature in galaxy clustering on large scales. On the other hands, in [14-16] the opposite claim is asserted that the observable quantities are not affected by these nonlinear relativistic effects of gravity, while the calculations are in general based on studying the special case (or the squeezed limit), where only the linear- * Electronic address: jyoo@physik.uzh.ch † Electronic address: jinn-ouk.gong@apctp.org order calculations are required. In this article, we present the ...