  Spectra of stochastic gravitational waves (GW) generated in cosmological
first-order phase transitions are computed within strongly correlated theories
with a dual holographic description. The theories are mostly used as models of
dark sectors. In particular, we consider the so-called Witten-Sakai-Sugimoto
model, a $SU(N)$ gauge theory coupled to different matter fields in both the
fundamental and the adjoint representations. The model has a well-known
top-down holographic dual description which allows us to perform reliable
calculations in the strongly coupled regime. We consider the GW spectra from
bubble collisions and sound waves arising from two different kinds of
first-order phase transitions: a confinement/deconfinement one and a chiral
symmetry breaking/restoration one. Depending on the model parameters, we find
that the GW spectra may fall within the sensibility region of ground-based and
space-based interferometers, as well as of Pulsar Timing Arrays. In the latter
case, the signal could be compatible with the recent potential observation by
NANOGrav. When the two phase transitions happen at different critical
temperatures, characteristic spectra with double frequency peaks show up.
Moreover, in this case we explicitly show how to correct the redshift factors
appearing in the formulae for the GW power spectra to account for the fact that
adiabatic expansion from the first transition to the present times cannot be
assumed anymore.
