Enhancing transcriptome analysis in medicinal plants: multiple unigene sets in Astragalus membranaceus.

Astragalus membranaceus is a medicinal plant mainly used in East Asia and contains abundant secondary metabolites. Despite the importance of this plant, the available genomic and genetic information is still limited. De novo transcriptome construction is recognized as an essential method for transcriptome research when reference genome information is incomplete. In this study, we constructed three individual transcriptome sets (unigene sets) for detailed analysis of the phenylpropanoid biosynthesis pathway, a major metabolite of A. membranaceus. Set-1 was a circular consensus sequence (CCS) generated using PacBio sequencing (PacBio-seq). Set-2 consisted of hybridized assembled unigenes with Illumina sequencing (Illumina-seq) reads and PacBio CCS using rnaSPAdes. Set-3 unigenes were assembled from Illumina-seq reads using the Trinity software. Construction of multiple unigene sets provides several advantages for transcriptome analysis. First, it provides an appropriate expression filtering threshold for assembly-based unigenes: a threshold transcripts per million (TPM) ≥ 5 removed more than 88% of assembly-based unigenes, which were mostly short and low-expressing unigenes. Second, assembly-based unigenes compensated for the incomplete length of PacBio CCSs: the ends of the 5`/3` untranslated regions of phenylpropanoid-related unigenes derived from set-1 were incomplete, which suggests that PacBio CCSs are unlikely to be full-length transcripts. Third, more isoform unigenes could be obtained from multiple unigene sets; isoform unigenes missing in Set-1 were detected in set-2 and set-3. Finally, gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses showed that phenylpropanoid biosynthesis and carbohydrate metabolism were highly activated in A. membranaceus roots. Various sequencing technologies and assemblers have been developed for de novo transcriptome analysis. However, no technique is perfect for de novo transcriptome analysis, suggesting the need to construct multiple unigene sets. This method enables efficient transcript filtering and detection of longer and more diverse transcripts.

Chromosome-Scale Genome Assembly and Triterpenoid Saponin Biosynthesis in Korean Bellflower (Platycodon grandiflorum).

Platycodon grandiflorum belongs to the Campanulaceae family and is an important medicinal and food plant in East Asia. However, on the whole, the genome evolution of P. grandiflorum and the molecular basis of its major biochemical pathways are poorly understood. We reported a chromosome-scale genome assembly of P. grandiflorum based on a hybrid method using Oxford Nanopore Technologies, Illumina sequences, and high-throughput chromosome conformation capture (Hi-C) analysis. The assembled genome was finalized as 574 Mb, containing 41,355 protein-coding genes, and the genome completeness was assessed as 97.6% using a Benchmarking Universal Single-Copy Orthologs analysis. The P. grandiflorum genome comprises nine pseudo-chromosomes with 56.9% repeat sequences, and the transcriptome analysis revealed an expansion of the 14 beta-amylin genes related to triterpenoid saponin biosynthesis. Our findings provide an understanding of P. grandiflorum genome evolution and enable genomic-assisted breeding for the mass production of important components such as triterpenoid saponins.

The chromosome-level genome assembly of lance asiabell (Codonopsis lanceolata), a medicinal and vegetable plant of the Campanulaceae family.

Codonopsis lanceolata (2n = 2x = 16) belongs to the Campanulaceae family and is a valuable medicinal and vegetable plant primarily found in East Asia. Several studies have demonstrated its excellent pharmacological effects, for example in bronchial treatment. However, genomic information of C. lanceolata is scarce, hindering studies on crop improvement of the species. Here, we report a high-quality chromosome-level genome assembly of C. lanceolata based on a hybrid method using Nanopore long-read, Illumina short-read, and Hi-C data. The assembled genome was completed as 1,273 Mb (84.5% of the estimated genome size), containing eight pseudo-chromosomes, ranging from 101.3 to 184.3 Mb. The genome comprised of 71.3% repeat sequences and 46,005 protein-coding genes, of which 85.7% genes were functionally annotated. Completeness of the assembled genome and genes was assessed to be 97.5% and 90.4%, respectively, by Benchmarking Universal Single-Copy Orthologs analysis. Phylogenetic and synteny analysis revealed that C. lanceolata was closely related to Platycodon grandiflorus in the Campanulaceae family. Gene family evolution revealed significant expansion of related genes involved in saponin biosynthesis in the C. lanceolata genome. This is the first reference genome reported for C. lanceolata. The genomic data produced in this study will provide essential information for further research to improve this medicinal plant and will broaden the understanding of the Campanulaceae family.

First Contiguous Genome Assembly of Japanese Lady Bell (Adenophora triphylla) and Insights into Development of Different Leaf Types.

Adenophora triphylla is an important medicinal and food plant found in East Asia. This plant is rich in secondary metabolites such as triterpenoid saponin, and its leaves can develop into different types, such as round and linear, depending on the origin of germination even within the same species. Despite this, few studies have comprehensively characterized the development processes of different leaf types and triterpenoid saponin pathways in this plant. Herein, we provide the first report of a high-quality genome assembly of A. triphylla based on a combination of Oxford Nanopore Technologies and Illumina sequencing methods. Its genome size was estimated to be 2.6 Gb, and the assembled genome finalized as 2.48 Gb, containing 57,729 protein-coding genes. Genome completeness was assessed as 95.6% using the Benchmarking Universal Single-Copy Orthologs score. The evolutionary divergence of A. triphylla was investigated using the genomes of five plant species, including two other species in the Campanulaceae family. The species A. triphylla diverged approximately 51-118 million years ago from the other four plants, and 579 expanded/contracted gene families were clustered in the Gene Ontology terms. The expansion of the ß-amyrin synthase (bAS) gene, a key enzyme in the triterpenoid saponin pathway, was identified in the A. triphylla genome. Furthermore, transcriptome analysis of the two leaf types revealed differences in the activity of starch, sucrose, unsaturated fatty acid pathways, and oxidoreductase enzymes. The heat and endoplasmic reticulum pathways related to plant stress were active in the development of round type leaf, while an enhancement of pyrimidine metabolism related to cell development was confirmed in the development of the linear type leaf. This study provides insight into the evolution of bAS genes and the development of different leaf types in A. triphylla.

Co-Expression Analysis Reveals Differential Expression of Homologous Genes Associated with Specific Terpenoid Biosynthesis in Rehmannia glutinosa.

Terpenoids are naturally occurring compounds involved in respiration, photosynthesis, membrane fluidity, and pathogen interactions and are classified according to the structure of their carbon skeleton. Although most terpenoids possess pharmacological activity, knowledge about terpenoid metabolism in medicinal plants is insufficient. Rehmannia glutinosa (R. glutinosa) is a traditional herb that is widely used in East Asia and has been reported to contain various terpenoids. In this study, we performed a comprehensive transcriptome analysis of terpenoid metabolism in R. glutinosa using two RNA sequencing platforms: Illumina and PacBio. The results show that the sterol, saponin, iridoid, and carotenoid pathways are active in R. glutinosa. Sterol and saponin biosynthesis were mevalonate pathway dependent, whereas iridoid and carotenoid biosynthesis were methylerythritol 4-phosphate pathway dependent. In addition, we found that the homologous genes of key enzymes involved in terpenoid metabolism were expressed differentially and that the differential expression of these genes was associated with specific terpenoid biosynthesis. The different expression of homologous genes encoding acetyl-CoA acetyltransferase, 3-hydroxy-3-methylglutaryl-CoA reductase, mevalonate kinase, mevalonate diphosphate decarboxylase, farnesyl pyrophosphate synthase, squalene synthase, and squalene epoxidase was associated with sterol and saponin biosynthesis. Homologous genes encoding 1-deoxy-D-xylulose 5-phosphate synthase were also differentially expressed and were associated with carotenoid and iridoid biosynthesis. These results suggest that the biosynthesis of specific terpenoids can be regulated by the homologous of key enzymes involved in plant terpenoid metabolism.

Analysis of Phenotypic Characteristics and Sucrose Metabolism in the Roots of Raphanus sativus L.

The taproot of radish (Raphanus sativus L.) is an important sink organ; it is morphologically diverse and contains large amounts of secondary metabolites. Sucrose metabolism is believed to be important in the development of sink organs. We measured the amounts of glucose, fructose, and sucrose in the roots of sixty three radish accessions and analyzed the association between the sugar content and the root phenotype. Fructose content correlated with the root color and length characteristics, glucose was the most abundant sugar in the roots, and the sucrose content was very low, compared to that of the hexoses in most of the accessions. Expression analysis of the genes involved in sucrose metabolism, transportation, starch synthesis, and cell wall synthesis was performed through RNA sequencing. The genes encoding sucrose synthases (SUSY) and the enzymes involved in the synthesis of cellulose were highly expressed, indicating that SUSY is involved in cell wall synthesis in radish roots. The positive correlation coefficient (R) between the sucrose content and the expression of cell wall invertase and sugar transporter proteins suggest that hexose accumulation could occur through the apoplastic pathway in radish roots. A positive R score was also obtained when comparing the expression of genes encoding SUSY and fructokinase (FK), suggesting that the fructose produced by SUSY is mostly phosphorylated by FK. In addition, we concluded that sucrose was the most metabolized sugar in radish roots.

Systemic Expression of Genes Involved in the Plant Defense Response Induced by Wounding in Senna tora.

Wounds in tissues provide a pathway of entry for pathogenic fungi and bacteria in plants. Plants respond to wounding by regulating the expression of genes involved in their defense mechanisms. To analyze this response, we investigated the defense-related genes induced by wounding in the leaves of Senna tora using RNA sequencing. The genes involved in jasmonate and ethylene biosynthesis were strongly induced by wounding, as were a large number of genes encoding transcription factors such as ERFs, WRKYs, MYBs, bHLHs, and NACs. Wounding induced the expression of genes encoding pathogenesis-related (PR) proteins, such as PR-1, chitinase, thaumatin-like protein, cysteine proteinase inhibitor, PR-10, and plant defensin. Furthermore, wounding led to the induction of genes involved in flavonoid biosynthesis and the accumulation of kaempferol and quercetin in S. tora leaves. All these genes were expressed systemically in leaves distant from the wound site. These results demonstrate that mechanical wounding can lead to a systemic defense response in the Caesalpinioideae, a subfamily of the Leguminosae. In addition, a co-expression analysis of genes induced by wounding provides important information about the interactions between genes involved in plant defense responses.

Machine learning, transcriptome, and genotyping chip analyses provide insights into SNP markers identifying flower color in Platycodon grandiflorus.

Bellflower is an edible ornamental gardening plant in Asia. For predicting the flower color in bellflower plants, a transcriptome-wide approach based on machine learning, transcriptome, and genotyping chip analyses was used to identify SNP markers. Six machine learning methods were deployed to explore the classification potential of the selected SNPs as features in two datasets, namely training (60 RNA-Seq samples) and validation (480 Fluidigm chip samples). SNP selection was performed in sequential order. Firstly, 96 SNPs were selected from the transcriptome-wide SNPs using the principal compound analysis (PCA). Then, 9 among 96 SNPs were later identified using the Random forest based feature selection method from the Fluidigm chip dataset. Among six machines, the random forest (RF) model produced higher classification performance than the other models. The 9 SNP marker candidates selected for classifying the flower color classification were verified using the genomic DNA PCR with Sanger sequencing. Our results suggest that this methodology could be used for future selection of breeding traits even though the plant accessions are highly heterogeneous.

Induction of Glucoraphasatin Biosynthesis Genes by MYB29 in Radish (Raphanus sativus L.) Roots.

Glucoraphasatin (GRH) is a specific aliphatic glucosinolate (GSL) that is only abundant in radish (Raphanus sativus L.). The gene expression regulating GRH biosynthesis in radish is still poorly understood. We employed a total of 59 radish accessions to analyze GSL profiles and showed that GRH was specific and predominant among the aliphatic GSLs in radish roots. We selected five accessions roots with high, moderate and low GSL biosynthesis, respectively, to conduct a comparative transcriptome analysis and the qRT-PCR of the biosynthesis genes for aliphatic GSLs. In this study, among all the accessions tested, roots with the accession RA157-74 had a high GRH content and showed a significant expression of the aliphatic GSL biosynthesis genes. We defined the genes involved in the GRH biosynthesis process and found that they were regulated by a transcription factor (RSG00789) at the MYB29 locus in radish roots. We found 13 aliphatic GSL biosynthesis genes regulated by the RSG00789 gene in the GRH biosynthesis pathway.

Compositional analyses of diverse phytochemicals and polar metabolites from different-colored potato (Solanum tubersum L.) tubers.

Lipophilic bioactive compounds and hydrophilic primary metabolites from potato (solanum tubersum L.) tubers with different-colored flesh (white-, yellow-, red-, and purple) were characterized. The carotenoid content was relatively higher in red-colored potatoes, in which lutein was most plentiful. Among the other lipophilic compounds analyzed, including policosanols, tocopherols, and phytosterols, octacosanol was measured in the largest amount, followed by ß-sitosterol, irrespective of color variations. Forty-three hydrophilics consisting of amino acids, organic acids, sugars, and sugar alcohols and 18 lipophilics were subjected to data-mining processes. The results of multivariate statistical analyses clearly distincted the different varieties and separated red-fleshed potatoes from other color-fleshed potatoes according to abundance of amino acids, sugars, and carotenoids. This study confirmed the metabolic association-related biochemical pathway between metabolite characteristic and color differences in potato tubers. These results can facilitate understanding the metabolic differences among diverse colored potatoes and provide fruitful information for genetic engineering of potato cultivars.

Comparative evaluation of nutritional compositions between transgenic rice harboring the CaMsrB2 gene and the conventional counterpart.

As a part of a safety assessment of new transgenic crops, compositional equivalence studies between transgenic crops with non-transgenic comparators are almost universally required. This study was conducted to compare nutritional profiles of proximate composition, and fatty acid, amino acid, mineral, and vitamin contents, and anti-nutrients, between transgenic drought-tolerant Agb0103 rice harboring the pepper methionine sulfoxide reductase B2 gene CaMsrB2 and the parental rice cultivar, 'Ilmi' as a non-transgenic control. Both transgenic and non-transgenic rice were grown and harvested in 2 different locations. Proximate compositions of moisture, starch, protein, lipid, and ash content of Agb0103 rice were similar to parental non-transgenic rice. There were no differences between transgenic and non-transgenic rice with respect to the whole nutritional composition, except for minor locality differences for a few nutritional components. Agb0103 rice with improved resistance to drought is nutritionally equivalent to the parental rice cultivar.

Comparative analysis of phytochemicals and polar metabolites from colored sweet potato (Ipomoea batatas L.) tubers.

We determined the phytochemical diversity, including carotenoids, flavonoids, anthocyanins, and phenolic acids, in sweet potatoes (Ipomoea batatas L.) with distinctive flesh colors (white, orange, and purple) and identified hydrophilic primary metabolites. Carotenoid content was considerably higher in orange-fleshed sweet potatoes, wherein ß-carotene was the most plentiful, and anthocyanins were detected only in purple-fleshed sweet potatoes. The levels of phenolic acids and flavonoids were relatively higher in purple-fleshed sweet potatoes than those in the other two varieties. Forty-one primary and 18 secondary metabolite profiles were subjected to multivariate statistical analyses, which fully distinguished among the varieties and separated orange- and purple-fleshed sweet potatoes from white-fleshed sweet potatoes based on the high levels of sugars, sugar alcohols, and secondary metabolites. This is the first study to determine comprehensive metabolic differences among different color-fleshed sweet potatoes and provides useful information for genetic manipulation of sweet potatoes to influence primary and secondary metabolism.

Proteomic identification of a novel toxin protein (Txp40) from Xenorhabdus nematophila and its insecticidal activity against larvae of Plutella xylostella.

For the identification of a novel insecticidal protein, a two-dimensional liquid chromatography (PF-2D) system was used in a quantitative proteomic analysis of Xenorhabdus nematophila CBNU strain isolated from entomophagous nematode Steinernema carpocapsae . Protein patterns obtained from minimum and maximum insecticidal activities during cultivation were contrasted, and a novel toxin protein (Txp40) was identified by MALDI-TOF/MS. The DNA sequence of the cloned toxin gene (1089 bp) has an open reading frame encoding 363 amino acids with a predicted molecular mass of 41162 Da. The txp40 identified in this study is most closely related to the known txp40 cloned from X. nematophila EB (ADQ92844) with 94.4% identical sequence residues. Following the expression of the newly identified toxin gene in Escherichia coli , the insecticidal activity of the recombinant toxin protein was determined against Plutella xylostella larvae; a 56.7% mortality rate was observed within 24 h.

Exploiting leaf starch synthesis as a transient sink to elevate photosynthesis, plant productivity and yields.

Improvements in plant productivity (biomass) and yield have centered on increasing the efficiency of leaf CO(2) fixation and utilization of products by non-photosynthetic sink organs. We had previously demonstrated a correlation between photosynthetic capacity, plant growth, and the extent of leaf starch synthesis utilizing starch-deficient mutants. This finding suggested that leaf starch is used as a transient photosynthetic sink to recycle inorganic phosphate and, in turn, maximize photosynthesis. To test this hypothesis, Arabidopsis thaliana and rice (Oryza sativa L.) lines were generated with enhanced capacity to make leaf starch with minimal impact on carbon partitioning to sucrose. The Arabidopsis engineered plants exhibited enhanced photosynthetic capacity; this translated into increased growth and biomass. These enhanced phenotypes were displayed by similarly engineered rice lines. Manipulation of leaf starch is a viable alternative strategy to increase photosynthesis and, in turn, the growth and yields of crop and bioenergy plants.

Expression of PAT and NPT II proteins during the developmental stages of a genetically modified pepper developed in Korea.

Estimation of the protein levels introduced in a biotechnology-derived product is conducted as part of an overall safety assessment. An enzyme-linked immunosorbent assay (ELISA) was used to analyze phosphinothricin acetyltransferase (PAT) and neomycin phosphotransferase II (NPT II) protein expression in a genetically modified (GM) pepper plant developed in Korea. PAT and NPT II expression levels, based on both dry weight and fresh weight, were variable among different plant generations and plant sections from isolated genetically modified organism (GMO) fields at four developmental stages. PAT expression was highest in leaves at anthesis (11.44 µg/gdw and 2.17 µg/gfw) and lowest in roots (0.12 µg/gdw and 0.01 µg/gfw). NPT II expression was also highest in leaves at anthesis (17.31 µg/gdw and 3.41 µg/gfw) and lowest in red pepper (0.65 µg/gdw and 0.12 µg/gfw). In pollen, PAT expression was 0.59-0.62 µg/gdw, while NPT II was not detected. Both PAT and NPT II showed a general pattern of decreased expression with progression of the growing season. As expected, PAT and NPT II protein expression was not detectable in control pepper plants.

Auto-excision of selectable marker genes from transgenic tobacco via a stress inducible FLP/FRT site-specific recombination system.

Antibiotic resistance marker genes are powerful selection tools for use in plant transformation processes. However, once transformation is accomplished, the presence of these resistance genes is no longer necessary and can even be undesirable. We herein describe the successful excision of antibiotic resistance genes from transgenic plants via the use of an oxidative stress-inducible FLP gene. FLP encodes a recombinase that can eliminate FLP and hpt selection genes flanked by two FRT sites. During a transformation procedure in tobacco, transformants were obtained by selection on hygromycin media. Regenerants of the initial transformants were screened for selective marker excision in hydrogen peroxide supplemented media and both the FLP and hpt genes were found to have been eliminated. About 13-41% of regenerated shoots on hydrogen peroxide media were marker-free. This auto-excision system, mediated by the oxidative stress-inducible FLP/FRT system to eliminate a selectable marker gene can be very readily adopted and used to efficiently generate marker-free transgenic plants.

A sepal-expressed ADP-glucose pyrophosphorylase gene (NtAGP) is required for petal expansion growth in 'Xanthi' tobacco.

In this study, a tobacco (Nicotiana tabacum 'Xanthi') ADP-glucose pyrophosphorylase cDNA (NtAGP) was isolated from a flower bud cDNA library and the role of NtAGP in the growth of the floral organ was characterized. The expression of NtAGP was high in the sepal, moderate in the carpel and stamen, and low in the petal tissues. NtAGP-antisense plants produced flowers with abnormal petal limbs due to the early termination of the expansion growth of the petal limbs between the corolla lobes. Microscopic observation of the limb region revealed that cell expansion was limited in NtAGP-antisense plants but that cell numbers remained unchanged. mRNA levels of NtAGP, ADP-glucose pyrophosphorylase activity, and starch content in the sepal tissues of NtAGP-antisense plants were reduced, resulting in significantly lower levels of sugars (sucrose, glucose, and fructose) in the petal limbs. The feeding of these sugars to flower buds of the NtAGP-antisense plants restored the expansion growth in the limb area between the corolla lobes. Expansion growth of the petal limb between the corolla lobes was severely arrested in 'Xanthi' flowers from which sepals were removed, indicating that sepal carbohydrates are essential for petal limb expansion growth. These results demonstrate that NtAGP plays a crucial role in the morphogenesis of petal limbs in 'Xanthi' through the synthesis of starch, which is the main carbohydrate source for expansion growth of petal limbs, in sepal tissues.

Ribosome-inactivating activity and cDNA cloning of antiviral protein isoforms of Chenopodium album.

We have characterized a novel type I ribosome-inactivating protein (CAP30) from the leaves of Chenopodium album. Purified native CAP30 depurinated the ribosomes of Chenopodium, tomato, and tobacco leaves in vitro. To further characterize this protein, cDNA clones were isolated from a leaf cDNA library using a DNA probe derived from the N-terminal amino acid sequence. Two full-length cDNA clones, CAP30A and CAP30B, were isolated. The two clones were highly homologous (91.4% identity over 280 amino acids) at the deduced amino acid level. Both contain a putative signal peptide of 25 amino acid and a conserved domain commonly found in ribosome-inactivating proteins. This suggests that CAP30 is a single-chain ribosome-inactivating protein. Expression of CAP30 mRNA peaked twice, at 12 and 72 h, after tobacco mosaic virus (TMV) infection or wounding. Transformed Escherichia coli cells expressing pre- or mature CAP had greatly reduced growth rates. These results suggest that CAP30 functions as a broad-spectrum defense-related protein with both antiviral and anti-microbial activity.