We re-examine the contentious question of constraints on anisotropic expansion from Type Ia supernovae (SNIa) in the light of a novel determination of peculiar velocities, which are crucial to test isotropy with SNe out to distances ∼ < 200/ℎ Mpc. We re-analyze the Joint Light-Curve Analysis (JLA) Supernovae (SNe) data, improving on previous treatments of peculiar velocity corrections and their uncertainties (both statistical and systematic) by adopting state-of-the-art flow models constrained independently via the 2M++ galaxy redshift compilation. We also introduce a novel procedure to account for colour-based selection effects, and adjust the redshift of low-𝑧 SNe self-consistently in the light of our improved peculiar velocity model.We adopt the Bayesian hierarchical model BAHAMAS to constrain a dipole in the distance modulus in the context of the ΛCDM model and the deceleration parameter in a phenomenological Cosmographic expansion. We do not find any evidence for anisotropic expansion, and place a tight upper bound on the amplitude of a dipole, |𝐷 𝜇 |< 5.93 × 10 −4 (95% credible interval) in a ΛCDM setting, and |𝐷 𝑞 0 |< 6.29 × 10 −2 in the Cosmographic expansion approach. Using Bayesian model comparison, we obtain posterior odds in excess of 900:1 (640:1) against a constant-in-redshift dipole for ΛCDM (the Cosmographic expansion). In the isotropic case, an accelerating universe is favoured with odds of ∼ 1100 : 1 with respect to a decelerating one.