Let $(a_n)_{n \ge 0}$ be a linear recurrence sequence taking only integer values. Then $a_n$ satisfies a recurrence with integer coefficients.
Notes:
This follows up an older question on this site. The title of the question mentioned "the defining relation". The counterexample was a constant sequence, that can satisfy many other linear recurrences ( multiply the polynomial $(T-1)$ by other poly). However, if a sequence with integral values satisfies an integral linear recurrence, then the defining relation is also integral ( follows from Gauss lemma, and a bit of linear algebra).
What progress I've made: I am able to show the following:
If we have $\alpha_1$, $\ldots$, $\alpha_{\ell}$ , $x_1$, $\ldots$, $x_{\ell}$ distinct complex numbers, such that
$$\sum \alpha_i x_i^n$$
are integral for all $n\ge 0$, then the $x_i$'s are all algebraic integers. From here we can show that the above sums (depending on $n$, and forming a recurrent sequence) satisfy an integral linear recurrence.
I haven't considered yet the case of sequences of the form :
$$\sum_{i=1}^{\ell} P_i(n) x_i^n$$
where $x_i$ are distinct, and $P_i(n)$ are polynomial in $n$, not $0$.
Any feedback would be appreciated!
$\bf{Added:}$
Some definitions for clarity:
A linear recurrence sequence is a sequence satisfying a linear recurrence, that is, there exists a $d$, and $c_k$, $1\le k \le d$ (apriori $c_k \in \mathbb{C}$) such that for all $n \ge d$ we have
$$a_n = \sum_{k=1}^d c_k a_{n-k} \ \ \ (*)$$
that is, $a_n$ dependes linearly on the previous $d$ terms, for all $n\ge d$.
We say that the recurrence has integers coefficients if all $c_0$, $\ldots$, $c_d$ are integers.
$\bf{Added:}$ A related question ( a particular case).
$\bf{Added:}$ Ewan Delanoy showed that $a_n$ satisfies a recurrence with rational coefficients. A big step forward.
I have a proof along this lines: say $K\subset L$ field and and $c_1$, $\ldots$ $c_d$ in $L$. Consider $V$ the $K$ span of $1, c_j$ in $L$ ( a finite dimensional $K$ vector subspace). There exists a $K$-linear projection $\pi$ from $V$ to $K$ with $\pi(1) = 1$ ( basic linear algebra).
Now consider an equality
$$a = \sum c_j a_j$$
with $a_j$, $a \in K$. Apply the map $\pi$ to it and get
$$a= \sum \pi(c_j) a_j$$
Moral: any linear dependence over $L$ of elements in $K$ produces a linear dependence over $K$ of said elements.