We propose a new method of calculating entanglement entropy of a many-body interacting Bosonic system (open or closed) in a field theoretic approach without replica methods. The Wigner function and Renyi entropy of a Bosonic system undergoing arbitrary non-equilibrium dynamics can be obtained from its Wigner characteristic function, which we identify with the Schwinger Keldysh partition function in presence of quantum sources turned on at the time of measurement. For non-interacting many body systems, starting from arbitrary density matrices, we provide exact analytic formulae for Wigner function and entanglement entropy in terms of the single particle Green's functions. For interacting systems, we relate the Wigner characteristic to the connected multi-particle correlators of the system. We use this formalism to study the evolution of an open quantum system from a Fock state with negative Wigner function and zero entropy, to a thermal state with positive Wigner function and finite entropy. The evolution of the Renyi entropy is non-monotonic in time for both Markovian and non-Markovian dynamics. The entropy is also found to be anti-correlated with negativity of the Wigner function of a
2 -mode open quantum system.