Gelsemium elegans is a traditional Chinese herb of medicinal importance, with indole terpene alkaloids as its main active components. To study the expression of the most suitable housekeeping reference genes in G. elegans, the root bark, stem segments, leaves and inflorescences of four different parts of G. elegans were used as materials in this study. The expression stability of 10 candidate housekeeping reference genes (18S, GAPDH, Actin, TUA, TUB, SAND, EF-1α, UBC, UBQ, and cdc25) was assessed through real-time fluorescence quantitative PCR, GeNorm, NormFinder, BestKeeper, ΔCT, and RefFinder. The results showed that EF-1α was stably expressed in all four parts of G. elegans and was the most suitable housekeeping gene. Based on the coexpression pattern of genome, full-length transcriptome and metabolome, the key candidate targets of 18 related genes (AS, AnPRT, PRAI, IGPS, TSA, TSB, TDC, GES, G8H, 8-HGO, IS, 7-DLS, 7-DLGT, 7-DLH, LAMT, SLS, STR, and SGD) involved in the Gelsemium alkaloid biosynthesis were obtained. The expression of 18 related enzyme genes were analyzed by qRT-PCR using the housekeeping gene EF-1α as a reference. The results showed that these genes’ expression and gelsenicine content trends were correlated and were likely to be involved in the biosynthesis of the Gelsemium alkaloid, gelsenicine.