Search

Amaryllidaceae alkaloids (AAs) are complex plant secondary metabolites possessing a wide range of biological activities. 4′-O-methylnorbelladine (4OMN) is the branchpoint intermediate for the entire AAs, and was the last common intermediate before AA pathway branches diverge. The cyclization of 4OMN by C-C oxidative coupling, which can afford para-para ′, ortho-para ′, and para-ortho ′ scaffold, was catalyzed by cytochrome P450 96T (CYP96T) family enzymes. To clarify the mechanisms involved in this controversial step, four CYP96T homologs (LauCYP96T1, LauCYP96T1-like-1, LauCYP96T1-like-2 and LauCYP96T1-like-3) were cloned from the full-length transcriptome of Lycoris aurea. All the four LauCYP96T are localized to endoplasmic reticulum. Functional analysis reveals that LauCYP96T1 and LauCYP96T1-like proteins display inverted regioselectivity for oxidative coupling of 4OMN, in which LauCYP96T1 and LauCYP96T1-like-2 dominantly afford para-para ′ scaffold, and LauCYP96T1-like-1 and LauCYP96T1-like-3 are responsible for para-ortho ′ scaffold formation. Using molecular homology modeling and docking studies, we predicted models for the binding of 4OMN to LauCYP96T, and identified two amino acid residues that might be responsible for the dominant changes in generated products of para-ortho ′ and para-para ′ oxidative coupling. Our results highlight the functional diversity and promiscuity of LauCYP96T enzymes and might provide valuable information for Amaryllidaceae alkaloid production. [ABSTRACT FROM AUTHOR]

Glenea cantor Fabricius (Cerambycidae: Lamiinae) is an important pest that damages kapok trees in Southeast Asia with a wide adaptability to temperature. Glycerol is a protectant and energy source for insects in low-temperature environments. However, glycerol biosynthesis-related genes at the molecular level are limited in G. cantor. In this study, the supercooling points and freezing points at different stages were measured, and the cold hardiness of male and female pupae significantly differed. Moreover, a full-length transcriptome of G. cantor was established; glycerol kinase (GK) and glycerol-3-phosphate dehydrogenase (GPDH) genes, which are related to glycerol metabolism, were identified, with a special focus on their expression profiles. A total of 24,476 isoforms stemmed from the full-length transcriptome, along with 568 lncRNAs, 56 transcription factor (TF) families, and 1467 alternative splicing (AS) events. The KEGG pathway enrichment analysis revealed that the isoforms associated with AS were enriched primarily in glycerolipid and glycerophospholipid metabolism. In total, three GK genes and one GPDH gene were identified, and GcGK1 and GcGK3 presented differential sex expression during the pupal stage, which may play a role in thermal adaptability. This study provides a valuable transcriptional database of G. cantor and helps to elucidate the function of glycerol in the thermal adaptation mechanism of longhorn beetles. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Amaryllidaceae alkaloids (AAs) are complex plant secondary metabolites possessing a wide range of biological activities. 4'-O-methylnorbelladine (4OMN) is the branchpoint intermediate for the entire AAs, and was the last common intermediate before AA pathway branches diverge. The cyclization of 4OMN by C-C oxidative coupling, which can afford para-para', ortho-para', and para-ortho' scaffold, was catalyzed by cytochrome P450 96T (CYP96T) family enzymes. To clarify the mechanisms involved in this controversial step, four CYP96T homologs (LauCYP96T1, LauCYP96T1-like-1, LauCYP96T1-like-2 and LauCYP96T1-like-3) were cloned from the full-length transcriptome of Lycoris aurea. All the four LauCYP96T are localized to endoplasmic reticulum. Functional analysis reveals that LauCYP96T1 and LauCYP96T1-like proteins display inverted regioselectivity for oxidative coupling of 4OMN, in which LauCYP96T1 and LauCYP96T1-like-2 dominantly afford para-para' scaffold, and LauCYP96T1-like-1 and LauCYP96T1-like-3 are responsible for para-ortho' scaffold formation. Using molecular homology modeling and docking studies, we predicted models for the binding of 4OMN to LauCYP96T, and identified two amino acid residues that might be responsible for the dominant changes in generated products of para-ortho' and para-para' oxidative coupling. Our results highlight the functional diversity and promiscuity of LauCYP96T enzymes and might provide valuable information for Amaryllidaceae alkaloid production. [ABSTRACT FROM AUTHOR]

The Asian honey bee (Apis cerana) plays a crucial role in providing abundant bee products and in maintaining ecological balance. Despite the availability of the genomic sequence of the Asian honey bee, its transcriptomic information remains largely incomplete. To address this issue, here we constructed three pooled RNA samples from the queen, drone, and worker bees of A. cerana and performed full-length RNA sequencing using Nanopore single-molecule sequencing technology. Ultimately, we obtained 160,811 full-length transcript sequences from 19,859 genes, with 141,189 being novel transcripts, of which 130,367 were functionally annotated. We detected 520, 324, and 1823 specifically expressed transcripts in the queen, worker, and drone bees, respectively. Furthermore, we identified 38,799 alternative splicing (AS) events from 5710 genes, 44,243 alternative polyadenylation (APA) sites from 1649 gene loci, 88,187 simple sequence repeats (SSRs), and 17,387 long noncoding RNAs (lncRNAs). Leveraging these transcripts as references, we identified 6672, 7795, and 6804 differentially expressed transcripts (DETs) in comparisons of queen ovaries vs drone testes, worker ovaries vs drone testes, and worker ovaries vs queen ovaries, respectively. Our research results provide a comprehensive set of reference transcript datasets for Apis cerana, offering important sequence information for further exploration of its gene functions. [ABSTRACT FROM AUTHOR]

Simple Summary: A total of 37 mitochondrial genes appears in a certain arrangement in most insects, and the transcription of this conserved mitochondrial genome (mitogenome) has been well-studied in species like Drosophila melanogaster, Erthesina fullo and Coridius chinensis. However, transcription of the mitogenome with extensive gene rearrangement has rarely been studied. In this research, we sequenced the mitogenome and mitochondrial transcriptome of Sericothrips houjii (Thysanoptera: Thripidae: Sericothripinae). The mitogenome of S. houjii exhibited extensive gene rearrangement and contained two control regions (CRs) with sequence repeats. Compared to the insect mitogenome with a typical gene order, the mitogenomic feature, relative gene expression level, transcriptional model and post-transcriptional cleavage of S. houjii were quite different. Unlike other insects where ribosomal RNA (rRNA) is typically highly expressed, ND4/ND4L in S. houjii exhibits the highest expression. Both strands of this mitogenome were entirely transcribed and the bicistronic messenger RNA (mRNA) COI/ND3 was reported for the first time in insects. Our study provides new insights into the transcriptional and post-transcriptional regulation processes in the insect mitogenome with extensive gene rearrangement and duplicated CRs. The mitochondrial genome (mitogenome) of Thysanoptera has extensive gene rearrangement, and some species have repeatable control regions. To investigate the characteristics of the gene expression, transcription and post-transcriptional processes in such extensively gene-rearranged mitogenomes, we sequenced the mitogenome and mitochondrial transcriptome of Sericothrips houjii to analyze. The mitogenome was 14,965 bp in length and included two CRs contains 140 bp repeats between COIII-trnN (CR1) and trnT-trnP (CR2). Unlike the putative ancestral arrangement of insects, S. houjii exhibited only six conserved gene blocks encompassing 14 genes (trnL2-COII, trnD-trnK, ND2-trnW, ATP8-ATP6, ND5-trnH-ND4-ND4L and trnV-lrRNA). A quantitative transcription map showed the gene with the highest relative expression in the mitogenome was ND4-ND4L. Based on analyses of polycistronic transcripts, non-coding RNAs (ncRNAs) and antisense transcripts, we proposed a transcriptional model of this mitogenome. Both CRs contained the transcription initiation sites (TISs) and transcription termination sites (TTSs) of both strands, and an additional TIS for the majority strand (J-strand) was found within antisense lrRNA. The post-transcriptional cleavage processes followed the "tRNA punctuation" model. After the cleavage of transfer RNAs (tRNAs), COI and ND3 matured as bicistronic mRNA COI/ND3 due to the translocation of intervening tRNAs, and the 3′ untranslated region (UTR) remained in the mRNAs for COII, COIII, CYTB and ND5. Additionally, isoform RNAs of ND2, srRNA and lrRNA were identified. In summary, the relative mitochondrial gene expression levels, transcriptional model and post-transcriptional cleavage process of S. houjii are notably different from those insects with typical mitochondrial gene arrangements. In addition, the phylogenetic tree of Thripidae including S. houjii was reconstructed. Our study provides insights into the phylogenetic status of Sericothripinae and the transcriptional and post-transcriptional regulation processes of extensively gene-rearranged insect mitogenomes. [ABSTRACT FROM AUTHOR]

Vairimorpha (Nosema) ceranae is a single-cellular fungus that obligately infects the midgut epithelial cells of adult honeybees, causing bee microsporidiosis and jeopardizing bee health and production. This work aims to construct the full-length transcriptome of V. ceranae and conduct a relevant investigation using PacBio single-molecule real-time (SMRT) sequencing technology. Following PacBio SMRT sequencing, 41,950 circular consensus (CCS) were generated, and 25,068 full-length non-chimeric (FLNC) reads were then detected. After polishing, 4387 high-quality, full-length transcripts were gained. There are 778, 2083, 1202, 1559, 1457, 1232, 1702, and 3896 full-length transcripts that could be annotated to COG, GO, KEGG, KOG, Pfam, Swiss-Prot, eggNOG, and Nr databases, respectively. Additionally, 11 alternative splicing (AS) events occurred in 6 genes were identified, including 1 alternative 5′ splice-site and 10 intron retention. The structures of 225 annotated genes in the V. ceranae reference genome were optimized, of which 29 genes were extended at both 5′ UTR and 3′ UTR, while 90 and 106 genes were, respectively, extended at the 5′ UTR as well as 3′ UTR. Furthermore, a total of 29 high-confidence lncRNAs were obtained, including 12 sense-lncRNAs, 10 lincRNAs, and 7 antisense-lncRNAs. Taken together, the high-quality, full-length transcriptome of V. ceranae was constructed and annotated, the structures of annotated genes in the V. ceranae reference genome were improved, and abundant new genes, transcripts, and lncRNAs were discovered. Findings from this current work offer a valuable resource and a crucial foundation for molecular and omics research on V. ceranae. [ABSTRACT FROM AUTHOR]

Next-generation sequencing has revolutionized cancer genomics by enabling high-throughput mutation screening yet detecting fusion genes reliably remains challenging. Long-read sequencing offers potential for accurate fusion transcript identification, though challenges persist. In this study, we present an optimized workflow using nanopore sequencing technology to precisely identify fusion transcripts. Our approach encompasses a tailored library preparation protocol, data processing, and fusion gene analysis pipeline. We evaluated the performance using Universal Human Reference RNA and human adenocarcinoma cell lines. Our optimized nanopore sequencing workflow generated high-quality full-length transcriptome data characterized by an extended length distribution and comprehensive transcript coverage. Validation experiments confirmed novel fusion events with potential clinical relevance. Our protocol aims to mitigate biases and enhance accuracy, facilitating increased adoption in clinical diagnostics. Continued advancements in long-read sequencing promise deeper insights into fusion gene biology and improved cancer diagnostics.

Flavonoid glycosides are widespread in plants, and are of great interest owing to their diverse biological activities and effectiveness in preventing chronic diseases. Periploca forrestii, a renowned medicinal plant of the Apocynaceae family, contains diverse flavonoid glycosides and is clinically used to treat rheumatoid arthritis and traumatic injuries. However, the mechanisms underlying the biosynthesis of these flavonoid glycosides have not yet been elucidated. In this study, we used widely targeted metabolomics and full-length transcriptome sequencing to identify flavonoid diversity and biosynthetic genes in P. forrestii. A total of 120 flavonoid glycosides, including 21 C-, 96 O-, and 3 C/O-glycosides, were identified and annotated. Based on 24,123 full-length coding sequences, 99 uridine diphosphate sugar-utilizing glycosyltransferases (UGTs) were identified and classified into 14 groups. Biochemical assays revealed that four UGTs exhibited O-glycosyltransferase activity toward apigenin and luteolin. Among them, PfUGT74B4 and PfUGT92A8 were highly promiscuous and exhibited multisite O-glycosylation or consecutive glycosylation activities toward various flavonoid aglycones. These four glycosyltransferases may significantly contribute to the diversity of flavonoid glycosides in P. forrestii. Our findings provide a valuable genetic resource for further studies on P. forrestii and insights into the metabolic engineering of bioactive flavonoid glycosides.

Abstract Background Daphnia galeata is a suitable model organism for investigating predator-induced defense. Genes and pathways exhibiting differential expression between fish kairomone-treated and untreated groups in D. galeata have been identified. However, understanding of the significance of alternative splicing, a crucial process of the regulation of gene expression in eukaryotes, to this mechanism remains limited. This study measured life-history traits and conducted short-read RNA sequencing and long-read isoform sequencing of two Korean D. galeata genotypes (KB1 and KE2) to uncover the genetic mechanism underlying their phenotypic plasticity under predation stress. Results KB1 exhibited strategies to enhance fertility and decrease body length when exposed to fish kairomones, while KE2 deployed an adaptive strategy to increase body length. Full-length transcriptomes from KB1 and KE2 yielded 65,736 and 57,437 transcripts, respectively, of which 32 differentially expressed transcripts (DETs) were shared under predation stress across both genotypes. Prominent DETs common to both genotypes were related to energy metabolism and the immune system. Additionally, differential alternative splicing (DAS) events were detected in both genotypes in response to fish kairomones. DAS genes shared between both genotypes may indicate their significant role in the post-transcriptional stress response to fish predation. Calpain-3, involved in digestion and nutrient absorption, was identified as a DAS gene in both genotypes when exposed to fish kairomones. In addition, the gene encoding thymosin beta, which is related to growth, was found to be a statistically significant DAS only in KB1, while that encoding ultraspiracle protein, also associated with growth, was only identified in KE2. Moreover, transcripts encoding proteins such as EGF-like domain-containing protein, vitellogenin fused with superoxide dismutase, and others were identified overlapping between DAS events and DETs and potentially elucidating their association with the observed phenotypic variation in each genotype. Conclusions Our findings highlight the crucial role of alternative splicing in modulating transcriptome landscape under predation stress in D. galeata, emphasizing the requirement for integrating gene expression and splicing analyses in evolutionary adaptation studies.

Auxin response factor (ARF) is a transcription factor family that mediates auxin signaling and regulates various biological processes. To investigate the ARF gene family members and their roles in response to high temperature and drought stress, the phylogenetic and expression patterns of the ARF gene family members in Huperzia serrata were analyzed using full-length transcriptome and RNA-seq data. Based on bioinformatics analysis, the physicochemical properties, domains, conserved motifs, phylogeny, tissue expression patterns, and expression patterns of the ARF gene family under high temperature and drought stress were analyzed. The results were as follows:(1) A total of 24 ARF family members were screened in the full-length transcriptome of H. serrata, all of which were acidic proteins and hydrophilic proteins. (2) Subcellular localization prediction revealed that 24 HsARF were all localized in nucleus. (3) Phylogenetic analysis revealed that HsARF had a distant genetic relationship with angiosperms Arabidopsis and rice, and only share two common ARF ancestors with higher flowering plants. (4) Domain analysis revealed that, except for HsARF18/23/24, most HsARF had B3 domains. Analysis of secondary structure found that the highest proportion of HsARF protein was random curling, followed by elongated chains and α-helix. Three-dimensional (3D) protein structure prediction model used in the 24 HsARF proteins was only four. (5) RNA-seq analysis showed that the expression levels of 7 HsARF were high in all detected organs. The expression levels of 10 HsARF in stems were higher than those in roots and leaves. Otherwise, the expression levels of HsARF13 and HsARF14 in leaves were lower than those in roots and stems. (6) The expression levels of HsARF undergo significant changes under high temperature and drought stress, with 18 HsARF being induced to varying degrees by high temperature stress, 7 HsARF responded to drought stress, with 3 HsARF induced by drought, while 4 HsARF inhibited by drought. The results provide theoretical references for the functional and biological breeding study of ARF gene family in H. serrata.

In order to clarify the molecular bioinformatics characteristics and tissue expression patterns of the COBRA gene family members of Huperzia serrata, the physicochemical properties, domains, conserved motifs, cis-acting elements, and genes expression of the family members (HsCOBRAs) were analyzed by bioinformatics techniques, based on the full-length transcriptome data of the H. serrata. The results were as follows: (1) A total of 24 HsCOBRAs family members were screened in the full-length transcriptome of H. serrata, including 9 acidic proteins, 11 stabilizing proteins, 5 hydrophobic proteins, 7 proteins with transmembrane structures, and 3 proteins with signal peptides. (2) Subcellular localization was found in the cell wall, chloroplast, nucleus, and cell membrane. (3) Structural analysis revealed that HsCOBRAs had 7 domains and 6 conserved motifs, and partial members had a highly conserved CCVS structure. (4) HsCOBRAs had 45 cis-acting elements such as CAAT-box and TATA-box. (5) HsCOBRA2 had the highest expressions in leaves, spores, stems and gemma. The study results can provide theoretical basis for further research and biological function verification of HsCOBRAs.

Abstract 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.

Abstract 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.

[Objective] The preliminary screening of bHLH genes related to C. reticulata flower color by bioinformatics method provides support for further study of C. reticulata flower color. [Method] Based on the full-length transcriptome data of C. reticulata, the members of the bHLH gene family were identified by bioinformatics method, and their physical and chemical properties, secondary structure, phylogenetic relationships and domains were analyzed. The second-generation transcriptome and metabolome data were used to mine the bHLH genes related to the flower color of C. reticulata. Real-time quantitative polymerase chain reaction (RT-qPCR) was used to verify the expression of bHLH gene during flowering of C. reticulata among different varieties. [Result] Total 71 CrbHLHs genes were identified in the full-length transcriptome of C. reticulata. Their molecular weight ranged from 22 994.58 to 102 996.96 Da, and their equipotential ranged from 4.98 to 9.51. All of them were hydrophilic proteins. The cluster analysis with Arabidopsis thaliana showed that the bHLH gene family members of C. reticulata could be distributed into 14 subfamilies. Multiple omics analysis showed that CrbHLH8, CrbHLH61 and CrbHLH63 were positively correlated with cyanidin, while CrbHLH59 was negatively correlated with cyanidin. CrbHLH13 and CrbHLH71 were positively correlated with proanthocyanidins, and CrbHLH8 and CrbHLH61 were members of the IIIf subfamily. [Conclusion] The 71 CrbHLHs family members were identified from the full-length transcriptome of C. reticulata, of which 6 CrbHLHs were related to C. reticulata flower color. [ABSTRACT FROM AUTHOR]

Platycladus orientalis (P. orientalis) is a common tree used for vegetation restoration in northern China, and its large area propagation helps to improve site conditions. However, under harsh conditions such as poor land, the survival rate of P. orientalis is very low. Numerous studies have shown that root pruning can promote the formation of lateral roots in seedlings, enhancing the roots' capacity to absorb soil nutrients and water, and thereby improving the survival rate of seedlings. In this study, a one-third root pruning treatment was applied to P. orientalis seedlings, and the whole transcriptome of seedlings subjected to both control (CK) and root pruning treatments was sequenced to analyze their gene expression profiles. This study investigated the regulatory mechanisms of lateral root development in response to root pruning damage at the molecular level. Using nine cells, 15.28 Gb of clean data were obtained, which yielded 101,688 high-quality full-length transcript sequences and 22,955 low-quality full-length transcript sequences after clustering. Redundancy was then removed using CD-HIT, and Illumina RNA-seq sequencing produced 139.26 Gb of clean data. A total of 2025 differentially expressed genes (DEGs) were identified at three time points following root pruning treatment. Enrichment analysis revealed that the peroxidase gene family plays a significant role in lateral root proliferation. Furthermore, the expression levels of the peroxidase gene family were notably upregulated in comparison to the control group. Pathway enrichment analysis identified 22 relevant genes, which appeared to be highly associated with root growth and resilience to stress. Through examining the expression patterns and correlations of these genes, five central genes emerged as key players. The findings of this research suggest that the peroxidase gene family plays a crucial role in the stress response and root development of P. orientalis, providing reference and guidance for root development in other plant species. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Simple Summary: This study described a full-length transcriptome of Neocaridina denticulata sinensis that was generated using third-generation sequencing technology. The resulting transcriptome assembly comprised 5831 transcripts, revealing 90.5% novel isoforms of known genes. A functional annotation showed that 24.8% of the transcripts can be searched across seven key databases. Additionally, 1236 alternative splicing events, 344 transcription factors, and 124 long non-coding RNAs were predicted. Notably, a RING finger protein NHL-1 coding gene was characterized with 15 transcripts, and a phylogenetic analysis underscored its evolutionary connection to NHL-1 in other crustaceans. This full-length transcriptome dataset constitutes a valuable genetic resource in Neocaridina denticulata sinensis and will facilitate functional genomic studies and investigations into environmental adaptations in this species. 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. [ABSTRACT FROM AUTHOR]

Skeletal muscle grows in response to a combination of genetic and environmental factors, and its growth and development influence the quality of pork. Elucidating the molecular mechanisms regulating the growth and development of skeletal muscle is of great significance to both animal husbandry and farm management. The Jiangquan black pig is an excellent pig breed based on the original Yimeng black pig, importing the genes of the Duroc pig for meat traits, and cultivated through years of scientific selection and breeding. In this study, full-length transcriptome sequencing was performed on three growth stages of Jiangquan black pigs, aiming to study the developmental changes in Jiangquan black pigs at different developmental stages at the molecular level and to screen the key genes affecting the growth of skeletal muscle in Jiangquan black pigs. We performed an enrichment analysis of genes showing differential expression and constructed a protein–protein interaction network with the aim of identifying core genes involved in the development of Jiangquan black pigs. Notably, genes such as TNNI2, TMOD4, PLDIM3, MYOZ1, and MYH1 may be potential regulators of muscle development in Jiangquan black pigs. Our results contribute to the understanding of the molecular mechanisms of skeletal muscle development in this pig breed, which will facilitate molecular breeding efforts and the development of pig breeds to meet the needs of the livestock industry. [ABSTRACT FROM AUTHOR]

In order to clarify the molecular bioinformatics characteristics and tissue expression patterns of the COBRA gene family members of Huperzia serrata, the physicochemical properties, domains, conserved motifs, cis-acting elements, and genes expression of the family members(HsCOBRAs)were analyzed by bioinformatics techniques, based on the full-length transcriptome data of the H. serrata. The results were as follows:(1)A total of 24 HsCOBRAs family members were screened in the full-length transcriptome of H. serrata, including 9 acidic proteins, 11 stabilizing proteins, 5 hydrophobic proteins, 7 proteins with transmembrane structures, and 3 proteins with signal peptides.(2) Subcellular localization was found in the cell wall, chloroplast, nucleus, and cell membrane.(3)Structural analysis revealed that HsCOBRAs had 7 domains and 6 conserved motifs, and partial members had a highly conserved CCVS structure.(4)HsCOBRAs had 45 cis-acting elements such as CAAT-box and TATA-box.(5)HsCOBRA2 had the highest expressions in leaves, spores, stems and gemma. The study results can provide theoretical basis for further research and biological function verification of HsCOBRAs. [ABSTRACT FROM AUTHOR]

Auxin response factor(ARF)is a transcription factor family that mediates auxin signaling and regulates various biological processes. To investigate the ARF gene family members and their roles in response to high temperature and drought stress, the phylogenetic and expression patterns of the ARF gene family members in Huperzia serrata were analyzed using full-length transcriptome and RNA-seq data. Based on bioinformatics analysis, the physicochemical properties, domains, conserved motifs, phylogeny, tissue expression patterns, and expression patterns of the ARF gene family under high temperature and drought stress were analyzed. The results were as follows:(1)A total of 24 ARF family members were screened in the full-length transcriptome of H. serrata, all of which were acidic proteins and hydrophilic proteins.(2)Subcellular localization prediction revealed that 24 HsARF were all localized in nucleus.(3)Phylogenetic analysis revealed that HsARF had a distant genetic relationship with angiosperms Arabidopsis and rice, and only share two common ARF ancestors with higher flowering plants.(4)Domain analysis revealed that, except for HsARF18/23/24, most HsARF had B3 domains. Analysis of secondary structure found that the highest proportion of HsARF protein was random curling, followed by elongated chains and α-helix. Three-dimensional(3D)protein structure prediction model used in the 24 HsARF proteins was only four.(5)RNA-seq analysis showed that the expression levels of 7 HsARF were high in all detected organs. The expression levels of 10 HsARF in stems were higher than those in roots and leaves. Otherwise, the expression levels of HsARF13 and HsARF14 in leaves were lower than those in roots and stems.(6)The expression levels of HsARF undergo significant changes under high temperature and drought stress, with 18 HsARF being induced to varying degrees by high temperature stress, 7 HsARF responded to drought stress, with 3 HsARF induced by drought, while 4 HsARF inhibited by drought. The results provide theoretical references for the functional and biological breeding study of ARF gene family in H. serrata. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Panax notoginseng is a traditional Chinese medicine rich in many pharmacological components. The root of the ‘Miaoxiang Sanqi No. 2’ is yellow or greenish yellow, while a novel cultivar-‘Wenyuan Ziqi No. 1’ shows purple root and is thought to have high medicinal value. Little information is available about the anthocyanin biosynthesis in P. notoginseng root. In this study, we compared the ‘Miaoxiang Sanqi No. 2’ and ‘Wenyuan Ziqi No. 1’ in morphological, transcriptional and metabolic levels. The results showed that purple rich in the periderm, rhizome and phloem around cambium of the ‘Wenyuan Ziqi No. 1’ root and cyanidin 3-O-galactoside was the main anthocyanin causing purple. Moreover, ‘Wenyuan Ziqi No. 1’ highly accumulated in 155 metabolites, including flavones, phenylpropanoids and lipids. Transcriptome data showed that phenylpropanoid biosynthesis pathway genes are highly expressed in ‘Wenyuan Ziqi No. 1’. Conjoint analysis showed that anthocyanin biosynthesis pathway substances were highly accumulated in ‘Wenyuan Ziqi No. 1’, and the expression level of structural genes involved in anthocyanin biosynthesis pathway was higher in ‘Wenyuan Ziqi No. 1’. Meanwhile, eight R2R3-MYB genes that might be involved in anthocyanin biosynthesis were identified. The comprehensive analysis of two cultivars provides new insights into the understanding of root coloration in P. notoginseng.

Abstract Tulips (Tulipa gesneriana) are one of the most widely cultivated bulbous plants with substantial ornamental value. However, the lack of well-documented reference genomes has limited the research progress and molecular breeding of tulips. In the present study, a full-length transcriptome of a commercial tulip cultivar was obtained using single-molecule long-read sequencing (PacBio Iso-Seq). In total, 244,357 full-length transcripts were identified, which had an average length of 2,044 bp and an N50 value of 3,861; 67,350 of these were annotated to databases. An inaugural integrated analysis of the transcriptome and phytohormone profiles during flower opening and petal senescence was performed using Illumina RNA-seq, coupled with Mfuzz (an R pakage, http://mfuzz.sysbiolab.eu ) and weighted gene coexpression network analysis (WGCNA). A total of 16 gene coexpression and six transcription factor (TF) modules were constructed. Additionally, 26 hormone analogs were comprehensively profiled. Finally, a prominently novel gene, Tulipa gesneriana Homeobox12-like (TgHB12-like), which encodes an homeodomain–leucine zipper (HD-zip) TF, was identified as a pivotal regulator of petal senescence. Overall, this work facilitates the identification of hormones and TFs in plants related to flower opening and senescence in tulips. It also provides an important and valuable genetic basis for further research in them.

Abstract Background The popularity of Muscovy ducks is attributed not only to their conformation traits but also to their slightly higher content of breast and leg meat, as well as their stronger-tasting meat compared to that of typical domestic ducks. However, there is a lack of comprehensive systematic research on the development of breast muscle in Muscovy ducks. In addition, since the number of skeletal muscle myofibers is established during the embryonic period, this study conducted a full-length transcriptome sequencing and microRNA sequencing of the breast muscle. Muscovy ducks at four developmental stages, namely Embryonic Day 21 (E21), Embryonic Day 27 (E27), Hatching Day (D0), and Post-hatching Day 7 (D7), were used to isolate total RNA for analysis. Results A total of 68,161 genes and 472 mature microRNAs were identified. In order to uncover deeper insights into the regulation of mRNA by miRNAs, we conducted an integration of the differentially expressed miRNAs (known as DEMs) with the differentially expressed genes (referred to as DEGs) across various developmental stages. This integration allowed us to make predictions regarding the interactions between miRNAs and mRNA. Through this analysis, we identified a total of 274 DEGs that may serve as potential targets for the 68 DEMs. In the predicted miRNA‒mRNA interaction networks, let-7b, miR-133a-3p, miR-301a-3p, and miR-338-3p were the hub miRNAs. In addition, multiple DEMs also showed predicted target relationships with the DEGs associated with skeletal system development. These identified DEGs and DEMs as well as their predicted interaction networks involved in the regulation of energy homeostasis and muscle development were most likely to play critical roles in facilitating the embryo-to-hatchling transition. A candidate miRNA, miR-301a-3p, exhibited increased expression during the differentiation of satellite cells and was downregulated in the breast muscle tissues of Muscovy ducks at E21 compared to E27. A dual-luciferase reporter assay suggested that the ANKRD1 gene, which encodes a transcription factor, is a direct target of miR-301a-3p. Conclusions miR-301a-3p suppressed the posttranscriptional activity of ANKRD1, which is an activator of satellite cell proliferation, as determined with gain- and loss-of-function experiments. miR-301a-3p functions as an inducer of myogenesis by targeting the ANKRD1 gene in Muscovy ducks. These results provide novel insights into the early developmental process of black Muscovy breast muscles and will improve understanding of the underlying molecular mechanisms.

Abstract 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.

Background: The popularity of Muscovy ducks is attributed not only to their conformation traits but also to their slightly higher content of breast and leg meat, as well as their stronger-tasting meat compared to that of typical domestic ducks. However, there is a lack of comprehensive systematic research on the development of breast muscle in Muscovy ducks. In addition, since the number of skeletal muscle myofibers is established during the embryonic period, this study conducted a full-length transcriptome sequencing and microRNA sequencing of the breast muscle. Muscovy ducks at four developmental stages, namely Embryonic Day 21 (E21), Embryonic Day 27 (E27), Hatching Day (D0), and Post-hatching Day 7 (D7), were used to isolate total RNA for analysis. Results: A total of 68,161 genes and 472 mature microRNAs were identified. In order to uncover deeper insights into the regulation of mRNA by miRNAs, we conducted an integration of the differentially expressed miRNAs (known as DEMs) with the differentially expressed genes (referred to as DEGs) across various developmental stages. This integration allowed us to make predictions regarding the interactions between miRNAs and mRNA. Through this analysis, we identified a total of 274 DEGs that may serve as potential targets for the 68 DEMs. In the predicted miRNA‒mRNA interaction networks, let-7b, miR-133a-3p, miR-301a-3p, and miR-338-3p were the hub miRNAs. In addition, multiple DEMs also showed predicted target relationships with the DEGs associated with skeletal system development. These identified DEGs and DEMs as well as their predicted interaction networks involved in the regulation of energy homeostasis and muscle development were most likely to play critical roles in facilitating the embryo-to-hatchling transition. A candidate miRNA, miR-301a-3p, exhibited increased expression during the differentiation of satellite cells and was downregulated in the breast muscle tissues of Muscovy ducks at E21 compared to E27. A dual-luciferase reporter assay suggested that the ANKRD1 gene, which encodes a transcription factor, is a direct target of miR-301a-3p. Conclusions: miR-301a-3p suppressed the posttranscriptional activity of ANKRD1, which is an activator of satellite cell proliferation, as determined with gain- and loss-of-function experiments. miR-301a-3p functions as an inducer of myogenesis by targeting the ANKRD1 gene in Muscovy ducks. These results provide novel insights into the early developmental process of black Muscovy breast muscles and will improve understanding of the underlying molecular mechanisms. [ABSTRACT FROM AUTHOR]

Reticulitermes chinensis Snyder is an important pest in forestry and construction and is widely distributed in China. We found that Serratia marcescens Bizio strain SM1 has insecticidal activity to R. chinensis , but the pathogenic mechanism of SM1 to R. chinensis is not clear. Therefore, full-length transcriptome sequencing was performed on R. chinensis infected with SM1 and the control group. A total of 230 differentially expressed genes were identified by comparing SM1 infection group and the control group, among which 103 were downregulated and 127 were upregulated. We found downregulated genes in nine metabolic pathway categories, among which carbohydrate metabolism had the most downregulated genes, followed by energy metabolism and amino acid metabolism. We also found that some downregulated genes were related to pattern recognition receptors, cellular immunity, and humoral immunity, indicating that R. chinensis immunity was negatively affected by SM1 infection. In addition, some genes in signal transduction and genetic information processing pathways were downregulated. In this study, high-throughput full-length transcriptome analysis was used to analyse the pathogenic mechanism of SM1 to R. chinensi s. The results of this study provide useful information for exploring the relationship between SM1 and R. chinensis , and provide theoretical support for the future application of SM1 and the prevention and treatment of R. chinensis. [ABSTRACT FROM AUTHOR]

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.Key message: The clustering of chloroplast genes in Ophioglossum vulgatum demonstrates a tendency towards elevated expression levels in sporangium, leaf, and stem, while exhibiting a decline in evolutionary rates and selection pressure. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Introduction: Peanut (Arachis hypogaea L.), also called groundnut is an important oil and cash crop grown widely in the world. The annual global production of groundnuts has increased to approximately 50 million tons, which provides a rich source of vegetable oils and proteins for humans. Low temperature (non-freezing) is one of the major factors restricting peanut growth, yield, and geographic distribution. Since the complexity of cold-resistance trait, the molecular mechanism of cold tolerance and related gene networks were largely unknown in peanut. Methods: In this study, comparative transcriptomic analysis of two peanut cultivars (SLH vs. ZH12) with differential cold tolerance under low temperature (10°C) was performed using Oxford Nanopore Technology (ONT) platform. Results and discussion: As a result, we identified 8,949 novel gene loci and 95,291 new/novel isoforms compared with the reference database. More differentially expressed genes (DEGs) were discovered in cold-sensitive cultivar (ZH12) than cold-tolerant cultivar (SLH), while more alternative splicing events were found in SLH compared to ZH12. Gene Ontology (GO) analyses of the common DEGs showed that the “response to stress”, “chloroplast part”, and “transcription factor activity” were the most enriched GO terms, indicating that photosynthesis process and transcription factors play crucial roles in cold stress response in peanut. We also detected a total of 708 differential alternative splicing genes (DASGs) under cold stress compared to normal condition. Intron retention (IR) and exon skipping (ES) were the most prevalent alternative splicing (AS) events. In total, 4,993 transcription factors and 292 splicing factors were detected, many of them had differential expression levels and/or underwent AS events in response to cold stress. Overexpression of two candidate genes (encoding trehalose-6-phosphatephosphatases, AhTPPs) in yeast improves cold tolerance. This study not only provides valuable resources for the study of cold resistance in peanut but also lay a foundation for genetic modification of cold regulators to enhance stress tolerance in crops [ABSTRACT FROM AUTHOR]

Summary: The advent of full‐length transcriptome sequencing technologies has accelerated the discovery of novel splicing isoforms. However, existing alternative splicing (AS) tools are either tailored for short‐read RNA‐Seq data or designed for human and animal studies. The disparities in AS patterns between plants and animals still pose a challenge to the reliable identification and functional exploration of novel isoforms in plants.Here, we developed integrated full‐length alternative splicing analysis (iFLAS), a plant‐optimized AS toolkit that introduced a semi‐supervised machine learning method known as positive‐unlabeled (PU) learning to accurately identify novel isoforms. iFLAS also enables the investigation of AS functions from various perspectives, such as differential AS, poly(A) tail length, and allele‐specific AS (ASAS) analyses.By applying iFLAS to three full‐length transcriptome sequencing datasets, we systematically identified and functionally characterized maize (Zea mays) AS patterns. We found intron retention not only introduces premature termination codons, resulting in lower expression levels of isoforms, but may also regulate the length of 3′UTR and poly(A) tail, thereby affecting the functional differentiation of isoforms. Moreover, we observed distinct ASAS patterns in two genes within heterosis offspring, highlighting their potential value in breeding.These results underscore the broad applicability of iFLAS in plant full‐length transcriptome‐based AS research. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Polyphenols play a crucial role in fruit flavor. To elucidate the mechanism of fruit polyphenol metabolism, we constructed a transcriptome atlas through PacBio single-molecule real-time (SMRT) sequencing and Illumina next-generation sequencing (NGS) using Canarium album (Lour.) Raeusch., which is a fantastic fruit rich in polyphenolic compounds. In this work, PacBio full-length transcriptome assembly generated 135,439 isoforms with an average length of all isoforms of 2687.94 bp and an N50 length of 3224 bp. To gain deeper insights into the molecular mechanisms of polyphenol biosynthesis in C. album, we constructed twelve RNA-Seq libraries from four developmental stages of the fruits. We identified a total of 28,658 differentially expressed genes (DEGs). We found that many DEGs were involved in metabolic pathways, biosynthesis of secondary metabolites, biosynthesis of antibiotics, starch and sucrose metabolism, and plant hormone signal transduction. Here, we report the expression profiles of 215 DEGs encoding 27 enzymes involved in the polyphenol biosynthesis pathway in C. album. In addition, 285 differentially expressed transcription factors (TFs) continuously down-regulated in four developmental periods of C. album fruit, which may indicate their potential role in the response to polyphenol metabolism and phenylpropanoid biosynthesis pathways. This report will help us understand polyphenol biosynthesis's functions and metabolic mechanism in C. album. The transcriptome data provide a valuable resource for genetic and genomics research. They will facilitate future work exploiting C. album and other fruits used as medicine and food. [ABSTRACT FROM AUTHOR]

Acanthopanax gracilistylus is a deciduous plant in the family Araliaceae, which is commonly used in Chinese herbal medicine, as the root bark has functions of nourishing the liver and kidneys, removing dampness and expelling wind, and strengthening the bones and tendons. Kaurenoic acid (KA) is the main effective substance in the root bark of A. gracilistylus with strong anti-inflammatory effects. To elucidate the KA biosynthesis pathway, second-generation (DNA nanoball) and third-generation (Pacific Biosciences) sequencing were performed to analyze the transcriptomes of the A. gracilistylus leaves, roots, and stems. Among the total 505,880 isoforms, 408,954 were annotated by seven major databases. Sixty isoforms with complete open reading frames encoding 11 key enzymes involved in the KA biosynthesis pathway were identified. Correlation analysis between isoform expression and KA content identified a total of eight key genes. Six key enzyme genes involved in KA biosynthesis were validated by real-time quantitative polymerase chain reaction. Based on the sequence analysis, the spatial structure of ent-kaurene oxidase was modeled, which plays roles in the three continuous oxidations steps of KA biosynthesis. This study greatly enriches the transcriptome data of A. gracilistylus and facilitates further analysis of the function and regulation mechanism of key enzymes in the KA biosynthesis pathway. [ABSTRACT FROM AUTHOR]

Cyprinus carpio var. Quanzhounensis, native to Quanzhou County, Guilin City, Guangxi, has a dark brown body color, translucent gill cover, and abdominal skin, and is an important farmed species in the local integrated rice-fishery industry. A comparative study on the skin of Cyprinus carpio var. Quanzhounensis revealed a lack of reflective guanine crystals on the body surface and a significantly higher melanin content than that of C. carpio var. Jian, which was tentatively considered the direct cause of variations in body color of this group. Iridocytes are a pigment cell species that contain regularly arranged guanine crystals, which are the key material basis for the metallic luster of the fish body surface. In species such as medaka (Oryzias latipes) and zebrafish (Danio rerio), the absence of guanine crystals is considered a manifestation of the absence of iridocyte differentiation and therefore is ideal for studying the mechanism of pigment cell differentiation. The pnp4a, Gbx2, sox10, tfec genes and other iridocyte-related genes have been mined using mutant materials, and we have uncovered the differentiation mechanism that regulates the formation of iridescent cells. C. carpio var. color has abundant genetic variation in body color and is a good material for studying the mechanism of body color determination in fish. The roles of ASIP and MC1R in the aggregation and distribution of melanin and formation of black spots in C. carpio var. color were verified. The guanine crystalline deletion trait of C. carpio var. Quanzhounensis may be loaded with regulatory mechanism diversity and mutation loci related to guanine crystal formation or iridescent cell differentiation. In addition, as a rice-fish culture species, the living environment of C. carpio var. Quanzhounensis harvestmen differ significantly from pond and net-pen culture species, requiring high resistance to disease, adversity, and transport. It is unclear whether the absence of guanine crystals in their skin leads to changes in their basal physiological state, and in-depth studies are beneficial for accurate assessment of culture performance. To reveal the structural basis and transcriptomic characteristics of guanine crystalline deficiency traits in the skin of C. carpio var. Quanzhounensis in this study, we selected the subscale tissues with the most significant differences in guanine crystalline distribution as the control material, and transmission electron microscopy was used to observe the tissue structure and full-length transcriptome sequencing to understand the structural and transcriptomic characteristics of guanine crystalline deficiency in the skin of C. carpio var. Quanzhounensis. The results of this study provide information for the analysis of body color traits, identification of economic traits, and utilization of germplasm resources. Transmission electron microscopy of the subscale tissue sections revealed two significant differences in the histological structure of C. carpio var. Quanzhounensis, and C. carpio var. Jian. First, guanine crystals were absent in C. carpio var. Quanzhounensis, whereas guanine crystals were widely present in the tissues of C. carpio var. Jian and cascading cavities were observed in the sections after guanine crystals were dislodged. Second, the number and density of melanin particles in the tissues of C. carpio var. Quanzhounensis harvestmen were significantly higher than those of C. carpio var. Jian, showing smaller, darker, and more numerous particles, which is consistent with the darker color and lack of silvery reflective material on the body surface of C. carpio var. Quanzhounensis harvestmen. The transcriptome characteristics were analyzed using Oxford Nanopore (ONT) sequencing technology, and 2.88~3.26 Gb of high-quality data were obtained for each sample. The number of full-length sequences after filtering ribosomal RNA for all sample data was 2 203 826~2 412 500, and the proportion of full-length sequences for each sample was 87.06%~88.57%. The comparison rate was 90.35%~92.46%. Variable splicing events in the transcripts were counted; 3 075 variable splicing events and 57 624 variable polyadenylation events were detected; and 15 615 new coding region sequences and 771 long-stranded non-coding RNAs were predicted. The number of exon jumps and intron retention in variable splicing events differed significantly (P

Abstract Background Sphaeropteris brunoniana and Alsophila latebrosa are both old relict and rare tree ferns, which have experienced the constant changes of climate and environment. However, little is known about their high-quality genetic information and related research on environmental adaptation mechanisms of them. In this study, combined with PacBio and Illumina platforms, transcriptomic analysis was conducted on the roots, rachis, and pinna of S. brunoniana and A. latebrosa to identify genes and pathways involved in environmental adaptation. Additionally, based on the transcriptomic data of tree ferns, chloroplast genes were mined to analyze their gene expression levels and RNA editing events. Results In the study, we obtained 11,625, 14,391 and 10,099 unigenes of S. brunoniana root, rachis, and pinna, respectively. Similarly, a total of 13,028, 11,431 and 12,144 unigenes were obtained of A. latebrosa root, rachis, and pinna, respectively. According to the enrichment results of differentially expressed genes, a large number of differentially expressed genes were enriched in photosynthesis and secondary metabolic pathways of S. brunoniana and A. latebrosa. Based on gene annotation results and phenylpropanoid synthesis pathways, two lignin synthesis pathways (H-lignin and G-lignin) were characterized of S. brunoniana. Among secondary metabolic pathways of A. latebrosa, three types of WRKY transcription factors were identified. Additionally, based on transcriptome data obtained in this study, reported transcriptome data, and laboratory available transcriptome data, positive selection sites were identified from 18 chloroplast protein-coding genes of four tree ferns. Among them, RNA editing was found in positive selection sites of four tree ferns. RNA editing affected the protein secondary structure of the rbcL gene. Furthermore, the expression level of chloroplast genes indicated high expression of genes related to the chloroplast photosynthetic system in all four species. Conclusions Overall, this work provides a comprehensive transcriptome resource of S. brunoniana and A. latebrosa, laying the foundation for future tree fern research.

Abstract Background Corynespora leaf spot is a common leaf disease occurring in sesame, and the disease causes leaf yellowing and even shedding, which affects the growth quality of sesame. At present, the mechanism of sesame resistance to this disease is still unclear. Understanding the resistance mechanism of sesame to Corynespora leaf spot is highly important for the control of infection. In this study, the leaves of the sesame resistant variety (R) and the sesame susceptible variety (S) were collected at 0–48 hpi for transcriptome sequencing, and used a combined third-generation long-read and next-generation short-read technology approach to identify some key genes and main pathways related to resistance. Results The gene expression levels of the two sesame varieties were significantly different at 0, 6, 12, 24, 36 and 48 hpi, indicating that the up-regulation of differentially expressed genes in the R might enhanced the resistance. Moreover, combined with the phenotypic observations of sesame leaves inoculated at different time points, we found that 12 hpi was the key time point leading to the resistance difference between the two sesame varieties at the molecular level. The WGCNA identified two modules significantly associated with disease resistance, and screened out 10 key genes that were highly expressed in R but low expressed in S, which belonged to transcription factors (WRKY, AP2/ERF-ERF, and NAC types) and protein kinases (RLK-Pelle_DLSV, RLK-Pelle_SD-2b, and RLK-Pelle_WAK types). These genes could be the key response factors in the response of sesame to infection by Corynespora cassiicola. GO and KEGG enrichment analysis showed that specific modules could be enriched, which manifested as enrichment in biologically important pathways, such as plant signalling hormone transduction, plant-pathogen interaction, carbon metabolism, phenylpropanoid biosynthesis, glutathione metabolism, MAPK and other stress-related pathways. Conclusions This study provides an important resource of genes contributing to disease resistance and will deepen our understanding of the regulation of disease resistance, paving the way for further molecular breeding of sesame.

The Asian honey bee (Apis cerana) plays a crucial role in providing abundant bee products and in maintaining ecological balance. Despite the availability of the genomic sequence of the Asian honey bee, its transcriptomic information remains largely incomplete. To address this issue, here we constructed three pooled RNA samples from the queen, drone, and worker bees of A. cerana and performed full-length RNA sequencing using Nanopore single-molecule sequencing technology. Ultimately, we obtained 160,811 full-length transcript sequences from 19,859 genes, with 141,189 being novel transcripts, of which 130,367 were functionally annotated. We detected 520, 324, and 1823 specifically expressed transcripts in the queen, worker, and drone bees, respectively. Furthermore, we identified 38,799 alternative splicing (AS) events from 5710 genes, 44,243 alternative polyadenylation (APA) sites from 1649 gene loci, 88,187 simple sequence repeats (SSRs), and 17,387 long noncoding RNAs (lncRNAs). Leveraging these transcripts as references, we identified 6672, 7795, and 6804 differentially expressed transcripts (DETs) in comparisons of queen ovaries vs drone testes, worker ovaries vs drone testes, and worker ovaries vs queen ovaries, respectively. Our research results provide a comprehensive set of reference transcript datasets for Apis cerana, offering important sequence information for further exploration of its gene functions.

The mitochondrial genome (mitogenome) of Thysanoptera has extensive gene rearrangement, and some species have repeatable control regions. To investigate the characteristics of the gene expression, transcription and post-transcriptional processes in such extensively gene-rearranged mitogenomes, we sequenced the mitogenome and mitochondrial transcriptome of Sericothrips houjii to analyze. The mitogenome was 14,965 bp in length and included two CRs contains 140 bp repeats between COIII-trnN (CR1) and trnT-trnP (CR2). Unlike the putative ancestral arrangement of insects, S. houjii exhibited only six conserved gene blocks encompassing 14 genes (trnL2-COII, trnD-trnK, ND2-trnW, ATP8-ATP6, ND5-trnH-ND4-ND4L and trnV-lrRNA). A quantitative transcription map showed the gene with the highest relative expression in the mitogenome was ND4-ND4L. Based on analyses of polycistronic transcripts, non-coding RNAs (ncRNAs) and antisense transcripts, we proposed a transcriptional model of this mitogenome. Both CRs contained the transcription initiation sites (TISs) and transcription termination sites (TTSs) of both strands, and an additional TIS for the majority strand (J-strand) was found within antisense lrRNA. The post-transcriptional cleavage processes followed the “tRNA punctuation” model. After the cleavage of transfer RNAs (tRNAs), COI and ND3 matured as bicistronic mRNA COI/ND3 due to the translocation of intervening tRNAs, and the 3′ untranslated region (UTR) remained in the mRNAs for COII, COIII, CYTB and ND5. Additionally, isoform RNAs of ND2, srRNA and lrRNA were identified. In summary, the relative mitochondrial gene expression levels, transcriptional model and post-transcriptional cleavage process of S. houjii are notably different from those insects with typical mitochondrial gene arrangements. In addition, the phylogenetic tree of Thripidae including S. houjii was reconstructed. Our study provides insights into the phylogenetic status of Sericothripinae and the transcriptional and post-transcriptional regulation processes of extensively gene-rearranged insect mitogenomes.

Tetrastigma hemsleyanum is a valuable herb widely used in Chinese traditional and modern medicine. Winter cold severely limits the artificial cultivation of this plant, but the physiological and molecular mechanisms upon exposure to cold stress in T. hemsleyanum are unclear. T. hemsleyanum plants with different geographical origins exhibit large differences in response to cold stress. In this research study, using T. hemsleyanum ecotypes that exhibit frost tolerance (FR) and frost sensitivity (FS), we analyzed the response of cottage seedlings to a simulated frost treatment; plant hormones were induced with both short (2 h) and long (9 h) frost treatments, which were used to construct the full-length transcriptome and obtained 76,750 transcripts with all transcripts mapped to 28,805 genes, and 27,215 genes, respectively, annotated to databases. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed enrichment in plant hormone signaling pathways. Further analysis shows that differently expressed genes (DEGs) concentrated on calcium signaling, ABA biosynthesis and signal transduction, and ethylene in response to cold stress. We also found that endogenous ABA and ethylene content were increased after cold treatment, and exogenous ABA and ethylene significantly improved cold tolerance in both ecotypes. Our results elucidated the pivotal role of ABA and ethylene in response to cold stress in T. hemsleyanum and identified key genes. [ABSTRACT FROM AUTHOR]

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 noncoding 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. [ABSTRACT FROM AUTHOR]

Background: Sphaeropteris brunoniana and Alsophila latebrosa are both old relict and rare tree ferns, which have experienced the constant changes of climate and environment. However, little is known about their high-quality genetic information and related research on environmental adaptation mechanisms of them. In this study, combined with PacBio and Illumina platforms, transcriptomic analysis was conducted on the roots, rachis, and pinna of S. brunoniana and A. latebrosa to identify genes and pathways involved in environmental adaptation. Additionally, based on the transcriptomic data of tree ferns, chloroplast genes were mined to analyze their gene expression levels and RNA editing events. Results: In the study, we obtained 11,625, 14,391 and 10,099 unigenes of S. brunoniana root, rachis, and pinna, respectively. Similarly, a total of 13,028, 11,431 and 12,144 unigenes were obtained of A. latebrosa root, rachis, and pinna, respectively. According to the enrichment results of differentially expressed genes, a large number of differentially expressed genes were enriched in photosynthesis and secondary metabolic pathways of S. brunoniana and A. latebrosa. Based on gene annotation results and phenylpropanoid synthesis pathways, two lignin synthesis pathways (H-lignin and G-lignin) were characterized of S. brunoniana. Among secondary metabolic pathways of A. latebrosa, three types of WRKY transcription factors were identified. Additionally, based on transcriptome data obtained in this study, reported transcriptome data, and laboratory available transcriptome data, positive selection sites were identified from 18 chloroplast protein-coding genes of four tree ferns. Among them, RNA editing was found in positive selection sites of four tree ferns. RNA editing affected the protein secondary structure of the rbcL gene. Furthermore, the expression level of chloroplast genes indicated high expression of genes related to the chloroplast photosynthetic system in all four species. Conclusions: Overall, this work provides a comprehensive transcriptome resource of S. brunoniana and A. latebrosa, laying the foundation for future tree fern research. [ABSTRACT FROM AUTHOR]

Background: Corynespora leaf spot is a common leaf disease occurring in sesame, and the disease causes leaf yellowing and even shedding, which affects the growth quality of sesame. At present, the mechanism of sesame resistance to this disease is still unclear. Understanding the resistance mechanism of sesame to Corynespora leaf spot is highly important for the control of infection. In this study, the leaves of the sesame resistant variety (R) and the sesame susceptible variety (S) were collected at 0–48 hpi for transcriptome sequencing, and used a combined third-generation long-read and next-generation short-read technology approach to identify some key genes and main pathways related to resistance. Results: The gene expression levels of the two sesame varieties were significantly different at 0, 6, 12, 24, 36 and 48 hpi, indicating that the up-regulation of differentially expressed genes in the R might enhanced the resistance. Moreover, combined with the phenotypic observations of sesame leaves inoculated at different time points, we found that 12 hpi was the key time point leading to the resistance difference between the two sesame varieties at the molecular level. The WGCNA identified two modules significantly associated with disease resistance, and screened out 10 key genes that were highly expressed in R but low expressed in S, which belonged to transcription factors (WRKY, AP2/ERF-ERF, and NAC types) and protein kinases (RLK-Pelle_DLSV, RLK-Pelle_SD-2b, and RLK-Pelle_WAK types). These genes could be the key response factors in the response of sesame to infection by Corynespora cassiicola. GO and KEGG enrichment analysis showed that specific modules could be enriched, which manifested as enrichment in biologically important pathways, such as plant signalling hormone transduction, plant-pathogen interaction, carbon metabolism, phenylpropanoid biosynthesis, glutathione metabolism, MAPK and other stress-related pathways. Conclusions: This study provides an important resource of genes contributing to disease resistance and will deepen our understanding of the regulation of disease resistance, paving the way for further molecular breeding of sesame. [ABSTRACT FROM AUTHOR]

Simple Summary: The substantial accumulation of fat in the bodies of sheep will reduce the lean meat percentage to a certain extent, affect the carcass quality, and reduce the production efficiency. Therefore, there is a growing focus on investigating variations in adipose tissue formation across different regions, deciphering genetic regulatory mechanisms, reducing unnecessary fat hoarding, and assessing the utility of fat in current research. In this study, the expression of subcutaneous, tail, and visceral fat transcripts of Dorper × Hu sheep was quantified at the transcriptome level by combining the second- and third-generation sequencing technologies. The results revealed 185 (3198), 394 (3592), and 83 (3286) differentially expressed genes (transcripts) between tail and subcutaneous, tail and visceral, and subcutaneous and visceral adipose tissues, respectively. The results revealed that certain differentially expressed genes and transcripts were related to fat metabolism. These results provide a comprehensive dataset for further verification of the regulatory pathway linked to fat metabolism in sheep. Sheep have historically been bred globally as a vital food source. To explore the transcriptome of adipose tissue and investigate key genes regulating adipose metabolism in sheep, adipose tissue samples were obtained from F1 Dorper × Hu sheep. High-throughput sequencing libraries for second- and third-generation sequencing were constructed using extracted total RNA. Functional annotation of differentially expressed genes and isoforms facilitated the identification of key regulatory genes and isoforms associated with sheep fat metabolism. SMRT-seq generated 919,259 high-accuracy cDNA sequences after filtering. Full-length sequences were corrected using RNA-seq sequences, and 699,680 high-quality full-length non-chimeric (FLNC) reads were obtained. Upon evaluating the ratio of total lengths based on FLNC sequencing, it was determined that 36,909 out of 56,316 multiple-exon isoforms met the criteria for full-length status. This indicates the identification of 330,375 full-length FLNC transcripts among the 370,114 multiple-exon FLNC transcripts. By comparing the reference genomes, 60,276 loci and 111,302 isoforms were identified. In addition, 43,423 new genes and 44,563 new isoforms were identified. The results identified 185 (3198), 394 (3592), and 83 (3286) differentially expressed genes (transcripts) between tail and subcutaneous, tail and visceral, and subcutaneous and visceral adipose tissues, respectively. Functional annotation and pathway analysis revealed the following observations. (1) Among the differentially expressed genes (DEGs) of TF and SF tissues, the downregulation of ACADL, ACSL6, and NC_056060.1.2536 was observed in SF, while FFAR4 exhibited upregulation. (2) Among the DEGs of TF and VF tissues, expressions of ACADL, ACSL6, COL1A1, COL1A2, and SCD were downregulated in VF, with upregulation of FFAR4. (3) Among SF and VF expressions of COL1A1, COL1A2, and NC_056060.1.2536 were downregulated in VF. Specific differentially expressed genes (ACADL, ACSL6, COL1A1, COL1A2, FFAR4, NC_056060.1.2536, and SCD) and transcripts (NC_056066.1.1866.16 and NC_056066.1.1866.22) were identified as relevant to fat metabolism. These results provide a dataset for further verification of the regulatory pathway associated with fat metabolism in sheep. [ABSTRACT FROM AUTHOR]

SUMMARY: Bougainvillea is a typical tropical flower of great ornamental value due to its colorful bracts. The molecular mechanism behind color formation is not well‐understood. Therefore, this research conducted metabolome analysis, transcriptome analysis, and multi‐flux full‐length sequencing in two color bracts of Bougainvillea × buttiana 'Chitra' to investigate the significantly different metabolites (SDMs) and differentially expressed genes (DEGs). Overall, 261 SDMs, including 62 flavonoids and 26 alkaloids, were detected, and flavonols and betalains were significantly differentially accumulated among the two bracts. Furthermore, the complete‐length transcriptome of Bougainvillea × buttiana was also developed, which contained 512 493 non‐redundant isoforms. Among them, 341 210 (66.58%) displayed multiple annotations in the KOG, GO, NR, KEGG, Pfam, Swissprot, and NT databases. RNA‐seq findings revealed that 3610 DEGs were identified between two bracts. Co‐expression analysis demonstrated that the DEGs and SDMs involved in flavonol metabolism (such as CHS, CHI, F3H, FLS, CYP75B1, kaempferol, and quercetin) and betacyanin metabolism (DODA, betanidin, and betacyanins) were the main contributors for the canary yellow and red bract formation, respectively. Further investigation revealed that several putative transcription factors (TFs) might interact with the promoters of the genes mentioned above. The expression profiles of the putative TFs displayed that they may positively and negatively regulate the structural genes' expression profiles. The data revealed a potential regulatory network between important genes, putative TFs, and metabolites in the flavonol and betacyanin biosynthesis of Bougainvillea × buttiana 'Chitra' bracts. These findings will serve as a rich genetic resource for future studies that could create new color bracts. [ABSTRACT FROM AUTHOR]

Chinese cabbage is an economically important Brassica vegetable worldwide, and clubroot, which is caused by the soil-borne protist plant pathogen Plasmodiophora brassicae is regarded as a destructive disease to Brassica crops. Previous studies on the gene transcripts related to Chinese cabbage resistance to clubroot mainly employed RNA-seq technology, although it cannot provide accurate transcript assembly and structural information. In this study, PacBio RS II SMRT sequencing was used to generate full-length transcriptomes of mixed roots at 0, 2, 5, 8, 13, and 22 days after P. brassicae infection in the clubroot-resistant line DH40R. Overall, 39 376 high-quality isoforms and 26 270 open reading frames (ORFs) were identified from the SMRT sequencing data. Additionally, 426 annotated long noncoding RNAs (lncRNAs), 56 transcription factor (TF) families, 1 883 genes with poly(A) sites and 1 691 alternative splicing (AS) events were identified. Furthermore, 1 201 of the genes had at least one AS event in DH40R. A comparison with RNA-seq data revealed six differentially expressed AS genes (one for disease resistance and five for defensive response) that are potentially involved in P. brassicae resistance. The results of this study provide valuable resources for basic research on clubroot resistance in Chinese cabbage.

Liquidambar formosana is an excellent landscape ecological tree species because its beautiful tree shape and red or orange leaves in autumn. In order to understand the genetic basis of discoloration and secondary metabolism of L. formosana leaves, the mixed samples of L. formosana leaves in five leaf discoloration stages were used for full-length transcriptome sequencing using single-molecule real-time sequencing technology (PacBio platform). The results were as follows: (1) High-quality 41.04 Gb data were obtained by full-length transcriptome sequencing, from which 563 180 full-length non-chimeric sequences were identified, and 27 269 high-quality full-length transcripts were obtained by clustering and de-redundancy. In 27 269 full-length transcripts, 2 035 long-chain non-coding RNA (lncRNA) were predicted, and 14 892 simple sequence repeat (SSR) sites and 1 856 transcription factors were detected. (2) The results of gene annotation showed that a total of 24 857 transcripts were annotated in eight databases such as NR, GO, COG and KEGG, etc., 124 metabolic pathways were obtained in KEGG database, including ribosome, carbon metabolism, amino acid biosynthesis and so on, and 49 and 71 transcripts were involved in flavonoid and chlorophyll metabolism respectively. The above results preliminarily reveal the transcriptome information and functional characteristics of L. formosana leaves during the leaf discoloration stage, and provide basic data for the follow-up study of the molecular mechanism of leaf discoloration, the pathway and regulation of pigment metabolism and synthesis, the cloning of related functional genes and the improvement of leaf color.

Platycladus orientalis (P. orientalis) is a common tree used for vegetation restoration in northern China, and its large area propagation helps to improve site conditions. However, under harsh conditions such as poor land, the survival rate of P. orientalis is very low. Numerous studies have shown that root pruning can promote the formation of lateral roots in seedlings, enhancing the roots’ capacity to absorb soil nutrients and water, and thereby improving the survival rate of seedlings. In this study, a one-third root pruning treatment was applied to P. orientalis seedlings, and the whole transcriptome of seedlings subjected to both control (CK) and root pruning treatments was sequenced to analyze their gene expression profiles. This study investigated the regulatory mechanisms of lateral root development in response to root pruning damage at the molecular level. Using nine cells, 15.28 Gb of clean data were obtained, which yielded 101,688 high-quality full-length transcript sequences and 22,955 low-quality full-length transcript sequences after clustering. Redundancy was then removed using CD-HIT, and Illumina RNA-seq sequencing produced 139.26 Gb of clean data. A total of 2025 differentially expressed genes (DEGs) were identified at three time points following root pruning treatment. Enrichment analysis revealed that the peroxidase gene family plays a significant role in lateral root proliferation. Furthermore, the expression levels of the peroxidase gene family were notably upregulated in comparison to the control group. Pathway enrichment analysis identified 22 relevant genes, which appeared to be highly associated with root growth and resilience to stress. Through examining the expression patterns and correlations of these genes, five central genes emerged as key players. The findings of this research suggest that the peroxidase gene family plays a crucial role in the stress response and root development of P. orientalis, providing reference and guidance for root development in other plant species.

Comastoma polycladum is one of the original plant of Tibetan medicine “Zangyinchen”, which contains abundant medical components. To further know the transcriptome of C. polycladum and enrich its genetic information of gene annotation and metabolic pathway, the Pacbio sequencing platform was used to perform full-length transcriptome sequencing of C. polycladum leaves. The results were as follows: (1) A total of 17 Gb of sequencing data was collected, and 87 814 high-quality full-length transcripts were obtained by clustering and de-redundancy of 795 698 circular consistency sequences (CCSs) sequences. (2) Comparing with seven databases, 277 451 transcripts were annotated successfully, and in NR database with 39 214 transcripts annotated the most transcripts. A total of 26 396 transcripts were annotated to the KOG database, with 26 subcategories, and a total of 39 104 transcripts with six major pathways and 40 secondary pathways to the KEGG database. A total of 39 102 transcripts were annotated to the GO database, which were divided into three major categories: molecular function, biological process and cellular component. (3) SSR analysis yielded 22 861 SSRs, with single-base repeats being the most abundant. A total of 1 874 transcription factors and 15 166 long non-coding RNAs (LncRNAs) were identified, and the C3H transcription factor family had the most annotated transcripts. (4) A total of 55 transcripts involved in the synthesis of monoterpenes and flavonoids were screened out. These results enrich the transcriptome information of C. polycladum, and provide significant genetic resources for further screening of candidate genes related to the synthesis of medicinal components of C. polycladum.