An asymptotic preserving automatic domain decomposition method for the Vlasov-Poisson-BGK system with applications to plasmas