We present lattice QCD results for masses and magnetic polarizabilities of light and strange pseudoscalar mesons, chiral condensates, decay constants of neutral pion, and neutral kaon in the presence of background magnetic fields with e B ranging up to around 3.35 GeV2( ~ 70 M2π) in the vacuum. The computations were carried out in (2 + 1)-flavor QCD mostly on 323 × 96 lattices using the highly improved staggered quark action with Mπ ≈ 220 MeV at zero temperature. We find that the masses of neutral pseudoscalar mesons monotonously decrease as the magnetic field strength grows and then saturate at a nonzero value, while there exists a nonmonotonous behavior of charged pion and kaon masses in the magnetic field. We observe a qB scaling of the up and down quark flavor components of neutral pion mass, neutral pion decay constant as well as the quark chiral condensates at 0.05 ≲ eB ≲ 3.35 GeV2. We show that the correction to the Gell-Mann-Oakes-Renner relation involving the neutral pion is less than 6% and the correction for the relation involving neutral kaon is less than 30% at e B ≲ 3.35 GeV2 . We also derive the Ward-Takahashi identities for QCD in the magnetic field in the continuum formulation including the relation between integrated neutral pseudoscalar meson correlators and chiral condensates.