Although the existing answer is complete and there is nothing more to be said, I believe a new answer can at best clarify notation. So I will re-write the expression that the OP has to make the appearance of the boundary conditions manifest:
Here I will call $q_0=q_i$, $q_N=q_f$
$$\int Dq(t) \exp\Big[\frac{i}{\hbar}\int_0^Tdt(\frac{1}{2}m\dot{q}^2-V(q)\Big]$$$$=\lim_{N\to\infty}\Big(\frac{-2\pi imN}{T}\Big)^{N/2}\prod_{k=1}^{N-1}\int_{-\infty}^{+\infty}dq_k.\exp\Big[\frac{i}{\hbar}\sum_{k=1}^N \epsilon\left(\frac{1}{2}m\left(\frac{q_k-q_{k-1}}{\epsilon}\right)^2-V(q_k)\right)\Big]$$
Here $\epsilon=\frac{T}{N}$.
In particular, the boundary conditions always influences the time derivative term at the ends of the path integral and the potential term at one boundary(which one is a matter of convention, and is deemed unimportant in the large $N$ limit):
$$=\lim_{N\to\infty}\Big(\frac{-2\pi imN}{T}\Big)^{N/2}\prod_{k=1}^{N-1}\int_{-\infty}^{+\infty}dq_k.\exp\Big[\frac{i}{\hbar}\sum_{k=2}^{N-1} \epsilon\left(\frac{1}{2}m\left(\frac{q_k-q_{k-1}}{\epsilon}\right)^2-V(q_k)\right)+\frac{i}{\hbar}\epsilon\left(\frac{1}{2}m\left(\frac{q_1-q_{0}}{\epsilon}\right)^2-V(q_1)+\frac{1}{2}m\left(\frac{q_N-q_{N-1}}{\epsilon}\right)^2-V(q_N)\right)\Big]$$
Where in the last line, I have separated the $k=1,N$ terms to explicitly see the boundary conditions in the path integral.
The wisdom in Qmechanic's answer indicates, correctly, that the continuum representation of the path integral, is entirely formal and has little tangible meaning. This limiting expression is the closest one can get to explicitly writing down a path integral. And in this expression, the boundary terms are apparent.