Skip to main content

REVIEW article

Front. Mol. Med, 22 February 2022
Sec. Molecular Medicine for Cardiology
This article is part of the Research Topic Epigenetic Mechanisms in Cardiovascular Disease View all 4 articles

Perspectives on Bulk-Tissue RNA Sequencing and Single-Cell RNA Sequencing for Cardiac Transcriptomics

  • 1Department of Molecular Genetics, Faculty of Science and Engineering, Faculty of Health, Medicine and Life Sciences, Maastricht University, Maastricht, Netherlands
  • 2Department of Biochemistry, CARIM School for Cardiovascular Diseases, Maastricht University, Maastricht, Netherlands
  • 3Department of Genetic Epidemiology, Institute of Human Genetics, University Hospital Münster, Münster, Germany

The heart has been the center of numerous transcriptomic studies in the past decade. Even though our knowledge of the key organ in our cardiovascular system has significantly increased over the last years, it is still not fully understood yet. In recent years, extensive efforts were made to understand the genetic and transcriptomic contribution to cardiac function and failure in more detail. The advent of Next Generation Sequencing (NGS) technologies has brought many discoveries but it is unable to comprehend the finely orchestrated interactions between and within the various cell types of the heart. With the emergence of single-cell sequencing more than 10 years ago, researchers gained a valuable new tool to enable the exploration of new subpopulations of cells, cell-cell interactions, and integration of multi-omic approaches at a single-cell resolution. Despite this innovation, it is essential to make an informed choice regarding the appropriate technique for transcriptomic studies, especially when working with myocardial tissue. Here, we provide a primer for researchers interested in transcriptomics using NGS technologies.

Introduction

The heart is the first fully functional organ to be formed during embryonic development. In the past decades, conventional molecular biology studies have revealed important physiological and pathological mechanisms within the cardiovascular system, however it was unable to comprehend the finely orchestrated interactions between the various cells and cell types within the heart. Thus, our understanding of cardiac cell diversity remains limited and cardiovascular diseases persist as the main cause for morbidity and mortality worldwide.

The functional phenotype of each cellular unity is largely determined by its underlying gene expression leading to the recent increase in publications addressing the cardiac transcriptome. RNA is essential to biological processes in cells and cell-cell communication, providing critical information directly associated with cell phenotypes. Consequently, the transcriptome portrays a representation of the current biological pathways and processes in the examined material. Next generation sequencing (NGS) technologies have served as powerful tools to study genomic traits and provided key insights in various research or clinical fields (Wang et al., 2009; Garraway, 2013). To date, the most commonly used technique to decipher the transcriptional landscape is high-throughput RNA sequencing (RNA-Seq), which offers a quantitative and open system for profiling transcriptional expression at genome scale and hence provides a variety of applications. Several robust RNA-Seq protocols have since been developed, each with its distinctive purpose and (dis)advantages and rendering RNA-Seq almost mandatory for every molecular biology study.

The introduction of single-cell RNA-Seq (scRNA-Seq) has revolutionized genomic research and is yet another milestone to add to the list of accomplishments that include the completion of the Human Genome Project in 2003 (International Human Genome Sequencing Consortium, 2004). Furthermore, since its invention more than 10 years ago scRNA-Seq has made tremendous progress driven by the rapid development of innovative technologies and computational analysis methods, and more importantly the in-depth knowledge of biological processes. Its biggest advantage lies in its capability to look at the transcriptome of individual cells compared to conventional bulk techniques, which measure the average gene expression across cells in a sample. It is therefore not surprising that the cardiovascular field quickly started to integrate transcriptomic techniques into their research. The most recent studies already identified a large heterogeneity among cardiac cell types and during cellular differentiation, allowing for the discovery of novel genes involved in the complex connectivity network.

Despite all these innovations, it is essential to make an informed choice regarding the appropriate technique for the study of interest, especially considering myocardial tissue. However, many researchers in the field do not yet have a complete apprehension of the various technologies and their benefits and pitfalls. Therefore, here we provide a primer for researchers interested in transcriptomics using NGS technologies.

Sequencing Methodologies

Bulk RNA Sequencing

For over a decade, researchers all over the world have used the conventional bulk sequencing methods on RNA extracted from a population of cells to study gene expression changes in different tissues (Cloonan et al., 2008; Lister et al., 2008; Mortazavi et al., 2008; Nagalakshmi et al., 2008; Wilhelm et al., 2008). Since then, the system has been optimized for different types of RNAs and qualities of starting material. In general, bulk RNA-Seq refers to every sequencing approach that relies on averaged gene expression from a population of cells to reveal RNA presence as well as quantity in a sample of cells during the time of measurement. Therefore, bulk-based approaches can identify differences between sample conditions. Bulk RNA-Seq is not particularly limited by technical applications to the heart, nevertheless there are many criteria to consider when choosing bulk sequencing to ensure high-quality data.

Sample and Library Preparation

Several steps can make a difference on data quality during bulk RNA-Seq pipelines. For instance, sample and library preparation will have a direct effect on the outcome of the analysis. Their workflow can be further subdivided into RNA isolation, RNA depletion and cDNA synthesis. Unfortunately, their single stranded nature makes RNA very unstable, and susceptible to hydrolysis and heat degradation. To ensure the optimal conditions before sequencing RNA quality must be assessed, which is commonly done using the RNA Integrity Number (RIN) with a value between 1 (low quality) and 10 (high quality) (Schroeder et al., 2006). A RIN value over six is considered good enough for sequencing (Kukurba and Montgomery, 2015), confirming the availability of high-quality RNA. Unfortunately, samples obtained from human biopsies or paraffin embedded tissues can have an adverse effect on the quality of RNA retained. Of note, even frozen RNA will lose quality over the years (Seelenfreund et al., 2014) and therefore, the RIN should always be assessed right before library preparation. In general, bulk RNA-Seq requires a minimal amount of RNA as input, but certain methodologies require more.

“Bulk” refers to the total source of RNA in a cell population allowing in depth analysis and therefore all molecules of the transcriptome can be evaluated using bulk sequencing. Interestingly, total RNA can be sequenced, or specific types of RNA can be isolated beforehand from the total RNA pool, which is composed of ribosomal RNA (rRNA), pre-mRNA and the different classes of non-coding RNA (ncRNA). Various methodologies have been developed to selectively deplete or enrich a specific type of RNA molecule before or during library preparation (Accerbi et al., 2010). It is recommended to remove rRNA transcripts before library construction due to its over presentation in the cell. Otherwise, rRNAs will overwhelm most sequencing reads and leading to an overall reduction in sequencing depth and detection of less-abundant RNAs, such as many ncRNAs (Chu and Corey, 2012). When focusing on protein-coding RNA molecules, many protocols aim to enrich for polyadenylated RNA by using poly(T) oligos targeting the poly(A)-tail of mRNA instead of depleting rRNA. In projects focusing on ncRNA, rRNA-depletion seems to be a more appropriate choice, which allows also quantification of pre-mRNA that has not been post-transcriptionally modified. Of note, slight differences exist between rRNA-depletion protocols in terms of rRNA removal efficiency and differential coverage of small genes, which should be investigated before selecting a method (Huang et al., 2011).

As mentioned above, selective protocols have been developed to target specific RNA molecules, such as small RNA, which are key regulators of gene expression. As small RNAs are lowly abundant, short in length (15–30 nt) and lacks polyadenylation, a separate strategy is often preferred to isolate and profile these RNAs using commercially available extraction kits similar to total RNA isolation protocols. Most kits involve isolation of small RNAs by size fractionation using gel electrophoresis or capture using silica spin columns. After isolation of small or other preferred RNA species from the total RNA, the sample is ready for library generation, which is universal for most RNA-Seq preparations. This step contains by fragmentation, reverse transcription into double stranded cDNA and adapter ligation. Even though there might be small differences in library preparation depending on the used platform/platform provider, but overall, these steps are generally applied. The fragmentation of reads, or in simple terms cutting the total RNA into smaller pieces, can be achieved by physical (e.g., sonication), enzymatic (e.g., RNAse II, transposase) or chemical (e.g., heat) means. The subsequent cDNA synthesis is essential for stability and improves confidence of base calling, which decreases with read length. Adapter ligation is necessary for sequencing, but also determines the next step, single-end or paired-end (PE) sequencing. In brief, single-end reads will only be sequenced from on end, either 3′ or 5′. PE sequencing starts to sequence at one end, and after a predetermined read length (∼100 bp), it stops and starts to read from the opposite end of the fragment, thus generating two reads per transcript (Cheng et al., 2020). Hence PE keeps strand information and is therefore more suited for studies of isoforms (Piskol et al., 2013).

Short fragmented sequencing is the most commonly used method, but involves a higher false-discovery rate in terms of reconstruction and read counting. To overcome this, long-read technologies have been developed to enable sequencing of entire transcripts from 5 ́end to 3 ́end, thus providing improved coverage. Companies such as PacBio and Oxford Nanopore technologies have provided direct sequencing of RNA platforms that belong to the Third Generation of sequencing and are capable to generate long-reads of around 10 kb. These long reads allow coverage of entire transcripts and improve the identification of new splicing events and eliminate amplification bias. However, a downside to these technologies is a lower sensitivity due to the high number of discarded reads during the pre-processing step, due to the high RNA degradation rate, which ameliorates the integrity of the transcripts and jeopardizes the accuracy of the reads and data analysis. Furthermore, RNA sequencing lacks the proofreading exonuclease activity as well as strand directionality information leading to elevated error rates, which are around 8–10%, compared to <1% of Illumina sequencers (Kovaka et al., 2019) direct (Levin et al., 2010; Venkataraman et al., 2018).

Conclusively, there are various sequencing methodologies available that focus on sequencing specific RNA molecules or targeted regions. This is further attributed by specialized computational pipelines focusing on specific RNA classes. Therefore, prior knowledge of available and current sequencing technologies as well as study designing with utmost care, will greatly benefit the impact and quality of produced results.

Single-Cell Sequencing

In 2013, Nature Methods titled scRNA-Seq as one of the most anticipated technologies of the year (Method of the Year 2013, 2014) with follow up nomination as technology of the year in 2019 due to its key role in the identification of cell types and functions, in addition to possible simultaneously multi-omics approaches in the future, highlighting its extensive role in genomic research.

Over the last couple of years, scRNA-Seq has had a massive effect on research. The reason is simple - while bulk RNA-Seq can measure the average gene expression across cells in a sample, identify differences between sample conditions and give a representation of highly regulated pathways, it fails to demonstrate the individual complexity of each cell and heterogeneity of tissues. Furthermore, some cell populations have high degrees of cellular and transcriptomic heterogeneity due to different cell types or indiscriminate states. The advent of scRNA-Seq technologies has addressed most of these limitations by facilitating the analysis of the transcriptome of every cell in each sample (Shapiro et al., 2013; Islam et al., 2014). This has enabled the creation of cell atlases (Regev et al., 2017; Cusanovich et al., 2018; Han et al., 2018; Tabula Muris Consortium., 2018; Cui et al., 2019; Litviňuková et al., 2020; Suryawanshi et al., 2020) at unprecedented resolution, the analysis of thousands of cells in parallel and even allows integration of chromatin status and multimodal analysis.

As a result, scRNA-Seq has shown to be a helpful technique for identifying cell subpopulations and elucidating dynamic cellular transitions during development and differentiation with an unparalleled level of detail and accuracy (Lafzi et al., 2018; Potter, 2018).

Several scRNA-Seq methods and technologies have been developed in recent years, all having different criteria regarding RNA transcript length, number of captured cells and read depth per cell (Svensson et al., 2018). Therefore, a specific scRNA-seq technology may be more beneficial for certain types of material than others, but all of them share similar workflows: samples preparation, single-cell capture, transcription and amplification, library preparation, sequencing and analysis (Potter, 2018). In general, there are three main technologies of scRNA-Seq: plate-, microfluidic- and droplet-based. In this section we will first explain the basic concepts of these currently used platforms and then their applications in cardiovascular science.

Single-Cell Dissociation and Capture

Sample preparation, similar to bulk RNA-Seq, is critical to enable single-cell capture and high-quality data. It starts with dissociation of material or tissue into a single cell form in order to extract cellular RNA. The main challenges included in this step are the fragility of the starting sample, physical stress, the choice of buffers, the duration of cell dissociation and the yield of individual and viable cells as most protocols require living cells (Lafzi et al., 2018; Nguyen et al., 2018). Cell isolation is a delicate process and obtaining a high-quality sample is of immense importance for a study to be successful. Therefore, single-cell dissociation is mostly achieved enzymatically using optimized protocols to limit cell lysis and loss of valuable RNAs. However, for all techniques, it is crucial to keep the processing time to a minimum to avoid cell damage and to prevent the unnecessary expression of stress-related genes, thereby altering original cellular transcriptomic profiles. It is noteworthy that cell isolation and single cell solution creation highly depends on the tissue in question. Primary samples, particularly when obtained from patients during an intervention, tend to be snap frozen, which is mostly incompatible with downstream analysis such as scRNA-Seq preparation from live cells. Creating single-cell solutions from snap frozen samples normally yield a lower overall number of cells and less viable cells (Reichard and Asosingh, 2019). Therefore, the most recommended for a single molecule approach is single nuclei sequencing (more information below) (Mimpen et al., 2021). Even if fresh material is available, certain tissue characteristics, e.g., as observed in disease depended fibrosis or scaring, can severely affect the digestion procedure. Therefore, we suggest thoughtful optimization to ensure high-quality scRNA-Seq data from primary tissue (Litvinukova et al., 2018; Vieira Braga and Miragaia, 2019).

After successful dissociation, single cells must be captured. As mentioned above, there are currently multiple techniques, and the choice largely depends on the sample of interest. However, capturing a high number of cells in the range of ≥10,000 cells will ensure significant improvements in data quality. One of the most critical properties determining the capture method is cell size. Plate-based methods incorporate either the use of micropipettes or fluorescence-activated cell sorting (FACS), guiding a single cell per well of a 96-well or 384-well plate. This method is mostly unaffected by size and has the option of long-term storage, as each cell is directly lysed within a well (Gladka et al., 2018; Hwang et al., 2018). However, since most steps have to be performed per well, the number of cells that can be processed at once is limited.

Microfluidic-based methods on the other hand separate the cells using narrow parallel microfluidic channels, where cell capture, cell lysis, reverse transcription, and multiplexing take place within an integrated fluidic circuit chip. In addition, a nice feature of this technique is the possibility to view capture cells before reverse transcription. An advantage of these methods is the high recovery/capture rate; however, they do require homogeneity in cell size. Furthermore, most microfluidic platforms require an input of >10.000 cells and are limited due to their restricted number of capture sites per microfluidic array.

Droplet-based approaches require the encapsulating of single cells into oil droplets with cell-specific barcoding. Due to its massive parallelization, ∼10.000 cells per sample can be captured within each run. For instance, 10X Genomics recommends a starting point of at least ∼1,600 cells per reaction for 3’ analyses, resulting in a recovery of ∼1.000 cells, and a multiplet rate of ∼0.8% (XGenomics, 2019). However, these platforms are limited to cell sizes less than 30 μm in diameter, and cells larger than that will clog the nozzle of the droplet system. Moreover, droplet- and microfluidic-based approaches have the need for living populations of single cells as input. Additionally, it is advised to remove cell aggregates, dead cells, and cell debris, before capturing to ensure a high percentage of viability among the selected population.

ICELL8 (Goldstein et al., 2017) is another promising platform, which is nanowell-based and therefore mostly similar to plate-based methods but with certain advantages. It utilizes a large-bore nozzle dispenser to distribute 1,000–1,800 single cells from diluted cell suspensions into 5.184 nanowells. Nanowell-specific barcoded are used to track sequencing data to its originator cell. Furthermore, ICELL8 has an imaging system visualizing all nanowells and by using fluorescence signals, it allows the differentiation of viable from dead cells and single-from multi-cells. A summary of main features of available scRNA-seq platforms can be found in Table 1.

TABLE 1
www.frontiersin.org

TABLE 1. Main features of available single-cell methodologies.

Sequencing of Single-Cells

Post cell capture, the individual cells/droplets are lysed, converted into cDNA via reverse transcription and finally sequenced. These steps depend highly on the capturing method chosen and show differences on multiple sequencing aspects. In terms of sequencing, plate-based methods rely on individual reverse transcription in each well, which can limit throughput and increase noise in downstream analyses. However, these methods allow full-length transcript sequencing, preferable for the identification of isoform splicing. Overall, plate-based platforms generally have high sensitivity and can reliably quantify up to 10,000 genes per cell.

Microfluidic-base methods are high in sensitivity and can use full-length transcript sequencing as well. Automatic systems like the Fluidigm C1 originally using Smart-Seq (Ramsköld et al., 2012; Picelli, 2019) before moving on the CEL-Seq method, were among the first scRNA-Seq platforms, increasing the transcript sensitivity and gene detection (Hashimshony et al., 2016). Smart-Seq methods can capture full-length transcripts (Ramsköld et al., 2012; Picelli et al., 2014) and therefore allowing the possibility to analysis alternate isoform splicing, which is further developed as a primary focus in Smart-Seq3 (Hagemann-Jensen et al., 2020). However, in comparison CEL-seq2 is limited to 3′-end reading and, therefore, cannot detect alternatively spliced isoforms, microRNAs or other non-polyadenylated transcripts.

Droplet-based methods have various advantages such as their massive parallelization, but they are size limited. Some commercial systems, such as the 10x Genomics Chromium, enable high-throughput processing by 3′- or 5′-end sequencing, but unfortunately also show reduced transcript recovery rates compared to other methods (Papalexi and Satija, 2018) with a read-depth of 104 to 105 per cell (Hwang et al., 2018). Even though this seems to be a drawback compared to the other techniques, droplet-based methods remain sufficient for the large-scale profiling of complex heterogeneous samples.

The actual sequencing of scRNA-Seq material can be performed on commonly used machines, and each single cell platform has particular features and come with individual recommendations to ensure best sequencing results and demultiplexing compatibilities (matching single cell tags with individual cells and their RNA content). However, one current limitation of high throughput scRNA-Seq, is that it generates either 3′ or 5′ sequence information via short length sequencing. This restricts analysis on splicing and sequence heterogeneity for most of the transcripts. As mentioned briefly above, short read sequencing involves fragmentation, so originally the transcripts are indeed full length and could theoretically undergo direct RNA sequencing using Third Generation Sequencing. i.e., Oxford Nanopore which is compatible with certain 10x Genomics protocols. Conversely, protocols for full or long length sequencing in single cell approaches are basically absent. This is largely due to technical and predominantly downstream analysis hardships, including the high error rate of these platforms. Nevertheless, there are methods, which attempt to correct errors and integrate single cell long read sequencing (Volden et al., 2018; Singh et al., 2019; Lebrigand et al., 2020; Zhang et al., 2021). For example, the rolling circle to concatemeric consensus (R2C2) method can produce full-length cDNA sequences, achieving ∼98% (Volden et al., 2018; Cole et al., 2020; Volden and Christopher, 2021) or Single-cell Nanopore sequencing with UMIs (ScNaUmi-seq) with ∼99.8% accuracy (Lebrigand et al., 2020), compared to former ∼50%. However, some of these methods require sufficient sequencing coverage to call consensus reads or high sequencing depth.

Single-Cell Application on Cardiac Tissue

At least 50% of the cells in an adult human heart are cardiomyocytes (CM), providing the contraction needed to keep our heart beating. The human adult CM are >100 μm in diameter along the major axis which can cause problems for several steps within the creation of single cell data. Options for overcoming the size limitation are the use of fetal mammalian CM, which are smaller in size. Another more recent approach includes the sequencing of human induced pluripotent stem cells (hiPSCs), which have been differentiated into CM. While the verdict on their comparability to adult CM due to their immaturity and heterogeneity (Shi et al., 2017; Paik et al., 2020) is still out, these cells have the advantage of smaller size and their non-invasive human origin.

In general, plate-based methods are not particularly limited by size, making them seem like a preferable platform for CMs. However, due to their expensive single-cell selection and limit in number of processed cells at a time, it has to be thoroughly evaluated. Furthermore, studies performed suggested the use of a FACS nozzle of 500 μm in diameter (Kannan et al., 2019) rather than the conventional 70–130 μm (Gladka et al., 2018) to ensure no terminal damage to the live CM, which may affect laser delays.

Even so, microfluidic- and plate-base systems are also commonly used for large and fragile cells, such as CMs. However, microfluidic-based systems would require a prior sorting via FACS into different cell types due to their vulnerability to heterogeneity and the difference in cell sizes within the heart. Furthermore, the elongated shape of CMs, in addition to unfortunate positioning could cause clogging of the system.

Unfortunately, droplet-based methods like the 10X Genomics Chromium, cannot cope with cells in the size range of mature adult CMs at present. However, as all these systems tend to evolve very quickly, continued protocol optimization and reductions in cost are to be expected soon.

Single Nuclei Sequencing

Per definition single-nuclei RNA sequencing (snRNA-Seq) is not a technique by itself, but more a modification of the scRNA-Seq methodology. However, due to its unique contributions and applications to the cardiac field, the technique has earned an exclusive overview of its own.

Single-cell sequencing is a powerful tool allowing in-depth characterization of cell populations within complex tissues. However, as discussed in the previous sub-section, scRNA-Seq has two major limitations when applied to cardiovascular tissues: dissociation of material and cell size. ScRNA-Seq systems require dissociation of tissue material and especially gaining high-quality single-cell resolution of the adult mammalian heart is rather difficult. Secondly, technical limitations regarding capturing techniques, leading to an under representation of individual cell types (i.e., CM) due to their large cell size and irregular shape (Ackers-Johnson et al., 2018). There are few platforms for sequencing individual adult CMs and given that all of them are plate-based systems, the number of cells they can analyze is limited.

As the name suggests, snRNA-Seq utilizes only the nuclear component of single cells as input material, which is most appropriate when viable intact cells cannot be harvested from fresh tissue. Furthermore, there are other advantages of snRNA-Seq, when scRNA-Seq capturing methodologies cannot be applied. For instance, it allows the extraction of nuclei from frozen samples, which is extremely valuable for rare materials like the heart. Moreover, scRNA-Seq involves a dissociation step, which might damage sensitive cells while failing to release cells more tightly bound or those surrounded by a matrix. Although sampling a much larger number of cells can partially overcome this limitation, it may not always be financially feasible. Additionally, as mentioned before current enzymatic and mechanical methods for single-cell dissociation tend to introduce stress-induced transcriptional artefacts. Even with commercial nuclei extraction kits, the solution should be examined to ensure no big cells are present anymore. Afterwards, single-nuclei can be prepared to be captured with any of the above-described methodologies.

Nevertheless, there are some strategic limitations to be considered. While some studies have shown similar comparable transcriptomic results acquired from sc- and snRNA-Seq (Bakken et al., 2018), by utilizing RNA from within the nucleus only, a fair amount of the total RNA is excluded from the analysis, which could constitute a limiting factor in identifying dynamic cellular states.

Furthermore, taking single nuclei from frozen samples might result in false transcriptomic representation due to RNA degradation by freezing. Interestingly, CM can have a diploid or even polyploid as well as multinucleate nature (Landim-Vieira et al., 2020), which is not fully understood yet and therefore very valuable for future studies. In addition, extracting the nuclei removes the highly abundant mitochondria in CMs leading to possible alterations in gene expression or missing information when looking into network-based approaches.

Integrational Approaches

While scRNA-Seq is currently the most widespread RNA-Seq technology, its prowess in deciphering the individual make-up of a single cell is just beginning to evolve. In recent years, many suppliers were able to integrate multi-omic approaches on the single cell level. Immunological studies have benefitted from paired B-cell or T-cell receptors, surface protein expression together with gene expression from a single cell (DeBerge et al., 2021). Single nuclei sequencing can be combined with ATAC (Assay for Transposase Accessible Chromatin) allowing for the analysis of chromatin accessibility at the single cell level, thus providing insights into cell types and states, and deeper understanding of gene regulatory mechanisms (Zhang et al., 2021). Further applications include proteomics (Cheung et al., 2021) and DNA methylation (Galvão and Kelsey, 2021). However, as most of these are still in their early stages including protocol optimization and integrational data analysis tools, it needs to be individually assessed whether specific research questions can be answered using these approaches (Jansen et al., 2019; Argelaguet et al., 2020). On the other side, computational approaches tried to leverage cell-type specific gene expression profiles from multiple scRNA-seq reference datasets to bulk RNA-Seq with the intention to apply cell compositions even to bulk RNA-Seq data. To date, multiple tools have been suggested, such as DESCEND (Wang et al., 2018), SCDC (Dong et al., 2021) and some others, but to our knowledge not one has been yet considered to be “state-of-art”.

Spatial Transcriptomics

Identifying gene expression profiles at single-cell resolution has revolutionized the field of transcriptomics. Using single-cell or single-nuclei based approaches has made it possible to look at individual cells within complex tissues, gaining valuable biological insights into rare cellular properties, cell-to-cell variability, and tissue identity. Despite such advantages, these approaches are limited by technical criteria to a certain degree. In this review, it has already been demonstrated how dissociating tissue can be a difficulty. In addition, this process causes the loss of information regarding the spatial organization of cells and cell populations throughout a tissue causing limited or incomplete interpretation. In in vivo situations though, spatial-specific expression has an enormous impact on biological networks. Gaining knowledge of spatial information in addition to timing, and level of developmental gene expression allows the description of interactive biological networks, where each element is influenced by its surrounding microenvironment (Zheng et al., 2017). Therefore, the persevation of spatial information holds promising approaches of combining genomic, transcriptomic, and proteomic features while retaining positional information.

Various scRNA-Seq methodologies have been discussed in the section above, but spatial approaches also have made extraordinary progress in the recent years (Ståhl et al., 2016; Qian et al., 2020) and gained increasing popularity. In general, spatial techniques vary from in situ hybridization (e.g., smFISH, seqFISH) over in situ sequencing (e.g., BaristaSeq, STARmap) to in situ capturing technologies (e.g., ST, Slide-Seq). All of them have their designated place in research, in situ spatial methodologies hold the prospect of possible combination with scRNA-Seq technologies. While some of these technologies and their drawbacks and advantages have been reviewed before (Asp et al., 2020), this review will highlight a few of the spatial methodologies and their applicability to cardiovascular research.

In situ Capturing Technologies

The idea behind in situ capturing technologies is to capture RNA in situ and then perform sequencing ex situ, avoiding typical limitations of direct visualization (e.g., limited marker amount, fluorescence exposure) as well as allowing an unbiased analysis of the complete transcriptome. The first of such technology to be developed was the Spatial Transcriptomics (ST), published 2016 (Ståhl et al., 2016), which was the basis for all following techniques. The basic idea behind this technology is transferring thin tissue section onto glass slides coated with positional molecular barcodes and synthesized cDNA and therefore able to capture mRNA while still maintaining positional representation. In detail, these slides are coded with the barcoded RT primers which the tissue is fixed, stained, imaged and permeabilized upon. The mRNA molecules diffuse down to the slide surface during the permeabilization phase and hybridize locally to the RT primers. Afterwards, the tissue is removed, the RNA is captured and retrotranscribed in situ and the cDNA-mRNA complexes are used for library preparation in the following ex situ NGS. By overlaying the different images, the barcoded reads are superimposed back onto the tissue image.

This approach provided spatial information on the transcriptome of the tissue, but the barcoded regions were only 10 μm in diameter. Depending on the tissue or cells of interest, this could be enough to reach a single-cell resolution. To tackle this, 10X Genomics acquired ST and started to further develop it into “10X Visium” by the end of 2018 (Tirado-Lee, 2020). 10X achieved further improvements in terms of resolution (55 μm in diameter and smaller distance between barcoded regions) and protocol running time. In addition, 10X has further expanded by adding Visium Gateway slides to their Visium repertoire (Tirado-Lee, 2020). Interestingly, 10X Genomics announced a new update “Visium HD”, which supposedly offers visibility 400 times higher than the current Visium solution with improvements addressing RNA degradation problems during formalin fixed paraffin embeddings.

Another method using a similar approach, Slide-Seq, was introduced in 2019 (Rodriques et al., 2019), promising a higher spatial-resolution technique by using smaller barcoded beads (Vickovic et al., 2019). Recently, Stickels et al. (2020b) published Slide-SeqV2, combining improvements in library generation, bead synthesis and array indexing resulting ∼10-fold greater RNA capture rate than its predecessor.

In March 2019, NanoString announced the commercial launch of their GeoMx Digital Spatial Profiling instrument (NanoString). In contrast, to previously mentioned platforms, GeoMx is capable of the highly multiplexed detection of mRNA targets in FFPE tissues (Merritt et al., 2020). Using multiple sections of a sample will help to acquire not only the spatial transcriptomic profiles but in addition also the spatial protein profiles. GeoMx can utilize regions in the range of 10–600 μm in diameter, which have to be selected manually. After selecting the regions of interest, fluorescence labeled antibodies used as morphology markers are then excited with UV light, triggering the release of either RNA target probes or antibodies for protein targeting, coupled with barcoded tags. Conveniently, using the NanoString’s nCounter instrument will assist to collect and quantify the tags from the regions of interest, resulting in a digital quantification of RNA expression with spatial context (Merritt et al., 2020). The workflow of selecting regions of interest makes it possible to analyze almost whole tissue sections, but the selection is still manual and therefore a biased process. However, despite increased advances leading to capture sensitivity and efficiency (Stickels et al., 2020a), scRNA-Seq data is still required to help mapping of cell types.

Data Analysis and Computational Approaches

All above mentioned procedures have a tremendous effect on data quality. However, in order to draw meaningful (biological) conclusions, downstream data procession and analysis is necessary. In fact, data analysis should be taken into consideration when designing a study in the first place, to allow for a qualitative interpretation of the data. Therefore, here we discuss the basic steps of data processing, analysis and opportunities for downstream analysis (Figure 1). Of note, this section will not be separated into bulk and single-cell sequencing as both follow similar pipelines. In addition, many single-cell advances are actually adapted from existing bulk approaches. At this stage, we would like to refer to other publications addressing data analysis steps and challenges (Koch et al., 2018; Dal Molin and Di Camillo, 2019; Lähnemann et al., 2020) in single cell vs. bulk RNA approaches.

FIGURE 1
www.frontiersin.org

FIGURE 1. Data Processing and Analysis pipeline. After sequencing files need to be pre-, processed and analysis (individually colored) in order to interpret the data and put it into biological context. These steps are very similar for single-cell (left) and bulk (right) RNA-Seq approaches. Italics represent most common used tools for the specific tasks.

After sequencing the libraries, the raw files need to be pre-processed. This step includes demultiplexing, quality control, trimming and alignment/de novo assembly. Even though existing in both bulk and single cell analysis, demultiplexing refers to slightly different processes. In bulk RNA-Seq, a sequencing lane typically will contain a pool of barcoded samples, which therefore require demultiplexing. In scRNA-Seq, each cell is considered a sample and, in most protocols, defined by a cell-specific barcode tag to each read. Therefore, each library normally contains multiple cell barcodes, which need to be distinguished as a sample. This process is called demultiplexing and gives each cell in each sample its specific RNA content. From here on, the steps of processing are in common between approaches, which are meant to cut off primer/adapters or poor-quality nucleotides, alignment of reads to a reference genome and counting reads per region/gene, creating a gene count matrix. Due to the overall high input amount for bulk sequencing, a high coverage can be achieved, but unfortunately scRNA-Seq normally deals with to the low concentration of input RNA and results in a high number of zero reads counts. However, these “false” zero counts, such as dropout events or sparsity, may mask the biological zero counts due to the non-expression of a group of genes in specific cells (Eraslan et al., 2019). Moreover, this high number of zero read counts in the gene expression matrix ultimately results in the reconstruction of an incomplete transcriptome, which renders the study of the different cellular phenotypes very challenging (Wang J. et al., 2019). Sparsity can impede downstream analysis and remains a challenge. To overcome this loss of information, computational imputation or denoising was developed, i.e., MAGIC and scIMPUTE (Lafzi et al., 2018; van Dijk et al., 2018). The last step is normalization. There are different approaches to normalization depending on the software used, but most of them are equally accepted and their individual usage mostly depends on the intented downstream analysis. However, the methods developed for bulk RNA-Seq such as RPKM, FPKM and TMM were not able to successfully remove the effects of different sampling from scRNA-Seq datasets, due to the considerable technical noise present in such data. Therefore, new normalization methods, which are able to handle noisy data and take advantage of the presence of UMIs were developed (Townes and Irizarry, 2020).

Afterwards, the data can be processed in various ways. Dimensional reduction is one of the most commonly used features and helps to visualize data and to identify relevant genes. For bulk-based studies, principal component analysis (PCA) is the most common used feature. Although this also applies to single cell approaches, because of the high dimensionality of the data other methodologies, such as t-SNE and UMAP, are superior to PCA. Other limitations include errors introduced by the inevitable technical variability during sample processing steps, sequencing depth or pipetting accuracy, also known as batch effect.

These are not unique to one methodology and can be found in the datasets generated by both, bulk and scRNA-Seq. However, these confounding factors are more amplified in scRNA-Seq, which may lead to a mix of technical variations with biological variables. In bulk RNA-Seq, instead, these batch effects are smaller and do not result as systematic errors, as bigger sample sizes usually lead to large-scale sample preparation and sequencing in parallel. Therefore, batch effect correction in scRNA-Seq requires other features capable of handling these variations, for example nearest neighbor integration. Further typical analyses include clustering and differential expression gene (DEG) detection. Unsupervised clustering further helps to identify groups and subgroups within single cell datasets and to pinpoint relevant cell groups for further analysis. DEG analysis has been established as a powerful tool to discover different regulation between groups, which can reveal regulatory effects. Although DEG can be also used to make between-sample comparisons, it is usually applied to identify transcriptome signatures that differentiate between cell groups. However, cell group assignment is still constrained by ambiguities, which could mask biological differences. Therefore, an improved reference set for cell type assignment is necessary to make DEG analysis in single-cells more reliable.

However, scRNA-Seq also provided new bioinformatics tools, with RNA velocity and trajectory analysis being the most prominent features. Both aim to describe the dynamic changes among captured cell signatures. RNA Velocity takes advantage of the ratio between spliced and unspliced RNA to develop dynamical patterns with possible cell fate prediction, while the trajectory attempts to predict a path that certain cells undergo during specific conditions (i.e., differentiation). Especially, deconvolution of RNA seq data has gained increased recognition in the last year (Baron et al., 2016; Wang X. et al., 2019; Newman et al., 2019). Deconvolution attempts to use single cell datasets to deconvolute bulk datasets by leveraging cell-type specific gene expression profiles from multiple scRNA-seq reference datasets, and thus synergistically combines specific information gained from both, bulk and single cell datasets.

Sequencing the Heart

Bulk Sequencing in Cardiovascular Research

While extracting RNA from heart tissue has no particular limiting factor, human heart samples are highly valuable due to the many ethical restrictions, making them hard to obtain. Many cardiac studies are therefore performed using other mammalian models such as mice or human tissue-derived cell lines. Nevertheless, patients with heart disease can undergo surgery, allowing for ascertainment of cardiac biopsies, but the amount of tissue retrieved can be insufficient or compromised. However, healthy individuals normally do not usually tend to undergo biopsy, rendering the collection of healthy living heart tissue as controls close to impossible.

Nevertheless, RNA-based studies have provided important insights to cardiac research, especially helping to decipher ncRNA and their potential role as regulators of numerous cellular processes in the progression of cardiovascular diseases Transcriptomic studies using NGS enabled clear distinguishment between various cardiac cell types, such as CMs, endothelial cells, fibroblasts, and immune cells. Arguably one of the first studies that compared accuracy and speed of NGS with microarray analyses demonstrated that in the Gαq transgenic mouse cardiomyopathy model RNA sequencing was accurate and sensitive enough to detect abundant and even rare cardiac transcripts (Matkovich et al., 2010). Subsequent studies highlighted that even atrial and ventricular cells show clear significant differences in their transcriptomic profiles (Sehnert et al., 2002; Ng et al., 2010). While atrial CM are mostly described with Myh4/ALC1 and Myh6/αMHC, the ventricle expresses more Myh3/ELC, Myh7/βMHC and Myl2 myosin genes (Ng et al., 2010). In depth reports of the current knowledge on the different molecules in cardiac cell subtypes have been reviewed before (Zhou et al., 2018; Colpaert and Calore, 2019; Siasos et al., 2020). NGS was essential to reveal that >80% of the genome is transcribed in various classes of RNA with protein-coding transcripts only accounting for max. 3% of the genome, and the remainder of transcripts lacking coding ability (Dunham et al., 2012). Many of these noncoding transcripts are functionally active RNA molecules that can be subdivided into small noncoding RNAs (<200 nt), such as microRNAs (miRs), and longer noncoding RNAs (>200 nt) that include ribosomal RNAs, natural antisense transcripts and long noncoding RNAs (lncRNAs) (Engreitz et al., 2016). NGS has helped to discover cardiac microRNAs (Matkovich et al., 2010; Leptidis et al., 2013) and lncRNAs (Matkovich et al., 2014; Yang et al., 2014) in the developing and diseased human heart and how they can potentially function as therapeutic target molecules. While these studies contributed significantly to our current understanding of cardiac (patho)physiology, they failed to address the importance of cardiac cellular heterogeneity and cell-cell interactions.

Single-Molecule Sequencing and Cardiovascular Research

ScRNA-Seq is able to look at cell compositions in an unbiased fashion. In 2016, two labs first reported the use of scRNA-Seq on embryonic mouse hearts by investigating lineage-specific gene programs underling early cardiac development (DeLaughter et al., 2016; Li et al., 2016). To date, multiple researchers have used scRNA-based techniques on various cardiac materials allowing to look at cell populations at a single-cell resolution, and started a new age of transcriptomics including the generation of transcriptomic and epigenetic cell atlases of adult mice (Cusanovich et al., 2018; Han et al., 2018; Tabula Muris Consortium, 2018) and humans (Regev et al., 2017; Cui et al., 2019; Litviňuková et al., 2020; Suryawanshi et al., 2020). Investigations deciphering the heterogeneity and subpopulations of cell types define region- or disease-specific gene expression profiles not only benefit cell atlases but also help to identify novel molecular mechanisms relevant for cardiac disease and new therapies. For instance, Li et al. (2019) investigated cardiac endothelial cells following ischemic injury and were able to show clonal proliferation of resident endothelial cells post-myocardial infarction.

Single-cell approaches, however, do not stop here. The heart is a multi-cellular organ and the interaction between cell of ligands and receptors on target gene expression is not yet fully understood. Browaeys et al. recently developed a tool for modelling intercellular communication on the single-cell level by linking ligands to target genes. This as well as other promising tools enable investigation of possible signaling mediators in addition to gene interaction by expression profiling (Browaeys et al., 2020; Cang and Nie, 2020; Efremova et al., 2020; Gladka, 2020; Jin et al., 2021). Furthermore, network and trajectory analyses on single-cell data led to the identification of fibroblasts as a critical constituent in promoting cardiomyocyte maturity (See et al., 2017) and regulatory interactions between transcription factors and target genes (Jackson et al., 2020). Another recent study by Litviňuková et al. (2020) provided comprehensive transcriptomic data on six distinct cardiac regions of the adult human heart using both sc- and snRNA-Seq technologies. Similar to other reports, they observed the cellular heterogeneity of cardiomyocytes, pericytes, and fibroblasts, and additionally, they analyzed cell-to-cell interactions suggesting a direct interaction through NOTCH2. Another study using microfluidics approaches in healthy human donors was able to identify and study the cellular and transcriptional diversity in a healthy heart (Tucker et al., 2020). The ICELL8 platforms claims to be able to deal with a wide range of sizes, while still making high throughput sequencing possible. Wang et al. used this platform on 21,422 single cells-including CMs and non-CMs from normal, failed and partially recovered adult human hearts to reveal inter- and intra-compartmental heterogeneity in response to stress as well as CM contractility and metabolism as the source of changes in heart function (Wang et al., 2020).

Taken together, an increasing number of studies have used methodologies on single cell level to shed light on the roles of various cell populations within in the mammalian heart helping future cardiac research and advances with their findings.

Spatial Transcriptomics in Cardiovascular Research

Interpreting single-cell properties can be difficult (Asp et al., 2020; Stickels et al., 2020a), but by adding spatial context to the transcriptome it will greatly benefit the current understanding of many biological networks. However, the main hurdle for these spatial methods is restricted RNA capture efficiency, which becomes increasingly more challenging with higher resolution. Even though spatial transcriptomics has proved to be a valuable source of transcriptomic information, spatial transcriptomics in the heart, cardiac development, and regional changes in gene expression during heart maturation are strong areas of great interest for the future.

Two studies utilized spatial transcription offered by 10X Genomics and subsequently combined spatial transcriptomic data with scRNA-Seq through computational deconvolution methods in order to reconstitute detailed spatial transcriptomic maps of human myocardial infarction (Kuppe et al., 2020) and the development of the human heart (Asp et al., 2019). In 2017 the group of Asp et al. (2017) already used the ST method to study fetal markers and their localization within adult human cardiac biopsies and later on created a spatiotemporal atlas of the developing human heart in the first trimester (Asp et al., 2019) by combining scRNA-seq data with spatial transcriptomics and in situ sequencing. In 2020, Asp et al. (2020) used their approach of combing scRNA-Seq along with ST and in situ sequencing to uncover novel cell types, such as clusters of fibrosis-associated fibroblast-like cells and a subpopulation of cardiac muscle cells. Furthermore, they also visualized their result by integrating the spatial information into 3D transcriptional maps. A major caveat of this combined approach is that it maximizes the efficiency of using data information by providing both cell-type and regional information from tissues, e.g., human embryo tissue, that are not easily obtained.

Conclusion

Over the last decades, studying the transcriptome became standard when investigating biological processes of physiological and pathological mechanisms. The increased interest in studying the RNA and its properties led to development of more advanced technology, which started from whole tissue, to cell population and finally reached single-cell level.

Sequencing techniques evolved quickly and increased our understanding of molecular signaling, intercellular communication networks and rare cell sub-populations. However, there is no one-size fits all approach; each technique has method-specific advantages and limitations. Hence researchers need to design their study with utmost care and choose a technique that aligns with their desired goals and sample availability. This furthermore includes the evaluation of multiple factors such as sample availability, sample preparation, platform advantages and disadvantages and sequencing attributes. This especially holds true for applications in cardiovascular research (Figure 2). Regardless, single-cell analysis has revolutionized our understanding of cardiovascular development and disease. The applications and insights that single-cell analysis enable in the investigation of cell subpopulation variants and differentiation in cardiovascular disease and heart development are invaluable to both basic and medical research. Furthermore, numerous cardiac single-cell atlases have been developed in multiple organisms and in multiple contexts including genetic variation, biological sex, and cardiac injury, all of which will provide useful resources for future work. Considering the deluge of recent advances in single cell analysis, it is becoming essential to develop and optimize methods for the integration of data of different types and from different sources. For instance, the development of antibody-based cell hashing, enabling multiplexing by indexing samples, holds future prospects, by lowering costs and allowing more complex experimental designs. In addition, multi-omic approaches trying to combine epigenetics, transcriptomic and proteomic aspects on single-cell level hold promising future innovations. Ultimately, scRNA-Seq technologies are expected to continue to develop rapidly with continuous improvement of experimental and analytical methods.

FIGURE 2
www.frontiersin.org

FIGURE 2. Schematic figure of single cell sequencing workflow on cardiac tissue. A summary of the main steps during scRNA-Seq experiments (Dissociation, Capture, Sequencing and Data Analysis). The main cell capturing techniques, FACS-based (left), Microfluidic (middle) and droplet-based (right), are schematically displayed including a comparison of their main advantages and disadvantages. Created with BioRender.com.

Taken together, single-cell analysis of the heart has revealed previously underappreciated cellular heterogeneity and the importance of intercellular communication. This diversity of cardiac cell types (and cell subtypes) acting together likely contributes to the homeostatic maintenance of cardiac tissue and is integral in the complex biological processes driving cell differentiation, cardiovascular development, disease, and regeneration. Finally, the use of single-cell level analytics will enable the definition of a healthy cardiac cell system and thereby better equip therapeutic pursuit toward the maintenance of this healthy cell system during physiological stress.

Author Contributions

J-CH wrote the manuscript, GL provided preliminary data. LdW and MS acquired funding for the study.

Funding

LdW was funded by Dutch CardioVascular Alliance (ARENA-PRIME). LdW was further supported by a VICI award 918-156-47 from the Dutch Research Council and Marie Sklodowska-Curie grant agreements no. 813716 (TRAIN-HEART) and Marie Sklodowska-Curie grant agreements no. 765274 (iPlacenta). MS is funded by the DFG (RTG2220, Project 281125614) and Marie Skłodowska-Curie Grant Agreement 81371.

Conflict of Interest

LdW is co-founder and stockholder of Mirabilis Therapeutics BV, a spin-off company of Maastricht University.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Accerbi, M., Schmidt, S. A., De Paoli, E., Park, S., Jeong, D.-H., and Green, P. J. (2010). Methods for Isolation of Total RNA to Recover miRNAs and Other Small RNAs from Diverse Species. Methods Mol. Biol. 592, 31–50. doi:10.1007/978-1-60327-005-2_3

PubMed Abstract | CrossRef Full Text | Google Scholar

Ackers-Johnson, M., Tan, W. L. W., and Foo, R. S.-Y. (2018). Following Hearts, One Cell at a Time: Recent Applications of Single-Cell RNA Sequencing to the Understanding of Heart Disease. Nat. Commun. 9, 4434. doi:10.1038/s41467-018-06894-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Argelaguet, R., Arnol, D., Bredikhin, D., Deloro, Y., Velten, B., and Marioni, J. C. (2020). MOFA+: a Statistical Framework for Comprehensive Integration of Multi-Modal Single-Cell Data. Genome Biol. 21, 111. doi:10.1186/s13059-020-02015-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Asp, M., Bergenstråhle, J., and Lundeberg, J. (2020). Spatially Resolved Transcriptomes—Next Generation Tools for Tissue Exploration. BioEssays 42, 1900221. doi:10.1002/bies.201900221

PubMed Abstract | CrossRef Full Text | Google Scholar

Asp, M., Giacomello, S., Larsson, L., Wu, C., Fürth, D., Qian, X., et al. (2019). A Spatiotemporal Organ-wide Gene Expression and Cell Atlas of the Developing Human Heart. Cell 179, 1647–1660. e19. doi:10.1016/j.cell.2019.11.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Asp, M., Salmén, F., Ståhl, P. L., Vickovic, S., Felldin, U., Löfling, M., et al. (2017). Spatial Detection of Fetal Marker Genes Expressed at Low Level in Adult Human Heart Tissue. Sci. Rep. 7, 12941. doi:10.1038/s41598-017-13462-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Bakken, T. E., Hodge, R. D., Miller, J. A., Yao, Z., Nguyen, T. N., Aevermann, B., et al. (2018). Single-nucleus and Single-Cell Transcriptomes Compared in Matched Cortical Cell Types. PLoS One 13, e0209648. doi:10.1371/journal.pone.0209648

PubMed Abstract | CrossRef Full Text | Google Scholar

Baron, M., Veres, A., Wolock, S. L., Faust, A. L., Gaujoux, R., Vetere, A., et al. (2016). A Single-Cell Transcriptomic Map of the Human and Mouse Pancreas Reveals Inter- and Intra-cell Population Structure. Cell Syst 3, 346–360. e4. doi:10.1016/j.cels.2016.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Browaeys, R., Saelens, W., and Saeys, Y. (2020). NicheNet: Modeling Intercellular Communication by Linking Ligands to Target Genes. Nat. Methods 17, 159–162. doi:10.1038/s41592-019-0667-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Cang, Z., and Nie, Q. (2020). Inferring Spatial and Signaling Relationships between Cells from Single Cell Transcriptomic Data. Nat. Commun. 11, 2084. doi:10.1038/s41467-020-15968-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, C. W., Beech, D. J., and Wheatcroft, S. B. (2020). Advantages of CEMiTool for Gene Co-expression Analysis of RNA-Seq Data. Comput. Biol. Med. 125, 103975. doi:10.1016/j.compbiomed.2020.103975

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheung, T. K., Lee, C.-Y., Bayer, F. P., McCoy, A., Kuster, B., and Rose, C. M. (2021). Defining the Carrier Proteome Limit for Single-Cell Proteomics. Nat. Methods 18, 76–83. doi:10.1038/s41592-020-01002-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Chu, Y., and Corey, D. R. (2012). RNA Sequencing: Platform Selection, Experimental Design, and Data Interpretation. Nucleic Acid Ther. 22, 271–274. doi:10.1089/nat.2012.0367

PubMed Abstract | CrossRef Full Text | Google Scholar

Cloonan, N., Forrest, A. R. R., Kolle, G., Gardiner, B. B. A., Faulkner, G. J., Brown, M. K., et al. (2008). Stem Cell Transcriptome Profiling via Massive-Scale mRNA Sequencing. Nat. Methods 5, 613–619. doi:10.1038/nmeth.1223

PubMed Abstract | CrossRef Full Text | Google Scholar

Cole, C., Byrne, A., Adams, M., Volden, R., and Vollmers, C. (2020). Complete Characterization of the Human Immune Cell Transcriptome Using Accurate Full-Length cDNA Sequencing. Genome Res. 30, 589–601. doi:10.1101/gr.257188.119

PubMed Abstract | CrossRef Full Text | Google Scholar

Colpaert, R. M. W., and Calore, M. (2019). MicroRNAs in Cardiac Diseases. Cells 8, 737. doi:10.3390/cells8070737

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, Y., Zheng, Y., Liu, X., Yan, L., Fan, X., Yong, J., et al. (2019). Single-Cell Transcriptome Analysis Maps the Developmental Track of the Human Heart. Cell Rep 26, 1934–1950. e5. doi:10.1016/j.celrep.2019.01.079

PubMed Abstract | CrossRef Full Text | Google Scholar

Cusanovich, D. A., Hill, A. J., Aghamirzaie, D., Daza, R. M., Pliner, H. A., Berletch, J. B., et al. (2018). A Single-Cell Atlas of In Vivo Mammalian Chromatin Accessibility. Cell 174, 1309–1324. e18. doi:10.1016/j.cell.2018.06.052

PubMed Abstract | CrossRef Full Text | Google Scholar

Dal Molin, A., and Di Camillo, B. (2019). How to Design a Single-Cell RNA-Sequencing experiment: Pitfalls, Challenges and Perspectives. Brief. Bioinform. 20, 1384–1394. doi:10.1093/bib/bby007

PubMed Abstract | CrossRef Full Text | Google Scholar

DeBerge, M., Lantz, C., Dehn, S., Sullivan, D. P., van der Laan, A. M., Niessen, H. W. M., et al. (2021). Hypoxia-inducible Factors Individually Facilitate Inflammatory Myeloid Metabolism and Inefficient Cardiac Repair. J. Exp. Med. 218. doi:10.1084/jem.20200667

CrossRef Full Text | Google Scholar

DeLaughter, D. M., Bick, A. G., Wakimoto, H., McKean, D., Gorham, J. M., Kathiriya, I. S., et al. (2016). Single-Cell Resolution of Temporal Gene Expression during Heart Development. Dev. Cel 39, 480–490. doi:10.1016/j.devcel.2016.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, M., Thennavan, A., Urrutia, E., Li, Y., Perou, C. M., Zou, F., et al. (2021). SCDC: Bulk Gene Expression Deconvolution by Multiple Single-Cell RNA Sequencing References. Brief. Bioinform. 22, 416–427. doi:10.1093/bib/bbz166

PubMed Abstract | CrossRef Full Text | Google Scholar

Dunham, I., Kundaje, A., Aldred, S. F., Collins, P. J., Davis, C. A., Doyle, F., et al. (2012). An Integrated Encyclopedia of DNA Elements in the Human Genome. Nature 489, 57–74. doi:10.1038/nature11247

PubMed Abstract | CrossRef Full Text | Google Scholar

Efremova, M., Vento-Tormo, M., Teichmann, S. A., and Vento-Tormo, R. (2020). CellPhoneDB: Inferring Cell–Cell Communication from Combined Expression of Multi-Subunit Ligand–Receptor Complexes. Nat. Protoc. 15, 1484–1506. doi:10.1038/s41596-020-0292-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Engreitz, J. M., Ollikainen, N., and Guttman, M. (2016). Long Non-coding RNAs: Spatial Amplifiers that Control Nuclear Structure and Gene Expression. Nat. Rev. Mol. Cel Biol. 17, 756–770. doi:10.1038/nrm.2016.126

CrossRef Full Text | Google Scholar

Eraslan, G., Simon, L. M., Mircea, M., Mueller, N. S., and Theis, F. J. (2019). Single-cell RNA-Seq Denoising Using a Deep Count Autoencoder. Nat. Commun. 10, 390. doi:10.1038/s41467-018-07931-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Galvão, A., and Kelsey, G. (2021). Profiling Genome-Wide in Single Cells. Method. Mole. Biol. 2214, 221–240. doi:10.1007/978-1-0716-0958-3_15

CrossRef Full Text | Google Scholar

Garraway, L. A. (2013). Genomics-driven Oncology: Framework for an Emerging Paradigm. J. Clin. Oncol. 31, 1806–1814. doi:10.1200/JCO.2012.46.8934

CrossRef Full Text | Google Scholar

Gladka, M. M. (2020). Cellular Communication in a ‘virtual Lab’: Going beyond the Classical Ligand-Receptor Interaction. Cardiovasc. Res. 116, e67–e69. doi:10.1093/cvr/cvaa076

PubMed Abstract | CrossRef Full Text | Google Scholar

Gladka, M. M., Molenaar, B., de Ruiter, H., van der Elst, S., Tsui, H., Versteeg, D., et al. (2018). Single-Cell Sequencing of the Healthy and Diseased Heart Reveals Cytoskeleton-Associated Protein 4 as a New Modulator of Fibroblasts Activation. Circulation 138, 166–180. doi:10.1161/CIRCULATIONAHA.117.030742

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldstein, L. D., Chen, Y.-J. J., Dunne, J., Mir, A., Hubschle, H., Guillory, J., et al. (2017). Massively Parallel Nanowell-Based Single-Cell Gene Expression Profiling. BMC Genomics 18, 519. doi:10.1186/s12864-017-3893-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Hagemann-Jensen, M., Ziegenhain, C., Chen, P., Ramsköld, D., Hendriks, G.-J., Larsson, A. J. M., et al. (2020). Single-cell RNA Counting at Allele and Isoform Resolution Using Smart-Seq3. Nat. Biotechnol. 38, 708–714. doi:10.1038/s41587-020-0497-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, X., Wang, R., Zhou, Y., Fei, L., Sun, H., Lai, S., et al. (2018). Mapping the Mouse Cell Atlas by Microwell-Seq. Cell 172, 1091–1107. e17. doi:10.1016/j.cell.2018.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Hashimshony, T., Senderovich, N., Avital, G., Klochendler, A., de Leeuw, Y., Anavy, L., et al. (2016). CEL-Seq2: Sensitive Highly-Multiplexed Single-Cell RNA-Seq. Genome Biol. 17, 77. doi:10.1186/s13059-016-0938-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, R., Jaritz, M., Guenzl, P., Vlatkovic, I., Sommer, A., Tamir, I. M., et al. (2011). An RNA-Seq Strategy to Detect the Complete Coding and Non-coding Transcriptome Including Full-Length Imprinted Macro ncRNAs. PLoS One 6, e27288. doi:10.1371/journal.pone.0027288

PubMed Abstract | CrossRef Full Text | Google Scholar

Hwang, B., Lee, J. H., and Bang, D. (2018). Single-cell RNA Sequencing Technologies and Bioinformatics Pipelines. Exp. Mol. Med. 50, 96. doi:10.1038/s12276-018-0071-8

CrossRef Full Text | Google Scholar

International Human Genome Sequencing Consortium (2004). Finishing the Euchromatic Sequence of the Human Genome. Nature 431, 931–945. doi:10.1038/nature03001

PubMed Abstract | CrossRef Full Text | Google Scholar

Islam, S., Zeisel, A., Joost, S., La Manno, G., Zajac, P., Kasper, M., et al. (2014). Quantitative Single-Cell RNA-Seq with Unique Molecular Identifiers. Nat. Methods 11, 163–166. doi:10.1038/nmeth.2772

PubMed Abstract | CrossRef Full Text | Google Scholar

Jackson, C. A., Castro, D. M., Saldi, G.-A., Bonneau, R., and Gresham, D. (2020). Gene Regulatory Network Reconstruction Using Single-Cell RNA Sequencing of Barcoded Genotypes in Diverse Environments. Elife 9. doi:10.7554/eLife.51254

PubMed Abstract | CrossRef Full Text | Google Scholar

Jansen, C., Ramirez, R. N., El-Ali, N. C., Gomez-Cabrero, D., Tegner, J., Merkenschlager, M., et al. (2019). Building Gene Regulatory Networks from scATAC-Seq and scRNA-Seq Using Linked Self Organizing Maps. PLOS Comput. Biol. 15, e1006555. doi:10.1371/journal.pcbi.1006555

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, S., Guerrero-Juarez, C. F., Zhang, L., Chang, I., Ramos, R., Kuan, C.-H., et al. (2021). Inference and Analysis of Cell-Cell Communication Using CellChat. Nat. Commun. 12, 1088. doi:10.1038/s41467-021-21246-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Kannan, S., Miyamoto, M., Lin, B. L., Zhu, R., Murphy, S., Kass, D. A., et al. (2019). Large Particle Fluorescence-Activated Cell Sorting Enables High-Quality Single-Cell RNA Sequencing and Functional Analysis of Adult Cardiomyocytes. Circ. Res. 125, 567–569. doi:10.1161/CIRCRESAHA.119.315493

PubMed Abstract | CrossRef Full Text | Google Scholar

Koch, C. M., Chiu, S. F., Akbarpour, M., Bharat, A., Ridge, K. M., Bartom, E. T., et al. (2018). A Beginner’s Guide to Analysis of RNA Sequencing Data. Am. J. Respir. Cel Mol. Biol. 59, 145–157. doi:10.1165/rcmb.2017-0430TR

PubMed Abstract | CrossRef Full Text | Google Scholar

Kovaka, S., Zimin, A. V., Pertea, G. M., Razaghi, R., Salzberg, S. L., and Pertea, M. (2019). Transcriptome Assembly from Long-Read RNA-Seq Alignments with StringTie2. Genome Biol. 20, 278. doi:10.1186/s13059-019-1910-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kukurba, K. R., and Montgomery, S. B. (20152015). RNA Sequencing and Analysis. Cold Spring Harb. Protoc., 084970. doi:10.1101/pdb.top084970

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuppe, C., Flores, R. O. R., Li, Z., Hannani, M., Tanevski, J., Halder, M., et al. (2020). Spatial Multi-Omic Map of Human Myocardial Infarction. bioRxiv. doi:10.1101/2020.12.08.411686

CrossRef Full Text | Google Scholar

Lafzi, A., Moutinho, C., Picelli, S., and Heyn, H. (2018). Tutorial: Guidelines for the Experimental Design of Single-Cell RNA Sequencing Studies. Nat. Protoc. 13, 2742–2757. doi:10.1038/s41596-018-0073-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Lähnemann, D., Köster, J., Szczurek, E., McCarthy, D. J., Hicks, S. C., Robinson, M. D., et al. (2020). Eleven Grand Challenges in Single-Cell Data Science. Genome Biol. 21, 31. doi:10.1186/s13059-020-1926-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Landim-Vieira, M., Schipper, J. M., Pinto, J. R., and Chase, P. B. (2020). Cardiomyocyte Nuclearity and Ploidy: when Is Double Trouble. J. Muscle Res. Cel Motil. 41, 329–340. doi:10.1007/s10974-019-09545-7

CrossRef Full Text | Google Scholar

Lebrigand, K., Magnone, V., Barbry, P., and Waldmann, R. (2020). High Throughput Error Corrected Nanopore Single Cell Transcriptome Sequencing. Nat. Commun. 11, 4025. doi:10.1038/s41467-020-17800-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Leptidis, S., El Azzouzi, H., Lok, S. I., de Weger, R., Olieslagers, S., Kisters, N., et al. (2013). A Deep Sequencing Approach to Uncover the miRNOME in the Human Heart. PLoS One 8, e57800. doi:10.1371/journal.pone.0057800

PubMed Abstract | CrossRef Full Text | Google Scholar

Levin, J. Z., Yassour, M., Adiconis, X., Nusbaum, C., Thompson, D. A., Friedman, N., et al. (2010). Comprehensive Comparative Analysis of Strand-specific RNA Sequencing Methods. Nat. Methods 7, 709–715. doi:10.1038/nmeth.1491

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, G., Xu, A., Sim, S., Tsao, P. S., Quake, S. R., and Wu, S. M. (2016). Transcriptomic Profiling Maps Anatomically Patterned Subpopulations Among Single Embryonic Cardiac Cells Resource Transcriptomic Profiling Maps Anatomically Patterned Subpopulations Among Single Embryonic Cardiac Cells. Dev. Cel 39, 491–507. doi:10.1016/j.devcel.2016.10.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Solomonidis, E. G., Meloni, M., Taylor, R. S., Duffin, R., Dobie, R., et al. (2019). Single-cell Transcriptome Analyses Reveal Novel Targets Modulating Cardiac Neovascularization by Resident Endothelial Cells Following Myocardial Infarction. Eur. Heart J. 40, 2507–2520. doi:10.1093/eurheartj/ehz305

PubMed Abstract | CrossRef Full Text | Google Scholar

Lister, R., O’Malley, R. C., Tonti-Filippini, J., Gregory, B. D., Berry, C. C., Millar, A. H., et al. (2008). Highly Integrated Single-Base Resolution Maps of the Epigenome in Arabidopsis. Cell 133, 523–536. doi:10.1016/j.cell.2008.03.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Litvinukova, M., Lindberg, E., Maatz, H., Zhang, H., Radke, M., Gotthardt, M., et al. (2018). Single Cell and Single Nuclei Analysis Human Heart Tissue. protocols 10, eae3ae. doi:10.17504/protocols.io.veae3ae

CrossRef Full Text | Google Scholar

Litviňuková, M., Talavera-López, C., Maatz, H., Reichart, D., Worth, C. L., Lindberg, E. L., et al. (2020). Cells of the Adult Human Heart. Nature 588, 466–472. doi:10.1038/s41586-020-2797-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Matkovich, S. J., Edwards, J. R., Grossenheider, T. C., de Guzman Strong, C., and Dorn, G. W. (2014). Epigenetic Coordination of Embryonic Heart Transcription by Dynamically Regulated Long Noncoding RNAs. Proc. Natl. Acad. Sci. 111, 12264–12269. doi:10.1073/pnas.1410622111

PubMed Abstract | CrossRef Full Text | Google Scholar

Matkovich, S. J., Zhang, Y., Van Booven, D. J., and Dorn, G. W. (2010). Deep mRNA Sequencing for In Vivo Functional Analysis of Cardiac Transcriptional Regulators: Application to Galphaq. Circ. Res. 106, 1459–1467. doi:10.1161/CIRCRESAHA.110.217513

PubMed Abstract | CrossRef Full Text | Google Scholar

Merritt, C. R., Ong, G. T., Church, S. E., Barker, K., Danaher, P., Geiss, G., et al. (2020). Multiplex Digital Spatial Profiling of Proteins and RNA in Fixed Tissue. Nat. Biotechnol. 38, 586–599. doi:10.1038/s41587-020-0472-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Method of the Year 2013 (2014). Nat. Methods 11, 1. doi:10.1038/nmeth.2801

PubMed Abstract | CrossRef Full Text

Mimpen, J. Y., Paul, C., Network, T. S., Cribbs, A., and Snelling, S. (2021). Nuclei Isolation from Snap-Frozen Tendon Tissue for Single Nucleus RNA Sequencing. protocol.io. doi:10.17504/protocols.io.bc6xizfn

CrossRef Full Text | Google Scholar

Mortazavi, A., Williams, B. A., McCue, K., Schaeffer, L., and Wold, B. (2008). Mapping and Quantifying Mammalian Transcriptomes by RNA-Seq. Nat. Methods 5, 621–628. doi:10.1038/nmeth.1226

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagalakshmi, U., Wang, Z., Waern, K., Shou, C., Raha, D., Gerstein, M., et al. (2008). The Transcriptional Landscape of the Yeast Genome Defined by RNA Sequencing. Science (80-. ) 320, 1344–1349. doi:10.1126/science.1158441

PubMed Abstract | CrossRef Full Text | Google Scholar

NanoString Geomx Protein Assays, . Available at: https://www.nanostring.com/products/geomx-digital-spatial-profiler/geomx-protein-assays/(Accessed February 11, 2021).

Newman, A. M., Steen, C. B., Liu, C. L., Gentles, A. J., Chaudhuri, A. A., Scherer, F., et al. (2019). Determining Cell Type Abundance and Expression from Bulk Tissues with Digital Cytometry. Nat. Biotechnol. 37, 773–782. doi:10.1038/s41587-019-0114-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Ng, S. Y., Wong, C. K., and Tsang, S. Y. (2010). Differential Gene Expressions in Atrial and Ventricular Myocytes: Insights into the Road of Applying Embryonic Stem Cell-Derived Cardiomyocytes for Future Therapies. Am. J. Physiol. Physiol. 299, C1234–C1249. doi:10.1152/ajpcell.00402.2009

CrossRef Full Text | Google Scholar

Nguyen, Q. H., Pervolarakis, N., Nee, K., and Kessenbrock, K. (2018). Experimental Considerations for Single-Cell RNA Sequencing Approaches. Front. Cell Dev. Biol. 6, 108. doi:10.3389/fcell.2018.00108

PubMed Abstract | CrossRef Full Text | Google Scholar

Paik, D. T., Chandy, M., and Wu, J. C. (2020). Patient and Disease-specific Induced Pluripotent Stem Cells for Discovery of Personalized Cardiovascular Drugs and Therapeutics. Pharmacol. Rev. 72, 320–342. doi:10.1124/pr.116.013003

PubMed Abstract | CrossRef Full Text | Google Scholar

Papalexi, E., and Satija, R. (2018). Single-cell RNA Sequencing to Explore Immune Cell Heterogeneity. Nat. Rev. Immunol. 18, 35–45. doi:10.1038/nri.2017.76

PubMed Abstract | CrossRef Full Text | Google Scholar

Picelli, S., Faridani, O. R., Björklund, A. K., Winberg, G., Sagasser, S., and Sandberg, R. (2014). Full-length RNA-Seq from Single Cells Using Smart-Seq2. Nat. Protoc. 9, 171–181. doi:10.1038/nprot.2014.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Picelli, S. (2019). Full-Length Single-Cell RNA Sequencing with Smart-Seq2. Methods Mol. Biol., 25–44. doi:10.1007/978-1-4939-9240-9_3

PubMed Abstract | CrossRef Full Text | Google Scholar

Piskol, R., Ramaswami, G., and Li, J. B. (2013). Reliable Identification of Genomic Variants from RNA-Seq Data. Am. J. Hum. Genet. 93, 641–651. doi:10.1016/j.ajhg.2013.08.008

CrossRef Full Text | Google Scholar

Potter, S. S. (2018). Single-cell RNA Sequencing for the Study of Development, Physiology and Disease. Nat. Rev. Nephrol. 14, 479–492. doi:10.1038/s41581-018-0021-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Qian, X., Harris, K. D., Hauling, T., Nicoloutsopoulos, D., Muñoz-Manchado, A. B., Skene, N., et al. (2020). Probabilistic Cell Typing Enables fine Mapping of Closely Related Cell Types In Situ. Nat. Methods 17, 101–106. doi:10.1038/s41592-019-0631-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramsköld, D., Luo, S., Wang, Y.-C., Li, R., Deng, Q., Faridani, O. R., et al. (2012). Full-length mRNA-Seq from Single-Cell Levels of RNA and Individual Circulating Tumor Cells. Nat. Biotechnol. 30, 777–782. doi:10.1038/nbt.2282

PubMed Abstract | CrossRef Full Text | Google Scholar

Regev, A., Teichmann, S. A., Lander, E. S., Amit, I., Benoist, C., Birney, E., et al. (2017). The Human Cell Atlas. Elife 6, 27041. doi:10.7554/eLife.27041

PubMed Abstract | CrossRef Full Text | Google Scholar

Reichard, A., and Asosingh, K. (2019). Best Practices for Preparing a Single Cell Suspension from Solid Tissues for Flow Cytometry. Cytom. Part. A. 95, 219–226. doi:10.1002/cyto.a.23690

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodriques, S. G., Stickels, R. R., Goeva, A., Martin, C. A., Murray, E., Vanderburg, C. R., et al. (2019). Slide-seq: A Scalable Technology for Measuring Genome-wide Expression at High Spatial Resolution. Science 363, 1463–1467. doi:10.1126/science.aaw1219

PubMed Abstract | CrossRef Full Text | Google Scholar

Schroeder, A., Mueller, O., Stocker, S., Salowsky, R., Leiber, M., Gassmann, M., et al. (2006). The RIN: an RNA Integrity Number for Assigning Integrity Values to RNA Measurements. BMC Mol. Biol. 7, 3. doi:10.1186/1471-2199-7-3

PubMed Abstract | CrossRef Full Text | Google Scholar

See, K., Tan, W. L. W., Lim, E. H., Tiang, Z., Lee, L. T., Li, P. Y. Q., et al. (2017). Single Cardiomyocyte Nuclear Transcriptomes Reveal a lincRNA-Regulated De-differentiation and Cell Cycle Stress-Response In Vivo. Nat. Commun. 8, 225. doi:10.1038/s41467-017-00319-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Seelenfreund, E., Robinson, W. A., Amato, C. M., Tan, A.-C., Kim, J., and Robinson, S. E. (2014). Long Term Storage of Dry versus Frozen RNA for Next Generation Molecular Studies. PLoS One 9, e111827. doi:10.1371/journal.pone.0111827

PubMed Abstract | CrossRef Full Text | Google Scholar

Sehnert, A. J., Huq, A., Weinstein, B. M., Walker, C., Fishman, M., and Stainier, D. Y. R. (2002). Cardiac Troponin T Is Essential in Sarcomere Assembly and Cardiac Contractility. Nat. Genet. 31, 106–110. doi:10.1038/ng875

PubMed Abstract | CrossRef Full Text | Google Scholar

Shapiro, E., Biezuner, T., and Linnarsson, S. (2013). Single-cell Sequencing-Based Technologies Will Revolutionize Whole-Organism Science. Nat. Rev. Genet. 14, 618–630. doi:10.1038/nrg3542

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, Y., Inoue, H., Wu, J. C., and Yamanaka, S. (2017). Induced Pluripotent Stem Cell Technology: a Decade of Progress. Nat. Rev. Drug Discov. 16, 115–130. doi:10.1038/nrd.2016.245

PubMed Abstract | CrossRef Full Text | Google Scholar

Siasos, G., Bletsa, E., Stampouloglou, P. K., Oikonomou, E., Tsigkou, V., Paschou, S. A., et al. (2020). MicroRNAs in Cardiovascular Disease. Hell. J. Cardiol. 61, 165–173. doi:10.1016/j.hjc.2020.03.003

CrossRef Full Text | Google Scholar

Singh, M., Al-Eryani, G., Carswell, S., Ferguson, J. M., Blackburn, J., Barton, K., et al. (2019). High-throughput Targeted Long-Read Single Cell Sequencing Reveals the Clonal and Transcriptional Landscape of Lymphocytes. Nat. Commun. 10, 3120. doi:10.1038/s41467-019-11049-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Ståhl, P. L., Salmén, F., Vickovic, S., Lundmark, A., Navarro, J. F., Magnusson, J., et al. (2016). Visualization and Analysis of Gene Expression in Tissue Sections by Spatial Transcriptomics. Science 353, 78–82. doi:10.1126/science.aaf2403

PubMed Abstract | CrossRef Full Text | Google Scholar

Stickels, R. R., Murray, E., Kumar, P., Li, J., Marshall, J. L., Bella, D. Di., et al. (2020a). Sensitive Spatial Genome Wide Expression Profiling at Cellular Resolution. bioRxiv. doi:10.1101/2020.03.12.989806

CrossRef Full Text | Google Scholar

Stickels, R. R., Murray, E., Kumar, P., Li, J., Marshall, J. L., Di Bella, D. J., et al. (2020b). Highly Sensitive Spatial Transcriptomics at Near-Cellular Resolution with Slide-seqV2. Nat. Biotechnol. doi:10.1038/s41587-020-0739-1

CrossRef Full Text | Google Scholar

Suryawanshi, H., Clancy, R., Morozov, P., Halushka, M. K., Buyon, J. P., and Tuschl, T. (2020). Cell Atlas of the Foetal Human Heart and Implications for Autoimmune-Mediated Congenital Heart Block. Cardiovasc. Res. 116, 1446–1457. doi:10.1093/cvr/cvz257

PubMed Abstract | CrossRef Full Text | Google Scholar

Svensson, V., Vento-Tormo, R., and Teichmann, S. A. (2018). Exponential Scaling of Single-Cell RNA-Seq in the Past Decade. Nat. Protoc. 13, 599–604. doi:10.1038/nprot.2017.149

PubMed Abstract | CrossRef Full Text | Google Scholar

Tabula Muris Consortium (2018). Overall coordination, Logistical coordination, Organ collection and processing, Library preparation and sequencing, Computational data analysis, et alSingle-cell transcriptomics of 20 mouse organs creates a Tabula Muris. Nature 562, 367–372. doi:10.1038/s41586-018-0590-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Tirado-Lee, L. (2020). Introducing Visium Gateway: Your Entry to the World of Spatial Gene Expression. Available at: https://www.10xgenomics.com/blog/introducing-visium-gateway-your-entry-to-the-world-of-spatial-gene-expression (Accessed February 11, 2021).

Google Scholar

Townes, F. W., and Irizarry, R. A. (2020). Quantile Normalization of Single-Cell RNA-Seq Read Counts without Unique Molecular Identifiers. Genome Biol. 21, 160. doi:10.1186/s13059-020-02078-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Tucker, N. R., Chaffin, M., Fleming, S. J., Hall, A. W., Parsons, V. A., Bedi, K. C., et al. (2020). Transcriptional and Cellular Diversity of the Human Heart. Circulation 142, 466–482. doi:10.1161/CIRCULATIONAHA.119.045401

PubMed Abstract | CrossRef Full Text | Google Scholar

van Dijk, D., Sharma, R., Nainys, J., Yim, K., Kathail, P., Carr, A. J., et al. (2018). Recovering Gene Interactions from Single-Cell Data Using Data Diffusion. Cell 174, 716–729. e27. doi:10.1016/j.cell.2018.05.061

PubMed Abstract | CrossRef Full Text | Google Scholar

Venkataraman, S., Prasad, B., and Selvarajan, R. (2018). RNA Dependent RNA Polymerases: Insights from Structure, Function and Evolution. Viruses 10, 76. doi:10.3390/v10020076

PubMed Abstract | CrossRef Full Text | Google Scholar

Vickovic, S., Eraslan, G., Salmén, F., Klughammer, J., Stenbeck, L., Schapiro, D., et al. (2019). High-definition Spatial Transcriptomics for In Situ Tissue Profiling. Nat. Methods 16, 987–990. doi:10.1038/s41592-019-0548-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Vieira Braga, F. A., and Miragaia, R. J. (2019). Tissue Handling and Dissociation for Single-Cell RNA-Seq. Methods Mol. Biol. 1979, 9–21. doi:10.1007/978-1-4939-9240-9_2

PubMed Abstract | CrossRef Full Text | Google Scholar

Volden, R., and Christopher, V. (2021). Highly Multiplexed Single-Cell Full-Length cDNA Sequencing of Human Immune Cells with 10X Genomics and R2C2. bioRxiv. doi:10.1101/2020.01.10.902361

CrossRef Full Text | Google Scholar

Volden, R., Palmer, T., Byrne, A., Cole, C., Schmitz, R. J., Green, R. E., et al. (2018). Improving Nanopore Read Accuracy with the R2C2 Method Enables the Sequencing of Highly Multiplexed Full-Length Single-Cell cDNA. Proc. Natl. Acad. Sci. 115, 9726–9731. doi:10.1073/pnas.1806447115

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Agarwal, D., Huang, M., Hu, G., Zhou, Z., Ye, C., et al. (2019a). Data Denoising with Transfer Learning in Single-Cell Transcriptomics. Nat. Methods 16, 875–878. doi:10.1038/s41592-019-0537-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Huang, M., Torre, E., Dueck, H., Shaffer, S., Murray, J., et al. (2018). Gene Expression Distribution Deconvolution in Single-Cell RNA Sequencing. Proc. Natl. Acad. Sci. 115, E6437–E6446. doi:10.1073/pnas.1721085115

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Yu, P., Zhou, B., Song, J., Li, Z., Zhang, M., et al. (2020). Single-cell Reconstruction of the Adult Human Heart during Heart Failure and Recovery Reveals the Cellular Landscape Underlying Cardiac Function. Nat. Cel Biol. 22, 108–119. doi:10.1038/s41556-019-0446-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Park, J., Susztak, K., Zhang, N. R., and Li, M. (2019b). Bulk Tissue Cell Type Deconvolution with Multi-Subject Single-Cell Expression Reference. Nat. Commun. 10, 380. doi:10.1038/s41467-018-08023-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Gerstein, M., and Snyder, M. (2009). RNA-seq: a Revolutionary Tool for Transcriptomics. Nat. Rev. Genet. 10, 57–63. doi:10.1038/nrg2484

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilhelm, B. T., Marguerat, S., Watt, S., Schubert, F., Wood, V., Goodhead, I., et al. (2008). Dynamic Repertoire of a Eukaryotic Transcriptome Surveyed at Single-Nucleotide Resolution. Nature 453, 1239–1243. doi:10.1038/nature07002

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, K.-C., Yamada, K. A., Patel, A. Y., Topkara, V. K., George, I., Cheema, F. H., et al. (2014). Deep RNA Sequencing Reveals Dynamic Regulation of Myocardial Noncoding RNAs in Failing Human Heart and Remodeling with Mechanical Circulatory Support. Circulation 129, 1009–1021. doi:10.1161/CIRCULATIONAHA.113.003863

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, K., Hocker, J. D., Miller, M., Hou, X., Chiou, J., Poirion, O. B., et al. (2021). A Single-Cell Atlas of Chromatin Accessibility in the Human Genome. Cell 184, 5985–6001. doi:10.1016/j.cell.2021.10.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, G., Li, S., and Szekely, G. (2017). Statistical Shape and Deformation Analysis. 1st Edition. London: Academic Press. doi:10.1016/C2015-0-06799-5

CrossRef Full Text | Google Scholar

Zhou, S., Jin, J., Wang, J., Zhang, Z., Freedman, J. H., Zheng, Y., et al. (2018). miRNAS in Cardiovascular Diseases: Potential Biomarkers, Therapeutic Targets and Challenges. Acta Pharmacol. Sin. 39, 1073–1084. doi:10.1038/aps.2018.30

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: sequencing, single-cell, bulk-seq, transcriptomics, heart

Citation: Hegenbarth J-C, Lezzoche G, De Windt LJ and Stoll M (2022) Perspectives on Bulk-Tissue RNA Sequencing and Single-Cell RNA Sequencing for Cardiac Transcriptomics. Front. Mol. Med. 2:839338. doi: 10.3389/fmmed.2022.839338

Received: 19 December 2021; Accepted: 31 January 2022;
Published: 22 February 2022.

Edited by:

Joerg Heineke, University of Heidelberg, Germany

Reviewed by:

Izhak Kehat, Technion Israel Institute of Technology, Israel
Lukas Tombor, Goethe University Frankfurt, Germany
Christoph Dieterich, Heidelberg University, Germany

Copyright © 2022 Hegenbarth, Lezzoche, De Windt and Stoll. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Monika Stoll, bS5zdG9sbEBtYWFzdHJpY2h0dW5pdmVyc2l0eS5ubA==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.