Estimating fatigue damage of marine structures using spectral methods under broadband random loadings is a complex and challenging task. In the presence of non-Gaussian loading, the spectral method typically employs the Hermite polynomial model (HPM) to address the non-Gaussianity. However, the applicability of HPM is limited in certain non-Gaussian cases due to the requirement of monotonicity. This paper adopts the Johnson transformation model (JTM) instead of HPM and integrates it with the non-Gaussian TB method (Ding and Chen, 2015) to overcome the constraints in broadband non-Gaussian fatigue estimation. To verify the performance of the new approach, this paper analyses the transformation accuracy of the two models. Additionally, taking the time-domain fatigue damage obtained by the rainflow counting method and the linear damage accumulation rule as a benchmark, this paper compares the fatigue damage estimates obtained from the spectral methods. The comparative assessment encompasses various spectra, S-N curves, and non-Gaussian characteristics. The results indicate the performance of JTM within the non-Gaussian TB method is comparable to HPM in most cases and even surpasses it in hardening non-Gaussian cases. Given the broader applicability of JTM across ranges of skewness and kurtosis, it emerges a superior option for the non-Gaussian TB method in estimating broadband non-Gaussian fatigue damage. Furthermore, the limitations of the proposed approach are discussed.