Matt Zaremsky commented above that there should be a CAT(0) cubical proof that's "easier" modulo some hard result; I think that's exactly right, and one such proof uses the rank-rigidity theorem for CAT(0) cube complexes.
The Salvetti complex $S_X$ is formed from the presentation complex of the given presentation of $A_X$ by adding an $n$--cube in the natural way whenever its $2$--skeleton appears (more precisely, see Section 2.6 in [A].)
By construction, $\pi_1S_X \cong A_X$, and the universal cover $\widetilde S_X$ is a CAT(0) cube complex on which $A_X$ acts freely, cocompactly, and without an invariant proper convex subcomplex.
Suppose that $A_X$ splits nontrivially as a direct product $G\times H$. Since, for example, finite groups of isometries of CAT(0) spaces have fixed points, $A_X$ is torsion-free and hence $G$ and $H$ are infinite. Hence no asymptotic cone of $A_X$, and hence no asymptotic cone of $\widetilde S_X$, has a cut-point (this follows from a general property of asymptotic cones, namely that any asymptotic cone of $G\times H$ is bilipschitz equivalent to the product of an asymptotic cone of $G$ with one of $H$).
By Corollary I in [B], $\widetilde S_X$ is reducible, i.e., there are CAT(0) cube complexes $C,D$ such that $\widetilde S_X$ is isomorphic to the complex $C\times D$, and neither $C$ nor $D$ is a single vertex. (This step is the application of rank rigidity.)
Label the (oriented) edges of $S_X$ by the vertices of $X$ (generators of $A_X$) in the usual way, and pull back the labels to $\widetilde S_X$. From the presentation, two edges lying on opposite sides of a square have the same label, and from one of the characterisations of hyperplanes in a CAT(0) cube complex, this gives a well-defined labelling of the hyperplanes of $\widetilde S_X$: the label of a hyperplane $h$ is the label of some (hence any) edge $e$ such that $h\cap e$ is the midpoint of $e$.
If $h_1, h_2$ are distinct hyperplanes that intersect, then $h_1\cap h_2$ contains the barycentre of a square whose edges are labelled by the label $v_1$ of $h_1$ and the label $v_2$ of $h_2$. Projecting this square down to $S_X$ tells you that $[v_1,v_2]=1$. Moreover, intersecting hyperplanes have different labels, since $S_X$ contains no square whose edges all have the same labels.
Fix vertices $c\in C,d\in D$ and consider the subcomplexes $C\times\{d\}$ and $\{c\}\times D$. By Lemma 2.5 in [B], the hyperplanes of $\widetilde S_X$ are partitioned into two subsets: $H_C$ --- those that intersect $C\times\{d\}$ --- and $H_D$ --- those that intersect $\{c\}\times D$. By the same lemma, every hyperplane in $H_C$ intersects every hyperplane in $H_D$.
In summary, we have a surjective map $\ell:H_C\sqcup H_D\to V(X)$ that sends each hyperplane to its label. Since intersecting hyperplanes have different labels, and hyperplanes in $H_C$ intersect hyperplanes in $H_D$, we get a partition $V(X) = \ell(H_C) \sqcup \ell(H_D)$ with all elements of the first subset commuting with all elements of the second. It remains to observe that $H_C$ and $H_D$ are nonempty, since $C$ and $D$ are CAT(0) cube complexes each having at least one edge.
(So from that point of view, the statement about product RAAGs is a manifestation of the fact that product decompositions of cube complexes correspond to partitions of the hyperplanes into two disjoint crossing sets, and the role of "RAAGness" is to translate that into an algebraic statement.)
[A]: Charney, Ruth, An introduction to right-angled Artin groups., Geom. Dedicata 125, 141-158 (2007).
[B]: Caprace, Pierre-Emmanuel; Sageev, Michah, Rank rigidity for CAT(0) cube complexes., Geom. Funct. Anal. 21, No. 4, 851-891 (2011).