In these proceedings, we summarize the main results from the first-ever calculations of the chiraleven and chiral-odd twist-3 parton distributions, π π (π₯) and β πΏ (π₯), of the proton from lattice QCD. We use an π π = 2 + 1 + 1 ensemble of maximally twisted mass fermions with a clover improvement. The lattice has a spatial extent of 3 fm, the lattice spacing is 0.093 fm, and the pion mass is 260 MeV. The matrix elements are obtained with a source-sink time separation of 1.12 fm to control contamination from excited states. The calculation is based on the quasi-PDF approach and employs three values for the proton momentum: 0.83 GeV, 1.25 GeV, and 1.67 GeV. The lattice data are renormalized non-perturbatively using the RI scheme, and the final results are presented in the MS scheme at the scale of 2 GeV. Furthermore, we compute in the same setup the helicity, π 1 (π₯), and transversity, β 1 (π₯), distributions, which are used to compare π π (π₯) and β πΏ (π₯) to their Wandzura-Wilczek approximations. For β πΏ (π₯), we combine results for the isovector and isoscalar flavor combinations to disentangle the individual up-and down-quark contributions.