The first motivation of this paper is to study stationarity and ergodic properties for a general class of time series models defined conditional on an exogenous covariates process. The dynamic of these models is given by an autoregressive latent process which forms a Markov chain in random environments. Contrarily to existing contributions in the field of Markov chains in random environments, the state space is not discrete and we do not use small set type assumptions or uniform contraction conditions for the random Markov kernels. Our assumptions are quite general and allow us to deal with models that are not fully contractive, such as threshold autoregressive processes. Using a coupling approach, we study the existence of a limit, in Wasserstein metric, for the backward iterations of the chain. We also derive ergodic properties for the corresponding skew-product Markov chain. Our results are illustrated with many examples of autoregressive processes widely used in statistics or in econometrics, including GARCH type processes, count autoregressions and categorical time series.