The crystal structures of 16 boron subphthalocyanines (BsubPcs) with structurally diverse axial groups were analyzed and compared to elucidate the impact of the axial group on the intermolecular π-π interactions, axial-group interactions, axial bond length and BsubPc bowl depth. π-π interactions between the isoindole units of adjacent BsubPc molecules most often involve concave-concave packing, whereas axial-group interactions with adjacent BsubPc molecules tend to favour the convex side of the BsubPc bowl. Furthermore, axial groups that contain O and/or F atoms tend to have significant hydrogen-bonding interactions, while axial groups containing arene site(s) can participate in π-π interactions with the BsubPc bowl, both of which can strongly influence the crystal packing. Bulky axial groups did tend to disrupt the π-π interactions and/or axial-group interactions, preventing some of the close packing that is seen in BsubPcs with less bulky axial groups. The atomic radius of the heteroatom bonded to boron directly influences the axial bond length, whereas the axial group has minimal impact on the BsubPc bowl depth. Finally, the crystal growth method did not generally appear to have a significant impact on the solid-state arrangement, with the exception of water occasionally being incorporated into crystal structures when hygroscopic solvents were used. These insights can help with the design and fine-tuning of the solid-state structures of BsubPcs as they continue to be developed as functional materials in organic electronics.