Co-integration is a theoretical, and unconditional relationship between stochastic processes, not a relationship between two finite series of specific realizations of these processes. Therefore when samples of different sizes give different co-integration tests results, the result from the largest sample size prevails.
Using OLS to estimate the co-integrating vector, has the following issues: to apply OLS we have to normalize the co-integrating vector and treat one variable as the "dependent variable" of the specification. If the theory behind the model does not specify which should this variable be, then the choice will be arbitrary - and the OLS estimates are affected by the normalization -and worse, the tests for the existence of co-integration are also affected (see for example Hamilton p. 589).
Second, the usual OLS may have large bias, as mentioned by Hayashi p. 651 . Hayashi calls this estimator "SOLS" (Static OLS) and advocates the "DOLS" estimator (Dynamic OLS), where in the RHS of the specification we include in addition a number of first-differences of the regressors (past and future), which result in the the regressors being strictly exogenous. Studies indicate that the DOLS estimator dominates in root mean squared error both the SOLS estimator, but also the Full-Information Maximum Likelihood estimator of Johansen (1988 and 1995).