I will present our recent results analyzing HSC weak lensing data using non-Gaussian statistics, where 30% improvement in cosmological parameter constraints was achieved. I will also summarize our ongoing effort in preparing for future non-Gaussian statistics analysis with Rubin LSST.