Fractional Brownian motion is a non-Markovian Gaussian process indexed by the Hurst exponent $Hin [0,1]$, generalising standard Brownian motion to account for anomalous diffusion. Functionals of this process are important for practical applications as a standard reference point for non-equilibrium dynamics. We describe a perturbation expansion allowing us to evaluate many non-trivial observables analytically: We generalize the celebrated three arcsine-laws of standard Brownian motion. The functionals are: (i) the fraction of time the process remains positive, (ii) the time when the process last visits the origin, and (iii) the time when it achieves its maximum (or minimum). We derive expressions for the probability of these three functionals as an expansion in $epsilon = H-tfrac{1}{2}$, up to second order. We find that the three probabilities are different, except for $H=tfrac{1}{2}$ where they coincide. Our results are confirmed to high precision by numerical simulations.