Effective spin-orbit (SO) Hamiltonians for conduction electrons in wurtzite heterostructures are lacking in the literature, in contrast to zincblende structures. Here we address this issue by deriving such an effective Hamiltonian valid for quantum wells, wires, and dots with arbitrary confining potentials and external magnetic fields. We start from an 8×8 Kane model accounting for the s-p z orbital mixing important to wurtzite structures, but absent in zincblende, and apply both quasi-degenerate perturbation theory (Löwdin partitioning) and the folding down approach to derive an effective 2×2 electron Hamiltonian. Focusing on wurtzite quantum wells, we later on also perform a self-consistent Poisson-Schrödinger calculation in the Hartree approximation to determine the relevant SO couplings. We obtain the usual k-linear Rashba term arising from the structural inversion asymmetry of the wells and, differently from zincblende structures, a bulk Rashba-type term induced by the inversion asymmetry of the wurtzite lattice. Our results show this latter term to be the main contribution to the Rashba coupling in wurtzite wells. We also find linear-and cubic-in-momentum Dresselhaus contributions. Both the bulk Rashba-type term and the Dresselhaus terms originate exclusively from the admixture of s-and p z -like states in wurtzites structures. Interestingly, in these systems the linear Rashba and the Dresselhaus terms have the same symmetry and can in principle cancel each other out completely, thus making the spin a conserved quantity. We determine the intrasubband (intersubband) Rashba α ν (η) and linear Dresselhaus β ν (Γ) SO strengths of GaN/AlGaN single and double wells with one and two occupied subbands (ν = 1, 2). For the GaN/Al 0.3 Ga 0.7 N single well with one occupied subband, we obtain the total spin splitting coefficient α eff 1 = α 1 + β 1 ∼ 7.16 meV·Å, in agreement with weak antilocalization measurements. In the case of two occupied subbands, we observe that the intersubband Rashba η is much weaker than the intrasubband coupling α ν . For double wells even in the presence of strong built-in electric fields (spontaneous and piezoelectric, crucial in GaN/AlGaN wells), we find a seemingly symmetric potential configuration at which both the Rashba η and Dresselhaus Γ intersubband couplings exhibit their highest strengths. On the other hand, we observe that the intrasubband Dresselhaus coefficients β 1 and β 2 interchange their values as the gate voltage V g varies across zero; a similar behavior, though less pronounced, is seen for the Rashba couplings α 1 and α 2 . We believe our general effective Hamiltonian for electrons in wurtzite heterostructures put forward here, should stimulate additional theoretical works on wurtzite quantum wells, wires, and dots with variously defined geometries and external magnetic fields.