We discuss implementation of the LHC experimental data sets in the new CT18 global analysis of quantum chromodynamics (QCD) at the next-to-next-leading order of the QCD coupling strength. New methodological developments in the fitting methodology are discussed. Behavior of the CT18 NNLO PDFs for the conventional and saturation-inspired factorization scales in deep-inelastic scattering is reviewed. Four new families of (N)NLO CTEQ-TEA PDFs are presented: CT18, A, X, and Z.