Impact of high predation risk on genome-wide hippocampal gene expression in snowshoe hares

The population dynamics of snowshoe hares (Lepus americanus) are fundamental to the ecosystem dynamics of Canada’s boreal forest. During the 8- to 11-year population cycle, hare densities can fluctuate up to 40-fold. Predators in this system (lynx, coyotes, greathorned owls) affect population numbers not only through direct mortality but also through sublethal effects. The chronic stress hypothesis posits that high predation risk during the decline severely stresses hares, leading to greater stress responses, heightened ability to mobilize cortisol and energy, and a poorer body condition. These effects may result in, or be mediated by, differential gene expression. We used an oligonucleotide microarray designed for a closely-related species, the European rabbit (Oryctolagus cuniculus), to characterize differences in genome-wide hippocampal RNA transcript abundance in wild hares from the Yukon during peak and decline phases of a single cycle. A total of 106 genes were differentially regulated between phases. Array results were validated with quantitative realtime PCR, and mammalian protein sequence similarity was used to infer gene function. In comparison to hares from the peak, decline phase hares showed increased expression of genes involved in metabolic processes and hormone response, and decreased expression of immune response and blood cell formation genes. We found evidence for predation risk effects on the expression of genes whose putative functions correspond with physiological impacts known to be induced by predation risk in snowshoe hares. This study shows, for the first time, a link between changes in demography and alterations in neural RNA transcript abundance in a natural population.