ANGLE: a sequencing errors resistant program for predicting protein coding regions in unfinished cDNA

Integrated Transcriptome and Proteome Analysis Provides New Insights into Starch and Sucrose Metabolism and Regulation of Corm Expansion Process in Colocasia esculenta

The global significance of Colocasia esculenta, a tuber crop rich in nutritional value and starch, prompts further investigation into its corm development.

BACKGROUND

Previous studies have focused on starch accumulation within the tubers, yet the genetic and proteomic basis of corm expansion remains largely unexplored. This study aims to elucidate the key genes and proteins involved in this process.

METHODS

We selected 'Lipu Taro No.1' and conducted a longitudinal starch content analysis, full-length transcriptome sequencing, and a proteomic analysis during three distinct stages of corm development.

RESULTS

Our findings reveal a significant increase in both amylose and amylopectin contents as the corm develops, indicating the temporal regulation of starch biosynthesis. The integration of transcriptome and proteomic data identified differentially expressed genes and proteins associated with starch and sucrose metabolism, as well as plant hormone signal transduction.

CONCLUSIONS

This study delineates a temporal gene expression pattern that is crucial for starch synthesis and provides insights into the regulatory mechanisms controlling corm expansion and starch deposition, offering valuable references for future molecular breeding strategies to enhance taro yield and quality.

Full-length transcriptome sequencing of the longissimus dorsi muscle of yak and cattle-yak using nanopore technology

Short-read RNA sequencing has been used to sequence the transcriptome of the skeletal muscle of yak and cattle-yak; however, full-length transcripts cannot be obtained and alternative splicing (AS) events cannot be inferred using this sequencing approach. Here, we used Oxford Nanopore Technologies (ONT) full-length sequencing to sequence the transcriptome of the longissimus dorsi of yak and cattle-yak. A total of 20,323 novel genes and 172,870 novel transcripts were identified, and 159,700 novel transcripts were successfully annotated. A total of 157,812 AS events, 58,073 simple sequence repeats, 57,468 complete open reading frames, 2296 transcription factors, and 20,404 lncRNAs were detected. Differentially expressed transcripts (DETs) in the longissimus dorsi muscle of yak and cattle-yak were involved in the MAPK and JAK-STAT signaling pathways related to muscle development and growth. Protein-protein interaction analysis of DETs suggested that TNNI2 might make a major contribution to differences in muscle growth and meat quality traits between yak and cattle-yak. The results have enriched the transcriptome data of dorsal muscles, providing new ideas for the study of transcriptional regulation processes, and also providing useful information for the production of higher yields of yak meat.

Construction of a Full-Length transcriptome resource for the African sharptooth catfish (Clarias gariepinus), a prototypical air-breathing Fish, based on isoform sequencing (Iso-Seq)

The African sharptooth catfish (Clarias gariepinus) assumes significance in aquaculture, given its role as a farmed freshwater species with modified gill structures functioning as an air-breathing organ (ABO). To provide a scientific basis for further elucidating the air-breathing formation mechanism and deeply utilizing the genetic resources of Clarias gariepinus, we utilized the PacBio sequencing platform to acquire a comprehensive full-length transcriptome from five juvenile developmental stages and various adult tissues, including the ABO, gills, liver, skin, and muscle. We generated 25,766,688 high-quality reads, with an average length of 2,006 bp and an N50 of 2,241 bp. Following rigorous quality control, 34,890 (97.7 %) of the high-quality isoforms were mapped to the reference genome for gene and transcript annotation, yielding 387 novel isoforms and 14,614 new isoforms. Additionally, we identified 28,582 open reading frames, 48 SNPs, 5,464 variable splices, and 6,141 variable polyadenylation sites, along with 475 long non-coding RNAs. Many DEGs were involved with low oxygen GO terms and KEGG pathways, such as response to stimulus, biological regulation and catalytic activities. Furthermore, it was found that transcription factors such as zf-C2H2, Homeobox, bHLH, and MYB could underpin the African sharptooth catfish's developmental plasticity and its capacity to adapt its morphology and function to its environment. Through the comprehensive analysis of its genomic characteristics, it was found that the African sharptooth catfish has developed a series of unique respiratory adaptive mechanisms during the evolutionary process, These results not only advances the understanding of genetic adaptations to hypoxia in Clarias fish but also provides a valuable framework for future studies aimed at improving aquaculture practices,besides provide important references and inspirations for the evolution of aquatic organisms.

Unraveling Pinus massoniana's Defense Mechanisms Against Bursaphelenchus xylophilus Under Aseptic Conditions: A Transcriptomic Analysis

Pine wilt disease (PWD) is caused by the pine wood nematode (PWN, Bursaphelenchus xylophilus) and significantly impacts pine forest ecosystems globally. This study focuses on Pinus massoniana, an important timber and oleoresin resource in China, which is highly susceptible to PWN. However, the defense mechanism of pine trees in response to PWN remains unclear. Addressing the complexities of PWD, influenced by diverse factors such as bacteria, fungi, and environment, we established a reciprocal system between PWN and P. massoniana seedlings under aseptic conditions. Utilizing combined second- and third-generation sequencing technologies, we identified 3,718 differentially expressed genes post PWN infection. Transcript analysis highlighted the activation of defense mechanisms via stilbenes, salicylic acid and jasmonic acid pathways, terpene synthesis, and induction of pathogenesis-related proteins and resistance genes, predominantly at 72 h postinfection. Notably, terpene synthesis pathways, particularly the mevalonate pathway, were crucial in defense, suggesting their significance in P. massoniana's response to PWN. This comprehensive transcriptome profiling offers insights into P. massoniana's intricate defense strategies against PWN under aseptic conditions, laying a foundation for future functional analyses of key resistance genes. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Exploring the innate immune system of Urechis unicinctus: Insights from full-length transcriptome analysis

The Echiura worm Urechis unicinctus refers to a common benthic invertebrate found in the intertidal zone of Huanghai as well as Bohai Bay. U. unicinctus is known to contain various physiologically active substances, making it highly valuable in terms of its edibility, medicinal properties, and economic potential. Nonetheless, the limited study on the immune system of U. unicinctus poses difficulties for its aquaculture and artificial reproduction. Marine invertebrates, including shellfish and U. unicinctus, are thought to primarily depend on their innate immune system for disease protection, owing to the severalinnate immune molecules they possess. Herein, we employed PacBio single-molecule real-time (SMRT) sequencing technology to perform the full-length transcriptome analysis of U. unicinctus individuals under five different conditions (room temperature (RT), low temperature (LT), high temperature (HT), without water (DRY), ultraviolet irradiation (UV)). Concequently, we identified 59,371 unigenes that had a 2,779 bp average length, 2,613 long non-coding RNAs (lncRNAs), 59,190 coding sequences (CDSs), 35,166 simple sequence repeats (SSRs), and 1,733 transcription factors (TFs), successfully annotating 90.58 % (53,778) of the unigenes. Subsequently, key factors associated with immune-related processes, such as non-self-recognition, cellular immune defenses, and humoral immune defenses, were searched. Our study also identified pattern recognition receptors (PRRs) that included 17 peptidoglycan recognition proteins (PGRPs), 13 Gram-negative binding proteins (GNBPs), 18 scavenger receptors (SRs), 74 toll-like receptors (TLRs), and 89 C-type lectins (CLTs). Altogether, the high-quality transcriptome obtained data will offer valuable insights for further investigations into U. unicinctus innate immune response, laying the foundation for subsequent molecular biology studies and aquaculture.

Full-length transcriptome characterization and analysis of Carrizo Citrange and molecular insights into pathogen defense

Citrus huanglongbing (HLB) is a major challenge that impacts the flourishing of the citrus industry. Therefore, analyzing the genomic information of HLB-resistant or tolerant citrus resources is crucial for breeding HLB-resistant citrus varieties. The Carrizo citrange, a hybrid of Citrus sinensis and Poncirus trifoliata, plays a pivotal role in citrus cultivation. However, its genetic explorations are difficult due to the absence of a reference genome or full-length transcriptome. In order to enhance our understanding of the genetic information of citrange, we conducted a full-length transcriptomic sequencing of multiple tissues from the Carrizo citrange using the PacBio Sequel II platform. Moreover, we performed gene ontology (GO) annotation, gene functional annotation, simple sequence repeats (SSR) types analysis, as well as identification of lncRNAs, alternative splicing events, and analysis of pathogen defense-related genes. Results showed that a total of 43,452 isoforms were generated, with 43,307 of them being annotated. GO annotation indicated the involvement of these isoforms in various biological processes, cellular components, and molecular functions. The coding sequence length of the isoforms ranged from 1,000 to 4,000 base pairs (bp). Moreover, we have discovered 54 varieties of transcription factors and regulators, along with 16 classifications of genes associated with resistance. Among all types of SSRs, trimer type SSRs were the most abundant. 130 lncRNAs were predicted to be highly reliable in the isoforms of the Carrizo citrange, with alternative splicing events identified, and the most frequent being retained intron. The analysis of gene family expansion and contraction revealed a significant increase in pathogen defense-related genes within the Carrizo citrange. The results of this study will be of great value for future investigations into gene function in citrange and for expanding the genetic pool for breeding citrus varieties resistant or tolerant to HLB.

Identifying resistance molecules in TiO2 nanoparticle-tolerant strains to facilitate the development of strategies for combating TiO2 nanoparticle pollution

The severity of environmental pollution caused by TiO2 nanoparticles (nTiO2) is increasing, highlighting the urgent need for the development of strategies to combat nTiO2 pollution. Insights into resistance molecules from nTiO2-tolerant strains may facilitate such development. In this study, we utilized multi-omics, genetic manipulation, physiological and biochemical experiments to identify relevant resistance molecules in two strains (Physarum polycephalum Z259 and T83) tolerated to mixed-phase nTiO2 (MPnTiO2). We discovered that a competing endogenous RNA (ceRNA) network, comprising one long non-coding RNA (lncRNA), four microRNAs, and nine mRNAs, influenced metabolic rearrangement and was associated with significant resistance in these strains. Additionally, we found that the lncRNA in the ceRNAs network and certain small-weight metabolites associated with the ceRNA exhibited notable mitigation effects not only against MPnTiO2 but also against other types of nTiO2 with broad species applicability (they significantly improved the resistance of several non-nTiO2-tolerant cells/organisms in the laboratory and reduced cell damage of non-nTiO2-tolerant cells/organisms in highly suspected nTiO2-polluted areas of the real world). In summary, this study deepens our understanding of nTiO2-tolerant strains, provides valuable insights into resistance molecules in these strains, and facilitates the development of strategies to combat nTiO2 pollution.

Analysis of the Spatholobus suberectus full-length transcriptome identified an R2R3-MYB transcription factor-encoding gene SsMYB158 that regulates flavonoid biosynthesis

Spatholobus suberectus Dunn (Leguminosae) has been used for medicinal purposes for a long period. Flavonoids are the major bioactive components of S. suberectus. However, there is still limited knowledge of the exact method via which transcription factors (TFs) regulate flavonoid biosynthesis. The full-length transcriptome of S. suberectus was analyzed using SMRT sequencing; 61,548 transcripts were identified, including 12,311 new gene loci, 53,336 novel transcripts, 44,636 simple sequence repeats, 36,414 complete coding sequences, 871 long non-coding RNAs and 6781 TFs. The SsMYB158 TF, which is associated with flavonoid biosynthesis, belongs to the R2R3-MYB class and is localized subcellularly to the nucleus. The overexpression of SsMYB158 in Nicotiana benthamiana and the transient overexpression of SsMYB158 in S. suberectus resulted in a substantial enhancement in both flavonoids and catechin levels. In addition, there was a remarkable upregulation in the expression of essential enzyme-coding genes associated with the flavonoid biosynthesis pathways. Our study revealed SsMYB158 as a critical regulator of flavonoid biosynthesis in S. suberectus and laying the foundation for its molecular breeding.

Transcriptome data for an ancient 'living-fossil' mollusc, Entemnotrochus rumphii

The Pleurotomarioidea, commonly referred to as slit shells, constitute one of the most ancient and enduring lineages within the phylum Mollusca, with its fossil record tracing back to the Upper Cambrian epoch. Its rareness and evolutionary antiquity surpass even that of the nautilus. In this study, we conducted the first transcriptome sequencing and analyses of Entemnotrochus rumphii (Schepman, 1879), a representative species of Pleurotomarioidea. Full-length transcriptome sequencing of E. rumphii was performed using the PacBio Sequel II platform with SMRT technology. A total of 64.38 gigabytes of data and 964,550 polymerase reads were generated, resulting in 28,068,998 subreads after data filtering. After de-duplication, correction, and clustering, we identified 19,273 genes. Additionally, next-generation sequencing was performed on 11 tissues of E. rumphii. This investigation provides a detailed portrayal and analytical scrutiny of its transcriptomic landscape.

Full-length agave transcriptome reveals candidate glycosyltransferase genes involved in hemicellulose biosynthesis

Agave species are typical crassulacean acid metabolism (CAM) plants commonly cultivated to produce beverages, fibers, and medicines. To date, few studies have examined hemicellulose biosynthesis in Agave H11648, which is the primary cultivar used for fiber production. We conducted PacBio sequencing to obtain full-length transcriptome of five agave tissues: leaves, shoots, roots, flowers, and fruits. A total of 41,807 genes were generated, with a mean length of 2394 bp and an annotation rate of 97.12 % using public databases. We identified 42 glycosyltransferase genes related to hemicellulose biosynthesis, including mixed-linkage glucan (1), glucomannan (5), xyloglucan (16), and xylan (20). Their expression patterns were examined during leaf development and fungal infection, together with hemicellulose content. The results revealed four candidate glycosyltransferase genes involved in xyloglucan and xylan biosynthesis, including glucan synthase (CSLC), xylosyl transferase (XXT), xylan glucuronyltransferase (GUX), and xylan α-1,3-arabinosyltransferase (XAT). These genes can be potential targets for manipulating xyloglucan and xylan traits in agaves, and can also be used as candidate enzymatic tools for enzyme engineering. We have provided the first full-length transcriptome of agave, which will be a useful resource for gene identification and characterization in agave species. We also elucidated the hemicellulose biosynthesis machinery, which will benefit future studies on hemicellulose traits in agave.

Unraveling the molecular landscape of breast muscle development in domestic Yuzhong pigeons and European meat pigeon: Insights from Iso-seq and RNA-seq analysis

The mechanisms governing gene regulation in domestic Yuzhong pigeon breast muscle development remain largely elusive. Here, we conducted a comparative analysis using Iso-seq and RNA-seq data from domestic Yuzhong pigeons and European meat pigeons to uncover signaling pathways and genes involved in breast muscle development. The Iso-seq data from domestic Yuzhong pigeons yielded 131,377,075 subreads, resulting in 16,587 non-redundant high-quality full-length transcripts post-correction. Furthermore, utilizing pfam, CPC, PLEK, and CPAT, we predicted 5575, 4973, 2333, and 4336 lncRNAs, respectively. Notably, several genes potentially implicated in breast muscle development were identified, including tropomyosin beta chain, myosin regulatory light chain 2, and myosin binding protein C. KEGG enrichment analysis revealed critical signaling pathways in breast muscle development, spanning carbon metabolism, biosynthesis of amino acids, glycolysis/gluconeogenesis, estrogen signaling, PI3K-AKT signaling, protein processing in the endoplasmic reticulum, oxidative phosphorylation, pentose phosphate pathway, fructose and mannose metabolism, and tight junctions. These findings offer insights into the biological processes driving breast muscle development in domestic Yuzhong pigeon, contributing to our understanding of this complex phenomenon.

Identification of key genes for triacylglycerol biosynthesis and storage in herbaceous peony (Paeonia lactifolra Pall.) seeds based on full-length transcriptome

BACKGROUND

The herbaceous peony (Paeonia lactiflora Pall.) is extensively cultivated in China due to its root being used as a traditional Chinese medicine known as 'Radix Paeoniae Alba'. In recent years, it has been discovered that its seeds incorporate abundant unsaturated fatty acids, thereby presenting a potential new oilseed plant. Surprisingly, little is known about the full-length transcriptome sequencing of Paeonia lactiflora, limiting research into its gene function and molecular mechanisms.

RESULTS

A total of 484,931 Reads of Inserts (ROI) sequences and 1,455,771 full-Length non-chimeric reads (FLNC) sequences were obtained for CDS prediction, TF analysis, SSR analysis and lncRNA identification. In addition, gene function annotation and gene structure analysis were performed. A total of 4905 transcripts were related to lipid metabolism biosynthesis pathway, belonging to 28 enzymes. We use these data to identify 10 oleosin (OLE) and 5 diacylglycerol acyltransferase (DGAT) gene members after de-redundancy. The analysis of physicochemical properties and secondary structure showed them similarity in gene family respectively. The phylogenetic analysis showed that the distribution of OLE and DGAT family members was roughly the same as that of Arabidopsis. Quantitative real-time polymerase chain reaction (qRT-PCR) analyses revealed expression changes in different seed development stages, and showed a trend of increasing and then decreasing.

CONCLUSION

In summary, these results provide new insights into the molecular mechanism of triacylglycerol (TAG) biosynthesis and storage during the seedling stage in Paeonia lactiflora. It provides theoretical references for selecting and breeding oil varieties and understanding the functions of oil storage as well as lipid synthesis related genes in Paeonia lactiflora.

Transcriptome sequencing and expression analysis in peanut reveal the potential mechanism response to Ralstonia solanacearum infection

BACKGROUND

Bacterial wilt caused by Ralstonia solanacearum severely affects peanut (Arachis hypogaea L.) yields. The breeding of resistant cultivars is an efficient means of controlling plant diseases. Therefore, identification of resistance genes effective against bacterial wilt is a matter of urgency. The lack of a reference genome for a resistant genotype severely hinders the process of identification of resistance genes in peanut. In addition, limited information is available on disease resistance-related pathways in peanut.

RESULTS

Full-length transcriptome data were used to generate wilt-resistant and -susceptible transcript pools. In total, 253,869 transcripts were retained to form a reference transcriptome for RNA-sequencing data analysis. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis of differentially expressed genes revealed the plant-pathogen interaction pathway to be the main resistance-related pathway for peanut to prevent bacterial invasion and calcium plays an important role in this pathway. Glutathione metabolism was enriched in wilt-susceptible genotypes, which would promote glutathione synthesis in the early stages of pathogen invasion. Based on our previous quantitative trait locus (QTL) mapping results, the genes arahy.V6I7WA and arahy.MXY2PU, which encode nucleotide-binding site-leucine-rich repeat receptor proteins, were indicated to be associated with resistance to bacterial wilt.

CONCLUSIONS

This study identified several pathways associated with resistance to bacterial wilt and identified candidate genes for bacterial wilt resistance in a major QTL region. These findings lay a foundation for investigation of the mechanism of resistance to bacterial wilt in peanut.

Full-length transcriptome characterization of Platycladus orientalis based on the PacBio platform

As a unique and native conifer in China, Platycladus orientalis is widely used in soil erosion control, garden landscapes, timber, and traditional Chinese medicine. However, due to the lack of reference genome and transcriptome, it is limited to the further molecular mechanism research and gene function mining. To develop a full-length reference transcriptome, tissues from five different parts of P. orientalis and four cone developmental stages were sequenced and analyzed by single-molecule real-time (SMRT) sequencing through the PacBio platform in this study. Overall, 37,111 isoforms were detected by PacBio with an N50 length of 2,317 nt, an average length of 1,999 bp, and the GC content of 41.81%. Meanwhile, 36,120 coding sequences, 5,645 simple sequence repeats (SSRs), 1,201 non-coding RNAs (lncRNAs), and 182 alternative splicing (AS) events with five types were identified using the results obtained from the PacBio transcript isoforms. Furthermore, 1,659 transcription factors (TFs) were detected and belonged to 51 TF families. A total of 35,689 transcripts (96.17%) were annotated through the NCBI nr, KOG, Swiss-Prot and KEGG databases, and 385 transcript isoforms related to 8 types of hormones were identified incorporated into plant hormone signal transduction pathways. The assembly and revelation of the full-length transcriptome of P. orientalis offer a pioneering insight for future investigations into gene function and genetic breeding within Platycladus species.

Postmortem Muscle Proteome Characteristics of Silver Carp (Hypophthalmichthys molitrix): Insights from Full-Length Transcriptome and Deep 4D Label-Free Proteomic

The coverage of the protein database directly determines the results of shotgun proteomics. In this study, PacBio single-molecule real-time sequencing technology was performed on postmortem silver carp muscle transcripts. A total of 42.43 Gb clean data, 35,834 nonredundant transcripts, and 15,413 unigenes were obtained. In total, 99.32% of the unigenes were successfully annotated and assigned specific functions. PacBio long-read isoform sequencing (Iso-Seq) analysis can provide more accurate protein information with a higher proportion of complete coding sequences and longer lengths. Subsequently, 2671 proteins were identified in deep 4D proteomics informed by a full-length transcriptomics technique, which has been shown to improve the identification of low-abundance muscle proteins and potential protein isoforms. The feature of the sarcomeric protein profile and information on more than 30 major proteins in the white dorsal muscle of silver carp were reported here for the first time. Overall, this study provides valuable transcriptome data resources and the comprehensive muscle protein information detected to date for further study into the processing characteristic of early postmortem fish muscle, as well as a spectral library for data-independent acquisition and data processing. This batch of muscle-specific dependent acquisition data is available via PRIDE with identifier PXD043702.

Comparative transcriptome analysis reveals deep molecular landscapes in stony coral Montipora clade

Introduction: Coral reefs, among the most invaluable ecosystems in the world, face escalating threats from climate change and anthropogenic activities. To decipher the genetic underpinnings of coral adaptation and resilience, we undertook comprehensive transcriptome profiling of two emblematic coral species, Montipora foliosa and Montipora capricornis, leveraging PacBio Iso-Seq technology. These species were strategically selected for their ecological significance and their taxonomic proximity within the Anthozoa class. Methods: Our study encompassed the generation of pristine transcriptomes, followed by thorough functional annotation via diverse databases. Subsequently, we quantified transcript abundance and scrutinized gene expression patterns, revealing notable distinctions between the two species. Results: Intriguingly, shared orthologous genes were identified across a spectrum of coral species, highlighting a substantial genetic conservation within scleractinian corals. Importantly, a subset of genes, integral to biomineralization processes, emerged as exclusive to scleractinian corals, shedding light on their intricate evolutionary history. Furthermore, we discerned pronounced upregulation of genes linked to immunity, stress response, and oxidative-reduction processes in M. foliosa relative to M. capricornis. These findings hint at the presence of more robust mechanisms in M. foliosa for maintaining internal equilibrium and effectively navigating external challenges, underpinning its potential ecological advantage. Beyond elucidating genetic adaptation in corals, our research underscores the urgency of preserving genetic diversity within coral populations. Discussion: These insights hold promise for informed conservation strategies aimed at safeguarding these imperiled ecosystems, bearing ecological and economic significance. In synthesis, our study seamlessly integrates genomic inquiry with ecological relevance, bridging the gap between molecular insights and the imperative to conserve coral reefs in the face of mounting threats.

Morphological, functional, compositional and transcriptional constraints shape the distinct venom profiles of the assassin bug Sycanus croceovittatus

Predatory bugs employ a salivary venom apparatus to generate complex venoms for capturing and digesting prey. The venom apparatus consists of different glands for the production of distinct venom sets, but the underlying mechanisms behind this process remain poorly understood. Here we present a comprehensive analysis of the morphological, functional, compositional and transcriptional characteristics of venoms derived from posterior main gland (PMG), anterior main gland (AMG), and accessory gland (AG) of the assassin bug Sycanus croceovittatus. Structural observations revealed the intricate constructions of the venom apparatus, enabling the production and storage of three distinct venom sets in anatomically varied glands and allowing them to be modulated in a context-dependent manner upon utilization. There were remarkable differences in the biological activities exhibited by PMG, AMG, and AG venoms. Proteotranscriptomic analysis demonstrated that these venoms displayed compositional heterogeneity at both the quantity and variety levels of proteins. Transcriptional profiles of the identified venom proteins revealed gland-specific or biased expression patterns. These findings indicate that the divergence in venom profiles among different glands arises from morphological, functional, compositional and transcriptional constraints on the venom apparatus, reflecting remarkable morphogenesis and regulatory gene networks responsible for the compartmentalized production of venom proteins in different glands.

Comprehensive transcripts analysis based on single-molecule real-time sequencing and Illumina sequencing provides insights into the mining of Toll-like receptor family in Schizothorax lissolabiatus

Schizothorax lissolabiatus is an economically important cold-water fish species in southwestern China. Because of water pollution and habitat destruction, the number of wild populations has dramatically decreased. In this study, we used PacBio single-molecule real-time (SMRT) sequencing and Illumina sequencing to generate the first full-length transcriptome and transcriptome, respectively. A total of 19 310 polished consensus reads (PC) were obtained, with an average length of 1379 bp and an N50 length of 1485 bp. Meanwhile, 12 253 transcripts were successfully annotated as known homologous genes. The pathway annotation indicated that the enrichment and expression of most genes were mainly related to membrane, signal transduction and binding, and immune response. Furthermore, we identified 16 Toll-like receptors (TLRs) by mining the data from the transcripts. Phylogeny analysis showed that S. lissolabiatus TLR genes (slTLRs) supported the classification of TLRs into six families as in other vertebrates. Selection pressure analyses showed that 16 slTLRs revealed purification selection at the overall evolutionary selection. Further, positive selection signals were still detected in eight slTLRs, and most of the positive selection sites were located in the leucine-rich repeat region (LRR domain) associated with the recognition of pathogenic microorganisms, indicating that the function of these slTLR genes may be affected. Tissue specific expression analysis showed all slTLRs are present in kidney, spleen and liver but the relative expression varied among tissues. In conclusion, this study not only provided a valuable resource of transcripts for further research on S. lissolabiatus, but also contributed to improve the current understanding of the evolutionary history of immune-related genes and the TLR gene family in S. lissolabiatus.

Integrating Full-Length Transcriptome and RNA Sequencing of Siberian Wildrye (Elymus sibiricus) to Reveal Molecular Mechanisms in Response to Drought Stress

Drought is one of the most significant limiting factors affecting plant growth and development on the Qinghai-Tibet Plateau (QTP). Mining the drought-tolerant genes of the endemic perennial grass of the QTP, Siberian wildrye (Elymus sibiricus), is of great significance to creating new drought-resistant varieties which can be used in the development of grassland livestock and restoring natural grassland projects in the QTP. To investigate the transcriptomic responsiveness of E. sibiricus to drought stress, PEG-induced short- and long-term drought stress was applied to two Siberian wildrye genotypes (drought-tolerant and drought-sensitive accessions), followed by third- and second-generation transcriptome sequencing analysis. A total of 40,708 isoforms were detected, of which 10,659 differentially expressed genes (DEGs) were common to both genotypes. There were 2107 and 2498 unique DEGs in the drought-tolerant and drought-sensitive genotypes, respectively. Additionally, 2798 and 1850 DEGs were identified in the drought-tolerant genotype only under short- and long-term conditions, respectively. DEGs numbering 1641 and 1330 were identified in the drought-sensitive genotype only under short- and long-term conditions, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that all the DEGs responding to drought stress in E. sibiricus were mainly associated with the mitogen-activated protein kinase (MAKP) signaling pathway, plant hormone signal transduction, the linoleic acid metabolism pathway, the ribosome pathway, and plant circadian rhythms. In addition, Nitrate transporter 1/Peptide transporter family protein 3.1 (NPF3.1) and Auxin/Indole-3-Acetic Acid (Aux/IAA) family protein 31(IAA31) also played an important role in helping E. sibiricus resist drought. This study used transcriptomics to investigate how E. sibiricus responds to drought stress, and may provide genetic resources and references for research into the molecular mechanisms of drought resistance in native perennial grasses and for breeding drought-tolerant varieties.

Molecular Mechanism Underlying the Sorghum sudanense (Piper) Stapf. Response to Osmotic Stress Determined via Single-Molecule Real-Time Sequencing and Next-Generation Sequencing

Drought, as a widespread environmental factor in nature, has become one of the most critical factors restricting the yield of forage grass. Sudangrass (Sorghum sudanense (Piper) Stapf.), as a tall and large grass, has a large biomass and is widely used as forage and biofuel. However, its growth and development are limited by drought stress. To obtain novel insight into the molecular mechanisms underlying the drought response and excavate drought tolerance genes in sudangrass, the first full-length transcriptome database of sudangrass under drought stress at different time points was constructed by combining single-molecule real-time sequencing (SMRT) and next-generation transcriptome sequencing (NGS). A total of 32.3 Gb of raw data was obtained, including 20,199 full-length transcripts with an average length of 1628 bp after assembly and correction. In total, 11,921 and 8559 up- and down-regulated differentially expressed genes were identified between the control group and plants subjected to drought stress. Additionally, 951 transcription factors belonging to 50 families and 358 alternative splicing events were found. A KEGG analysis of 158 core genes exhibiting continuous changes over time revealed that 'galactose metabolism' is a hub pathway and raffinose synthase 2 and β-fructofuranosidase are key genes in the response to drought stress. This study revealed the molecular mechanism underlying drought tolerance in sudangrass. Furthermore, the genes identified in this study provide valuable resources for further research into the response to drought stress.

SMRT Sequencing Technology Was Used to Construct the Batocera horsfieldi (Hope) Transcriptome and Reveal Its Features

Batocera horsfieldi (Hope) (Coleoptera: Cerambycidae) is an important forest pest in China that mainly infests timber and economic forests. This pest primarily causes plant tissue to necrotize, rot, and eventually die by feeding on the woody parts of tree trunks. To gain a deeper understanding of the genetic mechanism of B. horsfieldi, this study employed single-molecule real-time sequencing (SMRT) and Illumina RNA-seq technologies to conduct full-length transcriptome sequencing of the insect. Total RNA extracted from male and female adults was mixed and subjected to SMRT sequencing, generating a complete transcriptome. Transcriptome analysis, prediction of long non-coding RNA (lncRNA), coding sequences (CDs), analysis of simple sequence repeats (SSR), prediction of transcription factors, and functional annotation of transcripts were performed in this study. The collective 20,356,793 subreads (38.26 G, clean reads) were generated, including 432,091 circular consensus sequences and 395,851 full-length non-chimera reads. The full-length non-chimera reads (FLNC) were clustered and redundancies were removed, resulting in 39,912 consensus reads. SSR and ANGEL software v3.0 were used for predicting SSR and CDs. In addition, four tools were used for annotating 6058 lncRNAs, identifying 636 transcription factors. Furthermore, a total of 84,650 transcripts were functionally annotated in seven different databases. This is the first time that the full-length transcriptome of B. horsfieldi has been obtained using SMRT sequencing. This provides an important foundation for investigating the gene regulation underlying the interaction between B. horsfieldi and its host plants through gene editing in the future and provides a scientific basis for the prevention and control of B. horsfieldi.

Full-Length Transcriptome of the Great Himalayan Leaf-Nosed Bats (Hipposideros armiger) Optimized Genome Annotation and Revealed the Expression of Novel Genes

The Great Himalayan Leaf-nosed bat (Hipposideros armiger) is one of the most representative species of all echolocating bats and is an ideal model for studying the echolocation system of bats. An incomplete reference genome and limited availability of full-length cDNAs have hindered the identification of alternatively spliced transcripts, which slowed down related basic studies on bats' echolocation and evolution. In this study, we analyzed five organs from H. armiger for the first time using PacBio single-molecule real-time sequencing (SMRT). There were 120 GB of subreads generated, including 1,472,058 full-length non-chimeric (FLNC) sequences. A total of 34,611 alternative splicing (AS) events and 66,010 Alternative Polyadenylation (APA) sites were detected by transcriptome structural analysis. Moreover, a total of 110,611 isoforms were identified, consisting of 52% new isoforms of known genes and 5% of novel gene loci, as well as 2112 novel genes that have not been annotated before in the current reference genome of H. armiger. Furthermore, several key novel genes, including Pol, RAS, NFKB1, and CAMK4, were identified as being associated with nervous, signal transduction, and immune system processes, which may be involved in regulating the auditory nervous perception and immune system that helps bats to regulate in echolocation. In conclusion, the full-length transcriptome results optimized and replenished existing H. armiger genome annotation in multiple ways and offer advantages for newly discovered or previously unrecognized protein-coding genes and isoforms, which can be used as a reference resource.

Integrated metabolome, transcriptome analysis, and multi-flux full-length sequencing offer novel insights into the function of lignin biosynthesis as a Sesuvium portulacastrum response to salt stress

Sesuvium portulacastrum is a typical halophyte. However, few studies have investigated its salt-tolerant molecular mechanism. In this study, metabolome, transcriptome, and multi-flux full-length sequencing analysis were conducted to investigate the significantly different metabolites (SDMs) and differentially expressed genes (DEGs) of S. portulacastrum samples under salinity. The complete-length transcriptome of S. portulacastrum was developed, which contained 39,659 non-redundant unigenes. RNA-seq results showed that 52 DEGs involved in lignin biosynthesis may be responsible for S. portulacastrum salt tolerance. Furthermore, 130 SDMs were identified, and the salt response could be attributed to the p-coumaryl alcohol-rich in lignin biosynthesis. The co-expression network that was constructed after comparing the different salt treatment processes showed that the p-Coumaryl alcohol was linked to 30 DEGs. Herein, 8 structures genes, i.e., Sp4CL, SpCAD, SpCCR, SpCOMT, SpF5H, SpCYP73A, SpCCoAOMT, and SpC3'H were identified as significant factors in regulating lignin biosynthesis. Further investigation revealed that 64 putative transcription factors (TFs) may interact with the promoters of the above-mentioned genes. Together, the data revealed a potential regulatory network comprising important genes, putative TFs, and metabolites involved in the lignin biosynthesis of S. portulacastrum roots under salt stress, which could serve as a rich useful genetic resource for breeding excellent salt-tolerant plants.

Mechanisms underlying the formation of complex color patterns on Nigella orientalis (Ranunculaceae) petals

Complex color patterns on petals are widespread in flowering plants, yet the mechanisms underlying their formation remain largely unclear. Here, by conducting detailed morphological, anatomical, biochemical, optical, transcriptomic, and functional studies, we investigated the cellular bases, chromogenic substances, reflectance spectra, developmental processes, and underlying mechanisms of complex color pattern formation on Nigella orientalis petals. We found that the complexity of the N. orientalis petals in color pattern is reflected at multiple levels, with the amount and arrangement of different pigmented cells being the key. We also found that biosynthesis of the chromogenic substances of different colors is sequential, so that one color/pattern is superimposed on another. Expression and functional studies further revealed that a pair of R2R3-MYB genes function cooperatively to specify the formation of the eyebrow-like horizontal stripe and the Mohawk haircut-like splatters. Specifically, while NiorMYB113-1 functions to draw a large splatter region, NiorMYB113-2 functions to suppress the production of anthocyanins from the region where a gap will form, thereby forming the highly specialized pattern. Our results provide a detailed portrait for the spatiotemporal dynamics of the coloration of N. orientalis petals and help better understand the mechanisms underlying complex color pattern formation in plants.

Abscisic acid mediated strawberry receptacle ripening involves the interplay of multiple phytohormone signaling networks

As a canonical non-climacteric fruit, strawberry (Fragaria spp.) ripening is mainly mediated by abscisic acid (ABA), which involves multiple other phytohormone signalings. Many details of these complex associations are not well understood. We present an coexpression network, involving ABA and other phytohormone signalings, based on weighted gene coexpression network analysis of spatiotemporally resolved transcriptome data and phenotypic changes of strawberry receptacles during development and following various treatments. This coexpression network consists of 18,998 transcripts and includes transcripts related to phytohormone signaling pathways, MADS and NAC family transcription factors and biosynthetic pathways associated with fruit quality. Members of eight phytohormone signaling pathways are predicted to participate in ripening and fruit quality attributes mediated by ABA, of which 43 transcripts were screened to consist of the hub phytohormone signalings. In addition to using several genes reported from previous studies to verify the reliability and accuracy of this network, we explored the role of two hub signalings, small auxin up-regulated RNA 1 and 2 in receptacle ripening mediated by ABA, which are also predicted to contribute to fruit quality. These results and publicly accessible datasets provide a valuable resource to elucidate ripening and quality formation mediated by ABA and involves multiple other phytohormone signalings in strawberry receptacle and serve as a model for other non-climacteric fruits.

Transcriptome profiling reveals the roles of pigment formation mechanisms in yellow Paeonia delavayi flowers

The yellow colour of ornamental varieties of tree peony originated from Paeonia delavayi. However, but P. delavayi and Paeonia suffruticosa belong to different subgroups, so hybridization is difficult and results in a long breeding cycle. However, no comprehensive transcriptomic profiling has focused on the colour formation mechanisms of yellow tree peony petals. Analysing the colour formation mechanism of yellow petals in P. delavayi is very important for directional molecular breeding. In this study, the transcriptional map of yellow pigment development in petals was used to analyse the mechanism of petal colour formation. We analysed the genes related to the metabolism of flavonoids and carotenoids and the transcription factors (TFs) involved in P. delavayi var. lutea (pure yellow individual) yellow pigment development using transcriptome sequence profiling. Transcriptome sequence profiles revealed three and four differentially expressed transcripts (DETs) involved in flavonoid biosynthesis and carotenoid biosynthesis, respectively. An analysis of DETs in the flavonoid pathway showed that chalcone synthase (CHS) and chalcone 2´-glucosyltransferases (THC2'GT) act in synergy to synthesize isosalipurposide (ISP). CHS and flavonol synthase (FLS) synergistically synthesize quercetin and kaempferol. DEG analysis of the carotenoid pathway revealed that phytoene synthase (PSY), carotenoid isomerase (CRTISO) and β-carotene hydroxylases (CHYB) play a key role in regulating lutein formation, and carotenoid cleavage dioxygenase (CCD) plays an important role in the degradation of carotenoids. These two pathways may be regulated by TF families such as bHLH, ARF, and MYB. The results of the transient overexpression of genes showed that CHS and CHI are regulated by PdMYB2. In this study, the molecular mechanism of ISP synthesis was analysed in depth, and the complete metabolic pathway of carotenoids in Paeonia L. was reported for the first time. By studying the formation mechanism of yellow pigment in P. delavayi petals, a breeding strategy for improving flavonol and carotenoid contents and reducing anthocyanin synthesis by genetic engineering was suggested.

Full-Length Transcriptome of Myotis pilosus as a Reference Resource and Mining of Auditory and Immune Related Genes

Rickett's big-footed bat, Myotis pilosus, which belongs to the family Vespertilionida, is the only known piscivorous bat in East Asia. Accurate whole genome and transcriptome annotations are essential for the study of bat biological evolution. The lack of a whole genome for M. pilosus has limited our understanding of the molecular mechanisms underlying the species' evolution, echolocation, and immune response. In the present work, we sequenced the entire transcriptome using error-corrected PacBio single-molecule real-time (SMRT) data. Then, a total of 40 GB of subreads were generated, including 29,991 full-length non-chimeric (FLNC) sequences. After correction by Illumina short reads and de-redundancy, we obtained 26,717 error-corrected isoforms with an average length of 3018.91 bp and an N50 length of 3447 bp. A total of 1528 alternative splicing (AS) events were detected by transcriptome structural analysis. Furthermore, 1032 putative transcription factors (TFs) were identified, with additional identification of several long non-coding RNAs (lncRNAs) with high confidence. Moreover, several key genes, including PRL-2, DPP4, Glul, and ND1 were also identified as being associated with metabolism, immunity, nervous system processes, and auditory perception. A multitude of pattern recognition receptors was identified, including NLR, RLR, SRCR, the antiviral molecule IRF3, and the IFN receptor subunit IFNAR1. High-quality reference genomes at the transcriptome level may be used to quantify gene or transcript expression, evaluate alternative splicing levels, identify novel transcripts, and enhance genome annotation in bats.

The gastrodin biosynthetic pathway in Pholidota chinensis Lindl. revealed by transcriptome and metabolome profiling

Pholidota chinensis Lindl. is an epiphytic or lithophytic perennial herb of Orchidaceae family used as a garden flower or medicinal plant to treat high blood pressure, dizziness and headache in traditional Chinese medicine. Gastrodin (GAS) is considered as a main bioactive ingredient of this herb but the biosynthetic pathway remains unclear in P. chinensis. To elucidate the GAS biosynthesis and identify the related genes in P. chinensis, a comprehensive analysis of transcriptome and metabolome of roots, rhizomes, pseudobulbs and leaves were performed by using PacBio SMART, Illumina Hiseq and Ultra Performance Liquid Chromatography Tandem Mass Spectrometry (UPLC-MS/MS). A total of 1,156 metabolites were identified by UPLC-MS/MS, of which 345 differential metabolites were mainly enriched in phenylpropanoid/phenylalanine, flavone and flavonol biosynthesis. The pseudobulbs make up nearly half of the fresh weight of the whole plant, and the GAS content in the pseudobulbs was also the highest in four tissues. Up to 23,105 Unigenes were obtained and 22,029 transcripts were annotated in the transcriptome analysis. Compared to roots, 7,787, 8,376 and 9,146 differentially expressed genes (DEGs) were identified in rhizomes, pseudobulbs and leaves, respectively. And in total, 80 Unigenes encoding eight key enzymes for GAS biosynthesis, were identified. Particularly, glycosyltransferase, the key enzyme of the last step in the GAS biosynthetic pathway had 39 Unigenes candidates, of which, transcript28360/f2p0/1592, was putatively identified as the most likely candidate based on analysis of co-expression, phylogenetic analysis, and homologous searching. The metabolomics and transcriptomics of pseudobulbs versus roots showed that 8,376 DEGs and 345 DEMs had a substantial association based on the Pearson's correlation. This study notably enriched the metabolomic and transcriptomic data of P. chinensis, and it provides valuable information for GAS biosynthesis in the plant.

Full-length transcriptome analysis of the bloom-forming dinoflagellate Akashiwo sanguinea by single-molecule real-time sequencing

The dinoflagellate Akashiwo sanguinea is a harmful algal species and commonly observed in estuarine and coastal waters around the world. Harmful algal blooms (HABs) caused by this species lead to serious environmental impacts in the coastal waters of China since 1998 followed by huge economic losses. However, the full-length transcriptome information of A. sanguinea is still not fully explored, which hampers basic genetic and functional studies. Herein, single-molecule real-time (SMRT) sequencing technology was performed to characterize the full-length transcript in A. sanguinea. Totally, 83.03 Gb SMRT sequencing clean reads were generated, 983,960 circular consensus sequences (CCS) with average lengths of 3,061 bp were obtained, and 81.71% (804,016) of CCS were full-length non-chimeric reads (FLNC). Furthermore, 26,461 contigs were obtained after being corrected with Illumina library sequencing, with 20,037 (75.72%) successfully annotated in the five public databases. A total of 13,441 long non-coding RNA (lncRNA) transcripts, 3,137 alternative splicing (AS) events, 514 putative transcription factors (TFs) members from 23 TF families, and 4,397 simple sequence repeats (SSRs) were predicted, respectively. Our findings provided a sizable insights into gene sequence characteristics of A. sanguinea, which can be used as a reference sequence resource for A. sanguinea draft genome annotation, and will contribute to further molecular biology research on this harmful bloom algae.

Full-Length Transcriptome Sequencing Provides Insights into Flavonoid Biosynthesis in Camellia nitidissima Petals

Flavonoids are the main medicinal ingredients in Camellia nitidissima, but the regulatory mechanism of flavonoid biosynthesis in flowers is unclear; therefore, the flavonoids in C. nitidissima have not been effectively used. The present study performed full-length transcriptome sequencing of C. nitidissima flower. Furthermore, the reported RNA-sequencing data of C. nitidissima petals were reanalyzed using the full-length transcriptome as a reference, and the regulatory mechanism of flavonoid synthesis in petals was elucidated. The analysis identified 43,350 isoforms annotated in non-redundant protein (Nr), Kyoto Encyclopedia of Genes and Genomes (KEGG), EuKaryotic Orthologous Groups (KOG), and Swiss-Prot databases, among which 34,602 aligned to Camellia sinensis genes. A total of 11,857 differentially expressed genes (DEGs), including 112 related to flavonoid synthesis, were identified by pairwise comparison. Subsequently, analysis of the phylogeny and the conserved motifs of R2R3-MYB using the proteins sequences identified three R2R3-MYB transcription factors that regulated flavonoid biosynthesis. Weighted gene co-expression network analysis (WGCNA) identified phenylalanine ammonia-lyase (PAL) and 4-coumarate: CoA ligase(4CL) as the hub genes and showed that bHLH79 interacted with PAL. Finally, validated the expression of seven DEGs involved in flavonoid biosynthesis using real-time quantitative PCR (qRT-PCR). Thus, the present study generated and used the full-length transcriptome as the reference to analyze the transcriptome of petals and proposed a possible regulatory mechanism of flavonoid synthesis in C. nitidissima. The study's findings unravel the genetic mechanisms underlying flavonoid synthesis and suggest candidate genes for the genetic improvement of C. nitidissima.

Full-Length Transcriptome Maps of Reef-Building Coral Illuminate the Molecular Basis of Calcification, Symbiosis, and Circa-Dian Genes

Coral transcriptomic data largely rely on short-read sequencing, which severely limits the understanding of coral molecular mechanisms and leaves many important biological questions unresolved. Here, we sequence the full-length transcriptomes of four common and frequently dominant reef-building corals using the PacBio Sequel II platform. We obtain information on reported gene functions, structures, and expression profiles. Among them, a comparative analysis of biomineralization-related genes provides insights into the molecular basis of coral skeletal density. The gene expression profiles of the symbiont Symbiodiniaceae are also isolated and annotated from the holobiont sequence data. Finally, a phylogenetic analysis of key circadian clock genes among 40 evolutionarily representative species indicates that there are four key members in early metazoans, including cry genes; Clock or Npas2; cyc or Arntl; and tim, while per, as the fifth member, occurs in Bilateria. In summary, this work provides a foundation for further work on the manipulation of skeleton production or symbiosis to promote the survival of these important organisms.

Comparative transcriptomics of Pinus massoniana organs provides insights on terpene biosynthesis regulation

Terpenoids are the most important natural products collected from conifer species. However, the molecular mechanisms and core factors underlying terpenoid biosynthesis in Pinus massoniana remain unclear. To clarify these mechanisms, this study aimed to identify potential genes that might participate in the terpenoid biosynthesis of P. massoniana. In this study, single molecule real-time (SMRT) sequencing and expression analysis were used to confirm the expression patterns of genes involved in the cones, immature needles, mature needles, immature branches, and mature branches of P. massoniana. A total of 31,331 lncRNAs and 71,240 mRNAs were identified from these organs, and the greatest number of differentially expressed genes (DEGs) was discovered between needles and branches. Weighted gene coexpression network analysis (WGCNA) classified all expressed genes into nine typical modules with 11 kinds of transcription factors (TFs), namely, AP2-ERF, ARF, AUX-IAA, C2H2, Dof, F-box, SBP, WRKY, bHLH, bZIP, and GRAS, and seven kinds of functional genes, namely, ABC transporter, cellulose synthase (CesA), leucine-rich repeats (LRR), cytochrome P450 (CYT P450), pathogenesis-related protein (PR), terpene synthase (TPS), and chlorophyllase enzyme. A molecular network was constructed for hub genes, TFs, and functional genes in three modules. The potential function of eight candidate genes, including PmbHLH2, PmERF1, PmRGA, PmGAI, PmbZIP1, PmLOB1, PmMADS1, and PmMYB1, was validated through correlation analysis between terpenoid contents and expression levels, subcellular localization, and transcriptional activation activity, which provides us with probable regulators of terpenoid biosynthesis in conifers.

Characterization of the gene expression profile response to drought stress in Haloxylon using PacBio single-molecule real-time and Illumina sequencing

Haloxylon ammodendron and Haloxylon persicum are important drought-tolerant plants in northwest China. The whole-genome sequencing of H. ammodendron and H. persicum grown in their natural environment is incomplete, and their transcriptional regulatory network in response to drought environment remains unclear. To reveal the transcriptional responses of H. ammodendron and H. persicum to an arid environment, we performed single-molecule real-time (SMRT) and Illumina RNA sequencing. In total, 20,246,576 and 908,053 subreads and 435,938 and 210,334 circular consensus sequencing (CCS) reads were identified by SMRT sequencing of H. ammodendron and H. persicum, and 15,238 and 10,135 unigenes, respectively, were successfully obtained. In addition, 9,794 and 7,330 simple sequence repeats (SSRs) and 838 and 71 long non-coding RNAs were identified. In an arid environment, the growth of H. ammodendron was restricted; plant height decreased significantly; basal and branch diameters became thinner and hydrogen peroxide (H₂O₂) content and peroxidase (POD) activity were increased. Under dry and wet conditions, 11,803 and 15,217 differentially expressed genes (DEGs) were identified in H. ammodendron and H. persicum, respectively. There were 319 and 415 DEGs in the signal transduction pathways related to drought stress signal perception and transmission, including the Ca²⁺ signal pathway, the ABA signal pathway, and the MAPK signal cascade. In addition, 217 transcription factors (TFs) and 398 TFs of H. ammodendron and H. persicum were differentially expressed, including FAR1, MYB, and AP2/ERF. Bioinformatic analysis showed that under drought stress, the expression patterns of genes related to active oxygen [reactive oxygen species (ROS)] scavenging, functional proteins, lignin biosynthesis, and glucose metabolism pathways were altered. Thisis the first full-length transcriptome report concerning the responses of H. ammodendron and H. persicum to drought stress. The results provide a foundation for further study of the adaptation to drought stress. The full-length transcriptome can be used in genetic engineering research.

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.

SMRT Sequencing Reveals Candidate Genes and Pathways With Medicinal Value in Cipangopaludina chinensis

Cipangopaludina chinensis is an economically important aquatic snail with high medicinal value. However, molecular biology research on C. chinensis is limited by the lack of a reference genome, so the analysis of its transcripts is an important step to study the regulatory genes of various substances in C. chinensis. Herein, we conducted the first full-length transcriptome analysis of C. chinensis using PacBio single-molecule real-time (SMRT) sequencing technology. We identified a total of 26,312 unigenes with an average length of 2,572 bp, of which the largest number of zf-c2h2 transcription factor families (120,18.24%) were found, and also observed that the majority of the 8,058 SSRs contained 4-7 repeat units, which provided data for subsequent work on snail genetics Subsequently, 91.86% (24,169) of the genes were successfully annotated to the four major databases, while the highest homology was observed with Pomacea canaliculata. Functional annotation revealed that the majority of transcripts were enriched in metabolism, signal transduction and Immune-related pathways, and several candidate genes involved in drug metabolism and immune response were identified (e.g., CYP1A1, CYP2J, CYP2U1, GST, ,PIK3, PDE3A, PRKAG). This study lays a foundation for future molecular biology research and provides a reference for studying genes associated with the medicinal value of C. chinensis.

Full-Length Transcriptomics Reveals Complex Molecular Mechanism of Salt Tolerance in Bromus inermis L

Bromus inermis L. (commonly known as smooth bromegrass) is a grass species with high nutritional value, great palatability, cold tolerance, and grazing resistance, which has been widely cultivated for pasture and sand fixation in northern and northwestern China. Salt stress is a main environmental factor limiting growth and production of smooth bromegrass. In this study, we performed PacBio Iso-Seq to construct the first full-length transcriptome database for smooth bromegrass under 300 mM NaCl treatment at different time points. Third-generation full-length transcriptome sequencing yielded 19.67 G polymerase read bases, which were assembled into 355,836 full-length transcripts with an average length of 2,542 bp. A total of 116,578 differentially expressed genes were obtained by comparing the results of third-generation sequencing and second-generation sequencing. GO and KEGG enrichment analyses revealed that multiple pathways were differently activated in leaves and roots. In particular, a number of genes participating in the molecular network of plant signal perception, signal transduction, transcription regulation, antioxidant defense, and ion regulation were affected by NaCl treatment. In particular, the CBL-CIPK, MAPK, ABA signaling network, and SOS core regulatory pathways of Ca²⁺ signal transduction were activated to regulate salt stress response. In addition, the expression patterns of 10 salt-responsive genes were validated by quantitative real-time PCR, which were consistent with those detected by RNA-Seq. Our results reveal the molecular regulation of smooth bromegrass in response to salt stress, which are important for further investigation of critical salt responsive genes and molecular breeding of salt-tolerant smooth bromegrass.

Detection of Structural Variations and Fusion Genes in Breast Cancer Samples Using Third-Generation Sequencing

Background: Structural variations (SVs) are common genetic alterations in the human genome that could cause different phenotypes and diseases, including cancer. However, the detection of structural variations using the second-generation sequencing was limited by its short read length, which restrained our understanding of structural variations.

Methods: In this study, we developed a 28-gene panel for long-read sequencing and employed it to Oxford Nanopore Technologies and Pacific Biosciences platforms. We analyzed structural variations in the 28 breast cancer-related genes through long-read genomic and transcriptomic sequencing of tumor, para-tumor, and blood samples in 19 breast cancer patients.

Results: Our results showed that some somatic SVs were recurring among the selected genes, though the majority of them occurred in the non-exonic region. We found evidence supporting the existence of hotspot regions for SVs, which extended our previous understanding that they exist only for single nucleotide variations.

Conclusion: In conclusion, we employed long-read genomic and transcriptomic sequencing to identify SVs from breast cancer patients and proved that this approach holds great potential in clinical application.

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.

Identification and characterization of the key enzyme in the biosynthesis of the neurotoxin β-ODAP in grass pea

Grass pea (Lathyrus sativus L.) is a grain legume commonly grown in Asia and Africa for food and forage. It is a highly nutritious and robust crop, capable of surviving both droughts and floods. However, it produces a neurotoxic compound, β-N-oxalyl-L-α,β-diaminopropionic acid (β-ODAP), which can cause a severe neurological disorder when consumed as a primary diet component. While the catalytic activity associated with β-ODAP formation was demonstrated more than 50 years ago, the enzyme responsible for this activity has not been identified. Here, we report on the identity, activity, 3D structure, and phylogenesis of this enzyme—β-ODAP synthase (BOS). We show that BOS belongs to the benzylalcohol O-acetyltransferase, anthocyanin O-hydroxycinnamoyltransferase, anthranilate N-hydroxycinnamoyl/benzoyltransferase, deacetylvindoline 4-O-acetyltransferase superfamily of acyltransferases and is structurally similar to hydroxycinnamoyl transferase. Using molecular docking, we propose a mechanism for its catalytic activity, and using heterologous expression in tobacco leaves (Nicotiana benthamiana), we demonstrate that expression of BOS in the presence of its substrates is sufficient for β-ODAP production in vivo. The identification of BOS may pave the way toward engineering β-ODAP–free grass pea cultivars, which are safe for human and animal consumption.

Identification of essential genes involved in metabolism-based resistance mechanism to fenoxaprop-P-ethyl in Polypogon fugax

BACKGROUND

Metabolic resistance is a worldwide concern for weed control but has not yet been well-characterized at the genetic level. Previously, we have identified an Asia minor bluegrass (Polypogon fugax Nees ex Steud.) population AHHY exhibiting cytochrome P450 (P450)-involved metabolic resistance to fenoxaprop-P-ethyl. In this study, we aimed to confirm the metabolic fenoxaprop-P-ethyl resistance in AHHY and uncover the potential herbicide metabolism-related genes in this economically damaging weed species.

RESULTS

Liquid chromatography-tandem mass spectrometry (LC-MS/MS) assays indicated the metabolic rates of fenoxaprop-P-ethyl were significantly faster in resistant (R, AHHY) than in susceptible (S, SDTS) plants. The amount of phytotoxic fenoxaprop-P peaked at 12 h after herbicide treatment (HAT) and started to decrease at 24 HAT in both biotypes. R and S plants at 24 HAT were sampled to conduct isoform-sequencing (Iso-Seq) and RNA-sequencing (RNA-Seq). A reference transcriptome containing 24 972 full-length isoforms was obtained, of which 24 329 unigenes were successfully annotated. Transcriptomic profiling identified 28 detoxifying enzyme genes constitutively and/or herbicide-induced up-regulated in R than in S plants. Real-time quantitative polymerase chain reaction (RT-qPCR) confirmed 17 genes were consistently up-regulated in R and its F1 generation plants. They were selected as potential fenoxaprop-P-ethyl metabolism-related genes, including ten P450s, one glutathione-S-transferase, one UDP-glucosyltransferase, and five adenosine triphosphate (ATP)-binding cassette transporters.

CONCLUSION

This study revealed that the enhanced rates of fenoxaprop-P-ethyl metabolism in P. fugax were very likely driven by the herbicide metabolism-related genes. The transcriptome data generated by Iso-Seq combined with RNA-Seq will provide abundant gene resources for understanding the molecular mechanisms of resistance in P. fugax.

First Multi-Organ Full-Length Transcriptome of Tree Fern Alsophila spinulosa Highlights the Stress-Resistant and Light-Adapted Genes

Alsophila spinulosa, a relict tree fern, is a valuable plant for investigating environmental adaptations. Its genetic resources, however, are scarce. We used the PacBio and Illumina platforms to sequence the polyadenylated RNA of A. spinulosa root, rachis, and pinna, yielding 125,758, 89,107, and 89,332 unigenes, respectively. Combining the unigenes from three organs yielded a non-redundant reference transcriptome with 278,357 unigenes and N50 of 4141 bp, which were further reconstructed into 38,470 UniTransModels. According to functional annotation, pentatricopeptide repeat genes and retrotransposon-encoded polyprotein genes are the most abundant unigenes. Clean reads mapping to the full-length transcriptome is used to assess the expression of unigenes. The stress-induced ASR genes are highly expressed in all three organs, which is validated by qRT-PCR. The organ-specific upregulated genes are enriched for pathways involved in stress response, secondary metabolites, and photosynthesis. Genes for five types of photoreceptors, CRY signaling pathway, ABA biosynthesis and transduction pathway, and stomatal movement-related ion channel/transporter are profiled using the high-quality unigenes. The gene expression pattern coincides with the previously identified stomatal characteristics of fern. This study is the first multi-organ full-length transcriptome report of a tree fern species, the abundant genetic resources and comprehensive analysis of A. spinulosa, which provides the groundwork for future tree fern research.

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.

Binding Pattern Reconstructions of FGF-FGFR Budding-Inducing Signaling in Reef-Building Corals

Reef-building corals play an important role in marine ecosystems. However, owing to climate change, ocean acidification, and predation by invasive crown-of-thorns starfish, these corals are declining. As marine animals comprise polyps, reproduction by asexual budding is pivotal in scleractinian coral growth. The fibroblast growth factor (FGF) signaling pathway is essential in coral budding morphogenesis. Here, we sequenced the full-length transcriptomes of four common and frequently dominant reef-building corals and screened out the budding-related FGF and FGFR genes. Thereafter, three-dimensional (3D) models of FGF and FGFR proteins as well as FGF-FGFR binding models were reconstructed. Based on our findings, the FGF8-FGFR3 binding models in Pocillopora damicornis, Montipora capricornis, and Acropora muricata are typical receptor tyrosine kinase-signaling pathways that are similar to the Kringelchen (FGFR) in hydra. However, in P. verrucosa, FGF8 is not the FGFR3 ligand, which is found in other hydrozoan animals, and its FGFR3 must be activated by other tyrosine kinase-type ligands. Overall, this study provides background on the potentially budding propagation signaling pathway activated by the applications of biological agents in reef-building coral culture that could aid in the future restoration of coral reefs.

Full-length transcriptome sequencing of Heliocidaris crassispina using PacBio single-molecule real-time sequencing

The lack of high-throughput sequencing data makes the research progress of Heliocidaris crassispina slow. Therefore, we used PacBio single-molecule real-time sequencing to generate the first full-length transcriptome. Here, 31,181 isoforms were obtained, with an average length of 2383.20 and a N50 length of 2732 bp. Meanwhile, 764 alternative splicing (AS) events, 5098 long-noncoding RNAs (LncRNAs), 6978 simple sequence repeats (SSRs), and 950 hypothetical transcript factors (TFs) were identified. Moreover, five key innate immune pattern recognition receptors (PRRs), including toll-like receptor (TLR), NACHT domain and leucine-rich repeat (NLR), scavenger receptor cysteine-rich (SRCR), peptidoglycan recognition proteins (PGRP), and gram-negative binding proteins (GNBP), were searched in the transcriptome. In addition, 37 isoforms enriched in KEGG and GO immune systems were also detected. The study provid abundant data support for the current research on H. crassispina.

Full-Length Transcriptome: A Reliable Alternative for Single-Cell RNA-Seq Analysis in the Spleen of Teleost Without Reference Genome

Fish is considered as a supreme model for clarifying the evolution and regulatory mechanism of vertebrate immunity. However, the knowledge of distinct immune cell populations in fish is still limited, and further development of techniques advancing the identification of fish immune cell populations and their functions are required. Single cell RNA-seq (scRNA-seq) has provided a new approach for effective in-depth identification and characterization of cell subpopulations. Current approaches for scRNA-seq data analysis usually rely on comparison with a reference genome and hence are not suited for samples without any reference genome, which is currently very common in fish research. Here, we present an alternative, i.e. scRNA-seq data analysis with a full-length transcriptome as a reference, and evaluate this approach on samples from Epinephelus coioides-a teleost without any published genome. We show that it reconstructs well most of the present transcripts in the scRNA-seq data achieving a sensitivity equivalent to approaches relying on genome alignments of related species. Based on cell heterogeneity and known markers, we characterized four cell types: T cells, B cells, monocytes/macrophages (Mo/MΦ) and NCC (non-specific cytotoxic cells). Further analysis indicated the presence of two subsets of Mo/MΦ including M1 and M2 type, as well as four subsets in B cells, i.e. mature B cells, immature B cells, pre B cells and early-pre B cells. Our research will provide new clues for understanding biological characteristics, development and function of immune cell populations of teleost. Furthermore, our approach provides a reliable alternative for scRNA-seq data analysis in teleost for which no reference genome is currently available.

Full-length transcriptome assembly of Andrias davidianus (Amphibia: Caudata) skin via hybrid sequencing

The Chinese giant salamander, Andrias davidianus, is the largest amphibian species in the world; it is thus an economically and ecologically important species. The skin of A. davidianus exhibits complex adaptive structural and functional adaptations to facilitate survival in aquatic and terrestrial ecosystems. Here, we report the first full-length amphibian transcriptome from the dorsal skin of A. davidianus, which was assembled using hybrid sequencing and the PacBio and Illumina platforms. A total of 153,038 transcripts were hybrid assembled (mean length of 2039 bp and N50 of 2172 bp), and 133,794 were annotated in at least one database (nr, Swiss-Prot, KEGG, KOGs, GO, and nt). A total of 58,732, 68,742, and 115,876 transcripts were classified into 24 KOG categories, 1903 GO term categories, and 46 KEGG pathways (level 2), respectively. A total of 207,627 protein-coding regions, 785 transcription factors, 27,237 potential long non-coding RNAs, and 8299 simple sequence repeats were also identified. The hybrid-assembled transcriptome recovered more full-length transcripts, had a higher N50 contig length, and a higher annotation rate of unique genes compared with that assembled in previous studies using next-generation sequencing. The high-quality full-length reference gene set generated in this study will help elucidate the genetic characteristics of A. davidianus skin and aid the identification of functional skin proteins.

Full-Length Transcriptome of Thalassiosira weissflogii as a Reference Resource and Mining of Chitin-Related Genes

β-Chitin produced by diatoms is expected to have significant economic and ecological value due to its structure, which consists of parallel chains of chitin, its properties and the high abundance of diatoms. Nevertheless, few studies have functionally characterised chitin-related genes in diatoms owing to the lack of omics-based information. In this study, we first compared the chitin content of three representative Thalassiosira species. Cell wall glycosidic linkage analysis and chitin/chitosan staining assays showed that Thalassiosira weissflogii was an appropriate candidate chitin producer. A full-length (FL) transcriptome of T. weissflogii was obtained via PacBio sequencing. In total, the FL transcriptome comprised 23,362 annotated unigenes, 710 long non-coding RNAs (lncRNAs), 363 transcription factors (TFs), 3113 alternative splicing (AS) events and 3295 simple sequence repeats (SSRs). More specifically, 234 genes related to chitin metabolism were identified and the complete biosynthetic pathways of chitin and chitosan were explored. The information presented here will facilitate T. weissflogii molecular research and the exploitation of β-chitin-derived high-value enzymes and products.