Role of cytochrome P450 2B sequence variation and gene copy number in facilitating dietary specialization in mammalian herbivores

Theory postulates that dietary specialization in mammalian herbivores is enabled by a specialized set of liver enzymes that process the high concentrations of similar plant secondary metabolites (PSMs) in the diets of specialists. To investigate whether qualitative and quantitative differences in detoxification mechanisms distinguish dietary specialists from generalists, we compared the sequence diversity and gene copy number of detoxification enzymes in two woodrat species: a generalist, the white‐throated woodrat (Neotoma albigula) and a juniper specialist, Stephens’ woodrat (N. stephensi). We focused on enzymes in the cytochrome P450 subfamily 2B (CYP2B), because previous research suggests this subfamily plays a key role in the processing of PSMs. For both woodrat species, we obtained and sequenced CYP2B cDNA, generated CYP2B phylogenies, estimated CYP2B gene copy number and created a homology model of the active site. We found that the specialist possessed on average ~5 more CYP2B gene copies than the generalist, but the specialist's CYP2B sequences were less diverse. Phylogenetic analysis of putative CYP2B homologs resolved woodrat species as reciprocally monophyletic and suggested evolutionary convergence of distinct homologs on similar key amino acid residues in both species. Homology modelling of the CYP2B enzyme suggests that interspecific differences in substrate preference and function likely result from amino acid differences in the enzyme active site. The characteristics of CYP2B in the specialist, that is greater gene copy number coupled with less sequence variation, are consistent with specialization to a narrow range of dietary toxins.