© 2018 Institute of Physics Publishing. All rights reserved. Seismic anisotropy in the lowermost mantle has been attributed to texture development during mantle convection. This study models texture evolution in a subducting slab impinging on the core-mantle boundary. Using a 3-dimensional geodynamic model with tracers recording the deformation history, a visco-plastic self-consistent (VPSC) model for polycrystal deformation, and relying on experimentally determined slip systems, aggregate grains with volume fractions of 60% orthorhombic silicate perovskite (MgSiO3), 20% cubic calcium perovskite (CaSiO3), and 20% cubic ferropericlase ((Mg,Fe)O) were deformed plastically and developed crystal preferred orientation. Forward and reverse perovskite (Pv)-postperovskite (PPv) phase transitions were included by allowing for likely orientation variant selections. Grain orientations, P (compression) and S (shear) wave velocity pole figures were calculated for each phase as well as the aggregate. The results show that dominant (001) slip on PPv can produce strong texture and shear wave anisotropy of 0.01-3.07% with VSH >VSV which agrees with seismological observations in selected areas of the D" layer.