A clean discussion involves the assumption and discussion of time reversal symmetry. If there are no spin orbit coupling effects, time reversal will be represented by complex conjugation and the single particle wavefunctions in a periodic potential can always be chosen real as then E(k)=E(-k) so that instead of the solutions \psi_k(x)=u_k(x)\exp(ikx) and \psi_{-k}=(\psi_k(x))^* real valued combinations can be chosen. For k=0, only one real function will be obtained.
If spin orbit coupling is taken into account, time reversal is no longer just complex conjugation so that it does not always guarantee real valuedness. This is known as Kramers degeneracy.