Focus on Weed Control: Dual Function of the Cytochrome P450 CYP76 Family from Arabidopsis thaliana in the Metabolism of Monoterpenols and Phenylurea Herbicides

Comparative genomics analysis unravels lineage-specific bursts of gene duplications related to the emergence of specialized pathways. The CYP76C subfamily of cytochrome P450 enzymes is specific to Brassicaceae. Two of its members were recently associated with monoterpenol metabolism. This prompted us to investigate the CYP76C subfamily genetic and functional diversification. Our study revealed high rates of CYP76C gene duplication and loss in Brassicaceae, suggesting the association of the CYP76C subfamily with species-specific adaptive functions. Gene differential expression and enzyme functional specialization in Arabidopsis thaliana, including metabolism of different monoterpenols and formation of different products, support this hypothesis. In addition to linalool metabolism, CYP76C1, CYP76C2, and CYP76C4 metabolized herbicides belonging to the class of phenylurea. Their ectopic expression in the whole plant conferred herbicide tolerance. CYP76Cs from A. thaliana. thus provide a first example of promiscuous cytochrome P450 enzymes endowing effective metabolism of both natural and xenobiotic compounds. Our data also suggest that the CYP76C gene family provides a suitable genetic background for a quick evolution of herbicide resistance.Although extensive monoterpenol (especially linalool) oxidative metabolism has been described in many plant species, leading to fragrant and bioactive compounds as diverse as alcohols, aldehydes, acids, and epoxides (Williams et al., 1982; Matich et al., 2003, 2011; Luan et al., 2005, 2006; Ginglinger et al., 2013), pyranoid or furanoid linalool derivatives (Pichersky et al., 1994; Raguso and Pichersky, 1999), and geraniol-derived iridoids and secoiridoids (Dinda et al., 2007a, 2007b, 2011; Tundis et al., 2008), limited information is available on the enzymes generating these oxygenated compounds. Involvement of a cytochrome P450 (P450) enzyme extracted from Vinca rosea (now renamed Catharanthus roseus) in the hydroxylation of geraniol and nerol was suggested as early as 1976 (Madyastha et al., 1976). The first plant P450 gene to be isolated, CYP71A1 from avocado (Persea americana) fruit, was later shown to encode an enzyme with geraniol/nerol epoxidase activity (Hallahan et al., 1992, 1994). To our knowledge, a connection with compounds formed in the fruit has not yet been established. The geraniol 8-hydroxylase (often named geraniol 10-hydroxylase) CYP76B6, involved in the biosynthesis of secoiridoids and monoterpene indole alkaloid anticancer drugs in C. roseus, was found to belong to the CYP76 family in 2001 (Collu et al., 2001). The catalytic function of this enzyme was recently revised, and was shown to include a second oxidation activity, the conversion of 8-hydroxygeraniol into 8-oxogeraniol (Höfer et al., 2013). The same work also revealed a geraniol 8- and 9-hydroxylase activity of CYP76C4 from Arabidopsis thaliana. More recently, another CYP76 enzyme (CYP76A226) from C. roseus was found to metabolize oxidized geraniol derivatives and to have an iridoid oxidase activity, catalyzing the triple oxygenation of cis-trans-nepetalactol into 7-deoxyloganetic acid for the biosynthesis of secoiridoids and terpene indole alkaloids (Miettinen et al., 2014; Salim et al., 2014). Not all CYP76 enzymes seem to be devoted to the metabolism of monoterpenols. In most cases, however, CYP76s seem to be involved in terpenoid metabolism. CYP76Ms from monocots were found to metabolize diterpenoids for the synthesis of antifungal phytocassanes (Swaminathan et al., 2009; Wang et al., 2012; Wu et al., 2013), CYP76AH1 from Salvia miltiorhizza and its ortholog CYP76AH4 from rosemary (Rosmarinus officinalis) were shown to hydroxylate the norditerpene abietatriene in the pathway to labdane-related compounds (Zi and Peters, 2013), whereas CYP76Fs from sandalwood (Santalum album) were found to hydroxylate the sesquiterpenes santalene and bergamotene (Diaz-Chavez et al., 2013). CYP76B1 from Helianthus tuberosus was, however, found to metabolize herbicides belonging to the class of phenylurea (Robineau et al., 1998; Didierjean et al., 2002), but its physiological function was not reported. Other P450s from soybean (Glycine max; CYP71A10; Siminszky et al., 1999) or tobacco (Nicotiana tabacum; CYP71A11 and CYP81B1; Yamada et al., 2000) were also reported to metabolize phenylurea, but their physiological function was not investigated.A. thaliana ecotype Columbia-0 (Col-0) emits no geraniol and only tiny amounts of linalool, and extensive volatile profiling of different tissues detected only minor amounts of lilac aldehydes (oxygenated linalool derivatives; Rohloff and Bones, 2005). However, ectopic expression of a linalool/nerolidol synthase of strawberry (Fragaria × anannasa cv Elsanta) revealed a potentially efficient oxidative linalool metabolism in A. thaliana rosette leaves (Aharoni et al., 2003). Only recent work started to explore linalool metabolism in A. thaliana, which was found mainly localized in the flowers (Ginglinger et al., 2013). This work demonstrated the existence of two linalool synthases producing different enantiomers, and the concomitant involvement of two P450 enzymes, CYP76C3 and CYP71B31, with predominance of CYP76C3, in linalool oxidation. It also suggested the presence of partially redundant enzymes that may contribute to floral linalool metabolism.A family of eight CYP76 genes is detected in the A. thaliana genome. We report here an evolutionary and functional analysis of this family. We show that members of the CYP76C subfamily, when successfully expressed in yeast (Saccharomyces cerevisiae), all metabolize monoterpenols with different substrate specificities. Although CYP76Cs seem specific to Brassicaceae, they share common functions with CYP76s from other plants, such as CYP76B1 from H. tuberosus and CYP76B6 from C. roseus. These functions include not only monoterpenol oxidation, but also metabolism and detoxification of herbicides belonging to the class of phenylurea. Because of this property, CYP76Cs can be used simultaneously for monoterpenol oxidation and as selectable markers for plant transformation.     

Distinct Amino Acids in the C-Linker Domain of the Arabidopsis K+ Channel KAT2 Determine Its Subcellular Localization and Activity at the Plasma Membrane

Shaker K⁺ channels form the major K⁺ conductance of the plasma membrane in plants. They are composed of four subunits arranged around a central ion-conducting pore. The intracellular carboxy-terminal region of each subunit contains several regulatory elements, including a C-linker region and a cyclic nucleotide-binding domain (CNBD). The C-linker is the first domain present downstream of the sixth transmembrane segment and connects the CNBD to the transmembrane core. With the aim of identifying the role of the C-linker in the Shaker channel properties, we performed subdomain swapping between the C-linker of two Arabidopsis (Arabidopsis thaliana) Shaker subunits, K⁺ channel in Arabidopsis thaliana2 (KAT2) and Arabidopsis thaliana K⁺ rectifying channel1 (AtKC1). These two subunits contribute to K⁺ transport in planta by forming heteromeric channels with other Shaker subunits. However, they display contrasting behavior when expressed in tobacco mesophyll protoplasts: KAT2 forms homotetrameric channels active at the plasma membrane, whereas AtKC1 is retained in the endoplasmic reticulum when expressed alone. The resulting chimeric/mutated constructs were analyzed for subcellular localization and functionally characterized. We identified two contiguous amino acids, valine-381 and serine-382, located in the C-linker carboxy-terminal end, which prevent KAT2 surface expression when mutated into the equivalent residues from AtKC1. Moreover, we demonstrated that the nine-amino acid stretch ₃₁₂TVRAASEFA₃₂₀ that composes the first C-linker α-helix located just below the pore is a crucial determinant of KAT2 channel activity. A KAT2 C-linker/CNBD three-dimensional model, based on animal HCN (for Hyperpolarization-activated, cyclic nucleotide-gated K⁺) channels as structure templates, has been built and used to discuss the role of the C-linker in plant Shaker inward channel structure and function.In plants, potassium channels from the Shaker family dominate the plasma membrane (PM) conductance to K⁺ in most cell types and play crucial roles in sustained K⁺ transport (Blatt et al., 2012; Hedrich, 2012; Sharma et al., 2013). Plant Shaker channels, like their homologs in animals (Craven and Zagotta, 2006; Wahl-Schott and Biel, 2009), belong to the six transmembrane-one pore (6TM-1P) cation channel superfamily. Functional channels are tetrameric proteins arranged around a central pore (Daram et al., 1997; Urbach et al., 2000; Dreyer et al., 2004). These channels can result from the assembly of Shaker subunits encoded by the same gene (homotetramers) or by different Shaker genes (heterotetramers). Heterotetramerization has been extensively reported within the inwardly rectifying Shaker channel group (five members in Arabidopsis [Arabidopsis thaliana]) and increased channel functional diversity (Jeanguenin et al., 2008; Lebaudy et al., 2008a).Based on in silico sequence analyses, plant Shaker subunits display a short cytosolic N-terminal domain, followed by the 6TM-1P hydrophobic core, and a long C-terminal cytosolic region in which several domains can be identified. The first one, named C-linker (about 80 amino acids in length), is followed by a cyclic nucleotide-binding domain (CNBD), an ankyrin domain (involved in protein-protein interaction; Lee et al., 2007, Grefen and Blatt, 2012), and a domain named K_HA (Ehrhardt et al., 1997) rich in hydrophobic and acidic residues. Sequence analysis of plant Shaker channels indicates that, among these cytosolic domains, the highest levels of similarity are displayed by the C-linker and the CNBD domains. Interestingly, both domains are also highly conserved in some members from the animal K⁺ channel superfamily, like Hyperpolarization-activated, cyclic nucleotide-gated K⁺ channel (HCN), K⁺ voltage-gated channel, subfamily H (KCNH), and Cyclic-nucleotide-gated ion channel (CNGC). In these animal 6TM-1P channels, the roles of C-linker and CNBD domains have been extensively investigated via crystal structure analyses (Zagotta et al., 2003; Brelidze et al., 2012), whereas plant Shaker channels are still poorly characterized at the structural level (Dreyer et al., 2004; Gajdanowicz et al., 2009; Naso et al., 2009; Garcia-Mata et al., 2010).Aiming at investigating the structure-function relationship of plant Shaker channels, we have used the Arabidopsis Shaker subunit K⁺ channel in Arabidopsis thaliana2 (KAT2) as a model. We developed a subdomain-swapping strategy between KAT2 and another Shaker subunit displaying distinctive features, Arabidopsis thaliana K⁺ rectifying channel1 (AtKC1). The KAT2 subunit can form homomeric or heteromeric inwardly rectifying K⁺ channels at the PM and has been shown to be strongly expressed in guard cells, where it provides a major contribution to the membrane conductance to K⁺ (Pilot et al., 2001; Lebaudy et al., 2008b). In contrast, the behavior of AtKC1 is more complex. In planta, this subunit is coexpressed with other inwardly rectifying Shaker subunits, including KAT2, in different plant tissues (Jeanguenin et al., 2011), and in roots, direct evidence has been obtained that AtKC1 is involved in functional heterotetrameric channel formation with AKT1 (Reintanz et al., 2002; Honsbein et al., 2009). However, experiments performed in tobacco (Nicotiana tabacum) mesophyll protoplasts have revealed that when expressed alone, AtKC1 is entrapped in the endoplasmic reticulum (ER). However, in tobacco protoplasts and Xenopus laevis oocytes, coexpression of AtKC1 with KAT2 or other inwardly rectifying Shaker subunits (AKT1, KAT1, or AKT2) gives rise to functional heteromeric channels (Duby et al., 2008; Jeanguenin et al., 2011). In Arabidopsis, it is interesting that evidence of the AtKC1 retention in the ER compartment, in the absence of other Shaker subunits, is lacking, since in the native tissues, AtKC1 is always expressed with its inward partners, with which it is able to form heteromeric channels.Here, we took advantage of the unique behavior of AtKC1 when expressed in heterologous systems to investigate the structure-function relationship of the C-linker of KAT2 by sequence exchange between these two channel subunits and by site-directed mutagenesis. The C-linker domain, which, to our knowledge, had never been studied as such in plant Shaker channels before, could be predicted to play crucial roles in channel properties due to its strategic location between the channel transmembrane core and the cytoplasmic CNBD domain. The resulting KAT2-AtKC1 chimeras were expressed in tobacco mesophyll protoplasts and in X. laevis oocytes for investigating their subcellular localization and measuring their activity at the cell membrane. Here, we show that two amino acids present in the C-linker are important for channel subcellular location and that a stretch of nine amino acids forming a short helix just below the membrane, downstream of the sixth transmembrane segment of the channel hydrophobic core, is involved in channel gating. The obtained experimental results are discussed in relation with a KAT2 C-linker/CNBD three-dimensional (3D) model based on animal HCN channels as structure templates.     

Identification of germinal disk region derived genes potentially involved in hen fertility

Despite the regular decrease in fertility observed in hens, especially in “meat” lines, little is known about genes affecting fertility. We used the Affymetrix microarray to search for oocyte genes whose expression would vary in relation to fertility rate in both “laying” and “meat” line hens. We focused on oocyte genes because several of them have been found to be involved in fertility in other species. Based on microarray analysis, 54 and 84 genes were differentially expressed between germinal disc regions (GDR) of F1 maturation stage oocytes from hens exhibiting either high (100%) or low (from 22% to 80%) fertility rate from laying and meat lines respectively. Most of these differentially expressed genes were distributed between “laying” and “meat” lines indicating that mechanisms involved in the decrease in fertility rates in these two cases were independent. Real time RT‐PCR performed on the same samples which were used for microarray confirmed in several cases differences in gene expression levels detected by microarray. Moreover the correlations between gene expression levels and fertility rates were evaluated for the 10 most interesting genes at different stages of follicular maturation and early embryo development on individual GDR samples from hens exhibiting different fertility rates. In total, we identified five genes whose expression levels correlated with fertility rate in accordance with findings of microarray analysis and real time RT‐PCR: VWC2, CR407412, TAPA, FGL2, and TRAP6. The biological significance of these genes sheds light on potential mechanisms influencing fertility and could provide candidates for fertility markers in the hen. Mol. Reprod. Dev. 76: 1043–1055, 2009. © 2009 Wiley‐Liss, Inc.     

Direct Detection of Diverse Metabolic Changes in Virally Transformed and Tax-Expressing Cells by Mass Spectrometry

Background

Viral transformation of a cell starts at the genetic level, followed by changes in the proteome and the metabolome of the host. There is limited information on the broad metabolic changes in HTLV transformed cells.

Methods and Principal Findings

Here, we report the detection of key changes in metabolites and lipids directly from human T-lymphotropic virus type 1 and type 3 (HTLV1 and HTLV3) transformed, as well as Tax1 and Tax3 expressing cell lines by laser ablation electrospray ionization (LAESI) mass spectrometry (MS). Comparing LAESI-MS spectra of non-HTLV1 transformed and HTLV1 transformed cells revealed that glycerophosphocholine (PC) lipid components were dominant in the non-HTLV1 transformed cells, and PC(O-32∶1) and PC(O-34∶1) plasmalogens were displaced by PC(30∶0) and PC(32∶0) species in the HTLV1 transformed cells. In HTLV1 transformed cells, choline, phosphocholine, spermine and glutathione, among others, were downregulated, whereas creatine, dopamine, arginine and AMP were present at higher levels. When comparing metabolite levels between HTLV3 and Tax3 transfected 293T cells, there were a number of common changes observed, including decreased choline, phosphocholine, spermine, homovanillic acid, and glycerophosphocholine and increased spermidine and N-acetyl aspartic acid. These results indicate that the lipid metabolism pathway as well as the creatine and polyamine biosynthesis pathways are commonly deregulated after expression of HTLV3 and Tax3, indicating that the noted changes are likely due to Tax3 expression. N-acetyl aspartic acid is a novel metabolite that is upregulated in all cell types and all conditions tested.

Conclusions and Significance

We demonstrate the high throughput in situ metabolite profiling of HTLV transformed and Tax expressing cells, which facilitates the identification of virus-induced perturbations in the biochemical processes of the host cells. We found virus type-specific (HTLV1 vs. HTLV3), expression-specific (Tax1 vs. Tax3) and cell-type–specific (T lymphocytes vs. kidney epithelial cells) changes in the metabolite profiles. The new insight on the affected metabolic pathways can be used to better understand the molecular mechanisms of HTLV induced transformation, which in turn can result in new treatment strategies.     

Mutations Abrogating VP35 Interaction with Double-Stranded RNA Render Ebola Virus Avirulent in Guinea Pigs

Ebola virus (EBOV) protein VP35 is a double-stranded RNA (dsRNA) binding inhibitor of host interferon (IFN)-α/β responses that also functions as a viral polymerase cofactor. Recent structural studies identified key features, including a central basic patch, required for VP35 dsRNA binding activity. To address the functional significance of these VP35 structural features for EBOV replication and pathogenesis, two point mutations, K319A/R322A, that abrogate VP35 dsRNA binding activity and severely impair its suppression of IFN-α/β production were identified. Solution nuclear magnetic resonance (NMR) spectroscopy and X-ray crystallography reveal minimal structural perturbations in the K319A/R322A VP35 double mutant and suggest that loss of basic charge leads to altered function. Recombinant EBOVs encoding the mutant VP35 exhibit, relative to wild-type VP35 viruses, minimal growth attenuation in IFN-defective Vero cells but severe impairment in IFN-competent cells. In guinea pigs, the VP35 mutant virus revealed a complete loss of virulence. Strikingly, the VP35 mutant virus effectively immunized animals against subsequent wild-type EBOV challenge. These in vivo studies, using recombinant EBOV viruses, combined with the accompanying biochemical and structural analyses directly correlate VP35 dsRNA binding and IFN inhibition functions with viral pathogenesis. Moreover, these studies provide a framework for the development of antivirals targeting this critical EBOV virulence factor.Ebola viruses (EBOVs) are zoonotic, enveloped negative-strand RNA viruses belonging to the family Filoviridae which cause lethal viral hemorrhagic fever in humans and nonhuman primates (47). Currently, information regarding EBOV-encoded virulence determinants remains limited. This, coupled with our lack of understanding of biochemical and structural properties of virulence factors, limits efforts to develop novel prophylactic or therapeutic approaches toward these infections.It has been proposed that EBOV-encoded mechanisms to counter innate immune responses, particularly interferon (IFN) responses, are critical to EBOV pathogenesis (7). However, a role for viral immune evasion functions in the pathogenesis of lethal EBOV infection has yet to be demonstrated. Of the eight major EBOV gene products, two viral proteins have been demonstrated to counter host IFN responses. The VP35 protein is a viral polymerase cofactor and structural protein that also inhibits IFN-α/β production by preventing the activation of interferon regulatory factor (IRF)-3 and -7 (3, 4, 8, 24, 27, 34, 41). VP35 also inhibits the activation of PKR, an IFN-induced, double-stranded RNA (dsRNA)-activated kinase with antiviral activity, and inhibits RNA silencing (17, 20, 48). The VP24 protein is a minor structural protein implicated in virus assembly and regulation of viral RNA synthesis, and changes in VP24 coding sequences are also associated with adaptation of EBOVs to mice and guinea pigs (2, 13, 14, 27, 32, 37, 50, 52). Further, VP24 inhibits cellular responses to both IFN-α/β and IFN-γ by preventing the nuclear accumulation of tyrosine-phosphorylated STAT1 (44, 45). The functions of VP35 and VP24 proteins are manifested in EBOV-infected cells by the absence of IRF-3 activation, impaired production of IFN-α/β, and severely reduced expression of IFN-induced genes, even after treatment of infected cells with IFN-α (3, 19, 21, 22, 24, 25, 28).Previous studies proposed that VP35 basic residues 305, 309, and 312 are required for VP35 dsRNA binding activity (26). VP35 residues K309 and R312 were subsequently identified as critical for binding to dsRNA, and mutation of these residues impaired VP35 suppression of IFN-α/β production (8). In vivo, an EBOV engineered to carry a VP35 R312A point mutation exhibited reduced replication in mice (23). However, because the parental recombinant EBOV into which the mutation was built did not cause disease in these animals, the impact of the mutation on viral pathogenesis could not be fully evaluated. Further, the lack of available structural and biochemical data to explain how the R312A mutation affects VP35 function limited avenues for the therapeutic targeting of critical VP35 functions. Recent structural analyses of the VP35 carboxy-terminal interferon inhibitory domain (IID) suggested that additional residues from the central basic patch may contribute to VP35 dsRNA binding activity and IFN-antagonist function (30). However, a direct correlation between dsRNA and IFN inhibitory functions of VP35 with viral pathogenesis is currently lacking.In order to further define the molecular basis for VP35 dsRNA binding and IFN-antagonist function and to define the contribution of these functions to EBOV pathogenesis, an integrated molecular, structural, and virological approach was taken. The data presented below identify two VP35 carboxy-terminal basic amino acids, K319 and R322, as required for its dsRNA binding and IFN-antagonist functions. Interestingly, these residues are outside the region originally identified as being important for dsRNA binding and IFN inhibition (26). However, they lie within the central basic patch identified by prior structural studies (26, 30). Introduction of these mutations (VP35 with these mutations is designated KRA) into recombinant EBOV renders this otherwise fully lethal virus avirulent in guinea pigs. KRA-infected animals also develop EBOV-specific antibodies and become fully resistant to subsequent challenge with wild-type (WT) virus. Our data further reveal that the KRA EBOV is immunogenic and likely replicates to low levels early after infection in vivo. However, the mutant virus is subsequently cleared by host immune responses. These data demonstrate that the VP35 central basic patch is important not only for IFN-antagonist function but also for EBOV immune evasion and pathogenesis in vivo. High-resolution structural analysis, coupled with our in vitro and in vivo analyses of the recombinant Ebola viruses, provides the molecular basis for loss of function by the VP35 mutant and highlights the therapeutic potential of targeting the central basic patch with small-molecule inhibitors and for future vaccine development efforts.