The reason Julia is a good fit is it has both good numerics support (arrays, matrices, autodiff, samplers, etc) and good metaprogramming support. Probabilistic programming systems really need both to feel natural and well-integrated while still being able to have nice inference algorithms.
A PPS in an environment without those two things is at a disadvantage.
A PPS in an environment without those two things is at a disadvantage.