Here, we present a protocol describing cell nuclei preparation. After microdissection and enzymatic dissociation of cardiac tissue into single cells, the progenitor cells were frozen, followed by isolation of pure viable cells, which were used for single-nucleus RNA sequencing and the single-nucleus assay for transposase-accessible chromatin with high-throughput sequencing analyses.
The developing heart is a complex structure containing various progenitor cells controlled by complex regulatory mechanisms. The examination of the gene expression and chromatin state of individual cells allows the identification of the cell type and state. Single-cell sequencing approaches have revealed a number of important characteristics of cardiac progenitor cell heterogeneity. However, these methods are generally restricted to fresh tissue, which limits studies with diverse experimental conditions, as the fresh tissue must be processed at once in the same run to reduce the technical variability. Therefore, easy and flexible procedures to produce data from methods such as single-nucleus RNA sequencing (snRNA-seq) and the single-nucleus assay for transposase-accessible chromatin with high-throughput sequencing (snATAC-seq) are needed in this area. Here, we present a protocol to rapidly isolate nuclei for subsequent single-nuclei dual-omics (combined snRNA-seq and snATAC-seq). This method allows the isolation of nuclei from frozen samples of cardiac progenitor cells and can be combined with platforms that use microfluidic chambers.
Among birth defects, congenital heart defects (CHDs) are the most common, occurring in about 1% of live births each year1,2. Genetic mutations are identified in only a minority of cases, implying that other causes, such as abnormalities in gene regulation, are involved in the etiology of CHD2,3. Cardiac development is a complex process of diverse and interacting cell types, making the identification of causal noncoding mutations and their effects on gene regulation challenging. Organogenesis of the heart begins with cellular progenitors that give rise to different subtypes of cardiac cells, including myocardial, fibroblast, epicardial, and endocardial cells4,5. Single-cell genomics is emerging as a key method for studying heart development and assessing the impact of cellular heterogeneity in health and disease6. The development of multi-omics methods for the simultaneous measurement of different parameters and the expansion of computational pipelines has facilitated the discovery of cell types and subtypes in the normal and diseased heart6. This article describes a reliable single-nucleus isolation protocol for frozen cardiac progenitor cells obtained from mouse embryos that is compatible with downstream snRNA-seq and snATAC-seq (as well as snRNA-seq and snATAC-seq combined)7,8,9.
ATAC-seq is a robust method that allows the identification of regulatory open chromatin regions and the positioning of nucleosomes10,11. This information is used to draw conclusions about the location, identity, and activity of transcription factors. The activity of chromatin factors, including remodelers, as well as the transcriptional activity of RNA polymerase, can, thus, be analyzed since the method is highly sensitive for measuring quantitative changes in chromatin structure1,2. Thus, ATAC-seq provides a robust and impartial approach to uncovering the mechanisms controlling transcriptional regulation in a specific cell type. ATAC-seq protocols have also been validated to measure chromatin accessibility in single cells, revealing variability in the chromatin architecture within cell populations10,12,13.
Although there have been notable advances in the field of single cells in recent years, the main difficulty is the processing of the fresh samples needed to perform these experiments14. To circumvent this difficulty, various tests have been carried out with the aim of conducting analyses such as snRNA-seq and snATAC-seq with frozen cardiac tissue or cells15,16.
Several platforms have been used to analyze single-cell genomics data17. The widely used platforms for single-cell gene expression and ATAC profiling are platforms for multiple microfluidic droplet encapsulation17. As these platforms use microfluidic chambers, debris or aggregates can clog the system, resulting in non-usable data. Thus, the success of single-cell studies depends on the accurate isolation of individual cells/nuclei.
The protocol presented here uses a similar approach to recent studies using snRNA-seq and snATAC-seq to understand congenital heart defects18,19,20,21,22,23. This procedure utilizes the enzymatic dissociation of freshly microdissected cardiac tissue followed by the cryopreservation of mouse cardiac progenitor cells. After thawing, the viable cells are purified and processed for nuclear isolation. In this work, this protocol was successfully used to obtain snRNA-seq and snATAC-seq data from the same nuclear preparation of mouse cardiac progenitor cells.
The animal procedure adopted in this study was approved by the animal ethics committees of the Aix-Marseille University (C2EA-14) and was carried out according to protocols approved by the appointed national ethical committee for animal experimentation (Ministère de l'Education Nationale, de l'Enseignement Supérieur et de la Recherche; Authorization Apafis N°33927-2021111715507212).
1. Setting up the timed mating prior to dissection
2. Tissue preparation and cell isolation
3. Cell counting and viability assessment
NOTE: The cell count and viability are critical parameters for the success of single-cell experiments. The snATAC-seq approach is sensitive to minor variations in cell number. Too few cells lead to over-digestion of the chromatin, which results in a higher number of reads and the mapping of inaccessible chromatin regions (noise). Similarly, having too many cells generates high-molecular-weight fragments that are hard to sequence. For the accurate determination of the cell and nuclei concentrations, the cell count and viability were assessed by two different methods.
4. Cell freezing
5. Cell thawing and removal of dead cells
6. Nuclei isolation
NOTE: snATAC and snRNA-seq combined are performed with a suspension of clean and intact nuclei. The optimization of the lysis conditions (lysis time and NP40 concentration) for the used cell type is recommended. The optimal lysis timepoint and detergent concentration are those that result in the maximum number of cells being lysed without disrupting the nuclear morphology. The loss of cells/nuclei can be reduced by using a swinging bucket centrifuge instead of a fixed-angle centrifuge. To minimize nuclei retention on plastics, it is recommended to use low-retention pipette tips and centrifuge tubes. This can increase the nuclei recovery. The coating of the pipette tips and centrifuge tubes with 5% BSA is a less expensive but more time-consuming alternative.
7. Quality and quantity assessments of the isolated nuclei
NOTE: Based on the initial cell count and estimating approximately 50% nuclei loss during cell lysis, resuspend the appropriate number of cells in cold diluted nuclei buffer. The calculation of the resuspension volume with chilled diluted nuclei buffer is based on the targeted nuclei recovery and the corresponding nuclei stock concentrations recommended by the manufacturer's user guide. See an example of this calculation in Supplementary Protocol 1.
8. Quality analysis of the snRNA-seq and snATAC-seq libraries
Compared to the preparation of single-cell suspensions for single-cell approaches, the preparation of single-nuclei suspensions is much more challenging and requires a higher degree of resolution and processing. The key factor for successful combined snRNA-seq and snATAC-seq is a clean and intact nuclei suspension. The protocol for efficient nuclei isolation must be adapted to each tissue type and condition (fresh or frozen). Here, an optimized protocol is described for the isolation of nuclei from frozen mouse embryonic cardiac cells. All the steps from embryonic cardiac region dissection and cardiac cell dissociation to nuclei isolation are summarized in Figure 1. For the starting material, a cell count of 1 x 105-1 x 106 cells and a cell viability >70% are recommended for efficient nuclei isolation. As shown in Figure 2, the additional step performed in this protocol to sort and remove the dead cells after thawing allowed for obtaining a cleaner cell suspension with a significantly higher cell viability and for improving the quality of the starting sample. In addition to nuclei isolation, this protocol provides detailed information on how to assess the quality and quantity of the isolated nuclei (Figure 3). The quality of the isolated nuclei greatly impacts the quality of the snRNA-seq and snATAC-seq combined data; here, the quality of the isolated nuclei was assessed by the evaluation of the nuclear membrane morphology. The isolated nuclei appeared smooth and uniformly round, as shown in Figure 3C, under brightfield microscopy, indicating an effective isolation process. For more accuracy, two different counting methods were used, including trypan blue and a kit for the fluorescence-based assessment of viability, to quantify the cells and nuclei. The fluorescence assay was used to determine the cell viability based on cellular integrity and intracellular esterase activity. This dye solution allows for distinguishing live cells from dead cells by simultaneously visualizing green fluorescent calcein AM, which indicates the intracellular esterase activity, and red fluorescent ethidium homodimer 1, which reflects the loss of cellular integrity. Using a fluorescent dye for counting helps distinguish nuclei from cell clumps and debris. Using this protocol, the cell viability passed from 80%-90% before nuclei isolation to <5% after nuclei isolation, as shown in Figure 3A,B.
Our procedure was compared to the existing method, which involves freezing fresh tissue directly in liquid nitrogen and performing the lysing step without prior enzymatic dissociation (Supplementary Protocol 2 and Supplementary Figure 1). Both approaches were tested to determine which method gives a better-quality nuclei suspension. Snap-freezing whole embryonic cardiac tissue gave rise to a high percentage of non-lysed cells detected as live, indicating inappropriate cell lysis (Supplementary Figure 1A). Aggregates and non-individual cells were observed despite filtration through filters (Supplementary Figure 1A,B).
The isolated nuclei were used to generate libraries for joint snATAC-seq and snRNA-seq. Figure 4 illustrates the quality control results of the cDNA used for the construction of the gene expression (GEX) library, and the ATAC library. The quality control was performed using automated electrophoresis (Figure 5A). The mRNA-derived cDNAs were quantified, and the yield was sufficient for subsequent use in GEX library construction. A high-quality GEX library ranging from ~300 bp to 600 bp and a high-quality ATAC library with periodic nucleosome winding ranging from 200 bp to 7 kbp were obtained. These libraries were sequenced by high-throughput sequencing and successfully generated single-nucleus gene expression and chromatin accessibility (Figure 5A). These raw data were ready to be demultiplexed and bioinformatically analyzed to generate transcriptional and epigenetic profiling of mouse E9.5 cardiac progenitor cells. SnRNA-seq and snATAC-seq were clustered, identified, and labeled based on known marker genes using unsupervised clustering with the Seurat toolkit for single genomics24 . This initial analysis confirmed the presence of all the expected cardiac cell types at E9.5 (Figure 5B).
All of the above results support the success of this protocol in isolating nuclei that are suitable for downstream single nuclei approaches from frozen embryonic cardiac cells.
Figure 1: Representation of the main steps in the sample preparation for combined snRNA-seq and snATAC-seq profiling. This flowchart recapitulates all the steps, including cardiac tissue dissection and cardiac cell dissociation, cell freezing and thawing, live cell purification by removing dead cells, and nuclei isolation, that are carried out for the simultaneous detection of mRNA and chromatin accessibility from the same cell. Abbreviations: PBS = phosphate-buffered saline; BSA = bovine serum albumin; RT = room temperature. Created with BioRender.com Please click here to view a larger version of this figure.
Figure 2: Viability of the optimized thawed cells after dead cell sorting. Images showing the cell count and viability using an automated cell counter before (top) and after (bottom) dead cell removal using the trypan blue exclusion dye. Please click here to view a larger version of this figure.
Figure 3: Quality control of the nuclei preparation. (A,B) Images showing the cell count and viability using an automated cell counter before (top) and after (bottom) nuclei isolation. Two counting methods were used: (A) the fluorescent mortality assay and (B)the trypan blue assay. (C) Images showing the quality of the nuclei isolation. Images were taken at 20x (scale bar = 50 µm) using brightfield and fluorescence microscopy. The brightfield image (left) shows nuclear morphology; the RFP (middle) shows cells that lost their membrane integrity, in other words, isolated nuclei; and the GFP (right) that normally indicates the esterase activity of live cells here confirms the absence of live cells in the nuclei suspension. Please click here to view a larger version of this figure.
Figure 4: Quality control of the single-cell gene expression and ATAC libraries. DNA analysis with automated electrophoresis showing the quality and size distribution of the cDNA (top), GEX library (middle), and ATAC library (bottom). For the cDNA quantification, the zone ranging from 200 bp to 8,900 bp was selected, and the bioanalyzer estimated the cDNA concentration (in pg/µL; red circle). A sufficient cDNA yield was obtained to proceed with GEX library construction. Please click here to view a larger version of this figure.
Figure 5: Sample performance assessed with the Cell Ranger software25. (A) Cross-sensitivity plot showing the transposition events in peaks per barcode on the x-axis and RNA unique molecular identifiers (UMIs) per barcode on the y-axis. As expected, the cells are mainly located in the upper-right-hand corner, with high RNA and ATAC data. A clear separation between the cells and non-cells (empty GEMs, background signal) was observed. (B) Representative uniform manifold approximation and projection (UMAP) plot of all the captured cells in scRNA-seq (left) and scATAC-seq (right) colored by cluster identity and grouped based on cell types. A good separation of clusters was observed. Joint GEX and ATAC raw data were obtained from sequencing 7,000 nuclei from the dissected E9.5 mouse embryonic heart region. Please click here to view a larger version of this figure.
Supplementary Figure 1: Quality control of nuclei isolation from the existing method involving freezing fresh tissue directly in liquid nitrogen. (A) Images showing the isolated nuclei count using an automated cell counter and trypan blue as the exclusion dye before (top) and after (bottom) filtration through a 40 µm strainer. A high percentage of non-lysed cells were detected as live. The detection of 20% live cells indicates inappropriate cell lysis. The black arrows indicate aggregates. (B) Images showing the quality of the isolated nuclei. The images were taken at 20x (top and middle with a scale bar of 50 µm) and 40x (bottom with a scale bar of 25 µm) using brightfield and fluorescence microscopy. The brightfield image (left) shows the nuclear morphology; the RFP (right) shows cells that lost their membrane integrity, in other words, isolated nuclei. The black arrows indicate multiplets of nuclei attached by debris. Please click here to download this File.
Table 1: Table of buffers and solutions used in the procedure. Please click here to download this Table.
Supplementary Protocol 1: Assessment of the quantity of the isolated nuclei. Please click here to download this File.
Supplementary Protocol 2: Isolation of nuclei from flash-frozen tissue. Please click here to download this File.
The analysis of the cellular composition of the developing heart by combined snRNA-seq and snATAC-seq studies provides a deeper understanding of the origin of congenital heart disease26. Several research laboratories have studied the effects of cardiac tissue cryopreservation on snRNA-seq27. Conducting snRNA-seq and snATAC-seq using fresh micro-dissected tissue from mouse models of human disease can be logistically challenging when comparing different experimental conditions. For a given stage of development, one difficulty is that the control and experimental groups must be processed in the same run to reduce the technical variability introduced by sample processing. When processing fresh tissue is not possible, it is recommended, based on our results, to dissociate and then cryopreserve rather than snap-freeze the whole embryonic cardiac tissue. The choice of one method over the other is based on the fact that if the cardiac tissue is frozen whole, isolating viable single cells after thawing is much more challenging. However, if the tissue dissociation protocol is not optimized, then it may be better to snap-freeze the tissue whole. Dissociation followed by cryopreservation is preferred given that this tissue dissociation protocol yields >90% fresh viable cells with no cell-type bias.
This approach has several benefits, such as the use of frozen cells, which eliminates the requirement for fresh tissue, and the use of a single-cell dissociation step, which allows for the isolation of intact single nuclei; indeed, this is critical for downstream single-cell analysis. However, the enzymatic dissociation may induce transcriptional bias in cell populations. Enzymatic cell dissociation interferes with the cellular microenvironment and, thus, induces transcriptional changes in most cell types28. The primary determinant of these changes is the dissociation time, not the cell type29. To moderate dissociation-induced artifacts it is important to follow the recommended digestion time. The published literature in diverse tissues indicates that a digestion time of 30 min or less induces minor transcriptional changes28,30. In addition, bioinformatics interventions could be used to minimize data distortions31,32.
Obtaining pure and intact single nuclei is the most important factor in single-nucleus genomics. The clogging of the microfluidic chambers of single-cell platforms due to the presence of cell clusters or organic debris leads to poor data quality or, more problematically, to experimental failure. In this procedure, existing protocols were modified by combining enzymatic dissociation, cryopreservation, and purification of viable cardiac progenitor cells. Our procedure allows the isolation of pure nuclei without the need for FACS sorting. Using this procedure, we obtained a high-quality nuclear suspension without clumping, almost no debris, and intact single nuclei. The cell isolation step is crucial to avoid potential nuclei clumping; indeed, this can be avoided by gentle pipetting using large pipette tips or by filtering the nuclei. The use of filters at different steps and subsequent observation under a light microscope can prevent such issues.
The protocol described here was successfully used for snRNA-seq and snATAC-seq profiling from the same nuclei preparation obtained from E9.5 cardiac progenitor cells.This method can be used in pregnant wild-type mice versus mouse models for congenital heart defects. Congenital heart defects likely occur due to the disruption of pathways that control the differentiation, multiplication, migration, and survival of cardiac progenitor cells33,34. To better understand the origin of these defects, this study focused specifically on embryonic stage 9.5, in which the progenitors of the major areas affected in congenital heart defects are present, in particular the cardiac progenitors of the second heart field11 and the neural crest35. The application of this procedure to cells of the adult mouse heart and human heart requires the optimization of the protocol, particularly in terms of the cell dissociation and lysis steps.
The main considerations between manual and automated cell counting are the costs, labor, and accuracy. Automated counters are faster than most skilled scientists performing a manual count. Especially when dealing with a large number of samples, it is of great benefit to be able to simply press count and view the result. This is an important point because once the nuclei are isolated, they must be processed quickly to avoid the loss of their integrity. The primary accuracy benefit of automated systems is that these systems eliminate variability between users or labs. Another advantage is that automated systems typically have a larger field of view than hemocytometers. When the number of cells/nuclei or the concentration of the cells/nuclei are lower, this larger field of view is certainly an advantage over the traditional manual method.
A high percentage of mitochondrial DNA contamination can be observed when the concentration of the non-ionic detergent is too high. By decreasing its concentration in the lysis buffer, this contamination can be reduced without decreasing the yield of the nuclei. The most satisfactory results were obtained in our experiments with a lysis buffer containing 0.1% detergent. Spinning the nuclei in a swing bucket can also reduce the number of mitochondria in the pellet.
Having an appropriate ratio of Tn5 transposase to nuclei is crucial for successful snATAC-seq13,14. For example, having too much Tn5 may lead to a high background due to closed chromatin and low complexity of the sequencing libraries, whereas sub-lysis may not result in a complete PCR amplified library13. To determine the number of nuclei accurately, we used two complementary methods.
Multimodal omics datasets provide the opportunity to investigate multiple levels of genomic organization simultaneously9,36. The joint RNA and ATAC multimodal approach enables the study of upstream regulators and downstream metabolic genes and provides a comprehensive approach for the study of transcriptional networks and chromatin architecture in the context of heart development at the single-cell resolution. This protocol for single-nuclei isolation is compatible with both the individual and joint evaluation of expression and chromatin datasets.Chromatin accessibility is an important epigenetic level of regulation because it enables the activity of transcriptional regulators at specific genomic sites8,9,10,11,12,13. A multitude of information on the identity and target sites of dozens of possible transcription factors are, thus, provided by the DNA sequence in the chromatin profile. The comparison of the information obtained by snATAC-seq with snRNA expression profiling can help to identify the most relevant transcription factors, which can then be further analyzed by Cut&Run37. We hope that this procedure will help interested researchers and encourage them to use this powerful method for their research.
The authors have nothing to disclose.
This research was supported by ERA-CVD-2019 and ANR-JCJC-2020 to SS. We thank the Genomics and Bioinformatics facility (GBiM) from the U 1251/Marseille Medical Genetics lab and the anonymous reviewers for providing valuable comments.
2100 Bioanalyzer Instrument | Agilent | No catalog number | |
5M Sodium chloride (NaCl) | Sigma | 59222C-500ML | |
BSA 10% | Sigma | A1595-50ML | |
Chromium Next GEM Chip J Single Cell Kit, 16 rxns | 10X Genomics | 1000230 | |
Chromium Next GEM Single Cell Multiome ATAC + Gene Expression Reagent Bundle, 4 rxns (including Nuclei Buffer 20X) | 10X Genomics | 1000285 | |
Countess cell counting chamber slides | Invitrogen | C10283 | |
Countess III FL | Thermofisher | No catalog number | |
Digitonin (5%) | Thermofisher | BN2006 | |
DMSO | Sigma | D2650-5x5ML | |
DNA LoBind Tubes | Eppendorf | 22431021 | |
D-PBS | Thermofisher | 14190094 | Sterile and RNase-free |
Dual Index Kit TT Set A 96 rxns | 10X Genomics | 1000215 | |
Falcon 15 mL Conical Centrifuge Tubes | Fisher Scientific | 352096 | |
Falcon 50 mL Conical Centrifuge Tubes | Fisher Scientific | 10788561 | |
HI-FBS | Thermofisher | A3840001 | Heat inactivated |
High sensitivity DNA kit | Agilent | 5067-4626 | |
Igepal CA-630 | Sigma | I8896-50ML | |
LIVE/DEAD Viability/Cytotoxicity Kit | Thermofisher | L3224 | |
MACS Dead Cell Removal kit: Dead Cell Romoval MicroBeads, Binding Buffer 20X | Miltenyi Biotec | 130-090-101 | |
MACS SmartStrainers (30 µm) | Miltenyi Biotec | 130-098-458 | |
Magnesium chloride (MgCl2) | Sigma | M1028-100ML | |
Milieu McCoy 5A | Thermofisher | 16600082 | |
MS Columns | Miltenyi Biotec | 130-042-201 | |
NovaSeq 6000 S2 | Illumina | No catalog number | |
Penicillin Streptomycin (Pen/Strep) | Thermofisher | 15070063 | |
PluriStrainer Mini 40µm | PluriSelect | V-PM15-2021-12 | |
Rock inhibitor | Enzo Life Sciences | ALX-270-333-M005 | |
Single Index Kit N Set A, 96 rxn | 10X Genomics | 1000212 | |
Standard 90mm Petri dish Sterilin | Thermofisher | 101R20 | |
Sterile double-distilled water | Thermofisher | R0582 | |
Trizma Hydrochloride solution (HCl) | Sigma | T2194-100ML | |
Trypan Blue stain (0.4%) | Invitrogen | T10282 | |
Trypsin 0.05% – EDTA 1X | Thermofisher | 25300054 | |
Tween20 | Sigma | P9416-50ML | |
Wide orifice filtered pipette tips 200 μl | Labcon | 1152-965-008-9 | |
ZEISS SteREO Discovery.V8 | ZEISS | No catalog number |