We treat the problem of option pricing under a stochastic volatility model that exhibits long-range dependence. We model the price process as a Geometric Brownian Motion with volatility evolving as a fractional Ornstein-Uhlenbeck process. We assume that the model has long-memory, thus the memory parameter H in the volatility is greater than 0.5. Although the price process evolves in continuous time, the reality is that observations can only be collected in discrete time. Using historical stock price information we adapt an interacting particle stochastic filtering algorithm to estimate the stochastic volatility empirical distribution. In order to deal with the pricing problem we construct a multinomial recombining tree using sampled values of the volatility from the stochastic volatility empirical measure. Moreover, we describe how to estimate the parameters of our model, including the long-memory parameter of the fractional Brownian motion that drives the volatility process using an implied method. Finally, we compute option prices on the S&P 500 index and we compare our estimated prices with the market option prices.