Outflows and rotation are two ubiquitous kinematic features in the gas kinematics of galaxies. Here we introduce a semi-analytic model to quantify how rotating outflows impact the morphology of the Lyman-$\alpha$ emission line. The model is contrasted against Monte Carlo radiative transfer simulations of outflowing gas with additional solid body rotation. We explore a range of neutral Hydrogen optical depth of $10^5 \leq \tau_{\mathrm{H}} \leq10^7$, rotational velocity $0 \leq v_{\mathrm{rot}}/\mathrm{km\ s}^{-1} \leq 100$ and outflow velocity $0\leq v_{\mathrm{out}}/\mathrm{km\ s}^{-1} \leq 50$. We find three consequences of rotation. First, it introduces a dependency with viewing angle; second it broadens the line and third it increases the flux at the line's center. For all viewing angles, the semi-analytic model reproduces the radiative transfer results for the line width and flux change at the line's center within a $7\%$ and $50\%$ precision for an optical depth of $\tau_{\mathrm{H}}=10^5$, respectively, and within $2\%$ and $1\%$ for an optical depth of $\tau_{\mathrm{H}}=10^7$. Using this model we also show that the peaks of integrated spectra taken from opposite sides of an edge-on rotating gas distribution should have a separation of $\frac{1}{2}v_{\mathrm{rot}}$. The semi-analytic model presented here is a convenient tool to introduce rotational kinematics as a post-processing step of idealized Monte Carlo simulations; it provides a framework to interpret \lya spectra in systems where rotation is expected or directly measured through kinematic maps.

Star-forming Compact Dwarf Galaxies (CDGs) resemble the expected pristine conditions of the first galaxies in the Universe and are the best systems to test models on primordial galaxy formation and evolution. Here we report on one of such CDGs, Tololo 1214-277, which presents a broad, single peaked, highly symmetric Ly$\alpha$ emission line that had evaded theoretical interpretation so far. In this paper we reproduce for the first time these line features with two different physically motivated kinematic models: an interstellar medium composed by outflowing clumps with random motions and an homogeneous gaseous sphere undergoing solid body rotation. The multiphase model requires a clump velocity dispersion of $54.3\pm 0.6$ km s$^{-1}$ with outflows of $54.3\pm 5.1$ km s$^{-1}$, while the bulk rotation velocity is constrained to be $348^{+75}_{-48}$ km s$^{-1}$. We argue that the results from the multiphase model provide a correct interpretation of the data. In that case the clump velocity dispersion implies a dynamical mass of $2\times 10^{9}$ M$_{\odot}$, ten times its baryonic mass. If future kinematic maps of Tololo 1214-277 confirm the velocities suggested by the multiphase model, it would provide additional support to expect such kinematic state in primordial galaxies, opening the opportunity to use the models and methods presented in this paper to constrain the physics of star formation and feedback in the early generation of Ly-$\alpha$ emitting galaxies.