Xenobiotic Metabolism of Plant Secondary Compounds in Oak (Quercus Agrifolia) by Specialist and Generalist Woodrat Herbivores, Genus Neotoma

Where Does All the Poison Go? Investigating Toxicokinetics of Newt (Taricha) Tetrodotoxin (TTX) in Garter Snakes (Thamnophis)

Animals that consume toxic diets provide models for understanding the molecular and physiological adaptations to ecological challenges. Garter snakes (Thamnophis) in western North America prey on Pacific newts (Taricha), which employ tetrodotoxin (TTX) as an antipredator defense. These snakes possess mutations in voltage-gated sodium channels (Nav), the molecular targets of TTX, that decrease the binding ability of TTX to sodium channels (target-site resistance). However, genetic variation at these loci that cannot explain all the phenotypic variation in TTX resistance in Thamnophis. We explored a separate means of resistance, toxin metabolism, to determine if TTX-resistant snakes either rapidly remove TTX or sequester TTX. We examined the metabolism and distribution of TTX in the body (toxicokinetics), to determine differences between TTX-resistant and TTX-sensitive snakes in the rates at which TTX is eliminated from organs and the whole body (using TTX half-life as our metric). We assayed TTX half-life in snakes from TTX-resistant and TTX-sensitive populations of three garter snake species with a coevolutionary history with newts (T. atratus, T. couchii, T. sirtalis), as well as two non-resistant "outgroup" species (T. elegans, Pituophis catenifer) that seldom (if ever) engage newts. We found TTX half-life varied across species, populations, and tissues. Interestingly, TTX half-life was shortest in T. elegans and P. catenifer compared to all other snakes. Furthermore, TTX-resistant populations of T. couchii and T. sirtalis eliminated TTX faster (shorter TTX half-life) than their TTX-sensitive counterparts, while populations of TTX-resistant and TTX-sensitive T. atratus showed no difference rates of TTX removal (same TTX half-life). The ability to rapidly eliminate TTX may have permitted increased prey consumption, which may have promoted the evolution of additional resistance mechanisms. Finally, snakes still retain substantial amounts of TTX, and we projected that snakes could be dangerous to their own predators days to weeks following the ingestion of a single newt. Thus, aspects of toxin metabolism may have been key in driving predator-prey relationships, and important in determining other ecological interactions.

Climatic factors and population demography in big-eared woodrat, Neotoma macrotis

Changes in temperature and rainfall patterns can have marked impacts on small mammal populations that inhabit environments with highly fluctuating water availability. With projected increases in droughts and fewer but more intense rainfall events in the Southwestern United States, the persistence of many wildlife populations may be threatened. Our goal was to assess how temperature and rainfall during distinct dry and wet seasons influenced the dynamics of a population of big-eared woodrats (Neotoma macrotis) in a mixed oak woodland of coastal central California. We applied Pradel’s temporal symmetry models to our 21-year biannual capture–mark–recapture data set (1993–2014) to determine the effects of climatic factors on the woodrats’ apparent survival (Φ) and recruitment rate (f). Monthly Φ averaged 0.945 ± 0.001 and varied with season. Monthly f was 0.064 ± 0.001 in the wet season (f was fixed to 0 in the dry season). Monthly population growth rate (λ) varied from 0.996 ± 0.001 during the dry season to 1.001 ± 0.001 during the wet season, which indicated a stable population (0.999 ± 0.001). Total rainfall from the previous season and mean temperature during the same season positively influenced Φ and f. By contrast, Φ and f were negatively influenced by mean temperature from the previous season and total rainfall in the same season. The resulting λ fluctuated with total rainfall, particularly in the wet season. Our results suggest that the big-eared woodrat population may not be substantially affected by warm temperatures per se, potentially because of the microclimate provided by its stick houses. We also discuss its adaptability to local food resources and relatively slow life history relative to other cricetids, and propose that the big-eared woodrat population may be equipped to cope with future climate change.

Strategies in herbivory by mammals revisited: The role of liver metabolism in a juniper specialist (Neotoma stephensi) and a generalist (Neotoma albigula)

Although herbivory is widespread among mammals, few species have adopted a strategy of dietary specialization. Feeding on a single plant species often exposes herbivores to high doses of plant secondary metabolites (PSMs), which may exceed the animal's detoxification capacities. Theory predicts that specialists will have unique detoxification mechanisms to process high levels of dietary toxins. To evaluate this hypothesis, we compared liver microsomal metabolism of a juniper specialist, Neotoma stephensi (diet >85% juniper), to a generalist, N. albigula (diet ≤30% juniper). Specifically, we quantified the concentration of a key detoxification enzyme, cytochrome P450 2B (CYP2B) in liver microsomes, and the metabolism of α-pinene, the most abundant terpene in the juniper species consumed by the specialist woodrat. In both species, a 30% juniper diet increased the total CYP2B concentration (2-3×) in microsomes and microsomal α-pinene metabolism rates (4-fold). In N. stephensi, higher levels of dietary juniper (60% and 100%) further induced CYP2B and increased metabolism rates of α-pinene. Although no species-specific differences in metabolism rates were observed at 30% dietary juniper, total microsomal CYP2B concentration was 1.7× higher in N. stephensi than in N. albigula (p < .01), suggesting N. stephensi produces one or more variant of CYP2B that is less efficient at processing α-pinene. In N. stephensi, the rates of α-pinene metabolism increased with dietary juniper and were positively correlated with CYP2B concentration. The ability of N. stephensi to elevate CYP2B concentration and rate of α-pinene metabolism with increasing levels of juniper in the diet may facilitate juniper specialization in this species.

An association between differential expression and genetic divergence in the Patagonian olive mouse (Abrothrix olivacea)

Recent molecular studies have found striking differences between desert-adapted species and model mammals regarding water conservation. In particular, aquaporin 4, a classical gene involved in water regulation of model species, is absent or not expressed in the kidneys of desert-adapted species. To further understand the molecular response to water availability, we studied the Patagonian olive mouse Abrothrix olivacea, a species with an unusually broad ecological tolerance that exhibits a great urine concentration capability. The species is able to occupy both the arid Patagonian steppe and the Valdivian and Magellanic forests. We sampled 95 olive mouse specimens from four localities (two in the steppe and two in the forests) and analysed both phenotypic variables and transcriptomic data to investigate the response of this species to the contrasting environmental conditions. The relative size of the kidney and the ratio of urine to plasma concentrations were, as expected, negatively correlated with annual rainfall. Expression analyses uncovered nearly 3,000 genes that were differentially expressed between steppe and forest samples and indicated that this species resorts to the "classical" gene pathways for water regulation. Differential expression across biomes also involves genes that involved in immune and detoxification functions. Overall, genes that were differentially expressed showed a slight tendency to be more divergent and to display an excess of intermediate allele frequencies, relative to the remaining loci. Our results indicate that both differential expression in pathways involved in water conservation and geographical allelic variation are important in the occupation of contrasting habitats by the Patagonian olive mouse.

Evidence for functional convergence in genes upregulated by herbivores ingesting plant secondary compounds

BACKGROUND

Nearly 40 years ago, Freeland and Janzen predicted that liver biotransformation enzymes dictated diet selection by herbivores. Despite decades of research on model species and humans, little is known about the biotransformation mechanisms used by mammalian herbivores to metabolize plant secondary compounds (PSCs). We investigated the independent evolution of PSC biotransformation mechanisms by capitalizing on a dramatic diet change event-the dietary inclusion of creosote bush (Larrea tridentata)-that occurred in the recent evolutionary history of two species of woodrats (Neotoma lepida and N. bryanti).

RESULTS

By comparing gene expression profiles of two populations of woodrats with evolutionary experience to creosote and one population naïve to creosote, we identified genes either induced by a diet containing creosote PSCs or constitutively higher in populations with evolutionary experience of creosote. Although only one detoxification gene (an aldo-keto reductase) was induced by both experienced populations, these populations converged upon functionally equivalent strategies to biotransform the PSCs of creosote bush by constitutively expressing aldehyde and alcohol dehydrogenases, Cytochromes P450s, methyltransferases, glutathione S-transferases and sulfotransferases. The response of the naïve woodrat population to creosote bush was indicative of extreme physiological stress.

CONCLUSIONS

The hepatic detoxification system of mammals is notoriously complex, with hundreds of known biotransformation enzymes. The comparison herein of woodrat taxa that differ in evolutionary and ecological experience with toxins in creosote bush reveals convergence in the overall strategies used by independent species after a historical shift in diet. In addition, remarkably few genes seemed to be important in this dietary shift. The research lays the requisite groundwork for future studies of specific biotransformation pathways used by woodrats to metabolize the toxins in creosote and the evolution of diet switching in woodrats. On a larger level, this work advances our understanding of the mechanisms used by mammalian herbivores to process toxic diets and illustrates the importance of the selective relationship of PSCs in shaping herbivore diversity.

Functional characterization of cytochromes P450 2B from the desert woodrat Neotoma lepida

Mammalian detoxification processes have been the focus of intense research, but little is known about how wild herbivores process plant secondary compounds, many of which have medicinal value or are drugs. cDNA sequences that code for three enzymes of the cytochrome P450 (CYP) 2B subfamily, here termed 2B35, 2B36, and 2B37 have been recently identified from a wild rodent, the desert woodrat (Malenke et al., 2012). Two variant clones of each enzyme were engineered to increase protein solubility and to facilitate purification, as reported for CYP2B enzymes from multiple species. When expressed in Escherichia coli each of the woodrat proteins gave the characteristic maximum at 450nm in a reduced carbon monoxide difference spectrum but generally expressed at lower levels than rat CYP2B1. Two enzymes, 2B36 and 2B37, showed dealkylation activity with the model substrates 7-ethoxy-4-(trifluoromethyl)coumarin and 7-benzyloxyresorufin, whereas 2B35 was inactive. Binding of the monoterpene (+)-α-pinene produced a Type I shift in the absorbance spectrum of each enzyme. Mutation of 2B37 at residues 114, 262, or 480, key residues governing ligand interactions with other CYP2B enzymes, did not significantly change expression levels or produce the expected functional changes. In summary, two catalytic and one ligand-binding assay are sufficient to distinguish among CYP2B35, 2B36, and 2B37. Differences in functional profiles between 2B36 and 2B37 are partially explained by changes in substrate recognition site residue 114, but not 480. The results advance our understanding of the mechanisms of detoxification in wild mammalian herbivores and highlight the complexity of this system.

Xenobiotic metabolism of plant secondary compounds in the English grain aphid, Sitobion avenae (F.) (Hemiptera: Aphididae)

Plant secondary compounds have been documented to be deleterious to insects and other herbivores in diverse ways. In this study, the effect of catechol (phenolics), gramine (alkaloid) and L-ornithine-HCI (non-protein amino acid) on the activities of xenobiotic metabolizing enzymes in English grain aphid, Sitobion avenae, was evaluated. Phase I enzymes investigated in this study included carboxylesterase (CarE), and oxidoreductase, whereas Phase II enzymes were represented by glutathione S-transferase (GST). In general, CarE and GST activities in S. avenae were positively correlated with the concentration of plant secondary compounds in artificial diets. Oxidoreductase activity, however, displayed a different profile. Specifically, peroxidase (POD) and polyphenol oxidase (PPO) activities in S. avenae were positively correlated with concentrations of dietary catechol and gramine, respectively, whereas catalase (CAT) activity was significantly suppressed by the higher concentration of catechol, gramine and L-ornithine-HCl. These combined results suggest that CarE and GST in S. avenae are key enzymes to breakdown a broad spectrum of plant secondary compounds, whereas oxidoreductase, including PPO and POD, degrades specific groups of plant secondary compounds.

Transcriptome sequencing and microarray development for the woodrat (Neotoma spp.): custom genetic tools for exploring herbivore ecology

Massively parallel sequencing has enabled the creation of novel, in-depth genetic tools for nonmodel, ecologically important organisms. We present the de novo transcriptome sequencing, analysis and microarray development for a vertebrate herbivore, the woodrat (Neotoma spp.). This genus is of ecological and evolutionary interest, especially with respect to ingestion and hepatic metabolism of potentially toxic plant secondary compounds. We generated a liver transcriptome of the desert woodrat (Neotoma lepida) using the Roche 454 platform. The assembled contigs were well annotated using rodent references (99.7% annotation), and biotransformation function was reflected in the gene ontology. The transcriptome was used to develop a custom microarray (eArray, Agilent). We tested the microarray with three experiments: one across species with similar habitat (thus, dietary) niches, one across species with different habitat niches and one across populations within a species. The resulting one-colour arrays had high technical and biological quality. Probes designed from the woodrat transcriptome performed significantly better than functionally similar probes from the Norway rat (Rattus norvegicus). There were a multitude of expression differences across the woodrat treatments, many of which related to biotransformation processes and activities. The pattern and function of the differences indicate shared ecological pressures, and not merely phylogenetic distance, play an important role in shaping gene expression profiles of woodrat species and populations. The quality and functionality of the woodrat transcriptome and custom microarray suggest these tools will be valuable for expanding the scope of herbivore biology, as well as the exploration of conceptual topics in ecology.

Biotransformation enzyme expression in the nasal epithelium of woodrats

When herbivores come in contact with volatile plant secondary compounds (PSC) that enter the nasal passages the only barrier between the nasal cavity and the brain is the nasal epithelium and the biotransformation enzymes present there. The expression of two biotransformation enzymes Cytochrome P450 2B (CYP2B) and glutathione-S-transferase (GST) was investigated in the nasal epithelia and livers of three populations of woodrats. One population of Neotoma albigula was fed juniper that contains volatile terpenes. Juniper caused upregulation of CYP2B and GST in the nasal epithelium and the expression of CYP2B and GST in the nasal epithelium was correlated to liver expression, showing that the nasal epithelia responds to PSC and the response is similar to the liver. Two populations of Neotoma bryanti were fed creosote that contains less volatile phenolics. The creosote naive animals upregulated CYP2B in their nasal epithelia while the creosote experienced animals upregulated GST. There was no correlation between CYP2B and GST expression in the nasal epithelia and livers of either population. The response of the nasal epithelium to PSC seems to be an evolved response that is PSC and experience dependent.

Mechanisms for eliminating monoterpenes of sagebrush by specialist and generalist rabbits

Pygmy rabbits (Brachylagus idahoensis) are one of only three vertebrates that subsist virtually exclusively on sagebrush (Artemisia spp.), which contains high levels of monoterpenes that can be toxic. We examined the mechanisms used by specialist pygmy rabbits to eliminate 1,8-cineole, a monoterpene of sagebrush, and compared them with those of cottontail rabbits (Sylvilagus nuttalli), a generalist herbivore. Rabbits were offered food pellets with increasing concentrations of cineole, and we measured voluntary intake and excretion of cineole metabolites in feces and urine. We expected pygmy rabbits to consume more, but excrete cineole more rapidly by using less-energetically expensive methods of detoxification than cottontails. Pygmy rabbits consumed 3-5 times more cineole than cottontails relative to their metabolic body mass, and excreted up to 2 times more cineole metabolites in their urine than did cottontails. Urinary metabolites excreted by pygmy rabbits were 20 % more highly-oxidized and 6 times less-conjugated than those of cottontails. Twenty percent of all cineole metabolites recovered from pygmy rabbits were in feces, whereas cottontails did not excrete fecal metabolites. When compared to other mammals that consume cineole, pygmy rabbits voluntarily consumed more, and excreted more cineole metabolites in feces, but they excreted less oxidized and more conjugated cineole metabolites in urine. Pygmy rabbits seem to have a greater capacity to minimize systemic exposure to cineole than do cottontails, and other cineole-consumers, by minimizing absorption and maximizing detoxification of ingested cineole. However, mechanisms that lower systemic exposure to cineole may come with a higher energetic cost in pygmy rabbits than in other mammalian herbivores.

Significant interspecies differences in induction profiles of hepatic CYP enzymes by TCDD in bank and field voles

The gene expression and induction of cytochrome P450 (CYP)-enzymes following 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) peroral administration was studied in the livers of two wild vole species--the bank vole (Myodes glareolus) and the field vole (Microtus agrestis). The dioxin-sensitive C57BL/6 mouse was used as a reference. Doses of 0.05, 0.5, 5.0, and 50 µg/kg were applied to ascertain a dose-response relationship, and the dose of 50 µg/kg was applied to the study time course for up to 96 h. The cytochrome P450 1A1 (CYP1A1) mRNA expression showed an expected dose-dependent increase equally in both vole species. Bank voles expressed notably higher CYP2A mRNA levels as compared with field voles. Both species exhibited dose-dependent increases in putative CYP1A-, CYP2B-, and CYP2A-associated activities as measured by fluorometric assays for ethoxyresorufin-O-deethylase (EROD), penthoxyresorufin-O-depenthylase (PROD), and 7-ethoxycoumarin-O-deethylase (ECOD), respectively. Putative CYP2A-associated coumarin-7-hydroxylase (COH) activity showed a slight increase at the two highest doses of TCDD in field voles but not in bank voles, and their basal COH activity was only one-fourth or less of that in field voles. Overall, however, bank voles tended to exhibit higher CYP-associated enzyme activities measured at the two largest doses of TCDD than field voles. A western blot analysis of aryl hydrocarbon receptor (AhR) revealed that the two vole species had differential band patterns, suggesting dissimilar structures for their AhRs.