Comprehensive assessment of mRNA isoform detection methods for long-read sequencing data.

The advancement of Long-Read Sequencing (LRS) techniques has significantly increased the length of sequencing to several kilobases, thereby facilitating the identification of alternative splicing events and isoform expressions. Recently, numerous computational tools for isoform detection using long-read sequencing data have been developed. Nevertheless, there remains a deficiency in comparative studies that systemically evaluate the performance of these tools, which are implemented with different algorithms, under various simulations that encompass potential influencing factors. In this study, we conducted a benchmark analysis of thirteen methods implemented in nine tools capable of identifying isoform structures from long-read RNA-seq data. We evaluated their performances using simulated data, which represented diverse sequencing platforms generated by an in-house simulator, RNA sequins (sequencing spike-ins) data, as well as experimental data. Our findings demonstrate IsoQuant as a highly effective tool for isoform detection with LRS, with Bambu and StringTie2 also exhibiting strong performance. These results offer valuable guidance for future research on alternative splicing analysis and the ongoing improvement of tools for isoform detection using LRS data.

The StuA Transcription Factor and Alternative Splicing Mechanisms Drive the Levels of MAPK Hog1 Transcripts in the Dermatophyte Trichophyton rubrum.

Trichophyton rubrum is a human fungal pathogen that causes dermatophytosis, an infection that affects keratinized tissues. Integrated molecular signals coordinate mechanisms that control pathogenicity. Transcriptional regulation is a core regulation of relevant fungal processes. Previous RNA sequencing data revealed that the absence of the transcription factor StuA resulted in the differential expression of the MAPK-related high glycerol osmolarity gene (hog1) in T. rubrum. Here we validated the role of StuA in regulating the transcript levels of hog1. We showed through RT-qPCR that transcriptional regulation controls hog1 levels in response to glucose, keratin, and co-culture with human keratinocytes. In addition, we also detected hog1 pre-mRNA transcripts that underwent alternative splicing, presenting intron retention in a StuA-dependent mechanism. Our findings suggest that StuA and alternative splicing simultaneously, but not dependently, coordinate hog1 transcript levels in T. rubrum. As a means of preventing and treating dermatophytosis, our results contribute to the search for new potential drug therapies based on the molecular aspects of signaling pathways in T. rubrum.

Alternative Splicing Plays a Crucial Role in the Salt Tolerance of Foxtail Millet.

Foxtail millet is an important cereal crop that is relatively sensitive to salt stress, with its yield significantly affected by such stress. Alternative splicing (AS) widely affects plant growth, development, and adaptability to stressful environments. Through RNA-seq analysis of foxtail millet under different salt treatment periods, 2078 AS events were identified, and analyses were conducted on differential gene (DEG), differential alternative splicing gene (DASG), and overlapping gene. To investigate the regulatory mechanism of AS in response to salt stress in foxtail millet, the foxtail millet AS genes SiCYP19, with two AS variants (SiCYP19-a and SiCYP19-b), were identified and cloned. Yeast overexpression experiments indicated that SiCYP19 may be linked to the response to salt stress. Subsequently, we conducted overexpression experiments of both alternative splicing variants in foxtail millet roots to validate them experimentally. The results showed that, under salt stress, both SiCYP19-a and SiCYP19-b jointly regulated the salt tolerance of foxtail millet. Specifically, overexpression of SiCYP19-b significantly increased the proline content and reduced the accumulation of reactive oxygen species (ROS) in foxtail millet, compared to that in SiCYP19-a. This shows that SiCYP19-b plays an important role in increasing the content of proline and promoting the clearance of ROS, thus improving the salt tolerance of foxtail millet.

High expression of SRSF1 facilitates osteosarcoma progression and unveils its potential mechanisms.

BACKGROUND: SRSF1, a member of Serine/Arginine-Rich Splicing Factors (SRSFs), has been observed to significantly influence cancer progression. However, the precise role of SRSF1 in osteosarcoma (OS) remains unclear. This study aims to investigate the functions of SRSF1 and its underlying mechanism in OS. METHODS: SRSF1 expression level in OS was evaluated on the TCGA dataset, TAGET-OS database. qRT-PCR and Western blotting were employed to assess SRSF1 expression in human OS cell lines as well as the interfered ectopic expression states. The effect of SRSF1 on cell migration, invasion, proliferation, and apoptosis of OS cells were measured by transwell assay and flow cytometry. RNA sequence and bioinformatic analyses were conducted to elucidate the targeted genes, relevant biological pathways, and alternative splicing (AS) events regulated by SRSF1. RESULTS: SRSF1 expression was consistently upregulated in both OS samples and OS cell lines. Diminishing SRSF1 resulted in reduced proliferation, migration, and invasion and increased apoptosis in OS cells while overexpressing SRSF1 led to enhanced growth, migration, invasion, and decreased apoptosis. Mechanistically, Gene Ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, and Gene Set Enrichment Analysis (GSEA) revealed that the biological functions of SRSF1 were closely associated with the dysregulation of the protein targeting processes, location of the cytosolic ribosome, extracellular matrix (ECM), and proteinaceous extracellular matrix, along with the PI3K-AKT pathway, Wnt pathway, and HIPPO pathway. Transcriptome analysis identified AS events modulated by SRSF1, especially (Skipped Exon) SE events and (Mutually exclusive Exons) MXE events, revealing potential roles of targeted molecules in mRNA surveillance, RNA degradation, and RNA transport during OS development. qRT-PCR confirmed that SRSF1 knockdown resulted in the occurrence of alternative splicing of SRRM2, DMKN, and SCAT1 in OS. CONCLUSIONS: Our results highlight the oncogenic role of high SRSF1 expression in promoting OS progression, and further explore the potential mechanisms of action. The significant involvement of SRSF1 in OS development suggests its potential utility as a therapeutic target in OS.

Deciphering the role of alternative splicing in neoplastic diseases for immune-oncological therapies.

Alternative splicing (AS) is an important molecular biological mechanism regulated by complex mechanisms involving a plethora of cis and trans-acting elements. Furthermore, AS is tissue specific and altered in various pathologies, including infectious, inflammatory, and neoplastic diseases. Recently developed immuno-oncological therapies include monoclonal antibodies (mAbs) and chimeric antigen receptor (CAR) T cells targeting, among others, immune checkpoint (ICP) molecules. Despite therapeutic successes have been demonstrated, only a limited number of patients showed long-term benefit from these therapies with tumor entity-related differential response rates were observed. Interestingly, splice variants of common immunotherapeutic targets generated by AS are able to completely escape and/or reduce the efficacy of mAb- and/or CAR-based tumor immunotherapies. Therefore, the analyses of splicing patterns of targeted molecules in tumor specimens prior to therapy might help correct stratification, thereby increasing therapy success by antibody panel selection and antibody dosages. In addition, the expression of certain splicing factors has been linked with the patients' outcome, thereby highlighting their putative prognostic potential. Outstanding questions are addressed to translate the findings into clinical application. This review article provides an overview of the role of AS in (tumor) diseases, its molecular mechanisms, clinical relevance, and therapy response.

Temporal landscape and translational regulation of A-to-I RNA editing in mouse retina development.

BACKGROUND: The significance of A-to-I RNA editing in nervous system development is widely recognized; however, its influence on retina development remains to be thoroughly understood. RESULTS: In this study, we performed RNA sequencing and ribosome profiling experiments on developing mouse retinas to characterize the temporal landscape of A-to-I editing. Our findings revealed temporal changes in A-to-I editing, with distinct editing patterns observed across different developmental stages. Further analysis showed the interplay between A-to-I editing and alternative splicing, with A-to-I editing influencing splicing efficiency and the quantity of splicing events. A-to-I editing held the potential to enhance translation diversity, but this came at the expense of reduced translational efficiency. When coupled with splicing, it could produce a coordinated effect on gene translation. CONCLUSIONS: Overall, this study presents a temporally resolved atlas of A-to-I editing, connecting its changes with the impact on alternative splicing and gene translation in retina development.

Role of small nucleolar RNAs in alternative splicing of the doublesex gene in the silkworm, Bombyx mori.

More and more evidence shows that small noncoding RNAs (ncRNAs) play diverse roles in development, stress response and other cellular processes, but functional study of intermediate-size ncRNAs is still rare. Here, the expression profile of 16 intermediate-size ncRNAs in ovary and testis of silkworm Bombyx mori were analyzed. Twelve ncRNAs, including 5 small nucleolar RNAs (snoRNAs) and 7 unclassified ncRNAs, accumulated more in the testis than in the ovary of silkworm, especially Bm-163, Bm-51 and Bm-68. Four ncRNAs (including three orphan snoRNAs and one unclassified ncRNA) had higher expression level in the ovary than in the testis, especially Bm-86. Overexpression of the testis-enriched snoRNA Bm-68 in the female led to the accumulation of male-specific isoform of doublesex (BmdsxM) and increased the expression ratio of BmdsxM: BmdsxF. While overexpression of ovary-enriched snoRNA Bm-86 in the male decreased the expression ratio of BmdsxM: BmdsxF, indicating the roles of the two snoRNAs played in the alternative splicing of Bmdsx of silkworm, which will provide new clues for the functional study of snoRNAs in insects.

Alternative splicing of a chromatin modifier alters the transcriptional regulatory programs of stem cell maintenance and neuronal differentiation.

Development of embryonic stem cells (ESCs) into neurons requires intricate regulation of transcription, splicing, and translation, but how these processes interconnect is not understood. We found that polypyrimidine tract binding protein 1 (PTBP1) controls splicing of DPF2, a subunit of BRG1/BRM-associated factor (BAF) chromatin remodeling complexes. Dpf2 exon 7 splicing is inhibited by PTBP1 to produce the DPF2-S isoform early in development. During neuronal differentiation, loss of PTBP1 allows exon 7 inclusion and DPF2-L expression. Different cellular phenotypes and gene expression programs were induced by these alternative DPF2 isoforms. We identified chromatin binding sites enriched for each DPF2 isoform, as well as sites bound by both. In ESC, DPF2-S preferential sites were bound by pluripotency factors. In neuronal progenitors, DPF2-S sites were bound by nuclear factor I (NFI), while DPF2-L sites were bound by CCCTC-binding factor (CTCF). DPF2-S sites exhibited enhancer modifications, while DPF2-L sites showed promoter modifications. Thus, alternative splicing redirects BAF complex targeting to impact chromatin organization during neuronal development.

Discovery of immunotherapy targets for pediatric solid and brain tumors by exon-level expression.

Immunotherapy with chimeric antigen receptor T cells for pediatric solid and brain tumors is constrained by available targetable antigens. Cancer-specific exons present a promising reservoir of targets; however, these have not been explored and validated systematically in a pan-cancer fashion. To identify cancer specific exon targets, here we analyze 1532 RNA-seq datasets from 16 types of pediatric solid and brain tumors for comparison with normal tissues using a newly developed workflow. We find 2933 exons in 157 genes encoding proteins of the surfaceome or matrisome with high cancer specificity either at the gene (n = 148) or the alternatively spliced isoform (n = 9) level. Expression of selected alternatively spliced targets, including the EDB domain of fibronectin 1, and gene targets, such as COL11A1, are validated in pediatric patient derived xenograft tumors. We generate T cells expressing chimeric antigen receptors specific for the EDB domain or COL11A1 and demonstrate that these have antitumor activity. The full target list, explorable via an interactive web portal ( https://cseminer.stjude.org/ ), provides a rich resource for developing immunotherapy of pediatric solid and brain tumors using gene or AS targets with high expression specificity in cancer.

Transcriptomic analysis of differentially alternative splicing patterns in mice with inflammatory and neuropathic pain.

Although the molecular mechanisms of chronic pain have been extensively studied, a global picture of alternatively spliced genes and events in the peripheral and central nervous systems of chronic pain is poorly understood. The current study analyzed the changing pattern of alternative splicing (AS) in mouse brain, dorsal root ganglion, and spinal cord tissue under inflammatory and neuropathic pain. In total, we identified 6495 differentially alternatively spliced (DAS) genes. The molecular functions of shared DAS genes between these two models are mainly enriched in calcium signaling pathways, synapse organization, axon regeneration, and neurodegeneration disease. Additionally, we identified 509 DAS in differentially expressed genes (DEGs) shared by these two models, accounting for a small proportion of total DEGs. Our findings supported the hypothesis that the AS has an independent regulation pattern different from transcriptional regulation. Taken together, these findings indicate that AS is one of the important molecular mechanisms of chronic pain in mammals. This study presents a global description of AS profile changes in the full path of neuropathic and inflammatory pain models, providing new insights into the underlying mechanisms of chronic pain and guiding genomic clinical diagnosis methods and rational medication.

The six brain-specific TAU isoforms and their role in Alzheimer's disease and related neurodegenerative dementia syndromes.

INTRODUCTION: Alternative splicing of the human MAPT gene generates six brain-specific TAU isoforms. Imbalances in the TAU isoform ratio can lead to neurodegenerative diseases, underscoring the need for precise control over TAU isoform balance. Tauopathies, characterized by intracellular aggregates of hyperphosphorylated TAU, exhibit extensive neurodegeneration and can be classified by the TAU isoforms present in pathological accumulations. METHODS: A comprehensive review of TAU and related dementia syndromes literature was conducted using PubMed, Google Scholar, and preprint server. RESULTS: While TAU is recognized as key driver of neurodegeneration in specific tauopathies, the contribution of the isoforms to neuronal function and disease development remains largely elusive. DISCUSSION: In this review we describe the role of TAU isoforms in health and disease, and stress the importance of comprehending and studying TAU isoforms in both, physiological and pathological context, in order to develop targeted therapeutic interventions for TAU-associated diseases. HIGHLIGHTS: MAPT splicing is tightly regulated during neuronal maturation and throughout life. TAU isoform expression is development-, cell-type and brain region specific. The contribution of TAU to neurodegeneration might be isoform-specific. Ineffective TAU-based therapies highlight the need for specific targeting strategies.

The landscape of abiotic and biotic stress-responsive splice variants with deep RNA-seq datasets in hot pepper.

Alternative splicing (AS) is a widely observed phenomenon in eukaryotes that plays a critical role in development and stress responses. In plants, the large number of RNA-seq datasets in response to different environmental stressors can provide clues for identification of condition-specific and/or common AS variants for preferred agronomic traits. We report RNA-seq datasets (350.7 Gb) from Capsicum annuum inoculated with one of three bacteria, one virus, or one oomycete and obtained additional existing transcriptome datasets. In this study, we investigated the landscape of AS in response to environmental stressors, signaling molecules, and tissues from 425 total samples comprising 841.49 Gb. In addition, we identified genes that undergo AS under specific and shared stress conditions to obtain potential genes that may be involved in enhancing tolerance to stressors. We uncovered 1,642,007 AS events and identified 4,354 differential alternative splicing genes related to environmental stressors, tissues, and signaling molecules. This information and approach provide useful data for basic-research focused on enhancing tolerance to environmental stressors in hot pepper or establishing breeding programs.

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

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

The phosphatase inhibitor LB-100 creates neoantigens in colon cancer cells through perturbation of mRNA splicing.

Perturbation of protein phosphorylation represents an attractive approach to cancer treatment. Besides kinase inhibitors, protein phosphatase inhibitors have been shown to have anti-cancer activity. A prime example is the small molecule LB-100, an inhibitor of protein phosphatases 2A/5 (PP2A/PP5), enzymes that affect cellular physiology. LB-100 has proven effective in pre-clinical models in combination with immunotherapy, but the molecular underpinnings of this synergy remain understood poorly. We report here a sensitivity of the mRNA splicing machinery to phosphorylation changes in response to LB-100 in colorectal adenocarcinoma. We observe enrichment for differentially phosphorylated sites within cancer-critical splicing nodes of U2 snRNP, SRSF and hnRNP proteins. Altered phosphorylation endows LB-100-treated colorectal adenocarcinoma cells with differential splicing patterns. In PP2A-inhibited cells, over 1000 events of exon skipping and intron retention affect regulators of genomic integrity. Finally, we show that LB-100-evoked alternative splicing leads to neoantigens that are presented by MHC class 1 at the cell surface. Our findings provide a potential explanation for the pre-clinical and clinical observations that LB-100 sensitizes cancer cells to immune checkpoint blockade.

Alternative splicing analysis of lignocellulose-degrading enzyme genes and enzyme variants in Aspergillus niger.

Alternative splicing (AS) greatly expands the protein diversity in eukaryotes. Although AS variants have been frequently reported existing in filamentous fungi, it remains unclear whether lignocellulose-degrading enzyme genes in industrially important fungi undergo AS events. In this work, AS events of lignocellulose-degrading enzymes genes in Aspergillus niger under two carbon sources (glucose and wheat straw) were investigated by RNA-Seq. The results showed that a total of 23 out of the 56 lignocellulose-degrading enzyme genes had AS events and intron retention was the main type of these AS events. The AS variant enzymes from the annotated endo-ß-1,4-xylanase F1 gene (xynF1) and the endo-ß-1,4-glucanase D gene (eglD), noted as XYNF1-AS and EGLD-AS, were characterized compared to their normal splicing products XYNF1 and EGLD, respectively. The AS variant XYNF1-AS displayed xylanase activity whereas XYNF1 did not. As for EGLD-AS and EGLD, neither of them showed annotated endo-ß-1,4-glucanase activity. Instead, both showed lytic polysaccharide monooxygenase (LPMO) activity with some differences in catalytic properties. Our work demonstrated that the AS variants in A. niger were good sources for discovering novel lignocellulose-degrading enzymes. KEY POINTS: • AS events were identified in the lignocellulose-degrading enzyme genes of A. niger. • New ß-1,4-xylanase and LPMO derived from AS events were characterized.

Alternative splicing for tuneable expression of protein subunits at desired ratios.

The controlled expression of two or more proteins at a defined and stable ratio remains a substantial challenge, particularly in the bi- and multispecific antibody field. Achieving an optimal ratio of protein subunits can facilitate the assembly of multimeric proteins with high efficiency and minimize the production of by-products. In this study, we propose a solution based on alternative splicing, enabling the expression of a tunable and predefined ratio of two distinct polypeptide chains from the same pre-mRNA under the control of a single promoter. The pre-mRNA used in this study contains two open reading frames situated on separate exons. The first exon is flanked by two copies of the chicken troponin intron 4 (cTNT-I4) and is susceptible to excision from the pre-mRNA by means of alternative splicing. This specific design enables the modulation of the splice ratio by adjusting the strength of the splice acceptor. To illustrate this approach, we developed constructs expressing varying ratios of GFP and dsRED and extended their application to multimeric proteins such as monoclonal antibodies, achieving industrially relevant expression levels (>1 g/L) in a 14-day fed-batch process. The stability of the splice ratio was confirmed by droplet digital PCR in a stable pool cultivated over a 28-day period, while product quality was assessed via intact mass analysis, demonstrating absence of product-related impurities resulting from undesired splice events. Furthermore, we showcased the versatility of the construct by expressing two subunits of a bispecific antibody of the BEAT® type, which contains three distinct subunits in total.

[SRSF2 promotes glioblastoma cell proliferation by inducing alternative splicing of FSP1 and inhibiting ferroptosis].

Objective: To investigate the effect of serine/arginine-rich splicing factor 2 (SRSF2) on ferroptosis and its possible mechanism in glioblastoma cells. Methods: The online database of gene expression profiling interactive analysis 2 (GEPIA 2) and Chinese Glioma Genome Atlas were used to analyze the expression of SRSF2 in glioblastoma tissue and its association with patients prognosis. To validate the findings of the online databases, the pathological sections of glioblastoma and non-tumor brain tissues from Tianjin Medical University General Hospital, Tianjin, China were collected and analyzed by using immunohistochemistry. Silencing SRSF2 gene expression in glioblastoma cells by siRNA was analyzed with Western blot. The proliferation index was detected by using CCK8 assay. The rescued experiment was conducted by using expression plasmid of pcDNA3.1(+)-SRSF2. The activity of ferroptosis was assessed by using the levels of iron ions and malondialdehyde in glioblastoma cells and the changes in the ratio of glutathione to oxidized glutathione. The changes of gene expression and differential pre-mRNA alternative splicing (PMAS) induced by SRSF2 were monitored by using the third-generation sequencing technology analysis, namely Oxford nanopore technologies (ONT) sequencing analysis. Results: SRSF2 expression was higher in glioblastoma tissues than non-tumor brain tissues. Immunohistochemistry also showed a positive rate of 88.48%±4.60% in glioblastoma tissue which was much higher than the 9.97%±4.57% in non-tumor brain tissue. The expression of SRSF2 was inversely correlated with overall and disease-free disease survivals (P<0.01). The proliferation index of glioblastoma cells was significantly reduced by silencing with SRSF2 siRNA (P<0.01) and could be reversed with transfection of exogenous SRSF2. The levels of intracellulariron ions and malondialdehyde increased (P<0.05), but the glutathione/oxidized glutathione ratio and the expression of key proteins in the glutathione pathway remained unchanged (P>0.05). ONT sequencing results showed that silencing SRSF2 in glioblastoma cells could induce a significant alternative 3' splice site change on ferroptosis suppressor protein 1 (FSP1). Conclusion: SRSF2 inhibits the ferroptosis in glioblastoma cells and promotes their proliferation, which may be achieved by regulating FSP1 PMAS.

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.