For flavor $SU(2)$ (Isospin) we have two $\frac{1}{2}^+$ baryons, the nucleons. For flavor $SU(3)$ we have the eight baryons in the octet. In a world with $N$ light quarks we would see a baryon multiplet of dimension $\frac{N}{3}(N^2-1)$.
Such a theory would see the chiral symmetry breaking $SU(N)\times SU(N) \to SU(N)$ creating $N^2-1$ Goldstone bosons $\phi^a, a=1,\dots,N^2-1$. These mesons are usually parametrized in the $N\times N$ matrix $U=\exp\left( i T^a \phi^a \right)$, where $T^a$ are the generators of $SU(N)$. These fields than transform under $(L,R)\in SU(N)\times SU(N)$ as $U \mapsto R U L^\dagger$.
So far so good (please correct me if I already made any mistakes, I think this is correct). My question now is:
How can we include the $\frac{N}{3}(N^2-1)$ baryons in a Lagrangian?
We need to find a parametrization for these baryon fields and we have to find out how they transform under $SU(N) \times SU(N)$.
It feels like we got a bit lucky in reality, since for $N=2$: $\frac{N}{3}(N^2-1)=N=2$, and for $N=3$: $\frac{N}{3}(N^2-1)=N^2-1=8$ is the number of generators of $SU(3)$. So in these cases we can use the isospin doublet $N = \begin{pmatrix} p\\n\end{pmatrix}$ and for octet baryons we can use $B=\sum_{a=1}^{8} \frac{B^{a} \lambda^{a}}{\sqrt{2}}=\left[\begin{array}{ccc}{\Sigma^{0} / \sqrt{2}+\Lambda / \sqrt{6}} & {\Sigma^{+}} & {p} \\ {\Sigma^{-}} & {-\Sigma^{0} / \sqrt{2}+\Lambda / \sqrt{6}} & {n} \\ {\Xi^{-}} & {\Xi^{0}} & {-\sqrt{\frac{2}{3}} \Lambda}\end{array}\right].$
Is there literature for chiral perturbation theory for $N$ light flavors? How would one go about including baryons?