Gaining insight into the behavior of bacteria at the single-cell level is important given that heterogeneous microenvironments strongly influence microbial physiology. The hybridization chain reaction (HCR) is a technique that provides in situ molecular signal amplification, enabling simultaneous mapping of multiple target RNAs at small spatial scales. To refine this method for biofilm applications, we designed and validated new probes to visualize the expression of key catabolic genes in Pseudomonas aeruginosa aggregates. In addition to using existing probes for the dissimilatory nitrate reductase (narG), we developed probes for a terminal oxidase (ccoN1), nitrite reductase (nirS), nitrous oxide reductase (nosZ), and acetate kinase (ackA). These probes can be used to determine gene expression levels across heterogeneous populations such as biofilms. Using these probes, we quantified gene expression across oxygen gradients in aggregate populations grown using the agar block biofilm assay (ABBA). We observed distinct patterns of catabolic gene expression, with upregulation occurring in particular ABBA regions both within individual aggregates and over the aggregate population. Aerobic respiration (ccoN1) showed peak expression under oxic conditions, whereas fermentation (ackA) showed peak expression in the anoxic cores of high metabolic activity aggregates near the air-agar interface. Denitrification genes narG, nirS, and nosZ showed peak expression in hypoxic and anoxic regions, although nirS expression remained at peak levels deeper into anoxic environments than other denitrification genes. These results reveal that the microenvironment correlates with catabolic gene expression in aggregates, and they demonstrate the utility of HCR in unveiling cellular activities at the microscale level in heterogeneous populations.