The Linear Boltzmann Transport (LBT) model for jet propagation and interaction in quark-gluon plasma (QGP) has been used to study jet quenching in high-energy heavy-lion collisions. The suppression of single inclusive jet production, medium modification of $gamma$-jet correlation, jet profiles and fragmentation functions as observed in experiments at Large Hadron Collider (LHC) can be described well by LBT in which jet-induced medium response is shown to play an essential role. In event-by-event simulations of jet quenching within LBT, jet azimuthal anisotropies are found to correlate linearly with the anisotropic flows of bulk hadrons from the underlying hydrodynamic events.