I was pointed to this very nice paper On the CNOT-cost of TOFFOLI gates
. The issue appears to be much more complicated than I thought. Apparently there are theses yet to be written on the decomposition of Toffoli into $CNOT$s.
First, due to identity $HXH=Z$ decompositions into $CNOT$s and $CZ$s can be translated into each others. Similarly, decomposing $n$-qubit Toffoli can be reduced to decomposing $n$-qubit controlled-$Z$ gate, which might be technically simpler as it is symmetric and diagonal.
The main result derived in the paper is that the cost of a $Z$ gate with $n-1$ controls satisfies $|C^{(n-1)}Z|_{CZ}\ge 2n$, which implies the same bound for decomposition of an $n$-qubit Toffoli into $CNOT$s: $|C^{(n-1)}X|_{CX}\ge 2n$. For $n=3$ this gives 6 $CNOT$ gates and there is a known decomposition with this amount of gates (see OP), so for $n=3$ this is also the optimal result. Already for $n=4$ the bound appears quite weak saying that one needs at least $8$ gates, but the best known decomposition (which is probably optimal) requires 14 gates.
The paper also proves that for $n=3$ adding any amount of ancilla qubits can not lower the cost. For $n\ge4$ ancilla probably can help, but coming up with explicit examples seems to be surprisingly difficult.
I do not know if there were improvements to the results of the paper, I would be interested to learn that.