ORF Capture-Seq as a versatile method for targeted identification of full-length isoforms

Cell Type- and Age-Specific Expression of lncRNAs across Kidney Cell Types

Key Points

Background

Accumulated evidence demonstrates that long noncoding RNAs (lncRNAs) regulate cell differentiation and homeostasis, influencing kidney aging and disease. Despite their versatility, the function of lncRNA remains poorly understood because of the lack of a reference map of lncRNA transcriptome in various cell types.

Methods

In this study, we used a targeted single-cell RNA sequencing method to enrich and characterize lncRNAs in individual cells. We applied this method to various mouse tissues, including normal and aged kidneys.

Results

Through tissue-specific clustering analysis, we identified cell type–specific lncRNAs that showed a high correlation with known cell-type marker genes. Furthermore, we constructed gene regulatory networks to explore the functional roles of differentially expressed lncRNAs in each cell type. In the kidney, we observed dynamic expression changes of lncRNAs during aging, with specific changes in glomerular cells. These cell type– and age-specific expression patterns of lncRNAs suggest that lncRNAs may have a potential role in regulating cellular processes, such as immune response and energy metabolism, during kidney aging.

Conclusions

Our study sheds light on the comprehensive landscape of lncRNA expression and function and provides a valuable resource for future analysis of lncRNAs (https://gist-fgl.github.io/sc-lncrna-atlas/).

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

Chikungunya virus infection disrupts lymph node lymphatic endothelial cell composition and function via MARCO

Infection with chikungunya virus (CHIKV) causes disruption of draining lymph node (dLN) organization, including paracortical relocalization of B cells, loss of the B cell-T cell border, and lymphocyte depletion that is associated with infiltration of the LN with inflammatory myeloid cells. Here, we found that, during the first 24 hours of infection, CHIKV RNA accumulated in MARCO-expressing lymphatic endothelial cells (LECs) in both the floor and medullary LN sinuses. The accumulation of viral RNA in the LN was associated with a switch to an antiviral and inflammatory gene expression program across LN stromal cells, and this inflammatory response - including recruitment of myeloid cells to the LN - was accelerated by CHIKV-MARCO interactions. As CHIKV infection progressed, both floor and medullary LECs diminished in number, suggesting further functional impairment of the LN by infection. Consistent with this idea, antigen acquisition by LECs, a key function of LN LECs during infection and immunization, was reduced during pathogenic CHIKV infection.

Chikungunya virus infection disrupts lymph node lymphatic endothelial cell composition and function via MARCO

Infection with chikungunya virus (CHIKV) causes disruption of draining lymph node (dLN) organization, including paracortical relocalization of B cells, loss of the B cell-T cell border, and lymphocyte depletion that is associated with infiltration of the LN with inflammatory myeloid cells. Here, we find that during the first 24 h of infection, CHIKV RNA accumulates in MARCO-expressing lymphatic endothelial cells (LECs) in both the floor and medullary LN sinuses. The accumulation of viral RNA in the LN was associated with a switch to an antiviral and inflammatory gene expression program across LN stromal cells, and this inflammatory response, including recruitment of myeloid cells to the LN, was accelerated by CHIKV-MARCO interactions. As CHIKV infection progressed, both floor and medullary LECs diminished in number, suggesting further functional impairment of the LN by infection. Consistent with this idea, we find that antigen acquisition by LECs, a key function of LN LECs during infection and immunization, was reduced during pathogenic CHIKV infection.

Coordination of alternative splicing and alternative polyadenylation revealed by targeted long read sequencing

Nervous system development is associated with extensive regulation of alternative splicing (AS) and alternative polyadenylation (APA). AS and APA have been extensively studied in isolation, but little is known about how these processes are coordinated. Here, the coordination of cassette exon (CE) splicing and APA in Drosophila was investigated using a targeted long-read sequencing approach we call Pull-a-Long-Seq (PL-Seq). This cost-effective method uses cDNA pulldown and Nanopore sequencing combined with an analysis pipeline to quantify inclusion of alternative exons in connection with alternative 3' ends. Using PL-Seq, we identified genes that exhibit significant differences in CE splicing depending on connectivity to short versus long 3'UTRs. Genomic long 3'UTR deletion was found to alter upstream CE splicing in short 3'UTR isoforms and ELAV loss differentially affected CE splicing depending on connectivity to alternative 3'UTRs. This work highlights the importance of considering connectivity to alternative 3'UTRs when monitoring AS events.

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.

Complexities of mammalian transcriptome revealed by targeted RNA enrichment techniques

Studies using highly sensitive targeted RNA enrichment methods have shown that a large portion of the human transcriptome remains to be discovered and that most of the genome is transcribed in a complex, interleaved fashion characterized by a complex web of transcripts emanating from protein coding and noncoding loci. These results resonate with those from single-cell transcriptome profiling endeavors that reveal the existence of multiple novel, cell type-specific transcripts and clearly demonstrate that our understanding of the complexities of the human transcriptome is far from being complete. Here, we review the current status of the targeted RNA enrichment techniques, their application to the discovery of novel cell type-specific transcripts, and their impact on our understanding of the human genome and transcriptome.

Structure and Alternative Splicing of the Antisense FMR1 (ASFMR1) Gene

Fragile X-associated tremor/ataxia syndrome (FXTAS) is a neurodegenerative disorder caused by an expansion of 55-200 CGG repeats (premutation) in the 5'-UTR of the FMR1 gene. Bidirectional transcription at FMR1 locus has been demonstrated and specific alternative splicing of the Antisense FMR1 (ASFMR1) gene has been proposed to have a contributing role in the pathogenesis of FXTAS. The structure of ASFMR1 gene is still uncharacterized and it is currently unknown how many isoforms of the gene are expressed and at what level in premutation carriers (PM) and if they may contribute to the premutation pathology. In this study, we characterized the ASFMR1 gene structure and the transcriptional landscape by using PacBio SMRT sequencing with target enrichment (IDT customized probe panel). We identified 45 ASFMR1 isoforms ranging in sizes from 523 bp to 6 Kb, spanning approximately 59 kb of genomic DNA. Multiplexing and sequencing of six human brain samples from PM samples and normal control (HC) were carried out on the PacBio Sequel platform. We validated the presence of these isoforms by qRT-PCR and Sanger sequencing and characterized the acceptor and donor splicing site consensus sequences. Consistent with previous studies conducted in other tissue types, we found a high expression of ASFMR1 isoform Iso131bp in brain samples of PM as compared to HC, while no differences in expression levels were observed for the newly identified isoforms IsoAS1 and IsoAS2. We investigated the role of the splicing regulatory protein Sam68 which we did not observe in the alternative splicing of the ASFMR1 gene. Our study provides a useful insight into the structure of ASFMR1 gene and transcriptional landscape along with the expression pattern of various newly identified novel isoforms and on their potential role in premutation pathology.

Coordination of Alternative Splicing and Alternative Polyadenylation revealed by Targeted Long-Read Sequencing

Nervous system development is associated with extensive regulation of alternative splicing (AS) and alternative polyadenylation (APA). AS and APA have been extensively studied in isolation, but little is known about how these processes are coordinated. Here, the coordination of cassette exon (CE) splicing and APA in Drosophila was investigated using a targeted long-read sequencing approach we call Pull-a-Long-Seq (PL-Seq). This cost-effective method uses cDNA pulldown and Nanopore sequencing combined with an analysis pipeline to resolve the connectivity of alternative exons to alternative 3' ends. Using PL-Seq, we identified genes that exhibit significant differences in CE splicing depending on connectivity to short versus long 3'UTRs. Genomic long 3'UTR deletion was found to alter upstream CE splicing in short 3'UTR isoforms and ELAV loss differentially affected CE splicing depending on connectivity to alternative 3'UTRs. This work highlights the importance of considering connectivity to alternative 3'UTRs when monitoring AS events.

Bridging the splicing gap in human genetics with long-read RNA sequencing: finding the protein isoform drivers of disease

Aberrant splicing underlies many human diseases, including cancer, cardiovascular diseases and neurological disorders. Genome-wide mapping of splicing quantitative trait loci (sQTLs) has shown that genetic regulation of alternative splicing is widespread. However, identification of the corresponding isoform or protein products associated with disease-associated sQTLs is challenging with short-read RNA-seq, which cannot precisely characterize full-length transcript isoforms. Furthermore, contemporary sQTL interpretation often relies on reference transcript annotations, which are incomplete. Solutions to these issues may be found through integration of newly emerging long-read sequencing technologies. Long-read sequencing offers the capability to sequence full-length mRNA transcripts and, in some cases, to link sQTLs to transcript isoforms containing disease-relevant protein alterations. Here, we provide an overview of sQTL mapping approaches, the use of long-read sequencing to characterize sQTL effects on isoforms, the linkage of RNA isoforms to protein-level functions and comment on future directions in the field. Based on recent progress, long-read RNA sequencing promises to be part of the human disease genetics toolkit to discover and treat protein isoforms causing rare and complex diseases.

A roadmap for the functional annotation of protein families: a community perspective

Over the last 25 years, biology has entered the genomic era and is becoming a science of ‘big data’. Most interpretations of genomic analyses rely on accurate functional annotations of the proteins encoded by more than 500 000 genomes sequenced to date. By different estimates, only half the predicted sequenced proteins carry an accurate functional annotation, and this percentage varies drastically between different organismal lineages. Such a large gap in knowledge hampers all aspects of biological enterprise and, thereby, is standing in the way of genomic biology reaching its full potential. A brainstorming meeting to address this issue funded by the National Science Foundation was held during 3–4 February 2022. Bringing together data scientists, biocurators, computational biologists and experimentalists within the same venue allowed for a comprehensive assessment of the current state of functional annotations of protein families. Further, major issues that were obstructing the field were identified and discussed, which ultimately allowed for the proposal of solutions on how to move forward.

Enhanced protein isoform characterization through long-read proteogenomics

BACKGROUND

The detection of physiologically relevant protein isoforms encoded by the human genome is critical to biomedicine. Mass spectrometry (MS)-based proteomics is the preeminent method for protein detection, but isoform-resolved proteomic analysis relies on accurate reference databases that match the sample; neither a subset nor a superset database is ideal. Long-read RNA sequencing (e.g., PacBio or Oxford Nanopore) provides full-length transcripts which can be used to predict full-length protein isoforms.

RESULTS

We describe here a long-read proteogenomics approach for integrating sample-matched long-read RNA-seq and MS-based proteomics data to enhance isoform characterization. We introduce a classification scheme for protein isoforms, discover novel protein isoforms, and present the first protein inference algorithm for the direct incorporation of long-read transcriptome data to enable detection of protein isoforms previously intractable to MS-based detection. We have released an open-source Nextflow pipeline that integrates long-read sequencing in a proteomic workflow for isoform-resolved analysis.

CONCLUSIONS

Our work suggests that the incorporation of long-read sequencing and proteomic data can facilitate improved characterization of human protein isoform diversity. Our first-generation pipeline provides a strong foundation for future development of long-read proteogenomics and its adoption for both basic and translational research.

A comprehensive long-read isoform analysis platform and sequencing resource for breast cancer

Tumors display widespread transcriptome alterations, but the full repertoire of isoform-level alternative splicing in cancer is unknown. We developed a long-read (LR) RNA sequencing and analytical platform that identifies and annotates full-length isoforms and infers tumor-specific splicing events. Application of this platform to breast cancer samples identifies thousands of previously unannotated isoforms; ~30% affect protein coding exons and are predicted to alter protein localization and function. We performed extensive cross-validation with -omics datasets to support transcription and translation of novel isoforms. We identified 3059 breast tumor–specific splicing events, including 35 that are significantly associated with patient survival. Of these, 21 are absent from GENCODE and 10 are enriched in specific breast cancer subtypes. Together, our results demonstrate the complexity, cancer subtype specificity, and clinical relevance of previously unidentified isoforms and splicing events in breast cancer that are only annotatable by LR-seq and provide a rich resource of immuno-oncology therapeutic targets.

CD8+ T-cell Immune Surveillance against a Tumor Antigen Encoded by the Oncogenic Long Noncoding RNA PVT1

CD8⁺ T cells recognize peptides displayed by HLA class I molecules on cell surfaces, monitoring pathologic conditions such as cancer. Advances in proteogenomic analysis of HLA ligandomes have demonstrated that cells present a subset of cryptic peptides derived from noncoding regions of the genome; however, the roles of cryptic HLA ligands in tumor immunity remain unknown. In the current study, we comprehensively and quantitatively investigated the HLA class I ligandome of a set of human colorectal cancer and matched normal tissues, showing that cryptic translation products accounted for approximately 5% of the HLA class I ligandome. We also found that a peptide encoded by the long noncoding RNA (lncRNA) PVT1 was predominantly enriched in multiple colorectal cancer tissues. The PVT1 gene is located downstream of the MYC gene in the genome and is aberrantly overexpressed across a variety of cancers, reflecting its oncogenic property. The PVT1 peptide was recognized by patient CD8⁺ tumor-infiltrating lymphocytes, as well as peripheral blood mononuclear cells, suggesting the presence of patient immune surveillance. Our findings show that peptides can be translated from lncRNAs and presented by HLA class I and that cancer patient T cells are capable of sensing aberrations in noncoding regions of the genome.

Perspectives and Benefits of High-Throughput Long-Read Sequencing in Microbial Ecology

Short-read, high-throughput sequencing (HTS) methods have yielded numerous important insights into microbial ecology and function. Yet, in many instances short-read HTS techniques are suboptimal, for example, by providing insufficient phylogenetic resolution or low integrity of assembled genomes. Single-molecule and synthetic long-read (SLR) HTS methods have successfully ameliorated these limitations. In addition, nanopore sequencing has generated a number of unique analysis opportunities, such as rapid molecular diagnostics and direct RNA sequencing, and both Pacific Biosciences (PacBio) and nanopore sequencing support detection of epigenetic modifications. Although initially suffering from relatively low sequence quality, recent advances have greatly improved the accuracy of long-read sequencing technologies. In spite of great technological progress in recent years, the long-read HTS methods (PacBio and nanopore sequencing) are still relatively costly, require large amounts of high-quality starting material, and commonly need specific solutions in various analysis steps. Despite these challenges, long-read sequencing technologies offer high-quality, cutting-edge alternatives for testing hypotheses about microbiome structure and functioning as well as assembly of eukaryote genomes from complex environmental DNA samples.

Isoform Age - Splice Isoform Profiling Using Long-Read Technologies

Alternative splicing (AS) of RNA is a key mechanism that results in the expression of multiple transcript isoforms from single genes and leads to an increase in the complexity of both the transcriptome and proteome. Regulation of AS is critical for the correct functioning of many biological pathways, while disruption of AS can be directly pathogenic in diseases such as cancer or cause risk for complex disorders. Current short-read sequencing technologies achieve high read depth but are limited in their ability to resolve complex isoforms. In this review we examine how long-read sequencing (LRS) technologies can address this challenge by covering the entire RNA sequence in a single read and thereby distinguish isoform changes that could impact RNA regulation or protein function. Coupling LRS with technologies such as single cell sequencing, targeted sequencing and spatial transcriptomics is producing a rapidly expanding suite of technological approaches to profile alternative splicing at the isoform level with unprecedented detail. In addition, integrating LRS with genotype now allows the impact of genetic variation on isoform expression to be determined. Recent results demonstrate the potential of these techniques to elucidate the landscape of splicing, including in tissues such as the brain where AS is particularly prevalent. Finally, we also discuss how AS can impact protein function, potentially leading to novel therapeutic targets for a range of diseases.

Best practices on the differential expression analysis of multi-species RNA-seq

Advances in transcriptome sequencing allow for simultaneous interrogation of differentially expressed genes from multiple species originating from a single RNA sample, termed dual or multi-species transcriptomics. Compared to single-species differential expression analysis, the design of multi-species differential expression experiments must account for the relative abundances of each organism of interest within the sample, often requiring enrichment methods and yielding differences in total read counts across samples. The analysis of multi-species transcriptomics datasets requires modifications to the alignment, quantification, and downstream analysis steps compared to the single-species analysis pipelines. We describe best practices for multi-species transcriptomics and differential gene expression.

Top-ranked expressed gene transcripts of human protein-coding genes investigated with GTEx dataset

With considerable accumulation of RNA-Seq transcriptome data, we have extended our understanding about protein-coding gene transcript compositions. However, alternatively compounded patterns of human protein-coding gene transcripts would complicate gene expression data processing and interpretation. It is essential to exhaustively interrogate complex mRNA isoforms of protein-coding genes with an unified data resource. In order to investigate representative mRNA transcript isoforms to be utilized as transcriptome analysis references, we utilized GTEx data to establish a top-ranked transcript isoform expression data resource for human protein-coding genes. Distinctive tissue specific expression profiles and modulations could be observed for individual top-ranked transcripts of protein-coding genes. Protein-coding transcripts or genes do occupy much higher expression fraction in transcriptome data. In addition, top-ranked transcripts are the dominantly expressed ones in various normal tissues. Intriguingly, some of the top-ranked transcripts are noncoding splicing isoforms, which imply diverse gene regulation mechanisms. Comprehensive investigation on the tissue expression patterns of top-ranked transcript isoforms is crucial. Thus, we established a web tool to examine top-ranked transcript isoforms in various human normal tissue types, which provides concise transcript information and easy-to-use graphical user interfaces. Investigation of top-ranked transcript isoforms would contribute understanding on the functional significance of distinctive alternatively spliced transcript isoforms.