In recent papers [1,2], two new ways have been proposed to probe the linear polarization of gluons in unpolarized proton: using the azimuthal asymmetries and Callan-Gross ratio in heavy-quark pair leptoproduction, lN → l ′ QQX. In this talk, we discuss in details the sensitivity of the QCD predictions for the azimuthal cos ϕ and cos 2ϕ asymmetries to the contribution of linearly polarized gluons inside unpolarized proton, where the azimuth ϕ is the angle between the lepton scattering plane (l, l ′ ) and the heavy quark production plane (N, Q). Our analysis shows that the azimuthal distributions under consideration vary from 0 to 1 depending on the transverse-momentum dependent gluonic counterpart of the Boer-Mulders function, h ⊥g 1 . We conclude that the cos ϕ and cos 2ϕ asymmetries in heavy-quark pair production in DIS processes are predicted to be large in wide kinematic ranges and sensitive to the contribution of linearly polarized gluons. * e-mail: nikiv@yerphi.am 1 The opportunities to measure the function h ⊥g 1 in unpolarized hadron-hadron collisions are discussed in Ref. [9].