Ultra-high field scanners allow recording BOLD images with unprecedented spatio-temporal resolution. These highly precise measurements permit studying the human brain at the mesoscale level, investigating the roles of some of the most fundamental units of neural computations, namely, layers and columns. These submillimeter recordings are however SNR-starved and the ideal choice of sequence remains to be determined. Here, we show that by using a multislab multiband 3D EPI approach we can achieve simultaneously high spatial and temporal resolution for sub-millimeter fMRI applications at 7T. We compare this approach to the more common strategies of 3D EPI and 2D MB EPI.