A plane wall-jet flow is numerically investigated and compared to experiments. The measured base flow is matched to a boundary-layer solution developing from a coupled Blasius boundary layer and Blasius shear layer. Linear stability analysis is performed, revealing high instability of two-dimensional eigenmodes and non-modal streaks. The nonlinear stage of laminar-flow breakdown is studied with three-dimensional direct numerical simulations and experimentally visualized. In the direct numerical simulation, an investigation of the nonlinear interaction between two-dimensional waves and streaks is made. The role of subharmonic waves and pairing of vortex rollers is also investigated. It is demonstrated that the streaks play an important role in the breakdown process, where their growth is transformed from algebraic to exponential as they become part of the secondary instability of the two-dimensional waves. In the presence of streaks, pairing is suppressed and breakdown to turbulence is enhanced.