In this paper, we establish the two weight commutator of Calderón-Zygmund operators in the sense of Coifman-Weiss on spaces of homogeneous type, by studying the weighted Hardy and BMO space for A 2 weight and by proving the sparse operator domination of commutators. The main tool here is the Haar basis and the adjacent dyadic systems on spaces of homogeneous type, and the construction of a suitable version of a sparse operator on spaces of homogeneous type. As applications, we provide a two weight commutator theorem (including the high order commutator) for the following Calderón-Zygmund operators: Cauchy integral operator on R, Cauchy-Szegö projection operator on Heisenberg groups, Szegö projection operators on a family of unbounded weakly pseudoconvex domains, Riesz transform associated with the sub-Laplacian on stratified Lie groups, as well as the Bessel Riesz transforms (one-dimension and high dimension).to the subspace of functions {F b } that are boundary values of functions F ∈ H 2 (U n ). The associated Cauchy-Szegö kernel is as follows.Then it is natural to study the following question: is there a setting, by which the characterisation of two weight commutators and the related BMO space for Calderón-Zygmund operators T can be obtained, that can be applied to Calderón-Zygmund operators such as the Bessel Riesz transform, the Cauchy-Szegö projection operator on Heisenberg groups, and many other examples?To address this question we work in a general setting: spaces of homogeneous type introduced by Coifman and Weiss in the early 1970s, in [9], see also [10]. We say that (X, d, µ) is a space of homogeneous type in the sense of Coifman and Weiss if d is a quasi-metric on X and µ is a nonzero measure satisfying the doubling condition. A quasi-metric d on a set X is a function d : X × X −→ [0, ∞) satisfying (i) d(x, y) = d(y, x) ≥ 0 for all x, y ∈ X; (ii) d(x, y) = 0 if and only if x = y; and (iii) the quasi-triangle inequality: there is a constant A 0 ∈ [1, ∞) such that for all x, y, z ∈ X,