Ethnopharmacological relevanceHaizao Yuhu decoction (HYD) is a classic Chinese herbal formula described in the surgical monographs of the Ming Dynasty “Waikezhengzong.” It has been widely used to treat goiter for approximately 500 years and found to be particularly effective. HYD contains glycyrrhiza and sargassum. This pair of herbs belongs to “18 incompatible medicaments” of traditional Chinese medicine theory. Although these two herbs are opposite, our preliminary study proved that they have superior effect when added into HYD at 2 times the dose of Chinese Pharmacopoeia. However, the species of glycyrrhiza in HYD that are the most effective have not been recorded in ancient Chinese medical texts. According to the Chinese Pharmacopoeia, glycyrrhiza is divided into the following three species: Glycyrrhiza uralensis Fish., G. glabra L., and G. inflata Bat. The effect of HYD containing different species of glycyrrhiza and their mechanisms remain to be further explored. Aim of the studyTo investigate the effect of HYD containing three species of glycyrrhiza on goiter, and to elucidate the molecular mechanism using network pharmacology combined with RNA sequencing (RNA-seq). Materials and methodsA rat model of goiter was established by 14 days of intragastric gavage of propylthiouracil (PTU), and the rats were treated for 4 weeks with HYD containing three different species of glycyrrhiza. The body weight and rectal temperature of rats were tested weekly. At the end of the experiment, the serum and thyroid tissues of rats were collected. The effect of the three HYDs was assessed based on general observations (including body weight, rectal temperature, and living status of rats), absolute/relative thyroid weight, thyroid function (including triiodothyronine, thyroxine, free triiodothyronine, free thyroxine, and thyroid-stimulating hormone levels), and thyroid tissue pathology. Next, we explored their pharmacological mechanisms using network pharmacology combined with RNA-seq and validated key targets using real-time quantitative reverse-transcription polymerase chain reaction (RT-qPCR), western blotting (WB), and immunofluorescence (IF) assays. ResultsThe three HYDs reduced the absolute/relative weights of thyroid tissues and improved the pathological structure, thyroid function, and general findings of rats with goiter. Overall, the effect of HYD-G. uralensis Fish. (HYD-U) was better. Results from network pharmacology and RNA-seq jointly suggested that both the pathogenesis of goiter and the mechanism of action of HYD for goiter were related to the phosphatidylinositol 3-kinase-protein kinase B (PI3K-Akt) pathway. We validated the key targets in the pathway, namely, vascular endothelial growth factor (VEGF) A, VEGF receptor 2, phosphoinositide-3-kinase regulatory subunit 1 (PIK3R1) and its encoded protein PI3K (p85), AKT serine/threonine kinase 1 (AKT1), phospho-AKT and cyclin D1 using RT-qPCR, WB, and IF assays. The PI3K-Akt pathway was hyperactivated in rats with PTU-induced goiter, whereas the three HYDs could inhibit the pathway. ConclusionThis study confirmed the definite effect of the three HYDs in the treatment of goiter, and HYD-U was found to be more effective. The three HYDs inhibited angiogenesis and cell proliferation in goiter tissue by inhibiting the PI3K-Akt signaling pathway.