Two-dimensional (2D) covariance NMR spectroscopy, which has originally been established to extract homonuclear correlations (HOMCOR), is extended to include heteronuclear correlations (HETCOR). In a C-13/N-15 2D chemical shift correlation experiment, C-13 and N-15 signals of a polycrystalline sample of C-13, N-15-labeled amino acid are acquired simultaneously using a dual-receiver NMR system. The data sets are rearranged for the covariance data processing, and the C-13-N-15 heteronuclear correlations are obtained together with the C-13-C-13 and N-15-N-15 homonuclear correlations. The present approach retains the favorable feature of the original covariance HOMCOR that the spectral resolution along the indirect dimension is given by that of the detection dimension. As a result, much fewer amounts of data are required to obtain a well-resolved 2D spectrum compared to the case of the conventional 2D Fourier-Transformation (FT) scheme. Hence, one can significantly save the experimental time, or enhance the sensitivity by increasing the number of signal averaging within a given measurement time.