  This talk provides an update on the calculation of matrix elements of flavor
diagonal axial, scalar and tensor quark bilinear operators between the nucleon
ground state. The simulations are done using Wilson-clover fermions on a sea of
eight 2+1+1-flavor HISQ ensembles generated by the MILC collaboration. We
discuss the signal in the connected and disconnected contributions, calculation
of the renormalization constants and mixing in the RI-sMOM scheme, and control
over the simultaneous chiral-continuum-finite-volume fit used to extract the
final charges.
