Medicinal plants contain valuable compounds that have attracted worldwide interest for their use in the production of natural drugs. The presence of compounds such as rosmarinic acid, carnosic acid, and carnosol in Rosmarinus officinalis has made it a plant with unique therapeutic effects. The identification and regulation of the biosynthetic pathways and genes will enable the large-scale production of these compounds. Hence, we studied the correlation between the genes involved in biosynthesis of the secondary metabolites in R. officinalis using proteomics and metabolomics data by WGCNA. We identified three modules as having the highest potential for the metabolite engineering. Moreover, the hub genes highly connected to particular modules, TFs, PKs, and transporters were identified. The TFs of MYB, C3H, HB, and C2H2 were the most likely candidates associated with the target metabolic pathways. The results indicated that the hub genes including Copalyl diphosphate synthase (CDS), Phenylalanine ammonia lyase (PAL), Cineole synthase (CIN), Rosmarinic acid synthase (RAS), Tyrosine aminotransferase (TAT), Cinnamate 4-hydroxylase (C4H), and MYB58 are responsible for biosynthesis of important secondary metabolites. Thus, we confirmed these results using qRT-PCR after treating R. officinalis seedlings with methyl jasmonate. These candidate genes may be employed for genetic and metabolic engineering research to increase R. officinalis metabolite production.