The neocortex contains diverse populations of excitatory neurons segregated by layer and further definable by their specific cortical and subcortical projection targets. The current study describes a systematic approach to identify molecular correlates of specific projection neuron classes in mouse primary somatosensory cortex (S1), using a combination of in situ hybridization (ISH) data mining, marker gene colocalization, and combined retrograde labeling with ISH for layer-specific marker genes. First, we identified a large set of genes with specificity for each cortical layer, and that display heterogeneous patterns within those layers. Using these genes as markers, we find extensive evidence for the covariation of gene expression and projection target specificity in layer 2/3, 5, and 6, with individual genes labeling neurons projecting to specific subsets of target structures. The combination of gene expression and target specificity imply a great diversity of projection neuron classes that is similar to or greater than that of GABAergic interneurons. The covariance of these 2 phenotypic modalities suggests that these classes are both discrete and genetically specified.