Siegel modular forms as sections of line bundles over the period domain
Dear Samuel,
I don't know if the holomorphic line bundles on the quotients you ask about are classified, but one point of view that is adopted when thinking of (Siegel, or other) modular forms as sections of bundles is to consider $G$-equivariant holomorphic bundles on the symmetric space, so in your case, this would be $Sp_{2g}(\mathbb R)$-equivariant bundles on Siegel upper half-space. Such equivariant bundles can be used to induce a holomorphic bundle on the quotient of the upper half-space by any discrete group (in a way that is compatible with pull-back when one discrete group is contained in another), and so these are the bundles that are normally the most interesting in relation to automorphic forms.
The recipe for making these bundles (which works in some generality, say in an Shimura variety context; but I will stick to your Siegel situation) is as follows:
We first note that, just as the usual complex upper half-plane embeds into $\mathbb P^1(\mathbb C)$, the Siegel upper half-space embeds as an open subset of the complex points of a certain partial flag variety of $Sp_{2g},$ namely $Sp_{2g}/P$, where $P$ is the Siegel parabolic of $Sp_{2g}$ (the stablizer of a maximal isotropic subspace in the $2g$-dimensional symplectic space on which $Sp_{2g}$ acts). The Levi of $P$ is equal to $GL_g$, and so the irreducible $Sp_{2g}$-equivariant vector bundles on $Sp_{2g}/P$ correspond to irreducible reps. of $GL_g$. (If $V$ is a rep. of $GL_g$, we think of it as a rep. of $P$ through the quotient map $P \to GL_g$, and then make a bundle $\mathcal V$ on $Sp_{2g}/P$ via the formula $\mathcal V = (Sp_{2g} \times V)/P.$ In particular, the rank of $\mathcal V$ equals the dimension of $V$.
Restricting these bundles $\mathcal V$ to the Siegel upper half-space, one obtains a family of $Sp_{2g}(\mathbb R)$-equivariant holomorphic bundles on the Siegel upper half-space. As I already mentioned, because of their equivariance, they induce holomorphic bundles on the quotient by any discrete subgroup of the symplectic group.
The line bundles whose sections are classical Siegel modular forms correspond to the one-dimensional representations of $GL_g$ given by taking powers of the determinant.
Other (higher than one-dimensional) rep's of $GL_g$ give rise to bundles whose sections are vector-valued Siegel modular forms.
I found this kind of material a bit tricky to find in the literature, and pieced it together by reading Deligne's papers on Shimura varieties, as well as the early parts of various articles of Michael Harris and Jim Milne about automorphic vector bundles on Shimura varieties.
I wouldn't be surprised if there is a more classical literature that is complimentary to the Shimura-variety literature just mentioned, and which is perhaps more accessible. Unfortunately I don't know it!
Modulo torsion, the Picard group of the quotient of the Siegel upper half plane by a finite-index subgroup of $\text{Sp}_{2g}(\mathbb{Z})$ is just $\mathbb{Z}$. This result should probably be attributed to Borel.
For a calculation of the torsion and explicit line bundles corresponding to the various pieces (plus references), see my paper "The Picard group of the moduli space of curves with level structures", available here. The results you want are discussed starting on page 3.
In the case $g=1$ modular forms are sections of a line bundle on a (geometrically meaningful) compactification of the stack obtained by quotienting the upper half plane by $SL_2(\mathbb Z)$ (usually denoted $\overline{M}_{1,1}$), whose Picard group is indeed isomorphic to $\mathbb Z$. By contrast, the Picard group before compactifying is torsion of order 12.
Since such compactifications exists AFAIK for every $g$, I suspect that this is what you have to consider. Bundles on the stack quotient of $X$ by $G$ are the same as $G$-equivariant bundles on $X$, so this answer has nontrivial overlap with Emerton's.