Full-length transcriptome sequences and the identification of putative genes for flavonoid biosynthesis in safflower

The Carthamus tinctorius L. genome sequence provides insights into synthesis of unsaturated fatty acids

Domesticated safflower (Carthamus tinctorius L.) is a widely cultivated edible oil crop. However, despite its economic importance, the genetic basis underlying key traits such as oil content, resistance to biotic and abiotic stresses, and flowering time remains poorly understood. Here, we present the genome assembly for C. tinctorius variety Jihong01, which was obtained by integrating Oxford Nanopore Technologies (ONT) and BGI-SEQ500 sequencing results. The assembled genome was 1,061.1 Mb, and consisted of 32,379 protein-coding genes, 97.71% of which were functionally annotated. Safflower had a recent whole genome duplication (WGD) event in evolution history and diverged from sunflower approximately 37.3 million years ago. Through comparative genomic analysis at five seed development stages, we unveiled the pivotal roles of fatty acid desaturase 2 (FAD2) and fatty acid desaturase 6 (FAD6) in linoleic acid (LA) biosynthesis. Similarly, the differential gene expression analysis further reinforced the significance of these genes in regulating LA accumulation. Moreover, our investigation of seed fatty acid composition at different seed developmental stages unveiled the crucial roles of FAD2 and FAD6 in LA biosynthesis. These findings offer important insights into enhancing breeding programs for the improvement of quality traits and provide reference resource for further research on the natural properties of safflower.

Integrating Full-Length and Second-Generation Transcriptomics to Reveal Differentially Expressed Genes Associated with the Development of Corydalis yanhusuo Tuber

Corydalis yanhusuo is a medicinal herb in China that has been widely used to treat various kinds of pain. The tuber is the main organ of C. yanhusuo used for medicinal purposes, but changes in related genes during the development of the tuber have rarely been reported. To identify the differentially expressed genes during tuber development, C. yanhusuo full-length transcriptomic sequencing was performed using single-molecule real-time technology, and tubers at three development stages were selected for comparative transcriptome analysis. A total of 90,496 full-length non-chimeric transcripts were obtained, and 19,341 transcripts were annotated in at least one public database. A total of 9221 differentially expressed genes were identified during the swelling process of C. yanhusuo tuber. A Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis revealed that differentially expressed genes associated with a "starch and sucrose metabolism pathway", "phenylpropanoid biosynthesis pathway", "isoquinoline alkaloid biosynthesis pathway", "zeatin biosynthesis pathway", and "brassinosteroid biosynthesis pathway" were predominantly enriched. In addition, the genes involved in cell wall metabolism were potentially associated with tuber swelling. These processes regulated and were involved in C. yanhusuo tuber development. The results provide a foundation for further research on tuber formation in medicinal plants.

Identification and Characterization of CtUGT3 as the Key Player of Astragalin Biosynthesis in Carthamus tinctorius L

Safflower (Carthamus tinctorius L.) is a multipurpose economic crop that is distributed worldwide. Flavonoid glycosides are the main bioactive components in safflower, but only a few UDP-glycosyltransferases (UGT) have been identified. Three differentially expressed UGT genes related with the accumulation of 9 flavonoid O-glycosides were screened from metabolomics and transcriptome analysis. Safflower corolla protoplasts were used to confirm the glycosylation ability of UGT candidates in vivo for the first time. The astragalin content was significantly increased only when CtUGT3 was overexpressed. CtUGT3 also showed flavonoid 3-OH and 7-OH glycosylation activities in vitro. Molecular modeling and site-directed mutagenesis revealed that G15, T136, S276, and E384 were critical catalytic residues for the glycosylation ability of CtUGT3. These results demonstrate that CtUGT3 has a flavonoid 3-OH glycosylation function and is involved in the biosynthesis of astragalin in safflower. This study provides a reference for flavonoid biosynthesis genes research in nonmodel plants.

Full-Length Transcriptome Characterization and Functional Analysis of Pathogenesis-Related Proteins in Lilium Oriental Hybrid 'Sorbonne' Infected with Botrytis elliptica

Gray mold (Botrytis elliptica) causes a deleterious fungal disease that decreases the ornamental value and yield of lilies. Lilium oriental hybrid 'Sorbonne' is a variety that is resistant to gray mold. Understanding the mechanism of resistance against B. elliptica infection in 'Sorbonne' can provide a basis for the genetic improvement in lily plants. In this study, a PacBio Sequel II system was used to sequence the full-length transcriptome of Lilium 'Sorbonne' after inoculation with B. elliptica. A total of 46.64 Gb subreads and 19,102 isoforms with an average length of 1598 bp were obtained. A prediction analysis revealed 263 lncRNAs, and 805 transcription factors, 4478 simple sequence repeats, and 17,752 coding sequences were identified. Pathogenesis-related proteins (PR), which may play important roles in resistance against B. elliptica infection, were identified based on the full-length transcriptome data and previously obtained second-generation transcriptome data. Nine non-redundant potential LhSorPR proteins were identified and assigned to two groups that were composed of two LhSorPR4 and seven LhSorPR10 proteins based on their genetic relatedness. The real-time quantitative reverse transcription PCR (qRT-PCR) results showed that the patterns of expression of nine differentially expressed PR genes under B. elliptica stress were basically consistent with the results of transcriptome sequencing. The pattern of expression of LhSorPR4s and LhSorPR10s genes in different tissues was analyzed, and the expression of each gene varied. Furthermore, we verified the function of LhSorPR4-2 gene in Lilium. The expression of LhSorPR4-2 was induced by phytohormones such as methyl jasmonate, salicylic acid, and ethephon. Moreover, the promoter region of LhSorPR4-2 was characterized by several functional domains associated with phytohormones and stress response. The overexpression of LhSorPR4-2 gene in 'Sorbonne' increased the resistance of the lily plant to B. elliptica and correlated with high chitinase activity. This study provides a full-length transcript database and functionally analyzed the resistance of PR gene to B. elliptica in Lilium, thereby introducing the candidate gene LhSorPR4-2 to breed resistance in Lilium.

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.

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis

BACKGROUND

Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis.

RESULTS

We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts-twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage.

CONCLUSIONS

AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

Genome-wide transcriptome analysis and characterization of the cytochrome P450 flavonoid biosynthesis genes in pigeon pea (Cajanus cajan)

226 CcCYP450 genes were identified at the genomic level and were classified into 45 clades based on phylogenetic analysis. CcCYP75B165 gene was found that might play important roles in the biosynthesis of flavonoids in pigeon pea, and was significantly induced by methyl jasmonate (MeJA). The cytochrome P450 mono-oxygenase (CYP450) superfamily plays a key role in the flavonoid biosynthesis pathway and resists different kinds of stresses. Several CYP450 genes have been identified to be involved in the biosynthesis of crop protection agents. However, the CcCYP450 genes from pigeon pea have not been identified. Here, 226 CcCYP450 genes were identified at the genomic level by analysing the gene structure, distribution on chromosomes, gene duplication, and conserved motifs and were classified into 45 clades based on phylogenetic analysis. RNA-seq analysis revealed clear details of CcCYP450 genes that varied with time of MeJA (methyl jasmonate) induction. Among them, six CcCYP450 subfamily genes were found that might play important roles in the biosynthesis of flavonoids in pigeon pea. The overexpression of CcCYP75B165 in pigeon pea significantly induced the accumulation of genistin and downregulated the contents of cajaninstilbene acid, apigenin, isovitexin, and genistein and the expression of flavonoid synthase genes. This study provides theoretical guidance and plant genetic resources for cultivating new pigeon pea varieties with high flavonoid contents and exploring the molecular mechanisms of the biosynthesis of flavonoids under MeJA treatment.

Pangenomes as a Resource to Accelerate Breeding of Under-Utilised Crop Species

Pangenomes are a rich resource to examine the genomic variation observed within a species or genera, supporting population genetics studies, with applications for the improvement of crop traits. Major crop species such as maize (Zea mays), rice (Oryza sativa), Brassica (Brassica spp.), and soybean (Glycine max) have had pangenomes constructed and released, and this has led to the discovery of valuable genes associated with disease resistance and yield components. However, pangenome data are not available for many less prominent crop species that are currently under-utilised. Despite many under-utilised species being important food sources in regional populations, the scarcity of genomic data for these species hinders their improvement. Here, we assess several under-utilised crops and review the pangenome approaches that could be used to build resources for their improvement. Many of these under-utilised crops are cultivated in arid or semi-arid environments, suggesting that novel genes related to drought tolerance may be identified and used for introgression into related major crop species. In addition, we discuss how previously collected data could be used to enrich pangenome functional analysis in genome-wide association studies (GWAS) based on studies in major crops. Considering the technological advances in genome sequencing, pangenome references for under-utilised species are becoming more obtainable, offering the opportunity to identify novel genes related to agro-morphological traits in these species.

Unveiling the transcriptomic complexity of Miscanthus sinensis using a combination of PacBio long read- and Illumina short read sequencing platforms

Background

Miscanthus sinensis Andersson is a perennial grass that exhibits remarkable lignocellulose characteristics suitable for sustainable bioenergy production. However, knowledge of the genetic resources of this species is relatively limited, which considerably hampers further work on its biology and genetic improvement.

Results

In this study, through analyzing the transcriptome of mixed samples of leaves and stems using the latest PacBio Iso-Seq sequencing technology combined with Illumina HiSeq, we report the first full-length transcriptome dataset of M. sinensis with a total of 58.21 Gb clean data. An average of 15.75 Gb clean reads of each sample were obtained from the PacBio Iso-Seq system, which doubled the data size (6.68 Gb) obtained from the Illumina HiSeq platform. The integrated analyses of PacBio- and Illumina-based transcriptomic data uncovered 408,801 non-redundant transcripts with an average length of 1,685 bp. Of those, 189,406 transcripts were commonly identified by both methods, 169,149 transcripts with an average length of 619 bp were uniquely identified by Illumina HiSeq, and 51,246 transcripts with an average length of 2,535 bp were uniquely identified by PacBio Iso-Seq. Approximately 96 % of the final combined transcripts were mapped back to the Miscanthus genome, reflecting the high quality and coverage of our sequencing results. When comparing our data with genomes of four species of Andropogoneae, M. sinensis showed the closest relationship with sugarcane with up to 93 % mapping ratios, followed by sorghum with up to 80 % mapping ratios, indicating a high conservation of orthologs in these three genomes. Furthermore, 306,228 transcripts were successfully annotated against public databases including cell wall related genes and transcript factor families, thus providing many new insights into gene functions. The PacBio Iso-Seq data also helped identify 3,898 alternative splicing events and 2,963 annotated AS isoforms within 10 function categories.

Conclusions

Taken together, the present study provides a rich data set of full-length transcripts that greatly enriches our understanding of M. sinensis transcriptomic resources, thus facilitating further genetic improvement and molecular studies of the Miscanthus species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07971-x.

The chromosome-scale reference genome of safflower (Carthamus tinctorius) provides insights into linoleic acid and flavonoid biosynthesis

Safflower (Carthamus tinctorius L.), a member of the Asteraceae, is a popular crop due to its high linoleic acid (LA) and flavonoid (such as hydroxysafflor yellow A) contents. Here, we report the first high-quality genome assembly (contig N50 of 21.23 Mb) for the 12 pseudochromosomes of safflower using single-molecule real-time sequencing, Hi-C mapping technologies and a genetic linkage map. Phyloge nomic analysis showed that safflower diverged from artichoke (Cynara cardunculus) and sunflower (Helianthus annuus) approximately 30.7 and 60.5 million years ago, respectively. Comparative genomic analyses revealed that uniquely expanded gene families in safflower were enriched for those predicted to be involved in lipid metabolism and transport and abscisic acid signalling. Notably, the fatty acid desaturase 2 (FAD2) and chalcone synthase (CHS) families, which function in the LA and flavonoid biosynthesis pathways, respectively, were expanded via tandem duplications in safflower. CarFAD2-12 was specifically expressed in seeds and was vital for high-LA content in seeds, while tandemly duplicated CarFAD2 genes were up-regulated in ovaries compared to CarFAD2-12, which indicates regulatory divergence of FAD2 in seeds and ovaries. CarCHS1, CarCHS4 and tandem-duplicated CarCHS5˜CarCHS6, which were up-regulated compared to other CarCHS members at early stages, contribute to the accumulation of major flavonoids in flowers. In addition, our data reveal multiple alternative splicing events in gene families related to fatty acid and flavonoid biosynthesis. Together, these results provide a high-quality reference genome and evolutionary insights into the molecular basis of fatty acid and flavonoid biosynthesis in safflower.

The Full-Length Transcriptome Sequencing and Identification of Na+/H+ Antiporter Genes in Halophyte Nitraria tangutorum Bobrov

Nitraria tangutorum Bobrov is a halophyte that is resistant to salt and alkali and is widely distributed in northwestern China. However, its genome has not been sequenced, thereby limiting studies on this particular species. For species without a reference genome, the full-length transcriptome is a convenient and rapid way to obtain reference gene information. To better study N. tangutorum, we used PacBio single-molecule real-time technology to perform full-length transcriptome analysis of this halophyte. In this study, a total of 21.83 Gb of data were obtained, and 198,300 transcripts, 51,875 SSRs (simple sequence repeats), 55,574 CDS (coding sequence), and 74,913 lncRNAs (long non-coding RNA) were identified. In addition, using this full-length transcriptome, we identified the key Na⁺/H⁺ antiporter (NHX) genes that maintain ion balance in plants and found that these are induced to express under salt stress. The results indicate that the full-length transcriptome of N. tangutorum can be used as a database and be utilized in elucidating the salt tolerance mechanism of N. tangutorum.

Transcriptomic and metabolomic analyses provide insights into the biosynthesis of chlorogenic acids in Lonicera macranthoides Hand.-Mazz

Lonicera macranthoides Hand.-Mazz (L. macranthoides) is a medicinal herb that is widely distributed in South China. The developmental stage and corolla dehiscence of the flower are the important factors affecting the quality of medicinal ingredients. However, neither the regulatory mechanism controlling chlorogenic acids biosynthesis in L. macranthoides nor the molecular basis of effect of corolla dehiscence on the quality of medicinal materials is fully understood. In this study, metabolomics and transcriptomics were used to analyze the metabolic and transcriptional differences of two different cultivars closed bud type (Bt), and flowering type (Ft), as well as the effect of jasmonic acid methyl ester (MeJA) on chlorogenic acids (CGAs) biosynthesis. In total, large number of differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) were filtered among three lines of samples. Gene metabolite correlation analyses revealed a ‘core set’ of 30 genes and 54 genes that were strongly correlated with CGAs biosynthesis and regulating the flowering, respectively. Quantitative real-time polymerase chain reaction results proved the alterations in the expression levels of genes encoding the pathways involved in CGAs biosynthesis. The ion abundances of CGAs were most significantly increased, while some of the CGAs derived and Caffeoyl-CoA-derived substances showed the most largely reduced abundances in the closed bud type (Bt) compared to the flowering type (Ft). MeJA may leads to the activation of downstream genes in CGAs biosynthesis pathway. Overall, there were significant differences in the transcriptional and metabolic levels of CGAs biosynthesis pathway in flower buds of different flowering cultivars. The redirection of metabolic flux may contribute to increased accumulation of CGAs. However, whether MeJA and flowering have direct effects on the accumulation of CGAs needs further studied. These researches effectively expanded the functional genomic library and provide new insights into CGAs biosynthesis in L. macranthoides.

PacBio Single-Molecule Long-Read Sequencing Reveals Genes Tolerating Manganese Stress in Schima superba Saplings

Schima superba (Theaceae) is a subtropical evergreen tree and is used widely for forest firebreaks and gardening. It is a plant that tolerates salt and typically accumulates elevated amounts of manganese in the leaves. With large ecological amplitude, this tree species grows quickly. Due to its substantial biomass, it has a great potential for soil remediation. To evaluate the thorough framework of the mRNA, we employed PacBio sequencing technology for the first time to generate S. Superba transcriptome. In this analysis, overall, 511,759 full length non-chimeric reads were acquired, and 163,834 high-quality full-length reads were obtained. Overall, 93,362 open reading frames were obtained, of which 78,255 were complete. In gene annotation analyses, the Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Genes (COG), Gene Ontology (GO), and Non-Redundant (Nr) databases were allocated 91,082, 71,839, 38,914, and 38,376 transcripts, respectively. To identify long non-coding RNAs (lncRNAs), we utilized four computational methods associated with protein families (Pfam), Cooperative Data Classification (CPC), Coding Assessing Potential Tool (CPAT), and Coding Non-Coding Index (CNCI) databases and observed 8,551, 9,174, 20,720, and 18,669 lncRNAs, respectively. Moreover, nine genes were randomly selected for the expression analysis, which showed the highest expression of Gene 6 (Na_Ca_ex gene), and CAX (CAX-interacting protein 4) was higher in manganese (Mn)-treated group. This work provided significant number of full-length transcripts and refined the annotation of the reference genome, which will ease advanced genetic analyses of S. superba.

Full-Length Transcriptome Sequencing Provides Insights into Flavonoid Biosynthesis in Fritillaria hupehensis

One of the most commonly utilized medicinal plants in China is Fritillaria hupehensis (Hsiao et K.C. Hsia). However, due to a lack of genomic resources, little is known about the biosynthesis of relevant compounds, particularly the flavonoid biosynthesis pathway. A PacBio RS II sequencing generated a total of 342,044 reads from the bulb, leaf, root, and stem, of which 316,438 were full-length (FL) non-redundant reads with an average length of 1365 bp and a N50 of 1888 bp. There were also 38,607 long non-coding RNAs and 7914 simple sequence repeats detected. To improve our understanding of processes implicated in regulating secondary metabolite biosynthesis in F. hupehensis tissues, we evaluated potential metabolic pathways. Overall, this study provides a repertoire of FL transcripts in F. hupehensis for the first time, and it will be a valuable resource for marker-assisted breeding and research into bioactive compounds for medicinal and pharmacological applications.

Integrated Metabolome and Transcriptome Analysis Unveils Novel Pathway Involved in the Formation of Yellow Peel in Cucumber

Yellow peel will adversely affect the appearance quality of cucumber fruit, but the metabolites and the molecular mechanism of pigment accumulation in cucumber peel remain unclear. Flavonoid metabolome and transcriptome analyses were carried out on the young peel and old peel of the color mutant L19 and the near-isogenic line L14. The results showed that there were 165 differential flavonoid metabolites in the old peel between L14 and L19. The total content of representative flavonoid metabolites in the old peel of L14 was 95 times that of L19, and 35 times that of young peel of L14, respectively. This might explain the difference of pigment accumulation in yellow peel. Furthermore, transcriptome analysis showed that there were 3396 and 1115 differentially expressed genes in the yellow color difference group (Young L14 vs. Old L14 and Old L14 vs. Old L19), respectively. These differentially expressed genes were significantly enriched in the MAPK signaling pathway-plant, plant-pathogen interaction, flavonoid biosynthesis and cutin, suberine and wax biosynthesis pathways. By analyzing the correlation between differential metabolites and differentially expressed genes, six candidate genes related to the synthesis of glycitein, kaempferol and homoeriodictyol are potentially important. In addition, four key transcription factors that belong to R2R3-MYB, bHLH51 and WRKY23 might be the major drivers of transcriptional changes in the peel between L14 and L19. Then, the expression patterns of these important genes were confirmed by qRT-PCR. These results suggested that the biosynthesis pathway of homoeriodictyol was a novel way to affect the yellowing of cucumber peel. Together, the results of this study provide a research basis for the biosynthesis and regulation of flavonoids in cucumber peel and form a significant step towards identifying the molecular mechanism of cucumber peel yellowing.

Integrated Metabolomics and Transcriptome Analysis of Flavonoid Biosynthesis in Safflower (Carthamus tinctorius L.) With Different Colors

Safflower is widely used in dying and in traditional medicine, and C-glucosylquinochalcones are the main metabolic species in the red color of safflower. Various safflower cultivars have flowers with different colors. However, the metabolic and transcriptional differences among safflower cultivars with different-colored flowers and the genes participating in C-glucosylquinochalcone biosynthesis are largely unknown. To provide insights on this issue, we performed integrated metabolomics and transcriptome analyses on the flavonoid biosynthesis of flowers of different colors in safflower (white-W, yellow-Y, light red-LR, and deep red-DR). The metabolic analysis showed that flavonoid metabolites showed great differences among the different colors of safflower. More flavonoid metabolic species were detected in Y and W, while C-glucosylquinochalcones were not detected in W. The content of C-glucosylquinochalcones increased with increasing color. Transcriptional analysis showed that most of the annotated flavonoid biosynthesis genes were significantly increased in W. The expression of genes related to flavonoid biosynthesis decreased with increasing color. We analyzed the candidate genes associated with C-glucosylquinochalcones, and an integration of the metabolic and transcriptional analyses indicated that the differential expression of the chalcone synthase (CHS) gene is one of the main reasons for the difference in flavonoid species and content among the different colors of safflower. Combined with the expression pattern analysis, these results indicated that HH_035319, HH_032689, and HH_018025 are likely involved in C-glucosylquinochalcones biosynthesis. In addition, we found that their expression showed greatly increased after the methyl jasmonate (MeJA) treatment. Therefore, HH_035319, HH_032689, and HH_018025 might participate in C-glucosylquinochalcone biosynthesis, which ultimately leads to the red color in safflower.

Full-length transcriptome analysis and identification of genes involved in asarinin and aristolochic acid biosynthesis in medicinal plant Asarum sieboldii

Asarum sieboldii, a well-known traditional Chinese medicinal herb, is used for curing inflammation and ache. It contains both the bioactive ingredient asarinin and the toxic compound aristolochic acid. To address further breeding demand, genes involved in the biosynthetic pathways of asarinin and aristolochic acid should be explored. Therefore, the full-length transcriptome of A. sieboldii was sequenced using PacBio Iso-Seq to determine the candidate transcripts that encode the biosynthetic enzymes of asarinin and aristolochic acid. In this study, 63 023 full-length transcripts were generated with an average length of 1371 bp from roots, stems, and leaves, of which 49 593 transcripts (78.69%) were annotated against public databases. Furthermore, 555 alternative splicing (AS) events, 10 869 long noncoding RNAs (lncRNAs) as well as their 11 291 target genes, and 17 909 simple sequence repeats (SSRs) were identified. The data also revealed 97 candidate transcripts related to asarinin metabolism, of which six novel genes that encoded enzymes involved in asarinin biosynthesis were initially reported. In addition, 56 transcripts related to aristolochic acid biosynthesis were also identified, including CYP81B. In summary, these transcriptome data provide a useful resource to study gene function and genetic engineering in A. sieboldii.

High quality genome of Erigeron breviscapus provides a reference for herbal plants in Asteraceae

Erigeron breviscapus is an important medicinal plant in Compositae and the first species to realize the whole process from the decoding of the draft genome sequence to scutellarin biosynthesis in yeast. However, the previous low‐quality genome assembly has hindered the optimization of candidate genes involved in scutellarin synthesis and the development of molecular‐assisted breeding based on the genome. Here, the E. breviscapus genome was updated using PacBio RSII sequencing data and Hi‐C data, and increased in size from 1.2 Gb to 1.43 Gb, with a scaffold N50 of 156.82 Mb and contig N50 of 140.95 kb, and a total of 43,514 protein‐coding genes were obtained and oriented onto nine pseudo‐chromosomes, thus becoming the third plant species assembled to chromosome level after sunflower and lettuce in Compositae. Fourteen genes with evidence for positive selection were identified and found to be related to leaf morphology, flowering and secondary metabolism. The number of genes in some gene families involved in flavonoid biosynthesis in E. breviscapus have been significantly expanded. In particular, additional candidate genes involved in scutellarin biosynthesis, such as flavonoid‐7‐O‐glucuronosyltransferase genes (F7GATs) were identified using updated genome. In addition, three candidate genes encoding indole‐3‐pyruvate monooxygenase YUCCA2 (YUC2), serine carboxypeptidase‐like 18 (SCPL18), and F‐box protein (FBP), respectively, were identified to be probably related to leaf development and flowering by resequencing 99 individuals. These results provided a substantial genetic basis for improving agronomic and quality traits of E. breviscapus, and provided a platform for improving other draft genome assemblies to chromosome‐level.

Genome-wide identification and expression analysis of bZIP gene family in Carthamus tinctorius L

The basic leucine zipper (bZIP) is a widely known transcription factors family in eukaryotes. In plants, the role of bZIP proteins are crucial in various biological functions such as plant growth and development, seed maturation, response to light signal and environmental stress. To date, bZIP protein family has been comprehensively identified in Arabidopsis, castor, rice, ramie, soybean and other plant species, however, the complete genome-wide investigation of Carthamus tinctorius-bZIP family still remains unexplained. Here, we identified 52 putative bZIP genes from Carthamus tinctorius using a draft genome assembly and further analyzed their evolutionary classification, physicochemical properties, Conserved domain analysis, functional differentiation and the investigation of expression level in different tissues. Based on the common bZIP domain, CtbZIP family were clustered into 12 subfamilies renamed as (A-J, S, X), of which the X is a unique subfamily to Carthamus tinctorius. A total of 20 conserved protein motifs were found in CtbZIP proteins. The expression profiling of CtbZIP genes deciphered their tissue-specific pattern. Furthermore, the changes in CtbZIP transcript abundance suggested that their transcription regulation could be highly influenced by light intensity and hormones. Collectively, this study highlights all functional and regulatory elements of bZIP transcription factors family in Carthamus tinctorius which may serve as potential candidates for functional characterization in future.

Comparative transcriptome analysis of root, stem, and leaf tissues of Entada phaseoloides reveals potential genes involved in triterpenoid saponin biosynthesis

Background

Entada phaseoloides (L.) Merr. is an important traditional medicinal plant. The stem of Entada phaseoloides is popularly used as traditional medicine because of its significance in dispelling wind and dampness and remarkable anti-inflammatory activities. Triterpenoid saponins are the major bioactive compounds of Entada phaseoloides. However, genomic or transcriptomic technologies have not been used to study the triterpenoid saponin biosynthetic pathway in this plant.

Results

We performed comparative transcriptome analysis of the root, stem, and leaf tissues of Entada phaseoloides with three independent biological replicates and obtained a total of 53.26 Gb clean data and 116,910 unigenes, with an average N50 length of 1218 bp. Putative functions could be annotated to 42,191 unigenes (36.1%) based on BLASTx searches against the Non-redundant, Uniprot, KEGG, Pfam, GO, KEGG and COG databases. Most of the unigenes related to triterpenoid saponin backbone biosynthesis were specifically upregulated in the stem. A total of 26 cytochrome P450 and 17 uridine diphosphate glycosyltransferase candidate genes related to triterpenoid saponin biosynthesis were identified. The differential expressions of selected genes were further verified by qPT-PCR.

Conclusions

The dataset reported here will facilitate the research about the functional genomics of triterpenoid saponin biosynthesis and genetic engineering of Entada phaseoloides.

Integrated metabolomics and transcriptome analysis on flavonoid biosynthesis in safflower (Carthamus tinctorius L.) under MeJA treatment

BACKGROUND

Safflower (Carthamus tinctorius L.) is an important cash crop, of which the dried tube flower is not only an important raw material for dyes and cosmetics but also an important herb widely used in traditional Chinese medicine. The pigment and bioactive compounds are composed of flavonoids (mainly quinone chalcones), and studies have reported that MeJA can promote the biosynthesis of quinone chalcones, but the mechanism underlying the effect of MeJA in safflower remains unclear. Here, we attempt to use metabolomics and transcriptome technologies to analyse the molecular mechanism of flavonoid biosynthesis under MeJA treatment in safflower.

RESULTS

Based on a UHPLC-ESI-MS/MS detection platform and a self-built database (including hydroxysafflor yellow A, HSYA), a total of 209 flavonoid metabolites were detected, and 35 metabolites were significantly different after treatment with MeJA. Among them, 24 metabolites were upregulated upon MeJA treatment, especially HSYA. Eleven metabolites were downregulated after MeJA treatment. Integrated metabolomics and transcriptome analysis showed that MeJA might upregulate the expression of upstream genes in the flavonoid biosynthesis pathway (such as CHSs, CHIs and HCTs) and downregulate the expression of downstream genes (such as F3Ms, ANRs and ANSs), thus promoting the biosynthesis of quinone chalcones, such as HSYA. The transcription expressions of these genes were validated by real-time PCR. In addition, the promoters of two genes (CtCHI and CtHCT) that were significantly upregulated under MeJA treatment were cloned and analysed. 7 and 3 MeJA response elements were found in the promoters, respectively.

CONCLUSIONS

MeJA might upregulate the expression of the upstream genes in the flavonoid biosynthesis pathway and downregulate the expression of the downstream genes, thus promoting the biosynthesis of quinone chalcones. Our results provide insights and basic data for the molecular mechanism analysis of flavonoid synthesis in safflower under MeJA treatment.

Full-length sequencing of ginkgo transcriptomes for an in-depth understanding of flavonoid and terpenoid trilactone biosynthesis

Ginkgo biloba L. is regarded as the most ancient living tree, and its kernel has been used as a traditional Chinese medicine for more than 2,000 years. The leaf extracts of this tree have been among the bestselling herbal remedies in Western countries since the last century. To understand the biosynthesis of the pharmacologically active ingredients in G. biloba, flavonoids and terpenoid trilactones (TTLs), we sequenced the transcriptomes of G. biloba leaves, kernels and testae with Iso-Seq and RNA-Seq technologies and obtained 152,524 clean consensus reads. When these reads were used to improve the annotation of the G. biloba genome, 4,856 novel genes, 25,583 new isoforms of previously annotated genes and 4,363 lncRNAs were discovered. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses indicated that genes involved in growth, regulation and response to stress were more likely to be regulated by alternative splicing (AS) or alternative polyadenylation (APA), which represent the two most important posttranscriptional regulation mechanisms. It was found that some of the characterized genes involved in the biosynthesis of flavonoids and TTLs were also possibly regulated by AS and APA. Using phylogenetic and gene expression pattern analyses, some candidate genes for the biosynthesis of flavonoids and TTLs were screened. After qRT-PCR validation, the final candidate genes for flavonoid biosynthesis included three UDP-glycosyltransferases and one MYB transcription factor, while the candidate genes for TTL biosynthesis included two cytochrome P450 and one WRKY transcription factor. Our study suggested that Iso-Seq may play an important role in improving genome annotation, elucidating AS and APA mechanisms and discovering candidate genes involved in the biosynthesis of some secondary metabolites.

Integrated metabolomic and transcriptomic profiling reveals the tissue-specific flavonoid compositions and their biosynthesis pathways in Ziziphora bungeana

Background

Ziziphora bungeana Juz. is a folk medicine from the Xinjiang Uygur Autonomous Region. The herb or the aerial parts of it have been used to medicinally treat cardiovascular diseases. Flavonoids are the main pharmacologically active ingredients in Z. bungeana. Identification of the tissue-specific distribution of flavonoids in Z. bungeana is crucial for effective and sustainable medicinal use of the plant. Furthermore, understanding of the biosynthesis pathways of these flavonoids in Z. bungeana is of great biological significance.

Methods

The flavonoids from different tissues of Z. bungeana were identified using liquid chromatography-tandem mass spectrometry (LC–MS/MS). The full-length transcriptome of Z. bungeana was determined using a strategy based on a combination of Illumina and PacBio sequencing techniques. The functions of differentially expressed unigenes were predicted using bioinformatics methods and further investigated by real-time quantitative PCR and phylogenetic relationship analysis.

Results

Among the 12 major flavonoid components identified from Z. bungeana extracts, linarin was the most abundant component. Nine flavonoids were identified as characteristic components of specific tissues. Transcriptome profiling and bioinformatic analysis revealed that 18 genes were putatively involved in flavonoid biosynthesis. The gene expression and phylogenetic analysis results indicated that ZbPALs, Zb4CL3, ZbCHS1, and ZbCHI1 may be involved in the biosynthesis of the main flavonoid intermediate. ZbFNSII, ZbANS, and ZbFLS may be involved in the biosynthesis of flavones, anthocyanins, and flavonols, respectively. A map of the biosynthesis pathways of the 12 major flavonoids in Z. bungeana is proposed.

Conclusions

The chemical constituent analysis revealed the compositions of 9 characteristic flavonoids in different tissues of Z. bungeana. Linarin can be hydrolysed into acacetin to exert a pharmaceutical role. Apigenin-7-O-rutinoside is hypothesised to be the precursor of linarin in Z. bungeana. There was greater content of linarin in the aerial parts of the plant than in the whole herb, which provides a theoretical basis for using the aerial parts of Z. bungeana for medicine. These results provide a valuable reference for further research on the flavonoid biosynthesis pathways of Z. bungeana and will be significant for the effective utilisation and ecological protection of Z. bungeana.

A transcriptome atlas of silkworm silk glands revealed by PacBio single-molecule long-read sequencing

The silk gland of the silkworm Bombyx mori is a specialized organ where silk proteins are efficiently synthesized under precise regulation that largely determines the properties of silk fibers. To understand the genes involved in the regulation of silk protein synthesis, considerable research has focused on the transcripts expressed in silk glands; however, the complete transcriptome profile of this organ has yet to be elucidated. Here, we report a full-length silk gland transcriptome obtained by PacBio single-molecule long-read sequencing technology. In total, 11,697 non-redundant transcripts were identified in mixed samples of silk glands dissected from larvae at five developmental stages. When compared with the published reference, the full-length transcripts optimized the structures of 3002 known genes, and a total of 9061 novel transcripts with an average length of 2171 bp were detected. Among these, 1403 (15.5%) novel transcripts were computationally revealed to be lncRNAs, 8135 (89.8%) novel transcripts were annotated to different protein and nucleotide databases, and 5655 (62.4%) novel transcripts were predicted to have complete ORFs. Furthermore, we found 1867 alternative splicing events, 2529 alternative polyadenylation events, 784 fusion events and 6596 SSRs. This study provides a comprehensive set of reference transcripts and greatly revises and expands the available silkworm transcript data. In addition, these data will be very useful for studying the regulatory mechanisms of silk protein synthesis.

The pomegranate (Punica granatum L.) draft genome dissects genetic divergence between soft- and hard-seeded cultivars

Complete and highly accurate reference genomes and gene annotations are indispensable for basic biological research and trait improvement of woody tree species. In this study, we integrated single-molecule sequencing and high-throughput chromosome conformation capture techniques to produce a high-quality and long-range contiguity chromosome-scale genome assembly of the soft-seeded pomegranate cultivar 'Tunisia'. The genome covers 320.31 Mb (scaffold N50 = 39.96 Mb; contig N50 = 4.49 Mb) and includes 33 594 protein-coding genes. We also resequenced 26 pomegranate varieties that varied regarding seed hardness. Comparative genomic analyses revealed many genetic differences between soft- and hard-seeded pomegranate varieties. A set of selective loci containing SUC8-like, SUC6, FoxO and MAPK were identified by the selective sweep analysis between hard- and soft-seeded populations. An exceptionally large selective region (26.2 Mb) was identified on chromosome 1. Our assembled pomegranate genome is more complete than other currently available genome assemblies. Our results indicate that genomic variations and selective genes may have contributed to the genetic divergence between soft- and hard-seeded pomegranate varieties.

Transcriptomic and metabolomic analyses of Lycium ruthenicum and Lycium barbarum fruits during ripening

Red wolfberry (or goji berry, Lycium barbarum; LB) is an important agricultural product with a high content of pharmacologically important secondary metabolites such as phenylpropanoids. A close relative, black wolfberry (L. ruthenicum; LR), endemic to the salinized deserts of northwestern China, is used only locally. The two fruits exhibit many morphological and phytochemical differences, but genetic mechanisms underlying them remain poorly explored. In order to identify the genes of interest for further studies, we studied transcriptomic (Illumina HiSeq) and metabolomic (LC-MS) profiles of the two fruits during five developmental stages (young to ripe). As expected, we identified much higher numbers of significantly differentially regulated genes (DEGs) than metabolites. The highest numbers were identified in pairwise comparisons including the first stage for both species, but total numbers were consistently somewhat lower for the LR. The number of differentially regulated metabolites in pairwise comparisons of developmental stages varied from 66 (stages 3 vs 4) to 133 (stages 2 vs 5) in both species. We identified a number of genes (e.g. AAT1, metE, pip) and metabolites (e.g. rutin, raffinose, galactinol, trehalose, citrulline and DL-arginine) that may be of interest to future functional studies of stress adaptation in plants. As LB is also highly suitable for combating soil desertification and alleviating soil salinity/alkalinity/pollution, its potential for human use may be much wider than its current, highly localized, relevance.

Full-length transcriptome analysis of Coptis deltoidea and identification of putative genes involved in benzylisoquinoline alkaloids biosynthesis based on combined sequencing platforms

The study carry out comprehensive transcriptome analysis of C. deltoidea and exploration of BIAs biosynthesis and accumulation based on UHPLC-MS/MS and combined sequencing platforms. Coptis deltoidea is an important medicinal plant with a long history of medicinal use, which is rich in benzylisoquinoline alkaloids (BIAs). In this study, Ultra performance liquid chromatography-electrospray ionization tandem mass spectrometry (UHPLC-ESI-MS/MS) and combined sequencing platforms were performed for exploration of BIAs biosynthesis, accumulation and comprehensive transcriptome analysis of C. deltoidea. By metabolism profiling, the accumulation of ten BIAs was analyzed using UHPLC-MS/MS and different contents were observed in different organs. From transcriptome sequencing result, we applied single-molecule real-time (SMRT) sequencing to C. deltoidea and generated a total of 75,438 full-length transcripts. We proposed the candidate biosynthetic pathway of tyrosine, precursor of BIAs, and identified 64 full length-transcripts encoding enzymes putatively involved in BIAs biosynthesis. RNA-Seq data indicated that the majority of genes exhibited relatively high expression level in roots. Transport of BIAs was also important for their accumulation. Here, 9 ABC transporters and 2 MATE transporters highly homologous to known alkaloid transporters related with BIAs transport in roots and rhizomes were identified. These findings based on the combined sequencing platforms provide valuable genetic information for C. deltoidea and the results of transcriptome combined with metabolome analysis can help us better understand BIAs biosynthesis and transport in this medicinal plant. The information will be critical for further characterization of C. deltoidea transcriptome and molecular-assisted breeding for this medicinal plant with scarce resources.

Transcriptome analysis of flavonoid biosynthesis in safflower flowers grown under different light intensities

Background

Safflower (Carthamus tinctorius L.) is a domesticated species with a long history of cultivation and widespread distribution across the globe, and light plays an important role in controlling its distribution boundary. Flowers from safflower have been widely used in traditional Chinese medicine because of their ability to improve cerebral blood flow. Flavonoids are the main active compounds in safflower and have many pharmacological effects. In this study, we aimed to explore the relationship between different light intensities and flavonoid biosynthesis in safflower flowers cultivated in greenhouse.

Methods

The transcriptome of safflower flowers grown under different light intensities were sequenced through BGISEQ-500 platform. After assembled and filtered, Unigenes were annotated by aligning with seven functional databases. Differential expression analysis of two samples was performed with the DEseq2 package. Differentially expressed genes (DEGs) related with flavonoids biosynthesis were analyzed by Real-time PCR (RT-PCR). Flavonoids accumulation in flowers were determined by high performance liquid chromatography and spectrophotometer.

Results

Transcriptome analysis of safflower flowers cultivated under different light intensities was performed. A total of 99.16 Gb data were obtained, and 78,179 Unigenes were annotated. Among the DEGs, 13 genes were related to flavonoid biosynthesis. The differential expressions of seven key genes were confirmed by RT-PCR. In addition, the levels of some flavonoids were measured in safflower flowers grown under different light intensities. CtHCT3 gene expression showed a significantly negative correlation with kaempferol content in safflower grown under different light intensities.

Conclusion

Our results strongly suggested that the reduction in light intensity in a suitable range promoted flavonoid biosynthesis in safflower flowers. We suggest that the expressions of HCT genes played an important role in flavonoid accumulation in safflower flowers. Our study lays a foundation for further research on the effects of light on flavonoid biosynthesis in safflower.

Transcriptome and Phytochemical Analyses Provide New Insights Into Long Non-Coding RNAs Modulating Characteristic Secondary Metabolites of Oolong Tea (Camellia sinensis) in Solar-Withering

Oolong tea is a popular and semi-fermented beverage. During the processing of tea leaves, withering is the first indispensable process for improving flavor. However, the roles of long non-coding RNAs (lncRNAs) and the characteristic secondary metabolites during the withering of oolong tea leaves remain unknown. In this study, phytochemical analyses indicated that total polyphenols, flavonoids, catechins, epigallocatechin (EGC), catechin gallate (CG), gallocatechin gallate (GCG), epicatechin gallate (ECG), and epigallocatechin gallate (EGCG) were all less abundant in the solar-withered leaves (SW) than in the fresh leaves (FL) and indoor-withered leaves (IW). In contrast, terpenoid, jasmonic acid (JA), and methyl jasmonate (MeJA) contents were higher in the SW than in the FL and IW. By analyzing the transcriptome data, we detected 32,036 lncRNAs. On the basis of the Kyoto Encyclopedia of Genes and Genomes analysis, the flavonoid metabolic pathway, the terpenoid metabolic pathway, and the JA/MeJA biosynthesis and signal transduction pathway were enriched pathways. Additionally, 63 differentially expressed lncRNAs (DE-lncRNAs) and 23 target genes were identified related to the three pathways. A comparison of the expression profiles of the DE-lncRNAs and their target genes between the SW and IW revealed four up-regulated genes (FLS, CCR, CAD, and HCT), seven up-regulated lncRNAs, four down-regulated genes (4CL, CHI, F3H, and F3'H), and three down-regulated lncRNAs related to flavonoid metabolism; nine up-regulated genes (DXS, CMK, HDS, HDR, AACT, MVK, PMK, GGPPS, and TPS), three up-regulated lncRNAs, and six down-regulated lncRNAs related to terpenoid metabolism; as well as six up-regulated genes (LOX, AOS, AOC, OPR, ACX, and MFP2), four up-regulated lncRNAs, and three down-regulated lncRNAs related to JA/MeJA biosynthesis and signal transduction. These results suggested that the expression of DE-lncRNAs and their targets involved in the three pathways may be related to the low abundance of the total polyphenols, flavonoids, and catechins (EGC, CG, GCG, ECG, and EGCG) and the high abundance of terpenoids in the SW. Moreover, solar irradiation, high JA and MeJA contents, and the endogenous target mimic (eTM)-related regulatory mechanism in the SW were also crucial for increasing the terpenoid levels. These findings provide new insights into the greater contribution of solar-withering to the high-quality flavor of oolong tea compared with the effects of indoor-withering.

Full-length transcriptome sequences of Agropyron cristatum facilitate the prediction of putative genes for thousand-grain weight in a wheat-A. cristatum translocation line

Background

Agropyron cristatum (L.) Gaertn. (2n = 4x = 28; genomes PPPP) is a wild relative of common wheat (Triticum aestivum L.) and provides many desirable genetic resources for wheat improvement. However, there is still a lack of reference genome and transcriptome information for A. cristatum, which severely impedes functional and molecular breeding studies.

Results

Single-molecule long-read sequencing technology from Pacific Biosciences (PacBio) was used to sequence full-length cDNA from a mixture of leaves, roots, stems and caryopses and constructed the first full-length transcriptome dataset of A. cristatum, which comprised 44,372 transcripts. As expected, the PacBio transcripts were generally longer and more complete than the transcripts assembled via the Illumina sequencing platform in previous studies. By analyzing RNA-Seq data, we identified tissue-enriched transcripts and assessed their GO term enrichment; the results indicated that tissue-enriched transcripts were enriched for particular molecular functions that varied by tissue. We identified 3398 novel and 1352 A. cristatum-specific transcripts compared with the wheat gene model set. To better apply this A. cristatum transcriptome, the A. cristatum transcripts were integrated with the wheat genome as a reference sequence to try to identify candidate A. cristatum transcripts associated with thousand-grain weight in a wheat-A. cristatum translocation line, Pubing 3035.

Conclusions

Full-length transcriptome sequences were used in our study. The present study not only provides comprehensive transcriptomic insights and information for A. cristatum but also proposes a new method for exploring the functional genes of wheat relatives under a wheat genetic background. The sequence data have been deposited in the NCBI under BioProject accession number PRJNA534411.

Full-length transcriptome analysis of shade-induced promotion of tuber production in Pinellia ternata

BACKGROUND

Pinellia ternata is native to China and has been used as a traditional herb due to its antiemetic, antitussive, analgesic, and anxiolytic effects. When exposed to strong light intensity and high temperature during the reproductive growth process, P. ternata withers in a phenomenon known as "sprout tumble", which largely limits tuber production. Shade was previously found to delay sprout tumble formation (STF); however, no information exists regarding this process at the molecular level. Hence, we determined the genes involved in tuber development and STF in P. ternata.

RESULTS

Compared to that with natural sun-light (control), shade significantly induced chlorophyll accumulation, increased chlorophyll fluorescence parameters including initial fluorescence, maximal fluorescence, and qP, and dramatically repressed chlorophyll a:b and NPQ. Catalase (CAT) activity was largely induced by shade, and tuber products were largely increased in this environment. Transcriptome profiles of P. ternata grown in natural sun-light and shaded environments were analyzed by a combination of next generation sequencing (NGS) and third generation single-molecule real-time (SMRT) sequencing. Corrections of SMRT long reads based on NGS short reads yielded 136,163 non-redundant transcripts, with an average N50 length of 2578 bp. In total, 6738 deferentially-expressed genes (DEGs) were obtained from the comparisons, specifically D5S vs D5CK, D20S vs D20CK, D20S vs D5S, and D20CK vs D5CK, of which, 6384 DEGs (94.8%) were generated from the D20S vs D20CK comparison. Gene annotation and functional analyses revealed that these genes were related to auxin signal transduction, polysaccharide and sugar metabolism, phenylpropanoid biosynthesis, and photosynthesis. Moreover, the expression of genes enriched in photosynthesis appeared to be significantly altered by shade. The expression patterns of 16 candidate genes were consistent with changes in their transcript abundance as identified by RNA-Seq, and these might contribute to STF and tuber production.

CONCLUSION

The full-length transcripts identified in this study have provided a more accurate depiction of P. ternata gene transcription. Further, we identified potential genes involved in STF and tuber growth. Such data could serve as a genetic resource and a foundation for further research on this important traditional herb.

Single-molecule real-time transcript sequencing identified flowering regulatory genes in Crocus sativus

Background

Saffron crocus (Crocus sativus) is a valuable spice with medicinal uses in gynaecopathia and nervous system diseases. Identify flowering regulatory genes plays a vital role in increasing flower numbers, thereby resulting in high saffron yield.

Results

Two full length transcriptome gene sets of flowering and non-flowering saffron crocus were established separately using the single-molecule real-time (SMRT) sequencing method. A total of sixteen SMRT cells generated 22.85 GB data and 75,351 full-length saffron crocus unigenes on the PacBio RS II panel and further obtained 79,028 SSRs, 72,603 lncRNAs and 25,400 alternative splicing (AS) events. Using an Illumina RNA-seq platform, an additional fifteen corms with different flower numbers were sequenced. Many differential expression unigenes (DEGs) were screened separately between flowering and matched non-flowering top buds with cold treatment (1677), flowering top buds of 20 g corms and non-flowering top buds of 6 g corms (1086), and flowering and matched non-flowering lateral buds (267). A total of 62 putative flower-related genes that played important roles in vernalization (VRNs), gibberellins (G3OX, G2OX), photoperiod (PHYB, TEM1, PIF4), autonomous (FCA) and age (SPLs) pathways were identified and a schematic representation of the flowering gene regulatory network in saffron crocus was reported for the first time. After validation by real-time qPCR in 30 samples, two novel genes, PB.20221.2 (p = 0.004, r = 0.52) and PB.38952.1 (p = 0.023, r = 0.41), showed significantly higher expression levels in flowering plants. Tissue distribution showed specifically high expression in flower organs and time course expression analysis suggested that the transcripts increasingly accumulated during the flower development period.

Conclusions

Full-length transcriptomes of flowering and non-flowering saffron crocus were obtained using a combined NGS short-read and SMRT long-read sequencing approach. This report is the first to describe the flowering gene regulatory network of saffron crocus and establishes a reference full-length transcriptome for future studies on saffron crocus and other Iridaceae plants.

Overexpression of a Novel Cytochrome P450 Promotes Flavonoid Biosynthesis and Osmotic Stress Tolerance in Transgenic Arabidopsis

Flavonoids are mainly associated with growth, development, and responses to diverse abiotic stresses in plants. A growing amount of data have demonstrated the biosynthesis of flavonoids through multienzyme complexes of which the membrane-bounded cytochrome P450 supergene family shares a crucial part. However, the explicit regulation mechanism of Cytochrome P450s related to flavonoid biosynthesis largely remains elusive. In the present study, we reported the identification of a stress-tolerant flavonoid biosynthetic CtCYP82G24 gene from Carthamus tinctorius. The transient transformation of CtCYP82G24 determined the subcellular localization to the cytosol. Heterologously expressed CtCYP82G24 was effective to catalyze the substrate-specific conversion, promoting the de novo biosynthesis of flavonoids in vitro. Furthermore, a qRT-PCR assay and the accumulation of metabolites demonstrated that the expression of CtCYP82G24 was effectively induced by Polyethylene glycol stress in transgenic Arabidopsis. In addition, the overexpression of CtCYP82G24 could also trigger expression levels of several other flavonoid biosynthetic genes in transgenic plants. Taken together, our findings suggest that CtCYP82G24 overexpression plays a decisive regulatory role in PEG-induced osmotic stress tolerance and alleviates flavonoid accumulation in transgenic Arabidopsis.

Integrating molecular characterization and metabolites profile revealed CtCHI1's significant role in Carthamus tinctorius L

BACKGROUND

As a traditional Chinese herb, safflower (Carthamus tinctorius L.) is valued for its florets to prevent cardiovascular and cerebrovascular diseases. Basing on previous chemical analysis, the main active compounds are flavonoids in its florets. Although flavonoid biosynthetic pathway has been well-documented in many model species, unique biosynthetic pathway remains to be explored in safflower. Of note, as an important class of transitional enzymes, chalcone isomerase (CHI) has not been characterized in safflower.

RESULTS

According to our previous research, CHIs were identified in a safflower transcriptome library built by our lab. To characterize CHI in safflower, a CHI gene named CtCHI1 was identified. A multiple sequences alignment and phylogenetic tree demonstrate that CtCHI1 shares 92% amino acid identity and close relationship with CHI to Saussurea medusa. Additionally, subcellular localization analysis indicated CtCHI1-GFP fusion protein was mainly in the cell nucleus. Further, we purified CtCHI1 protein from E. coli which can effectively catalyze isomerization of 2',4',4,6'-tetrahydroxychalcone into naringenin in vitro. Via genetic engineer technology, we successfully obtained transgenic tobacco and safflower lines. In transgenic tobacco, overexpression of CtCHI1 significantly inhibited main secondary metabolites accumulation, including quercetin (~ 79.63% for ovx-5 line) and anthocyanins (~ 64.55% for ovx-15 line). As shown in transgenic safflower, overexpression of CtCHI1 resulted in upstream genes CtPAL3 and CtC4H1 increasing dramatically (up to ~ 3.9fold) while Ct4CL3, CtF3H and CtDFR2 were inhibited. Also, comparing the whole metabolomics database by PCA and PLS-DA between transgenic and control group, 788 potential differential metabolites were marked and most of them displayed up-regulated trends. In parallel, some isolated secondary metabolites, such as hydroxysafflor yellow A (HSYA), rutin, kaempferol-3-O-β-rutinoside and dihydrokaempferol, accumulated in transgenic safflower plants.

CONCLUSIONS

In this study, we found that CtCHI1 is an active, functional, catalytic protein. Moreover, CtCHI1 can negatively and competitively regulate anthocyanins and quercetin pathway branches in tobacco. By contrast, CtCHI1 can positively regulate flavonol and chalcone metabolic flow in safflower. This research provides some clues to understand CHI's differential biochemical functional characterization involving in flavonoid pathway. More molecular mechanisms of CHI remain to be explored in the near future.

Transcriptome Profiling Provides Insight into the Genes in Carotenoid Biosynthesis during the Mesocarp and Seed Developmental Stages of Avocado (Persea americana)

Avocado (Persea americana Mill.) is an economically important crop because of its high nutritional value. However, the absence of a sequenced avocado reference genome has hindered investigations of secondary metabolism. For next-generation high-throughput transcriptome sequencing, we obtained 365,615,152 and 348,623,402 clean reads as well as 109.13 and 104.10 Gb of sequencing data for avocado mesocarp and seed, respectively, during five developmental stages. High-quality reads were assembled into 100,837 unigenes with an average length of 847.40 bp (N50 = 1725 bp). Additionally, 16,903 differentially expressed genes (DEGs) were detected, 17 of which were related to carotenoid biosynthesis. The expression levels of most of these 17 DEGs were higher in the mesocarp than in the seed during five developmental stages. In this study, the avocado mesocarp and seed transcriptome were also sequenced using single-molecule long-read sequencing to acquired 25.79 and 17.67 Gb clean data, respectively. We identified 233,014 and 238,219 consensus isoforms in avocado mesocarp and seed, respectively. Furthermore, 104 and 59 isoforms were found to correspond to the putative 11 carotenoid biosynthetic-related genes in the avocado mesocarp and seed, respectively. The isoform numbers of 10 out of the putative 11 genes involved in the carotenoid biosynthetic pathway were higher in the mesocarp than those in the seed. Besides, alpha- and beta-carotene contents in the avocado mesocarp and seed during five developmental stages were also measured, and they were higher in the mesocarp than in the seed, which validated the results of transcriptome profiling. Gene expression changes and the associated variations in gene dosage could influence carotenoid biosynthesis. These results will help to further elucidate carotenoid biosynthesis in avocado.

Identification of putative genes for polyphenol biosynthesis in olive fruits and leaves using full-length transcriptome sequencing

Olive (Olea europaea) is a rich source of valuable bioactive polyphenols, which has attracted widespread interest. In this study, we combined targeted metabolome, Pacbio ISOseq transcriptome, and Illumina RNA-seq transcriptome to investigate the association between polyphenols and gene expression in the developing olive fruits and leaves. A total of 12 main polyphenols were measured, and 122 transcripts of 17 gene families, 101 transcripts of 9 gene families, and 106 transcripts of 6 gene families that encode for enzymes involved in flavonoid, oleuropein, and hydroxytyrosol biosynthesis were separately identified. Additionally, 232 alternative splicing events of 18 genes related to polyphenol synthesis were analyzed. This is the first time that the third generations of full-length transcriptome technology were used to study the gene expression pattern of olive fruits and leaves. The results of transcriptome combined with targeted metabolome can help us better understand the polyphenol biosynthesis pathways in the olive.

Full-length transcriptome sequencing and methyl jasmonate-induced expression profile analysis of genes related to patchoulol biosynthesis and regulation in Pogostemon cablin

BACKGROUND

Pogostemon cablin (Blanco) Benth. (Patchouli) is an important aromatic and medicinal plant and widely used in traditional Chinese medicine as well as in the perfume industry. Patchoulol is the primary bioactive component in P. cablin, its biosynthesis has attracted widespread interests. Previous studies have surveyed the putative genes involved in patchoulol biosynthesis using next-generation sequencing method; however, technical limitations generated by short-read sequencing restrict the yield of full-length genes. Additionally, little is known about the expression pattern of genes especially patchoulol biosynthesis related genes in response to methyl jasmonate (MeJA). Our understanding of patchoulol biosynthetic pathway still remained largely incomplete to date.

RESULTS

In this study, we analyzed the morphological character and volatile chemical compounds of P. cablin cv. 'Zhanxiang', and 39 volatile chemical components were detected in the patchouli leaf using GC-MS, most of which were sesquiterpenes. Furthermore, high-quality RNA isolated from leaves and stems of P. cablin were used to generate the first full-length transcriptome of P. cablin using PacBio isoform sequencing (Iso-Seq). In total, 9.7 Gb clean data and 82,335 full-length UniTransModels were captured. 102 transcripts were annotated as 16 encoding enzymes involved in patchouli alcohol biosynthesis. Accorded with the uptrend of patchoulol content, the vast majority of genes related to the patchoulol biosynthesis were up-regulated after MeJA treatment, indicating that MeJA led to an increasing synthesis of patchoulol through activating the expression level of genes involved in biosynthesis pathway of patchoulol. Moreover, expression pattern analysis also revealed that transcription factors participated in JA regulation of patchoulol biosynthesis were differentially expressed.

CONCLUSIONS

The current study comprehensively reported the morphological specificity, volatile chemical compositions and transcriptome characterization of the Chinese-cultivated P. cablin cv. 'Zhanxiang', these results contribute to our better understanding of the physiological and molecular features of patchouli, especially the molecular mechanism of biosynthesis of patchoulol. Our full-length transcriptome data also provides a valuable genetic resource for further studies in patchouli.

Differential gene expression for carotenoid biosynthesis in a green alga Ulva prolifera based on transcriptome analysis

BACKGROUND

Carotenoids are widely distributed in plants and algae, and their biosynthesis has attracted widespread interest. Carotenoid-related research has mostly focused on model species, and there is a lack of data on the carotenoid biosynthetic pathway in U. prolifera that is the main species leading to green tide, a harmful plague of floating green algae.

RESULTS

The carotenoid content of U. prolifera samples, that is the main species leading to green tide, a harmful plague of floating green algae at different temperatures revealed that its terpenoid was highest in the samples subjected to high temperature at 28 °C (H), followed by the samples subjected to low temperature at 12 °C (L). Its terpenoid was lowest in the samples subjected to medium temperature at 20 °C (M). We conducted transcriptome sequencing (148.5 million raw reads and 49,676 unigenes in total) of samples that were subjected to different temperatures to study the carotenoid biosynthesis of U. prolifera. There were 1125-3164 significant differentially expressed genes between L, M and H incubation temperatures, of which 11-672 genes were upregulated and 453-3102 genes were downregulated. A total of 3164 genes were significantly differentially expressed between H and M, of which 62 genes were upregulated and 3102 genes were downregulated. A total of 2669 significant differentially expressed genes were observed between L and H, of which 11 genes were upregulated and 2658 genes were downregulated. A total of 13 genes were identified to be involved in carotenoid biosynthesis in U. prolifera, and the expression levels of the majority were highest at H and lowest at M of incubation temperature. Both the carotenoid concentrations and the expression of the analysed genes were lowest in the normal temperature group, while low temperature and high temperature seemed to activate the biosynthesis of carotenoids in U. prolifera.

CONCLUSIONS

In this study, transcriptome sequencing provided critical information for understanding the accumulation of carotenoids and will serve as an important reference for the study of other metabolic pathways in U. prolifera.

Flavonoid biosynthetic pathways in plants: Versatile targets for metabolic engineering

Plants, fungi, and microorganisms are equipped with biosynthesis machinery for producing thousands of secondary metabolites. These compounds have important functions in nature as a defence against predators or competitors as well as other ecological significances. The full utilization of these compounds for food, medicine, and other purposes requires a thorough understanding of their structures and the distinct biochemical pathways of their production in cellular systems. In this review, flavonoids as classical examples of secondary metabolites are employed to highlight recent advances in understanding how valuable compounds can be regulated at various levels. With extensive diversity in their chemistry and pharmacology, understanding the metabolic engineering of flavonoids now allows us to fine-tune the eliciting of their production, accumulation, and extraction from living systems. More specifically, recent advances in the shikimic acid and acetate biosynthetic pathways of flavonoids production from metabolic engineering point of view, from genes expression to multiple principles of regulation, are addressed. Specific examples of plants and microorganisms as the sources of flavonoids-based compounds with particular emphasis on therapeutic applications are also discussed.