We investigate the effects of the nonminimal coupling between the scalar field dark energy (quintessence) and the dark matter on the two-point correlation function. It is well known that this coupling shifts the turnover scale as well as suppresses the amplitude of the matter power spectrum. However, these effects are too small to be observed when we limit the coupling strength to be consistent with observations. Since the coupling of quintessence to baryons is strongly constrained, species-dependent coupling may arise. This results in a baryon bias that is different from unity. Thus, we investigate the correlation function in this coupled model. We are able to observe the enhancement of the baryon acoustic oscillation (BAO) peak due to the increasing bias factor of baryon from this species-dependent coupling. In order to avoid the damping effect of the BAO signature in the matter power spectrum due to nonlinear clustering, we consider the coupling effect on the BAO bump in the linear regime. This provides an alternative method to constrain the coupling of dark energy to dark matter. PACS numbers: 95.36.+x, 95.35.+d, 98.80.Cq Because of the strong constraint on the coupling of the scalar field dark energy (quintessence) to baryons from the local gravity, we investigate the effect of the speciesdependent coupling [1,2] by considering a model in which the quintessence Q is only coupled to the cold dark matter (CDM). We assume a Yukawa-type coupling, m c ¼ e ncQ m à c , where m à c is the bare mass of the CDM [3,4]. This specific choice of coupling requires that the present value of the scalar field vanishes in order to satisfy m c ¼ m à c at present. Then, we are able to write the general action including this interaction asis the reduced Planck mass, VðQÞ is the potential of Q, and L r anddenote the Lagrangian of radiation, CDM, and baryons, respectively. We adopt VðQÞ ¼ V 0 expðQ 2 =2Þ with ¼ 5 in the following [4,5]. However, the main conclusions are independent of the form of the scalar field potential (see below). Because of the coupling, the scalings of the CDM and the quintessence energy densities are changed, respectively, to [5] c ðaÞ ¼ 0 c a À3þ ; where lnðaÞ ¼ n c ½QðaÞ À Qð1Þ;where H ðda=dÞ=a, ! Q is the equation of state of the quintessential dark energy (DE), 0 c denotes the present value of the CDM energy density, the present value of scale factor a 0 ¼ 1, and primes mean the differentiation with respect to the conformal time . Generally, the sign of depends on both the model and the form of the coupling. The linear perturbation equations for the CDM and the scalar field Q in the synchronous gauge are [6] 0 c ¼ À c À 1 2 h 0 þ n c Q 0 ; (4) 0 c ¼ ÀH c þ n c ðk 2 Q À Q 0 c Þ;Q 00 þ 2H Q 0 þ k 2 Q þ a 2 " M 2