The phase behavior of hard boardlike biaxial particles of relative dimensions close to the clamitic to discotic crossover is explored by means of Monte Carlo molecular simulations. Transitions between two distinct biaxial nematic phases as well as transitions from a biaxial nematic to a uniaxial Sm-A phase are obtained. The formation of anisotropic supramolecular assemblies is demonstrated and is quantified by means of rotationally invariant pair correlation functions.