A survey of the sorghum transcriptome using single-molecule long reads

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

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

Comparative transcriptome analysis of Cyperus esculentus and C. rotundus with contrasting oil contents in tubers defines genes and regulatory networks involved in oil accumulation

Plant vegetative organs present great potential for lipid storage, with tubers of Cyperus esculentus as a unique example. To investigate the genome and transcriptomic features of C. esculentus and related species, we sequenced and assembled the C. esculentus genome at the contig level. Through a comparative study of high-quality transcriptomes across 36 tissues from high-oil and intermediate-oil C. esculentus and low-oil Cyperus rotundus, we identified potential genes and regulatory networks related to tuber oil accumulation. First, we identified tuber-specific genes in two C. esculentus cultivars. Second, genes involved in fatty acid (FA) biosynthesis, triacylglycerol synthesis, and TAG packaging presented increased activity in the later stages of tuber development. Notably, tubers with high oil contents presented higher levels of these genes than those with intermediate oil contents did, whereas tubers with low oil contents presented minimal gene expression. Notably, a large fragment of the FA biosynthesis rate-limiting enzyme-encoding gene BCCP1 was missing from the C. rotundus transcript, which might be responsible for blocking FA biosynthesis in its tubers. WGCNA pinpointed a gene module linked to tuber oil accumulation, with a coexpression network involving the transcription factors WRI1, MYB4, and bHLH68. The ethylene-related genes in this module suggest a role for ethylene signaling in oil accumulation, which is supported by the finding that ethylene (ETH) treatment increases the oil content in C. esculentus tubers. This study identified potential genes and networks associated with tuber oil accumulation in C. esculentus, highlighting the role of specific genes, transcription factors, and ethylene signaling in this process.

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

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

The full-length transcriptional of the multiple spatiotemporal embryo-gonad tissues in chicken (Gallus gallus)

OBJECTIVES

Chicken (Gallus gallus), as the most economically important poultry, is a classical and ideal model for studying the mechanism of vertebrate developmental biology and embryology. However, the sex determination and differentiation in chicken is still elusive, which limited the application and slowed down many basic studies in chicken.

DATA DESCRIPTION

We applied PacBio Iso-seq to multiple spatiotemporal embryo-gonad tissues in the male and female chicken, which contain the blastoderm (E0, un-differentiation stage), genital ridge (E3.5-6.5, sex-differentiation stage) and gonads (E18.5, full-sex-differentiation stage). We obtained 51,479 and 48,356 full-length transcripts in male and female chicken embryo, respectively. The comprehensive annotated and evaluated these transcripts. The 1,293 and 1,556 candidate lncRNAs, 5,766 and 4,211 AS events in male and female. Collectively, our data constitutes a grand increase in the known number of lncRNA, AS (Alternative splicing) and Poly(A) during chicken embryo sex-differentiation and plays an important role in improving current genome annotation. In the meantime, the data will be enriched the functional studies in other birds.

Full-Length Transcriptome Profile of Apis cerana Revealed by Nanopore Sequencing

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.

Integrative analyses of long and short-read RNA sequencing reveal the spliced isoform regulatory network of seedling growth dynamics in upland cotton

The polyploid genome of cotton has significantly increased the transcript complexity. Recent advances in full-length transcript sequencing are now widely used to characterize the complete landscape of transcriptional events. Such studies in cotton can help us to explore the genetic mechanisms of the cotton seedling growth. Through long-read single-molecule RNA sequencing, this study compared the transcriptomes of three yield contrasting genotypes of upland cotton. Our analysis identified different numbers of spliced isoforms from 31,166, 28,716, and 28,713 genes in SJ48, Z98, and DT8 cotton genotypes, respectively, most of which were novel compared to previous cotton reference transcriptomes, and showed significant differences in the number of exon structures and coding sequence length due to intron retention. Quantification of isoform expression revealed significant differences in expression in the root and leaf of each genotype. An array of key isoform target genes showed protein kinase or phosphorylation functions, and their protein interaction network contained most of the circadian oscillator proteins. Spliced isoforms from the GIGANTEA (GI) protien were differentially regulated in each genotype and might be expected to regulate translational activities, including the sequence and function of target proteins. In addition, these spliced isoforms generate diurnal expression profiles in cotton leaves, which may alter the transcriptional regulatory network of seedling growth. Silencing of the novel spliced GI isoform Gh_A02G0645_N17 significantly affected biomass traits, contributed to variable growth, and increased transcription of the early flowering pathway gene ELF in cotton. Our high-throughput hybrid sequencing results will be useful to dissect functional differences among spliced isoforms in the polyploid cotton genome.

A chromosome-level reference genome assembly and a full-length transcriptome assembly of the giant freshwater prawn (Macrobrachium rosenbergii)

The giant freshwater prawn (Macrobrachium rosenbergii) is a key species in the aquaculture industry in several Asian, African, and South American countries. Despite a considerable growth in its production worldwide, the genetic complexities of M. rosenbergii various morphotypes pose challenges in cultivation. This study reports the first chromosome-scale reference genome and a high-quality full-length transcriptome assembly for M. rosenbergii. We employed the PacBio High Fidelity (HiFi) sequencing to obtain an initial draft assembly and further scaffolded it with the chromatin contact mapping (Hi-C) technique to achieve a final assembly of 3.73-Gb with an N50 scaffold length of 33.6 Mb. Repetitive elements constituted nearly 60% of the genome assembly, with simple sequence repeats and retrotransposons being the most abundant. The availability of both the chromosome-scale assembly and the full-length transcriptome assembly enabled us to thoroughly probe alternative splicing events in M. rosenbergii. Among the 2,041 events investigated, exon skipping represented the most prevalent class, followed by intron retention. Interestingly, specific isoforms were observed across multiple tissues. Additionally, within a single tissue type, transcripts could undergo alternative splicing, yielding multiple isoforms. We believe that the availability of a chromosome-level reference genome for M. rosenbergii, along with its full-length transcriptome, will be instrumental in advancing our understanding of the giant freshwater prawn biology and enhancing its molecular breeding programs, paving the way for the development of M. rosenbergii with valuable traits in commercial aquaculture.

A crosswalk on the genetic and conventional strategies for enhancing astaxanthin production in Haematococcus pluvialis

Astaxanthin is a naturally occurring xanthophyll with powerful: antioxidant, antitumor, and antibacterial properties that are widely employed in food, feed, medicinal and nutraceutical industries. Currently, chemical synthesis dominates the world's astaxanthin market, but the increasing demand for natural products is shifting the market for natural astaxanthin. Haematococcus pluvialis (H. pluvialis) is the factory source of natural astaxanthin when grown in optimal conditions. Currently, various strategies for the production of astaxanthin have been proposed or are being developed in order to meet its market demand. This up-to-date review scrutinized the current approaches or strategies that aim to increase astaxanthin yield from H. pluvialis. We have emphasized the genetic and environmental parameters that increase astaxanthin yield. We also looked at the transcriptomic dynamics caused by environmental factors (phytohormones induction, light, salt, temperature, and nutrient starvation) on astaxanthin synthesizing genes and other metabolic changes. Genetic engineering and culture optimization (environmental factors) are effective approaches to producing more astaxanthin for commercial purposes. Genetic engineering, in particular, is accurate, specific, potent, and safer than conventional random mutagenesis approaches. New technologies, such as CRISPR-Cas9 coupled with omics and emerging computational tools, may be the principal strategies in the future to attain strains that can produce more astaxanthin. This review provides accessible data on the strategies to increase astaxanthin accumulation natively. Also, this review can be a starting point for new scholars interested in H. pluvialis research.

High-quality chromosome scale genome assemblies of two important Sorghum inbred lines, Tx2783 and RTx436

Sorghum bicolor (L.) Moench is a significant grass crop globally, known for its genetic diversity. High quality genome sequences are needed to capture the diversity. We constructed high-quality, chromosome-level genome assemblies for two vital sorghum inbred lines, Tx2783 and RTx436. Through advanced single-molecule techniques, long-read sequencing and optical maps, we improved average sequence continuity 19-fold and 11-fold higher compared to existing Btx623 v3.0 reference genome and obtained 19 and 18 scaffolds (N50 of 25.6 and 14.4) for Tx2783 and RTx436, respectively. Our gene annotation efforts resulted in 29 612 protein-coding genes for the Tx2783 genome and 29 265 protein-coding genes for the RTx436 genome. Comparative analyses with 26 plant genomes which included 18 sorghum genomes and 8 outgroup species identified around 31 210 protein-coding gene families, with about 13 956 specific to sorghum. Using representative models from gene trees across the 18 sorghum genomes, a total of 72 579 pan-genes were identified, with 14% core, 60% softcore and 26% shell genes. We identified 99 genes in Tx2783 and 107 genes in RTx436 that showed functional enrichment specifically in binding and metabolic processes, as revealed by the GO enrichment Pearson Chi-Square test. We detected 36 potential large inversions in the comparison between the BTx623 Bionano map and the BTx623 v3.1 reference sequence. Strikingly, these inversions were notably absent when comparing Tx2783 or RTx436 with the BTx623 Bionano map. These inversion were mostly in the pericentromeric region which is known to have low complexity regions and harder to assemble and suggests the presence of potential artifacts in the public BTx623 reference assembly. Furthermore, in comparison to Tx2783, RTx436 exhibited 324 883 additional Single Nucleotide Polymorphisms (SNPs) and 16 506 more Insertions/Deletions (INDELs) when using BTx623 as the reference genome. We also characterized approximately 348 nucleotide-binding leucine-rich repeat (NLR) disease resistance genes in the two genomes. These high-quality genomes serve as valuable resources for discovering agronomic traits and structural variation studies.

Source leaves are regulated by sink strengths through non-coding RNAs and alternative polyadenylation in cucumber (Cucumis sativus L.)

BACKGROUND

The yield of major crops is generally limited by sink capacity and source strength. Cucumber is a typical raffinose family oligosaccharides (RFOs)-transporting crop. Non-coding RNAs and alternative polyadenylation (APA) play important roles in the regulation of growth process in plants. However, their roles on the sink‒source regulation have not been demonstrated in RFOs-translocating species.

RESULTS

Here, whole-transcriptome sequencing was applied to compare the leaves of cucumber under different sink strength, that is, no fruit-carrying leaves (NFNLs) and fruit-carrying leaves (FNLs) at 12th node from the bottom. The results show that 1101 differentially expressed (DE) mRNAs, 79 DE long non-coding RNAs (lncRNAs) and 23 DE miRNAs were identified, which were enriched in photosynthesis, energy production and conversion, plant hormone signal transduction, starch and carbohydrate metabolism and protein synthesis pathways. Potential co-expression networks like, DE lncRNAs-DE mRNAs/ DE miRNAs-DE mRNAs, and competing endogenous RNA (ceRNA) regulation models (DE lncRNAs-DE miRNAs-DE mRNAs) associated with sink‒source allocation, were constructed. Furthermore, 37 and 48 DE genes, which enriched in MAPK signaling and plant hormone signal transduction pathway, exist differentially APA, and SPS (CsaV3_2G033300), GBSS1 (CsaV3_5G001560), ERS1 (CsaV3_7G029600), PNO1 (CsaV3_3G003950) and Myb (CsaV3_3G022290) may be regulated by both ncRNAs and APA between FNLs and NFNLs, speculating that ncRNAs and APA are involved in the regulation of gene expression of cucumber sink‒source carbon partitioning.

CONCLUSIONS

These results reveal a comprehensive network among mRNAs, ncRNAs, and APA in cucumber sink-source relationships. Our findings also provide valuable information for further research on the molecular mechanism of ncRNA and APA to enhance cucumber yield.

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

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

Maize splicing-mediated mRNA surveillance impeded by sugarcane mosaic virus-coded pathogenic protein NIa-Pro

The eukaryotic mRNA surveillance pathway, a pivotal guardian of mRNA fidelity, stands at the nexus of diverse biological processes, including antiviral immunity. Despite the recognized function of splicing factors on mRNA fate, the intricate interplay shaping the mRNA surveillance pathway remains elusive. We illustrate that the conserved splicing factor U2 snRNP auxiliary factor large subunit B (U2AF65B) modulates splicing of mRNA surveillance complex, contributing to transcriptomic homeostasis in maize. The functionality of the mRNA surveillance pathway requires ZmU2AF65B-mediated normal splicing of upstream frameshift 3 (ZmUPF3) pre-mRNA, encoding a core factor in this pathway. Intriguingly, sugarcane mosaic virus (SCMV)-coded nuclear inclusion protein a protease (NIa-Pro) hinders the splicing function of ZmU2AF65B. Furthermore, NIa-Pro disrupts ZmU2AF65B binding to ZmUPF3 pre-mRNA, leading to dysregulated splicing of ZmUPF3 transcripts and, consequently, impairing mRNA surveillance, thus facilitating viral infection. Together, this study establishes that splicing governs the mRNA surveillance pathway and identifies a pathogenic protein capable of disrupting this regulation to compromise RNA immunity.

Direct RNA sequencing in plants: Practical applications and future perspectives

The transcriptome serves as a bridge that links genomic variation to phenotypic diversity. A vast number of studies using next-generation RNA sequencing (RNA-seq) over the last 2 decades have emphasized the essential roles of the plant transcriptome in response to developmental and environmental conditions, providing numerous insights into the dynamic changes, evolutionary traces, and elaborate regulation of the plant transcriptome. With substantial improvement in accuracy and throughput, direct RNA sequencing (DRS) has emerged as a new and powerful sequencing platform for precise detection of native and full-length transcripts, overcoming many limitations such as read length and PCR bias that are inherent to short-read RNA-seq. Here, we review recent advances in dissecting the complexity and diversity of plant transcriptomes using DRS as the main technological approach, covering many aspects of RNA metabolism, including novel isoforms, poly(A) tails, and RNA modification, and we propose a comprehensive workflow for processing of plant DRS data. Many challenges to the application of DRS in plants, such as the need for machine learning tools tailored to plant transcriptomes, remain to be overcome, and together we outline future biological questions that can be addressed by DRS, such as allele-specific RNA modification. This technology provides convenient support on which the connection of distinct RNA features is tightly built, sustainably refining our understanding of the biological functions of the plant transcriptome.

SMRT sequencing of a full-length transcriptome reveals cold induced alternative splicing in Vitis amurensis root

Alternative splicing enhances diversity at the transcriptional and protein levels that widely involved in plant response to biotic and abiotic stresses. V. amurensis is an extremely cold-tolerant wild grape variety, however, studies on alternative splicing (AS) in amur grape at low temperatures are currently poorly understood. In this study, we analyzed full-length transcriptome and RNA seq data at 0, 2, and 24 h after cold stress in V. amurensis roots. Following quality control and correction, 221,170 high-quality full-length non-concatemer (FLNC) reads were identified. A total of 16,181 loci and 30,733 isoforms were identified. These included 22,868 novel isoforms from annotated genes and 2815 isoforms from 2389 novel genes. Among the distinguished novel isoforms, 673 Long non-coding RNAs (LncRNAs) and 18,164 novel isoforms open reading frame (ORF) region were found. A total of 2958 genes produced 8797 AS events, of which 189 genes were involved in the low-temperature response. Twelve transcription factors show AS during cold treatment and VaMYB108 was selected for initial exploration. Two transcripts, Chr05.63.1 (VaMYB108short) and Chr05.63.2 (VaMYB108normal) of VaMYB108, display up-regulated expression after cold treatment in amur grape roots and are both localized in the nucleus. Only VaMYB108normal exhibits transcriptional activation activity. Overexpression of either VaMYB108short or VaMYB108normal in grape roots leads to increased expression of the other transcript and both increased chilling resistance of amur grape roots. The results improve and supplement the genome annotations and provide insights for further investigation into AS mechanisms during cold stress in V. amurensis.

High-throughput single-molecule long-read RNA sequencing analysis of tissue-specific genes and isoforms in lettuce (Lactuca sativa L.)

Lettuce is one of the most widely cultivated and consumed dicotyledonous vegetables globally. Despite the availability of its reference genome sequence, lettuce gene annotation remains incomplete, impeding comprehensive research and the broad application of genomic resources. Long-read RNA isoform sequencing (Iso-Seq) offers substantial advantages for analyzing RNA alternative splicing and aiding gene annotation, yet it faces throughput limitations. We present the HIT-ISOseq method tailored for bulk sample analysis, significantly enhancing RNA sequencing throughput on the PacBio platform by concatenating cDNA. Here we show, HIT-ISOseq generates 3-4 cDNA molecules per CCS read in lettuce, yielding 15.7 million long reads per PacBio Sequel II SMRT Cell 8 M. We validate its effectiveness in analyzing six lettuce tissue samples, including roots, stems, and leaves, revealing tissue-specific gene expression patterns and RNA isoforms. Leveraging diverse tissue long-read RNA sequencing, we refine the transcript annotation of the lettuce reference genome, expanding its GO and KEGG annotation repertoire. Collectively, this study serves as a foundational reference for genome annotation and the analysis of multi-sample isoform expression, utilizing high-throughput long-read transcriptome sequencing.

Accurate isoform quantification by joint short- and long-read RNA-sequencing

Accurate quantification of transcript isoforms is crucial for understanding gene regulation, functional diversity, and cellular behavior. Existing RNA sequencing methods have significant limitations: short-read (SR) sequencing provides high depth but struggles with isoform deconvolution, whereas long-read (LR) sequencing offers isoform resolution at the cost of lower depth, higher noise, and technical biases. Addressing this gap, we introduce Multi-Platform Aggregation and Quantification of Transcripts (MPAQT), a generative model that combines the complementary strengths of different sequencing platforms to achieve state-of-the-art isoform-resolved transcript quantification, as demonstrated by extensive simulations and experimental benchmarks. By applying MPAQT to an in vitro model of human embryonic stem cell differentiation into cortical neurons, followed by machine learning-based modeling of transcript abundances, we show that untranslated regions (UTRs) are major determinants of isoform proportion and exon usage; this effect is mediated through isoform-specific sequence features embedded in UTRs, which likely interact with RNA-binding proteins that modulate mRNA stability. These findings highlight MPAQT's potential to enhance our understanding of transcriptomic complexity and underline the role of splicing-independent post-transcriptional mechanisms in shaping the isoform and exon usage landscape of the cell.

Metabolism-Based Herbicide Resistance to Mesosulfuron-methyl and Identification of Candidate Genes in Bromus japonicus

The evolved resistance of Bromus japonicus Houtt. to ALS-inhibiting herbicides is well established. Previous studies have primarily focused on target-site resistance; however, non-target-site resistance has not been well characterized. This investigation demonstrated that ALS gene sequencing did not detect any previously known resistance mutations in a mesosulfuron-methyl-resistant (MR) population, and notably, treatment with the P450 monooxygenase (P450) inhibitor malathion markedly heightened susceptibility to mesosulfuron-methyl. Utilizing UPLC-MS/MS analysis confirmed elevated mesosulfuron-methyl metabolism in MR plants. The integration of Isoform Sequencing (Iso-Seq) and RNA Sequencing (RNA-Seq) facilitated the identification of candidate genes associated with non-target sites in a subpopulation with two generations of herbicide selection. Through qRT-PCR analysis, 21 differentially expressed genes were characterized, and among these, 10 genes (comprising three P450s, two glutathione S-transferases, one glycosyltransferase, two ATP-binding cassette transporters, one oxidase, and one hydrolase) exhibited constitutive upregulation in resistant plants. Our findings substantiated that increased herbicide metabolism is a driving force behind mesosulfuron-methyl resistance in this B. japonicus population.

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

FUS reads histone H3K36me3 to regulate alternative polyadenylation

Complex organisms generate differential gene expression through the same set of DNA sequences in distinct cells. The communication between chromatin and RNA regulates cellular behavior in tissues. However, little is known about how chromatin, especially histone modifications, regulates RNA polyadenylation. In this study, we found that FUS was recruited to chromatin by H3K36me3 at gene bodies. The H3K36me3 recognition of FUS was mediated by the proline residues in the ZNF domain. After these proline residues were mutated or H3K36me3 was abolished, FUS dissociated from chromatin and bound more to RNA, resulting in an increase in polyadenylation sites far from stop codons genome-wide. A proline mutation corresponding to a mutation in amyotrophic lateral sclerosis contributed to the hyperactivation of mitochondria and hyperdifferentiation in mouse embryonic stem cells. These findings reveal that FUS is an H3K36me3 reader protein that links chromatin-mediated alternative polyadenylation to human disease.

Comprehensive transcriptome sequencing of silkworm Midguts: Uncovering extensive isoform diversity and alternative splicing in BmNPV-Sensitive and BmNPV-resistant strains

The silkworm, Bombyx mori, stands out as one of the few economically valuable insects within the realm of model organisms. However, Bombyx mori nucleopolyhedrovirus (BmNPV) poses a significant threat, decreasing the quality and quantity of silkworm cocoons. Over the past few decades, a multitude of researchers has delved into the mechanisms that underlie silkworm resistance to BmNPV, employing diverse methodologies and approaching the problem from various angles. Despite this extensive research, the role of alternative splicing (AS) in the silkworm's response to BmNPV infection has been largely unexplored. This study leveraged both third-generation (Oxford Nanopore Technologies) and second-generation (Illumina) high-throughput sequencing technologies to meticulously identify and analyze AS patterns in the context of BmNPV response, utilizing two distinct silkworm strains-the susceptible strain 306 and the resistant strain NB. Consequently, we identified five crucial genes (Dsclp, LOC692903, LOC101743583, LOC101742498, LOC101743809) that are linked to the response to BmNPV infection through AS and differential expression. Additionally, a thorough comparative analysis was conducted on their diverse transcriptomic expression profiles, including alternative polyadenylation, simple sequence repeats, and transcription factors.

Chromosomal-level genome and metabolome analyses of highly heterozygous allohexaploid Dendrocalamus brandisii elucidate shoot quality and developmental characteristics

Dendrocalamus brandisii (Munro) Kurz is a sympodial bamboo species with inimitable taste and flavorful shoots. Its rapid growth and use as high-quality material make this bamboo species highly valued for both food processing and wood applications. However, genome information for D. brandisii is lacking, primarily due to its polyploidy and large genome size. Here, we assembled a high-quality genome for hexaploid D. brandisii, which comprises 70 chromosomes with a total size of 2,756 Mb, using long-read HiFi sequencing. Furthermore, we accurately separated the genome into its three constituent subgenomes. We used Oxford Nanopore Technologies long reads to construct a transcriptomic dataset covering 15 tissues for gene annotation to complement our genome assembly, revealing differential gene expression and post-transcriptional regulation. By integrating metabolome analysis, we unveiled that well-balanced lignin formation, as well as abundant flavonoid and fructose contents, contribute to the superior quality of D. brandisii shoots. Integrating genomic, transcriptomic, and metabolomic datasets provided a solid foundation for enhancing bamboo shoot quality and developing efficient gene-editing techniques. This study should facilitate research on D. brandisii and enhance its use as a food source and wood material by providing crucial genomic resources.

Multiple strategies, including 6mA methylation, affecting plant alternative splicing in allopolyploid peanut

Alternative splicing (AS), an important post-transcriptional regulation mechanism in eukaryotes, can significantly increase transcript diversity and contribute to gene expression regulation and many other complicated developmental processes. While plant gene AS events are well described, few studies have investigated the comprehensive regulation machinery of plant AS. Here, we use multi-omics to analyse peanut AS events. Using long-read isoform sequencing, 146 464 full-length non-chimeric transcripts were obtained, resulting in annotation corrections for 1782 genes and the identification of 4653 new loci. Using Iso-Seq RNA sequences, 271 776 unique splice junctions were identified, 82.49% of which were supported by transcriptome data. We characterized 50 977 polyadenylation sites for 23 262 genes, 12 369 of which had alternative polyadenylation sites. AS allows differential regulation of the same gene by miRNAs at the isoform level coupled with polyadenylation. In addition, we identified many long non-coding RNAs and fusion transcripts. There is a suppressed effect of 6mA on AS and gene expression. By analysis of chromatin structures, the genes located in the boundaries of topologically associated domains, proximal chromosomal telomere regions, inter- or intra-chromosomal loops were found to have more unique splice isoforms, higher expression, lower 6mA and more transposable elements (TEs) in their gene bodies than the other genes, indicating that chromatin interaction, 6mA and TEs play important roles in AS and gene expression. These results greatly refine the peanut genome annotation and contribute to the study of gene expression and regulation in peanuts. This work also showed AS is associated with multiple strategies for gene regulation.

A Comparative Full-Length Transcriptome Analysis Using Oxford Nanopore Technologies (ONT) in Four Tissues of Bovine Origin

The transcriptome complexity and splicing patterns in male and female cattle are ambiguous, presenting a substantial obstacle to genomic selection programs that seek to improve productivity, disease resistance, and reproduction in cattle. A comparative transcriptomic analysis using Oxford Nanopore Technologies (ONT) was conducted in bovine testes (TESTs), ovaries (OVAs), muscles (MUSCs), and livers (LIVs). An average of 5,144,769 full-length reads were obtained from each sample. The TESTs were found to have the greatest number of alternative polyadenylation (APA) events involved in processes such as sperm flagellum development and fertilization in male reproduction. In total, 438 differentially expressed transcripts (DETs) were identified in the LIVs in a comparison of females vs. males, and 214 DETs were identified in the MUSCs between females and males. Additionally, 14,735, 36,347, and 33,885 DETs were detected in MUSC vs. LIV, MUSC vs. TEST, and OVA vs. TEST comparisons, respectively, revealing the complexity of the TEST. Gene Set Enrichment Analysis (GSEA) showed that these DETs were mainly involved in the "spermatogenesis", "flagellated sperm motility", "spermatid development", "reproduction", "reproductive process", and "microtubule-based movement" KEGG pathways. Additional studies are necessary to further characterize the transcriptome in different cell types, developmental stages, and physiological conditions in bovines and ascertain the functions of the novel transcripts.

Pan-transcriptomic analysis reveals alternative splicing control of cold tolerance in rice

Plants have evolved complex mechanisms to adapt to harsh environmental conditions. Rice (Oryza sativa) is a staple food crop that is sensitive to low temperatures. However, its cold stress responses remain poorly understood, thus limiting possibilities for crop engineering to achieve greater cold tolerance. In this study, we constructed a rice pan-transcriptome and characterized its transcriptional regulatory landscape in response to cold stress. We performed Iso-Seq and RNA-Seq of 11 rice cultivars subjected to a time-course cold treatment. Our analyses revealed that alternative splicing-regulated gene expression plays a significant role in the cold stress response. Moreover, we identified CATALASE C (OsCATC) and Os03g0701200 as candidate genes for engineering enhanced cold tolerance. Importantly, we uncovered central roles for the 2 serine-arginine-rich proteins OsRS33 and OsRS2Z38 in cold tolerance. Our analysis of cold tolerance and resequencing data from a diverse collection of 165 rice cultivars suggested that OsRS2Z38 may be a key selection gene in japonica domestication for cold adaptation, associated with the adaptive evolution of rice. This study systematically investigated the distribution, dynamic changes, and regulatory mechanisms of alternative splicing in rice under cold stress. Overall, our work generates a rich resource with broad implications for understanding the genetic basis of cold response mechanisms in plants.

A Combinatorial Single-Molecule Real-Time and Illumina Sequencing Analysis of Postembryonic Gene Expression in the Asian Citrus Psyllid Diaphorina citri

Huanglongbing (HLB) is a systemic plant disease caused by 'Candidatus Liberibacter asiaticus (CLas)' and transmitted by Diaphorina citri. D. citri acquires the CLas bacteria in the nymph stage and transmits it in the adult stage, indicating that molting from the nymph to adult stages is crucial for HLB transmission. However, the available D. citri reference genomes are incomplete, and gene function studies have been limited to date. In the current research, PacBio single-molecule real-time (SMRT) and Illumina sequencing were performed to investigate the transcriptome of D. citri nymphs and adults. In total, 10,641 full-length, non-redundant transcripts (FLNRTs), 594 alternative splicing (AS) events, 4522 simple sequence repeats (SSRs), 1086 long-coding RNAs (lncRNAs), 281 transcription factors (TFs), and 4459 APA sites were identified. Furthermore, 3746 differentially expressed genes (DEGs) between nymphs and adults were identified, among which 30 DEGs involved in the Hippo signaling pathway were found. Reverse transcription-quantitative PCR (RT-qPCR) further validated the expression levels of 12 DEGs and showed a positive correlation with transcriptome data. Finally, the spatiotemporal expression pattern of genes involved in the Hippo signaling pathway exhibited high expression in the D. citri testis, ovary, and egg. Silencing of the D. citri transcriptional co-activator (DcYki) gene significantly increased D. citri mortality and decreased the cumulative molting. Our results provide useful information and a reliable data resource for gene function research of D. citri.

Integrative analysis of Iso-Seq and RNA-seq data reveals transcriptome complexity and differential isoform in skin tissues of different hair length Yak

BACKGROUND

The hair follicle development process is regulated by sophisticated genes and signaling networks, and the hair grows from the hair follicle. The Tianzhu white yak population exhibits differences in hair length, especially on the forehead and shoulder region. However, the genetic mechanism is still unclear. Isoform sequencing (Iso-seq) technology with advantages in long reads sequencing. Hence, we combined the Iso-seq and RNA-seq methods to investigate the transcript complexity and difference between long-haired yak (LHY) and normal-haired yak (NHY).

RESULTS

The hair length measurement result showed a significant difference between LHY and NHY on the forehead and the shoulder (P-value < 0.001). The skin samples from the forehead and the shoulder of LHY and NHY were pooled for isoform sequencing (Iso-seq). We obtained numerous long transcripts, including novel isoforms, long non-coding RNA, alternative splicing events, and alternative polyadenylation events. Combined with RNA-seq data, we performed differential isoforms (DEIs) analysis between LHY and NHY. We found that some hair follicle and skin development-related DEIs, like BMP4, KRT2, IGF2R, and COL1A2 in the forehead skin; BMP1, KRT1, FGF5, COL2A1, and IGFBP5 in the shoulder skin. Enrichment analysis revealed that DEIs in both two comparable groups significantly participated in skin and hair follicle development-related pathways, like ECM-receptor interaction, focal adhesion, and PI3K-Akt signaling pathways. The results indicated that the hair follicle development of Tianzhu white yak may influence the hair length difference. Besides, the protein-protein interaction (PPI) network of DEIs showed COL2A1 and COL3A1 exhibited a high degree of centrality, and these two genes were suggested as potential candidates for the hair length growth of Tianzhu white yak.

CONCLUSIONS

The results provided a comprehensive analysis of the transcriptome complexity and identified differential transcripts that enhance our understanding of the molecular mechanisms underlying the variation in hair length growth in Tianzhu white yak.

Full-length transcriptome analysis of a bloom-forming dinoflagellate Prorocentrum shikokuense (Dinophyceae)

Prorocentrum shikokuense (formerly P. donghaiense) is a pivotal dinoflagellate species associating with the HABs in the East China Sea. The complexity of its large nuclear genome hindered us from understanding its genomic characteristics. Full-length transcriptome sequencing offers a practical solution to decipher the physiological mechanisms of a species without the reference genome. In this study, we employed single-molecule real-time (SMRT) sequencing technology to sequence the full-length transcriptome of Prorocentrum shikokuense. We successfully generated 41.73 Gb of clean SMRT sequencing reads and isolated 105,249 non-redundant full-length non-chimeric reads. Our trial has led to the identification of 11,917 long non-coding RNA transcripts, 514 alternative splicing events, 437 putative transcription factor genes from 17 TF gene families, and 34,723 simple sequence repeats. Additionally, a total of 78,265 open reading frames were identified, of them 15,501 were the protein coding sequences. This dataset is valuable for annotating P. shikokuense genome, and will contribute significantly to the in-depth studies on the molecular mechanisms underlining the dinoflagellate bloom formation.

Genome-wide characterization of post-transcriptional processes related to wood formation in Dalbergia odorifera

BACKGROUND

Alternative polyadenylation (APA), alternative splicing (AS), and long non-coding RNAs (lncRNAs) play regulatory roles in post-transcriptional processes in plants. However, little is known about their involvement in xylem development in Dalbergia odorifera, a valuable rosewood species with medicinal and commercial significance. We addressed this by conducting Isoform Sequencing (Iso-Seq) using PacBio's SMRT technology and combined it with RNA-seq analysis (RNA sequencing on Illumina platform) after collecting xylem samples from the transition zone and the sapwood of D. odorifera.

RESULTS

We identified 14,938 full-length transcripts, including 9,830 novel isoforms, which has updated the D. odorifera genome annotation. Our analysis has revealed that 4,164 genes undergo APA, whereas 3,084 genes encounter AS. We have also annotated 118 lncRNAs. Furthermore, RNA-seq analysis identified 170 differential alternative splicing (DAS) events, 344 genes with differential APA site usage (DE-APA), and 6 differentially expressed lncRNAs in the transition zone when compared to the sapwood. AS, APA, and lncRNAs are differentially regulated during xylem development. Differentially expressed APA genes were enriched for terpenoid and flavonoid metabolism, indicating their role in the heartwood formation. Additionally, DE-APA genes were associated with cell wall biosynthesis and terpenoid metabolism, implying an APA's role in wood formation. A DAS gene (involved in chalcone accumulation) with a significantly greater inclusion of the last exon in the transition zone than in the sapwood was identified. We also found that differentially expressed lncRNAs targeted the genes related to terpene synthesis.

CONCLUSIONS

This study enhances our understanding of the molecular regulatory mechanisms underlying wood formation in D. odorifera, and provides valuable genetic resources and insights for its molecular-assisted breeding.

High-quality faba bean reference transcripts generated using PacBio and Illumina RNA-seq data

The genome of faba bean was first published in 2023. To promote future molecular breeding studies, we improved the quality of the faba genome based on high-density genetic maps and the Illumina and Pacbio RNA-seq datasets. Two high-density genetic maps were used to conduct the scaffold ordering and orientation of faba bean, culminating in an increased length (i.e., 14.28 Mbp) of chromosomes and a decrease in the number of scaffolds by 45. In gene model mining and optimisation, the PacBio and Illumina RNA-seq datasets from 37 samples allowed for the identification and correction 121,606 transcripts, and the data facilitated a prediction of 15,640 alternative splicing events, 2,148 lncRNAs, and 1,752 fusion transcripts, thus allowing for a clearer understanding of the gene structures underlying the faba genome. Moreover, a total of 38,850 new genes including 56,188 transcripts were identified compared with the reference genome. Finally, the genetic data of the reference genome was integrated and a comprehensive and complete faba bean transcriptome sequence of 103,267 transcripts derived from 54,753 uni-genes was formed.

ASAPA: a bioinformatic pipeline based on Iso-Seq that identifies the links among alternative splicing, alternative transcription initiation and alternative polyadenylation

BACKGROUND

Although the events associated with alternative splicing (AS), alternative polyadenylation (APA) and alternative transcription initiation (ATI) can be identified by many approaches based on isoform sequencing (Iso-Seq), these analyses are generally independent of each other and the links between these events are still rarely mentioned. However, an interdependency analysis can be achieved because the transcriptional start site, splice sites and polyA site could be simultaneously included in a long, full-length read from Iso-Seq.

RESULTS

We create ASAPA pipeline that enables streamlined analysis for a robust detection of potential links among AS, ATI and APA using Iso-Seq data. We tested this pipeline using Arabidopsis data and found some interesting results: some adjacent introns tend to be simultaneously spliced or retained; coupling between AS and ATI or APA is limited to the initial or terminal intron; and ATI and APA are potentially linked in some special cases.

CONCLUSION

Our pipeline enables streamlined analysis for a robust detection of potential links among AS, ATI and APA using Iso-Seq data, which is conducive to a better understanding of transcription landscape generation.

Bulk and single-cell alternative splicing analyses reveal roles of TRA2B in myogenic differentiation

Alternative splicing (AS) disruption has been linked to disorders of muscle development, as well as muscular atrophy. However, the precise changes in AS patterns that occur during myogenesis are not well understood. Here, we employed isoform long-reads RNA-seq (Iso-seq) and single-cell RNA-seq (scRNA-seq) to investigate the AS landscape during myogenesis. Our Iso-seq data identified 61,146 full-length isoforms representing 11,682 expressed genes, of which over 52% were novel. We identified 38,022 AS events, with most of these events altering coding sequences and exhibiting stage-specific splicing patterns. We identified AS dynamics in different types of muscle cells through scRNA-seq analysis, revealing genes essential for the contractile muscle system and cytoskeleton that undergo differential splicing across cell types. Gene-splicing analysis demonstrated that AS acts as a regulator, independent of changes in overall gene expression. Two isoforms of splicing factor TRA2B play distinct roles in myogenic differentiation by triggering AS of TGFBR2 to regulate canonical TGF-β signalling cascades differently. Our study provides a valuable transcriptome resource for myogenesis and reveals the complexity of AS and its regulation during myogenesis.

Nanopore Direct RNA Sequencing Reveals the Short-Term Salt Stress Response in Maize Roots

Transcriptome analysis, relying on the cutting-edge sequencing of cDNA libraries, has become increasingly prevalent within functional genome studies. However, the dependence on cDNA in most RNA sequencing technologies restricts their ability to detect RNA base modifications. To address this limitation, the latest Oxford Nanopore Direct RNA Sequencing (ONT DRS) technology was employed to investigate the transcriptome of maize seedling roots under salt stress. This approach aimed to unveil both the RNA transcriptional profiles and alterations in base modifications. The analysis of the differential expression revealed a total of 1398 genes and 2223 transcripts that exhibited significant variation within the maize root system following brief exposure to salt stress. Enrichment analyses, such as the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway assessments, highlighted the predominant involvement of these differentially expressed genes (DEGs) in regulating ion homeostasis, nitrogen metabolism, amino acid metabolism, and the phytohormone signaling pathways. The protein-protein interaction (PPI) analysis showed the participation of various proteins related to glycolytic metabolism, nitrogen metabolism, amino acid metabolism, abscisic acid signaling, and the jasmonate signaling pathways. It was through this intricate molecular network that these proteins collaborated to safeguard root cells against salt-induced damage. Moreover, under salt stress conditions, the occurrence of variable shear events (AS) in RNA modifications diminished, the average length of poly(A) tails underwent a slight decrease, and the number of genes at the majority of the variable polyadenylation (APA) sites decreased. Additionally, the levels of N5-methylcytosine (m5C) and N6-methyladenosine (m6A) showed a reduction. These results provide insights into the mechanisms of early salt tolerance in maize.

Full-length transcriptome and RNA-Seq analyses reveal the resistance mechanism of sesame in response to Corynespora cassiicola

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.

Full-Length Transcriptome and Gene Expression Analysis of Different Ovis aries Adipose Tissues Reveals Transcript Variants Involved in Lipid Biosynthesis

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.

The subgenome Saccharum spontaneum contributes to sugar accumulation in sugarcane as revealed by full-length transcriptomic analysis

INTRODUCTION

Modern sugarcane cultivars (Saccharum spp. hybrids) derived from crosses between S. officinarum and S. spontaneum, with high-sugar traits and excellent stress tolerance inherited respectively. However, the contribution of the S. spontaneum subgenome to sucrose accumulation is still unclear.

OBJECTIVE

To compensate for the absence of a high-quality reference genome, a transcriptome analysis method is needed to analyze the molecular basis of differential sucrose accumulation in sugarcane hybrids and to find clues to the contribution of the S. spontaneum subgenome to sucrose accumulation.

METHODS

PacBio full-length sequencing was used to complement genome annotation, followed by the identification of differential genes between the high and low sugar groups using differential alternative splicing analysis and differential expression analysis. At the subgenomic level, the factors responsible for differential sucrose accumulation were investigated from the perspective of transcriptional and post-transcriptional regulation.

RESULTS

A full-length transcriptome annotated at the subgenomic level was provided, complemented by 263,378 allele-defined transcript isoforms and 139,405 alternative splicing (AS) events. Differential alternative splicing (DA) analysis and differential expression (DE) analysis identified differential genes between high and low sugar groups and explained differential sucrose accumulation factors by the KEGG pathways. In some gene models, different or even opposite expression patterns of alleles from the same gene were observed, reflecting the potential evolution of these alleles toward novel functions in polyploid sugarcane. Among DA and DE genes in the sucrose source-sink complex pathway, we found some alleles encoding sucrose accumulation-related enzymes derived from the S. spontaneum subgenome were differentially expressed or had DA events between the two contrasting sugarcane hybrids.

CONCLUSION

Full-length transcriptomes annotated at the subgenomic level could better characterize sugarcane hybrids, and the S. spontaneum subgenome was found to contribute to sucrose accumulation.

Genomics and transcriptomics of the Chinese mitten crabs (Eriocheir sinensis)

To gain a deeper understanding of the genetic factors influencing the growth and development of Eriocheir sinensis, a well-known species of hairy crab found in Yangcheng Lake, this study focused on the de novo genome and full-length transcriptome information of the selected subjects. Specifically, Yangcheng Lake hairy crabs were chosen as the experimental samples. Initially, a genome analysis was performed, resulting in the identification of gene fragments with a combined length of 1266,092,319 bp. Subsequently, a transcriptome analysis was conducted on a mixture of tissues from four different sites, namely muscle, brain, eye, and heart, to further investigate the genetic characteristics at the transcriptome level. The Pacific Biosciences (Pacio) single-molecule real-time sequencing system generated a total of 36.93 G sub-fragments and 175,90041 effective inserts. This research contributes to the indirect comprehension of genetic variations underlying individual traits. Furthermore, a comparison of the obtained data with relevant literature emphasizes the advantages of this study and establishes a basis for further investigations on the Chinese mitten crab.

Alternative polyadenylation regulates acetyl-CoA carboxylase function in peanut

BACKGROUND

Polyadenylation is a crucial process that terminates mRNA molecules at their 3'-ends. It has been observed that alternative polyadenylation (APA) can generate multiple transcripts from a single gene locus, each with different polyadenylation sites (PASs). This leads to the formation of several 3' untranslated regions (UTRs) that vary in length and composition. APA has a significant impact on approximately 60-70% of eukaryotic genes and has far-reaching implications for cell proliferation, differentiation, and tumorigenesis.

RESULTS

In this study, we conducted long-read, single-molecule sequencing of mRNA from peanut seeds. Our findings revealed that over half of all peanut genes possess over two PASs, with older developing seeds containing more PASs. This suggesting that the PAS exhibits high tissue specificity and plays a crucial role in peanut seed maturation. For the peanut acetyl-CoA carboxylase A1 (AhACCA1) gene, we discovered four 3' UTRs referred to UTR1-4. RT-PCR analysis showed that UTR1-containing transcripts are predominantly expressed in roots, leaves, and early developing seeds. Transcripts containing UTR2/3 accumulated mainly in roots, flowers, and seeds, while those carrying UTR4 were constitutively expressed. In Nicotiana benthamiana leaves, we transiently expressed all four UTRs, revealing that each UTR impacted protein abundance but not subcellular location. For functional validation, we introduced each UTR into yeast cells and found UTR2 enhanced AhACCA1 expression compared to a yeast transcription terminator, whereas UTR3 did not. Furthermore, we determined ACC gene structures in seven plant species and identified 51 PASs for 15 ACC genes across four plant species, confirming that APA of the ACC gene family is universal phenomenon in plants.

CONCLUSION

Our data demonstrate that APA is widespread in peanut seeds and plays vital roles in peanut seed maturation. We have identified four 3' UTRs for AhACCA1 gene, each showing distinct tissue-specific expression patterns. Through subcellular location experiment and yeast transformation test, we have determined that UTR2 has a stronger impact on gene expression regulation compared to the other three UTRs.

Comprehensive Profiling of Alternative Splicing and Alternative Polyadenylation during Fruit Ripening in Watermelon (Citrullus lanatus)

Fruit ripening is a highly complicated process that is accompanied by the formation of fruit quality. In recent years, a series of studies have demonstrated post-transcriptional control play important roles in fruit ripening and fruit quality formation. Till now, the post-transcriptional mechanisms for watermelon fruit ripening have not been comprehensively studied. In this study, we conducted PacBio single-molecule long-read sequencing to identify genome-wide alternative splicing (AS), alternative polyadenylation (APA) and long non-coding RNAs (lncRNAs) in watermelon fruit. In total, 6,921,295 error-corrected and mapped full-length non-chimeric (FLNC) reads were obtained. Notably, more than 42,285 distinct splicing isoforms were derived from 5,891,183 intron-containing full-length FLNC reads, including a large number of AS events associated with fruit ripening. In addition, we characterized 21,506 polyadenylation sites from 11,611 genes, 8703 of which have APA sites. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that fructose and mannose metabolism, starch and sucrose metabolism and carotenoid biosynthesis were both enriched in genes undergoing AS and APA. These results suggest that post-transcriptional regulation might potentially have a key role in regulation of fruit ripening in watermelon. Taken together, our comprehensive PacBio long-read sequencing results offer a valuable resource for watermelon research, and provide new insights into the molecular mechanisms underlying the complex regulatory networks of watermelon fruit ripening.

Alternative Splicing Variation: Accessing and Exploiting in Crop Improvement Programs

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

Uncovering genomic and transcriptional variations facilitates utilization of wild resources in cotton disease resistance improvement

BACKGROUND

Upland cotton wild/landraces represent a valuable resource for disease resistance alleles. Genetic differentiation between genotypes, as well as variation in Verticillium wilt (VW) resistance, has been poorly characterized for upland cotton accessions on the domestication spectrum (from wild/landraces to elite lines).

RESULTS

To illustrate the effects of modern breeding on VW resistance in upland cotton, 37 wild/landraces were resequenced and phenotyped for VW resistance. Genomic patterns of differentiation were identified between wild/landraces and improved upland cotton, and a significant decline in VW resistance was observed in association with improvement. Four genotypes representing different degrees of improvement were used in a full-length transcriptome analysis to study the genetic basis of VW resistance. ROS signaling was highly conserved at the transcriptional level, likely providing the basis for VW resistance in upland cotton. ASN biosynthesis and HSP90-mediated resistance moderated the response to VW in wild/landraces, and loss of induction activity of these genes resulted in VW susceptibility. The observed genomic differentiation contributed to the loss of induction of some important VW resistance genes such as HSP90.4 and PR16.

CONCLUSIONS

Besides providing new insights into the evolution of upland cotton VW resistance, this study also identifies important resistance pathways and genes for both fundamental research and cotton breeding.

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

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

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

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

Mapping intron retention events contributing to complex traits using splice quantitative trait locus

BACKGROUND

Alternative splicing (AS) of mRNA plays an important roles in transcriptome diversity, involving regulation of plant growth and stress response. Understanding the variation of AS events underlying GWAS loci in a crop population can provide insight into the molecular mechanisms of complex agronomic traits. To date, genome-wide association studies relating AS events to agronomic traits have rarely been conducted at the population level in crops.

RESULTS

Here, a pipeline was constructed to identify candidate AS events related to complex traits. Firstly, ovule transcriptome data were used to characterize intron retention (IR), the predominant type of AS in plants, on a genome-wide scale. This was done in a natural population consisting of 279 upland cotton lines. Secondly, splice quantitative trait locus (sQTL) analysis was carried out, which yielded a total of 2295 sQTLs involving 1607 genes. Of these, 14.25% (n = 427) were cis-regulatory loci. Integration with expression quantitative trait loci (eQTL) revealed that 53 (21.4%) cis-sGenes were regulated by both cis-sQTLs and cis-eQTLs. Finally, co-localization analysis integrated with GWAS loci in this population showed 32 cis-QTLs to be co-located with genetic regulatory loci related to fiber yield and quality traits, indicating that sQTLs are likely to participate in regulating cotton fiber yield and quality. An in-depth evaluation confirmed that differences in the IR rates of sQTL-regulated candidate genes such as GhLRRK1 and GhGC1 are associated with lint percentage (LP), which has potential in breeding applications.

CONCLUSION

This study provides a clue that AS of mRNA has an impact on crop yield, along with functional sQTLs are new genetic resources for cotton precision breeding.

PacBio single-molecule long-read sequencing provides new insights into the complexity of full-length transcripts in oriental river prawn, macrobrachium nipponense

BACKGROUND

Oriental river prawn (Macrobrachium nipponense) is one of the most dominant species in shrimp farming in China, which is a rich source of protein and contributes to a significant impact on the quality of human life. Thus, more complete and accurate annotation of gene models are important for the breeding research of oriental river prawn.

RESULTS

A full-length transcriptome of oriental river prawn muscle was obtained using the PacBio Sequel platform. Then, 37.99 Gb of subreads were sequenced, including 584,498 circular consensus sequences, among which 512,216 were full length non-chimeric sequences. After Illumina-based correction of long PacBio reads, 6,599 error-corrected isoforms were identified. Transcriptome structural analysis revealed 2,263 and 2,555 alternative splicing (AS) events and alternative polyadenylation (APA) sites, respectively. In total, 620 novel genes (NGs), 197 putative transcription factors (TFs), and 291 novel long non-coding RNAs (lncRNAs) were identified.

CONCLUSIONS

In summary, this study offers novel insights into the transcriptome complexity and diversity of this prawn species, and provides valuable information for understanding the genomic structure and improving the draft genome annotation of oriental river prawn.

IsoTools: a flexible workflow for long-read transcriptome sequencing analysis

MOTIVATION

Long-read transcriptome sequencing (LRTS) has the potential to enhance our understanding of alternative splicing and the complexity of this process requires the use of versatile computational tools, with the ability to accommodate various stages of the workflow with maximum flexibility.

RESULTS

We introduce IsoTools, a Python-based LRTS analysis framework that offers a wide range of functionality for transcriptome reconstruction and quantification of transcripts. Furthermore, we integrate a graph-based method for identifying alternative splicing events and a statistical approach based on the beta-binomial distribution for detecting differential events. To demonstrate the effectiveness of our methods, we applied IsoTools to PacBio LRTS data of human hepatocytes treated with the histone deacetylase inhibitor valproic acid. Our results indicate that LRTS can provide valuable insights into alternative splicing, particularly in terms of complex and differential splicing patterns, in comparison to short-read RNA-seq.

AVAILABILITY AND IMPLEMENTATION

IsoTools is available on GitHub and PyPI, and its documentation, including tutorials, CLI, and API references, can be found at https://isotools.readthedocs.io/.

Cochlear transcript diversity and its role in auditory functions implied by an otoferlin short isoform

Isoforms of a gene may contribute to diverse biological functions. In the cochlea, the repertoire of alternative isoforms remains unexplored. We integrated single-cell short-read and long-read RNA sequencing techniques and identified 236,012 transcripts, 126,612 of which were unannotated in the GENCODE database. Then we analyzed and verified the unannotated transcripts using RNA-seq, RT-PCR, Sanger sequencing, and MS-based proteomics approaches. To illustrate the importance of identifying spliced isoforms, we investigated otoferlin, a key protein involved in synaptic transmission in inner hair cells (IHCs). Upon deletion of the canonical otoferlin isoform, the identified short isoform is able to support normal hearing thresholds but with reduced sustained exocytosis of IHCs, and further revealed otoferlin functions in endocytic membrane retrieval that was not well-addressed previously. Furthermore, we found that otoferlin isoforms are associated with IHC functions and auditory phenotypes. This work expands our mechanistic understanding of auditory functions at the level of isoform resolution.

Recent advancements in the breeding of sorghum crop: current status and future strategies for marker-assisted breeding

Sorghum is emerging as a model crop for functional genetics and genomics of tropical grasses with abundant uses, including food, feed, and fuel, among others. It is currently the fifth most significant primary cereal crop. Crops are subjected to various biotic and abiotic stresses, which negatively impact on agricultural production. Developing high-yielding, disease-resistant, and climate-resilient cultivars can be achieved through marker-assisted breeding. Such selection has considerably reduced the time to market new crop varieties adapted to challenging conditions. In the recent years, extensive knowledge was gained about genetic markers. We are providing an overview of current advances in sorghum breeding initiatives, with a special focus on early breeders who may not be familiar with DNA markers. Advancements in molecular plant breeding, genetics, genomics selection, and genome editing have contributed to a thorough understanding of DNA markers, provided various proofs of the genetic variety accessible in crop plants, and have substantially enhanced plant breeding technologies. Marker-assisted selection has accelerated and precised the plant breeding process, empowering plant breeders all around the world.

Alternative splicing affects synapses in the hippocampus of offspring after maternal fructose exposure during gestation and lactation

Increased fructose over-intake is a global issue. Maternal fructose exposure during gestation and lactation can impair brain development in offspring. However, the effect on synapses is still unknown. For the diversification of RNA and biological functions, alternative splicing (AS) and alternative polyadenylation (APA) are essential. We constructed a maternal high-fructose diet model by administering 13% and 40% fructose water. The student's t-test analyzed the results of RT-qPCR. All other results were analyzed by one-way analysis of variance. The animal behavior experiment results revealed that conditioning and associative memory had been damaged. The proteins that form synapses were consistently low-expressed. In addition, compared with the control group, the Oxford Nanopore Technologies platform's full-length RNA-sequencing identified 298 different spliced genes (DSGs) and 51 differentially expressed alternative splicing (DEAS) genes in the 13% fructose group. 313 DSGs and 74 DEAS genes were in the 40% fructose group. Enrichment analysis based on these altered genes revealed some enlightening items and pathways. Our findings demonstrated the transcriptome mechanism underlying maternal fructose exposure during gestation and lactation and impaired synapse function during the transcripts' editing.

Transcriptomic Complexity of Culm Growth and Development in Different Types of Moso Bamboo

Moso bamboo is capable of both sexual and asexual reproduction during natural growth, resulting in four distinct types of culms: the bamboo shoot-culm, the seedling stem, the leptomorph rhizome, and a long-ignored culm-the outward-rhizome. Sometimes, when the outward rhizomes break through the soil, they continue to grow longitudinally and develop into a new individual. However, the roles of alternative transcription start sites (aTSS) or termination sites (aTTS) as well as alternative splicing (AS) have not been comprehensively studied for their development. To re-annotate the moso bamboo genome and identify genome-wide aTSS, aTTS, and AS in growing culms, we utilized single-molecule long-read sequencing technology. In total, 169,433 non-redundant isoforms and 14,840 new gene loci were identified. Among 1311 lncRNAs, most of which showed a positive correlation with their target mRNAs, one-third of these IncRNAs were preferentially expressed in winter bamboo shoots. In addition, the predominant AS type observed in moso bamboo was intron retention, while aTSS and aTTS events occurred more frequently than AS. Notably, most genes with AS events were also accompanied by aTSS and aTTS events. Outward rhizome growth in moso bamboo was associated with a significant increase in intron retention, possibly due to changes in the growth environment. As different types of moso bamboo culms grow and develop, a significant number of isoforms undergo changes in their conserved domains due to the regulation of aTSS, aTTS, and AS. As a result, these isoforms may play different roles than their original functions. These isoforms then performed different functions from their original roles, contributing to the transcriptomic complexity of moso bamboo. Overall, this study provided a comprehensive overview of the transcriptomic changes underlying different types of moso bamboo culm growth and development.