Let $S$ be a countable state space and $(X_n)$ an $S$-valued discrete Markov chain. Then $$ \mathbb P [X_{n+1} = x_{n+1} |X_n = x_n, \ldots, X_0=x_0] = \mathbb P [X_{n+1} = x_{n+1} |X_n = x_n] $$ for all $x_0, \ldots, x_{n+1} \in S$. This is called the Markov property. Inspired by this thread, I have proved below generalization, i.e.,
Theorem $$ \begin{align} & \mathbb P [X_{n+1} \in S_{n+1} |X_{i_M} = x_{i_M}, X_{i_{M-1}} \in S_{i_{M-1}} \ldots, X_{i_m} \in S_{i_{m}}] \\ = &\mathbb P [X_{n+1} \in S_{n+1} |X_{i_M} = x_{i_M}] \end{align} $$ for all $i_m \le i_{m+1} \le \cdots \le i_M \le n$ such that $S_{n+1} \subset S, x_{i_M} \in S, S_{i_{M-1}} \subset S, \ldots, S_{i_m} \subset S$.
Can we generalize further by relacing $X_{i_M} = x_{i_M}$ with $X_{i_M} \in S_{i_{M}}$? If not, is there some intuition why it does not hold?
Thank you so much for your elaboration!
Proof Because $\mathbb P [\bigcup_n A_n | B] = \sum_n \mathbb P [A_n | B]$, it suffices to consider the case $S_{n+1}$ is a singleton, i.e., $S_{n+1} = \{x_{n+1}\}$. Let $\pmb S := S_{i_{M-1}} \times \cdots \times S_{i_{m}}$ and $\pmb Y := (X_{i_{M-1}}, \ldots, X_{i_m})$. Then $$ \begin{align*} & \mathbb P [X_{n+1} = x_{n+1} |X_{i_M} = x_{i_M}, X_{i_{M-1}} \in S_{i_{M-1}} \ldots, X_{i_m} \in S_{i_{m}}] \\ = &\mathbb P [X_{n+1} = x_{n+1} |X_{i_M} = x_{i_M}, \pmb Y \in \pmb S] \\ = & \frac{\mathbb P [X_{n+1} = x_{n+1} ,X_{i_M} = x_{i_M}, \pmb Y \in \pmb S]}{\mathbb P [X_{i_M} = x_{i_M}, \pmb Y \in \pmb S]} \\ = & \frac{\sum_{\pmb y \in \pmb S } \mathbb P [X_{n+1} = x_{n+1} ,X_{i_M} = x_{i_M}, \pmb Y = \pmb y]}{\sum_{\pmb y \in \pmb S }\mathbb P [X_{i_M} = x_{i_M}, \pmb Y = \pmb y]} \\ = & \frac{\sum_{\pmb y \in \pmb S } \mathbb P [X_{n+1} = x_{n+1} |X_{i_M} = x_{i_M}, \pmb Y = \pmb y] \mathbb P [ X_{i_M} = x_{i_M}, \pmb Y = \pmb y]}{\sum_{\pmb y \in \pmb S }\mathbb P [X_{i_M} = x_{i_M}, \pmb Y = \pmb y]} \\ = & \frac{\sum_{\pmb y \in \pmb S } \mathbb P [X_{n+1} = x_{n+1} |X_{i_M} = x_{i_M}] \mathbb P [ X_{i_M} = x_{i_M}, \pmb Y = \pmb y]}{\sum_{\pmb y \in \pmb S }\mathbb P [X_{i_M} = x_{i_M}, \pmb Y = \pmb y]} \\ = & \mathbb P [X_{n+1} = x_{n+1} |X_{i_M} = x_{i_M}], \end{align*} $$ where the fifth inequality follows from the following Lemma, i.e.,
Lemma $$ \mathbb P [X_{n+1} = x_{n+1} |X_{i_M} = x_{i_M}, \ldots, X_{i_m}=x_{i_m}] = \mathbb P [X_{n+1} = x_{n+1} |X_{i_M} = x_{i_M}] $$ for all $i_m \le i_{m+1} \le \cdots \le i_M \le n$ and $x_{i_m}, \ldots, x_{i_M}, x_{n+1} \in S$.