Well yes, I chose a frame of the tangent bundle coming from a coordinate chart, which implies that the Lie derivative vanishes, since the conventions I am familiar with reserve the symbol \Gamma^a_{bc} only for these types of bases. You can simply replace the special frame \left\{\frac{\partial}{\partial x^b} \right\} I chose in my previous post with an arbitrary local frame \{e_i \} of the tangent bundle and the identical argument works but now the symbols \Gamma^a_{bc} are usually not called Christoffel symbols and I've always seen them denoted by different symbols.
In any case I believe you are asking about the definition of the torsion tensor not proving anything about it. In fact the derivation you posted is already a complete general definition of the torsion tensor. Defining a tensor in every coordinate coordinate chart of some atlas for your manifold completely defines the tensor and you don't need to even consider the existence of nonholonomic coordinate bases when giving the definition. So we have a tensor defined in local charts by T^i_{jk}=\Gamma^i_{jk}-\Gamma^i_{kj} and your question is simply asking what form this tensor takes when you write it in terms of some basis which did not come from a chart.
From this definition, if we have some (possibly nonholonomic) local frame \{ e_1,\cdots, e_n\}, you can derive the equation T(e_b,e_c)=\left(\Gamma^a_{bc}-\Gamma^a_{cb}-\gamma^a_{bc}\right)e_a by choosing a chart (U,x^i) which overlaps the domain of our local frame around any given point and then expanding the basis elements as e_b=\sum_{\alpha=1}^n K_b^\alpha \frac{\partial}{\partial x^\alpha}. I don't have time to type up the complete computation since it is a bit tedious but there are no tricks, just direct computation.
Just to get you started, by definition T^a_{bc}e_a<br />
=T(e_b,e_c)<br />
=T\left( K_b^j \frac{\partial}{\partial x^j}, K_c^k \frac{\partial}{\partial x^k}\right)<br />
=K^j_b K_c^kT\left( \frac{\partial}{\partial x^j},\frac{\partial}{\partial x^k}\right)<br />
=K^j_b K_c^k\left( \Gamma^i_{jk}-\Gamma^i_{kj} \right) \frac{\partial}{\partial x^i}.<br />
Now do a similar thing using the inverse of the matrix K_b^\alpha and the definition of the \Gamma's in terms of covariant derivatives to write the terms of the form \Gamma^k_{ij} in terms of \Gamma^a_{bc}. Substituting this back into the equation for T(e_b,e_c) will result in the formula you are trying to derive. The Lie derivative of vector fields appears in the computation since you will need to use the Leibniz rule to compute the covariant derivative at some point which will provide the extra terms.