In the Next-to-Minimal Supersymmetric Standard Model (NMSSM) with extra heavy neutrino superfields, neutrino may acquire its mass via a seesaw mechanism and sneutrino may act as a viable dark matter (DM) candidate. Given the strong tension between the naturalness for Z boson mass and the DM direct detection experiments for customary neutralino DM candidate, we augment the NMSSM with Type-I seesaw mechanism, which is the simplest extension of the theory to predict neutrino mass, and study the scenarios of sneutrino DM. We construct likelihood function with LHC Higgs data, B-physics measurements, DM relic density and its direct and indirect search limits, and perform a comprehensive scan over the parameter space of the theory by Nested Sampling method. We adopt both Bayesian and frequentist statistical quantities to illustrate the favored parameter space of the scenarios, the DM annihilation mechanism as well as the features of DMnucleon scattering. We find that the scenarios are viable over broad parameter regions, especially the Higgsino mass µ can be below about 250GeV for a significant part of the region, which predicts Z boson mass in a natural way. We also find that the DM usually co-annihilated with the Higgsinos to get the measured relic density, and consequently the DM-nucleon scattering rate is naturally suppressed to coincide with the recent XENON-1T results even for light Higgsinos. Other issues, such as the LHC search for the Higgsinos, are also addressed.1 Numerically speaking, the mixing angle should satisfy sin θν ∼ 0.02 for mν = 100GeV to predict the right DM relic density by the Z boson mediated annihilation, which corresponds to the scattering rate at the order of 10 −45 cm −2 [34]. Such a rate has been excluded by the latest XENON-1T experiment, which, on the other side, limits sin θν < 0.01 by the recent calculation in [35]. Moreover, we find that the correlation between the relic density and the scattering rate is underestimated in FIG.1 of [36].