We report on our progress in studying exotic, heavy tetraquark states, qq′ Q̅Q̅′. Using publicly available dynamical nf = 2 + 1 Wilson-Clover gauge configurations, generated by the PACS-CS collaboration, with pion masses ≃ 164, 299 and 415 MeV, we extend our previous analysis to heavy quark components containing heavier than physical bottom quarks Q̅Q̅′ = b̅′b̅′ or Q̅Q̅′ = b̅b̅′, charm and bottom quarks c̅b̅ and also only charm quarks c̅ c̅. Throughout we employ NRQCD and relativistic heavy quarks for the heavier than bottom, bottom and charm quarks. Using our previously established diquark-antidiquark and meson-meson operator basis we comment in particular on the dependence of the binding energy on the mass of the heavy quark component Q̅Q̅, with heavy quarks ranging from mQ = 0:85… 6.3 · mb. In the heavy flavor non-degenerate case, Q̅Q̅′, and especially for the tetraquark channel udc̅d̅, we extend our work to utilize a 3 × 3 GEVP to study the ground and threshold states thereby enabling a clear identification of possible binding. Finally, we present initial work on the Q̅Q̅′ = c̅c̅ system where a much larger operator basis is available in comparison to flavor combinations with NRQCD quarks.