We study analytic smooth solutions of a general, strongly parabolic semilinear Cauchy problem of 2m-th order in R-N x (0, T) with analytic coefficients (in space and time variables) and analytic initial data (in space variables). They are expressed in terms of holomorphic continuation of global (weak) solutions to the system valued in a suitable Besov interpolation space of B-s;p,B-p-type at every time moment t is an element of [0, T]. Given 0 < T' < T <= infinity , it is proved that any B-s;p,B-p-type solution u : R-N x (0, T) -> C-M with analytic initial data possesses a bounded holomorphic continuation u(x + iy, sigma + i tau) into a complex domain in C-N x C defined by (x, sigma) is an element of R-N x (T', T), vertical bar y vertical bar < A' and vertical bar tau vertical bar < B', where A', B' > 0 are constants depending upon T'. The proof uses the extension of a weak solution to a B-s;p,B-p-type solution in a domain in C-N x C, such that this extension satisfies the Cauchy-Riemann equations. The holomorphic extension is obtained with a help from holomorphic semigroups and maximal regularity theory for parabolic problems in Besov interpolation spaces of B-s;p,B-p-type imbedded (densely and continuously) into an L-p-type Lebesgue space. Applications include risk models for European options in Mathematical Finance.