Transcript Isoforms Research Articles

Identification of cell-type specific alternative transcripts in the multicellular alga Volvox carteri

BackgroundCell type specialization is a hallmark of complex multicellular organisms and is usually established through implementation of cell-type-specific gene expression programs. The multicellular green alga Volvox carteri has just two cell types, germ and soma, that have previously been shown to have very different transcriptome compositions which match their specialized roles. Here we interrogated another potential mechanism for differentiation in V. carteri, cell type specific alternative transcript isoforms (CTSAI).MethodsWe used pre-existing predictions of alternative transcripts and de novo transcript assembly with HISAT2 and Ballgown software to compile a list of loci with two or more transcript isoforms, identified a small subset that were candidates for CTSAI, and manually curated this subset of genes to remove false positives. We experimentally verified three candidates using semi-quantitative RT-PCR to assess relative isoform abundance in each cell type.ResultsOf the 1978 loci with two or more predicted transcript isoforms 67 of these also showed cell type isoform expression biases. After curation 15 strong candidates for CTSAI were identified, three of which were experimentally verified, and their predicted gene product functions were evaluated in light of potential cell type specific roles. A comparison of genes with predicted alternative splicing from Chlamydomonas reinhardtii, a unicellular relative of V. carteri, identified little overlap between ortholog pairs with alternative splicing in both species. Finally, we interrogated cell type expression patterns of 126 V. carteri predicted RNA binding protein (RBP) encoding genes and found 40 that showed either somatic or germ cell expression bias. These RBPs are potential mediators of CTSAI in V. carteri and suggest possible pre-adaptation for cell type specific RNA processing and a potential path for generating CTSAI in the early ancestors of metazoans and plants.ConclusionsWe predicted numerous instances of alternative transcript isoforms in Volvox, only a small subset of which showed cell type specific isoform expression bias. However, the validated examples of CTSAI supported existing hypotheses about cell type specialization in V. carteri, and also suggested new hypotheses about mechanisms of functional specialization for their gene products. Our data imply that CTSAI operates as a minor but important component of V. carteri cellular differentiation and could be used as a model for how alternative isoforms emerge and co-evolve with cell type specialization.

Principal Component Analysis of Alternative Splicing Profiles Revealed by Long-Read ONT Sequencing in Human Liver Tissue and Hepatocyte-Derived HepG2 and Huh7 Cell Lines.

The long-read RNA sequencing developed by Oxford Nanopore Technology provides a direct quantification of transcript isoforms. That makes the number of transcript isoforms per gene an intrinsically suitable metric for alternative splicing (AS) profiling in the application to this particular type of RNA sequencing. By using this simple metric and recruiting principal component analysis (PCA) as a tool to visualize the high-dimensional transcriptomic data, we were able to group biospecimens of normal human liver tissue and hepatocyte-derived malignant HepG2 and Huh7 cells into clear clusters in a 2D space. For the transcriptome-wide analysis, the clustering was observed regardless whether all genes were included in analysis or only those expressed in all biospecimens tested. However, in the application to a particular set of genes known as pharmacogenes, which are involved in drug metabolism, the clustering worsened dramatically in the latter case. Based on PCA data, the subsets of genes most contributing to biospecimens' grouping into clusters were selected and subjected to gene ontology analysis that allowed us to determine the top 20 biological processes among which translation and processes related to its regulation dominate. The suggested metrics can be a useful addition to the existing metrics for describing AS profiles, especially in application to transcriptome studies with long-read sequencing.

Long read sequencing to reveal the full complexity of a plant transcriptome by targeting both standard and long workflows

BackgroundLong read sequencing allows the analysis of full-length transcripts in plants without the challenges of reliable transcriptome assembly. Long read sequencing of transcripts from plant genomes has often utilized sized transcript libraries. However, the value of including libraries of differing sizes has not been established.MethodsA comprehensive transcriptome of the leaves of Jojoba (Simmondsia chinensis) was generated from two different PacBio library preparations: standard workflow (SW) and long workflow (LW).ResultsThe importance of using both transcript groups in the analysis was demonstrated by the high proportion of unique sequences (74.6%) that were not shared between the groups. A total of 37.8% longer transcripts were only detected in the long dataset. The completeness of the combined transcriptome was indicated by the presence of 98.7% of genes predicted in the jojoba male reference genome. The high coverage of the transcriptome was further confirmed by BUSCO analysis showing the presence of 96.9% of the genes from the core viridiplantae_odb10 lineage. The high-quality isoforms post Cd-Hit merged dataset of the two workflows had a total of 167,866 isoforms. Most of the transcript isoforms were protein-coding sequences (71.7%) containing open reading frames (ORFs) ≥ 100 amino acids (aa). Alternative splicing and intron retention were the basis of most transcript diversity when analysed at the whole genome level and by specific analysis of the apetala2 gene families.ConclusionThis suggests the need to specifically target the capture of longer transcripts to provide more comprehensive genome coverage in plant transcriptome analysis and reveal the high level of alternative splicing.

ScTEA-db: a comprehensive database of novel terminal exon isoforms identified from human single cell transcriptomes.

The usage of alternative terminal exons results in messenger RNA (mRNA) isoforms that differ in their 3' untranslated regions (3' UTRs) and often also in their protein-coding sequences. Alternative 3' UTRs contain different sets of cis-regulatory elements known to regulate mRNA stability, translation and localization, all of which are vital to cell identity and function. In previous work, we revealed that ∼25 percent of the experimentally observed RNA 3' ends are located within regions currently annotated as intronic, indicating that many 3' end isoforms remain to be uncovered. Also, the inclusion of not yet annotated terminal exons is more tissue specific compared to the already annotated ones. Here, we present the single cell-based Terminal Exon Annotation database (scTEA-db, www.scTEA-db.org) that provides the community with 12063 so far not yet annotated terminal exons and associated transcript isoforms identified by analysing 53069 publicly available single cell transcriptomes. Our scTEA-db web portal offers an array of features to find and explore novel terminal exons belonging to 5538 human genes, 110 of which are known cancer drivers. In summary, scTEA-db provides the foundation for studying the biological role of large numbers of so far not annotated terminal exon isoforms in cell identity and function.

CanISO: a database of genomic and transcriptomic variations in domestic dog (Canis lupus familiaris)

BackgroundThe domestic dog, Canis lupus familiaris, is a companion animal for humans as well as an animal model in cancer research due to similar spontaneous occurrence of cancers as humans. Despite the social and biological importance of dogs, the catalogue of genomic variations and transcripts for dogs is relatively incomplete.ResultsWe developed CanISO, a new database to hold a large collection of transcriptome profiles and genomic variations for domestic dogs. CanISO provides 87,692 novel transcript isoforms and 60,992 known isoforms from whole transcriptome sequencing of canine tumors (N = 157) and their matched normal tissues (N = 64). CanISO also provides genomic variation information for 210,444 unique germline single nucleotide polymorphisms (SNPs) from the whole exome sequencing of 183 dogs, with a query system that searches gene- and transcript-level information as well as covered SNPs. Transcriptome profiles can be compared with corresponding human transcript isoforms at a tissue level, or between sample groups to identify tumor-specific gene expression and alternative splicing patterns.ConclusionsCanISO is expected to increase understanding of the dog genome and transcriptome, as well as its functional associations with humans, such as shared/distinct mechanisms of cancer. CanISO is publicly available at https://www.kobic.re.kr/caniso/.

Regulation of alternative splicing and polyadenylation in neurons.

Cell-type-specific gene expression is a fundamental feature of multicellular organisms and is achieved by combinations of regulatory strategies. Although cell-restricted transcription is perhaps the most widely studied mechanism, co-transcriptional and post-transcriptional processes are also central to the spatiotemporal control of gene functions. One general category of expression control involves the generation of multiple transcript isoforms from an individual gene, whose balance and cell specificity are frequently tightly regulated via diverse strategies. The nervous system makes particularly extensive use of cell-specific isoforms, specializing the neural function of genes that are expressed more broadly. Here, we review regulatory strategies and RNA-binding proteins that direct neural-specific isoform processing. These include various classes of alternative splicing and alternative polyadenylation events, both of which broadly diversify the neural transcriptome. Importantly, global alterations of splicing and alternative polyadenylation are characteristic of many neural pathologies, and recent genetic studies demonstrate how misregulation of individual neural isoforms can directly cause mutant phenotypes.

Gene Expression Profiling Reveals Potential Players of Sex Determination and Asymmetrical Development in Chicken Embryo Gonads.

Despite the notable progress made in recent years, the understanding of the genetic control of gonadal sex differentiation and asymmetrical ovariogenesis in chicken during embryonic development remains incomplete. This study aimed to identify potential key genes and speculate about the mechanisms associated with ovary and testis development via an analysis of the results of PacBio and Illumina transcriptome sequencing of embryonic chicken gonads at the initiation of sexual differentiation (E4.5, E5.5, and E6.5). PacBio sequencing detected 328 and 233 significantly up-regulated transcript isoforms in females and males at E4.5, respectively. Illumina sequencing detected 95, 296 and 445 DEGs at E4.5, E5.5, and E6.5, respectively. Moreover, both sexes showed asymmetrical expression in gonads, and more DEGs were detected on the left side. There were 12 DEGs involved in cell proliferation shared between males and females in the left gonads. GO analysis suggested that coagulation pathways may be involved in the degradation of the right gonad in females and that blood oxygen transport pathways may be involved in preventing the degradation of the right gonad in males. These results provide a comprehensive expression profile of chicken embryo gonads at the initiation of sexual differentiation, which can serve as a theoretical basis for further understanding the mechanism of bird sex determination and its evolutionary process.

Profiling the polyadenylated transcriptome of extracellular vesicles with long-read nanopore sequencing

BackgroundWhile numerous studies have described the transcriptomes of extracellular vesicles (EVs) in different cellular contexts, these efforts have typically relied on sequencing methods requiring RNA fragmentation, which limits interpretations on the integrity and isoform diversity of EV-targeted RNA populations. It has been assumed that mRNA signatures in EVs are likely to be fragmentation products of the cellular mRNA material, and the extent to which full-length mRNAs are present within EVs remains to be clarified.ResultsUsing long-read nanopore RNA sequencing, we sought to characterize the full-length polyadenylated (poly-A) transcriptome of EVs released by human chronic myelogenous leukemia K562 cells. We detected 443 and 280 RNAs that were respectively enriched or depleted in EVs. EV-enriched poly-A transcripts consist of a variety of biotypes, including mRNAs, long non-coding RNAs, and pseudogenes. Our analysis revealed that 10.58% of all EV reads, and 18.67% of all cellular (WC) reads, corresponded to known full-length transcripts, with mRNAs representing the largest biotype for each group (EV = 58.13%, WC = 43.93%). We also observed that for many well-represented coding and non-coding genes, diverse full-length transcript isoforms were present in EV specimens, and these isoforms were reflective-of but often in different ratio compared to cellular samples.ConclusionThis work provides novel insights into the compositional diversity of poly-A transcript isoforms enriched within EVs, while also underscoring the potential usefulness of nanopore sequencing to interrogate secreted RNA transcriptomes.

Large-scale RNAseq analysis provide a new insight into the critical genes and regulatory networks of tiller development mediated by gibberellin in sugarcane

Tillering plays a key role in sugarcane growth, development, and cane yield. Gibberellin (GA) and its biosynthesis inhibitor, paclobutrazol, have been used to manipulate the tiller production in sugarcane planting. However, the key genes and regulatory networks of sugarcane tiller development remain unclear. In this study, we applied GA3 and paclobutrazol to soak the sugarcane seeds for generating samples with varying tiller numbers. A total of 54 samples from multiple tissues of three tillering stages were employed for transcriptome sequencing using two full-length transcript isoform sequences datasets as reference. We observed significant downregulation of two genes within the GA metabolism pathway in response to the GA synthesis inhibitor. Additionally, we noted a substantial enrichment of genes associated with photosynthesis following treatment with GA3 or paclobutrazol. Using all the differentially expressed genes (DEGs) between each sample group, 41 gene modules were generated in the coexpression network. The associations between modules and traits were performed. The Gene Ontology (GO) enrichment, KEGG pathway enrichment, and hub genes analysis were implemented in the modules correlated with traits, especially, the modules correlated with GA3 content and tiller number were emphasized. Notably, the pathway of ubiquitin-mediated proteolysis were enriched in both the GA3 content positive correlated module and the tiller number negative correlated module, while an E3 ubiquitin ligase gene, RZFP34, was the top hub gene in these two modules, indicating the vital role of ubiquitin-mediated proteasomal degradation system in tillering regulation in sugarcane. Our study provides a comprehensive insight into regulatory networks involved in GA-mediated regulation of tiller development in sugarcane, which would lay the foundation of hormone regulation to promote sugarcane production.

The sequencing of the key genes and end products in the TLR4 signaling pathway from the kidney of Rana dybowskii exposed to Aeromonas hydrophila.

Infectious diseases caused by Aeromonas hydrophila (AH) have reduced the populations of Rana dybowskii). However, little is known about the immune response of R. dybowskii against AH infections. The toll-like receptor (TLR) signaling pathway has been identified as a critical component in innate immunity, responsible for identifying pathogen-associated molecular patterns in pathogens. Our study used the next-generation sequencing technique and single-molecule long-read sequencing to determine the structures of transcript isoforms and functions of genes in the kidneys of R. dybowskii, as well as identify and validate the related genes in the TLR4 signaling pathway. In total, 628,774 reads of inserts were identified, including 300,053 full-length non-chimeric reads and 233,592 non-full-length reads. Among the transcriptome sequences, 124 genes were identified as homologs of known genes in the TLR4 pathway especially inflammatory cytokines and receptors. Our findings shed light on the structures and functions of R. dybowskii genes exposed to AH and confirm the presence of both MyD88-dependent and independent pathways in R. dybowskii. Our work reveals how various functional proteins in amphibians at the initial stage of immune response are activated and complete their corresponding functions in a short time.

A new mechanism of flowering regulation by the competition of isoforms in Osmanthus fragrans.

The regulation of flowering time is typically governed by transcription factors or epigenetic modifications. Transcript isoforms can play important roles in flowering regulation. Recently, transcript isoforms were discovered in the key genes, OfAP1 and OfTFL1, of the flowering regulatory network in Osmanthus fragrans. OfAP1-b generates a full-length isoform of OfAP1-b1 as well as an isoform of OfAP1-b2 that lacks the C-terminal domain. Although OfAP1-b2 does not possess an activation domain, it has a complete K domain that allows it to form heterodimers. OfAP1-b2 competes with OfAP1-b1 by binding with OfAGL24 to create non-functional and functional heterodimers. As a result, OfAP1-b1 promotes flowering while OfAP1-b2 delays flowering. OfTFL1 produces two isoforms located in different areas: OfTFL1-1 in the cytoplasm and OfTFL1-2 in the nucleus. When combined with OfFD, OfTFL1-1 does not enter the nucleus to repress AP1 expression, leading to early flowering. Conversely, when combined with OfFD, OfTFL1-2 enters the nucleus to repress AP1 expression, resulting in later flowering. Tissue-specific expression and functional conservation testing of OfAP1 and OfTFL1 support the new model's effectiveness in regulating flowering. Overall, this study provides new insights into regulating flowering time by the competition of isoforms.

Genome-Wide Analysis of Alternative Splicing Events Responding to High Temperatures in Populus tomentosa Carr.

Through alternative splicing (AS) processes, eukaryotic genes can generate a variety of transcription isoforms that lower the expression levels of the normal transcripts or result in diversity in the genes’ activities. Then, AS plays a significant role in the control of plant development and stress tolerance. In this study, we analyzed Populus tomentosa Carr. TC1521’s AS episodes in response to high temperatures. The samples treated at 25 °C, 30 °C, 35 °C, and 40 °C produced a total of 10,418, 11,202, 9947, and 14,121 AS events, respectively, which responded to 4105, 4276, 4079, and 4915 genes, respectively, representing 9.84%, 10.25%, 9.78%, and 11.78% of the total number of transcribed genes, respectively. The most common AS pattern, accounting for 42.31% to 51.00% of all AS events, was intron retention (IR), followed by exon skipping (ES), which accounted for 9.14% to 10.23% of all AS events. respectively. According to sequence characterization, AS was negatively correlated with guanine-cytosine content (GC content) but favorably correlated with intron length, exon number, exon length, and gene transcription level. Compared to treatment at 25 °C, 2001 distinct AS genes were discovered at 40 °C. They were primarily enriched in the RNA degradation pathway and the valine, leucine, and isoleucine degradation route, according to (gene ontology) GO and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. These findings demonstrated how the AS process might be severely impacted by high temperatures. In addition, the information on AS isoforms helped us comprehend stress-resistance mechanisms in new ways and completed molecular design breeding.

Linc2function: A Comprehensive Pipeline and Webserver for Long Non-Coding RNA (lncRNA) Identification and Functional Predictions Using Deep Learning Approaches.

Long non-coding RNAs (lncRNAs), comprising a significant portion of the human transcriptome, serve as vital regulators of cellular processes and potential disease biomarkers. However, the function of most lncRNAs remains unknown, and furthermore, existing approaches have focused on gene-level investigation. Our work emphasizes the importance of transcript-level annotation to uncover the roles of specific transcript isoforms. We propose that understanding the mechanisms of lncRNA in pathological processes requires solving their structural motifs and interactomes. A complete lncRNA annotation first involves discriminating them from their coding counterparts and then predicting their functional motifs and target bio-molecules. Current in silico methods mainly perform primary-sequence-based discrimination using a reference model, limiting their comprehensiveness and generalizability. We demonstrate that integrating secondary structure and interactome information, in addition to using transcript sequence, enables a comprehensive functional annotation. Annotating lncRNA for newly sequenced species is challenging due to inconsistencies in functional annotations, specialized computational techniques, limited accessibility to source code, and the shortcomings of reference-based methods for cross-species predictions. To address these challenges, we developed a pipeline for identifying and annotating transcript sequences at the isoform level. We demonstrate the effectiveness of the pipeline by comprehensively annotating the lncRNA associated with two specific disease groups. The source code of our pipeline is available under the MIT licensefor local use by researchers to make new predictions using the pre-trained models or to re-train models on new sequence datasets. Non-technical users can access the pipeline through a web server setup.

P02.13.A THIRD-GENERATION SEQUENCING OF BRAIN TUMORS ON SINGLE-CELL LEVEL

Abstract BACKGROUND Intratumor heterogeneity (ITH) in brain tumors has been suggested to increase tumor adaptability and thus drive tumor progression and treatment resistance. Single cell sequencing has widened the range of possibilities in detecting clonal diversity and tumor evolution. As of now, most single cell experiments rely on next-generation sequencing (NGS) platforms. However, they detect only short RNA/DNA fragments that do not show the full picture of existing transcript isoforms or genetic variants. New, long-read sequencing methods such as Oxford Nanopore Technologies (ONT) or PacBio sequencing may now enable novel approaches to single cell experiments. We here set-out to complement existing single-cell analysis workflows and data pipelines with long-read sequencing approaches. We compare and integrate the data output in order to advance the understanding of ITH in brain tumors, paradigmatically demonstrated on meningioma. MATERIAL AND METHODS We set up a pipeline for single-cell/nuclei RNA long-read sequencing on ONT devices. First, libraries as derived from the 10x genomics workflow were employed, subsequently followed by independently generated libraries. To widen the insight into ITH on multiple levels, we next assessed single-cell/nuclei DNA for small somatic variants and methylation along a modified version of this platform. RESULTS Long-read RNA-sequencing yielded results comparable to short-read sequencing. Depending on sample quality, higher reads per cell and library were required for ONT sequencing to achieve similar discrimination of clusters by expression. On the other hand, cell-type specific isoforms could be identified. Specific small somatic variants known from bulk sequencing data were identified when the region was sufficiently covered in the RNA data, which was naturally more often achieved in long read than short-read data. Thus, cell-specific variant calling was enabled in a high-throughput pipeline. Maintaining the methylation information in DNA sequencing could so far only be set-up in a plate-based approach. CONCLUSION Using these approaches, we resolved features of ITH and clonal tumor evolution, aiding us in gaining a broader understanding of brain tumor biology. Furthermore, it showed novel approaches that could possibly be implemented in routine diagnostic procedures in the future. These novel techniques could be applied to other tumor entities and possibly further omics layers such as single cell proteomics and metabolomics could be added to gain a greater understanding of different tumor types and their clonal diversity.

Cross-species analysis of differential transcript usage in humans and chickens with fatty liver disease.

Fatty liver disease is a common condition, characterized by excess fat accumulation in the liver. It can contribute to more severe liver-related health issues, making it a critical concern in avian and human medicine. Apart from modifying the gene expression of liver cells, the disease also alters the expression of specific transcript isoforms, which might serve as new biological markers for both species. This study aimed to identify cross-species genes displaying differential expressions in their transcript isoforms in humans and chickens with fatty liver disease. We performed differential gene expression and differential transcript usage (DTU) analyses on messenger RNA datasets from the livers of both chickens and humans with fatty liver disease. Using appropriate cross-species gene identification methods, we reviewed the acquired candidate genes and their transcript isoforms to determine their potential role in fatty liver disease's pathogenesis. We identified seven genes - ALG5, BRD7, DIABLO, RSU1, SFXN5, STIMATE, TJP3, and VDAC2 - and their corresponding transcript isoforms as potential candidates (false discovery rate ≤0.05). Our findings showed that these genes most likely contribute to fatty disease development and progression. This study successfully identified novel human-chicken DTU genes in fatty liver disease. Further research is encouraged to verify the functions and regulations of these transcript isoforms as potential diagnostic markers for fatty liver disease in humans and chickens.

Regulation of epithelial cell differentiation by the Ubiquitous expressed transcript isoform 1 in ulcerative colitis.

Mucosal healing has emerged as a desirable treatment goal for patients with ulcerative colitis (UC). Healing of mucosal wounds involves epithelial cell proliferation and differentiation, and Y-box transcription factor ZONAB has recently been identified as the key modulator of intestinal epithelial restitution. We studied the characteristics of UXT-V1 expression in UC patients using immunohistochemistry and qPCR. The functional role of UXT-V1 in the colonic epithelium was investigated using lentivirus-mediated shRNA in vitro and ex vivo. Through endogenous Co-immunoprecipitation and LC-MS/MS, we identified ZONAB as a UXT-V1-interactive protein. Herein, we report that UXT-V1 promotes differentiation of intestinal epithelial cells by regulating the nuclear translocation of ZONAB. UXT-V1 was upregulated in the intestinal epithelia of UC patients compared with that of healthy controls. Knocking down UXT-V1 in NCM-460 cells led to the enrichment of pathways associated with proliferation and differentiation. Furthermore, the absence of UXT-V1 in cultured intestinal epithelial cells and colonic organoids inhibited differentiation to the goblet cell phenotype. Mechanistically, the loss of UXT-V1 in the intestinal epithelial cells allowed nuclear translocation of ZONAB, wherein it regulated the transcription of differentiation-related genes, including AML1 and KLF4. Taken together, our study reveals a potential role of UXT-V1 in regulating epithelial cell differentiation, proving a molecular basis for mucosal healing in UC.

Allele-specific quantification of human leukocyte antigen transcript isoforms by nanopore sequencing.

While tens of thousands of HLA alleles have been identified by DNA sequencing, the contribution of alternative splicing to HLA diversity is not well characterized. In this study, we sought to determine if long-read sequencing could be used to accurately quantify allele-specific HLA transcripts in primary human lymphocytes. cDNA libraries were prepared from peripheral blood lymphocytes from 12 donors and sequenced by nanopore long-read sequencing. HLA reads were aligned to donor-specific reference sequences based on the known type of each donor. Allele-specific exon utilization was calculated as the proportion of reads aligning to each allele containing known exons, and transcript isoforms were quantified based on patterns of exon utilization within individual reads. Splice variants were rare among class I HLA genes (median exon retention rate 99%-100%), except for several HLA-C alleles with exon 5 spliced out of up to 15% of reads. Splice variants were also rare among class II HLA genes (median exon retention rate 98%-100%), except for HLA-DQB1. Consistent with previous work, exon 5 of HLA-DQB1 was spliced out in alleles with a mutated splice acceptor site at rs28688207. Surprisingly, a 28% loss of exon 5 was also observed in HLA-DQB1 alleles with an intact splice acceptor site at rs28688207. We describe a simple bioinformatic workflow to quantify allele-specific expression of HLA transcript isoforms. Further studies are warranted to characterize the repertoire of HLA transcripts expressed in different cell types and tissues across diverse populations.

Hepatitis B virus serum RNA transcript isoform composition and proportion in chronic hepatitis B patients by nanopore long-read sequencing.

Serum hepatitis B virus (HBV) RNA is a promising new biomarker to manage and predict clinical outcomes of chronic hepatitis B (CHB) infection. However, the HBV serum transcriptome within encapsidated particles, which is the biomarker analyte measured in serum, remains poorly characterized. This study aimed to evaluate serum HBV RNA transcript composition and proportionality by PCR-cDNA nanopore sequencing of samples from CHB patients having varied HBV genotype (gt, A to F) and HBeAg status. Longitudinal specimens from 3 individuals during and following pregnancy (approximately 7 months between time points) were also investigated. HBV RNA extracted from 16 serum samples obtained from 13 patients (73.3% female, 84.6% Asian) was sequenced and serum HBV RNA isoform detection and quantification were performed using three bioinformatic workflows; FLAIR, RATTLE, and a GraphMap-based workflow within the Galaxy application. A spike-in RNA variant (SIRV) control mix was used to assess run quality and coverage. The proportionality of transcript isoforms was based on total HBV reads determined by each workflow. All chosen isoform detection workflows showed high agreement in transcript proportionality and composition for most samples. HBV pregenomic RNA (pgRNA) was the most frequently observed transcript isoform (93.8% of patient samples), while other detected transcripts included pgRNA spliced variants, 3' truncated variants and HBx mRNA, depending on the isoform detection method. Spliced variants of pgRNA were primarily observed in HBV gtB, C, E, or F-infected patients, with the Sp1 spliced variant detected most frequently. Twelve other pgRNA spliced variant transcripts were identified, including 3 previously unidentified transcripts, although spliced isoform identification was very dependent on the workflow used to analyze sequence data. Longitudinal sampling among pregnant and post-partum antiviral-treated individuals showed increasing proportions of 3' truncated pgRNA variants over time. This study demonstrated long-read sequencing as a promising tool for the characterization of the serum HBV transcriptome. However, further studies are needed to better understand how serum HBV RNA isoform type and proportion are linked to CHB disease progression and antiviral treatment response.

TEQUILA-seq: a versatile and low-cost method for targeted long-read RNA sequencing

Long-read RNA sequencing (RNA-seq) is a powerful technology for transcriptome analysis, but the relatively low throughput of current long-read sequencing platforms limits transcript coverage. One strategy for overcoming this bottleneck is targeted long-read RNA-seq for preselected gene panels. We present TEQUILA-seq, a versatile, easy-to-implement, and low-cost method for targeted long-read RNA-seq utilizing isothermally linear-amplified capture probes. When performed on the Oxford nanopore platform with multiple gene panels of varying sizes, TEQUILA-seq consistently and substantially enriches transcript coverage while preserving transcript quantification. We profile full-length transcript isoforms of 468 actionable cancer genes across 40 representative breast cancer cell lines. We identify transcript isoforms enriched in specific subtypes and discover novel transcript isoforms in extensively studied cancer genes such as TP53. Among cancer genes, tumor suppressor genes (TSGs) are significantly enriched for aberrant transcript isoforms targeted for degradation via mRNA nonsense-mediated decay, revealing a common RNA-associated mechanism for TSG inactivation. TEQUILA-seq reduces the per-reaction cost of targeted capture by 2-3 orders of magnitude, as compared to a standard commercial solution. TEQUILA-seq can be broadly used for targeted sequencing of full-length transcripts in diverse biomedical research settings.

KITENIN promotes aerobic glycolysis through PKM2 induction by upregulating the c-Myc/hnRNPs axis in colorectal cancer

PurposeThe oncoprotein KAI1 C-terminal interacting tetraspanin (KITENIN; vang-like 1) promotes cell metastasis, invasion, and angiogenesis, resulting in shorter survival times in cancer patients. Here, we aimed to determine the effects of KITENIN on the energy metabolism of human colorectal cancer cells.Experimental designThe effects of KITENIN on energy metabolism were evaluated using in vitro assays. The GEPIA web tool was used to extrapolate the clinical relevance of KITENIN in cancer cell metabolism. The bioavailability and effect of the disintegrator of KITENIN complex compounds were evaluated by LC–MS, in vivo animal assay.ResultsKITENIN markedly upregulated the glycolytic proton efflux rate and aerobic glycolysis by increasing the expression of GLUT1, HK2, PKM2, and LDHA. β-catenin, CD44, CyclinD1 and HIF-1A, including c-Myc, were upregulated by KITENIN expression. In addition, KITENIN promoted nuclear PKM2 and PKM2-induced transactivation, which in turn, increased the expression of downstream mediators. This was found to be mediated through an effect of c-Myc on the transcription of hnRNP isoforms and a switch to the M2 isoform of pyruvate kinase, which increased aerobic glycolysis. The disintegration of KITENIN complex by silencing the KITENIN or MYO1D downregulated aerobic glycolysis. The disintegrator of KITENIN complex compound DKC1125 and its optimized form, DKC-C14S, exhibited the inhibition activity of KITENIN-mediated aerobic glycolysis in vitro and in vivo.ConclusionsThe oncoprotein KITENIN induces PKM2-mediated aerobic glycolysis by upregulating the c-Myc/hnRNPs axis.