Phlomoides rotata is a traditional Chinese herbal medicine that grows in the Qinghai–Tibet Plateau region at a 3100–5000 m altitude. Iridoid compounds are the main active compounds of the P. rotata used as medical ingredients and display anti-inflammatory, analgesic, and hepatoprotective properties. To better understand the biological mechanisms of iridoid compounds in this species, we performed a comprehensive analysis of the transcriptome and metabolome of P. rotata leaves from four different regions (3540–4270 m). Global metabolome profiling detected 575 metabolites, and 455 differentially accumulated metabolites (DAMs) were detected in P. rotata leaves from the four regions. Eight major DAMs related to iridoid metabolism in P. rotata leaves were investigated: shanzhiside methyl ester, 8-epideoxyloganic acid, barlerin, shanzhiside, geniposide, agnuside, feretoside, and catalpin. In addition, five soil physical and chemical indicators in P. rotata rhizosphere soils were analyzed. Four significant positive correlations were observed between alkaline nitrogen and geniposide, exchangeable calcium and geniposide, available potassium and shanzhiside, and available phosphorus and shanzhiside methyl ester. The transcriptome data showed 12 P. rotata cDNA libraries with 74.46 Gb of clean data, which formed 29,833 unigenes. Moreover, 78.91% of the unigenes were annotated using the eight public databases. Forty-one candidate genes representing 23 enzymes involved in the biosynthesis of iridoid compounds were identified in P. rotata leaves. Moreover, the DXS1, IDI1, 8-HGO1, and G10H2 genes associated with iridoid biosynthesis were specifically expressed in P. rotata. The integration of transcriptome and metabolome analyses highlights the crucial role of soil physical and chemical indicators and major gene expression related to iridoid metabolism pathways in P. rotata from different areas. Our findings provide a theoretical foundation for exploring the molecular mechanisms underlying iridoid compound accumulation in P. rotata.

Lamiophlomis rotata (Benth.) Kudo is a perennial and unique medicinal plant of the Qinghai–Tibet Plateau. It has the effects of diminishing inflammation, activating blood circulation, removing blood stasis, reducing swelling, and relieving pain. However, thus far, reliable reference gene identifications have not been reported in wild L. rotata. In this study, we identified suitable reference genes for the analysis of gene expression related to the medicinal compound synthesis in wild L. rotata subjected to five different-altitude habitats. Based on the RNA-Seq data of wild L. rotata from five different regions, the stability of 15 candidate internal reference genes was analyzed using geNorm, NormFinder, BestKeeper, and RefFinder. TFIIS, EF-1α, and CYP22 were the most suitable internal reference genes in the leaves of L. rotata from different regions, while OBP, TFIIS, and CYP22 were the optimal reference genes in the roots of L. rotata. The reference genes identified here would be very useful for gene expression studies with different tissues in L. rotata from different habitats.