Cytochrome P450 2B diversity and dietary novelty in the herbivorous, desert woodrat (Neotoma lepida)

The desert woodrat (Neotoma lepida) induces a diversity of biotransformation genes in response to creosote bush resin

Liver biotransformation enzymes have long been thought to enable animals to feed on diets rich in xenobiotic compounds. However, despite decades of pharmacological research in humans and rodents, little is known about hepatic gene expression in specialized mammalian herbivores feeding on toxic diets. Leveraging a recently identified population of the desert woodrat (Neotoma lepida) found to be highly tolerant to toxic creosote bush (Larrea tridentata), we explored the expression changes of suites of biotransformation genes in response to diets enriched with varying amounts of creosote resin. Analysis of hepatic RNA-seq data indicated a dose-dependent response to these compounds, including the upregulation of several genes encoding transcription factors and numerous phase I, II, and III biotransformation families. Notably, elevated expression of five biotransformation families - carboxylesterases, cytochromes P450, aldo-keto reductases, epoxide hydrolases, and UDP-glucuronosyltransferases - corresponded to species-specific duplication events in the genome, suggesting that these genes play a prominent role in N. lepida's adaptation to creosote bush. Building on pharmaceutical studies in model rodents, we propose a hypothesis for how the differentially expressed genes are involved in the biotransformation of creosote xenobiotics. Our results provide some of the first details about how these processes likely operate in the liver of a specialized mammalian herbivore.

Toxin tolerance across landscapes: Ecological exposure not a prerequisite

Little is known about the tolerances of mammalian herbivores to plant specialized metabolites across landscapes.We investigated the tolerances of two species of herbivorous woodrats, Neotoma lepida (desert woodrat) and Neotoma bryanti (Bryant's woodrat) to creosote bush (Larrea tridentata), a widely distributed shrub with a highly toxic resin. Woodrats were sampled from 13 locations both with and without creosote bush across a 900 km transect in the US southwest. We tested whether these woodrat populations consume creosote bush using plant metabarcoding of feces and quantified their tolerance to creosote bush through feeding trials using chow amended with creosote resin.Toxin tolerance was analyzed in the context of population structure across collection sites with microsatellite analyses. Genetic differentiation among woodrats collected from different locations was minimal within either species. Tolerance differed substantially between the two species, with N. lepida persisting 20% longer than N. bryanti in feeding trials with creosote resin. Furthermore, in both species, tolerance to creosote resin was similar among woodrats near or within creosote bush habitat. In both species, woodrats collected greater than 25 km from creosote had markedly lower tolerances to creosote resin compared to animals from within the range of creosote bush.The results imply that mammalian herbivores are adapted to the specialized metabolites of plants in their diet, and that this tolerance can extend several kilometers outside of the range of dietary items. That is, direct ecological exposure to the specialized chemistry of particular plant species is not a prerequisite for tolerance to these compounds. These findings lay the groundwork for additional studies to investigate the genetic mechanisms underlying toxin tolerance and to identify how these mechanisms are maintained across landscape-level scales in mammalian herbivores.

Mammalian cytochrome P450 biodiversity: Physiological importance, function, and protein and genomic structures of cytochromes P4502B in multiple species of woodrats with different dietary preferences

The vast diversity of cytochrome P450 enzymes in mammals has been proposed to result in large measure from plant-animal warfare, whereby evolution of chemical defenses such as phenolics and terpenoids in plants led to duplication and divergence of P450 genes in herbivores. Over evolutionary time, natural selection is predicted to have produced P450s with high affinity and enhanced metabolism of substrates that are ingested regularly by herbivores. Interestingly, however, almost all knowledge of the interactions of mammalian P450 enzymes with substrates stems from studies of the metabolism of drugs and model compounds rather than studies on wild mammalian herbivores and their respective PSMs. A question of particular interest centers on the role of individual P450 enzymes in the ability of certain herbivores to specialize on plants that are lethal to most other species, including those from the same genus as the specialists. We tackled this intricate problem using a tractable natural system (herbivorous woodrats, genus Neotoma) focusing on comparisons of the specialist N. stephensi, the facultative specialist N. lepida, and the generalist N. albigula, and employing a cross-disciplinary approach involving ecology, biochemistry, pharmacology, structural biology, and genomics. Based on multiple findings suggesting the importance of CYP2B enzymes for ingestion of juniper and a major constituent, α-pinene, we characterized the structure, function and activity of several CYP2B enzymes in woodrats with different dietary habits. Results to date suggest that differences in CYP2B gene copy number may contribute to differential tolerance of PSMs among woodrat species, although additional work is warranted to firmly link gene copy number to juniper tolerance.

Multigenerational backcrossing and introgression between two woodrat species at an abrupt ecological transition

When organisms experience secondary contact after allopatric divergence, genomic regions can introgress differentially depending on their relationships with adaptation, reproductive isolation, recombination, and drift. Analyses of genome-wide patterns of divergence and introgression could provide insight into the outcomes of hybridization and the potential relationship between allopatric divergence and reproductive isolation. Here, we generate population genetic data (26,262 SNPs; 353 individuals) using a reduced-representation sequencing approach to quantify patterns of ancestry, differentiation, and introgression between a pair of ecologically distinct mammals-the desert woodrat (N. lepida) and Bryant's woodrat (N. bryanti)-that hybridize at a sharp ecotone in southern California. Individual ancestry estimates confirmed that hybrids were rare in this bimodal hybrid zone, and entirely consisted of a few F₁ individuals and a broad range of multigenerational backcrosses. Genomic cline analyses indicated more than half of loci had elevated introgression from one genomic background into the other. However, introgression was not associated with relative or absolute measures of divergence, and loci with extreme values for both were not typically found near detoxification enzymes previously implicated in dietary specialization for woodrats. The decoupling of differentiation and introgression suggests that processes other than adaptation, such as drift, may underlie the extreme clines at this contact zone.

Differences in dietary composition and preference maintained despite gene flow across a woodrat hybrid zone

Ecotones, characterized by adjacent yet distinct biotic communities, provide natural laboratories in which to investigate how environmental selection influences the ecology and evolution of organisms. For wild herbivores, differential plant availability across sharp ecotones may be an important source of dietary-based selection.We studied small herbivore diet composition across a sharp ecotone where two species of woodrat, Neotoma bryanti and N. lepida, come into secondary contact with one another and hybridize. We quantified woodrat dietary preference through trnL metabarcoding of field-collected fecal pellets and experimental choice trials. Despite gene flow, parental N. bryanti and N. lepida maintain distinct diets across this fine spatial scale, and across temporal scales that span both wet and dry conditions. Neotoma bryanti maintained a more diverse diet, with Frangula californica (California coffeeberry) making up a large portion of its diet. Neotoma lepida maintains a less diverse diet, with Prunus fasciculata (desert almond) comprising more than half of its diet. Both F. californica and P. fasciculata are known to produce potentially toxic plant secondary compounds (PSCs), which should deter herbivory, yet these plants have relatively high nutritional value as measured by crude protein content. Neotoma bryanti and N. lepida consumed F. californica and P. fasciculata, respectively, in greater abundance than these plants are available on the landscape-indicating dietary selection. Finally, experimental preference trials revealed that N. bryanti exhibited a preference for F. californica, while N. lepida exhibited a relatively stronger preference for P. fasciculata. We find that N. bryanti exhibit a generalist herbivore strategy relative to N. lepida, which exhibit a more specialized feeding strategy in this study system.Our results suggest that woodrats respond to fine-scale environmental differences in plant availability that may require different metabolic strategies in order to balance nutrient acquisition while minimizing exposure to potentially toxic PSCs.

Cytochrome P450 2B diversity in a dietary specialist—the red tree vole (Arborimus longicaudus)

Although herbivores rely on liver enzymes to biotransform plant secondary metabolites ingested in plant-based diets, only a few enzymes from a handful of species have been characterized at the genomic level. In this study, we examined cytochrome P450 2B (CYP2B) sequence diversity and gene copy number in a conifer specialist, the red tree vole (Arborimus longicaudus). We fed captive individuals exclusively Douglas-fir (Pseudotsuga menziesii) foliage, cloned and sequenced their liver CYP2B cDNA, and estimated CYP2B gene copy number. We identified 21 unique CYP2B nucleotide sequences, and 20 unique CYP2B amino acid sequences. Gene copy number of CYP2B was estimated at 7.7 copies per haploid genome. We compared red tree vole CYP2B with CYP2B sequences of a generalist, the prairie vole (Microtus ochrogaster), found in GenBank. Our study revealed that the CYP2B enzymes of red tree voles possess unique sequences compared to CYP2B enzymes of other herbivorous species. The unique combination of amino acid residues at key substrate recognition sites of CYP2B enzymes may underlie the ability of the red tree vole to specialize on a highly toxic diet of Douglas-fir.

Human cytochrome P450 2B6 genetic variability in Botswana: a case of haplotype diversity and convergent phenotypes

Identification of inter-individual variability for drug metabolism through cytochrome P450 2B6 (CYP2B6) enzyme is important for understanding the differences in clinical responses to malaria and HIV. This study evaluates the distribution of CYP2B6 alleles, haplotypes and inferred metabolic phenotypes among subjects with different ethnicity in Botswana. A total of 570 subjects were analyzed for CYP2B6 polymorphisms at position 516 G > T (rs3745274), 785 A > G (rs2279343) and 983 T > C (rs28399499). Samples were collected in three districts of Botswana where the population belongs to Bantu (Serowe/Palapye and Chobe) and San-related (Ghanzi) ethnicity. The three districts showed different haplotype composition according to the ethnic background but similar metabolic inferred phenotypes, with 59.12%, 34.56%, 2.10% and 4.21% of the subjects having, respectively, an extensive, intermediate, slow and rapid metabolic profile. The results hint at the possibility of a convergent adaptation of detoxifying metabolic phenotypes despite a different haplotype structure due to the different genetic background. The main implication is that, while there is substantial homogeneity of metabolic inferred phenotypes among the country, the response to drugs metabolized via CYP2B6 could be individually associated to an increased risk of treatment failure and toxicity. These are important facts since Botswana is facing malaria elimination and a very high HIV prevalence.

Structure-Function Analysis of Mammalian CYP2B Enzymes Using 7-Substituted Coumarin Derivatives as Probes: Utility of Crystal Structures and Molecular Modeling in Understanding Xenobiotic Metabolism

Crystal structures of CYP2B35 and CYP2B37 from the desert woodrat were solved in complex with 4-(4-chlorophenyl)imidazole (4-CPI). The closed conformation of CYP2B35 contained two molecules of 4-CPI within the active site, whereas the CYP2B37 structure demonstrated an open conformation with three 4-CPI molecules, one within the active site and the other two in the substrate access channel. To probe structure-function relationships of CYP2B35, CYP2B37, and the related CYP2B36, we tested the O-dealkylation of three series of related substrates-namely, 7-alkoxycoumarins, 7-alkoxy-4-(trifluoromethyl)coumarins, and 7-alkoxy-4-methylcoumarins-with a C1-C7 side chain. CYP2B35 showed the highest catalytic efficiency (kcat/KM) with 7-heptoxycoumarin as a substrate, followed by 7-hexoxycoumarin. In contrast, CYP2B37 showed the highest catalytic efficiency with 7-ethoxy-4-(trifluoromethyl)coumarin (7-EFC), followed by 7-methoxy-4-(trifluoromethyl)coumarin (7-MFC). CYP2B35 had no dealkylation activity with 7-MFC or 7-EFC. Furthermore, the new CYP2B-4-CPI-bound structures were used as templates for docking the 7-substituted coumarin derivatives, which revealed orientations consistent with the functional studies. In addition, the observation of multiple -Cl and -NH-π interactions of 4-CPI with the aromatic side chains in the CYP2B35 and CYP2B37 structures provides insight into the influence of such functional groups on CYP2B ligand binding affinity and specificity. To conclude, structural, computational, and functional analysis revealed striking differences between the active sites of CYP2B35 and CYP2B37 that will aid in the elucidation of new structure-activity relationships.

Using the Specialization Framework to Determine Degree of Dietary Specialization in a Herbivorous Woodrat

To be considered a dietary specialist, mammalian herbivores must consume large quantities of a plant species considered "difficult" with respect to nutrient or toxin content, and possess specialized adaptations to deal with plant defensive compounds or low nutritional content. Populations of Neotoma lepida in the Great Basin consume Juniperus osteosperma, a plant heavily defended by terpenes, but a detailed dietary analysis of this population is lacking. Therefore, we investigated the extent of dietary specialization in this species in comparison with the better-studied specialist species, N. stephensi. Microhistological analysis of feces from N. lepida revealed that greater than 90% of their diet in nature was comprised of juniper. In laboratory tolerance trials, N. lepida tolerated a diet of 80% J. osteosperma, similar to that observed for N. stephensi. There was no difference in the abilities of N. lepida and N. stephensi to metabolize hexobarbital, a proxy compound for terpene metabolism. In preference tests of native and non-native juniper species, N. lepida did not exhibit a preference for its native or co-occurring juniper, J. osteosperma, over the non-native species, J. monosperma, whereas N. stephensi preferred its native or co-occurring juniper J. monosperma over non-native J. osteosperma. Behavioral and habitat differences between these woodrat species lead to the categorization of N. stephensi as an obligate juniper specialist with a small range that overlaps that of its preferred food, J. monosperma, and N. lepida as a facultative juniper specialist with a large range, and only a portion of its distribution containing populations that feed extensively on J. osteosperma.

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.

Warmer ambient temperatures depress liver function in a mammalian herbivore

Diet selection in mammalian herbivores is thought to be mainly influenced by intrinsic factors such as nutrients and plant secondary compounds, yet extrinsic factors like ambient temperature may also play a role. In particular, warmer ambient temperatures could enhance the toxicity of plant defence compounds through decreased liver metabolism of herbivores. Temperature-dependent toxicity has been documented in pharmacology and agriculture science but not in wild mammalian herbivores. Here, we investigated how ambient temperature affects liver metabolism in the desert woodrat, Neotoma lepida. Woodrats (n = 21) were acclimated for 30 days to two ambient temperatures (cool = 21°C, warm = 29°C). In a second experiment, the temperature exposure was reduced to 3.5 h. After temperature treatments, animals were given a hypnotic agent and clearance time of the agent was estimated from the duration of the hypnotic state. The average clearance time of the agent in the long acclimation experiment was 45% longer for animals acclimated to 29°C compared with 21°C. Similarly, after the short exposure experiment, woodrats at 29°C had clearance times 26% longer compared with 21°C. Our results are consistent with the hypothesis that liver function is reduced at warmer environmental temperatures and may provide a physiological mechanism through which climate change affects herbivorous mammals.

An in vivo assay for elucidating the importance of cytochromes P450 for the ability of a wild mammalian herbivore (Neotoma lepida) to consume toxic plants

An in vivo assay using the cytochrome P450 (P450) suicide inhibitor 1-aminobenzotriazole (ABT) and 24-h food intake was developed to determine the relative importance of P450s in two populations of Neotoma lepida with respect to biotransformation of plant secondary compounds in the animals' natural diets. The efficacy of ABT as a P450 inhibitor was first validated using hypnotic-state assays with and without pretreatment with ABT. Pretreatment with 100 mg/kg ABT by gavage increased hexobarbital sleep times 3-4-fold in both populations, showing effective inhibition of P450s in woodrats. Next, the Great Basin population was fed a terpene-rich juniper diet, and the Mojave population was fed a phenolic-rich creosote diet, with rabbit chow serving as the control diet in each group. Treatment with ABT inhibited food intake in the Great Basin population fed the juniper diet to a greater extent (35%) than the Great Basin population fed the control diet (19%) or the Mojave population fed the creosote diet (16%). The food intake of the Mojave population fed the control diet was not significantly inhibited by ABT. The findings suggest that the biotransformation of terpenes in juniper relies more heavily on P450s than that of phenolics in creosote. This assay provides an inexpensive and noninvasive method to explore the relative importance of P450s in the biotransformation strategies of wild mammalian herbivores.