Coevolutionary Information Captures Catalytic Functions and Reveals Divergent Roles of Terpene Synthase Interdomain Connections.

Inferring the historical and biophysical causes of diversity within protein families is a complex puzzle. A key to unraveling this problem is characterizing the rugged topography of sequence-function adaptive landscapes. Using biochemical data from a 29 = 512 combinatorial library of tobacco 5-epi-aristolochene synthase (TEAS) mutants engineered to make the native major product of Egyptian henbane premnaspirodiene synthase (HPS) and a complementary 512 mutant HPS library, we address the question of how product specificity is controlled. These data sets reveal that HPS is far more robust and resistant to mutations than TEAS, where most mutants are promiscuous. We also combine experimental data with a sequence Potts Hamiltonian model and direct coupling analysis to quantify mutant fitness. Our results demonstrate that the Hamiltonian captures variation in product outputs across both libraries, clusters native family members based on their substrate specificities, and exposes the divergent catalytic roles of couplings between the catalytic and noncatalytic domains of TEAS versus HPS. Specifically, we found that the role of the interdomain connectivities in specifying product output is more important in TEAS than connectivities within the catalytic domain. Despite being 75% identical, this property is not shared by HPS, where connectivities within the catalytic domain are more important for specificity. By solving the X-ray crystal structure of HPS, we assessed structural bases for their interdomain network differences. Last, we calculate the product profile Shannon entropies of the two libraries, which showcases that site-site connectivities also play divergent roles in catalytic accuracy.

High-throughput discovery of plastid genes causing albino phenotypes in ornamental chimeric plants.

Chimeric plants composed of green and albino tissues have great ornamental value. To unveil the functional genes responsible for albino phenotypes in chimeric plants, we inspected the complete plastid genomes (plastomes) in green and albino leaf tissues from 23 ornamental chimeric plants belonging to 20 species, including monocots, dicots, and gymnosperms. In nine chimeric plants, plastomes were identical between green and albino tissues. Meanwhile, another 14 chimeric plants were heteroplasmic, showing a mutation between green and albino tissues. We identified 14 different point mutations in eight functional plastid genes related to plastid-encoded RNA polymerase (rpo) or photosystems which caused albinism in the chimeric plants. Among them, 12 were deleterious mutations in the target genes, in which early termination appeared due to small deletion-mediated frameshift or single nucleotide substitution. Another was single nucleotide substitution in an intron of the ycf3 and the other was a missense mutation in coding region of the rpoC2 gene. We inspected chlorophyll structure, protein functional model of the rpoC2, and expression levels of the related genes in green and albino tissues of Reynoutria japonica. A single amino acid change, histidine-to-proline substitution, in the rpoC2 protein may destabilize the peripheral helix of plastid-encoded RNA polymerase, impairing the biosynthesis of the photosynthesis system in the albino tissue of R. japonica chimera plant.

Comparative transcriptome and metabolome analyses of four Panax species explore the dynamics of metabolite biosynthesis.

Background: The genus Panax in the Araliaceae family has been used as traditional medicinal plants worldwide and is known to biosynthesize ginsenosides and phytosterols. However, genetic variation between Panax species has influenced their biosynthetic pathways is not fully understood. Methods: Simultaneous analysis of transcriptomes and metabolomes obtained from adventitious roots of two tetraploid species (Panax ginseng and P. quinquefolius) and two diploid species (P. notoginseng and P. vietnamensis) revealed the diversity of their metabolites and related gene expression profiles. Results: The transcriptome analysis showed that 2,3-OXIDOSQUALENE CYCLASEs (OSCs) involved in phytosterol biosynthesis are upregulated in the diploid species, while the expression of OSCs contributing to ginsenoside biosynthesis is higher in the tetraploid species. In agreement with these results, the contents of dammarenediol-type ginsenosides were higher in the tetraploid species relative to the diploid species. Conclusion: These results suggest that a whole-genome duplication event has influenced the triterpene biosynthesis pathway in tetraploid Panax species during their evolution or ecological adaptation. This study provides a basis for further efforts to explore the genetic variation of the Panax genus.

Complete chloroplast genome of Gypsophila oldhamiana Miq. (Caryophyllales: Caryophyllaceae).

The complete chloroplast genome sequence of Gypsophila oldhamiana Miq., a species of the Caryophyllaceae family, was assembled and analyzed from the de novo assembly of Illumina paired-end sequencing data. The total length of the chloroplast genome of G. oldhamiana was 152,675 bp including a large single-copy (LSC) region of 83,552 bp, a small single-copy (SSC) region of 17,349 bp, and a pair of identical inverted repeat regions (IRs) of 25,887 bp. The genome possessed a total of 130 genes, including 85 protein-coding genes, 37 transfer RNA (tRNA) genes, and 8 ribosomal RNA (rRNA) genes. The phylogenetic analysis of G. oldhamiana with 14 related species discovered the closest taxonomical relationship with Gypsophila vaccaria voucher in the Caryophyllaceae family.

RNA editing may stabilize membrane-embedded proteins by increasing phydrophobicity: A study of Zanthoxylum piperitum and Z. schinifolium chloroplast NdhG.

Most chloroplast genes in Zanthoxylum schinifolium (Korean pepper) and Z. piperitum (Japanese pepper) are subject to neutral or purifying selection (Ka/Ks values < 1); however, NAD(P)H dehydrogenase subunit G (ndhG) has a Ka/Ks value of 1.43, which may indicate positive selection. Here, we modeled the ZsNdhG and ZpNdhG structures by comparing them with the NuoJ subunit of respiratory complex I in Escherichia coli, revealing the locations of four amino acid differences between ZsNdhG and ZpNdhG. As these polymorphisms were located at the end of a membrane-spanning α-helix or in peptide loops external to the membrane, they are not expected to have major effects on the membrane-embedding properties of these proteins. However, we found that C-to-U RNA editing occurred at the ndhG-50 sites of ndhG (uCa to uUa, Ser to Leu) in both species, resulting in changes to an amino acid located in the middle of a membrane-spanning α-helix, which may maintain its hydrophobicity. RNA editing at the ndhG-50 site was conserved in many plant species, and the modeled structures of Anthoceros formosae NdhG and Spirodela polyrhiza NdhB provided further evidence that RNA editing increases the hydrophobicity of membrane-embedded proteins. We also speculate that the polar residues inside membrane-spanning α-helices serve to support the protein structure. This report represents the first RNA-editing site identified in Zanthoxylum and points to the importance of considering RNA editing when identifying positively selected genes based on Ka/Ks values.

Identification of candidate UDP-glycosyltransferases involved in protopanaxadiol-type ginsenoside biosynthesis in Panax ginseng.

Ginsenosides are dammarane-type or triterpenoidal saponins that contribute to the various pharmacological activities of the medicinal herb Panax ginseng. The putative biosynthetic pathway for ginsenoside biosynthesis is known in P. ginseng, as are some of the transcripts and enzyme-encoding genes. However, few genes related to the UDP-glycosyltransferases (UGTs), enzymes that mediate glycosylation processes in final saponin biosynthesis, have been identified. Here, we generated three replicated Illumina RNA-Seq datasets from the adventitious roots of P. ginseng cultivar Cheongsun (CS) after 0, 12, 24, and 48 h of treatment with methyl jasmonate (MeJA). Using the same CS cultivar, metabolomic data were also generated at 0 h and every 12-24 h thereafter until 120 h of MeJA treatment. Differential gene expression, phylogenetic analysis, and metabolic profiling were used to identify candidate UGTs. Eleven candidate UGTs likely to be involved in ginsenoside glycosylation were identified. Eight of these were considered novel UGTs, newly identified in this study, and three were matched to previously characterized UGTs in P. ginseng. Phylogenetic analysis further asserted their association with ginsenoside biosynthesis. Additionally, metabolomic analysis revealed that the newly identified UGTs might be involved in the elongation of glycosyl chains of ginsenosides, especially of protopanaxadiol (PPD)-type ginsenosides.

Genome and evolution of the shade-requiring medicinal herb Panax ginseng.

Panax ginseng C. A. Meyer, reputed as the king of medicinal herbs, has slow growth, long generation time, low seed production and complicated genome structure that hamper its study. Here, we unveil the genomic architecture of tetraploid P. ginseng by de novo genome assembly, representing 2.98 Gbp with 59 352 annotated genes. Resequencing data indicated that diploid Panax species diverged in association with global warming in Southern Asia, and two North American species evolved via two intercontinental migrations. Two whole genome duplications (WGD) occurred in the family Araliaceae (including Panax) after divergence with the Apiaceae, the more recent one contributing to the ability of P. ginseng to overwinter, enabling it to spread broadly through the Northern Hemisphere. Functional and evolutionary analyses suggest that production of pharmacologically important dammarane-type ginsenosides originated in Panax and are produced largely in shoot tissues and transported to roots; that newly evolved P. ginseng fatty acid desaturases increase freezing tolerance; and that unprecedented retention of chlorophyll a/b binding protein genes enables efficient photosynthesis under low light. A genome-scale metabolic network provides a holistic view of Panax ginsenoside biosynthesis. This study provides valuable resources for improving medicinal values of ginseng either through genomics-assisted breeding or metabolic engineering.

The complete chloroplast genome sequence of Korean Lonicera japonica and intra-species diversity.

Lonicera japonica is a traditional medicinal plant well known for its anti-inflammatory effect. The complete chloroplast genome sequence of L. japonica collected from Korea was obtained by de novo assembly using whole genome sequence data. The chloroplast genome is 155,060 bp in length, containing 88,853 bp in a large single copy (LSC), 18,653 bp in a small single copy (SSC) and 23,777 bp in a pair of inverted repeats (IRs). A total of 112 genes including 78 protein-coding genes and 34 structural RNA genes were identified. The sequence comparison of two L. japonica collected from Korea and China revealed 48 single nucleotide polymorphisms (SNPs) and 45 insertions/deletions (InDels). In addition, phylogenetic analysis represented intraspecific diversity within L. japonica species collected in Korea and China.

Authentication of Zanthoxylum Species Based on Integrated Analysis of Complete Chloroplast Genome Sequences and Metabolite Profiles.

We performed chloroplast genome sequencing and comparative analysis of two Rutaceae species, Zanthoxylum schinifolium (Korean pepper tree) and Z. piperitum (Japanese pepper tree), which are medicinal and culinary crops in Asia. We identified more than 837 single nucleotide polymorphisms and 103 insertions/deletions (InDels) based on a comparison of the two chloroplast genomes and developed seven DNA markers derived from five tandem repeats and two InDel variations that discriminated between Korean Zanthoxylum species. Metabolite profile analysis pointed to three metabolic groups, one with Korean Z. piperitum samples, one with Korean Z. schinifolium samples, and the last containing all the tested Chinese Zanthoxylum species samples, which are considered to be Z. bungeanum based on our results. Two markers were capable of distinguishing among these three groups. The chloroplast genome sequences identified in this study represent a valuable genomics resource for exploring diversity in Rutaceae, and the molecular markers will be useful for authenticating dried Zanthoxylum berries in the marketplace.

Integrated Transcriptomic and Metabolomic Analysis of Five Panax ginseng Cultivars Reveals the Dynamics of Ginsenoside Biosynthesis.

Panax ginseng C.A. Meyer is a traditional medicinal herb that produces bioactive compounds such as ginsenosides. Here, we investigated the diversity of ginsenosides and related genes among five genetically fixed inbred ginseng cultivars (Chunpoong [CP], Cheongsun [CS], Gopoong [GO], Sunhyang [SH], and Sunun [SU]). To focus on the genetic diversity related to ginsenoside biosynthesis, we utilized in vitro cultured adventitious roots from the five cultivars grown under controlled environmental conditions. PCA loading plots based on secondary metabolite composition classified the five cultivars into three groups. We selected three cultivars (CS, SH, and SU) to represent the three groups and conducted further transcriptome and gas chromatography-mass spectrometry analyses to identify genes and intermediates corresponding to the variation in ginsenosides among cultivars. We quantified ginsenoside contents from the three cultivars. SH had more than 12 times the total ginsenoside content of CS, with especially large differences in the levels of panaxadiol-type ginsenosides. The expression levels of genes encoding squalene epoxidase (SQE) and dammarenediol synthase (DDS) were also significantly lower in CS than SH and SU, which is consistent with the low levels of ginsenoside produced in this cultivar. Methyl jasmonate (MeJA) treatment increased the levels of panaxadiol-type ginsenosides up to 4-, 13-, and 31-fold in SH, SU, and CS, respectively. MeJA treatment also greatly increased the quantity of major intermediates and the expression of the underlying genes in the ginsenoside biosynthesis pathway; these intermediates included squalene, 2,3-oxidosqualene, and dammarenediol II, especially in CS, which had the lowest ginsenoside content under normal culture conditions. We conclude that SQE and DDS are the most important genetic factors for ginsenoside biosynthesis with diversity among ginseng cultivars.

Biosynthetic potential of sesquiterpene synthases: product profiles of Egyptian Henbane premnaspirodiene synthase and related mutants.

The plant terpene synthase (TPS) family is responsible for the biosynthesis of a variety of terpenoid natural products possessing diverse biological functions. TPSs catalyze the ionization and, most commonly, rearrangement and cyclization of prenyl diphosphate substrates, forming linear and cyclic hydrocarbons. Moreover, a single TPS often produces several minor products in addition to a dominant product. We characterized the catalytic profiles of Hyoscyamus muticus premnaspirodiene synthase (HPS) and compared it with the profile of a closely related TPS, Nicotiana tabacum 5-epi-aristolochene synthase (TEAS). The profiles of two previously studied HPS and TEAS mutants, each containing nine interconverting mutations, dubbed HPS-M9 and TEAS-M9, were also characterized. All four TPSs were compared under varying temperature and pH conditions. In addition, we solved the X-ray crystal structures of TEAS and a TEAS quadruple mutant complexed with substrate and products to gain insight into the enzymatic features modulating product formation. These informative structures, along with product profiles, provide new insight into plant TPS catalytic promiscuity.

Formation of a Novel Macrocyclic Alkaloid from the Unnatural Farnesyl Diphosphate Analogue Anilinogeranyl Diphosphate by 5-Epi-Aristolochene Synthase.

As part of an effort to identify substrate analogs suitable for helping to resolve structural features important for terpene synthases, the inhibition of 5-epi-aristolochene biosynthesis from farnesyl diphosphate (FPP) by the tobacco 5-epi-aristolochene synthase incubated with anilinogeranyl diphosphate (AGPP) was examined. The apparent noncompetitive nature of the inhibition supported further assessment of how AGPP might be bound to crystallographic forms of the enzyme. Surprisingly, the bound form of the inhibitor appeared to have undergone a cyclization event consistent with the native mechanism associated with FPP catalysis. Biocatalytic formation of a novel 13-membered macrocyclic paracyclophane alkaloid was confirmed by high-resolution GC-MS and NMR analysis. This work provides insights into new biosynthetic means for generating novel, functionally diversified, medium-sized terpene alkaloids.

Small heat shock proteins can release light dependence of tobacco seed during germination.

Small heat shock proteins (sHSPs) function as ATP-independent molecular chaperones, and although the production and function of sHSPs have often been described under heat stress, the expression and function of sHSPs in fundamental developmental processes, such as pollen and seed development, have also been confirmed. Seed germination involves the breaking of dormancy and the resumption of embryo growth that accompany global changes in transcription, translation, and metabolism. In many plants, germination is triggered simply by imbibition of water; however, different seeds require different conditions in addition to water. For small-seeded plants, like Arabidopsis (Arabidopsis thaliana), lettuce (Lactuca sativa), tomato (Solanum lycopersicum), and tobacco (Nicotiana tabacum), light is an important regulator of seed germination. The facts that sHSPs accumulate during seed development, sHSPs interact with various client proteins, and seed germination accompanies synthesis and/or activation of diverse proteins led us to investigate the role of sHSPs in seed germination, especially in the context of light dependence. In this study, we have built transgenic tobacco plants that ectopically express sHSP, and the effect was germination of the seeds in the dark. Administering heat shock to the seeds also resulted in the alleviation of light dependence during seed germination. Subcellular localization of ectopically expressed sHSP was mainly observed in the cytoplasm, whereas heat shock-induced sHSPs were transported to the nucleus. We hypothesize that ectopically expressed sHSPs in the cytoplasm led the status of cytoplasmic proteins involved in seed germination to function during germination without additional stimulus and that heat shock can be another signal that induces seed germination.

Tobacco class I cytosolic small heat shock proteins are under transcriptional and translational regulations in expression and heterocomplex prevails under the high-temperature stress condition in vitro.

Seven genomic clones of tobacco (Nicotiana tabacum W38) cytosolic class I small heat shock proteins (sHSPs), probably representing all members in the class, were isolated and found to have 66 to 92% homology between their nucleotide sequences. Even though all seven sHSP genes showed heat shock-responsive accumulation of their transcripts and proteins, each member showed discrepancies in abundance and timing of expression upon high-temperature stress. This was mainly the result of transcriptional regulation during mild stress conditions and transcriptional and translational regulation during strong stress conditions. Open reading frames (ORFs) of these genomic clones were expressed in Escherichia coli and the sHSPs were purified from E. coli. The purified tobacco sHSPs rendered citrate synthase and luciferase soluble under high temperatures. At room temperature, non-denaturing pore exclusion polyacrylamide gel electrophoresis on three sHSPs demonstrated that the sHSPs spontaneously formed homo-oligomeric complexes of 200 ∼ 240 kDa. However, under elevated temperatures, hetero-oligomeric complexes between the sHSPs gradually prevailed. Atomic force microscopy showed that the hetero-oligomer of NtHSP18.2/NtHSP18.3 formed a stable oligomeric particle similar to that of the NtHSP18.2 homo-oligomer. These hetero-oligomers positively influenced the revival of thermally inactivated luciferase. Amino acid residues mainly in the N-terminus are suggested for the exchange of the component sHSPs and the formation of dominant hetero-oligomers under high temperatures.

Ginger and turmeric expressed sequence tags identify signature genes for rhizome identity and development and the biosynthesis of curcuminoids, gingerols and terpenoids.

BACKGROUND: Ginger (Zingiber officinale) and turmeric (Curcuma longa) accumulate important pharmacologically active metabolites at high levels in their rhizomes. Despite their importance, relatively little is known regarding gene expression in the rhizomes of ginger and turmeric. RESULTS: In order to identify rhizome-enriched genes and genes encoding specialized metabolism enzymes and pathway regulators, we evaluated an assembled collection of expressed sequence tags (ESTs) from eight different ginger and turmeric tissues. Comparisons to publicly available sorghum rhizome ESTs revealed a total of 777 gene transcripts expressed in ginger/turmeric and sorghum rhizomes but apparently absent from other tissues. The list of rhizome-specific transcripts was enriched for genes associated with regulation of tissue growth, development, and transcription. In particular, transcripts for ethylene response factors and AUX/IAA proteins appeared to accumulate in patterns mirroring results from previous studies regarding rhizome growth responses to exogenous applications of auxin and ethylene. Thus, these genes may play important roles in defining rhizome growth and development. Additional associations were made for ginger and turmeric rhizome-enriched MADS box transcription factors, their putative rhizome-enriched homologs in sorghum, and rhizomatous QTLs in rice. Additionally, analysis of both primary and specialized metabolism genes indicates that ginger and turmeric rhizomes are primarily devoted to the utilization of leaf supplied sucrose for the production and/or storage of specialized metabolites associated with the phenylpropanoid pathway and putative type III polyketide synthase gene products. This finding reinforces earlier hypotheses predicting roles of this enzyme class in the production of curcuminoids and gingerols. CONCLUSION: A significant set of genes were found to be exclusively or preferentially expressed in the rhizome of ginger and turmeric. Specific transcription factors and other regulatory genes were found that were common to the two species and that are excellent candidates for involvement in rhizome growth, differentiation and development. Large classes of enzymes involved in specialized metabolism were also found to have apparent tissue-specific expression, suggesting that gene expression itself may play an important role in regulating metabolite production in these plants.

Suites of terpene synthases explain differential terpenoid production in ginger and turmeric tissues.

The essential oils of ginger (Zingiber officinale) and turmeric (Curcuma longa) contain a large variety of terpenoids, some of which possess anticancer, antiulcer, and antioxidant properties. Despite their importance, only four terpene synthases have been identified from the Zingiberaceae family: (+)-germacrene D synthase and (S)-ß-bisabolene synthase from ginger rhizome, and α-humulene synthase and ß-eudesmol synthase from shampoo ginger (Zingiber zerumbet) rhizome. We report the identification of 25 mono- and 18 sesquiterpene synthases from ginger and turmeric, with 13 and 11, respectively, being functionally characterized. Novel terpene synthases, (-)-caryolan-1-ol synthase and α-zingiberene/ß-sesquiphellandrene synthase, which is responsible for formation of the major sesquiterpenoids in ginger and turmeric rhizomes, were also discovered. These suites of enzymes are responsible for formation of the majority of the terpenoids present in these two plants. Structures of several were modeled, and a comparison of sets of paralogs suggests how the terpene synthases in ginger and turmeric evolved. The most abundant and most important sesquiterpenoids in turmeric rhizomes, (+)-α-turmerone and (+)-ß-turmerone, are produced from (-)-α-zingiberene and (-)-ß-sesquiphellandrene, respectively, via α-zingiberene/ß-sesquiphellandrene oxidase and a still unidentified dehydrogenase.

Tobacco mitochondrial small heat shock protein NtHSP24.6 adopts a dimeric configuration and has a broad range of substrates.

There is a broad range of different small heat shock proteins (sHSPs) that have diverse structural and functional characteristics. To better understand the functional role of mitochondrial sHSP, NtHSP24.6 was expressed in Escherichia coli with a hexahistidine tag and purified. The protein was analyzed by non-denaturing PAGE, chemical cross-linking and size exclusion chromatography and the H6NtHSP24.6 protein was found to form a dimer in solution. The in vitro functional analysis of H6NtHSP24.6 using firefly luciferase and citrate synthase demonstrated that this protein displays typical molecular chaperone activity. When cell lysates of E. coli were heated after the addition of H6NtHSP24.6, a broad range of proteins from 10 to 160 kD in size remained in the soluble state. These results suggest that NtHSP24.6 forms a dimer and can function as a molecular chaperone to protect a diverse range of proteins from thermal aggregation.

Metabolic profiling and phylogenetic analysis of medicinal Zingiber species: Tools for authentication of ginger (Zingiber officinale Rosc).

Phylogenetic analysis and metabolic profiling were used to investigate the diversity of plant material within the ginger species and between ginger and closely related species in the genus Zingiber (Zingiberaceae). In addition, anti-inflammatory data were obtained for the investigated species. Phylogenetic analysis demonstrated that all Zingiber officinale samples from different geographical origins were genetically indistinguishable. In contrast, other Zingiber species were significantly divergent, allowing all species to be clearly distinguished using this analysis. In the metabolic profiling analysis, the Z. officinale samples derived from different origins showed no qualitative differences in major volatile compounds, although they did show some significant quantitative differences in non-volatile composition, particularly regarding the content of [6]-, [8]-, and [10]-gingerols, the most active anti-inflammatory components in this species. The differences in gingerol content were verified by HPLC. The metabolic profiles of other Zingiber species were very different, both qualitatively and quantitatively, when compared to Z. officinale and to each other. Comparative DNA sequence/chemotaxonomic phylogenetic trees showed that the chemical characters of the investigated species were able to generate essentially the same phylogenetic relationships as the DNA sequences. This supports the contention that chemical characters can be used effectively to identify relationships between plant species. Anti-inflammatory in vitro assays to evaluate the ability of all extracts from the Zingiber species examined to inhibit LPS-induced PGE(2) and TNF-alpha production suggested that bioactivity may not be easily predicted by either phylogenetic analysis or gross metabolic profiling. Therefore, identification and quantification of the actual bioactive compounds are required to guarantee the bioactivity of a particular Zingiber sample even after performing authentication by molecular and/or chemical markers.