Single-molecule, full-length transcript isoform sequencing reveals disease-associated RNA isoforms in cardiomyocytes

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.

Single-cell and spatial transcriptomics: Bridging current technologies with long-read sequencing

Single-cell technologies have transformed biomedical research over the last decade, opening up new possibilities for understanding cellular heterogeneity, both at the genomic and transcriptomic level. In addition, more recent developments of spatial transcriptomics technologies have made it possible to profile cells in their tissue context. In parallel, there have been substantial advances in sequencing technologies, and the third generation of methods are able to produce reads that are tens of kilobases long, with error rates matching the second generation short reads. Long reads technologies make it possible to better map large genome rearrangements and quantify isoform specific abundances. This further improves our ability to characterize functionally relevant heterogeneity. Here, we show how researchers have begun to combine single-cell, spatial transcriptomics, and long-read technologies, and how this is resulting in powerful new approaches to profiling both the genome and the transcriptome. We discuss the achievements so far, and we highlight remaining challenges and opportunities.

Near-gapless and haplotype-resolved apple genomes provide insights into the genetic basis of rootstock-induced dwarfing

Dwarfing rootstocks have transformed the production of cultivated apples; however, the genetic basis of rootstock-induced dwarfing remains largely unclear. We have assembled chromosome-level, near-gapless and haplotype-resolved genomes for the popular dwarfing rootstock 'M9', the semi-vigorous rootstock 'MM106' and 'Fuji', one of the most commonly grown apple cultivars. The apple orthologue of auxin response factor 3 (MdARF3) is in the Dw1 region of 'M9', the major locus for rootstock-induced dwarfing. Comparing 'M9' and 'MM106' genomes revealed a 9,723-bp allele-specific long terminal repeat retrotransposon/gypsy insertion, DwTE, located upstream of MdARF3. DwTE is cosegregated with the dwarfing trait in two segregating populations, suggesting its prospective utility in future dwarfing rootstock breeding. In addition, our pipeline discovered mobile mRNAs that may contribute to the development of dwarfed scion architecture. Our research provides valuable genomic resources and applicable methodology, which have the potential to accelerate breeding dwarfing rootstocks for apple and other perennial woody fruit trees.

Mechanisms of RBM20 Cardiomyopathy: Insights From Model Systems

RBM20 (RNA-binding motif protein 20) is a vertebrate- and muscle-specific RNA-binding protein that belongs to the serine-arginine-rich family of splicing factors. The RBM20 gene was first identified as a dilated cardiomyopathy-linked gene over a decade ago. Early studies in Rbm20 knockout rodents implicated disrupted splicing of RBM20 target genes as a causative mechanism. Clinical studies show that pathogenic variants in RBM20 are linked to aggressive dilated cardiomyopathy with early onset heart failure and high mortality. Subsequent studies employing pathogenic variant knock-in animal models revealed that variants in a specific portion of the arginine-serine-rich domain in RBM20 not only disrupt splicing but also hinder nucleocytoplasmic transport and lead to the formation of RBM20 biomolecular condensates in the sarcoplasm. Conversely, mice harboring a disease-associated variant in the RRM (RNA recognition motif) do not show evidence of adverse remodeling or exhibit sudden death despite disrupted splicing of RBM20 target genes. Thus, whether disrupted splicing, biomolecular condensates, or both contribute to dilated cardiomyopathy is under debate. Beyond this, additional questions remain, such as whether there is sexual dimorphism in the presentation of RBM20 cardiomyopathy. What are the clinical features of RBM20 cardiomyopathy and why do some individuals develop more severe disease than others? In this review, we summarize the reported observations and discuss potential mechanisms of RBM20 cardiomyopathy derived from studies employing in vivo animal models and in vitro human-induced pluripotent stem cell-derived cardiomyocytes. Potential therapeutic strategies to treat RBM20 cardiomyopathy are also discussed.

Patterns of Gene Expression, Splicing, and Allele-Specific Expression Vary among Macular Tissues and Clinical Stages of Age-Related Macular Degeneration

Age-related macular degeneration (AMD) is a leading cause of blindness, and elucidating its underlying disease mechanisms is vital to the development of appropriate therapeutics. We identified differentially expressed genes (DEGs) and differentially spliced genes (DSGs) across the clinical stages of AMD in disease-affected tissue, the macular retina pigment epithelium (RPE)/choroid and the macular neural retina within the same eye. We utilized 27 deeply phenotyped donor eyes (recovered within a 6 h postmortem interval time) from Caucasian donors (60-94 years) using a standardized published protocol. Significant findings were then validated in an independent set of well-characterized donor eyes (n = 85). There was limited overlap between DEGs and DSGs, suggesting distinct mechanisms at play in AMD pathophysiology. A greater number of previously reported AMD loci overlapped with DSGs compared to DEGs between disease states, and no DEG overlap with previously reported loci was found in the macular retina between disease states. Additionally, we explored allele-specific expression (ASE) in coding regions of previously reported AMD risk loci, uncovering a significant imbalance in C3 rs2230199 and CFH rs1061170 in the macular RPE/choroid for normal eyes and intermediate AMD (iAMD), and for CFH rs1061147 in the macular RPE/choroid for normal eyes and iAMD, and separately neovascular AMD (NEO). Only significant DEGs/DSGs from the macular RPE/choroid were found to overlap between disease states. STAT1, validated between the iAMD vs. normal comparison, and AGTPBP1, BBS5, CERKL, FGFBP2, KIFC3, RORα, and ZNF292, validated between the NEO vs. normal comparison, revealed an intricate regulatory network with transcription factors and miRNAs identifying potential upstream and downstream regulators. Findings regarding the complement genes C3 and CFH suggest that coding variants at these loci may influence AMD development via an imbalance of gene expression in a tissue-specific manner. Our study provides crucial insights into the multifaceted genomic underpinnings of AMD (i.e., tissue-specific gene expression changes, potential splice variation, and allelic imbalance), which may open new avenues for AMD diagnostics and therapies specific to iAMD and NEO.

Identification and quantification of small exon-containing isoforms in long-read RNA sequencing data

Small exons are pervasive in transcriptomes across organisms, and their quantification in RNA isoforms is crucial for understanding gene functions. Although long-read RNA-seq based on Oxford Nanopore Technologies (ONT) offers the advantage of covering transcripts in full length, its lower base accuracy poses challenges for identifying individual exons, particularly microexons (≤ 30 nucleotides). Here, we systematically assess small exons quantification in synthetic and human ONT RNA-seq datasets. We demonstrate that reads containing small exons are often not properly aligned, affecting the quantification of relevant transcripts. Thus, we develop a local-realignment method for misaligned exons (MisER), which remaps reads with misaligned exons to the transcript references. Using synthetic and simulated datasets, we demonstrate the high sensitivity and specificity of MisER for the quantification of transcripts containing small exons. Moreover, MisER enabled us to identify small exons with a higher percent spliced-in index (PSI) in neural, particularly neural-regulated microexons, when comparing 14 neural to 16 non-neural tissues in humans. Our work introduces an improved quantification method for long-read RNA-seq and especially facilitates studies using ONT long-reads to elucidate the regulation of genes involving small exons.

Mislocalization of pathogenic RBM20 variants in dilated cardiomyopathy is caused by loss-of-interaction with Transportin-3

Severe forms of dilated cardiomyopathy (DCM) are associated with point mutations in the alternative splicing regulator RBM20 that are frequently located in the arginine/serine-rich domain (RS-domain). Such mutations can cause defective splicing and cytoplasmic mislocalization, which leads to the formation of detrimental cytoplasmic granules. Successful development of personalized therapies requires identifying the direct mechanisms of pathogenic RBM20 variants. Here, we decipher the molecular mechanism of RBM20 mislocalization and its specific role in DCM pathogenesis. We demonstrate that mislocalized RBM20 RS-domain variants retain their splice regulatory activity, which reveals that aberrant cellular localization is the main driver of their pathological phenotype. A genome-wide CRISPR knockout screen combined with image-enabled cell sorting identified Transportin-3 (TNPO3) as the main nuclear importer of RBM20. We show that the direct RBM20-TNPO3 interaction involves the RS-domain, and is disrupted by pathogenic variants. Relocalization of pathogenic RBM20 variants to the nucleus restores alternative splicing and dissolves cytoplasmic granules in cell culture and animal models. These findings provide proof-of-principle for developing therapeutic strategies to restore RBM20's nuclear localization in RBM20-DCM patients.

Striated muscle-specific base editing enables correction of mutations causing dilated cardiomyopathy

Dilated cardiomyopathy is the second most common cause for heart failure with no cure except a high-risk heart transplantation. Approximately 30% of patients harbor heritable mutations which are amenable to CRISPR-based gene therapy. However, challenges related to delivery of the editing complex and off-target concerns hamper the broad applicability of CRISPR agents in the heart. We employ a combination of the viral vector AAVMYO with superior targeting specificity of heart muscle tissue and CRISPR base editors to repair patient mutations in the cardiac splice factor Rbm20, which cause aggressive dilated cardiomyopathy. Using optimized conditions, we repair >70% of cardiomyocytes in two Rbm20 knock-in mouse models that we have generated to serve as an in vivo platform of our editing strategy. Treatment of juvenile mice restores the localization defect of RBM20 in 75% of cells and splicing of RBM20 targets including TTN. Three months after injection, cardiac dilation and ejection fraction reach wild-type levels. Single-nuclei RNA sequencing uncovers restoration of the transcriptional profile across all major cardiac cell types and whole-genome sequencing reveals no evidence for aberrant off-target editing. Our study highlights the potential of base editors combined with AAVMYO to achieve gene repair for treatment of hereditary cardiac diseases.

Cardiac splicing as a diagnostic and therapeutic target

Despite advances in therapeutics for heart failure and arrhythmias, a substantial proportion of patients with cardiomyopathy do not respond to interventions, indicating a need to identify novel modifiable myocardial pathobiology. Human genetic variation associated with severe forms of cardiomyopathy and arrhythmias has highlighted the crucial role of alternative splicing in myocardial health and disease, given that it determines which mature RNA transcripts drive the mechanical, structural, signalling and metabolic properties of the heart. In this Review, we discuss how the analysis of cardiac isoform expression has been facilitated by technical advances in multiomics and long-read and single-cell sequencing technologies. The resulting insights into the regulation of alternative splicing - including the identification of cardiac splice regulators as therapeutic targets and the development of a translational pipeline to evaluate splice modulators in human engineered heart tissue, animal models and clinical trials - provide a basis for improved diagnosis and therapy. Finally, we consider how the medical and scientific communities can benefit from facilitated acquisition and interpretation of splicing data towards improved clinical decision-making and patient care.

Double-sided plasmonic metasurface for simultaneous biomolecular separation and SERS detection

Porous membrane-based nanofiltration separation of small biomolecules is a widely used biotechnology for which size-based selectivity is a critical parameter of technological relevance. Efficient determination of size selectivity calls for an advanced detection method capable of performing sensitive, rapid, and on-membrane examination. Surface-enhanced Raman spectroscopy (SERS) is such a detection method that has been widely recognized as an ultrasensitive technique for trace-level detection with sensitivity down to the single-molecule level. In this work, we for the first time develop a double-sided hierarchical porous membrane-like plasmonic metasurface to realize high-selectivity bimolecular separation and simultaneous ultrasensitive SERS detection. This highly flexible device, consisting of subwavelength nanocone pairs surrounded by randomly orientated sub-5 nm nanogrooves, was prepared by combining customized "top-down" fabrication of conical nanopores in an ion-track registered polycarbonate membrane and self-assembly of nanogrooves on the membrane surface through physical vapor deposition. The unique tip-to-tip oriented conical nanopores in the device enables excellent size-based molecular selectivity; the hierarchical groove-pore structure supports a peculiar cascaded electromagnetic near-field enhancement mechanism, endowing the device with SERS-based molecular detection of ultrahigh sensitivity, uniformity, repeatability, and polarization independence. With such dual structural merits and performance enhancement, we demonstrate effective nanofiltration separation of small-sized adenine from big-sized ss-DNA and synergistic SERS determination of their species. We experimentally demonstrate an ultrasensitive detection of 4-mercaptopyridine down to 10 pM. Together with its unparalleled mechanical flexibility, this double-side-responsive plasmonic metasurface membrane can find great potential in real-world molecular filtration and detection under extremely complex working conditions.

Single-cell and long-read sequencing to enhance modelling of splicing and cell-fate determination

Single-cell sequencing technologies have revolutionised the life sciences and biomedical research. Single-cell sequencing provides high-resolution data on cell heterogeneity, allowing high-fidelity cell type identification, and lineage tracking. Computational algorithms and mathematical models have been developed to make sense of the data, compensate for errors and simulate the biological processes, which has led to breakthroughs in our understanding of cell differentiation, cell-fate determination and tissue cell composition. The development of long-read (a.k.a. third-generation) sequencing technologies has produced powerful tools for investigating alternative splicing, isoform expression (at the RNA level), genome assembly and the detection of complex structural variants (at the DNA level). In this review, we provide an overview of the recent advancements in single-cell and long-read sequencing technologies, with a particular focus on the computational algorithms that help in correcting, analysing, and interpreting the resulting data. Additionally, we review some mathematical models that use single-cell and long-read sequencing data to study cell-fate determination and alternative splicing, respectively. Moreover, we highlight the emerging opportunities in modelling cell-fate determination that result from the combination of single-cell and long-read sequencing technologies.

Genotype Complements the Phenotype: Identification of the Pathogenicity of an LMNA Splice Variant by Nanopore Long-Read Sequencing in a Large DCM Family

Dilated cardiomyopathy (DCM) is a common cause of heart failure (HF) and is of familial origin in 20−40% of cases. Genetic testing by next-generation sequencing (NGS) has yielded a definite diagnosis in many cases; however, some remain elusive. In this study, we used a combination of NGS, human-induced pluripotent-stem-cell-derived cardiomyocytes (iPSC-CMs) and nanopore long-read sequencing to identify the causal variant in a multi-generational pedigree of DCM. A four-generation family with familial DCM was investigated. Next-generation sequencing (NGS) was performed on 22 family members. Skin biopsies from two affected family members were used to generate iPSCs, which were then differentiated into iPSC-CMs. Short-read RNA sequencing was used for the evaluation of the target gene expression, and long-read RNA nanopore sequencing was used to evaluate the relevance of the splice variants. The pedigree suggested a highly penetrant, autosomal dominant mode of inheritance. The phenotype of the family was suggestive of laminopathy, but previous genetic testing using both Sanger and panel sequencing only yielded conflicting evidence for LMNA p.R644C (rs142000963), which was not fully segregated. By re-sequencing four additional affected family members, further non-coding LMNA variants could be detected: rs149339264, rs199686967, rs201379016, and rs794728589. To explore the roles of these variants, iPSC-CMs were generated. RNA sequencing showed the LMNA expression levels to be significantly lower in the iPSC-CMs of the LMNA variant carriers. We demonstrated a dysregulated sarcomeric structure and altered calcium homeostasis in the iPSC-CMs of the LMNA variant carriers. Using targeted nanopore long-read sequencing, we revealed the biological significance of the variant c.356+1G>A, which generates a novel 5′ splice site in exon 1 of the cardiac isomer of LMNA, causing a nonsense mRNA product with almost complete RNA decay and haploinsufficiency. Using novel molecular analysis and nanopore technology, we demonstrated the pathogenesis of the rs794728589 (c.356+1G>A) splice variant in LMNA. This study highlights the importance of precise diagnostics in the clinical management and workup of cardiomyopathies.

Structural insights into crista junction formation by the Mic60-Mic19 complex

Mitochondrial cristae membranes are the oxidative phosphorylation sites in cells. Crista junctions (CJs) form the highly curved neck regions of cristae and are thought to function as selective entry gates into the cristae space. Little is known about how CJs are generated and maintained. We show that the central coiled-coil (CC) domain of the mitochondrial contact site and cristae organizing system subunit Mic60 forms an elongated, bow tie-shaped tetrameric assembly. Mic19 promotes Mic60 tetramerization via a conserved interface between the Mic60 mitofilin and Mic19 CHCH (CC-helix-CC-helix) domains. Dimerization of mitofilin domains exposes a crescent-shaped membrane-binding site with convex curvature tailored to interact with the curved CJ neck. Our study suggests that the Mic60-Mic19 subcomplex traverses CJs as a molecular strut, thereby controlling CJ architecture and function.

Full-length transcriptomic analysis in murine and human heart reveals diversity of PGC-1α promoters and isoforms regulated distinctly in myocardial ischemia and obesity

Background

Peroxisome proliferator-activated receptor gamma coactivator-1 alpha (PGC-1α) acts as a transcriptional coactivator and regulates mitochondrial function. Various isoforms are generated by alternative splicing and differentially regulated promoters. In the heart, total PGC-1α deficiency knockout leads to dilatative cardiomyopathy, but knowledge on the complexity of cardiac isoform expression of PGC-1α remains sparse. Thus, this study aims to generate a reliable dataset on cardiac isoform expression pattern by long-read mRNA sequencing, followed by investigation of differential regulation of PGC-1α isoforms under metabolic and ischemic stress, using high-fat-high-sucrose-diet-induced obesity and a murine model of myocardial infarction.

Results

Murine (C57Bl/6J) or human heart tissue (obtained during LVAD-surgery) was used for long-read mRNA sequencing, resulting in full-length transcriptomes including 58,000 mRNA isoforms with 99% sequence accuracy. Automatic bioinformatic analysis as well as manual similarity search against exonic sequences leads to identification of putative coding PGC-1α isoforms, validated by PCR and Sanger sequencing. Thereby, 12 novel transcripts generated by hitherto unknown splicing events were detected. In addition, we postulate a novel promoter with homologous and strongly conserved sequence in human heart. High-fat diet as well as ischemia/reperfusion (I/R) injury transiently reduced cardiac expression of PGC-1α isoforms, with the most pronounced effect in the infarcted area. Recovery of PGC-1α-isoform expression was even more decelerated when I/R was performed in diet-induced obese mice.

Conclusions

We deciphered for the first time a complete full-length transcriptome of the murine and human heart, identifying novel putative PGC-1α coding transcripts including a novel promoter. These transcripts are differentially regulated in I/R and obesity suggesting transcriptional regulation and alternative splicing that may modulate PGC-1α function in the injured and metabolically challenged heart.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01360-w.

Single-nuclei isoform RNA sequencing unlocks barcoded exon connectivity in frozen brain tissue

Single-nuclei RNA sequencing characterizes cell types at the gene level. However, compared to single-cell approaches, many single-nuclei cDNAs are purely intronic, lack barcodes and hinder the study of isoforms. Here we present single-nuclei isoform RNA sequencing (SnISOr-Seq). Using microfluidics, PCR-based artifact removal, target enrichment and long-read sequencing, SnISOr-Seq increased barcoded, exon-spanning long reads 7.5-fold compared to naive long-read single-nuclei sequencing. We applied SnISOr-Seq to adult human frontal cortex and found that exons associated with autism exhibit coordinated and highly cell-type-specific inclusion. We found two distinct combination patterns: those distinguishing neural cell types, enriched in TSS-exon, exon-polyadenylation-site and non-adjacent exon pairs, and those with multiple configurations within one cell type, enriched in adjacent exon pairs. Finally, we observed that human-specific exons are almost as tightly coordinated as conserved exons, implying that coordination can be rapidly established during evolution. SnISOr-Seq enables cell-type-specific long-read isoform analysis in human brain and in any frozen or hard-to-dissociate sample.

RBM20-Related Cardiomyopathy: Current Understanding and Future Options

Splice regulators play an essential role in the transcriptomic diversity of all eukaryotic cell types and organ systems. Recent evidence suggests a contribution of splice-regulatory networks in many diseases, such as cardiomyopathies. Adaptive splice regulators, such as RNA-binding motif protein 20 (RBM20) determine the physiological mRNA landscape formation, and rare variants in the RBM20 gene explain up to 6% of genetic dilated cardiomyopathy (DCM) cases. With ample knowledge from RBM20-deficient mice, rats, swine and induced pluripotent stem cells (iPSCs), the downstream targets and quantitative effects on splicing are now well-defined and the prerequisites for corrective therapeutic approaches are set. This review article highlights some of the recent advances in the field, ranging from aspects of granule formation to 3D genome architectures underlying RBM20-related cardiomyopathy. Promising therapeutic strategies are presented and put into context with the pathophysiological characteristics of RBM20-related diseases.