Objectives In preterm infants, postnatal myocardial adaptation may be complicated by bronchopulmonary dysplasia (BPD). We aimed to describe the development of left ventricular function by serial 2D, Doppler, and speckle tracking echocardiography (2D-STE) in infants with and without BPD during the neonatal period and compare these to anthropometric and conventional hemodynamic parameters. Study Design Prospective echocardiography on day of life (DOL) 1, 7, 14, and 28 in 119 preterm infants <1500 g birth weight of whom 36 developed BPD (need for oxygen supplementation at 36 weeks gestational age). Non-BPD and BPD infants differed significantly in median (IQR) gestational age (25.5(24–26.5) weeks vs. 29(27–30) weeks, p<0.001) and birth weight (661(552–871) g vs. 1100(890–1290) g, p<0.001). Results The intra- and inter-observer variability of the 2D-STE parameters measured did not depend on time of measurement, although there were significant differences in the reproducibility of the parameters. Low intra- and inter-observer variability was seen for longitudinal systolic strain and strain rate mid septum with a median CV (coefficient of variation) of <4.6%. Much higher CVs (>10%) were seen for the apical segment. While anthropometric parameters show rapid development during the first 4 weeks of life, the speckle tracking parameters did not differ statistically significantly during the neonatal period. Infants with and without BPD differed significantly (p<0.001) in the development of anthropometric parameters, conventional hemodynamic parameters except for heart rate, and 2D-STE parameters: global longitudinal systolic strain rate (GLSSR) and longitudinal systolic strain for the mid left wall (LSSR). The largest differences were seen at DOL 1 and 7 in GLSSR (p<0.001) and in LSSR (p<0.01). Conclusions Reproducible 2D-STE measurements are possible in preterm infants <1500 g. Cardiac deformation reveals early (DOL 1 and 7) ventricular changes (GLSSR and LSSR) in very low birth weight infants who develop BPD.