In this work, we present a theoretical study for the collective oscillation modes, i.e. quadrupole, radial and axial mode, of a mixture of Bose and Fermi superfluids in the crossover from a Bardeen-Cooper-Schrieffer (BCS) superfluid to a molecular Bose-Einstein condensate (BEC) in harmonic trapping potentials with cylindrical symmetry of experimental interest. To this end, we start from the coupled superfluid hydrodynamic equations for the dynamics of Bose-Fermi superfluid mixtures and use the scaling theory that has been developed for a coupled system. The collective oscillation modes of Bose-Fermi superfluid mixtures are found to crucially depend on the overlap integrals of the spatial derivations of density profiles of the Bose and Fermi superfluids at equilibrium. We not only present the explicit expressions for the overlap density integrals, as well as the frequencies of the collective modes provided that the effective Bose-Fermi coupling is weak, but also test the valid regimes of the analytical approximations by numerical calculations in realistic experimental conditions. In the presence of a repulsive Bose-Fermi interaction, we find that the frequencies of the three collective modes of the Bose and Fermi superfluids are all upshifted, and the change speeds of the frequency shifts in the BCS-BEC crossover can characterize the different groundstate phases of the Bose-Fermi superfluid mixtures for different trap geometries. superfluid mixtures of 7 Li-6 Li [1, 2], 41 K-6 Li [6], and 174 Yb-6 Li [5], and the Bose-Fermi interaction gives rise to a rich behavior. In the presence of a repulsive Bose-Fermi interaction, the frequencies of the dipole oscillations of the Bose and Fermi superfluids are both downshifted in the weakly confined direction [1,5,6], and the frequency shifts increase monotonically from the BCS side to the BEC side [48,49]. In contrast, the frequency in the tight confinement for the Bose superfluid is upshifted, whereas the frequency for the Fermi superfluid is still downshifted [6]. The frequency shifts show non-monotonic and resonantlike behaviors in both directions around the BCS side [6], which may be originated from the effects of fermionic pairs breaking [50].It is naturally to ask how Bose-Fermi interaction affects the collective modes of Bose and Fermi superfluids in the BCS-BEC crossover, which is of great interest recently. However, a theoretical study for the collective oscillation modes of Bose-Fermi superfluid mixtures in a realistic experimental situation is highly nontrivial. In this work, to study quadrupole, radial and axial modes in cylindrically symmetric traps, we start from the coupled superfluid hydrodynamic equations describing the dynamics of the Bose-Fermi superfluid mixtures. The scaling method for coupled systems is then applied, and the eigenvalue equations for the coupled collective modes are obtained, which are crucially sensitive to the overlap integrals of the spatial derivations of the Bose and Fermi densities at groundstate. To present the explicit expressio...