The definition of Riemmannian manifold fairly straightforwardly means that every compact smooth (I think you only need C^2) manifold can be isometrically embedded in a high enough dimensional euclidean space (with the euclidean metric) just by using the fact that every riemannian manifold is locally euclidean, which means we can cover it with balls on which there are local isomorphisms, and from these we can pick a finite subcover and partition of unity that glues together for a global isometry. Basically you are gluing together local solutions to a PDE. This is something often covered in a first year course on geometry. For the details of this approach worked out, see this expository paper: https://www.math.mcgill.ca/gantumur/math580f12/siyuan.lu.pdf
To the best of my knowledge, this result doesn't even have a name attached to it (I could be wrong, but it was always presented as just a series of unnamed propositions whenever I saw this). That's not the hard part, the hard part is putting an explicit bound on the size of the target space other than "finite". So the genius of Nash was to put a strong bound on the dimension of the target space as well as the size of the target space. Rather than using general compactness and existence arguments, this requires careful arguments and innovative geometric ideas and is really a suprising result.
I tried to find a more detailed history, but in a few minutes' searching on the web (and flipping through the textbooks I can access from home) did not find one. Somewhat surprising given interest in the subject! Possibly this is covered in, e.g., Spivak.
[1] https://en.wikipedia.org/wiki/Embedding#Riemannian_and_pseud...
[2] https://en.wikipedia.org/wiki/Nash_embedding_theorem
[3] https://en.wikipedia.org/wiki/Whitney_embedding_theorem