A web-based resource for the Arabidopsis P450, cytochromes b5, NADPH-cytochrome P450 reductases, and family 1 glycosyltransferases (http://www.P450.kvl.dk)

Genome-wide characterization of cytochrome P450 genes reveals the potential roles in fruit ripening and response to cold stress in tomato

Plant cytochrome P450 (CYP) superfamily, the largest enzyme metabolism family, has been identified in many species and plays a vital role in plant development and stress response via secondary metabolite biosynthesis. A comprehensive identification and functional investigation of CYPs in tomato plants would contribute to deeper understanding of their biological significance. In this study, 268 tomato CYP genes were identified and found to be unevenly located on 12 chromosomes. Based on the phylogenetic analysis, these 268 SlCYPs were classed into two distinct clades (A-type and non-A-type) and nine clans, including 48 families. Moreover, 67 tandem and 22 WGD (whole genome duplication)/segmental duplication events were detected, of which 12 SlCYP genes experienced both WGD/segmental and tandem duplication events, indicating that tandem duplication plays a major role in the expansion of the SlCYP family. Besides, 48 pairs containing 41 SlCYP and 44 AtCYP genes were orthologous, while 216 orthologous pairs were obtained between tomato and potato. The expression level of all SlCYP genes in tomato tissues at different development stages was analyzed, and most expressed SlCYPs showed a tissue-specific pattern. Meanwhile, 143 differentially expressed SlCYPs were identified under cold stress. Furthermore, the RT-qPCR results indicated that SlCYPs may be involved in fruit ripening and cold tolerance in tomato seedlings. These findings provide valuable insights into the evolutionary relationships and functional characteristics of SlCYPs, which can be utilized for further investigation of fruit metabolic pathways and cold tolerance in tomato.

Genome-wide investigation of UDP-Glycosyltransferase family in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum) belongs to Polygonaceae family and has attracted increasing attention owing to its high nutritional value. UDP-glycosyltransferases (UGTs) glycosylate a variety of plant secondary metabolites to control many metabolic processes during plant growth and development. However, there have been no systematic reports of UGT superfamily in F. tataricum.

RESULTS

We identified 173 FtUGTs in F. tataricum based on their conserved UDPGT domain. Phylogenetic analysis of FtUGTs with 73 Arabidopsis UGTs clustered them into 21 families. FtUGTs from the same family usually had similar gene structure and motif compositions. Most of FtUGTs did not contain introns or had only one intron. Tandem repeats contributed more to FtUGTs amplification than segmental duplications. Expression analysis indicates that FtUGTs are widely expressed in various tissues and likely play important roles in plant growth and development. The gene expression analysis response to different abiotic stresses showed that some FtUGTs were involved in response to drought and cadmium stress. Our study provides useful information on the UGTs in F. tataricum, and will facilitate their further study to better understand their function.

CONCLUSIONS

Our results provide a theoretical basis for further exploration of the functional characteristics of FtUGTs and for understanding the growth, development, and metabolic model in F. tataricum.

Phylogenomic Analysis of Cytochrome P450 Gene Superfamily and Their Association with Flavonoids Biosynthesis in Peanut (Arachis hypogaea L.)

Cytochrome P450s (CYPs) constitute extensive enzyme superfamilies in the plants, playing pivotal roles in a multitude of biosynthetic and detoxification pathways essential for growth and development, such as the flavonoid biosynthesis pathway. However, CYPs have not yet been systematically studied in the cultivated peanuts (Arachis hypogaea L.), a globally significant cash crop. This study addresses this knowledge deficit through a comprehensive genome-wide analysis, leading to the identification of 589 AhCYP genes in peanuts. Through phylogenetic analysis, all AhCYPs were systematically classified into 9 clans, 43 gene families. The variability in the number of gene family members suggests specialization in biological functions. Intriguingly, both tandem duplication and fragment duplication events have emerged as pivotal drivers in the evolutionary expansion of the AhCYP superfamily. Ka/Ks analysis underscored the substantial influence of strong purifying selection on the evolution of AhCYPs. Furthermore, we selected 21 genes encoding 8 enzymes associated with the flavonoid pathway. The results of quantitative real-time PCR (qRT-PCR) experiments unveiled stage-specific expression patterns during the development of peanut testa, with discernible variations between pink and red testa. Importantly, we identified a direct correlation between gene expression levels and the accumulation of metabolites. These findings offer valuable insights into elucidating the comprehensive functions of AhCYPs and the underlying mechanisms governing the divergent accumulation of flavonoids in testa of different colors.

Genome-wide identification, comparative analysis and functional roles in flavonoid biosynthesis of cytochrome P450 superfamily in pear (Pyrus spp.)

BACKGROUND

The cytochrome P450 (CYP) superfamily is the largest enzyme metabolism family in plants identified to date, and it is involved in many biological processes, including secondary metabolite biosynthesis, hormone metabolism and stress resistance. However, the P450 gene superfamily has not been well studied in pear (Pyrus spp.).

RESULTS

Here, the comprehensive identification and a comparative analysis of P450 superfamily members were conducted in cultivated and wild pear genomes. In total, 338, 299 and 419 P450 genes were identified in Chinese white pear, European pear and the wild pear, respectively. Based on the phylogenetic analyses, pear P450 genes were divided into ten clans, comprising 48 families. The motif and gene structure analyses further supported this classification. The expansion of the pear P450 gene family was attributed to whole-genome and single-gene duplication events. Several P450 gene clusters were detected, which have resulted from tandem and proximal duplications. Purifying selection was the major force imposed on the long-term evolution of P450 genes. Gene dosage balance, subfunctionalization and neofunctionalization jointly drove the retention and functional diversification of P450 gene pairs. Based on the association analysis between transcriptome expression profiles and flavonoid content during fruit development, three candidate genes were identified as being closely associated with the flavonoid biosynthesis, and the expression of one gene was further verified using qRT-PCR and its function was validated through transient transformation in pear fruit.

CONCLUSIONS

The study results provide insights into the evolution and biological functions of P450 genes in pear.

Cytochrome P450 Gene Families: Role in Plant Secondary Metabolites Production and Plant Defense

Cytochrome P450s (CYPs) are the most prominent family of enzymes involved in NADPH- and O2-dependent hydroxylation processes throughout all spheres of life. CYPs are crucial for the detoxification of xenobiotics in plants, insects, and other organisms. In addition to performing this function, CYPs serve as flexible catalysts and are essential for producing secondary metabolites, antioxidants, and phytohormones in higher plants. Numerous biotic and abiotic stresses frequently affect the growth and development of plants. They cause a dramatic decrease in crop yield and a deterioration in crop quality. Plants protect themselves against these stresses through different mechanisms, which are accomplished by the active participation of CYPs in several biosynthetic and detoxifying pathways. There are immense potentialities for using CYPs as a candidate for developing agricultural crop species resistant to biotic and abiotic stressors. This review provides an overview of the plant CYP families and their functions to plant secondary metabolite production and defense against different biotic and abiotic stresses.

Molecular Characterization of an Isoflavone 2'-Hydroxylase Gene Revealed Positive Insights into Flavonoid Accumulation and Abiotic Stress Tolerance in Safflower

Flavonoids with significant therapeutic properties play an essential role in plant growth, development, and adaptation to various environments. The biosynthetic pathway of flavonoids has long been studied in plants; however, its regulatory mechanism in safflower largely remains unclear. Here, we carried out comprehensive genome-wide identification and functional characterization of a putative cytochrome P45081E8 gene encoding an isoflavone 2'-hydroxylase from safflower. A total of 15 CtCYP81E genes were identified from the safflower genome. Phylogenetic classification and conserved topology of CtCYP81E gene structures, protein motifs, and cis-elements elucidated crucial insights into plant growth, development, and stress responses. The diverse expression pattern of CtCYP81E genes in four different flowering stages suggested important clues into the regulation of secondary metabolites. Similarly, the variable expression of CtCYP81E8 during multiple flowering stages further highlighted a strong relationship with metabolite accumulation. Furthermore, the orchestrated link between transcriptional regulation of CtCYP81E8 and flavonoid accumulation was further validated in the yellow- and red-type safflower. The spatiotemporal expression of CtCYP81E8 under methyl jasmonate, polyethylene glycol, light, and dark conditions further highlighted its likely significance in abiotic stress adaption. Moreover, the over-expressed transgenic Arabidopsis lines showed enhanced transcript abundance in OE-13 line with approximately eight-fold increased expression. The upregulation of AtCHS, AtF3'H, and AtDFR genes and the detection of several types of flavonoids in the OE-13 transgenic line also provides crucial insights into the potential role of CtCYP81E8 during flavonoid accumulation. Together, our findings shed light on the fundamental role of CtCYP81E8 encoding a putative isoflavone 2'-hydroxylase via constitutive expression during flavonoid biosynthesis.

Functional and structural insight into the flexibility of cytochrome P450 reductases from Sorghum bicolor and its implications for lignin composition

Plant NADPH-dependent cytochrome P450 reductase (CPR) is a multidomain enzyme that donates electrons for hydroxylation reactions catalyzed by class II cytochrome P450 monooxygenases involved in the synthesis of many primary and secondary metabolites. These P450 enzymes include trans-cinnamate-4-hydroxylase, p-coumarate-3′-hydroxylase, and ferulate-5-hydroxylase involved in monolignol biosynthesis. Because of its role in monolignol biosynthesis, alterations in CPR activity could change the composition and overall output of lignin. Therefore, to understand the structure and function of three CPR subunits from sorghum, recombinant subunits SbCPR2a, SbCPR2b, and SbCPR2c were subjected to X-ray crystallography and kinetic assays. Steady-state kinetic analyses demonstrated that all three CPR subunits supported the oxidation reactions catalyzed by SbC4H1 (CYP73A33) and SbC3′H (CYP98A1). Furthermore, comparing the SbCPR2b structure with the well-investigated CPRs from mammals enabled us to identify critical residues of functional importance and suggested that the plant flavin mononucleotide–binding domain might be more flexible than mammalian homologs. In addition, the elucidated structure of SbCPR2b included the first observation of NADP⁺ in a native CPR. Overall, we conclude that the connecting domain of SbCPR2, especially its hinge region, could serve as a target to alter biomass composition in bioenergy and forage sorghums through protein engineering.

Metabolic Pathways for S-Metolachlor Detoxification Differ Between Tolerant Corn and Multiple-Resistant Waterhemp

Herbicide resistance in weeds can be conferred by target-site and/or non-target-site mechanisms, such as rapid metabolic detoxification. Resistance to the very-long-chain fatty acid-inhibiting herbicide, S-metolachlor, in multiple herbicide-resistant populations (CHR and SIR) of waterhemp (Amaranthus tuberculatus) is conferred by rapid metabolism compared with sensitive populations. However, enzymatic pathways for S-metolachlor metabolism in waterhemp are unknown. Enzyme assays using S-metolachlor were developed to determine the specific activities of glutathione S-transferases (GSTs) and cytochrome P450 monooxygenases (P450s) from CHR and SIR seedlings to compare with tolerant corn and sensitive waterhemp (WUS). GST activities were greater (∼2-fold) in CHR and SIR compared to WUS but much less than corn. In contrast, P450s in microsomal extracts from CHR and SIR formed O-demethylated S-metolachlor, and their NADPH-dependent specific activities were greater (>20-fold) than corn or WUS. Metabolite profiles of S-metolachlor generated via untargeted and targeted liquid chromatography-mass spectrometry from CHR and SIR differed from WUS, with greater relative abundances of O-demethylated S-metolachlor and O-demethylated S-metolachlor-glutathione conjugates formed by CHR and SIR. In summary, our results demonstrate that S-metolachlor metabolism in resistant waterhemp involves Phase I and Phase II metabolic activities acting in concert, but the initial O-demethylation reaction confers resistance.

Resistance to a nonselective 4-hydroxyphenylpyruvate dioxygenase-inhibiting herbicide via novel reduction-dehydration-glutathione conjugation in Amaranthus tuberculatus

Metabolic resistance to 4-hydroxyphenylpyruvate dioxygenase (HPPD)-inhibiting herbicides is a threat in controlling waterhemp (Amaranthus tuberculatus) in the USA. We investigated resistance mechanisms to syncarpic acid-3 (SA3), a nonselective, noncommercial HPPD-inhibiting herbicide metabolically robust to Phase I oxidation, in multiple-herbicide-resistant (MHR) waterhemp populations (SIR and NEB) and HPPD inhibitor-sensitive populations (ACR and SEN). Dose-response experiments with SA3 provided ED₅₀ -based resistant : sensitive ratios of at least 18-fold. Metabolism experiments quantifying parent SA3 remaining in excised leaves during a time course indicated MHR populations displayed faster rates of SA3 metabolism compared to HPPD inhibitor-sensitive populations. SA3 metabolites were identified via LC-MS-based untargeted metabolomics in whole plants. A Phase I metabolite, likely generated by cytochrome P450-mediated alkyl hydroxylation, was detected but was not associated with resistance. A Phase I metabolite consistent with ketone reduction followed by water elimination was detected, creating a putative α,β-unsaturated carbonyl resembling a Michael acceptor site. A Phase II glutathione-SA3 conjugate was associated with resistance. Our results revealed a novel reduction-dehydration-GSH conjugation detoxification mechanism. SA3 metabolism in MHR waterhemp is thus atypical compared to commercial HPPD-inhibiting herbicides. This previously uncharacterized detoxification mechanism presents a unique opportunity for future biorational design by blocking known sites of herbicide metabolism in weeds.

An independent evolutionary origin for insect deterrent cucurbitacins in Iberis amara

Pieris rapae and Phyllotreta nemorum are Brassicaceae specialists, but do not feed on Iberis amara spp. that contain cucurbitacins. The cucurbitacins are highly oxygenated triterpenoid, occurring widespread in cucurbitaceous species and in a few other plant families. Using de-novo assembled transcriptomics from I. amara, gene co-expression analysis and comparative genomics, we unraveled the evolutionary origin of the insect deterrent cucurbitacins in I. amara. Phylogenetic analysis of five oxidosqualene cyclases and heterologous expression allowed us to identify the first committed enzyme in cucurbitacin biosynthesis in I. amara, cucurbitadienol synthase (IaCPQ). In addition, two species-specific cytochrome P450s (CYP708A16 and CYP708A15) were identified that catalyse the unique C16 and C22 hydroxylation of the cucurbitadienol backbone, enzymatic steps that have not been reported before. Furthermore, the draft genome assembly of I. amara showed that the IaCPQ was localized to the same scaffold together with CYP708A15 but spanning over 100 kb, this contrasts with the highly organized cucurbitacin gene cluster in the cucurbits. These results reveal that cucurbitacin biosynthesis has evolved convergently via different biosynthetic routes in different families rather than through divergence from an ancestral pathway. This study thus provides new insight into the mechanism of recurrent evolution and diversification of a plant defensive chemical.

Leveraging computational genomics to understand the molecular basis of metal homeostasis

Genome-based data is helping to reveal the diverse strategies plants and algae use to maintain metal homeostasis. In addition to acquisition, distribution and storage of metals, acclimating to feast or famine can involve a wealth of genes that we are just now starting to understand. The fast-paced acquisition of genome-based data, however, is far outpacing our ability to experimentally characterize protein function. Computational genomic approaches are needed to fill the gap between what is known and unknown. To avoid misconstruing bioinformatically derived data, which is the root cause of the inaccurate functional annotations that plague databases, functional inferences from diverse sources and contextualization of that evidence with a robust understanding of protein family evolution is needed. Phylogenomic- and comparative-genomic-based studies can aid in the interpretation of experimental data or provide a spark for the discovery of a new function. These analyses not only lead to novel insight into a target protein's function but can generate thought-provoking insights across protein families.

Diversity and evolution of cytochrome P450s of Jacobaea vulgaris and Jacobaea aquatica

BACKGROUND

Collectively, plants produce a huge variety of secondary metabolites (SMs) which are involved in the adaptation of plants to biotic and abiotic stresses. The most characteristic feature of SMs is their striking inter- and intraspecific chemical diversity. Cytochrome P450 monooxygenases (CYPs) often play an important role in the biosynthesis of SMs and thus in the evolution of chemical diversity. Here we studied the diversity and evolution of CYPs of two Jacobaea species which contain a characteristic group of SMs namely the pyrrolizidine alkaloids (PAs).

RESULTS

We retrieved CYPs from RNA-seq data of J. vulgaris and J. aquatica, resulting in 221 and 157 full-length CYP genes, respectively. The analyses of conserved motifs confirmed that Jacobaea CYP proteins share conserved motifs including the heme-binding signature, the PERF motif, the K-helix and the I-helix. KEGG annotation revealed that the CYPs assigned as being SM metabolic pathway genes were all from the CYP71 clan but no CYPs were assigned as being involved in alkaloid pathways. Phylogenetic analyses of full-length CYPs were conducted for the six largest CYP families of Jacobaea (CYP71, CYP76, CYP706, CYP82, CYP93 and CYP72) and were compared with CYPs of two other members of the Asteraceae, Helianthus annuus and Lactuca sativa, and with Arabidopsis thaliana. The phylogenetic trees showed strong lineage specific diversification of CYPs, implying that the evolution of CYPs has been very fast even within the Asteraceae family. Only in the closely related species J. vulgaris and J. aquatica, CYPs were found often in pairs, confirming a close relationship in the evolutionary history.

CONCLUSIONS

This study discovered 378 full-length CYPs in Jacobaea species, which can be used for future exploration of their functions, including possible involvement in PA biosynthesis and PA diversity.

Discovery and modification of cytochrome P450 for plant natural products biosynthesis

Cytochrome P450s are widespread in nature and play key roles in the diversification and functional modification of plant natural products. Over the last few years, there has been remarkable progress in plant P450s identification with the rapid development of sequencing technology, "omics" analysis and synthetic biology. However, challenges still persist in respect of crystal structure, heterologous expression and enzyme engineering. Here, we reviewed several research hotspots of P450 enzymes involved in the biosynthesis of plant natural products, including P450 databases, gene mining, heterologous expression and protein engineering.

Role of Cytochrome P450 Enzymes in Plant Stress Response

Cytochrome P450s (CYPs) are the largest enzyme family involved in NADPH- and/or O₂-dependent hydroxylation reactions across all the domains of life. In plants and animals, CYPs play a central role in the detoxification of xenobiotics. In addition to this function, CYPs act as versatile catalysts and play a crucial role in the biosynthesis of secondary metabolites, antioxidants, and phytohormones in higher plants. The molecular and biochemical processes catalyzed by CYPs have been well characterized, however, the relationship between the biochemical process catalyzed by CYPs and its effect on several plant functions was not well established. The advent of next-generation sequencing opened new avenues to unravel the involvement of CYPs in several plant functions such as plant stress response. The expression of several CYP genes are regulated in response to environmental stresses, and they also play a prominent role in the crosstalk between abiotic and biotic stress responses. CYPs have an enormous potential to be used as a candidate for engineering crop species resilient to biotic and abiotic stresses. The objective of this review is to summarize the latest research on the role of CYPs in plant stress response.

Genome-wide characterization of the UDP-glycosyltransferase gene family in upland cotton

Uridine diphosphate (UDP)-glycosyltransferases (UGTs) involved in many metabolic processes are ubiquitous in plants, animals, microorganisms and other organisms and are essential for their growth and development. Upland cotton contains a large number of UGT genes. In this study, we aimed to identify UGT family members in the genome of upland cotton (Gossypium hirsutum L.) and analyze their expression patterns. Bioinformatics methods were used to identify UGT genes from the whole genome of upland cotton (Gossypium hirsutum L. acc. TM-1). Phylogenetic analysis was conducted based on alignment of UGT proteins from upland cotton, and the gene structure, motif and chromosome localization were analyzed for the H subgroup of the UGT family. And the physical and chemical properties and expressions of the genes in the H subgroup of this family were also analyzed. A total of 274 UGT genes were identified from the whole genome of upland cotton and were divided into nine subgroups based on phylogenetic analyses. In subgroup H, 36 genes were distributed on 18 chromosomes. The subfamily genes were simple in the structure, 19 of its members contained two introns, and the others contained only one intron. The qRT-PCR results and transcriptomic data indicated that most of the genes had a wide range of tissue expression characteristics. And the phylogenetic analysis results and expression profiles of these genes revealed tissues and different UGT genes from this crop. Taking RNA-seq, RT-qPCR, and quantitative trait locus (QTL) mapping together, our results suggested that GhUGT6 and GhUGT105 in subgroup H of the GhUGT gene family could be potential candidate genes for cotton yield, and GhUGT16, GhUGT103 might play a vital role in fiber development.

Transcriptome analyses of Paris polyphylla var. chinensis, Ypsilandra thibetica, and Polygonatum kingianum characterize their steroidal saponin biosynthesis pathway

Steroidal saponins, one of the most diverse groups of plant-derived natural products, elicit biological and pharmacological activities; however, the genes involved in their biosynthesis and the corresponding biosynthetic pathway in monocotyledon plants remain unclear. This study aimed to identify genes involved in the biosynthesis of steroidal saponins by performing a comparative analysis among transcriptomes of Paris polyphylla var. chinensis (PPC), Ypsilandra thibetica (YT), and Polygonatum kingianum (PK). De novo transcriptome assemblies generated 57,537, 140,420, and 151,773 unigenes from PPC, YT, and PK, respectively, of which 56.54, 47.81, and 44.30% were successfully annotated, respectively. Among the transcriptomes for PPC, YT, and PK, we identified 194, 169, and 131; 17, 14, and 26; and, 80, 122, and 113 unigenes corresponding to terpenoid backbone biosynthesis; sesquiterpenoid and triterpenoid biosynthesis; and, steroid biosynthesis pathways, respectively. These genes are putatively involved in the biosynthesis of cholesterol that is the primary precursor of steroidal saponins. Phylogenetic analyses indicated that lanosterol synthase may be exclusive to dicotyledon plant species, and the cytochrome P450 unigenes were closely related to clusters CYP90B1 and CYP734A1, which are UDP-glycosyltransferases unigenes homologous with the UGT73 family. Thus, unigenes of β-glucosidase may be candidate genes for catalysis of later period modifications of the steroidal saponin skeleton. Our data provide evidence to support the hypothesis that monocotyledons biosynthesize steroidal saponins from cholesterol via the cycloartenol pathway.

Genome-wide analysis of family-1 UDP glycosyltransferases (UGT) and identification of UGT genes for FHB resistance in wheat (Triticum aestivum L.)

BACKGROUND

Fusarium head blight (FHB), a devastating disease in wheat worldwide, results in yield loses and mycotoxin, such as deoxynivalenol (DON), accumulation in infected grains. DON also facilitates the pathogen colonization and spread of FHB symptoms during disease development. UDP-glycosyltransferase enzymes (UGTs) are known to contribute to detoxification and enhance FHB resistance by glycosylating DON into DON-3-glucoside (D3G) in wheat. However, a comprehensive investigation of wheat (Triticum aestivum) UGT genes is still lacking.

RESULTS

In this study, we carried out a genome-wide analysis of family-1 UDP glycosyltransferases in wheat based on the PSPG conserved box that resulted in the identification of 179 putative UGT genes. The identified genes were clustered into 16 major phylogenetic groups with a lack of phylogenetic group K. The UGT genes were invariably distributed among all the chromosomes of the 3 genomes. At least 10 intron insertion events were found in the UGT sequences, where intron 4 was observed as the most conserved intron. The expression analysis of the wheat UGT genes using both online microarray data and quantitative real-time PCR verification suggested the distinct role of UGT genes in different tissues and developmental stages. The expression of many UGT genes was up-regulated after Fusarium graminearum inoculation, and six of the genes were further verified by RT-qPCR.

CONCLUSION

We identified 179 UGT genes from wheat using the available sequenced wheat genome. This study provides useful insight into the phylogenetic structure, distribution, and expression patterns of family-1 UDP glycosyltransferases in wheat. The results also offer a foundation for future work aimed at elucidating the molecular mechanisms underlying the resistance to FHB and DON accumulation.

Identification and analysis of CYP450 genes from transcriptome of Lonicera japonica and expression analysis of chlorogenic acid biosynthesis related CYP450s

Background

Lonicera japonica is an important medicinal plant that has been widely used in traditional Chinese medicine for thousands of years. The pharmacological activities of L. japonica are mainly due to its rich natural active ingredients, most of which are secondary metabolites. CYP450s are a large, complex, and widespread superfamily of proteins that participate in many endogenous and exogenous metabolic reactions, especially secondary metabolism. Here, we identified CYP450s in L. japonica transcriptome and analyzed CYP450s that may be involved in chlorogenic acid (CGA) biosynthesis.

Methods

The recent availability of L. japonica transcriptome provided opportunity to identify CYP450s in this herb. BLAST based method and HMM based method were used to identify CYP450s in L. japonica transcriptome. Then, phylogenetic analysis, conserved motifs analysis, GO annotation, and KEGG annotation analyses were conducted to characterize the identified CYP450s. qRT-PCR was used to explore expression patterns of five CGA biosynthesis related CYP450s.

Results

In this study, 151 putative CYP450s with complete cytochrome P450 domain, which belonged to 10 clans, 45 families and 76 subfamilies, were identified in L. japonica transcriptome. Phylogenetic analysis classified these CYP450s into two major branches, A-type (47%) and non-A type (53%). Both types of CYP450s had conserved motifs in L. japonica. The differences of typical motif sequences between A-type and non-A type CYP450s in L. japonica were similar with other plants. GO classification indicated that non-A type CYP450s participated in more molecular functions and biological processes than A-type. KEGG pathway annotation totally assigned 47 CYP450s to 25 KEGG pathways. From these data, we cloned two LjC3Hs (CYP98A subfamily) and three LjC4Hs (CYP73A subfamily) that may be involved in biosynthesis of CGA, the major ingredient for pharmacological activities of L. japonica. qRT-PCR results indicated that two LjC3Hs exhibited oppositing expression patterns during the flower development and LjC3H2 exhibited a similar expression pattern with CGA concentration measured by HPLC. The expression patterns of three LjC4Hs were quite different and the expression pattern of LjC4H3 was quite similar with that of LjC3H1.

Discussion

Our results provide a comprehensive identification and characterization of CYP450s in L. japonica. Five CGA biosynthesis related CYP450s were cloned and their expression patterns were explored. The different expression patterns of two LjC3Hs and three LjC4Hs may be due to functional divergence of both substrate and catalytic specificity during plant evolution. The co-expression pattern of LjC3H1 and LjC4H3 strongly suggested that they were under coordinated regulation by the same transcription factors due to same cis elements in their promoters. In conclusion, this study provides insight into CYP450s and will effectively facilitate the research of biosynthesis of CGA in L. japonica.

High-throughput sequencing and de novo transcriptome assembly of Swertia japonica to identify genes involved in the biosynthesis of therapeutic metabolites

Here, we report potential transcripts involved in the biosynthesis of therapeutic metabolites in Swertia japonica , the first report of transcriptome assembly, and characterization of the medicinal plant from Swertia genus. Swertia genus, representing over 170 plant species including herbs such as S. chirata, S. hookeri, S. longifolia, S. japonica, among others, have been used as the traditional medicine in China, India, Korea, and Japan for thousands of years. Due to the lack of genomic and transcriptomic resources, little is known about the molecular basis involved in the biosynthesis of characteristic key bioactive metabolites. Here, we performed deep-transcriptome sequencing for the aerial tissues and the roots of S. japonica, generating over 2 billion raw reads with an average length of 101 bps. Using a combined approach of three popular assemblers, de novo transcriptome assembly for S. japonica was obtained, yielding 81,729 unigenes having an average length of 884 bps and N50 value of 1452 bps, of which 46,963 unigenes were annotated based on the sequence similarity against NCBI-nr protein database. Annotation of transcriptome assembly resulted in the identification of putative genes encoding all enzymes from the key therapeutic metabolite biosynthesis pathways. Transcript abundance analysis, gene ontology enrichment analysis, and KEGG pathway enrichment analysis revealed metabolic processes being up-regulated in the aerial tissues with respect to the roots of S. japonica. We also identified 37 unigenes as potential candidates involved in the glycosylation of bioactive metabolites. Being the first report of transcriptome assembly and annotation for any of the Swertia species, this study will be a valuable resource for future investigations on the biosynthetic pathways of therapeutic metabolites and their regulations.

RNA-seq Transcriptome Analysis of Panax japonicus, and Its Comparison with Other Panax Species to Identify Potential Genes Involved in the Saponins Biosynthesis

The Panax genus has been a source of natural medicine, benefitting human health over the ages, among which the Panax japonicus represents an important species. Our understanding of several key pathways and enzymes involved in the biosynthesis of ginsenosides, a pharmacologically active class of metabolites and a major chemical constituents of the rhizome extracts from the Panax species, are limited. Limited genomic information, and lack of studies on comparative transcriptomics across the Panax species have restricted our understanding of the biosynthetic mechanisms of these and many other important classes of phytochemicals. Herein, we describe Illumina based RNA sequencing analysis to characterize the transcriptome and expression profiles of genes expressed in the five tissues of P. japonicus, and its comparison with other Panax species. RNA sequencing and de novo transcriptome assembly for P. japonicus resulted in a total of 135,235 unigenes with 78,794 (58.24%) unigenes being annotated using NCBI-nr database. Transcriptome profiling, and gene ontology enrichment analysis for five tissues of P. japonicus showed that although overall processes were evenly conserved across all tissues. However, each tissue was characterized by several unique unigenes with the leaves showing the most unique unigenes among the tissues studied. A comparative analysis of the P. japonicus transcriptome assembly with publically available transcripts from other Panax species, namely, P. ginseng, P. notoginseng, and P. quinquefolius also displayed high sequence similarity across all Panax species, with P. japonicus showing highest similarity with P. ginseng. Annotation of P. japonicus transcriptome resulted in the identification of putative genes encoding all enzymes from the triterpene backbone biosynthetic pathways, and identified 24 and 48 unigenes annotated as cytochrome P450 (CYP) and glycosyltransferases (GT), respectively. These CYPs and GTs annotated unigenes were conserved across all Panax species and co-expressed with other the transcripts involved in the triterpenoid backbone biosynthesis pathways. Unigenes identified in this study represent strong candidates for being involved in the triterpenoid saponins biosynthesis, and can serve as a basis for future validation studies.

Evolution of NADPH-cytochrome P450 oxidoreductases (POR) in Apiales - POR 1 is missing

The NADPH-dependent cytochrome P450 oxidoreductase (POR) is the obligate electron donor to eukaryotic microsomal cytochromes P450 enzymes. The number of PORs within plant species is limited to one to four isoforms, with the most common being two PORs per plant. These enzymes provide electrons to a huge number of different cytochromes P450s (from 50 to several hundred within one plant). Within the eudicotyledons, PORs can be divided into two major clades, POR 1 and POR 2. Based on our own sequencing analysis and publicly available data, we have identified 45 PORs from the angiosperm order Apiales. These were subjected to a phylogenetic analysis along with 237 other publicly available (NCBI and oneKP) POR sequences found within the clade Asterids. Here, we show that the order Apiales only harbor members of the POR 2 clade, which are further divided into two distinct subclades. This is in contrast to most other eudicotyledon orders that have both POR 1 and POR 2. This suggests that through gene duplications and one gene deletion, Apiales only contain members of the POR 2 clade. Three POR 2 isoforms from Thapsia garganica L., Apiaceae, were all full-length in an Illumina root transcriptome dataset (available from the SRA at NCBI). All three genes were shown to be functional upon reconstitution into nanodiscs, confirming that none of the isoforms are pseudogenes.

Expanded investigations of the aglycon promiscuity and catalysis characteristic of flavonol 3-O-rhamnosyltransferase AtUGT78D1 from Arabidopsis thaliana

Rhamnosides usually possess better bioavailabilities and improved solubilities compared with their aglycons and are a major source of bioactive natural products.

Comparative interactions of withanolides and sterols with two members of sterol glycosyltransferases from Withania somnifera

Background

Sterol glycosyltransferases (SGTs) are ubiquitous but one of the most diverse group of enzymes of glycosyltransferases family. Members of this family modulate physical and chemical properties of secondary plant products important for various physiological processes. The role of SGTs has been demonstrated in the biosynthesis of pharmaceutically important molecules of medicinal plants like Withania somnifera.

Results

Analysis suggested conserved behaviour and high similarity in active sites of WsSGTs with other plant GTs. Substrate specificity of WsSGTs were analysed through docking performance of WsSGTs with different substrates (sterols and withanolides). Best docking results of WsSGTL1 in the form of stable enzyme-substrate complex having lowest binding energies were obtained with brassicasterol, transandrosteron and WsSGTL4 with solasodine, stigmasterol and 24-methylene cholesterol.

Conclusion

This study reveals topological characters and conserved nature of two SGTs from W. somnifera (WsSGTs) i.e. WsSGTL1 and WsSGTL4. However, besides being ubiquitous in nature and with broad substrate specificity, difference between WsSGTL1 and WsSGTL4 is briefly described by difference in stability (binding energy) of enzyme-substrate complexes through comparative docking.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0563-7) contains supplementary material, which is available to authorized users.

Computational Identification and Systematic Classification of Novel Cytochrome P450 Genes in Salvia miltiorrhiza

Salvia miltiorrhiza is one of the most economically important medicinal plants. Cytochrome P450 (CYP450) genes have been implicated in the biosynthesis of its active components. However, only a dozen full-length CYP450 genes have been described, and there is no systematic classification of CYP450 genes in S. miltiorrhiza. We obtained 77,549 unigenes from three tissue types of S. miltiorrhiza using RNA-Seq technology. Combining our data with previously identified CYP450 sequences and scanning with the CYP450 model from Pfam resulted in the identification of 116 full-length and 135 partial-length CYP450 genes. The 116 genes were classified into 9 clans and 38 families using standard criteria. The RNA-Seq results showed that 35 CYP450 genes were co-expressed with CYP76AH1, a marker gene for tanshinone biosynthesis, using r≥0.9 as a cutoff. The expression profiles for 16 of 19 randomly selected CYP450 obtained from RNA-Seq were validated by qRT-PCR. Comparing against the KEGG database, 10 CYP450 genes were found to be associated with diterpenoid biosynthesis. Considering all the evidence, 3 CYP450 genes were identified to be potentially involved in terpenoid biosynthesis. Moreover, we found that 15 CYP450 genes were possibly regulated by antisense transcripts (r≥0.9 or r≤–0.9). Lastly, a web resource (SMCYP450, http://www.herbalgenomics.org/samicyp450) was set up, which allows users to browse, search, retrieve and compare CYP450 genes and can serve as a centralized resource.

Genome-wide identification and expression analyses of cytochrome P450 genes in mulberry (Morus notabilis)

Cytochrome P450s play critical roles in the biosynthesis of physiologically important compounds in plants. These compounds often act as defense toxins to prevent herbivory. In the present study, a total of 174 P450 genes of mulberry (Morus notabilis C.K.Schn) were identified based on bioinformatics analyses. These mulberry P450 genes were divided into nine clans and 47 families and were found to be expressed in a tissue-preferential manner. These genes were compared to the P450 genes in Arabidopsis thaliana. Families CYP80, CYP92, CYP728, CYP733, CYP736, and CYP749 were found to exist in mulberry, and they may play important roles in the biosynthesis of mulberry secondary metabolites. Analyses of the functional and metabolic pathways of these genes indicated that mulberry P450 genes may participate in the metabolism of lipids, other secondary metabolites, xenobiotics, amino acids, cofactors, vitamins, terpenoids, and polyketides. These results provide a foundation for understanding of the structures and biological functions of mulberry P450 genes.

Exploiting natural variation of secondary metabolism identifies a gene controlling the glycosylation diversity of dihydroxybenzoic acids in Arabidopsis thaliana

Plant secondary metabolism is an active research area because of the unique and important roles the specialized metabolites have in the interaction of plants with their biotic and abiotic environment, the diversity and complexity of the compounds and their importance to human medicine. Thousands of natural accessions of Arabidopsis thaliana characterized with increasing genomic precision are available, providing new opportunities to explore the biochemical and genetic mechanisms affecting variation in secondary metabolism within this model species. In this study, we focused on four aromatic metabolites that were differentially accumulated among 96 Arabidopsis natural accessions as revealed by leaf metabolic profiling. Using UV, mass spectrometry, and NMR data, we identified these four compounds as different dihydroxybenzoic acid (DHBA) glycosides, namely 2,5-dihydroxybenzoic acid (gentisic acid) 5-O-β-D-glucoside, 2,3-dihydroxybenzoic acid 3-O-β-D-glucoside, 2,5-dihydroxybenzoic acid 5-O-β-D-xyloside, and 2,3-dihydroxybenzoic acid 3-O-β-D-xyloside. Quantitative trait locus (QTL) mapping using recombinant inbred lines generated from C24 and Col-0 revealed a major-effect QTL controlling the relative proportion of xylosides vs. glucosides. Association mapping identified markers linked to a gene encoding a UDP glycosyltransferase gene. Analysis of Transfer DNA (T-DNA) knockout lines verified that this gene is required for DHBA xylosylation in planta and recombinant protein was able to xylosylate DHBA in vitro. This study demonstrates that exploiting natural variation of secondary metabolism is a powerful approach for gene function discovery.

Cytochrome P450 genes induced by bispyribac-sodium treatment in a multiple-herbicide-resistant biotype of Echinochloa phyllopogon

BACKGROUND

Incremental herbicide metabolism by cytochrome P450 monooxygenases (P450s) has been proposed as the basis for resistance to bispyribac-sodium (bispyribac) in a multiple-herbicide-resistant biotype of Echinochloa phyllopogon. Upon exposure to bispyribac, strong induction of bispyribac-metabolising P450 activity has been reported in the resistant line, indicating that P450s induced by bispyribac are involved in the bispyribac resistance.

RESULTS

A polymerase chain reaction (PCR)-based cloning strategy was used to isolate 39 putative P450 genes from the bispyribac-resistant line of E. phyllopogon. Expression analysis by real-time PCR revealed that seven of the isolated genes were upregulated in response to bispyribac treatment of seedlings at the three-leaf stage. The transcript levels and protein sequences of the seven genes were compared between the bispyribac-resistant line and a susceptible line. CYP71AK2 and CYP72A254 were transcribed prominently in the bispyribac-resistant line. Amino acid polymorphisms were found in three genes, including CYP72A254.

CONCLUSION

Upregulated expression of these genes is consistent with the inducible herbicide-metabolising P450 activity under bispyribac stress that was reported in a previous study. This is the first study to compare P450 genes in arable weed species in order to elucidate the mechanism for P450-mediated herbicide resistance.

New insights into the targeting of a subset of tail-anchored proteins to the outer mitochondrial membrane

Tail-anchored (TA) proteins are a unique class of functionally diverse membrane proteins defined by their single C-terminal membrane-spanning domain and their ability to insert post-translationally into specific organelles with an Ncytoplasm-Corganelle interior orientation. The molecular mechanisms by which TA proteins are sorted to the proper organelles are not well-understood. Herein we present results indicating that a dibasic targeting motif (i.e., -R-R/K/H-X({X≠E})) identified previously in the C terminus of the mitochondrial isoform of the TA protein cytochrome b 5, also exists in many other A. thaliana outer mitochondrial membrane (OMM)-TA proteins. This motif is conspicuously absent, however, in all but one of the TA protein subunits of the translocon at the outer membrane of mitochondria (TOM), suggesting that these two groups of proteins utilize distinct biogenetic pathways. Consistent with this premise, we show that the TA sequences of the dibasic-containing proteins are both necessary and sufficient for targeting to mitochondria, and are interchangeable, while the TA regions of TOM proteins lacking a dibasic motif are necessary, but not sufficient for localization, and cannot be functionally exchanged. We also present results from a comprehensive mutational analysis of the dibasic motif and surrounding sequences that not only greatly expands the functional definition and context-dependent properties of this targeting signal, but also led to the identification of other novel putative OMM-TA proteins. Collectively, these results provide important insight to the complexity of the targeting pathways involved in the biogenesis of OMM-TA proteins and help define a consensus targeting motif that is utilized by at least a subset of these proteins.

A transcript profiling approach reveals an abscisic acid-specific glycosyltransferase (UGT73C14) induced in developing fiber of Ligon lintless-2 mutant of cotton (Gossypium hirsutum L.)

Ligon lintless-2, a monogenic dominant cotton (Gossypium hirsutum L.) fiber mutation, causing extreme reduction in lint fiber length with no pleiotropic effects on vegetative growth, represents an excellent model system to study fiber elongation. A UDP-glycosyltransferase that was highly expressed in developing fibers of the mutant Ligon lintless-2 was isolated. The predicted amino acid sequence showed ~53% similarity with Arabidopsis UGT73C sub-family members and the UDP-glycosyltransferase was designated as UGT73C14. When expressed in Escherichia coli as a recombinant protein with a maltose binding protein tag, UGT73C14 displayed enzymatic activity toward ABA and utilized UDP-glucose and UDP-galactose as the sugar donors. The recombinant UGT73C14 converted natural occurring isoform (+)-cis, trans-ABA better than (+)-trans, trans-ABA and (-)-cis, trans-ABA. Transgenic Arabidopsis plants constitutively overexpressing UGT73C14 did not show phenotypic changes under standard growth conditions. However, the increased glycosylation of ABA resulted in phenotypic changes in post-germinative growth and seedling establishment, confirming in vivo activity of UGT73C14 for ABA. This suggests that the expression level of UGT73C14 is regulated by the observed elevated levels of ABA in developing fibers of the Li₂ mutant line and may be involved in the regulation of ABA homeostasis.

Plant P450s as versatile drivers for evolution of species-specific chemical diversity

The irreversible nature of reactions catalysed by P450s makes these enzymes landmarks in the evolution of plant metabolic pathways. Founding members of P450 families are often associated with general (i.e. primary) metabolic pathways, restricted to single copy or very few representatives, indicative of purifying selection. Recruitment of those and subsequent blooms into multi-member gene families generates genetic raw material for functional diversification, which is an inherent characteristic of specialized (i.e. secondary) metabolism. However, a growing number of highly specialized P450s from not only the CYP71 clan indicate substantial contribution of convergent and divergent evolution to the observed general and specialized metabolite diversity. We will discuss examples of how the genetic and functional diversification of plant P450s drives chemical diversity in light of plant evolution. Even though it is difficult to predict the function or substrate of a P450 based on sequence similarity, grouping with a family or subfamily in phylogenetic trees can indicate association with metabolism of particular classes of compounds. Examples will be given that focus on multi-member gene families of P450s involved in the metabolic routes of four classes of specialized metabolites: cyanogenic glucosides, glucosinolates, mono- to triterpenoids and phenylpropanoids.

A combined biochemical screen and TILLING approach identifies mutations in Sorghum bicolor L. Moench resulting in acyanogenic forage production

Cyanogenic glucosides are present in several crop plants and can pose a significant problem for human and animal consumption, because of their ability to release toxic hydrogen cyanide. Sorghum bicolor L. contains the cyanogenic glucoside dhurrin. A qualitative biochemical screen of the M2 population derived from EMS treatment of sorghum seeds, followed by the reverse genetic technique of Targeted Induced Local Lesions in Genomes (TILLING), was employed to identify mutants with altered hydrogen cyanide potential (HCNp). Characterization of these plants identified mutations affecting the function or expression of dhurrin biosynthesis enzymes, and the ability of plants to catabolise dhurrin. The main focus in this study is on acyanogenic or low cyanide releasing lines that contain mutations in CYP79A1, the cytochrome P450 enzyme catalysing the first committed step in dhurrin synthesis. Molecular modelling supports the measured effects on CYP79A1 activity in the mutant lines. Plants harbouring a P414L mutation in CYP79A1 are acyanogenic when homozygous for this mutation and are phenotypically normal, except for slightly slower growth at early seedling stage. Detailed biochemical analyses demonstrate that the enzyme is present in wild-type amounts but is catalytically inactive. Additional mutants capable of producing dhurrin at normal levels in young seedlings but with negligible leaf dhurrin levels in mature plants were also identified. No mutations were detected in the coding sequence of dhurrin biosynthetic genes in this second group of mutants, which are as tall or taller, and leafier than nonmutated lines. These sorghum mutants with reduced or negligible dhurrin content may be ideally suited for forage production.

Two glycosyltransferases involved in anthocyanin modification delineated by transcriptome independent component analysis in Arabidopsis thaliana

To identify candidate genes involved in Arabidopsis flavonoid biosynthesis, we applied transcriptome coexpression analysis and independent component analyses with 1388 microarray data from publicly available databases. Two glycosyltransferases, UGT79B1 and UGT84A2 were found to cluster with anthocyanin biosynthetic genes. Anthocyanin was drastically reduced in ugt79b1 knockout mutants. Recombinant UGT79B1 protein converted cyanidin 3-O-glucoside to cyanidin 3-O-xylosyl(1→2)glucoside. UGT79B1 recognized 3-O-glucosylated anthocyanidins/flavonols and uridine diphosphate (UDP)-xylose, but not 3,5-O-diglucosylated anthocyanidins, indicating that UGT79B1 encodes anthocyanin 3-O-glucoside: 2''-O-xylosyltransferase. UGT84A2 is known to encode sinapic acid: UDP-glucosyltransferase. In ugt84a2 knockout mutants, a major sinapoylated anthocyanin was drastically reduced. A comparison of anthocyanin profiles in ugt84a knockout mutants indicated that UGT84A2 plays a major role in sinapoylation of anthocyanin, and that other UGT84As contribute the production of 1-O-sinapoylglucose to a lesser extent. These data suggest major routes from cyanidin 3-O-glucoside to the most highly modified cyanidin in the potential intricate anthocyanin modification pathways in Arabidopsis.

Cytochromes p450

There are 244 cytochrome P450 genes (and 28 pseudogenes) in the Arabidopsis genome. P450s thus form one of the largest gene families in plants. Contrary to what was initially thought, this family diversification results in very limited functional redundancy and seems to mirror the complexity of plant metabolism. P450s sometimes share less than 20% identity and catalyze extremely diverse reactions leading to the precursors of structural macromolecules such as lignin, cutin, suberin and sporopollenin, or are involved in biosynthesis or catabolism of all hormone and signaling molecules, of pigments, odorants, flavors, antioxidants, allelochemicals and defense compounds, and in the metabolism of xenobiotics. The mechanisms of gene duplication and diversification are getting better understood and together with co-expression data provide leads to functional characterization.

Natural diversity and adaptation in plant secondary metabolism

Technological advances in metabolomics, transcriptomics and genomics have facilitated the detection of genes that contribute to diversification in plant secondary metabolism. Statistical tools from molecular population genetics may help in evaluating whether the corresponding genes or genomic regions carry a signature of selection and answering the question of whether novel compounds are 'adaptive'. Gene duplication fuels diversification in plant secondary metabolism and the evolutionary mechanism for adaptation may follow a path of neofunctionalization subsequent to gene duplication, or gene duplication may occur subsequent to--and resolve--an adaptive conflict present in a single ancestral gene sequence.

Molecular activities, biosynthesis and evolution of triterpenoid saponins

Saponins are bioactive compounds generally considered to be produced by plants to counteract pathogens and herbivores. Besides their role in plant defense, saponins are of growing interest for drug research as they are active constituents of several folk medicines and provide valuable pharmacological properties. Accordingly, much effort has been put into unraveling the modes of action of saponins, as well as in exploration of their potential for industrial processes and pharmacology. However, the exploitation of saponins for bioengineering crop plants with improved resistances against pests as well as circumvention of laborious and uneconomical extraction procedures for industrial production from plants is hampered by the lack of knowledge and availability of genes in saponin biosynthesis. Although the ability to produce saponins is rather widespread among plants, a complete synthetic pathway has not been elucidated in any single species. Current conceptions consider saponins to be derived from intermediates of the phytosterol pathway, and predominantly enzymes belonging to the multigene families of oxidosqualene cyclases (OSCs), cytochromes P450 (P450s) and family 1 UDP-glycosyltransferases (UGTs) are thought to be involved in their biosynthesis. Formation of unique structural features involves additional biosynthetical enzymes of diverse phylogenetic background. As an example of this, a serine carboxypeptidase-like acyltransferase (SCPL) was recently found to be involved in synthesis of triterpenoid saponins in oats. However, the total number of identified genes in saponin biosynthesis remains low as the complexity and diversity of these multigene families impede gene discovery based on sequence analysis and phylogeny. This review summarizes current knowledge of triterpenoid saponin biosynthesis in plants, molecular activities, evolutionary aspects and perspectives for further gene discovery.

An evolutionary view of functional diversity in family 1 glycosyltransferases

Glycosyltransferases (GTs) (EC 2.4.x.y) catalyze the transfer of sugar moieties to a wide range of acceptor molecules, such as sugars, lipids, proteins, nucleic acids, antibiotics and other small molecules, including plant secondary metabolites. These enzymes can be classified into at least 92 families, of which family 1 glycosyltransferases (GT1), often referred to as UDP glycosyltransferases (UGTs), is the largest in the plant kingdom. To understand how UGTs expanded in both number and function during evolution of land plants, we screened genome sequences from six plants (Physcomitrella patens, Selaginella moellendorffii, Populus trichocarpa, Oryza sativa, Arabidopsis thaliana and Arabidopsis lyrata) for the presence of a conserved UGT protein domain. Phylogenetic analyses of the UGT genes revealed a significant expansion of UGTs, with lineage specificity and a higher duplication rate in vascular plants after the divergence of Physcomitrella. The UGTs from the six species fell into 24 orthologous groups that contained genes derived from the common ancestor of these six species. Some orthologous groups contained multiple UGT families with known functions, suggesting that UGTs discriminate compounds as substrates in a lineage-specific manner. Orthologous groups containing only a single UGT family tend to play a crucial role in plants, suggesting that such UGT families may have not expanded because of evolutionary constraints.

Tail-anchored membrane proteins: exploring the complex diversity of tail-anchored-protein targeting in plant cells

Tail-anchored (TA) proteins are special class of integral membrane proteins that in recent years have received a considerable amount of attention due to their diverse cellular functions and unique targeting and insertion mechanisms. Defined by the presence of a single, hydrophobic membrane-spanning domain at or near their C terminus, TA proteins must be inserted into membranes post-translationally and are orientated such that their larger N-terminal domain (most often the functional domain) faces the cytosol, while their shorter C-terminal domain faces the interior of the organelle. The C-terminal domain of TA proteins also usually contains the information responsible for their selective targeting to the proper subcellular membrane, a process that, based primarily on studies with yeasts and mammals, appears to be highly complex due to the presence of multiple pathways. Within this context, we discuss here the biogenesis of plant TA proteins and the potential for hundreds of new TA proteins identified via bioinformatics screens to contribute to the already remarkable number of roles that this class of membrane proteins participates in throughout plant growth and development.

Whole genome co-expression analysis of soybean cytochrome P450 genes identifies nodulation-specific P450 monooxygenases

BACKGROUND

Cytochrome P450 monooxygenases (P450s) catalyze oxidation of various substrates using oxygen and NAD(P)H. Plant P450s are involved in the biosynthesis of primary and secondary metabolites performing diverse biological functions. The recent availability of the soybean genome sequence allows us to identify and analyze soybean putative P450s at a genome scale. Co-expression analysis using an available soybean microarray and Illumina sequencing data provides clues for functional annotation of these enzymes. This approach is based on the assumption that genes that have similar expression patterns across a set of conditions may have a functional relationship.

RESULTS

We have identified a total number of 332 full-length P450 genes and 378 pseudogenes from the soybean genome. From the full-length sequences, 195 genes belong to A-type, which could be further divided into 20 families. The remaining 137 genes belong to non-A type P450s and are classified into 28 families. A total of 178 probe sets were found to correspond to P450 genes on the Affymetrix soybean array. Out of these probe sets, 108 represented single genes. Using the 28 publicly available microarray libraries that contain organ-specific information, some tissue-specific P450s were identified. Similarly, stress responsive soybean P450s were retrieved from 99 microarray soybean libraries. We also utilized Illumina transcriptome sequencing technology to analyze the expressions of all 332 soybean P450 genes. This dataset contains total RNAs isolated from nodules, roots, root tips, leaves, flowers, green pods, apical meristem, mock-inoculated and Bradyrhizobium japonicum-infected root hair cells. The tissue-specific expression patterns of these P450 genes were analyzed and the expression of a representative set of genes were confirmed by qRT-PCR. We performed the co-expression analysis on many of the 108 P450 genes on the Affymetrix arrays. First we confirmed that CYP93C5 (an isoflavone synthase gene) is co-expressed with several genes encoding isoflavonoid-related metabolic enzymes. We then focused on nodulation-induced P450s and found that CYP728H1 was co-expressed with the genes involved in phenylpropanoid metabolism. Similarly, CYP736A34 was highly co-expressed with lipoxygenase, lectin and CYP83D1, all of which are involved in root and nodule development.

CONCLUSIONS

The genome scale analysis of P450s in soybean reveals many unique features of these important enzymes in this crop although the functions of most of them are largely unknown. Gene co-expression analysis proves to be a useful tool to infer the function of uncharacterized genes. Our work presented here could provide important leads toward functional genomics studies of soybean P450s and their regulatory network through the integration of reverse genetics, biochemistry, and metabolic profiling tools. The identification of nodule-specific P450s and their further exploitation may help us to better understand the intriguing process of soybean and rhizobium interaction.

Flavocytochrome P450 BM3 mutant W1046A is a NADH-dependent fatty acid hydroxylase: implications for the mechanism of electron transfer in the P450 BM3 dimer

Bacillus megaterium P450 BM3 (BM3) is a P450/P450 reductase fusion enzyme, where the dimer is considered the active form in NADPH-dependent fatty acid hydroxylation. The BM3 W1046A mutant was generated, removing an aromatic "shield" from its FAD isoalloxazine ring. W1046A BM3 is a catalytically active NADH-dependent lauric acid hydroxylase, with product formation slightly superior to the NADPH-driven enzyme. The W1046A BM3 K(m) for NADH is 20-fold lower than wild-type BM3, and catalytic efficiency of W1046A BM3 with NADH and NADPH are similar in lauric acid oxidation. Wild-type BM3 also catalyzes NADH-dependent lauric acid hydroxylation, but less efficiently than W1046A BM3. A hypothesis that W1046A BM3 is inactive [15] helped underpin a model of electron transfer from FAD in one BM3 monomer to FMN in the other in order to drive fatty acid hydroxylation in native BM3. Our data showing W1046A BM3 is a functional fatty acid hydroxylase are consistent instead with a BM3 catalytic model involving electron transfer within a reductase monomer, and from FMN of one monomer to heme of the other [12]. W1046A BM3 is an efficient NADH-utilizing fatty acid hydroxylase with potential biotechnological applications.

Conformational changes of the NADPH-dependent cytochrome P450 reductase in the course of electron transfer to cytochromes P450

The NADPH-dependent cytochrome P450 reductase (CPR) is a key electron donor to eucaryotic cytochromes P450 (CYPs). CPR shuttles electrons from NADPH through the FAD and FMN-coenzymes into the iron of the prosthetic heme-group of the CYP. In the course of these electron transfer reactions, CPR undergoes large conformational changes. This mini-review discusses the new evidence provided for such conformational changes involving a combination of a "swinging" and "rotating" model and highlights the molecular mechanisms by which formation of these conformations are controlled and thereby enables CPR to serve as an effective electron transferring "nano-machine".

ATR3 encodes a diflavin reductase essential for Arabidopsis embryo development

*The Arabidopsis genome possesses two confirmed Cytochrome P450 Reductase (CPR) genes, ATR1 and ATR2, together with a third putative homologue, ATR3, which annotation is questionable. *Phylogenetic analysis classified ATR3 as a CPR-like protein sharing homologies with the animal cytosolic dual flavin reductases, NR1 and Fre-1, distinct from the microsomal CPRs, ATR1 and ATR2. Like NR1 and Fre-1, ATR3 lacks the N-terminal endoplasmic reticulum (ER) anchor domain of CPRs and is localized in the cytoplasm. Recombinant ATR3 in plant soluble extracts was able to reduce cytochrome c but failed to reduce the human P450 CYP1A2. *Loss of ATR3 function resulted in early embryo lethality indicating that this reductase activity is essential. A yeast 2-hybrid screen identified a unique interaction of ATR3 with the homologue of the human anti-apoptotic CIAPIN1 and the yeast Dre2 protein. *This interaction suggests two possible roles for ATR3 in the control of cell death and in chromosome segregation at mitosis. Consistent with these results, the promoter of ATR3 is activated during cell cycle progression. Together these results demonstrated that ATR3 belongs to the NR1 subfamily of diflavin reductases whose characterized members are involved in essential cellular functions.

Plant NADPH-cytochrome P450 oxidoreductases

NADPH-cytochrome P450 oxidoreductase (CPR) serves as the electron donor to almost all eukaryotic cytochromes P450. It belongs to a small family of diflavin proteins and is built of cofactor binding domains with high structural homology to those of bacterial flavodoxins and to ferredoxin-NADP(+) oxidoreductases. CPR shuttles electrons from NADPH through the FAD and FMN-cofactors into the central heme-group of the P450s. Mobile domains in CPR are essential for electron transfer between FAD and FMN and for P450 interaction. Blast searches identified 54 full-length gene sequences encoding CPR derived from a total of 35 different plant species. CPRs from vascular plants cluster into two major phylogenetic groups. Depending on the species, plants contain one, two or three paralogs of which one is inducible. The nature of the CPR-P450 interacting domains is well conserved as demonstrated by the ability of CPRs from different species or even from different kingdoms to at least partially complement each other functionally. This makes CPR an ideal bio-brick in synthetic biology approaches to re-design or develop entirely different combinations of existing biological systems to gain improved or completely altered functionalities based on the "share your parts" principle.

454 pyrosequencing based transcriptome analysis of Zygaena filipendulae with focus on genes involved in biosynthesis of cyanogenic glucosides

Background

An essential driving component in the co-evolution of plants and insects is the ability to produce and handle bioactive compounds. Plants produce bioactive natural products for defense, but some insects detoxify and/or sequester the compounds, opening up for new niches with fewer competitors. To study the molecular mechanism behind the co-adaption in plant-insect interactions, we have investigated the interactions between Lotus corniculatus and Zygaena filipendulae. They both contain cyanogenic glucosides which liberate toxic hydrogen cyanide upon breakdown. Moths belonging to the Zygaena family are the only insects known, able to carry out both de novo biosynthesis and sequestration of the same cyanogenic glucosides as those from their feed plants. The biosynthetic pathway for cyanogenic glucoside biosynthesis in Z. filipendulae proceeds using the same intermediates as in the well known pathway from plants, but none of the enzymes responsible have been identified. A genomics strategy founded on 454 pyrosequencing of the Z. filipendulae transcriptome was undertaken to identify some of these enzymes in Z. filipendulae.

Results

Comparisons of the Z. filipendulae transcriptome with the sequenced genomes of Bombyx mori, Drosophila melanogaster, Tribolium castaneum, Apis mellifera and Anopheles gambiae indicate a high coverage of the Z. filipendulae transcriptome. 11% of the Z. filipendulae transcriptome sequences were assigned to Gene Ontology categories. Candidate genes for enzymes functioning in the biosynthesis of cyanogenic glucosides (cytochrome P450 and family 1 glycosyltransferases) were identified based on sequence length, number of copies and presence/absence of close homologs in D. melanogaster, B. mori and the cyanogenic butterfly Heliconius. Examination of biased codon usage, GC content and selection on gene candidates support the notion of cyanogenesis as an "old" trait within Ditrysia, as well as its origins being convergent between plants and insects.

Conclusion

Pyrosequencing is an attractive approach to gain access to genes in the biosynthesis of bio-active natural products from insects and other organisms, for which the genome sequence is not known. Based on analysis of the Z. filipendulae transcriptome, promising gene candidates for biosynthesis of cyanogenic glucosides was identified, and the suitability of Z. filipendulae as a model system for cyanogenesis in insects is evident.