The manipulation of charge and lattice boundary conditions of atomically precise low-dimensional complex oxide heterostructures can stabilize ferroelectric topologies such as polar vortices and skyrmion-like bubble domains. Recent discovery the continuous rotation of clockwise/counterclockwise ferroelectric vortices were observed in ferroelectric/paraelectric (PbTiO₃)n/(SrTiO₃)n superlattices on the DyScO₃ substrate. In this work, we reduce the complexity of the ferroelectric/paraelectric superlattice to the study of a model system of a single confined PbTiO₃ layer on different substrates. We will present experimental-theoretical study of the emergence of complex topologies of electrical polarization. Phase-field modeling predicts a systematic evolution of the polarization states corresponding to different strain conditions of the substrates and second principle calculation predicts the polar distribution of electrical topologies. By changing strain conditions, we characterize from polar vortices to skyrmion-like bubble domains by high resolution scanning transmission electron microscopy in conjunction with piezoelectric force microscopy study of ferroelectric structure and reciprocal space mapping study of domain periodicity.