We report the discovery of the energy-dependent morphology of the GeV gamma-ray emission from HESS J1857+026 with more than 13 yr of Fermi Large Area Telescope data. The GeV gamma-ray emission from this region is composed of two extended components. The hard component with an index of 1.74 ± 0.07 in the energy range of 0.5–500 GeV is spatially coincident with HESS J1857+026, and its 68% containment radius varies from ∼0.°44 below 40 GeV to ∼0.°30 above 140 GeV. The hard GeV gamma-ray spectrum and the energy-dependent morphology of HESS J1857+026 make it favor a pulsar wind nebula origin, which is associated with the energetic pulsar PSR J1856+0245. The soft component with an index of 2.70 ± 0.16 and another extended gamma-ray source detected in this region, 4FGL J1857.9+0313e, with an index of 2.55 ± 0.07, are spatially coincident with two molecular clumps in the northeast and southwest of HESS J1857+026, which favors the hadronic process, and the protons could be accelerated by the hypothetical supernova remnant associated with PSR J1856+0245.