The remarkable properties and potential applications of Majorana fermions have led to considerable efforts in recent years to realize topological matters that host these excitations. For a number conserving system, there have been a few proposals, using either coupled chain models or a multi-component system with spin-orbit coupling (SOC), to create number fluctuation of fermion pairs in achieving Majorana fermions. In this work, we show that Majorana edge states can occur in a spinless Fermi gas in 1D lattices with tunable p-wave interaction. This is facilitated by the conversion between a pair of (open-channel) fermions and a (close channel) boson, thereby allowing the number fluctuation of fermion pairs in a single chain. This scheme requires neither SOC nor a multi-chain setup and can be implemented easily. Using the density-matrixrenormalization-group method, we have identified the Majorana phase in a wide range of parameter regimes as well as its associated phase transitions. The topological nature of the Majorana phase manifests itself in a strong edge-edge correlation in an open chain that is robust against disorder, as well as in a non-trivial winding number in the bulk generated by using twisted boundary condition. It is also shown that the Majorana phase in this system can be stable against atom losses due to few-body collisions on the same site, and can be easily identified from the fermion momentum distribution. These results pave the way for probing the intriguing Majorana physics in a simple and stable cold atoms system. † ]. It is a mode of excitation rather than a particle in the usual sense. In 2001, Kitaev showed that spinless fermions in a 1D chain coupled to a pairing field will have Majorana fermions at the ends [4]. Efforts to simulate this model in solid state matter have led to the proposal of using 1D semiconducting wires with spin-orbit coupling (SOC) in contact with a superconductor [5,6]. Similar proposals have also been made in the cold atom studies by engineering SOC on an attractive Fermi gas [7][8][9][10][11][12].Since Majorana fermions also emerge in the number conserving systems such as the Pfaffian quantum Hall state [13,14] and Kitaev's honeycomb spin model [15], there have been questions of whether proximity superconductivity (or lack of number conservation) is necessary for realizing Majaronas in 1D chains. While a single Majorana excitation cannot exist in a number conserving system, the correlation of two Majoranas at different locations (i i j 1 2 l l á ñ) is well defined. This provides a natural generalization of the presence of Majorana edge modes in a number conserving system, which is defined as a non-zero correlation of the Majoranas at the OPEN ACCESS RECEIVED