Branching Brownian motion, branching random walks, and the F-KPP equation have been the subject of intensive research during the last couple of decades. By means of Feynman-Kac and McKean formulas, the understanding of the maximal particles of the former two Markov processes is related to insights into the position of the front of the solution to the F-KPP equation. We will discuss some recent result on extensions of the above models to spatially random branching rates and random nonlinearities. Interestingly, the introduction of such inhomogeneities leads to a richer and much more nuanced picture when compared to the homogeneous setting.