Garnet is the second most abundant mineral phase in the upper mantle and transition zone settings. The crystal structure of garnet is quite flexible and hence it is able to accommodate various cations, including large incompatible cation such as sodium. We used first-principles simulation based on density functional theory and two widely used approximations—local density approximation (LDA) and generalized gradient approximation (GGA)—to explore the crystal structure, equation of state, and elasticity of sodium bearing majorite garnet with Na2MgSi5O12 stoichiometry at pressures relevant to the upper mantle and transition zone. We find that the pressure-volume results based on LDA can be explained by a Birch Murnaghan finite strain equation of state with V0 = 1447.6 (± 0.1) Å3, K0 = 177.4 (± 0.4) GPa, and K′0 = 3.93 (± 0.02). The results based on GGA can be explained by a Birch Murnaghan finite strain equation of state with V0 = 1525.8 (± 0.2) Å3, K0 = 160.2 (± 0.4) GPa, and K′0 = 3.96 (± 0.02). The full elastic moduli tensor for Na-majorite with tetragonal symmetry exhibits slight deviation from the cubic symmetry with C11 < C33, C12 ~ C13, and C44 ~ C66. The magnitude of the tetragonal strain also captures the slight deviation from the cubic symmetry. At pressures corresponding to the upper mantle and mantle transition zone, the compressional wave velocity, vP, and shear wave velocity, vS, for the Na-majorite garnet are fast compared to a wide variety of garnets such as pyrope, grossular, almandine, and majorite garnet. Although single-crystal anisotropy of Na-majorite is greater than pyrope, it is still low compared to the major mantle phases.