We review, compare and extend recent studies searching for evidence for a preferred cosmological axis. We start from the Union2 SnIa dataset and use the hemisphere comparison method to search for a preferred axis in the data. We find that the hemisphere of maximum accelerating expansion rate is in the direction $(l,b)=({309^circ}^{+23^circ}_{-3^circ}, {18^circ}^{+11^circ}_{-10^circ})$ ($omm=0.19$) while the hemisphere of minimum acceleration is in the opposite direction $(l,b)=({129^circ}^{+23^circ}_{-3^circ},{-18^circ}^{+10^circ}_{-11^circ})$ ($omm=0.30$). The level of anisotropy is described by the normalized difference of the best fit values of $omm$ between the two hemispheres in the context of lcdm fits. We find a maximum anisotropy level in the Union2 data of $frac{Delta ommax}{bomm}=0.43pm 0.06$. Such a level does not necessarily correspond to statistically significant anisotropy because it is reproduced by about $30%$ of simulated isotropic data mimicking the best fit Union2 dataset. However, when combined with the axes directions of other cosmological observations (bulk velocity flow axis, three axes of CMB low multipole moments and quasar optical polarization alignment axis), the statistical evidence for a cosmological anisotropy increases dramatically. We estimate the probability that the above independent six axes directions would be so close in the sky to be less than $1%$. Thus either the relative coincidence of these six axes is a very large statistical fluctuation or there is an underlying physical or systematic reason that leads to their correlation.