  Working within the deep-MOND limit (DML), I describe spherical,
self-gravitating systems governed by a polytropic equation of state,
$P=\mathcal{K}\rho^\gamma$. As self-consistent structures, such systems can
serve as heuristic models for DML, astronomical systems, such as dwarf
spheroidal galaxies, low-surface-density elliptical galaxies and star clusters,
and diffuse galaxy groups. They can also serve as testing ground for various
theoretical MOND inferences. In dimensionless form, the equation satisfied by
the radial density profile $\zeta(y)$ is (for $\gamma\not=1$) $[\int_0^y \zeta
\bar y^2 d\bar y]^{1/2}=-yd(\zeta^{\gamma-1})/dy$. Or,
$\theta^n(y)=y^{-2}[(y\theta')^2]'$, where $\theta=\zeta^{\gamma-1}$, and
$n\equiv (\gamma-1)^{-1}$. I discuss properties of the solutions, contrasting
them with those of their Newtonian analogues -- the Lane-Emden polytropes. Due
to the stronger MOND gravity, all DML polytropes have a finite mass, and for
$n<\infty$ ($\gamma>1$) all have a finite radius. (Lane-Emden spheres have a
finite mass only for $n\le 5$.) I use the DML polytropes to study DML scaling
relations. For example, they satisfy a universal relation (for all
$\mathcal{K}$ and $\gamma$) between the total mass, $M$, and the mass-average
velocity dispersion $\sigma$: $MGa_0=(9/4)\sigma^4$. However, the relation
between $M$ and other measures of the velocity dispersion, such as the central,
projected one, $\bar\sigma$, does depend on $n$ (but not $\mathcal{K}$),
defining a `fundamental surface' in the $[M,~\bar\sigma,~n]$ space. I also
describe the generalization to anisotropic polytropes, which also all have a
finite radius (for $\gamma>1$), and all satisfy the above universal $M-\sigma$
relation. This more extended class of models exhibits the
central-surface-densities relation: a tight relation between the baryonic and
the dynamical central surface densities predicted by MOND.
