Comprehensive identification of the full-length transcripts and alternative splicing related to the secondary metabolism pathways in the tea plant (Camellia sinensis)

Alternative Splicing Underpins the ALMT9 Transporter Function for Vacuolar Malic Acid Accumulation in Apple

Vacuolar malic acid accumulation largely determines fruit acidity, a key trait for the taste and flavor of apple and other fleshy fruits. Aluminum-activated malate transporter 9 (ALMT9/Ma1) underlies a major genetic locus, Ma, for fruit acidity in apple, but how the protein transports malate across the tonoplast is unclear. Here, it is shown that overexpression of the coding sequence of Ma1 (Ma1α) drastically decreases fruit acidity in "Royal Gala" apple, leading to uncovering alternative splicing underpins Ma1's function. Alternative splicing generates two isoforms: Ma1β is 68 amino acids shorter with much lower expression than the full-length protein Ma1α. Ma1β does not transport malate itself but interacts with the functional Ma1α to form heterodimers, creating synergy with Ma1α for malate transport in a threshold manner (When Ma1β/Ma1α ≥ 1/8). Overexpression of Ma1α triggers feedback inhibition on the native Ma1 expression via transcription factor MYB73, decreasing the Ma1β level well below the threshold that leads to significant reductions in Ma1 function and malic acid accumulation in fruit. Overexpression of Ma1α and Ma1β or genomic Ma1 increases both isoforms proportionally and enhances fruit malic acid accumulation. These findings reveal an essential role of alternative splicing in ALMT9-mediated malate transport underlying apple fruit acidity.

Molecular Mechanism of Resistance to Alternaria alternata Apple Pathotype in Apple by Alternative Splicing of Transcription Factor MdMYB6-like

As a fruit tree with great economic value, apple is widely cultivated in China. However, apple leaf spot disease causes significant damage to apple quality and economic value. In our study, we found that MdMYB6-like is a transcription factor without auto-activation activity and with three alternative spliced variants. Among them, MdMYB6-like-β responded positively to the pathogen infection. Overexpression of MdMYB6-like-β increased the lignin content of leaves and improved the pathogenic resistance of apple flesh callus. In addition, all three alternative spliced variants of MdMYB6-like could bind to the promoter of MdBGLU H. Therefore, we believe that MdMYB6-like plays an important role in the infection process of the pathogen and lays a solid foundation for breeding disease-resistant cultivars of apple in the future.

Tea-Derived Polyphenols Enhance Drought Resistance of Tea Plants (Camellia sinensis) by Alleviating Jasmonate-Isoleucine Pathway and Flavonoid Metabolism Flow

Extreme drought weather has occurred frequently in recent years, resulting in serious yield loss in tea plantations. The study of drought in tea plantations is becoming more and more intensive, but there are fewer studies on drought-resistant measures applied in actual production. Therefore, in this study, we investigated the effect of exogenous tea polyphenols on the drought resistance of tea plant by pouring 100 mg·L-1 of exogenous tea polyphenols into the root under drought. The exogenous tea polyphenols were able to promote the closure of stomata and reduce water loss from leaves under drought stress. Drought-induced malondialdehyde (MDA) accumulation in tea leaves and roots was also significantly reduced by exogenous tea polyphenols. Combined transcriptomic and metabolomic analyses showed that exogenous tea polyphenols regulated the abnormal responses of photosynthetic and energy metabolism in leaves under drought conditions and alleviated sphingolipid metabolism, arginine metabolism, and glutathione metabolism in the root system, which enhanced the drought resistance of tea seedlings. Exogenous tea polyphenols induced jasmonic acid-isoleucine (JA-ILE) accumulation in the root system, and the jasmonic acid-isoleucine synthetase gene (TEA028623), jasmonic acid ZIM structural domain proteins (JAMs) synthesis genes (novel.22237, TEA001821), and the transcription factor MYC2 (TEA014288, TEA005840) were significantly up-regulated. Meanwhile, the flavonoid metabolic flow was significantly altered in the root; for example, the content of EGCG, ECG, and EGC was significantly increased. Thus, exogenous tea polyphenols enhance the drought resistance of tea plants through multiple pathways.

Alternative splicing of CsWRKY21 positively regulates cold response in tea plant

Alternative splicing (AS) was an important post-transcriptional mechanism that involved in plant resistance to adversity stress. WRKY transcription factors function as transcriptional activators or repressors to modulate plant growth, development and stress response. However, the role of alternate splicing of WRKY in cold tolerance is poorly understood in tea plants. In this study, we found that the CsWRKY21 transcription factor, a member of the WRKY IId subfamily, was induced by low temperature. Subcellular localization and transcriptional activity assays showed that CsWRKY21 localized to the nucleus and had no transcriptional activation activity. Y1H and dual-luciferase reporter assays showed that CsWRKY21 suppressed expression of CsABA8H and CsUGT by binding with their promoters. Transient overexpression of CsABA8H and CsUGT reduced abscisic acid (ABA) content in tobacco leaves. Furthermore, we discovered that CsWRKY21 undergoes AS in the 5'UTR region. The AS transcript CsWRKY21-b was induced at low temperature, up to 6 folds compared to the control, while the full-length CsWRKY21-a transcript did not significantly change. Western blot analysis showed that the retention of introns in the 5'UTR region of CsWRKY21-b led to higher CsWRKY21 protein content. These results revealed that alternative splicing of CsWRKY21 involved in cold tolerance of tea plant by regulating the protein expression level and then regulating the content of ABA, and provide insights into molecular mechanisms of low temperature defense mediated by AS in tea plant.

Comparative Transcriptome Profiling of CMS-D2 and CMS-D8 Systems Characterizes Fertility Restoration Genes Network in Upland Cotton

Resolving the genetic basis of fertility restoration for cytoplasmic male sterility (CMS) can improve the efficiency of three-line hybrid breeding. However, the genetic determinants of male fertility restoration in cotton are still largely unknown. This study comprehensively compared the full-length transcripts of CMS-D2 and CMS-D8 systems to identify potential genes linked with fertility restorer genes Rf₁ or Rf₂. Target comparative analysis revealed a higher percentage of differential genes in each restorer line as compared to their corresponding sterile and maintainer lines. An array of genes with specific expression in the restorer line of CMS-D2 had functional annotations related to floral development and pathway enrichments in various secondary metabolites, while specifically expressed genes in the CMS-D8 restorer line showed functional annotations related to anther development and pathway enrichment in the biosynthesis of secondary metabolites. Further analysis identified potentially key genes located in the target region of fertility restorer genes Rf₁ or Rf₂. In particular, Ghir_D05G032450 can be the candidate gene related to restorer gene Rf₁, and Ghir_D05G035690 can be the candidate gene associated with restorer gene Rf₂. Further gene expression validation with qRT-PCR confirmed the accuracy of our results. Our findings provide useful insights into decoding the potential regulatory network that retrieves pollen fertility in cotton and will help to further reveal the differences in the genetic basis of fertility restoration for two CMS systems.

Alternative splicing in ABA signaling during seed germination

Seed germination is an essential step in a plant's life cycle. It is controlled by complex physiological, biochemical, and molecular mechanisms and external factors. Alternative splicing (AS) is a co-transcriptional mechanism that regulates gene expression and produces multiple mRNA variants from a single gene to modulate transcriptome diversity. However, little is known about the effect of AS on the function of generated protein isoforms. The latest reports indicate that alternative splicing (AS), the relevant mechanism controlling gene expression, plays a significant role in abscisic acid (ABA) signaling. In this review, we present the current state of the art about the identified AS regulators and the ABA-related changes in AS during seed germination. We show how they are connected with the ABA signaling and the seed germination process. We also discuss changes in the structure of the generated AS isoforms and their impact on the functionality of the generated proteins. Also, we point out that the advances in sequencing technology allow for a better explanation of the role of AS in gene regulation by more accurate detection of AS events and identification of full-length splicing isoforms.

Full-length transcriptome sequencing and comparative transcriptome analysis of Eriocheir sinensis in response to infection by the microsporidian Hepatospora eriocheir

As a new generation of high-throughput sequencing technology, PacBio Iso-Seq technology (Iso-Seq) provides a better alternative sequencing method for the acquisition of full-length unigenes. In this study, a total of 22.27 gigabyte (Gb) subread bases and 128,614 non-redundant unigenes (mean length: 2,324 bp) were obtained from six main tissues of Eriocheir sinensis including the heart, nerve, intestine, muscle, gills and hepatopancreas. In addition, 74,732 unigenes were mapped to at least one of the following databases: Non-Redundant Protein Sequence Database (NR), Gene Ontology (GO), Kyoto Encyclopaedia of Genes and Genomes (KEGG), KEGG Orthology (KO) and Protein family (Pfam). In addition, 6696 transcription factors (TFs), 28,458 long non-coding RNAs (lncRNAs) and 94,230 mRNA-miRNA pairs were identified. Hepatospora eriocheir is the primary pathogen of E. sinensis and can cause hepatopancreatic necrosis disease (HPND); the intestine is the main target tissue. Here, we attempted to identify the key genes related to H. eriocheir infection in the intestines of E. sinensis. By combining Iso-Seq and Illumina RNA-seq analysis, we identified a total of 12,708 differentially expressed unigenes (DEUs; 6,696 upregulated and 6,012 downregulated) in the crab intestine following infection with H. eriocheir. Based on the biological analysis of these DEUs, several key processes were identified, including energy metabolism-related pathways, cell apoptosis and innate immune-related pathways. Twelve selected genes from these DEUs were subsequently verified by quantitative real-time PCR (qRT-PCR) analysis. Our findings enhance our understanding of the E. sinensis transcriptome and the specific association between E. sinensis and H. eriocheir infection.

Combination of long-read and short-read sequencing provides comprehensive transcriptome and new insight for Chrysanthemum morifolium ray-floret colorization

Chrysanthemum morifolium is one of the most popular ornamental plants globally. Owing to its large and complex genome (around 10 Gb, segmental hexaploid), it has been difficult to obtain comprehensive transcriptome, which will promote to perform new breeding technique, such as genome editing, in C. morifolium. In this study, we used single-molecule real-time (SMRT) sequencing and RNA-seq technologies, combined them with an error-correcting process, and obtained high-coverage ray-floret transcriptome. The SMRT-seq data increased the ratio of long mRNAs containing complete open-reading frames, and the combined dataset provided a more complete transcriptomic data than those produced from either SMRT-seq or RNA-seq-derived transcripts. We finally obtained 'Sei Arabella' transcripts containing 928,645 non-redundant mRNA, which showed 96.6% Benchmarking Universal Single-Copy Orthologs (BUSCO) score. We also validated the reliability of the dataset by analyzing a mapping rate, annotation and transcript expression. Using the dataset, we searched anthocyanin biosynthesis gene orthologs and performed a qRT-PCR experiment to assess the usability of the dataset. The assessment of the dataset and the following analysis indicated that our dataset is reliable and useful for molecular biology. The combination of sequencing methods provided genetic information and a way to analyze the complicated C. morifolium transcriptome.

United States tea: A synopsis of ongoing tea research and solutions to United States tea production issues

Tea is a steeped beverage made from the leaves of Camellia sinensis. Globally, this healthy, caffeine-containing drink is one of the most widely consumed beverages. At least 50 countries produce tea and most of the production information and tea research is derived from international sources. Here, we discuss information related to tea production, genetics, and chemistry as well as production issues that affect or are likely to affect emerging tea production and research in the United States. With this review, we relay current knowledge on tea production, threats to tea production, and solutions to production problems to inform this emerging market in the United States.

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.

A Global Analysis of Alternative Splicing of Dichocarpum Medicinal Plants, Ranunculales

Background: The multiple isoforms are often generated from a single gene via Alternative Splicing (AS) in plants, and the functional diversity of the plant genome is significantly increased. Despite well-studied gene functions, the specific functions of isoforms are little known, therefore, the accurate prediction of isoform functions is exceedingly wanted. Methods: Here we perform the first global analysis of AS of Dichocarpum, a medicinal genus of Ranunculales, by utilizing full-length transcriptome datasets of five Chinese endemic Dichocarpum taxa. Multiple software were used to identify AS events, the gene function was annotated based on seven databases, and the protein-coding sequence of each AS isoform was translated into an amino acid sequence. The self-developed software DIFFUSE was used to predict the functions of AS isoforms. Results: Among 8,485 genes with AS events, the genes with two isoforms were the most (6,038), followed by those with three isoforms and four isoforms. Retained intron (RI, 551) was predominant among 1,037 AS events, and alternative 3' splice sites and alternative 5' splice sites were second. The software DIFFUSE was effective in predicting functions of Dichocarpum isoforms, which have not been unearthed. When compared with the sequence alignment-based database annotations, DIFFUSE performed better in differentiating isoform functions. The DIFFUSE predictions on the terms GO:0003677 (DNA binding) and GO: 0010333 (terpene synthase activity) agreed with the biological features of transcript isoforms. Conclusion: Numerous AS events were for the first time identified from full-length transcriptome datasets of five Dichocarpum taxa, and functions of AS isoforms were successfully predicted by the self-developed software DIFFUSE. The global analysis of Dichocarpum AS events and predicting isoform functions can help understand the metabolic regulations of medicinal taxa and their pharmaceutical explorations.

Alternative Splicing and Its Roles in Plant Metabolism

Plant metabolism, including primary metabolism such as tricarboxylic acid cycle, glycolysis, shikimate and amino acid pathways as well as specialized metabolism such as biosynthesis of phenolics, alkaloids and saponins, contributes to plant survival, growth, development and interactions with the environment. To this end, these metabolic processes are tightly and finely regulated transcriptionally, post-transcriptionally, translationally and post-translationally in response to different growth and developmental stages as well as the constantly changing environment. In this review, we summarize and describe the current knowledge of the regulation of plant metabolism by alternative splicing, a post-transcriptional regulatory mechanism that generates multiple protein isoforms from a single gene by using alternative splice sites during splicing. Numerous genes in plant metabolism have been shown to be alternatively spliced under different developmental stages and stress conditions. In particular, alternative splicing serves as a regulatory mechanism to fine-tune plant metabolism by altering biochemical activities, interaction and subcellular localization of proteins encoded by splice isoforms of various genes.

Regulation of biosynthesis of the main flavor-contributing metabolites in tea plant (Camellia sinensis): A review

In the process of adapting to the environment, tea plants (Camellia sinensis) endow tea with unique flavor and health functions, which should be attributed to secondary metabolites, including catechins, L-theanine, caffeine and terpene volatiles. Since the content of these flavor-contributing metabolites are mainly determined by the growth of tea plant, it is very important to understand their alteration and regulation mechanisms. In the present work, we first summarize the distribution, change characteristics of the main flavor-contributing metabolites in different cultivars, organs and under environmental stresses of tea plant. Subsequently, we discuss the regulating mechanisms involved in the biosynthesis of these metabolites based on the existing evidence. Finally, we propose the remarks and perspectives on the future study relating flavor-contributing metabolites. This review would contribute to the acceleration of research on the characteristic secondary metabolites and the breeding programs in tea plants.

Comparative transcriptome analysis of pollen and anther wall reveals novel insights into the regulatory mechanisms underlying anther wall development and its dehiscence in rice

To further understand the regulatory mechanism for anther dehiscence in rice, we carried out transcriptome analysis for the following two tissues: the anther wall and pollen at the anthesis stage. With the anatomical meta-expression data, in addition to these tissues, the differentially expressed genes (DEGs) between the two tissues were further refined to identify 1,717 pollen-preferred genes and 534 anther wall-preferred genes. A GUS transgenic line and RT-qPCR analysis for anther wall-preferred genes supported the fidelity of our gene candidates for further analysis. The refined DEGs were functionally classified through Gene Ontology (GO) enrichment and MapMan analyses. Through the analysis of cis-acting elements and alternative splicing variants, we also suggest the feature of regulatory sequences in promoter regions for anther wall-preferred expression and provide information of the unique splicing variants in anther wall. Subsequently, it was found that hormone signaling and the resulting transcriptional regulation pathways may play an important role in anther dehiscence and anther wall development. Our results could provide useful insights into future research to broaden the molecular mechanism of anther dehiscence or anther wall development in rice.

Theanine metabolism and transport in tea plants (Camellia sinensis L.): advances and perspectives

Theanine, a tea plant-specific non-proteinogenic amino acid, is the most abundant free amino acid in tea leaves. It is also one of the most important quality components of tea because it endows the "umami" taste, relaxation-promoting, and many other health benefits of tea infusion. Its content in tea leaves is directly correlated with the quality and price of green tea. Theanine biosynthesis primarily occurs in roots and is transported to new shoots in tea plants. Recently, great advances have been made in theanine metabolism and transport in tea plants. Along with the deciphering of the genomic sequences of tea plants, new genes in theanine metabolic pathway were discovered and functionally characterized. Theanine transporters were identified and were characterized on the affinity for: theanine, substrate specificity, spatiotemporal expression, and the role in theanine root-to-shoot transport. The mechanisms underlying the regulation of theanine accumulation by: cultivars, seasons, nutrients, and environmental factors are also being rapidly uncovered. Transcription factors were identified to be critical regulators of theanine biosynthesis. In this review, we summarize the progresses in theanine: biosynthesis, catabolism, and transport processes. We also discuss the future studies on theanine in tea plants, and application of the knowledge to crops to synthesize theanine to improve the health-promoting quality of non-tea crops.

Full-Length Transcriptome Characterization and Comparative Analysis of Chosenia arbutifolia

As a unique tree species in the Salicaceae family, Chosenia arbutifolia is used primarily for construction materials and landscape planting in China. Compared with other Salicaceae species members, the genomic resources of C. arbutifolia are extremely scarce. Thus, in the present study, the full-length transcriptome of C. arbutifolia was sequenced by single-molecular real-time sequencing (SMRT) technology based on the PacBio platform. Then, it was compared against those of other Salicaceae species. We generated 17,397,064 subreads and 95,940 polished reads with an average length of 1812 bp, which were acquired through calibration, clustering, and polishing. In total, 50,073 genes were reconstructed, of which 48,174 open reading frames, 4281 long non-coding RNAs, and 3121 transcription factors were discovered. Functional annotation revealed that 47,717 genes had a hit in at least one of five reference databases. Moreover, a set of 12,332 putative SSR markers were screened among the reconstructed genes. Single-copy and special orthogroups, and divergent and conserved genes, were identified and analyzed to find divergence among C. arbutifolia and the five Salicaceae species. To reveal genes involved in a specific function and pathway, enrichment analyses for GO and KEGG were also performed. In conclusion, the present study empirically confirmed that SMRT sequencing realistically depicted the C. arbutifolia transcriptome and provided a comprehensive reference for functional genomic research on Salicaceae species.

EST-SSR Marker Development and Full-Length Transcriptome Sequence Analysis of Tiger Lily (Lilium lancifolium Thunb)

The fast advancement and deployment of sequencing technologies after the Human Genome Project have greatly increased our knowledge of the eukaryotic genome sequences. However, due to technological concerns, high-quality genomic data has been confined to a few key organisms. Moreover, our understanding of which portions of genomes make up genes and which transcript isoforms synthesize these genes is scarce. Therefore, the current study has been designed to explore the reliability of the tiger lily (Lilium lancifolium Thunb) transcriptome. The PacBio-SMRT was used for attaining the complete transcriptomic profile. We obtained a total of 815,624 CCS (Circular Consensus Sequence) reads with an average length of 1295 bp. The tiger lily transcriptome has been sequenced for the first time using third-generation long-read technology. Furthermore, unigenes (38,707), lncRNAs (6852), and TF members (768) were determined based on the transcriptome data, followed by evaluating SSRs (3319). It has also been revealed that 105 out of 128 primer pairs effectively amplified PCR products. Around 15,608 transcripts were allocated to 25 distinct KOG Clusters, and 10,706 unigenes were grouped into 52 functional categories in the annotated transcripts. Until now, no tiger lily lncRNAs have been discovered. Results of this study may serve as an extensive set of reference transcripts and help us learn more about the transcriptomes of tiger lilies and pave the path for further research.

Transcriptomic Complexity in Strawberry Fruit Development and Maturation Revealed by Nanopore Sequencing

The use of alternative transcription start or termination sites (aTSS or aTTS) as well as alternative splicing (AS) produce diverse transcript isoforms, playing indispensable roles in the plant development and environmental adaptations. Despite the advances in the finding of the genome-wide alternatively spliced genes in strawberry, it remains unexplored how AS responds to the developmental cues and what relevance do these outcomes have to the gene function. In this study, we have systematically investigated the transcriptome complexity using long-read Oxford Nanopore Technologies along the four successive developmental stages. The full-length cDNA sequencing results unraveled thousands of previously unexplored transcript isoforms raised from aTSS, aTTS, and AS. The relative contributions of these three processes to the complexity of strawberry fruit transcripts were compared. The aTSS and aTTS were more abundant than the AS. Differentially expressed transcripts unraveled the key transitional role of the white fruit stage. Isoform switches of transcripts from 757 genes were observed. They were associated with protein-coding potential change and domain gain or loss as the main consequences. Those genes with switched isoforms take part in the key processes of maturation in the late stages. A case study using yeast two hybrid analysis supported the functional divergence of the two isoforms of the B-box protein 22. Our results provided a new comprehensive overview of the dynamic transcriptomic landscape during strawberry fruit development and maturation.

Full-length transcriptome analysis of pecan (Carya illinoinensis) kernels

Pecan is rich in bioactive components such as fatty acids (FAs) and flavonoids and is an important nut type worldwide. Therefore, the molecular mechanisms of phytochemical biosynthesis in pecan are a focus of research. Recently, a draft genome and several transcriptomes have been published. However, the full-length mRNA transcripts remain unclear, and the regulatory mechanisms behind the quality components biosynthesis and accumulation have not been fully investigated. In this study, single-molecule long-read sequencing technology was used to obtain full-length transcripts of pecan kernels. In total, 37,504 isoforms of 16,702 genes were mapped to the reference genome. The numbers of known isoforms, new isoforms, and novel isoforms were 9013 (24.03%), 26,080 (69.54%), and 2411 (6.51%), respectively. Over 80% of the transcripts (30,751, 81.99%) had functional annotations. A total of 15,465 alternative splicing (AS) events and 65,761 alternative polyadenylation events were detected; wherein, the retained intron was the predominant type (5652, 36.55%) of AS. Furthermore, 1894 long noncoding RNAs and 1643 transcription factors were predicted using bioinformatics methods. Finally, the structural genes associated with FA and flavonoid biosynthesis were characterized. A high frequency of AS accuracy (70.31%) was observed in FA synthesis-associated genes. This study provides a full-length transcriptome data set of pecan kernels, which will significantly enhance the understanding of the regulatory basis of phytochemical biosynthesis during pecan kernel maturation.

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.

Formation of volatiles in response to tea green leafhopper (Empoasca onukii Matsuda) herbivory in tea plants: a multi-omics study

Combined transcriptome and metabolome analysis of fresh leaf infestation by tea green leafhoppers (Empoasca (Matsumurasca) onukii Matsuda) suggests roles for alternative pre-mRNA splicing and mRNAs in the regulation of aroma formation in tea plants. Oriental Beauty is a high-grade, oolong tea with a pronounced honey-like aroma and rich ripe fruit flavor that develops primarily as a result of the infestation of the fresh leaves by tea green leafhoppers (Empoasca (Matsumurasca) onukii Matsuda). Here, we used PacBio Iso-Seq and RNA-seq analyses to determine the full-length transcripts and gene expression profiles of fresh tea leaves in response to E. (M.) onukii herbivory. We investigated the relationship between RNA-seq, tea metabolites, and aroma response mechanisms in leaves infested by leafhoppers. We found 3644 differentially expressed genes, of which 2552 were up- and 1092 were down-regulated. A total of 49,913 alternative splicing events were predicted, including 324 differential AS events. Moreover, 3105 differentially expressed transcripts were also identified, of which 2295 were up- and 810 were down-regulated. The characterization of expression patterns of the key gene transcript isoforms involved in the aroma formation pathways identified 130 differentially expressed metabolites, 97 of which were up- and 33 were down-regulated. Two key aroma compounds (phenylacetaldehyde and 4-hydroxybenzaldehyde) were highly correlated with genes of the aroma formation pathways. Our results revealed that pre-mRNA AS plays a crucial role in the metabolic regulation surrounding aroma formation under leafhopper herbivory in tea plants.

Molecular Basis of the Distinct Metabolic Features in Shoot Tips and Roots of Tea Plants (Camellia sinensis): Characterization of MYB Regulator for Root Theanine Synthesis

The physiological and metabolic differences between shoot tips and roots of tea plants are significant, and understanding them is required for improvement of tea quality and plant growth. A high-quality full-length transcriptome sequencing on tea plant roots and shoot tips by PacBio SMRT technology was done to gain a further understanding. Approximately 160699 and 166120 full-length transcripts were recovered in roots and shoots, respectively, including 31232 and 41068 novel isoforms and 16960 and 26029 alternative splicing (AS) isoforms. These supported 21699 full-length reads and 31232 and 41068 novel transcripts from root and shoot, respectively, including 1679 long noncoding RNAs (lncRNAs) and 2772 fusion transcripts, which significantly upgrade the Camellia sinensis genome annotation. The respective 6475 and 6981 transcripts in roots and shoots differ in 3'-untranslated regions. Meanwhile, extensive analyses of novel transcripts, ASs, and lncRNAs also revealed a large number of ASs and lincRNAs closely related to the regulation of characteristic secondary metabolites including catechins, theanine, and caffeine. Finally, a root-specific CsMYB6 was characterized to regulate theanine biosynthesis by genetic and molecular analyses. CsMYB6 directly bound to and activate the promoter of theanine synthetase gene (CsTSI). The study lays a foundation for the further investigation of metabolic genomics and regulation in tea plants.

Transcriptome profiling analysis of tea plant (Camellia sinensis) using Oxford Nanopore long-read RNA-Seq technology

Transcriptome profiles have been widely captured using short-read sequencing technology, but there are still limitations partially due to the read length. Here, we generated long reads using Oxford Nanopore PromethION™ technology and short reads using the Illumina sequencing platform to study the transcriptome of root, stem, and leaf of Camellia sinensis cv. Fudingdabai. We mapped the Nanopore reads to the Shuchazao of C. sinensis genome sequence, and the mapping rates ranged from 82.63% to 90.59% (average 86.44%); this is lower than that of the Illumina reads which was 87.83% to 91.14% (average 90.12%). Gene expression level was quantified using the Nanopore and Illumina data and we observed a good agreement. The same tea leaf flavor synthesis pathways were highlighted using both sequencing technologies when analyzing the differentially expressed genes between leaf and root. Alternative splicing was then analyzed, and the intron-retention was observed as the most common alternative splicing. Moreover Nanopore long reads could correct transcript isoform annotation for differential expression investigation purposes. Nanopore sequencing techniques can provide a novel reference basis for molecular analysis of tea plants.

Kirigami-Inspired Inflatables with Programmable Shapes

Kirigami, the Japanese art of paper cutting, has recently enabled the design of stretchable mechanical metamaterials that can be easily realized by embedding arrays of periodic cuts into an elastic sheet. Here, kirigami principles are exploited to design inflatables that can mimic target shapes upon pressurization. The system comprises a kirigami sheet embedded into an unstructured elastomeric membrane. First, it is shown that the inflated shape can be controlled by tuning the geometric parameters of the kirigami pattern. Then, by applying a simple optimization algorithm, the best parameters that enable the kirigami inflatables to transform into a family of target shapes at a given pressure are identified. Furthermore, thanks to the tessellated nature of the kirigami, it is shown that we can selectively manipulate the parameters of the single units to allow the reproduction of features at different scales and ultimately enable a more accurate mimicking of the target.

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.

Integrative expression network analysis of microRNA and gene isoforms in sacred lotus

Background

Gene expression is complex and regulated by multiple molecular mechanisms, such as miRNA-mediated gene inhibition and alternative-splicing of pre-mRNAs. However, the coordination of interaction between miRNAs with different splicing isoforms, and the change of splicing isoform in response to different cellular environments are largely unexplored in plants. In this study, we analyzed the miRNA and mRNA transcriptome from lotus (Nelumbo nucifera), an economically important flowering plant.

Results

Through RNA-seq analyses on miRNAs and their target genes (isoforms) among six lotus tissues, expression of most miRNAs seem to be negatively correlated with their targets and tend to be tissue-specific. Further, our results showed that preferential interactions between miRNAs and hub gene isoforms in one coexpression module which is highly correlated with leaf. Intriguingly, for many genes, their corresponding isoforms were assigned to different co-expressed modules, and they exhibited more divergent mRNA structures including presence and absence of miRNA binding sites, suggesting functional divergence for many isoforms is escalated by both structural and expression divergence. Further detailed functional enrichment analysis of miRNA targets revealed that miRNAs are involved in the regulation of lotus growth and development by regulating plant hormone-related pathway genes.

Conclusions

Taken together, our comprehensive analyses of miRNA and mRNA transcriptome elucidate the coordination of interaction between miRNAs and different splicing isoforms, and highlight the functional divergence of many transcript isoforms from the same locus in lotus.

Transcriptional regulation of amino acid metabolism in response to nitrogen deficiency and nitrogen forms in tea plant root (Camellia sinensis L.)

Free amino acids, including theanine, glutamine and glutamate, contribute greatly to the pleasant taste and multiple health benefits of tea. Amino acids in tea plants are mainly synthesized in roots and transported to new shoots, which are significantly affected by nitrogen (N) level and forms. However, the regulatory amino acid metabolism genes have not been systemically identified in tea plants. Here, we investigated the dynamic changes of free amino acid contents in response to N deficiency and forms in tea plant roots, and systemically identified the genes associated amino acid contents in individual metabolism pathways. Our results showed that glutamate-derived amino acids are the most dynamic in response to various forms of N and N deficiency. We then performed transcriptomic analyses of roots treated with N deficiency and various forms of N, and differentially expressed amino acid metabolic genes in each pathway were identified. The analyses on expression patterns and transcriptional responses of metabolic genes to N treatments provided novel insights for the molecular basis of high accumulation of theanine in tea plant root. These analyses also identified potential regulatory genes in dynamic amino acid metabolism in tea plant root. Furthermore, our findings indicated that the dynamic expression levels of CsGDH, CsAlaDC, CsAspAT, CsSDH, CsPAL, CsSHMT were highly correlated with changes of amino acid contents in their corresponding pathways. Herein, this study provides comprehensive insights into transcriptional regulation of amino acid metabolism in response to nitrogen deficiency and nitrogen forms in tea plant root.

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.

Single-molecule real-time sequencing facilitates the analysis of transcripts and splice isoforms of anthers in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

BACKGROUND

Anther development has been extensively studied at the transcriptional level, but a systematic analysis of full-length transcripts on a genome-wide scale has not yet been published. Here, the Pacific Biosciences (PacBio) Sequel platform and next-generation sequencing (NGS) technology were combined to generate full-length sequences and completed structures of transcripts in anthers of Chinese cabbage.

RESULTS

Using single-molecule real-time sequencing (SMRT), a total of 1,098,119 circular consensus sequences (CCSs) were generated with a mean length of 2664 bp. More than 75% of the CCSs were considered full-length non-chimeric (FLNC) reads. After error correction, 725,731 high-quality FLNC reads were estimated to carry 51,501 isoforms from 19,503 loci, consisting of 38,992 novel isoforms from known genes and 3691 novel isoforms from novel genes. Of the novel isoforms, we identified 407 long non-coding RNAs (lncRNAs) and 37,549 open reading frames (ORFs). Furthermore, a total of 453,270 alternative splicing (AS) events were identified and the majority of AS models in anther were determined to be approximate exon skipping (XSKIP) events. Of the key genes regulated during anther development, AS events were mainly identified in the genes SERK1, CALS5, NEF1, and CESA1/3. Additionally, we identified 104 fusion transcripts and 5806 genes that had alternative polyadenylation (APA).

CONCLUSIONS

Our work demonstrated the transcriptome diversity and complexity of anther development in Chinese cabbage. The findings provide a basis for further genome annotation and transcriptome research in Chinese cabbage.

Establishing a System for Functional Characterization of Full-Length cDNAs of Camellia sinensis

Tea (Camellia sinensis) is enriched with bioactive secondary metabolites, and is one of the most popular nonalcoholic beverages globally. Two tea reference genomes have been reported; however, the functional analysis of tea genes has lagged, mainly due to tea’s recalcitrance to genetic transformation and the absence of alternative high throughput heterologous expression systems. A full-length cDNA collection with a streamlined cloning system is needed in this economically important woody crop species. RNAs were isolated from nine different vegetative tea tissues, pooled, then used to construct a normalized full-length cDNA library. The titer of unamplified and amplified cDNA library was 6.89 × 10⁶ and 1.8 × 10¹⁰ cfu/mL, respectively; the library recombinant rate was 87.2%. Preliminary characterization demonstrated that this collection can complement existing tea reference genomes and facilitate rare gene discovery. In addition, to streamline tea cDNA cloning and functional analysis, a binary vector (pBIG2113SF) was reengineered, seven tea cDNAs isolated from this library were successfully cloned into this vector, then transformed into Arabidopsis. One FL-cDNA, which encodes a putative P_1B-type ATPase 5 (CsHMA5), was characterized further as a proof of concept. We demonstrated that overexpression of CsHMA5 in Arabidopsis resulted in copper hyposensitivity. Thus, our data demonstrated that this represents an efficient system for rare gene discovery and functional characterization of tea genes. The integration of a tea FL-cDNA collection with efficient cloning and a heterologous expression system would facilitate functional annotation and characterization of tea genes.

Optimized sequencing depth and de novo assembler for deeply reconstructing the transcriptome of the tea plant, an economically important plant species

Background

Tea is the oldest and among the world’s most popular non-alcoholic beverages, which has important economic, health and cultural values. Tea is commonly produced from the leaves of tea plants (Camellia sinensis), which belong to the genus Camellia of family Theaceae. In the last decade, many studies have generated the transcriptomes of tea plants at different developmental stages or under abiotic and/or biotic stresses to investigate the genetic basis of secondary metabolites that determine tea quality. However, these results exhibited large differences, particularly in the total number of reconstructed transcripts and the quality of the assembled transcriptomes. These differences largely result from limited knowledge regarding the optimized sequencing depth and assembler for transcriptome assembly of structurally complex plant species genomes.

Results

We employed different amounts of RNA-sequencing data, ranging from 4 to 84 Gb, to assemble the tea plant transcriptome using five well-known and representative transcript assemblers. Although the total number of assembled transcripts increased with increasing sequencing data, the proportion of unassembled transcripts became saturated as revealed by plant BUSCO datasets. Among the five representative assemblers, the Bridger package shows the best performance in both assembly completeness and accuracy as evaluated by the BUSCO datasets and genome alignment. In addition, we showed that Bridger and BinPacker harbored the shortest runtimes followed by SOAPdenovo and Trans-ABySS.

Conclusions

The present study compares the performance of five representative transcript assemblers and investigates the key factors that affect the assembly quality of the transcriptome of the tea plants. This study will be of significance in helping the tea research community obtain better sequencing and assembly of tea plant transcriptomes under conditions of interest and may thus help to answer major biological questions currently facing the tea industry.

The Impact of cDNA Normalization on Long-Read Sequencing of a Complex Transcriptome

Normalization of cDNA is widely used to improve the coverage of rare transcripts in analysis of transcriptomes employing next-generation sequencing. Recently, long-read technology has been emerging as a powerful tool for sequencing and construction of transcriptomes, especially for complex genomes containing highly similar transcripts and transcript-spliced isoforms. Here, we analyzed the transcriptome of sugarcane, a highly polyploidy plant genome, by PacBio isoform sequencing (Iso-Seq) of two different cDNA library preparations, with and without a normalization step. The results demonstrated that, while the two libraries included many of the same transcripts, many longer transcripts were removed, and many new generally shorter transcripts were detected by normalization. For the same input cDNA and data yield, the normalized library recovered more total transcript isoforms and number of predicted gene families and orthologous groups, resulting in a higher representation for the sugarcane transcriptome, compared to the non-normalized library. The non-normalized library, on the other hand, included a wider transcript length range with more longer transcripts above ∼1.25 kb and more transcript isoforms per gene family and gene ontology terms per transcript. A large proportion of the unique transcripts comprising ∼52% of the normalized library were expressed at a lower level than the unique transcripts from the non-normalized library, across three tissue types tested including leaf, stalk, and root. About 83% of the total 5,348 predicted long noncoding transcripts was derived from the normalized library, of which ∼80% was derived from the lowly expressed fraction. Functional annotation of the unique transcripts suggested that each library enriched different functional transcript fractions. This demonstrated the complementation of the two approaches in obtaining a complete transcriptome of a complex genome at the sequencing depth used in this study.

Expanding Alternative Splicing Identification by Integrating Multiple Sources of Transcription Data in Tomato

Tomato (Solanum lycopersicum) is an important vegetable and fruit crop. Its genome was completely sequenced and there are also a large amount of available expressed sequence tags (ESTs) and short reads generated by RNA sequencing (RNA-seq) technologies. Mapping transcripts including mRNA sequences, ESTs, and RNA-seq reads to the genome allows identifying pre-mRNA alternative splicing (AS), a post-transcriptional process generating two or more RNA isoforms from one pre-mRNA transcript. We comprehensively analyzed the AS landscape in tomato by integrating genome mapping information of all available mRNA and ESTs with mapping information of RNA-seq reads which were collected from 27 published projects. A total of 369,911 AS events were identified from 34,419 genomic loci involving 161,913 transcripts. Within the basic AS events, intron retention is the prevalent type (18.9%), followed by alternative acceptor site (12.9%) and alternative donor site (7.3%), with exon skipping as the least type (6.0%). Complex AS types having two or more basic event accounted for 54.9% of total AS events. Within 35,768 annotated protein-coding gene models, 23,233 gene models were found having pre-mRNAs generating AS isoform transcripts. Thus the estimated AS rate was 65.0% in tomato. The list of identified AS genes with their corresponding transcript isoforms serves as a catalog for further detailed examination of gene functions in tomato biology. The post-transcriptional information is also expected to be useful in improving the predicted gene models in tomato. The sequence and annotation information can be accessed at plant alternative splicing database (http://proteomics.ysu.edu/altsplice).