  We develop a Galois theory for systems of linear difference equations with an
action of an endomorphism {\sigma}. This provides a technique to test whether
solutions of such systems satisfy {\sigma}-polynomial equations and, if yes,
then characterize those. We also show how to apply our work to study
isomonodromic difference equations and difference algebraic properties of
meromorphic functions.
