Updated and enhanced pig cardiac transcriptome based on long-read RNA sequencing and proteomics

Transcriptomic Characterization of Genes Regulating the Stemness in Porcine Atrial Cardiomyocytes during Primary In Vitro Culture

Heart failure remains a major cause of death worldwide. There is a need to establish new management options as current treatment is frequently suboptimal. Clinical approaches based on autologous stem cell transplant is potentially a good alternative. The heart was long considered an organ unable to regenerate and renew. However, several reports imply that it may possess modest intrinsic regenerative potential. To allow for detailed characterization of cell cultures, whole transcriptome profiling was performed after 0, 7, 15, and 30 days of in vitro cell cultures (IVC) from the right atrial appendage and right atrial wall utilizing microarray technology. In total, 4239 differentially expressed genes (DEGs) with ratio > abs |2| and adjusted p-value ≤ 0.05 for the right atrial wall and 4662 DEGs for the right atrial appendage were identified. It was shown that a subset of DEGs, which have demonstrated some regulation of expression levels with the duration of the cell culture, were enriched in the following GO BP (Gene Ontology Biological Process) terms: "stem cell population maintenance" and "stem cell proliferation". The results were validated by RT-qPCR. The establishment and detailed characterization of in vitro culture of myocardial cells may be important for future applications of these cells in heart regeneration processes.

Full-Length RNA Sequencing Provides Insights into Goldfish Evolution under Artificial Selection

Goldfish Carassius auratus is an ideal model for exploring fish morphology evolution. Although genes underlying several ornamental traits have been identified, little is known about the effects of artificial selection on embryo gene expression. In the present study, hybrid transcriptome sequencing was conducted to reveal gene expression profiles of Celestial-Eye (CE) and Ryukin (RK) goldfish embryos. Full-length transcriptome sequencing on the PacBio platform identified 54,218 and 54,106 transcript isoforms in CE and RK goldfish, respectively. Of particular note was that thousands of alternative splicing (AS) and alternative polyadenylation (APA) events were identified in both goldfish breeds, and most of them were inter-breed specific. RT-PCR and Sanger sequencing showed that most of the predicted AS and APA were correct. Moreover, abundant long non-coding RNA and fusion genes were detected, and again most of them were inter-breed specific. Through RNA-seq, we detected thousands of differentially expressed genes (DEGs) in each embryonic stage between the two goldfish breeds. KEGG enrichment analysis on DEGs showed extensive differences between CE and RK goldfish in gene expression. Taken together, our results demonstrated that artificial selection has led to far-reaching influences on goldfish gene expression, which probably laid the genetic basis for hundreds of goldfish variations.

Transcriptome annotation of 17 porcine tissues using nanopore sequencing technology

The annotation of animal genomes plays an important role in elucidating molecular mechanisms behind the genetic control of economically important traits. Here, we employed long-read sequencing technology, Oxford Nanopore Technology, to annotate the pig transcriptome across 17 tissues from two Yorkshire littermate pigs. More than 9.8 million reads were obtained from a single flow cell, and 69 781 unique transcripts at 50 108 loci were identified. Of these transcripts, 16 255 were found to be novel isoforms, and 22 344 were found at loci that were novel and unannotated in the Ensembl (release 102) and NCBI (release 106) annotations. Novel transcripts were mostly expressed in cerebellum, followed by lung, liver, spleen, and hypothalamus. By comparing the unannotated transcripts to existing databases, there were 21 285 (95.3%) transcripts matched to the NT database (v5) and 13 676 (61.2%) matched to the NR database (v5). Moreover, there were 4324 (19.4%) transcripts matched to the SwissProt database (v5), corresponding to 11 356 proteins. Tissue-specific gene expression analyses showed that 9749 transcripts were highly tissue-specific, and cerebellum contained the most tissue-specific transcripts. As the same samples were used for the annotation of cis-regulatory elements in the pig genome, the transcriptome annotation generated by this study provides an additional and complementary annotation resource for the Functional Annotation of Animal Genomes effort to comprehensively annotate the pig genome.

Nematode gene annotation by machine-learning-assisted proteotranscriptomics enables proteome-wide evolutionary analysis

Nematodes encompass more than 24,000 described species, which were discovered in almost every ecological habitat, and make up >80% of metazoan taxonomic diversity in soils. The last common ancestor of nematodes is believed to date back to ∼650-750 million years, generating a large and phylogenetically diverse group to be explored. However, for most species high-quality gene annotations are incomprehensive or missing. Combining short-read RNA sequencing with mass spectrometry-based proteomics and machine-learning quality control in an approach called proteotranscriptomics, we improve gene annotations for nine genome-sequenced nematode species and provide new gene annotations for three additional species without genome assemblies. Emphasizing the sensitivity of our methodology, we provide evidence for two hitherto undescribed genes in the model organism Caenorhabditis elegans Extensive phylogenetic systems analysis using this comprehensive proteome annotation provides new insights into evolutionary processes of this metazoan group.

Proteotranscriptomics - A facilitator in omics research

Applications in omics research, such as comparative transcriptomics and proteomics, require the knowledge of the species-specific gene sequence and benefit from a comprehensive high-quality annotation of the coding genes to achieve high coverage. While protein-coding genes can in simple cases be detected by scanning the genome for open reading frames, in more complex genomes exonic sequences are separated by introns. Despite advances in sequencing technologies that allow for ever-growing numbers of genomes, the quality of many of the provided genome assemblies do not reach reference quality. These non-contiguous assemblies with gaps and the necessity to predict splice sites limit accurate gene annotation from solely genomic data. In contrast, the transcriptome only contains transcribed gene regions, is devoid of introns and thus provides the optimal basis for the identification of open reading frames. The additional integration of proteomics data to validate predicted protein-coding genes further enriches for accurate gene models. This review outlines the principles of the proteotranscriptomics approach, discusses common challenges and suggests methods for improvement.

Full-Length Spatial Transcriptomics Reveals the Unexplored Isoform Diversity of the Myocardium Post-MI

We introduce Single-cell Nanopore Spatial Transcriptomics (scNaST), a software suite to facilitate the analysis of spatial gene expression from second- and third-generation sequencing, allowing to generate a full-length near-single-cell transcriptional landscape of the tissue microenvironment. Taking advantage of the Visium Spatial platform, we adapted a strategy recently developed to assign barcodes to long-read single-cell sequencing data for spatial capture technology. Here, we demonstrate our workflow using four short axis sections of the mouse heart following myocardial infarction. We constructed a de novo transcriptome using long-read data, and successfully assigned 19,794 transcript isoforms in total, including clinically-relevant, but yet uncharacterized modes of transcription, such as intron retention or antisense overlapping transcription. We showed a higher transcriptome complexity in the healthy regions, and identified intron retention as a mode of transcription associated with the infarct area. Our data revealed a clear regional isoform switching among differentially used transcripts for genes involved in cardiac muscle contraction and tissue morphogenesis. Molecular signatures involved in cardiac remodeling integrated with morphological context may support the development of new therapeutics towards the treatment of heart failure and the reduction of cardiac complications.

Discovering functional motifs in long noncoding RNAs

Long noncoding RNAs (lncRNAs) are products of pervasive transcription that closely resemble messenger RNAs on the molecular level, yet function through largely unknown modes of action. The current model is that the function of lncRNAs often relies on specific, typically short, conserved elements, connected by linkers in which specific sequences and/or structures are less important. This notion has fueled the development of both computational and experimental methods focused on the discovery of functional elements within lncRNA genes, based on diverse signals such as evolutionary conservation, predicted structural elements, or the ability to rescue loss-of-function phenotypes. In this review, we outline the main challenges that the different methods need to overcome, describe the recently developed approaches, and discuss their respective limitations. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

Identification of Protein Isoforms Using Reference Databases Built from Long and Short Read RNA-Sequencing

Alternative splicing can lead to distinct protein isoforms. These can have different functions in specific cells and tissues or in different developmental stages. In this study, we explored whether transcripts assembled from long read, nanopore-based, direct RNA-sequencing (RNA-seq) could improve the identification of protein isoforms in human K562 cells. By comparing with Illumina-based short read RNA-seq, we showed that a large proportion of Ensembl transcripts (5949/14,326) and genes expressing alternatively spliced transcripts (486/2981) identified with long direct reads were missed by short paired-end reads. By co-analyzing proteomic and transcriptomic data, we also showed that some peptides (826/35,976), proteins (262/3215), and protein isoforms arising from distinct transcript variants (574/1212) identified with isoform-specific peptides via custom long-read-based databases were missed in Illumina-derived databases. Finally, we generated unequivocal peptide evidence for a set of protein isoforms and showed that long read, direct RNA-seq allows the discovery of novel protein isoforms not already in reference databases or custom databases built from short read RNA-seq data. Our analysis highlights the benefits of long read RNA-seq data in the generation of reference databases to increase tandem mass spectrometry (MS/MS) identification of protein isoforms.

Transcriptional atlas analysis from multiple tissues reveals the expression specificity patterns in beef cattle

Background

A comprehensive analysis of gene expression profiling across tissues can provide necessary information for an in-depth understanding of their biological functions. We performed a large-scale gene expression analysis and generated a high-resolution atlas of the transcriptome in beef cattle.

Results

Our transcriptome atlas was generated from 135 bovine tissues in adult beef cattle, covering 51 tissue types of major organ systems (e.g., muscular system, digestive system, immune system, reproductive system). Approximately 94.76% of sequencing reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We detected a total of 60,488 transcripts, and 32% of them were not reported before. We identified 2654 housekeeping genes (HKGs) and 477 tissue-specific genes (TSGs) across tissues. Using weighted gene co-expression network analysis, we obtained 24 modules with 237 hub genes (HUBGs). Functional enrichment analysis showed that HKGs mainly maintain the basic biological activities of cells, while TSGs were involved in tissue differentiation and specific physiological processes. HKGs in bovine tissues were more conserved in terms of expression pattern as compared to TSGs and HUBGs among multiple species. Finally, we obtained a subset of tissue-specific differentially expressed genes (DEGs) between beef and dairy cattle and several functional pathways, which may be involved in production and health traits.

Conclusions

We generated a large-scale gene expression atlas across the major tissues in beef cattle, providing valuable information for enhancing genome assembly and annotation. HKGs, TSGs, and HUBGs further contribute to better understanding the biology and evolution of multiple tissues in cattle. DEGs between beef and dairy cattle also fill in the knowledge gaps about differential transcriptome regulation of bovine tissues underlying economically important traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01269-4.

Full-Length Transcriptome Sequencing: An Insight Into the Dog Model of Heart Failure

Heart failure (HF) leads to a progressive increase in morbidity and mortality rates. This study aimed to explore the transcriptional landscape during HF and identify differentially expressed transcripts (DETs) and alternative splicing events associated with HF. We generated a dog model of HF (n = 3) using right ventricular pacemaker implantation. We performed full-length transcriptome sequencing (based on nanopore platform) on the myocardial tissues and analyzed the transcripts using differential expression analysis and functional annotation methods [Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses]. Additionally, we estimated the expression of the selected genes by quantitative real-time PCR (qRT-PCR) and detected the proportion of immune cells using flow cytometry. We found that increased B-type natriuretic peptide reduced ejection fraction, and apparent clinical signs were observed in the dog model of HF. We identified 67,458 transcripts using full-length transcriptome sequencing. A total of 785 DETs were obtained from the HF and control groups. These DETs were mainly enriched in the immune responses, especially Th1, Th2, and Th17 cell differentiation processes. Furthermore, flow cytometry results revealed that the proportion of Th1 and Th17 cells increased in patients with HF compared to controls, while the proportion of Th2 cells decreased. Differentially expressed genes in the HF and control groups associated with Th1, Th2, and Th17 cell differentiation were quantified using qRT-PCR. We also identified variable splicing events of sarcomere genes (e.g., MYBPC3, TNNT2, TTN, FLNC, and TTNI3). In addition, we detected 4,892 transcription factors and 406 lncRNAs associated with HF. Our analysis based on full-length transcript sequencing provided an analysis perspective in a dog model of HF, which is valuable for molecular research in an increasingly relevant large animal model of HF.

Improved nanopore direct RNA sequencing of cardiac myocyte samples by selective mt-RNA depletion

RNA sequencing is a powerful tool to analyze gene expression transcriptome wide. However, RNA sequencing in general and especially the recently developed methods of long read RNA sequencing are still low-throughput and cost-intensive. Here, one important design choice is to concentrate the sequencing capacity on specific parts of the transcriptome. Especially, abundant transcripts as ribosomal RNAs may dominate the available sequencing space, if not removed prior to sequencing. Several methods exist to reduce ribosomal RNA read numbers: either based on enrichment of the relevant fraction (polyA⁺ RNA) or depletion, respectively degradation of ribosomal RNAs. Furthermore, commercial kits are available to deplete globin transcripts from blood samples. However, so far, no solution exists to deal with other tissue-specific highly abundant transcripts. This is especially of interest in the heart and other muscle derived samples, where reads originating from mitochondrial RNAs make up to 30% of reads in polyA⁺ selected libraries and around 70% in single cell sequencing experiments. We present a simple method to diminish sequencing of mitochondrial RNAs in Oxford Nanopore direct RNA sequencing libraries by RNase H based clipping of the polyA tail. We show that mt-clipping enables enhanced detection of cytoplasmic mRNAs, among them genes involved in heart development and pathogenesis. Mt-clipping may be applied as well to other sequencing protocols that are based on oligo(dT) priming and can be easily adapted to other tissue-specific high-abundant transcripts.

Full-length transcriptome analysis reveals the mechanism of acupuncture at PC6 improves cardiac function in myocardial ischemia model

Background

The pathological process of myocardial ischemia (MI) is very complicated. Acupuncture at PC6 has been proved to be effective against MI injury, but the mechanism remains unclear. This study investigated the mechanism that underlies the effect of acupuncture on MI through full-length transcriptome.

Methods

Adult male C57/BL6 mice were randomly divided into control, MI, and PC6 groups. Mice in MI and PC6 group generated MI model by ligating the left anterior descending (LAD) coronary artery. The samples were collected 5 days after acupuncture treatment.

Results

The results showed that treatment by acupuncture improved cardiac function, decreased myocardial infraction area, and reduced the levels of cTnT and cTnI. Based on full-length transcriptome sequencing, 5083 differential expression genes (DEGs) and 324 DEGs were identified in the MI group and PC6 group, respectively. These genes regulated by acupuncture were mainly enriched in the inflammatory response pathway. Alternative splicing (AS) is a post-transcriptional action that contributes to the diversity of protein. In all samples, 8237 AS events associated with 1994 genes were found. Some differential AS-involved genes were enriched in the pathway related to heart disease. We also identified 602 new genes, 4 of which may the novel targets of acupuncture in MI.

Conclusions

Our findings suggest that the effect of acupuncture on MI may be based on the multi-level regulation of the transcriptome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13020-021-00465-8.

Large-Scale Multiplexing Permits Full-Length Transcriptome Annotation of 32 Bovine Tissues From a Single Nanopore Flow Cell

A comprehensive annotation of transcript isoforms in domesticated species is lacking. Especially considering that transcriptome complexity and splicing patterns are not well-conserved between species, this presents a substantial obstacle to genomic selection programs that seek to improve production, disease resistance, and reproduction. Recent advances in long-read sequencing technology have made it possible to directly extrapolate the structure of full-length transcripts without the need for transcript reconstruction. In this study, we demonstrate the power of long-read sequencing for transcriptome annotation by coupling Oxford Nanopore Technology (ONT) with large-scale multiplexing of 93 samples, comprising 32 tissues collected from adult male and female Hereford cattle. More than 30 million uniquely mapping full-length reads were obtained from a single ONT flow cell, and used to identify and characterize the expression dynamics of 99,044 transcript isoforms at 31,824 loci. Of these predicted transcripts, 21% exactly matched a reference transcript, and 61% were novel isoforms of reference genes, substantially increasing the ratio of transcript variants per gene, and suggesting that the complexity of the bovine transcriptome is comparable to that in humans. Over 7,000 transcript isoforms were extremely tissue-specific, and 61% of these were attributed to testis, which exhibited the most complex transcriptome of all interrogated tissues. Despite profiling over 30 tissues, transcription was only detected at about 60% of reference loci. Consequently, additional studies will be necessary to continue characterizing the bovine transcriptome in additional cell types, developmental stages, and physiological conditions. However, by here demonstrating the power of ONT sequencing coupled with large-scale multiplexing, the task of exhaustively annotating the bovine transcriptome - or any mammalian transcriptome - appears significantly more feasible.