We consider a class of SUSY models in which the MSSM gauge group is supplemented with a gauged U(1)B−L symmetry and a global U(1)R symmetry. This extension introduces only electrically neutral states, and the new SUSY partners effectively double the number of states in the neutralino sector that now includes a blino (from B − L) and singlino from a gauge singlet superfield. If the DM density is saturated by a LSP neutralino, the model yields quite a rich phenomenology depending on the DM composition. The LSP relic density constraint provides a lower bound on the stop and gluino masses of about 3 TeV and 4 TeV respectively, which is testable in the near future collider experiments such as HL-LHC. The chargino mass lies between 0.24 TeV and about 2.0 TeV, which can be tested based on the allowed decay channels. We also find $$ {m}_{\tilde{\tau}1}\gtrsim $$
m
τ
˜
1
≳
500 GeV, and $$ {m}_{\tilde{e}},{m}_{\tilde{\mu}},{m}_{{\tilde{v}}^{S,P}}\gtrsim $$
m
e
˜
,
m
μ
˜
,
m
v
˜
S
,
P
≳
1 TeV. We identify chargino-neutralino coannihilation processes in the mass region 0.24 TeV $$ \lesssim {m}_{{\tilde{\upchi}}_1^0}\approx {m}_{{\tilde{\upchi}}_1^{\pm }}\lesssim $$
≲
m
χ
˜
1
0
≈
m
χ
˜
1
±
≲
1.5 TeV, and also coannihilation processes involving stau, selectron, smuon and sneutrinos for masses around 1 TeV. In addition, A2 resonance solutions are found around 1 TeV, and H2 and H3 resonance solutions are also shown around 0.5 TeV and 1 TeV . Some of the A2 resonance solutions with tan β ≳ 20 may be tested by the A/H → τ+τ− LHC searches.. While the relic density constraint excludes the bino-like DM, it is still possible to realize higgsino, singlino and blino-like DM for various mass scales. We show that all these solutions will be tested in future direct detection experiments such as LUX-Zeplin and Xenon-nT.