Recent progress in gyrokinetic particle-in-cell simulations of turbulent plasmas using the gyrokinetic toroidal code (GTC) is surveyed. In particular, recent results for electron temperature gradient (ETG) modes and their resulting transport are presented. Also, turbulence spreading, and the effects of the parallel nonlinearity, are described. The GTC code has also been generalized for non-circular plasma cross-section, and initial results are presented. In addition, two distinct methods of generalizing the GTC code to be electromagnetic are described, along with preliminary results. Finally, a related code, GTC-Neo, for calculating neoclassical fluxes, electric fields, and velocities, are described.