I'll try to provide a proof of the general case $[y_0,\ldots,y_n;f]$. This proof has many similarities to that of timur, but I'll formulate things more in terms of limits rather than using $\epsilon$s and $\delta$s. My other answer was wrong, but I leave it since it had some useful bits, one of which I'll redo here.
Hope I didn't overlook anything grave this time. However, I'm sure I've taken a few technical short-cuts which could need filling in: the lemma provided at the end is an attempt at that.
Notation: Let $R^{n+1}_k\subset\mathbb{R}^{n+1}$ consist of the points $y=(y_0,\ldots,y_n)$ for which at most $k$ coefficients have the same value: i.e. $R^{n+1}_1$ are those where all $y_i$ are distinct, while $R^{n+1}_{n+1}=\mathbb{R}^{n+1}$. We then have $[\cdot;f]:R^{n+1}_1\rightarrow\mathbb{R}$ defined as $ [y_0,\ldots,y_n;f]=\frac{[y_1,\ldots,y_n;f]-[y_0,\ldots,y_{n-1};f]}{y_n-y_0} \tag{1} $ for $y\in R^{n+1}_{1}$: generalised from the cases $n=0,1,2$ stated in the original problem.
Definition: We say that a function $f$ is $\mathcal{F}^n$ if for any $y,y'\in R^{n+1}_1$, $\lim_{y,y'\rightarrow x^{[n+1]}}\Big|[y;f]-[y';f]\Big|=0\tag{2}$ where, for $x\in\mathbb{R}$, the point $x^{[n+1]}=(x,\ldots,x)\in\mathbb{R}^{n+1}$. This is just the equivalent of the $\epsilon$-$\delta$ formulation in terms of limits.
(i) Assume $f$ is $\mathcal{F}^n$. If we take any sequence $y^{(k)}\in R^{n+1}_1$ which converges to $x^{[n+1]}$, then $[y^{(k)};f]$ becomes a Cauchy sequence: i.e. for $p,q\ge k$, $\big|[y^{(p)};f]-[y^{(q)};f]\big|\rightarrow0$ as $n\rightarrow\infty$. Hence, $[f^{(k)};f]$ must converge. Let's denote the limit $[x^{[n+1]};f]=g_n(x)$.
(ii) Let's prove by induction that $[y_0,\ldots,y_n;f]$ is symmetric in the $y_i$. Obviously, it holds for $n=0,1$. Assume it holds for $n. Setting $n=N\ge2$, let $y=(u,v,A,w)\in R^{n+1}_1$ where $A=(a_1,\ldots,a_{n-2})$ (possibly empty). We then have $ \begin{split} [u,v,A,w;f]=&\frac{[u,v,A;f]-[v,A,w;f]}{u-w}\\ =&\frac{1}{u-w}\Big([u,A,v;f]-[v,A,w;f]\Big)\\ =&\frac{1}{u-w}\left(\frac{[u,A;f]-[A,v;f]}{u-v}-\frac{[v,A;f]-[A,w;f]}{v-w}\right)\\ =&\frac{[u,A;f]}{(u-v)(u-w)}+\frac{[v,A;f]}{(v-u)(v-w)}+\frac{[w,A;f]}{(w-u)(w-v)} \end{split} $ which is symmetric in $u,v,w$. From the first step, we also see that it must be symmetric in permutations of $v,A$. These symmetries generate the whole permutation group, so symmetry also holds for $n=N$.
(iii) Again, we use induction, assuming the following holds for $n: A function $f$ is $\mathcal{F}^n$ if and only if it is $C^n$, i.e. $\mathcal{F}^n=C^n$: in particular, it is also $\mathcal{F}^k$ for all $k. The definition of $[y;f]$ for $y\in R^{n+1}_1$ may be extended continuously to all $y\in\mathbb{R}^{n+1}$. Moreover, $g_n(x)=f^{(n)}(x)/n!$ where $f^{(n)}$ is the $n$th derivative.
For $n=0,1$, the asumption holds true. So we let $n=N\ge2$.
In order to utilise the induction hypothesis, we first need to show that when $f$ is $\mathcal{F}^n$, it is also $\mathcal{F}^{n-1}$. To do this, we rewrite equation (1) $(u-v)[u,A,v;f]=[u,A;f]-[A,v;f]$ for $A=(a_1,\ldots,a_{n-1})$. For $u=(u_1,\ldots,u_n)$ and $v=(v_1,\ldots,v_n)$, both in $R^n_1$, we then get $ \begin{split} (u_1-v_1)[u_1,\ldots,u_n,v_1;f]&+(u_2-v_2)[u_2,\ldots,u_n,v_1,v_2;f]\\ &+\cdots+(u_n-v_n)[u_n,v_1,\ldots,v_n;f]\\ =&[u_1,\ldots,u_n;f]-[v_1,\ldots,v_n;f]\\ \end{split}\tag{3} $ where $u,v\rightarrow x^{[n]}$ leads to $\big|[u;f]-[v;f]\big|\rightarrow0$, i.e. the definition of $f$ being $\mathcal{F}^{n-1}$.
For any $y\in R^{n+1}_n$, i.e. not all $y_i$ are equal, we can rearrange the $y_i$ to make $[y;f]=[u,A,v;f]$ where $u\not=v$ and use (1) to get $[y;f]$ expressed in terms of $[u,A;f]$ and $[A,v;f]$, both of which are defined on the whole $\mathbb{R}^n$. This provides us with an extension of $[y;f]$ to all $y\in R^{n+1}_n$.
Finally, in equation (3), we let $u=x^{[n]}$, $v=x'^{[n]}$: by continuity, we can do that, or alternatively we can take limits $u\rightarrow x^{[n]}$ and $v\rightarrow x'^{[n]}$ (and we take these limits before the later limit). Dividing by $x-x'$ on both sides then gives $ \frac{g_{n-1}(x)-g_{n-1}(x')}{x-x'} =[x,\ldots,x,x';f]+\cdots+[x,x',\ldots,x';f] $ where $x'\rightarrow x$ gives $g_{n-1}'(x)=n\cdot[g_n(x)]$. Hence, $g_n(x)=f^{(n)}/n!$.
The only remaining thing to point out is that the final extension of $[y;f]$ from $y\in R^{n+1}_n$ to all of $\mathbb{R}^{n+1}$ is continuous. This basically follows from the definition, but we have to use the following lemma (which I have implicitly used further up as well), to ensure the limit in (2) holds true even if the sequence even if $y\in R^{n+1}_n$ rather than $R^{n+1}_1$.
Lemma: Given $x,v_i,u^{(i)}_j\in\mathbb{R}$ so that $v_i\rightarrow x$ and $u^{(i)}_{j}\rightarrow v_i$, then we can pick $j_i$ so that $u^{(i)}_{j_i}\rightarrow x$.
One use of the lemma is that if $y^{(i)}\in R^{n+1}_n$, $z^{(i,j)}\in R^{n+1}_1$, we can set $v_i=[y^{(i)};f]$, $u^{(i)}_j=[z^{(i,j)};f]$, and replace the limit of $[y^{(i)};f]$, which is taken in $R^{n+1}_n$, with a limit taken in $R^{n+1}_1$.