Analysis of Transcriptome and Epitranscriptome in Plants Using PacBio Iso-Seq and Nanopore-Based Direct RNA Sequencing

NERD-seq: a novel approach of Nanopore direct RNA sequencing that expands representation of non-coding RNAs

Non-coding RNAs (ncRNAs) are frequently documented RNA modification substrates. Nanopore Technologies enables the direct sequencing of RNAs and the detection of modified nucleobases. Ordinarily, direct RNA sequencing uses polyadenylation selection, studying primarily mRNA gene expression. Here, we present NERD-seq, which enables detection of multiple non-coding RNAs, excluded by the standard approach, alongside natively polyadenylated transcripts. Using neural tissues as a proof of principle, we show that NERD-seq expands representation of frequently modified non-coding RNAs, such as snoRNAs, snRNAs, scRNAs, srpRNAs, tRNAs, and rRFs. NERD-seq represents an RNA-seq approach to simultaneously study mRNA and ncRNA epitranscriptomes in brain tissues and beyond.

Global analysis of key post-transcriptional regulation in early leaf development of Limonium bicolor identifies a long non-coding RNA that promotes salt gland development and salt resistance

Limonium bicolor, known horticulturally as sea lavender, is a typical recretohalophyte with salt glands in its leaf epidermis that secrete excess Na+ out of the plant. Although many genes have been proposed to contribute to salt gland initiation and development, a detailed analysis of alternative splicing, alternative polyadenylation patterns, and long non-coding RNAs (lncRNAs) has been lacking. Here, we applied single-molecule long-read mRNA isoform sequencing (Iso-seq) to explore the complexity of the L. bicolor transcriptome in leaves during salt gland initiation (stage A) and salt gland differentiation (stage B) based on the reference genome. We identified alternative splicing events and the use of alternative poly(A) sites unique to stage A or stage B, leading to the hypothesis that they might contribute to the differentiation of salt glands. Based on the Iso-seq data and RNA in situ hybridization of candidate genes, we selected the lncRNA Btranscript_153392 for gene editing and virus-induced gene silencing to dissect its function. In the absence of this transcript, we observed fewer salt glands on the leaf epidermis, leading to diminished salt secretion and salt tolerance. Our data provide transcriptome resources for unraveling the mechanisms behind salt gland development and furthering crop transformation efforts towards enhanced survivability in saline soils.

Quantitative profiling N1-methyladenosine (m1A) RNA methylation from Oxford nanopore direct RNA sequencing data

With the recent advanced direct RNA sequencing technique that proposed by the Oxford Nanopore Technologies, RNA modifications can be detected and profiled in a simple and straightforward manner. Majority nanopore-based modification studies were devoted to those popular types such as m6A and pseudouridine. To address current limitations on studying the crucial regulator, m1A modification, we conceived this study. We have developed an integrated computational workflow designed for the detection of m1A modifications from direct RNA sequencing data. This workflow comprises a feature extractor responsible for capturing signal characteristics (such as mean, standard deviations, and length of electric signals), a single molecule-level m1A predictor trained with features extracted from the IVT dataset using classical machine learning algorithms, a confident m1A site selector employing the binomial test to identify statistically significant m1A sites, and an m1A modification rate estimator. Our model achieved accurate molecule-level prediction (Average AUC = 0.9689) and reliable m1A site detection and quantification. To show the feasibility of our workflow, we conducted a study on in vivo transcribed human HEK293 cell line, and the results were carefully annotated and compared with other techniques (i.e., Illumina sequencing-based techniques). We believed that this tool will enabling a comprehensive understanding of the m1A modification and its functional mechanisms within cells and organisms.

Transcriptome Analysis Reveals the Immunosuppression in Tiger Pufferfish (Takifugu rubripes) under Cryptocaryon irritans Infection

The tiger pufferfish (Takifugu rubripes), also known as fugu, has recently suffered from severe C. irritans infections under aquaculture environment, yet the underlying immune mechanisms against the parasite remain poorly understood. In this study, we conducted a comprehensive transcriptome analysis of the gill tissue from infected and uninfected fish using PacBio long-read (one pooled sample each for seriously infected and healthy individuals, respectively) and Illumina short-read (three pools for mildly infected, seriously infected, and healthy individuals, respectively) RNA sequencing technologies. After aligning sequence data to fugu's reference genome, 47,307 and 34,413 known full-length transcripts were identified and profiled in healthy and infected fish, respectively. Similarly, we identified and profiled 1126 and 803 novel genes that were obtained from healthy and infected fish, respectively. Interestingly, we found a decrease in the number of alternative splicing (AS) events and long non-coding RNAs (lncRNAs) after infection with C. irritans, suggesting that they may be involved in the regulation of the immune response in fugu. There were 687 and 1535 differentially expressed genes (DEGs) in moderately and heavily infected fish, respectively, compared to uninfected fish. Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that immune-related DEGs in the two comparison groups were mainly enriched in cytokine-cytokine receptor interactions, ECM-receptor interactions, T-cell receptor signaling pathways, Th1 and Th2 cell differentiation, and Th17 cell differentiation pathways. Further analysis revealed that a large number of immune-related genes were downregulated in infected fish relative to uninfected ones, such as CCR7, IL7R, TNFRSF21, CD4, COL2A1, FOXP3B, and ITGA8. Our study suggests that C. irritans is potentially a highly efficient parasite that may disrupt the defense mechanisms of fugu against it. In addition, in combination of short-read RNA sequencing and previous genome-wide association analyses, we identified five key genes (NDUFB6, PRELID1, SMOX, SLC25A4, and DENND1B) that might be closely associated with C. irritans resistance. This study not only provides valuable resources of novel genic transcripts for further research, but also provides new insights into the immune mechanisms underlying C. irritans infection response in farmed fugu.

Full-length transcriptome sequencing of Arabidopsis plants provided new insights into the autophagic regulation of photosynthesis

Autophagy is a highly conserved eukaryotic pathway and plays a crucial role in cell survival under stress conditions. Here, we applied a full-length transcriptome approach to study an Arabidopsis autophagy mutant (atg5-1) subjected to nitrogen-starvation, using Oxford Nanopore Technologies. A total of 39,033 transcripts were identified, including 11,356 new transcripts. In addition, alternative splicing (AS) events and lncRNAs were also detected between Col-0 (WT) and atg5-1. Differentially expressed transcript enrichment showed that autophagy upregulates the expression of many stress-responsive genes and inhibits the transcription of photosynthesis-associated genes. The qRT-PCR results showed that the expression patterns of photosynthesis-related genes in the atg5-1 differed under the conditions of nitrogen starvation and carbon starvation. Under nitrogen starvation treatment, many genes related to photosynthesis also exhibited AS. Chlorophyll fluorescence images revealed that the Fv/Fm and ΦPSII of old atg5-1 leaves were significantly reduced after nitrogen starvation treatment, but the Y(NPQ) indices were significantly increased compared to those of the WT plants. The results of qRT-PCR suggest that autophagy appears to be involved in the degradation of genes related to photodamage repair in PSII. Taken together, the full-length transcriptiome sequencing provide new insights into how new transcripts, lncRNAs and alternative splicing (AS) are involved in plant autophagy through full-length transcriptome sequencing and suggest a new potential link between autophagy and photosynthesis.

Multi Omics Applications in Biological Systems

Traditional methodologies often fall short in addressing the complexity of biological systems. In this regard, system biology omics have brought invaluable tools for conducting comprehensive analysis. Current sequencing capabilities have revolutionized genetics and genomics studies, as well as the characterization of transcriptional profiling and dynamics of several species and sample types. Biological systems experience complex biochemical processes involving thousands of molecules. These processes occur at different levels that can be studied using mass spectrometry-based (MS-based) analysis, enabling high-throughput proteomics, glycoproteomics, glycomics, metabolomics, and lipidomics analysis. Here, we present the most up-to-date techniques utilized in the completion of omics analysis. Additionally, we include some interesting examples of the applicability of multi omics to a variety of biological systems.

Advances in single-cell long-read sequencing technologies

With an increase in accuracy and throughput of long-read sequencing technologies, they are rapidly being assimilated into the single-cell sequencing pipelines. For transcriptome sequencing, these techniques provide RNA isoform-level information in addition to the gene expression profiles. Long-read sequencing technologies not only help in uncovering complex patterns of cell-type specific splicing, but also offer unprecedented insights into the origin of cellular complexity and thus potentially new avenues for drug development. Additionally, single-cell long-read DNA sequencing enables high-quality assemblies, structural variant detection, haplotype phasing, resolving high-complexity regions, and characterization of epigenetic modifications. Given that significant progress has primarily occurred in single-cell RNA isoform sequencing (scRiso-seq), this review will delve into these advancements in depth and highlight the practical considerations and operational challenges, particularly pertaining to downstream analysis. We also aim to offer a concise introduction to complementary technologies for single-cell sequencing of the genome, epigenome and epitranscriptome. We conclude by identifying certain key areas of innovation that may drive these technologies further and foster more widespread application in biomedical science.

Construction of a Full-Length Transcriptome of Western Honeybee Midgut Tissue and Improved Genome Annotation

Honeybees are an indispensable pollinator in nature with pivotal ecological, economic, and scientific value. However, a full-length transcriptome for Apis mellifera, assembled with the advanced third-generation nanopore sequencing technology, has yet to be reported. Here, nanopore sequencing of the midgut tissues of uninoculated and Nosema ceranae-inoculated A. mellifera workers was conducted, and the full-length transcriptome was then constructed and annotated based on high-quality long reads. Next followed improvement of sequences and annotations of the current reference genome of A. mellifera. A total of 5,942,745 and 6,664,923 raw reads were produced from midguts of workers at 7 days post-inoculation (dpi) with N. ceranae and 10 dpi, while 7,100,161 and 6,506,665 raw reads were generated from the midguts of corresponding uninoculated workers. After strict quality control, 6,928,170, 6,353,066, 5,745,048, and 6,416,987 clean reads were obtained, with a length distribution ranging from 1 kb to 10 kb. Additionally, 16,824, 17,708, 15,744, and 18,246 full-length transcripts were respectively detected, including 28,019 nonredundant ones. Among these, 43,666, 30,945, 41,771, 26,442, and 24,532 full-length transcripts could be annotated to the Nr, KOG, eggNOG, GO, and KEGG databases, respectively. Additionally, 501 novel genes (20,326 novel transcripts) were identified for the first time, among which 401 (20,255), 193 (13,365), 414 (19,186), 228 (12,093), and 202 (11,703) were respectively annotated to each of the aforementioned five databases. The expression and sequences of three randomly selected novel transcripts were confirmed by RT-PCR and Sanger sequencing. The 5' UTR of 2082 genes, the 3' UTR of 2029 genes, and both the 5' and 3' UTRs of 730 genes were extended. Moreover, 17,345 SSRs, 14,789 complete ORFs, 1224 long non-coding RNAs (lncRNAs), and 650 transcription factors (TFs) from 37 families were detected. Findings from this work not only refine the annotation of the A. mellifera reference genome, but also provide a valuable resource and basis for relevant molecular and -omics studies.

A Full-Length Transcriptome and Analysis of the NHL-1 Gene Family in Neocaridina denticulata sinensis

Neocaridina denticulata sinensis has emerged as a promising model organism for basic studies in Decapod. However, the current transcriptome information on this species is based on next-generation sequencing technologies, which are limited by a short read length. Therefore, the present study aimed to generate a full-length transcriptome assembly of N. denticulata sinensis utilizing the PacBio Sequel Ⅱ platform. The resulting transcriptome assembly comprised 5831 transcripts with an N50 value of 3697 bp. Remarkably, 90.5% of these transcripts represented novel isoforms of known genes. The transcripts were further searched against the NR, SwissProt, KEGG, KOG, GO, NT, and Pfam databases. A total of 24.8% of the transcripts can be annotated across all seven databases. Additionally, 1236 alternative splicing events, 344 transcription factors, and 124 long non-coding RNAs (LncRNAs) were predicted. Based on the alternative splicing annotation results, a RING finger protein NHL-1 gene from N. denticulata sinensis (NdNHL-1) was identified. There are 15 transcripts in NdNHL-1. The longest transcript is 4995 bp in length and encodes a putative protein of 1665 amino acids. A phylogenetic analysis showed its close relationship with NHL-1 from other crustacean species. This report represents the full-length transcriptome of N. denticulata sinensis and will facilitate research on functional genomics and environmental adaptation in this species.

The Third-Generation Sequencing Challenge: Novel Insights for the Omic Sciences

The understanding of the human genome has been greatly improved by the advent of next-generation sequencing technologies (NGS). Despite the undeniable advantages responsible for their widespread diffusion, these methods have some constraints, mainly related to short read length and the need for PCR amplification. As a consequence, long-read sequencers, called third-generation sequencing (TGS), have been developed, promising to overcome NGS. Starting from the first prototype, TGS has progressively ameliorated its chemistries by improving both read length and base-calling accuracy, as well as simultaneously reducing the costs/base. Based on these premises, TGS is showing its potential in many fields, including the analysis of difficult-to-sequence genomic regions, structural variations detection, RNA expression profiling, DNA methylation study, and metagenomic analyses. Protocol standardization and the development of easy-to-use pipelines for data analysis will enhance TGS use, also opening the way for their routine applications in diagnostic contexts.

Long-read direct RNA sequencing of the mitochondrial transcriptome of Saccharomyces cerevisiae reveals condition-dependent intron abundance

Mitochondria fulfil many essential roles and have their own genome, which is expressed as polycistronic transcripts that undergo co- or posttranscriptional processing and splicing. Due to the inherent complexity and limited technical accessibility of the mitochondrial transcriptome, fundamental questions regarding mitochondrial gene expression and splicing remain unresolved, even in the model eukaryote Saccharomyces cerevisiae. Long-read sequencing could address these fundamental questions. Therefore, a method for the enrichment of mitochondrial RNA and sequencing using Nanopore technology was developed, enabling the resolution of splicing of polycistronic genes and the quantification of spliced RNA. This method successfully captured the full mitochondrial transcriptome and resolved RNA splicing patterns with single-base resolution and was applied to explore the transcriptome of S. cerevisiae grown with glucose or ethanol as the sole carbon source, revealing the impact of growth conditions on mitochondrial RNA expression and splicing. This study uncovered a remarkable difference in the turnover of Group II introns between yeast grown in either mostly fermentative or fully respiratory conditions. Whether this accumulation of introns in glucose medium has an impact on mitochondrial functions remains to be explored. Combined with the high tractability of the model yeast S. cerevisiae, the developed method enables to monitor mitochondrial transcriptome responses in a broad range of relevant contexts, including oxidative stress, apoptosis and mitochondrial diseases.

Epitranscriptomic orchestrations: Unveiling the regulatory paradigm of m6A, A-to-I editing, and m5C in breast cancer via long noncoding RNAs and microRNAs

Breast cancer (BC) poses a persistent global health challenge, particularly in countries with elevated human development indices linked to factors such as increased life expectancy, education, and wealth. Despite therapeutic progress, challenges persist, and the role of epitranscriptomic RNA modifications in BC remains inadequately understood. The epitranscriptome, comprising diverse posttranscriptional modifications on RNA molecules, holds the potential to intricately modulate RNA function and regulation, implicating dysregulation in various diseases, including BC. Noncoding RNAs (ncRNAs), acting as posttranscriptional regulators, influence physiological and pathological processes, including cancer. RNA modifications in long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) add an extra layer to gene expression control. This review delves into recent insights into epitranscriptomic RNA modifications, such as N-6-methyladenosine (m6A), adenine-to-inosine (A-to-I) editing, and 5-methylcytosine (m5C), specifically in the context of lncRNA and miRNAs in BC, highlighting their potential implications in BC development and progression. Understanding this intricate regulatory landscape is vital for deciphering the molecular mechanisms underlying BC and identifying potential therapeutic targets.

Full-length transcriptome analysis of Ophioglossum vulgatum: effects of experimentally identified chloroplast gene clusters on expression and evolutionary patterns

Genes with similar or related functions in chloroplasts are often arranged in close proximity, forming clusters on chromosomes. These clusters are transcribed coordinated to facilitate the expression of genes with specific function. Our previous study revealed a significant negative correlation between the chloroplast gene expression level of the rare medicinal fern Ophioglossum vulgatum and its evolutionary rates as well as selection pressure. Therefore, in this study, we employed a combination of SMRT and Illumina sequencing technology to analyze the full-length transcriptome sequencing of O. vulgatum for the first time. In particular, we experimentally identified gene clusters based on transcriptome data and investigated the effects of chloroplast gene clustering on expression and evolutionary patterns. The results revealed that the total sequenced data volume of the full-length transcriptome of O. vulgatum amounted to 71,950,652,163 bp, and 110 chloroplast genes received transcript coverage. Nine different types of gene clusters were experimentally identified in their transcripts. The chloroplast cluster genes may cause a decrease in non-synonymous substitution rate and selection pressure, as well as a reduction in transversion rate, transition rate, and their ratio. While expression levels of chloroplast cluster genes in leaf, sporangium, and stem would be relatively elevated. The Mann-Whitney U test indicated statistically significant in the selection pressure, sporangia and leaves groups (P < 0.05). We have contributed novel full-length transcriptome data resources for ferns, presenting new evidence on the effects of chloroplast gene clustering on expression land evolutionary patterns, and offering new theoretical support for transgenic research through gene clustering.

Full-Length Transcriptome Sequencing and Comparative Transcriptomic Analyses Provide Comprehensive Insight into Molecular Mechanisms of Flavonoid Metabolites Biosynthesis in Styphnolobium japonicum

Styphnolobium japonicum L. is a commonly consumed plant in China, known for its medicinal and nutritional benefits. This study focuses on the medicinal properties influenced by flavonoid metabolites, which vary during flower development. Utilizing full-length transcriptome sequencing on S. japonicum flowers, we observed changes in gene expression levels as the flowers progressed through growth stages. During stages S1 and S2, key genes related to flavonoid synthesis (PAL, 4CL, CHS, F3H, etc.) exhibited heightened expression. A weighted gene co-expression network analysis (WGCNA) identified regulatory genes (MYB, bHLH, WRKY) potentially involved in the regulatory network with flavonoid biosynthesis-related genes. Our findings propose a regulatory mechanism for flavonoid synthesis in S. japonicum flowers, elucidating the genetic underpinnings of this process. The identified candidate genes present opportunities for genetic enhancements in S. japonicum, offering insights into potential applications for improving its medicinal attributes.

Full-Length Transcriptome Sequencing and RNA-Seq Analysis Offer Insights into Terpenoid Biosynthesis in Blumea balsamifera (L.) DC

Blumea balsamifera (L.) DC., an important economic and medicinal herb, has a long history of being used as a traditional Chinese medicine. Its leaves have always been used as a raw material for the extraction of essential oils, comprising large amounts of terpenoids, which have good therapeutic effects on many diseases, such as eczema, bacterial infection, and hypertension. However, the genetic basis of terpenoid biosynthesis in this plant is virtually unknown on account of the lack of genomic data. Here, a combination of next-generation sequencing (NGS) and full-length transcriptome sequencing was applied to identify genes involved in terpenoid biosynthesis at five developmental stages. Then, the main components of essential oils in B. balsamifera were identified using GC-MS. Overall, 16 monoterpenoids and 20 sesquiterpenoids were identified and 333,860 CCS reads were generated, yielding 65,045 non-redundant transcripts. Among these highly accurate transcripts, 59,958 (92.18%) transcripts were successfully annotated using NR, eggNOG, Swissprot, KEGG, KOG, COG, Pfam, and GO databases. Finally, a total of 56 differently expressed genes (DEGs) involved in terpenoid biosynthesis were identified, including 38 terpenoid backbone genes and 18 TPSs, which provide a significant amount of genetic information for B. balsamifera. These results build a basis for resource protection, molecular breeding, and the metabolic engineering of this plant.

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

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

Chromosome-level genome assembly and annotation of the yellow grouper, Epinephelus awoara

Epinephelus awoara, as known as yellow grouper, is a significant economic marine fish that has been bred artificially in China. However, the genetic structure and evolutionary history of yellow grouper remains largely unknown. Here, this work presents the high-quality chromosome-level genome assembly of yellow grouper using PacBio single molecule sequencing technique (SMRT) and High-through chromosome conformation capture (Hi-C) technologies. The 984.48 Mb chromosome-level genome of yellow grouper was assembled, with a contig N50 length of 39.77 Mb and scaffold N50 length of 41.39 Mb. Approximately 99.76% of assembled sequences were anchored into 24 pseudo-chromosomes with the assistance of Hi-C reads. Furthermore, approximately 41.17% of the genome was composed of repetitive elements. In total, 24,541 protein-coding genes were predicted, of which 22,509 (91.72%) genes were functionally annotated. The highly accurate, chromosome-level reference genome assembly and annotation are crucial to the understanding of population genetic structure, adaptive evolution and speciation of the yellow grouper.

The Evolution, Assembly, and Dynamics of Marine Holobionts

The holobiont concept (i.e., multiple living beings in close symbiosis with one another and functioning as a unit) is revolutionizing our understanding of biology, especially in marine systems. The earliest marine holobiont was likely a syntrophic partnership of at least two prokaryotic members. Since then, symbiosis has enabled marine organisms to conquer all ocean habitats through the formation of holobionts with a wide spectrum of complexities. However, most scientific inquiries have focused on isolated organisms and their adaptations to specific environments. In this review, we attempt to illustrate why a holobiont perspective-specifically, the study of how numerous organisms form a discrete ecological unit through symbiosis-will be a more impactful strategy to advance our understanding of the ecology and evolution of marine life. We argue that this approach is instrumental in addressing the threats to marine biodiversity posed by the current global environmental crisis.

New insights into transcriptome variation during cattle adipocyte adipogenesis by direct RNA sequencing

We performed direct RNA sequencing (DRS) together with PCR-amplified cDNA long and short read sequencing for cattle adipocyte at different stages. We proved that the DRS was with advantages to avoid artificial transcripts and questionable exitrons. Totally, we obtained 68,124 transcripts with information of alternative splicing, poly (A) length and mRNA modification. The number of transcripts for adipogenesis was expanded by alternative splicing, which lead regulation mechanisms far more complex than ever known. We detected 891 differentially expressed genes (DEGs). However, 62.78% transcripts of DEGs were not significantly differentially expressed, and 248 transcripts showed opposite changing directions with their genes. The poly (A) tail became globally shorter in differentiated adipocyte than in primary adipocyte, and had a weak negative correlation with gene/transcript expression. Moreover, the study of different mRNA modifications implied their potential roles in gene expression and alternative splicing. Overall, our study promoted better understanding of adipogenesis mechanisms in cattle adipocytes.

Alternative Splicing Variation: Accessing and Exploiting in Crop Improvement Programs

Alternative splicing (AS) is a gene regulatory mechanism modulating gene expression in multiple ways. AS is prevalent in all eukaryotes including plants. AS generates two or more mRNAs from the precursor mRNA (pre-mRNA) to regulate transcriptome complexity and proteome diversity. Advances in next-generation sequencing, omics technology, bioinformatics tools, and computational methods provide new opportunities to quantify and visualize AS-based quantitative trait variation associated with plant growth, development, reproduction, and stress tolerance. Domestication, polyploidization, and environmental perturbation may evolve novel splicing variants associated with agronomically beneficial traits. To date, pre-mRNAs from many genes are spliced into multiple transcripts that cause phenotypic variation for complex traits, both in model plant Arabidopsis and field crops. Cataloguing and exploiting such variation may provide new paths to enhance climate resilience, resource-use efficiency, productivity, and nutritional quality of staple food crops. This review provides insights into AS variation alongside a gene expression analysis to select for novel phenotypic diversity for use in breeding programs. AS contributes to heterosis, enhances plant symbiosis (mycorrhiza and rhizobium), and provides a mechanistic link between the core clock genes and diverse environmental clues.

MkcDBGAS: a reference-free approach to identify comprehensive alternative splicing events in a transcriptome

Alternative splicing (AS) is an essential post-transcriptional mechanism that regulates many biological processes. However, identifying comprehensive types of AS events without guidance from a reference genome is still a challenge. Here, we proposed a novel method, MkcDBGAS, to identify all seven types of AS events using transcriptome alone, without a reference genome. MkcDBGAS, modeled by full-length transcripts of human and Arabidopsis thaliana, consists of three modules. In the first module, MkcDBGAS, for the first time, uses a colored de Bruijn graph with dynamic- and mixed- kmers to identify bubbles generated by AS with precision higher than 98.17% and detect AS types overlooked by other tools. In the second module, to further classify types of AS, MkcDBGAS added the motifs of exons to construct the feature matrix followed by the XGBoost-based classifier with the accuracy of classification greater than 93.40%, which outperformed other widely used machine learning models and the state-of-the-art methods. Highly scalable, MkcDBGAS performed well when applied to Iso-Seq data of Amborella and transcriptome of mouse. In the third module, MkcDBGAS provides the analysis of differential splicing across multiple biological conditions when RNA-sequencing data is available. MkcDBGAS is the first accurate and scalable method for detecting all seven types of AS events using the transcriptome alone, which will greatly empower the studies of AS in a wider field.

Development of Nanopore sequencing-based full-length transcriptome database toward functional genome annotation of the Pacific oyster, Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is a widely cultivated shellfish in the world, while its transcriptome diversity remains less unexplored due to the limitation of short reads. In this study, we used Oxford Nanopore sequencing to develop the full-length transcriptome database of C. gigas. We identified 77,920 full-length transcripts from 21,523 genes, and uncovered 9668 alternative splicing events and 87,468 alternative polyadenylation sites. Notably, a total of 16,721 novel transcripts were annotated in this work. Furthermore, integrative analysis of 25 publicly available RNA-seq datasets revealed the transcriptome diversity involved in post-transcriptional regulation in C. gigas. We further developed a Drupal based webserver, Cgtdb, which can be used for transcriptome visualization, sequence alignment, and functional genome annotation analyses. This work provides valuable resources and a useful tool for integrative analysis of various transcriptome datasets in C. gigas, which will serve as an essential reference for functional annotation of the oyster genome.

Detection of Transcription Factors Related to Axillary Bud Development after Exposure to Cold Conditions in Hexaploid Chrysanthemum morifolium Using Arabidopsis Information

Chrysanthemum is one of the most commercially used ornamental flowering plants in the world. As chrysanthemum is self-incompatible, the propagation of identical varieties is carried out through cuttings rather than through seed. Axillary bud development can be controlled by changing the temperature; for instance, axillary bud development in some varieties is suppressed at high temperatures. In this study, we focused on the simultaneous axillary bud growth from multiple lines of chrysanthemum upon changing conditions from low to normal temperature. Transcriptome analysis was conducted on the Chrysanthemum morifolium cultivar 'Jinba' to identify the important genes for axillary bud development seen when moved from low-temperature treatment to normal cultivation temperature. We performed RNA-Seq analysis on plants after cold conditions in two-day time-course experiments. Under these settings, we constructed a transcriptome of 415,923 C. morifolium and extracted 7357 differentially expressed genes. Our understanding of Arabidopsis axillary meristem development and growth showed that at least 101 genes in our dataset were homologous to transcription factors involved in the biological process. In addition, six genes exhibited statistically significant variations in expression throughout conditions. We hypothesized that these genes were involved in the formation of axillary buds in C. morifolium after cold conditions.

FASTdRNA: a workflow for the analysis of ONT direct RNA sequencing

Motivation

Direct RNA-seq (dRNA-seq) using Oxford Nanopore Technology (ONT) has revolutionized transcript mapping by offering enhanced precision due to its long-read length. Unlike traditional techniques, dRNA-seq eliminates the need for PCR amplification, reducing the impact of GC bias, and preserving valuable base physical information, such as RNA modification and poly(A) length estimation. However, the rapid advancement of ONT devices has set higher standards for analytical software, resulting in potential challenges of software incompatibility and reduced efficiency.

Results

We present a novel workflow, called FASTdRNA, to manipulate dRNA-seq data efficiently. This workflow comprises two modules: a data preprocessing module and a data analysis module. The preprocessing data module, dRNAmain, encompasses basecalling, mapping, and transcript counting, which are essential for subsequent analyses. The data analysis module consists of a range of downstream analyses that facilitate the estimation of poly(A) length, prediction of RNA modifications, and assessment of alternative splicing events across different conditions with duplication. The FASTdRNA workflow is designed for the Snakemake framework and can be efficiently executed locally or in the cloud. Comparative experiments have demonstrated its superior performance compared to previous methods. This innovative workflow enhances the research capabilities of dRNA-seq data analysis pipelines by optimizing existing processes and expanding the scope of analysis.

Availability and implementation

The workflow is freely available at https://github.com/Tomcxf/FASTdRNA under an MIT license. Detailed install and usage guidance can be found in the GitHub repository.

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.

An Efficient Method to Prepare Barcoded cDNA Libraries from Plant Callus for Long-Read Sequencing

Long-read sequencing methods allow a comprehensive analysis of transcriptomes in identifying full-length transcripts. This revolutionary method represents a considerable breakthrough for non-model species since it allows enhanced gene annotation and gene expression studies when compared to former sequencing methods. However, woody plant tissues are challenging to the successful preparation of cDNA libraries, thus, impairing further cutting-edge sequencing analyses. Here, a detailed protocol for preparing cDNA libraries suitable for high throughput RNA sequencing using Oxford Nanopore Technologies^® is described. This method was used to prepare eight barcoded cDNA libraries from two Solanum betaceum cell lines: one with compact morphology and embryogenic competency (EC) and another with friable and non-embryogenic (NEC). The libraries were successfully sequenced, and data quality assessment showed high mean quality scores. Using this method, long-read sequencing will allow a comprehensive analysis of plant transcriptomes.

Application of third-generation sequencing to herbal genomics

There is a long history of traditional medicine use. However, little genetic information is available for the plants used in traditional medicine, which limits the exploitation of these natural resources. Third-generation sequencing (TGS) techniques have made it possible to gather invaluable genetic information and develop herbal genomics. In this review, we introduce two main TGS techniques, PacBio SMRT technology and Oxford Nanopore technology, and compare the two techniques against Illumina, the predominant next-generation sequencing technique. In addition, we summarize the nuclear and organelle genome assemblies of commonly used medicinal plants, choose several examples from genomics, transcriptomics, and molecular identification studies to dissect the specific processes and summarize the advantages and disadvantages of the two TGS techniques when applied to medicinal organisms. Finally, we describe how we expect that TGS techniques will be widely utilized to assemble telomere-to-telomere (T2T) genomes and in epigenomics research involving medicinal plants.

Comprehensive profiling of epigenetic modifications in fast-growing Moso bamboo shoots

The fast growth of Moso bamboo (Phyllostachys edulis) shoots is caused by the rapid elongation of each internode. However, the key underlying cellular processes and epigenetic mechanisms remain largely unexplored. We used microscopy and multi-omics approaches to investigate two regions (bottom and middle) of the 18th internode from shoots of two different heights (2 and 4 m). We observed that internode cells become longer, and that lignin biosynthesis and glycosyltransferase family 43 (GT43) genes are substantially upregulated with shoot height. Nanopore direct RNA sequencing (DRS) revealed a higher N6-methyladenine (m6A) modification rate in 2-m shoots than in 4-m shoots. In addition, different specific m6A modification sites were enriched at different growth stages. Global DNA methylation profiling indicated that DNA methylation levels are higher in 4-m shoots than in 2-m shoots. We also detected shorter poly(A) tail lengths (PALs) in 4-m shoots compared with 2-m shoots. Genes showing differential PAL were mainly enriched in the functional terms of protein translation and vesicle fusion. An association analysis between PALs and DNA methylation strongly suggested that gene body CG methylation levels are positively associated with PAL. This study provides valuable information to better understand post-transcriptional regulations responsible for fast-growing shoots in Moso bamboo.

Epigenetic and transcriptional activation of the secretory kinase FAM20C as an oncogene in glioma

Gliomas are the most prevalent and aggressive malignancies of the nervous system. Previous bioinformatic studies have revealed the crucial role of the secretory pathway kinase FAM20C in the prediction of glioma invasion and malignancy. However, little is known about the pathogenesis of FAM20C in the regulation of glioma. Here, we construct the full-length transcriptome atlas in paired gliomas and observe that 22 genes are upregulated by full-length transcriptome and differential APA analysis. Analysis of ATAC-seq data reveals that both FAM20C and NPTN are the hub genes with chromatin openness and differential expression. Further, in vitro and in vivo studies suggest that FAM20C stimulates the proliferation and metastasis of glioma cells. Meanwhile, NPTN, a novel cancer suppressor gene, counteracts the function of FAM20C by inhibiting both the proliferation and migration of glioma. The blockade of FAM20C by neutralizing antibodies results in the regression of xenograft tumors. Moreover, MAX, BRD4, MYC, and REST are found to be the potential trans-active factors for the regulation of FAM20C. Taken together, our results uncover the oncogenic role of FAM20C in glioma and shed new light on the treatment of glioma by abolishing FAM20C.

MAEL gene contributes to bovine testicular development through the m5C-mediated splicing

Knowledge of RNA molecules regulating testicular development and spermatogenesis in bulls is essential for elite bull selection and an ideal breeding program. Herein, we performed direct RNA sequencing (DRS) to explore the functional characterization of RNA molecules produced in the testicles of 9 healthy Simmental bulls at three testicular development stages (prepuberty, puberty, and postpuberty). We identified 5,043 differentially expressed genes associated with testicular weight. These genes exhibited more alternative splicing at sexual maturity, particularly alternative 3' (A3) and 5' (A5) splice sites usage and exon skipping (SE). The expression of hub genes in testicular developmental stages was also affected by both m6A and m5C RNA modifications. We found m5C-mediated splicing events significantly (p < 0.05) increased MAEL gene expression at the isoform level, likely promoting spermatogenesis. Our findings highlight the complexity of RNA processing and expression as well as the regulation of transcripts involved in testicular development and spermatogenesis.

Measuring the tail: Methods for poly(A) tail profiling

The 3'-end poly(A) tail is an important and potent feature of most mRNA molecules that affects mRNA fate and translation efficiency. Polyadenylation is a posttranscriptional process that occurs in the nucleus by canonical poly(A) polymerases (PAPs). In some specific instances, the poly(A) tail can also be extended in the cytoplasm by noncanonical poly(A) polymerases (ncPAPs). This epitranscriptomic regulation of mRNA recently became one of the most interesting aspects in the field. Advances in RNA sequencing technologies and software development have allowed the precise measurement of poly(A) tails, identification of new ncPAPs, expansion of the function of known enzymes, discovery and a better understanding of the physiological role of tail heterogeneity, and recognition of a correlation between tail length and RNA translatability. Here, we summarize the development of polyadenylation research methods, including classic low-throughput approaches, Illumina-based genome-wide analysis, and advanced state-of-art techniques that utilize long-read third-generation sequencing with Pacific Biosciences and Oxford Nanopore Technologies platforms. A boost in technical opportunities over recent decades has allowed a better understanding of the regulation of gene expression at the mRNA level. This article is categorized under: RNA Methods > RNA Analyses In Vitro and In Silico.

The Transition from Cancer "omics" to "epi-omics" through Next- and Third-Generation Sequencing

Deciphering cancer etiopathogenesis has proven to be an especially challenging task since the mechanisms that drive tumor development and progression are far from simple. An astonishing amount of research has revealed a wide spectrum of defects, including genomic abnormalities, epigenomic alterations, disturbance of gene transcription, as well as post-translational protein modifications, which cooperatively promote carcinogenesis. These findings suggest that the adoption of a multidimensional approach can provide a much more precise and comprehensive picture of the tumor landscape, hence serving as a powerful tool in cancer research and precision oncology. The introduction of next- and third-generation sequencing technologies paved the way for the decoding of genetic information and the elucidation of cancer-related cellular compounds and mechanisms. In the present review, we discuss the current and emerging applications of both generations of sequencing technologies, also referred to as massive parallel sequencing (MPS), in the fields of cancer genomics, transcriptomics and proteomics, as well as in the progressing realms of epi-omics. Finally, we provide a brief insight into the expanding scope of sequencing applications in personalized cancer medicine and pharmacogenomics.

How do emerging long-read sequencing technologies function in transforming the plant pathology research landscape?

Long-read sequencing technologies are revolutionizing the sequencing and analysis of plant and pathogen genomes and transcriptomes, as well as contributing to emerging areas of interest in plant-pathogen interactions, disease management techniques, and the introduction of new plant varieties or cultivars. Long-read sequencing (LRS) technologies are progressively being implemented to study plants and pathogens of agricultural importance, which have substantial economic effects. The variability and complexity of the genome and transcriptome affect plant growth, development and pathogen responses. Overcoming the limitations of second-generation sequencing, LRS technology has significantly increased the length of a single contiguous read from a few hundred to millions of base pairs. Because of the longer read lengths, new analysis methods and tools have been developed for plant and pathogen genomics and transcriptomics. LRS technologies enable faster, more efficient, and high-throughput ultralong reads, allowing direct sequencing of genomes that would be impossible or difficult to investigate using short-read sequencing approaches. These benefits include genome assembly in repetitive areas, creating more comprehensive and exact genome determinations, assembling full-length transcripts, and detecting DNA and RNA alterations. Furthermore, these technologies allow for the identification of transcriptome diversity, significant structural variation analysis, and direct epigenetic mark detection in plant and pathogen genomic regions. LRS in plant pathology is found efficient for identifying and characterization of effectors in plants as well as known and unknown plant pathogens. In this review, we investigate how these technologies are transforming the landscape of determination and characterization of plant and pathogen genomes and transcriptomes efficiently and accurately. Moreover, we highlight potential areas of interest offered by LRS technologies for future study into plant-pathogen interactions, disease control strategies, and the development of new plant varieties or cultivars.

Identification of nucleoside monophosphates and their epigenetic modifications using an engineered nanopore

RNA modifications play critical roles in the regulation of various biological processes and are associated with many human diseases. Direct identification of RNA modifications by sequencing remains challenging, however. Nanopore sequencing is promising, but the current strategy is complicated by sequence decoding. Sequential nanopore identification of enzymatically cleaved nucleoside monophosphates may simultaneously provide accurate sequence and modification information. Here we show a phenylboronic acid-modified hetero-octameric Mycobacterium smegmatis porin A nanopore, with which direct distinguishing between monophosphates of canonical nucleosides, 5-methylcytidine, N⁶-methyladenosine, N⁷-methylguanosine, N¹-methyladenosine, inosine, pseudouridine and dihydrouridine was achieved. A custom machine learning algorithm, which reports an accuracy of 0.996, was also applied to the quantitative analysis of modifications in microRNA and natural transfer RNA. It is generally suitable for sensing of a variety of other nucleoside or nucleotide derivatives and may bring new insights to epigenetic RNA sequencing.

Epigenetic tumor heterogeneity in the era of single-cell profiling with nanopore sequencing

Nanopore sequencing has brought the technology to the next generation in the science of sequencing. This is achieved through research advancing on: pore efficiency, creating mechanisms to control DNA translocation, enhancing signal-to-noise ratio, and expanding to long-read ranges. Heterogeneity regarding epigenetics would be broad as mutations in the epigenome are sensitive to cause new challenges in cancer research. Epigenetic enzymes which catalyze DNA methylation and histone modification are dysregulated in cancer cells and cause numerous heterogeneous clones to evolve. Detection of this heterogeneity in these clones plays an indispensable role in the treatment of various cancer types. With single-cell profiling, the nanopore sequencing technology could provide a simple sequence at long reads and is expected to be used soon at the bedside or doctor's office. Here, we review the advancements of nanopore sequencing and its use in the detection of epigenetic heterogeneity in cancer.

Novel Method of Full-Length RNA-seq That Expands the Identification of Non-Polyadenylated RNAs Using Nanopore Sequencing

The occurrence of diseases displayed transcriptome alteration, including both coding and non-coding transcripts. The third-generation sequencing (TGS) technologies allow for intensive and comprehensive research of the transcriptome. However, the present standard TGS RNA sequencing method is unable to detect many of the non-polyadenylated [non-poly(A)] RNAs. To obtain more complete transcriptome information, we presented a new comprehensive sequencing approach by performing conventional poly(A) RNA-sequencing combined with the sequencing of non-poly(A) RNA fraction which was tailed by poly(U) on HepG2 and HL-7702 cell lines, enabling the detection of multiple categories of non-poly(A) RNAs excluded by the existing standard approach. Moreover, the length distribution of the full-splice match transcripts was longer than that assembled by short-reads, which contributed to characterizing alternative splicing events and provided a comprehensive portrait of transcriptional complexity. Besides the detection of genes with differential expression patterns in the poly(A) library between HepG2 and HL-7702, we also found a cancer-related non-coding gene in the poly(U) data, that is, growth arrest special 5 (GAS5). Collectively, our results suggested that the novel method effectively captured both poly(A) and non-poly(A) transcripts in the tested cell lines and allowed a deeper exploration of the transcriptome.

Alternative splicing as a key player in the fine-tuning of the immunity response in Arabidopsis

Plants, like animals, are constantly exposed to abiotic and biotic stresses, which often inhibit plant growth and development, and cause tissue damage, disease, and even plant death. Efficient and timely response to stress requires appropriate co- and posttranscriptional reprogramming of gene expression. Alternative pre-mRNA splicing provides an important layer of this regulation by controlling the level of factors involved in stress response and generating additional protein isoforms with specific features. Recent high-throughput studies have revealed that several defence genes undergo alternative splicing that is often affected by pathogen infection. Despite extensive work, the exact mechanisms underlying these relationships are still unclear, but the contribution of alternative protein isoforms to the defence response and the role of regulatory factors, including components of the splicing machinery, have been established. Modulation of gene expression in response to stress includes alternative splicing, chromatin remodelling, histone modifications, and nucleosome occupancy. How these processes affect plant immunity is mostly unknown, but these facets open new regulatory possibilities. Here we provide an overview of the current state of knowledge and recent findings regarding the growing importance of alternative splicing in plant response to biotic stress.

Full-length transcriptome sequencing analysis and characterization, development and validation of microsatellite markers in Kengyilia melanthera

As a typical psammophyte of the Triticeae, Kengyilia melanthera possesses high feeding potential and great utilization values in desertification control in the Qinghai-Tibet Plateau. However, few gene function and genetic studies have been performed in K. melanthera. In this study, single-molecule real-time sequencing technology was used to obtain the full-length transcriptome sequence of K. melanthera, following the functional annotation of transcripts and prediction of coding sequences (CDSs), transcription factors (TFs), and long noncoding RNA (lncRNA) sequences. Meanwhile, a total of 42,433 SSR loci were detected, with 5'-UTRs having the most SSR loci and trinucleotide being the most abundant type. In total, 108,399 SSR markers were designed, and 300 SSR markers were randomly selected for diversity verification of K. melanthera. A total of 49 polymorphic SSR markers were used to construct the genetic relationships of 56 K. melanthera accessions, among which 21 SSR markers showed good cross-species transferability among the related species. In conclusion, the full-length transcriptome sequence of the K. melanthera will assist gene prediction and promote molecular biology and genomics research, and the polymorphic SSR markers will promote molecular-assisted breeding and related research of K. melanthera and its relatives.

Exploring epitranscriptomics for crop improvement and environmental stress tolerance

Climate change and stressful environmental conditions severely hamper crop growth, development and yield. Plants respond to environmental perturbations, through their plasticity provided by key-genes, governed at post-/transcriptional levels. Gene-regulation in plants is a multilevel process controlled by diverse cellular entities that includes transcription factors (TF), epigenetic regulators and non-coding RNAs beside others. There are successful studies confirming the role of epigenetic modifications (DNA-methylation/histone-modifications) in gene expression. Recent years have witnessed emergence of a highly specialized field the "Epitranscriptomics". Epitranscriptomics deals with investigating post-transcriptional RNA chemical-modifications present across the life forms that change structural, functional and biological characters of RNA. However, deeper insights on of epitranscriptomic modifications, with >140 types known so far, are to be understood fully. Researchers have identified epitranscriptome marks (writers, erasers and readers) and mapped the site-specific RNA modifications (m6A, m⁵C, 3' uridylation, etc.) responsible for fine-tuning gene expression in plants. Simultaneous advancement in sequencing platforms, upgraded bioinformatic tools and pipelines along with conventional labelled techniques have further given a statistical picture of these epitranscriptomic modifications leading to their potential applicability in crop improvement and developing climate-smart crops. We present herein the insights on epitranscriptomic machinery in plants and how epitranscriptome and epitranscriptomic modifications underlying plant growth, development and environmental stress responses/adaptations. Third-generation sequencing technology, advanced bioinformatics tools and databases being used in plant epitranscriptomics are also discussed. Emphasis is given on potential exploration of epitranscriptome engineering for crop-improvement and developing environmental stress tolerant plants covering current status, challenges and future directions.

Comparative full-length transcriptome analysis by Oxford Nanopore Technologies reveals genes involved in anthocyanin accumulation in storage roots of sweet potatoes (Ipomoea batatas L.)

Background

Storage roots of sweet potatoes (Ipomoea batatas L.) with different colors vary in anthocyanin content, indicating different economically agronomic trait. As the newest DNA/RNA sequencing technology, Oxford Nanopore Technologies (ONT) have been applied in rapid transcriptome sequencing for investigation of genes related to nutrient metabolism. At present, few reports concern full-length transcriptome analysis based on ONT for study on the molecular mechanism of anthocyanin accumulation leading to color change of tuberous roots of sweet potato cultivars.

Results

The storage roots of purple-fleshed sweet potato (PFSP) and white-fleshed sweet potato (WFSP) at different developmental stages were subjected to anthocyanin content comparison by UV-visible spectroscopy as well as transcriptome analysis at ONT MinION platform. UV-visible spectrophotometric measurements demonstrated the anthocyanin content of PFSP was much higher than that of WFSP. ONT RNA-Seq results showed each sample generated average 2.75 GB clean data with Full-Length Percentage (FL%) over 70% and the length of N50 ranged from 1,192 to 1,395 bp, indicating reliable data for transcriptome analysis. Subsequent analysis illustrated intron retention was the most prominent splicing event present in the resulting transcripts. As compared PFSP with WFSP at the relative developmental stages with the highest (PH vs. WH) and the lowest (PL vs. WL) anthocyanin content, 282 and 216 genes were up-regulated and two and 11 genes were down-regulated respectively. The differential expression genes involved in flavonoid biosynthesis pathway include CCoAOMT, PpLDOX, DFR, Cytochrome P450, CHI, and CHS. The genes encoding oxygenase superfamily were significantly up-regulated when compared PFSP with WFSP at the relative developmental stages.

Conclusions

Comparative full-length transcriptome analysis based on ONT serves as an effective approach to detect the differences in anthocyanin accumulation in the storage roots of different sweet potato cultivars at transcript level, with noting that some key genes can now be closely related to flavonoids biosynthesis. This study helps to improve understanding of molecular mechanism for anthocyanin accumulation in sweet potatoes and also provides a theoretical basis for high-quality sweet potato breeding.

In-Depth Temporal Transcriptome Profiling of an Alphaherpesvirus Using Nanopore Sequencing

In this work, a long-read sequencing (LRS) technique based on the Oxford Nanopore Technology MinION platform was used for quantifying and kinetic characterization of the poly(A) fraction of bovine alphaherpesvirus type 1 (BoHV-1) lytic transcriptome across a 12-h infection period. Amplification-based LRS techniques frequently generate artefactual transcription reads and are biased towards the production of shorter amplicons. To avoid these undesired effects, we applied direct cDNA sequencing, an amplification-free technique. Here, we show that a single promoter can produce multiple transcription start sites whose distribution patterns differ among the viral genes but are similar in the same gene at different timepoints. Our investigations revealed that the circ gene is expressed with immediate–early (IE) kinetics by utilizing a special mechanism based on the use of the promoter of another IE gene (bicp4) for the transcriptional control. Furthermore, we detected an overlap between the initiation of DNA replication and the transcription from the bicp22 gene, which suggests an interaction between the two molecular machineries. This study developed a generally applicable LRS-based method for the time-course characterization of transcriptomes of any organism.

Full-Length Transcriptional Analysis of the Same Soybean Genotype With Compatible and Incompatible Reactions to Heterodera glycines Reveals Nematode Infection Activating Plant Defense Response

Full-length transcriptome sequencing with long reads is a powerful tool to analyze transcriptional and post-transcriptional events; however, it has not been applied on soybean (Glycine max). Here, a comparative full-length transcriptome analysis was performed on soybean genotype 09-138 infected with soybean cyst nematode (SCN, Heterodera glycines) race 4 (SCN4, incompatible reaction) and race 5 (SCN5, compatible reaction) using Oxford Nanopore Technology. Each of 9 full-length samples collected 8 days post inoculation with/without nematodes generated an average of 6.1 GB of clean data and a total of 65,038 transcript sequences. After redundant transcripts were removed, 1,117 novel genes and 41,096 novel transcripts were identified. By analyzing the sequence structure of the novel transcripts, a total of 28,759 complete open reading frame (ORF) sequences, 5,337 transcription factors, 288 long non-coding RNAs, and 40,090 novel transcripts with function annotation were predicted. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of differentially expressed genes (DEGs) revealed that growth hormone, auxin-activated signaling pathway and multidimensional cell growth, and phenylpropanoid biosynthesis pathway were enriched by infection with both nematode races. More DEGs associated with stress response elements, plant-hormone signaling transduction pathway, and plant-pathogen interaction pathway with more upregulation were found in the incompatible reaction with SCN4 infection, and more DEGs with more upregulation involved in cell wall modification and carbohydrate bioprocess were detected in the compatible reaction with SCN5 infection when compared with each other. Among them, overlapping DEGs with a quantitative difference was triggered. The combination of protein-protein interaction with DEGs for the first time indicated that nematode infection activated the interactions between transcription factor WRKY and VQ (valine-glutamine motif) to contribute to soybean defense. The knowledge of the SCN-soybean interaction mechanism as a model will present more understanding of other plant-nematode interactions.

Anthocyanin regulatory networks in Solanum tuberosum L. leaves elucidated via integrated metabolomics, transcriptomics, and StAN1 overexpression

BACKGROUND

Anthocyanins, which account for color variation and remove reactive oxygen species, are widely synthesized in plant tissues and organs. Using targeted metabolomics and nanopore full-length transcriptomics, including differential gene expression analysis, we aimed to reveal potato leaf anthocyanin biosynthetic pathways in different colored potato varieties.

RESULTS

Metabolomics analysis revealed 17 anthocyanins. Their levels varied significantly between the different colored varieties, explaining the leaf color differences. The leaves of the Purple Rose2 (PurpleR2) variety contained more petunidin 3-O-glucoside and malvidin 3-O-glucoside than the leaves of other varieties, whereas leaves of Red Rose3 (RedR3) contained more pelargonidin 3-O-glucoside than the leaves of other varieties. In total, 114 genes with significantly different expression were identified in the leaves of the three potato varieties. These included structural anthocyanin synthesis-regulating genes such as F3H, CHS, CHI, DFR, and anthocyanidin synthase and transcription factors belonging to multiple families such as C3H, MYB, ERF, NAC, bHLH, and WRKY. We selected an MYB family transcription factor to construct overexpression tobacco plants; overexpression of this factor promoted anthocyanin accumulation, turning the leaves purple and increasing their malvidin 3-o-glucoside and petunidin 3-o-glucoside content.

CONCLUSIONS

This study elucidates the effects of anthocyanin-related metabolites on potato leaves and identifies anthocyanin metabolic network candidate genes.

Full-Length Transcriptome Analysis of the Halophyte Nitraria sibirica Pall

Background: Nitraria sibirica Pall. is one of the pioneer tree species in saline–alkali areas due to its extreme salt tolerance. However, the lack of information on its genome limits the further exploration of the molecular mechanisms in N. sibirica under salt stress. Methods: In this study, we used single-molecule real-time (SMRT) technology based on the PacBio Iso-Seq platform to obtain transcriptome data from N. sibirica under salt treatment for the first time, which is helpful for our in-depth analysis of the salt tolerance and molecular characteristics of N. sibirica. Results: Our results suggested that a total of 234,508 circular consensus sequences (CCSs) with a mean read length of 2121 bp were obtained from the 19.26 Gb raw data. Furthermore, based on transcript cluster analysis, 93,713 consensus isoforms were obtained, including 92,116 high-quality isoforms. After removing redundant sequences, 49,240 non-redundant transcripts were obtained from high-quality isoforms. A total of 37,261 SSRs, 1816 LncRNAs and 47,314 CDSs, of which 40,160 carried complete ORFs, were obtained. Based on our transcriptome data, we also analyzed the coding genes of H+-PPase, and the results of both bioinformatics and functional analyses indicated that the gene prediction via full-length transcripts obtained by SMRT technology is reliable and effective. In summary, our research data obtained by SMRT technology provides more reliable and accurate information for the further analysis of the regulatory network and molecular mechanism of N. sibirica under salt stress.

Full-length transcriptome sequencing reveals the molecular mechanism of potato seedlings responding to low-temperature

BACKGROUND

Potato (Solanum tuberosum L.) is one of the world's most important crops, the cultivated potato is frost-sensitive, and low-temperature severely influences potato production. However, the mechanism by which potato responds to low-temperature stress is unclear. In this research, we apply a combination of second-generation sequencing and third-generation sequencing technologies to sequence full-length transcriptomes in low-temperature-sensitive cultivars to identify the important genes and main pathways related to low-temperature resistance.

RESULTS

In this study, we obtained 41,016 high-quality transcripts, which included 15,189 putative new transcripts. Amongst them, we identified 11,665 open reading frames, 6085 simple sequence repeats out of the potato dataset. We used public available genomic contigs to analyze the gene features, simple sequence repeat, and alternative splicing event of 24,658 non-redundant transcript sequences, predicted the coding sequence and identified the alternative polyadenylation. We performed cluster analysis, GO, and KEGG functional analysis of 4518 genes that were differentially expressed between the different low-temperature treatments. We examined 36 transcription factor families and identified 542 transcription factors in the differentially expressed genes, and 64 transcription factors were found in the AP2 transcription factor family which was the most. We measured the malondialdehyde, soluble sugar, and proline contents and the expression genes changed associated with low temperature resistance in the low-temperature treated leaves. We also tentatively speculate that StLPIN10369.5 and StCDPK16 may play a central coordinating role in the response of potatoes to low temperature stress.

CONCLUSIONS

Overall, this study provided the first large-scale full-length transcriptome sequencing of potato and will facilitate structure-function genetic and comparative genomics studies of this important crop.

Nanopore sequencing of RNA and cDNA molecules in Escherichia coli

High-throughput sequencing dramatically changed our view of transcriptome architectures and allowed for ground-breaking discoveries in RNA biology. Recently, sequencing of full-length transcripts based on the single-molecule sequencing platform from Oxford Nanopore Technologies (ONT) was introduced and is widely used to sequence eukaryotic and viral RNAs. However, experimental approaches implementing this technique for prokaryotic transcriptomes remain scarce. Here, we present an experimental and bioinformatic workflow for ONT RNA-seq in the bacterial model organism Escherichia coli, which can be applied to any microorganism. Our study highlights critical steps of library preparation and computational analysis and compares the results to gold standards in the field. Furthermore, we comprehensively evaluate the applicability and advantages of different ONT-based RNA sequencing protocols, including direct RNA, direct cDNA, and PCR-cDNA. We find that (PCR)-cDNA-seq offers improved yield and accuracy compared to direct RNA sequencing. Notably, (PCR)-cDNA-seq is suitable for quantitative measurements and can be readily used for simultaneous and accurate detection of transcript 5' and 3' boundaries, analysis of transcriptional units, and transcriptional heterogeneity. In summary, based on our comprehensive study, we show nanopore RNA-seq to be a ready-to-use tool allowing rapid, cost-effective, and accurate annotation of multiple transcriptomic features. Thereby nanopore RNA-seq holds the potential to become a valuable alternative method for RNA analysis in prokaryotes.

Allele-aware chromosome-scale assembly of the allopolyploid genome of hexaploid Ma bamboo (Dendrocalamus latiflorus Munro)

Dendrocalamus latiflorus Munro is a woody clumping bamboo with rapid shoot growth. Both genetic transformation and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) gene editing techniques are available for D. latiflorus, enabling reverse genetic approaches. Thus, D. latiflorus has the potential to be a model bamboo species. However, the genome sequence of D. latiflorus has remained unreported due to its polyploidy and large genome size. Here, we sequenced the D. latiflorus genome and assembled it into three allele-aware subgenomes (AABBCC), representing the largest genome of a major bamboo species. We assembled 70 allelic chromosomes (2, 737 Mb) for hexaploid D. latiflorus using both single-molecule sequencing from the Pacific Biosciences (PacBio) Sequel platform and chromosome conformation capture sequencing (Hi-C). Repetitive sequences comprised 52.65% of the D. latiflorus genome. We annotated 135 231 protein-coding genes in the genome based on transcriptomes from eight different tissues. Transcriptome sequencing using RNA-Seq and PacBio single-molecule real-time long-read isoform sequencing revealed highly differential alternative splicing (AS) between non-abortive and abortive shoots, suggesting that AS regulates the abortion rate of bamboo shoots. This high-quality hexaploid genome and comprehensive strand-specific transcriptome datasets for this Poaceae family member will pave the way for bamboo research using D. latiflorus as a model species.