Many efforts to detect Earth-like planets around low-mass stars are presently devoted in almost every extra-solar planetsearch. M dwarfs are considered ideal targets for Doppler radial velocity searches because their low masses and luminosities makelow-mass planets orbiting in their habitable zones more easily detectable than those around higher mass stars. Nonetheless, thestatistics of frequency of low-mass planets hosted by low mass stars remains poorly constrained.Our M-dwarf radial velocity monitoring with HARPS-N within the GAPS (Global architectures of Planetary Systems) - ICE(Institut de Ciències de l'Espai CSIC-IEEC) - IAC (Instituto de Astrofísica de Canarias) projectcan provide a major contributionto the widening of the current statistics through the in-depth analysis of accurate radial velocity observations in a narrow range ofspectral sub-types (79 stars, between dM0 to dM3). Spectral accuracy will enable us to reach the precision needed to detect smallplanets with a few earth masses. Our survey will bring a contribute to the surveys devoted to the search for planets around M-dwarfs, mainly focused on the M-dwarf population of the northern emisphere, for which we will provide an estimate of the planet occurrence.We present here a long duration radial velocity monitoring of theM1 dwarf star GJ 3998 with HARPS-N to identify periodicsignals in the data. Almost simultaneous photometric observations were carried out within the APACHE and EXORAP programs tocharacterize the stellar activity and to distinguish from the periodic signals those due to activity and to the presence of planetarycompanions. We run an MCMC simulation and use Bayesian model selection to determine the number of planets in this system, toestimate their orbital parameters and minimum masses and for a proper treatment of the activity noise.The radial velocities have a dispersion in excess of their internal errors due to at least four superimposed signals, with periodsof 30.7, 13.7, 42.5 and 2.65 days. Our data are well described by a 2-planet Keplerian (13.7 d and 2.65 d) and 2 sinusoidal functions(stellar activity, 30.7 d and 42.5 d) fit. The analysis of spectral indices based on Ca II H & K and Hα lines demonstrates that theperiods of 30.7 and 42.5 days are due to chromospheric inhomogeneities modulated by stellar rotation and differential rotation. Thisresult is supported by photometry and is consistent with the results on differential rotation of M stars obtained with Kepler. Theshorter periods of 13.74 ± 0.02 d and 2.6498 ± 0.0008 d are well explained with the presence of two planets, with minimum massesof 6.26 ± 0.79M ⊕ and 2.47 ± 0.27 M ⊕ and distances of 0.089 AU and 0.029 AU from the host, respectively.