Changes in gene expression patterns can reflect the adaptation of organisms to divergent environments. Quantitative real-time PCR (qRT-PCR) is an important tool for ecological adaptation studies at the gene expression level. The quality of the results of qRT-PCR analysis largely depends on the availability of reliable reference genes (RGs). To date, reliable RGs have not been determined for adaptive evolution studies in insects using a standard approach. Here, we evaluated the reliability of 17 candidate RGs for five populations inhabiting various altitudes of the Tibetan Plateau (TP) using four independent (geNorm, NormFinder, BestKeeper, and the deltaCt method) and one comprehensive (RefFinder) algorithms. Our results showed that α, , and were the top three most suitable RGs, and a combination of these three RGs was the most optimal for normalization. Conversely, ,, and were the most unstable RGs. The expression profiles of two target genes ( and ) were used to confirm the reliability of the chosen RGs. Additionally, the expression patterns of four other genes (,,, and ) associated with adaptation to extreme environments were assessed to explore the adaptive mechanisms of TP species to divergent environments. Each of these six target genes showed discrepant expression patterns among the five populations, suggesting that the observed expression differences may be associated with the local adaptation of to divergent altitudinal environments. This study is a useful resource for studying the adaptive evolution of TP to divergent environments using qRT-PCR, and it also acts as a guide for selecting suitable RGs for ecological and evolutionary studies in insects.