Isoform-resolved transcriptome of the human preimplantation embryo

Integrating Iso-seq and RNA-seq data for the reannotation of the greater amberjack genome

The greater amberjack is a very important fishery species with high commercial value, and it is distributed worldwide. Transcriptome-based studies on S. dumerili have been limited by an inadequate reference genome and a lack of well-annotated full-length transcripts. In this study, a total of 12 tissues from juvenile and adult fish both sexes were collected for next-generation RNA sequencing (RNA-seq) and full-length isoform sequencing (Iso-seq). For Iso-seq, a total of 163,218, 149,716, and 189,169 high-quality unique transcript sequences were obtained, with an N50 of 5,441, 5,255, and 5,939, from juvenile, adult male and adult female S. dumerili, respectively. We integrated the Iso-seq and RNA-seq data to construct a comprehensive gene annotation and systematically profiled the dynamics of gene expression across the 12 tissues. Our gene models had greater detail and accuracy than those from NCBI and Ensembl, with more precise polyA locations. These resources serve as a foundation for functional genomic studies and provide valuable insights into the molecular mechanisms underlying the development, reproduction and commercial traits of amberjack.

Long-read RNA sequencing identifies region- and sex-specific C57BL/6J mouse brain mRNA isoform expression and usage

Alternative splicing (AS) contributes to the biological heterogeneity between species, sexes, tissues, and cell types. Many diseases are either caused by alterations in AS or by alterations to AS. Therefore, measuring AS accurately and efficiently is critical for assessing molecular phenotypes, including those associated with disease. Long-read sequencing enables more accurate quantification of differentially spliced isoform expression than short-read sequencing approaches, and third-generation platforms facilitate high-throughput experiments. To assess differences in AS across the cerebellum, cortex, hippocampus, and striatum by sex, we generated and analyzed Oxford Nanopore Technologies (ONT) long-read RNA sequencing (lrRNA-Seq) C57BL/6J mouse brain cDNA libraries. From > 85 million reads that passed quality control metrics, we calculated differential gene expression (DGE), differential transcript expression (DTE), and differential transcript usage (DTU) across brain regions and by sex. We found significant DGE, DTE, and DTU across brain regions and that the cerebellum had the most differences compared to the other three regions. Additionally, we found region-specific differential splicing between sexes, with the most sex differences in DTU in the cortex and no DTU in the hippocampus. We also report on two distinct patterns of sex DTU we observed, sex-divergent and sex-specific, that could potentially help explain sex differences in the prevalence and prognosis of various neurological and psychiatric disorders in future studies. Finally, we built a Shiny web application for researchers to explore the data further. Our study provides a resource for the community; it underscores the importance of AS in biological heterogeneity and the utility of long-read sequencing to better understand AS in the brain.

Zygotic Splicing Activation of the Transcriptome is a Crucial Aspect of Maternal-to-Zygotic Transition and Required for the Conversion from Totipotency to Pluripotency

During maternal-to-zygotic transition (MZT) in the embryo, mRNA undergoes complex post-transcriptional regulatory processes. However, it is unclear whether and how alternative splicing plays a functional role in MZT. By analyzing transcriptome changes in mouse and human early embryos, dynamic changes in alternative splicing during MZT are observed and a previously unnoticed process of zygotic splicing activation (ZSA) following embryonic transcriptional activation is described. As the underlying mechanism of RNA splicing, splicing factors undergo dramatic maternal-to-zygotic conversion. This conversion relies on the key maternal factors BTG4 and PABPN1L and is zygotic-transcription-dependent. CDK11-dependent phosphorylation of the key splicing factor, SF3B1, and its aggregation with SRSF2 in the subnuclear domains of 2-cell embryos are prerequisites for ZSA. Isoforms generated by erroneous splicing, such as full-length Dppa4, hinder normal embryonic development. Moreover, alternative splicing regulates the conversion of early embryonic blastomeres from totipotency to pluripotency, thereby affecting embryonic lineage differentiation. ZSA is an essential post-transcriptional process of MZT and has physiological significance in generating new life. In addition to transcriptional activation, appropriate expression of transcript isoforms is also necessary for preimplantation embryonic development.

Long-read single-cell RNA sequencing enables the study of cancer subclone-specific genotype and phenotype in chronic lymphocytic leukemia

Bruton's tyrosine kinase (BTK) inhibitors are effective for the treatment of chronic lymphocytic leukemia (CLL) due to BTK's role in B cell survival and proliferation. Treatment resistance is most commonly caused by the emergence of the hallmark BTKC481S mutation that inhibits drug binding. In this study, we aimed to investigate whether the presence of additional CLL driver mutations in cancer subclones harboring a BTKC481S mutation accelerates subclone expansion. In addition, we sought to determine whether BTK-mutated subclones exhibit distinct transcriptomic behavior when compared to other cancer subclones. To achieve these goals, we employ our recently published method (Qiao et al. 2024) that combines bulk DNA sequencing and single-cell RNA sequencing (scRNA-seq) data to genotype individual cells for the presence or absence of subclone-defining mutations. While the most common approach for scRNA-seq includes short-read sequencing, transcript coverage is limited due to the vast majority of the reads being concentrated at the priming end of the transcript. Here, we utilized MAS-seq, a long-read scRNAseq technology, to substantially increase transcript coverage across the entire length of the transcripts and expand the set of informative mutations to link cells to cancer subclones in six CLL patients who acquired BTKC481S mutations during BTK inhibitor treatment. We found that BTK-mutated subclones often acquire additional mutations in CLL driver genes, leading to faster subclone proliferation. When examining subclone-specific gene expression, we found that in one patient, BTK-mutated subclones are transcriptionally distinct from the rest of the malignant B cell population with an overexpression of CLL-relevant genes.

Transcriptome Responses to Different Salinity Conditions in Litoditis marina, Revealed by Long-Read Sequencing

The marine nematode Litoditis marina is widely distributed in intertidal zones around the globe, yet the mechanisms underlying its broad adaptation to salinity remain elusive. In this study, we applied ONT long-read sequencing technology to unravel the transcriptome responses to different salinity conditions in L. marina. Through ONT sequencing under 3‱, 30‱ and 60‱ salinity environments, we obtained 131.78 G clean data and 26,647 non-redundant long-read transcripts, including 6464 novel transcripts. The DEGs obtained from the current ONT lrRNA-seq were highly correlated with those identified in our previously reported Illumina short-read RNA sequencing data. When we compared the 30‱ to the 3‱ salinity condition, we found that GO terms such as oxidoreductase activity, cation transmembrane transport and ion transmembrane transport were shared between the ONT lrRNA-seq and Illumina data. Similarly, GO terms including extracellular space, structural constituents of cuticle, substrate-specific channel activity, ion transport and substrate-specific transmembrane transporter activity were shared between the ONT and Illumina data under 60‱ compared to 30‱ salinity. In addition, we found that 79 genes significantly increased, while 119 genes significantly decreased, as the salinity increased. Furthermore, through the GO enrichment analysis of 214 genes containing DAS, in 30‱ compared to 3‱ salinity, we found that GO terms such as cellular component assembly and coenzyme biosynthetic process were enriched. Additionally, we observed that GO terms such as cellular component assembly and coenzyme biosynthetic process were also enriched in 60‱ compared to 30‱ salinity. Moreover, we found that 86, 125, and 81 genes that contained DAS were also DEGs, in comparisons between 30‱ and 3‱, 60‱ and 30‱, and 60‱ and 3‱ salinity, respectively. In addition, we demonstrated the landscape of alternative polyadenylation in marine nematode under different salinity conditions This report provides several novel insights for the further study of the mechanisms by which euryhalinity formed and evolved, and it might also contribute to the investigation of salinity dynamics induced by global climate change.

Long-read RNA sequencing identifies region- and sex-specific C57BL/6J mouse brain mRNA isoform expression and usage

Alternative splicing (AS) contributes to the biological heterogeneity between species, sexes, tissues, and cell types. Many diseases are either caused by alterations in AS or by alterations to AS. Therefore, measuring AS accurately and efficiently is critical for assessing molecular phenotypes, including those associated with disease. Long-read sequencing enables more accurate quantification of differentially spliced isoform expression than short-read sequencing approaches, and third-generation platforms facilitate high-throughput experiments. To assess differences in AS across the cerebellum, cortex, hippocampus, and striatum by sex, we generated and analyzed Oxford Nanopore Technologies (ONT) long-read RNA sequencing (lrRNA-Seq) C57BL/6J mouse brain cDNA libraries. From >85 million reads that passed quality control metrics, we calculated differential gene expression (DGE), differential transcript expression (DTE), and differential transcript usage (DTU) across brain regions and by sex. We found significant DGE, DTE, and DTU across brain regions and that the cerebellum had the most differences compared to the other three regions. Additionally, we found region-specific differential splicing between sexes, with the most sex differences in DTU in the cortex and no DTU in the hippocampus. We also report on two distinct patterns of sex DTU we observed, sex-divergent and sex-specific, that could potentially help explain sex differences in the prevalence and prognosis of various neurological and psychiatric disorders in future studies. Finally, we built a Shiny web application for researchers to explore the data further. Our study provides a resource for the community; it underscores the importance of AS in biological heterogeneity and the utility of long-read sequencing to better understand AS in the brain.