I wonder if the following upper bound holds:
$${}_2F_1[-m, -m; -(m+l); z]\leq 1,\tag{1}$$
where ${}_2F_1[a,b;c;z]$ is the Gauss hypergeometric function, $m,l=0,1,2,\ldots$, and $0<z<1$. Clearly, this holds for $l=0$ since (1) then reduces to the binomial identity in the answer to a related question. Thus, induction seems to be a promising path to proving this upper bound. However, direct application of the expansion:
$$\begin{align}{}_2F_1[-m, -m; -(m+l); z] &= {}_2F_1[-m, -m; -(m+l-1); z] \\ &\quad + \dfrac{m^2 z}{(m+l)(m+l-1)}{}_2F_1[-(m-1), -(m-1); -(m-1+l-1); z]\end{align}$$ derived for the aforementioned solution does not seem to work here. Numerical experiments seem to confirm this bound, though perhaps there is a counterexample. I tried induction on $m$ instead of $l$ to no avail. Any ideas?