Error with MCMCThreads()

[mcreel]

The following code runs fine with Turing v0.20.4, but errors with v0.21.1. The only package installed in the project is Turing

using Turing, LinearAlgebra


lb = [0.0, 0.0, 0.0, 0.0]
ub = [3.0, 3.0, 1.0, 3.0] 
m = [1.0, 1.0, 0.2, 1.8] 

k = size(m,1)
V = 1.0* Matrix(I, k, k)

macro Prior()
    return :( arraydist([Uniform(lb[i], ub[i]) for i = 1:size(lb,1)]) )
end

@model function MSM(m)
    θ ~ @Prior
    m ~ MvNormal(θ, V)
end

length = 1000
nchains = 4
# This does not work
chain = sample(MSM(m), NUTS(), MCMCThreads(), length, nchains; init_params=m)

## single thread, this works
chain = sample(MSM(m), NUTS(), length; init_params=m)

[mcreel]

The solution, pointed out by @sethaxen on Julia Discourse, is in the AbstractMCMC docs:
init_params (default: nothing): if set to init_params !== nothing, then the ith element of init_params is used as initial parameters of the ith chain. If one wants to use the same initial parameters x for every chain, one can specify e.g. init_params = Iterators.repeated(x) or init_params = FillArrays.Fill(x, N).